A Second Order Weighted Numerical Scheme on Nonuniform Meshes for Convection Diffusion Parabolic Problems

Authors

  • Lolugu Govindarao Department of Mathematics, National Institute of Technology Rourkela – 769008, India
  • Jugal Mohapatra Department of Mathematics, National Institute of Technology Rourkela – 769008, India

DOI:

https://doi.org/10.13052/ejcm2642-2085.2854

Keywords:

Singular perturbation, parabolic problems, boundary layer, monotone/weighted schemes, uniform convergence, implicit schemes

Abstract

In this article, a singularly perturbed parabolic convection-diffusion equation on a rectangular domain is considered. The solution of the problem possesses regular boundary layer which appears in the spatial variable. To discretize the time derivative, we use two type of schemes, first the implicit Euler scheme and second the implicit trapezoidal scheme on a uniform mesh. For approximating the spatial derivatives, we use the monotone hybrid scheme, which is a combination of midpoint upwind scheme and central difference scheme with variable weights on Shishkin-type meshes (standard Shishkin mesh, Bakhvalov-Shishkin mesh and modified Bakhvalov-Shishkin mesh). We prove that both numerical schemes converge uniformly with respect to the perturbation parameter and are of second order accurate. Thomas algorithm is used to solve the tri-diagonal system. Finally, to support the theoretical results, we present a numerical experiment by using the proposed methods.

Downloads

Download data is not yet available.

Author Biographies

Lolugu Govindarao, Department of Mathematics, National Institute of Technology Rourkela – 769008, India

Lolugu Govindarao received M.Sc. and M.Phil from Pondichery University, India. Currently he is pursuing Ph.D. at Department of Mathematics, National Institute of Technology Rourkela, India. His research interests include numerical methods for singular perturbation problems.

Jugal Mohapatra, Department of Mathematics, National Institute of Technology Rourkela – 769008, India

Jugal Mohapatra Ph.D. from Indian Institute of Technology Guwahati, India. He is currently working at Department of Mathematics, National Institute of Technology Rourkela. His research interest is Numerical Analysis mainly finite difference methods for singularly perturbed differential equations.

References

Andreev, V. B. & Kopteva, N. V. (1998). On the convergence, uniform with respect

to a small parameter, of monotone three- point finite difference approximation.

Journal of Difference Equations, 34, pp. 921–929.

Bansal, K., Rai, P. & Sharma, K. K. (2017). Numerical treatment for the class of time

dependent singularly perturbed parabolic problems with general shift arguments.

Differential Equations and Dynamical Systems, 25(2), pp. 327–346.

Brandt, A. & Yavneh, I. (1991). Inadequacy of first-order upwind difference schemes

for some recirculating flows. Journal of Computational Physics, 93, pp. 128–143.

Clavero, C., Gracia, J. L. & Jorge, J. C. (2005). High-order numerical methods

for one-dimensional parabolic singularly perturbed problems with regular layers.

Numerical Methods Partial Differential Equations, 21, pp. 148–169.

Clavero, C., Jorge, J. C. & Lisbona, F. (2003). A uniformly convergent scheme

on a nonuniform mesh for convection-diffusion parabolic problems. Journal of

Computational Methods in Applied Mathematics, 154, pp. 415–429.

Farrell, P. A., Hegarty, A. F., Miller, J. J. H., O’Riordan, E. & Shishkin, G. I. (2000)

Robust Computational Techniques for Boundary Layers. Chapman & Hall/CRC

Press, Boca Raton, FL.

Gartland, E. C. (1988). An analysis of a uniformly convergence finite difference or

finite element scheme for a model singular-perturbation problem. Mathematics of

Computation, 51(183), pp. 93–106.

Gartland, E. C. (1988a). Graded-mesh difference schemes for singularly perturbed

two-point boundary value problem. Mathematics of Computation, 51(184),

pp. 631–657.

Guo, W. & Stynes, M. (1993). Finite element analysis of exponentialy fitted

lumped schemes for time-dependent convection-diffusion problems. Numerische

Mathematik, 66, pp. 347–371.

Kay, B. D. & Elrick, D. E. (1967). Absorption and movement of lindane in soils.

Soil Science, 104, pp. 314–22.

Kellogg, R. B. & Tsan, A. (1978). Analysis of some diffenece approximations

for a singularly perturbed problem without turning points. Mathematics of

Computation, 32, pp. 1025–1039.

Kumar, M. & Sekhara Rao, S. C. (2010). High order parameter-robust numerical

method for time dependent singularly perturbed reaction-diffusion problems.

Computing, 90(1–2), pp. 15–38.

Linß, T., Roos, H.-G. & Vulanovi´c, R. (2000). Uniform pointwise convergence

on Shishkin-type meshes for quasi-linear convection-diffusion problems. SIAM

Journal on Numerical Analysis, 38, pp. 897–912.

Linß, T. (2001). Sufficient conditions for uniform convergence on layer-adapted

grids. Applied Numerical Mathematics, 37, pp. 241–255.

Linß, T. (2001). Uniform second-order pointwise convergence of a finite difference

discretization for a quasilinear problem. Zhurnal Vychislitel’noi Matematiki i

Matematicheskoi Fiziki, 41(6), pp. 947–958.

Lunardi, A. (1995). Analytic Semigroups and Optimal Regularity in Parabolic

Problems. Progress in nonlinear differential equations and their applications, 16.

Birkhäuser, Basel.

Kopteva, N. & Linß, T. (2001). Uniform second-order pointwise convergence of a

central difference approximation for a quasilinear convection-diffusion problem.

Journal of Computational Methods in Applied Mathematics, 137, pp. 257–267.

Melenk, J. M., Xenophontos, C. & Oberbroeckling, L. (2013). Robust exponential

convergence of hp FEM for singularly perturbed reaction-diffusion systems with

multiple scales. IMA Journal of Numerical Analysis, 33(2), pp. 609–628.

Miller, J. J. H., O’Riordan, E. & Shishkin, G. I. (2012). Fitted Numerical Methods

for Singular Perturbation Problems, Revised Edition.World Scientific, Singapore.

Misra, C. & Mishra, B. K. (1977). Miscible displacement of nitrate and chloride

under field conditions. Soil Science Society of America Journal, 41, pp. 496–499.

Mukherjee, K. & Natesan, S. (2009). Parameter-uniform hybrid numerical scheme

for time-dependent convection-dominated initial-boundary-value problems. Computing,

, pp. 209–230.

Ng-Stynes, M. J., O’Riordan, E. & Stynes, M. (1988). Numerical methods for timedependent

convection-diffusion equations. Journal of Computational Methods in

Applied Mathematics, 21, pp. 289–310.

O’Riordan, E. & Stynes, M. (1991). A globally uniformly convergent finite element

method for a singularly perturbed elliptic problem in two dimensions.

Mathematics of Computation, 51, pp. 47–62.

Palencia, C. (1993). A stability result for sectorical operators in Banach spaces.

SIAM Journal on Numerical Analysis, 30, 1373–1384.

Raji Reddy, N. & Mohapatra, J. (2015). An efficient numerical method for singularly

perturbed two point boundary value problems exhibiting boundary Layers.

National Academy Science Letters, 38(4), pp. 355–359.

Salama, A. A. & Al-Amerya, D. G. (2017). A higher order uniformly convergent

method for singularly perturbed delay parabolic partial differential equations.

International Journal of Computer Mathematics, 12(94), 2520–2546.

Schwab, C., Suri, M. & Xenophontos, C. A. (1996). Boundary layer approximation

by spectral/hp methods. Houston Journal of Mathematics Spec. Issue of

ICOSAHOM, 95, pp. 501–508.

Shishkin, G. I. & Shishkina, L. P. (2009). Difference Methods for Singular

Perturbation Problems. CRC Press, Boca Raton.

Stynes, M. & Roos, H. G. (1997). The midpoint upwind scheme. Applied Numerical

Mathematics, 23, pp. 361–374.

Vulanovi´c, R. (1989). A uniform method for quasilinear singular perturbation

problems without turning points. Computing, 4(1), pp. 97–106.

Vulanovi´c, R. (2001). A priori meshes for singularly perturbed quasilinear twopoint

boundary value problems. IMA Journal of Numerical Analysis, 21(1),

pp. 349–366.

Xenophontos, C. (1998). The h-p finite element method for singularly perturbed

boundary value problems on smooth domains. Math. Models and Methods in

Applied Sciences, 8, pp. 69–79.

Xenophontos, C. (1998). The h-p finite element method for singularly perturbed

problems on non-smooth domains. Numerical Methods for Partial Differential

Equations, 16, 1–28.

Downloads

Published

2020-01-14

How to Cite

Govindarao, L., & Mohapatra, J. (2020). A Second Order Weighted Numerical Scheme on Nonuniform Meshes for Convection Diffusion Parabolic Problems. European Journal of Computational Mechanics, 28(5), 467–498. https://doi.org/10.13052/ejcm2642-2085.2854

Issue

Section

Original Article