×

Studies of the accuracy of time integration methods for reaction-diffusion equations. (English) Zbl 1039.65069

Summary: We present numerical experiments of time integration methods applied to systems of reaction-diffusion equations. Our main interest is in evaluating the relative accuracy and asymptotic order of accuracy of the methods on problems which exhibit an approximate balance between the competing component time scales. Nearly balanced systems can produce a significant coupling of the physical mechanisms and introduce a slow dynamical time scale of interest. These problems provide a challenging test for this evaluation and tend to reveal subtle differences between the various methods.
The methods we consider include first- and second-order semi-implicit, fully implicit, and operator-splitting techniques. The test problems include a prototype propagating nonlinear reaction-diffusion wave, a non-equilibrium radiation-diffusion system, a Brusselator chemical dynamics system and a blow-up example. In this evaluation we demonstrate a “split personality” for the operator-splitting methods that we consider. While operator-splitting methods often obtain very good accuracy, they can also manifest a serious degradation in accuracy due to stability problems.

MSC:

65M60 Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs
35K57 Reaction-diffusion equations
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Bandle, C.; Brunner, H., Blowup in diffusion equation: a survey, Journal of Computational and Applied Mathematics, 97, 3-22 (1998) · Zbl 0932.65098
[2] (Bowers, R. L.; Wilson, J. R., Numerical Modeling in Applied Physics and Astrophysics (1991), Jones and Bartlett: Jones and Bartlett Boston) · Zbl 0786.76001
[3] (Brackbill, J.; Cohen, B., Multiple Time Scales (1985), Academic Press: Academic Press Orlando) · Zbl 0563.00023
[4] Brown, P. N.; Woodward, C. S., Preconditioning strategies or fully implicit radiation diffusion with material-energy transfer, SIAM Journal of Scientific Computing, 23, 2, 499-516 (2001) · Zbl 0992.65102
[5] Chung, J.; Hulbert, G. M., A time integration algorithm for structural dynamics with improved numerical dissipation: the generalized-\(α\) method, Journal of Applied Mechanics, 60, 371-375 (1993) · Zbl 0775.73337
[6] Cohen, S. D.; Hindmarsh, A. C., CVODE, a stiff/nonstiff ODE solver in C, Computers in Physics, 10, 2, 138-143 (1996)
[7] D.J. Estep, S. Tavener, Notes on determining stability of split-step integrators using green’s functions, personal communication; D.J. Estep, S. Tavener, Notes on determining stability of split-step integrators using green’s functions, personal communication
[8] Fletcher, C. A.J., Computational Techniques for Fluid Dynamics, vol. 1 (1988), Springer: Springer Berlin · Zbl 0706.76001
[9] Jansen, K. E.; Whiting, C. H.; Hulbert, G. M., A generalized-\(α\) method for integrating the filtered Navier-Stokes equations with a stabilized finite element method, Computer Methods in Applied Mechanics and Engineering, 190, 305-319 (2000) · Zbl 0973.76048
[10] Knio, O. M.; Najm, H. N.; Wyckoff, P. S., A semi-implicit numerical scheme of reacting flow: II. Stiff, operator-split formulation, Journal Computational Physics, 154, 428-467 (1999) · Zbl 0958.76061
[11] Knoll, D. A.; Chacon, L.; Margolin, L. G.; Mousseau, V. A., On balanced approximations for time integration of multiple time scale systems, Journal of Computational Physics, 185, 583-611 (2003) · Zbl 1047.76074
[12] Knoll, D. A.; Rider, W. J.; Olson, G. L., An efficient nonlinear solution methods for non-equilibrium radiation diffusion, Journal of Quantitative Spectroscopy and Radiative Transfer, 63, 15-29 (1999)
[13] Knoll, D. A.; Rider, W. J.; Olson, G. L., Nonlinear convergence, accuracy and time step control in non-equilibrium radiation diffusion, Journal of Quantitative Spectroscopy and Radiative Transfer, 65, 25-36 (2001)
[14] R.B. Lowrie, A comparison of time integration methods for nonlinear relaxation and diffusion, Journal of Computational Physics (accepted for publication); R.B. Lowrie, A comparison of time integration methods for nonlinear relaxation and diffusion, Journal of Computational Physics (accepted for publication) · Zbl 1053.65080
[15] G.I. Marchuk, On the theory of the splitting-up method, in: Proceedings of the 2nd Symposium on Numerical Solution of Partial Differential Equations, SVNSPADE, 1970, pp. 469-500; G.I. Marchuk, On the theory of the splitting-up method, in: Proceedings of the 2nd Symposium on Numerical Solution of Partial Differential Equations, SVNSPADE, 1970, pp. 469-500
[16] Mousseau, V. A.; Knoll, D. A.; Rider, W. J., Physics-based preconditioning and the Newton-Krylov method for non-equilibrium radiation diffusion, Journal of Computational Physics, 160, 743-765 (2000) · Zbl 0949.65092
[17] Nie, Y. Y.; Thomee, V., A lumped mass finite element method with quadrature for a non-linear parabolic problem, IMA Journal of Numerical Analysis, 5, 371-396 (1985) · Zbl 0591.65079
[18] C.C. Ober, J.N. Shadid, Studies on the accuracy of time-integration schemes for the radiation-diffusion equations, Journal of Computational Physics (2004), in press; C.C. Ober, J.N. Shadid, Studies on the accuracy of time-integration schemes for the radiation-diffusion equations, Journal of Computational Physics (2004), in press · Zbl 1053.65082
[19] Oran, E. S.; Boris, J. P., Numerical Simulation of Reactive Flow (2001), Cambridge University Press: Cambridge University Press Cambridge · Zbl 0762.76098
[20] Prigogine, I.; Lefever, R., Symmetry breaking instabilities in dissipative systems, ii, Journal of Chemical Physics, 48, 4, 1695-1700 (1968)
[21] L.A. Romero, On the accuracy of operator-splitting methods for problems with multiple time scales, Technical Report SAND2002-1448, Sandia National Laboratories, August 2002; L.A. Romero, On the accuracy of operator-splitting methods for problems with multiple time scales, Technical Report SAND2002-1448, Sandia National Laboratories, August 2002
[22] C. J. Roy, Grid convergence error analysis for mixed-order numerical schemes, AIAA Paper 2001-2606, 2001; C. J. Roy, Grid convergence error analysis for mixed-order numerical schemes, AIAA Paper 2001-2606, 2001
[23] Samarskii, A. A.; Galaktionov, V. A.; Kurdyumov, S. P.; Mikhailov, A. P., Blow-up in Quasilinear Parabolic Equations. Blow-up in Quasilinear Parabolic Equations, de Gruyter Expositions in Mathematics, vol. 19 (1995), Walter de Gruyter, Inc: Walter de Gruyter, Inc Berlin · Zbl 1020.35001
[24] J.N. Shadid, A.G. Salinger, R.C. Schmidt, T.M. Smith, S.A. Hutchinson, G.L. Hennigan, K.D. Devine, H.K. Moffat, A finite element computer program for reacting flow problems, Part 1: theoretical development, Technical Report SAND98-2864, Sandia National Laboratories, January 1999; J.N. Shadid, A.G. Salinger, R.C. Schmidt, T.M. Smith, S.A. Hutchinson, G.L. Hennigan, K.D. Devine, H.K. Moffat, A finite element computer program for reacting flow problems, Part 1: theoretical development, Technical Report SAND98-2864, Sandia National Laboratories, January 1999
[25] Shampine, L. F.; Reichelt, M. W., The matlab ode suite, SIAM Journal of Scientific Computing, 18, 1, 1-22 (1997) · Zbl 0868.65040
[26] Sportisee, B., An analysis of operator splitting techniques in the stiff case, Journal of Computational Physics, 161, 140-168 (2000)
[27] Strang, G., On the construction and comparison of difference schemes, SIAM Journal of Numerical Analysis, 5, 3, 506-517 (1968) · Zbl 0184.38503
[28] Szilard, R. H.; Pomraning, G. C., Numerical transport and diffusion methods in radiative transfer, Nuclear Science and Engineering, 112, 256-269 (1992)
[29] Thomée, V., Galerkin Finite Element Methods for Parabolic Problems (1997), Springer-Verlag: Springer-Verlag New York · Zbl 0884.65097
[30] Yanenko, N. N., The Method of Fractional Steps (1971), Springer-Verlag: Springer-Verlag New York, Translation Ed. M. Holt · Zbl 0209.47103
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.