A parameter-robust finite difference method for singularly perturbed delay parabolic partial differential equations

A Dirichlet boundary value problem for a delay parabolic differential equation is studied on a rectangular domain in the x-t plane. The second-order space derivative is multiplied by a small singular perturbation parameter, which gives rise to parabolic boundary layers on the two lateral sides of the rectangle. A numerical method comprising a standard finite difference operator (centred in space, implicit in time) on a rectangular piecewise uniform fitted mesh of N"xxN"t elements condensing in the boundary layers is proved to be robust with respect to the small parameter, or parameter-uniform, in the sense that its numerical solutions converge in the maximum norm to the exact solution uniformly well for all values of the parameter in the half-open interval (0,1]. More specifically, it is shown that the errors are bounded in the maximum norm by C(N"x^-^2ln^2N"x+N"t^-^1), where C is a constant independent not only of N"x and N"t but also of the small parameter. Numerical results are presented, which validate numerically this theoretical result and show that a numerical method consisting of the standard finite difference operator on a uniform mesh of N"xxN"t elements is not parameter-robust.

[1]  Shagi-Di Shih.  On a Class of Singularly Perturbed Parabolic Equations , 2001 .

[2]  Jianhong Wu.  Theory and Applications of Partial Functional Differential Equations , 1996 .

[3]  Ali R. Ansari,et al.  Numerical solution of a convection diffusion problem with Robin boundary conditions , 2003 .

[4]  Lutz Tobiska,et al.  Numerical Methods for Singularly Perturbed Differential Equations , 1996 .

[5]  Ali R. Ansari,et al.  A note on iterative methods for solving singularly perturbed problems using non-monotone methods on Shishkin meshes , 2003 .

[6]  L. P. Shishkina,et al.  FITTED MESH METHODS FOR PROBLEMS WITH PARABOLIC BOUNDARY LAYERS , 1998 .

[7]  J. C. Jorge,et al.  An alternating direction scheme on a nonuniform mesh for reaction-diffusion parabolic problems , 2000 .

[8]  G. I. Shishkin.  Approximation of the solutions of singularly perturbed boundary-value problems with a parabolic boundary layer , 1991 .

[9]  James D. Murray.  Mathematical Biology: I. An Introduction , 2007 .

[10]  Pieter W. Hemker,et al.  Novel Defect-correction High-order, in Space and Time, Accurate Schemes for Parabolic Singularly Perturbed Convection-diffusion Problems , 2003 .

[11]  JohnM . Miller,et al.  Robust Computational Techniques for Boundary Layers , 2000 .

[12]  Piet Hemker,et al.  ε-uniform schemes with high-order time-accuracy for parabolic singular perturbation problems , 2000 .

[13]  Ali R. Ansari,et al.  A note on fitted operator methods for a laminar jet problem , 2003 .

[14]  O. Ladyženskaja.  Linear and Quasilinear Equations of Parabolic Type , 1968 .

[15]  John J. H. Miller.  Fitted Numerical Methods for Singular Perturbation Problems: Error Estimates in the Maximum Norm for Linear Problems in One and Two Dimensions , 1996 .

[16]  Ulrich Behn,et al.  Inverted pendulum parametrically driven by dichotomous noise , 1996 .

[17]  Grigorii I. Shishkin,et al.  The Use of Defect Correction for the Solution of Parabolic Singular Perturbation Problems , 1995 .