The Midpoint Upwind Finite Difference Scheme for Time-Dependent Singularly Perturbed Convection-Diffusion Equations on Non-Uniform Mesh

被引:13
作者
Kadalbajoo, Mohan K. [1 ]
Awasthi, Ashish [2 ]
机构
[1] IIT Kanpur, Dept Math, Kanpur, Uttar Pradesh, India
[2] Natl Inst Technol, Calicut, Kerala, India
关键词
Singular perturbation; time dependent convection-diffusion equation; Backward Euler Method; midpoint upwind; Shishkin mesh; uniform convergence;
D O I
10.1080/15502287.2011.564264
中图分类号
O3 [力学];
学科分类号
08 ; 0801 ;
摘要
A numerical approach is proposed to examine the singularly perturbed time dependent convection-diffusion equation in one space dimension on a rectangular domain. The solution of considered problem exhibits a boundary layer on the right side of the domain. We semidiscretize the continuous problem by means of backward Euler finite difference method in the temporal direction. The semi-discretization process yields a set of ordinary differential equations at each time level. A resulting set of ordinary differential equations are discretized by using midpoint upwind finite difference scheme on a non-uniform mesh of Shishkin type. The resulting finite difference method is shown to be almost of second order accurate in the coarse mesh and almost of first order accurate in fine mesh in the spatial direction. First order accuracy is achieved in the temporal direction. An extensive amount of analysis has been carried out in order to obtain uniform convergence of the method. Finally, we have found that the method is uniformly convergent with respect to the singular perturbation parameter i.e. e-uniform. Some numerical experiments have been carried out to validate the predicted theory.
引用
收藏
页码:150 / 159
页数:10
相关论文
共 18 条
[1]   DIFFERENCE APPROXIMATIONS FOR SINGULAR PERTURBATIONS OF SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS [J].
ABRAHAMSSON, LR ;
KELLER, HB ;
KREISS, HO .
NUMERISCHE MATHEMATIK, 1974, 22 (05) :367-391
[2]  
Baumert H., 1981, LECT NOTES CONTROL I, V23, P482
[3]  
BOBISUD L, 1968, ARCH RATION MECH AN, V27, P385
[4]   A uniformly convergent scheme on a nonuniform mesh for convection-diffusion parabolic problems [J].
Clavero, C ;
Jorge, JC ;
Lisbona, F .
JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS, 2003, 154 (02) :415-429
[5]   An alternating direction scheme on a nonuniform mesh for reaction-diffusion parabolic problems [J].
Clavero, C ;
Jorge, JC ;
Lisbona, F ;
Shishkin, GI .
IMA JOURNAL OF NUMERICAL ANALYSIS, 2000, 20 (02) :263-280
[6]  
Ewing R.E, 1983, MATH RESERVOIR SIMUL
[7]  
Farrel PA, 2000, ROBUST COMPUTATIONAL
[8]   AN UPWIND FINITE-ELEMENT METHOD FOR ELECTROMAGNETIC-FIELD PROBLEMS IN MOVING-MEDIA [J].
HAHN, SY ;
BIGEON, J ;
SABONNADIERE, JC .
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, 1987, 24 (11) :2071-2086
[9]   ε-uniform schemes with high-order time-accuracy for parabolic singular perturbation problems [J].
Hemker, PW ;
Shishkin, GI ;
Shishkina, LP .
IMA JOURNAL OF NUMERICAL ANALYSIS, 2000, 20 (01) :99-121
[10]  
Hirsch C., 1988, NUMERICAL COMPUTATIO, V1