×

Positivity-preserving, flux-limited finite difference and finite element methods for reactive transport. (English) Zbl 1025.76029

Summary: A new class of positivity-preserving, flux-limited finite-difference and Petrov-Galerkin (PG) finite-element methods are devised for reactive transport problems.The methods are similar to classical TVD flux-limited schemes with the main difference being that the flux-limiter constraint is designed to preserve positivity for problems involving diffusion and reaction. In the finite-element formulation, we also consider the effect of numerical quadrature in the lumped and consistent mass matrix forms on the positivity-preserving property. Analysis of the latter scheme shows that positivity-preserving solutions of the resulting difference equations can only be guaranteed if the flux-limited scheme is both implicit and satisfies an additional lower-bound condition on time-step size. We show that this condition also applies to standard Galerkin linear finite-element approximations to the linear diffusion equation. Numerical experiments are provided to demonstrate the behavior of the methods and confirm the theoretical conditions on time-step size, mesh spacing, and flux limiting for transport problems with and without nonlinear reaction.

MSC:

76M20 Finite difference methods applied to problems in fluid mechanics
76M10 Finite element methods applied to problems in fluid mechanics
76V05 Reaction effects in flows
Full Text: DOI

References:

[1] Harten, Journal of Computational Physics 49 pp 357– (1983)
[2] Sweby, SIAM Journal on Numerical Analysis 21 pp 995– (1984)
[3] Numerical Methods for Conservation Laws. Birkhauser Verlag: Basal, 1992. · Zbl 0847.65053 · doi:10.1007/978-3-0348-8629-1
[4] Yee, Journal of Computational Physics 68 pp 151– (1987)
[5] Yee, Journal of Computational Physics 88 pp 31– (1990)
[6] Wang, Journal of Computational Physics 97 pp 53– (1991)
[7] Arora, Journal of Computational Physics 132 pp 3– (1997)
[8] Daru, Computers and Fluids 30 pp 89– (2001)
[9] Andrews, International Journal for Numerical Methods in Fluids 21 pp 205– (1995)
[10] Gupta, In-situ 15 pp 289– (1991)
[11] Blunt, Journal of Computational Physics 102 pp 194– (1992)
[12] Liu, Applied Mathematical Modelling 19 pp 162– (1995)
[13] Shi, International Journal for Numerical Methods in Fluids 23 pp 241– (1996)
[14] Hubbard, Journal of Computational Physics 155 pp 54– (1999)
[15] Suresh, SIAM Journal on Scientific Computing 22 pp 1184– (2000)
[16] Codina, Computer Methods in Applied Mechanics and Engineering 156 pp 185– (1998)
[17] Berzins, Communications in Numerical Methods in Engineering 17 pp 658– (2001)
[18] Christie, International Journal for Numerical Methods in Engineering 10 pp 1389– (1976)
[19] Biochemical Engineering Fundamentals (2nd edn). McGraw-Hill: New York, 1986.
[20] Chemical Kinetics and Reaction Mechanisms (2nd edn). McGraw-Hill: New York, 1995.
[21] Kindred, Water Resources Research 25 pp 1149– (1989)
[22] van Leer, Journal of Computational Physics 14 pp 361– (1974)
[23] Finite Elements Fluid Mechanics, vol. VI. Prentice-Hall: Englewood Cliffs, NJ, 1986.
[24] Wilcoxson, Journal of Computational Physics 115 pp 376– (1994)
[25] Groundwater Transport: Handbook of Mathematical Models.Water Resources Monograph Series 10, American Geophysical Union: Washington, DC, 1984. · doi:10.1029/WM010
[26] On the nonnegativity conservation of finite element solution of parabolic problems. In 3D Finite Element Conference Proceedings, et al. (eds). Gakkotosho Co.: Tokyo, 2001.
[27] Conduction of Heat in Solids (2nd edn). Oxford University Press: Oxford, UK, 1959.
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.