A High-Order Finite Difference Method for 1D Nonhomogeneous Heat Equations

被引:7
作者
Lin, Yuan [1 ]
Gao, Xuejun [2 ]
Xiao, MingQing [1 ]
机构
[1] So Illinois Univ, Dept Math, Carbondale, IL 62901 USA
[2] Guangdong Univ Technol, Fac Appl Math, Guangzhou, Guangdong, Peoples R China
关键词
heat equations; finite difference method; Bartels-Stewart method; boundary value method; unconditional stability; DIFFUSION EQUATION; SCHEMES; ACCURACY;
D O I
10.1002/num.20345
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
In this article a sixth-order approximation method (in both temporal and spatial variables) for solving non-homogeneous heat equations is proposed. We first develop a sixth-order finite difference approximation scheme for a two-point boundary Value problem, and then heat equation is approximated by a system of ODEs defined on spatial grid points. The ODE system is discretized to a Sylvester matrix equation via boundary value method. The obtained algebraic system is solved by a modified Bartels-Stewart method. The proposed approach is unconditionally stable. Numerical results are provided to illustrate the accuracy and efficiency of our approximation method along with comparisons with those generated by the standard second-order Crank-Nicolson scheme as well as Sun-Zhang's recent fourth-order method. (C) 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 25: 327-346, 2009
引用
收藏
页码:327 / 346
页数:20
相关论文
共 26 条
[1]   A TAYLOR-SERIES APPROACH TO NUMERICAL ACCURACY AND A 3RD-ORDER SCHEME FOR STRONG CONVECTIVE FLOWS [J].
BRADLEY, D ;
MISSAGHI, M ;
CHIN, SB .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 1988, 69 (02) :133-151
[2]  
Brugnano L., 1998, Solving Differential Equations By Multistep Initial and Boundary Value Methods
[3]   THE SENSITIVITY AND ACCURACY OF 4TH-ORDER FINITE-DIFFERENCE SCHEMES ON NONUNIFORM GRIDS IN ONE-DIMENSION [J].
CASTILLO, JE ;
HYMAN, JM ;
SHASHKOV, MJ ;
STEINBERG, S .
COMPUTERS & MATHEMATICS WITH APPLICATIONS, 1995, 30 (08) :41-55
[4]   Strang-type preconditioners for systems of LMF-based ODE codes [J].
Chan, RH ;
Ng, MK ;
Jin, XQ .
IMA JOURNAL OF NUMERICAL ANALYSIS, 2001, 21 (02) :451-462
[5]  
Chawla M. M., 1998, Neural, Parallel & Scientific Computations, V6, P523
[6]  
DATTAGUPTA A, 1991, IN SITU, V15, P289
[7]   HIGHER-ORDER MONOTONIC CONVECTIVE DIFFERENCE SCHEMES [J].
FORESTER, CK .
JOURNAL OF COMPUTATIONAL PHYSICS, 1977, 23 (01) :1-22
[8]   SOLUTION OF THE SYLVESTER MATRIX EQUATION AXB(T)+CXD(T)=E [J].
GARDINER, JD ;
LAUB, AJ ;
AMATO, JJ ;
MOLER, CB .
ACM TRANSACTIONS ON MATHEMATICAL SOFTWARE, 1992, 18 (02) :223-231
[9]  
GARTLAND EC, 1982, SIAM J SCI STAT COMP, V3, P460, DOI 10.1137/0903030
[10]  
HARTEN A, 1987, J COMPUT PHYS, V71, P231, DOI [10.1016/0021-9991(87)90031-3, 10.1006/jcph.1996.5632]