Abstract
In this paper, we construct a family of modified Patankar Runge–Kutta methods, which is conservative and unconditionally positivity-preserving, for production–destruction equations, and derive necessary and sufficient conditions to obtain second-order accuracy. This ordinary differential equation solver is then extended to solve a class of semi-discrete schemes for PDEs. Combining this time integration method with the positivity-preserving finite difference weighted essentially non-oscillatory (WENO) schemes, we successfully obtain a positivity-preserving WENO scheme for non-equilibrium flows. Various numerical tests are reported to demonstrate the effectiveness of the methods.
Similar content being viewed by others
References
Burchard, H., Deleersnijder, E., Meister, A.: A high-order conservative Patankar-type discretisation for stiff systems of production-destruction equations. Appl. Numer. Math. 47(1), 1–30 (2003)
Chertock, A., Cui, S., Kurganov, A., Wu, T.: Steady state and sign preserving semi-implicit Runge–Kutta methods for ODEs with stiff damping term. SIAM J. Numer. Anal. 53(4), 2008–2029 (2015)
Chertock, A., Cui, S., Kurganov, A., Wu, T.: Well-balanced positivity preserving central-upwind scheme for the shallow water system with friction terms. Int. J. Numer. Methods Fluids 78(6), 355–383 (2015)
Formaggia, L., Scotti, A.: Positivity and conservation properties of some integration schemes for mass action kinetics. SIAM J. Numer. Anal. 49(3), 1267–1288 (2011)
Gottlieb, S., Shu, C.-W., Tadmor, E.: Strong stability-preserving high-order time discretization methods. SIAM Rev. 43(1), 89–112 (2001)
Hairer, E., Wanner, G.: Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems. Springer, Heidelberg (1996)
Hu, J., Shu, R., Zhang, X.: Asymptotic-preserving and positivity-preserving implicit–explicit schemes for the stiff BGK equation. SIAM J. Numer. Anal. 56(2), 942–973 (2018)
Huang, J., Shu, C.-W.: A second-order asymptotic-preserving and positivity-preserving discontinuous Galerkin scheme for the Kerr–Debye model. Math. Models Methods Appl. Sci. 27(03), 549–579 (2017)
Huang, J., Shu, C.-W.: Bound-preserving modified exponential Runge–Kutta discontinuous Galerkin methods for scalar hyperbolic equations with stiff source terms. J. Comput. Phys. 361, 111–135 (2018)
Kopecz, S., Meister, A.: On order conditions for modified Patankar–Runge–Kutta schemes. Appl. Numer. Math. 123, 159–179 (2018)
Kopecz, S., Meister, A.: Unconditionally positive and conservative third order modified Patankar–Runge–Kutta discretizations of production–destruction systems. BIT Numer. Math. 58, 691–728 (2018). https://doi.org/10.1007/s10543-018-0705-1
Liu, X.-D., Osher, S., Chan, T.: Weighted essentially non-oscillatory schemes. J. Comput. Phys. 115(1), 200–212 (1994)
Meister, A., Ortleb, S.: On unconditionally positive implicit time integration for the DG scheme applied to shallow water flows. Int. J. Numer. Methods Fluids 76(2), 69–94 (2014)
Patankar, S.: Numerical Heat Transfer and Fluid Flow. CRC Press, London (1980)
Shu, C.-W.: Essentially Non-oscillatory and Weighted Essentially Non-oscillatory Schemes for Hyperbolic Conservation Laws, pp. 325–432. Springer, Berlin (1998)
Shu, C.-W., Osher, S.: Efficient implementation of essentially non-oscillatory shock-capturing schemes. J. Comput. Phys. 77(2), 439–471 (1989)
Strang, G.: On the construction and comparison of difference schemes. SIAM J. Numer. Anal. 5(3), 506–517 (1968)
Wang, C., Zhang, X., Shu, C.-W., Ning, J.: Robust high order discontinuous Galerkin schemes for two-dimensional gaseous detonations. J. Comput. Phys. 231(2), 653–665 (2012)
Wang, R., Spiteri, R.J.: Linear instability of the fifth-order WENO method. SIAM J. Numer. Anal. 45(5), 1871–1901 (2007)
Wang, W., Shu, C.-W., Yee, H., Sjögreen, B.: High-order well-balanced schemes and applications to non-equilibrium flow. J. Comput. Phys. 228(18), 6682–6702 (2009)
Xing, Y., Zhang, X.: Positivity-preserving well-balanced discontinuous Galerkin methods for the shallow water equations on unstructured triangular meshes. J. Sci. Comput. 57(1), 19–41 (2013)
Xing, Y., Zhang, X., Shu, C.-W.: Positivity-preserving high order well-balanced discontinuous Galerkin methods for the shallow water equations. Adv. Water Resour. 33(12), 1476–1493 (2010)
Zhang, X.: On positivity-preserving high order discontinuous Galerkin schemes for compressible navier-stokes equations. J. Comput. Phys. 328, 301–343 (2017)
Zhang, X., Shu, C.-W.: On maximum-principle-satisfying high order schemes for scalar conservation laws. J. Comput. Phys. 229(9), 3091–3120 (2010)
Zhang, X., Shu, C.-W.: On positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations on rectangular meshes. J. Comput. Phys. 229(23), 8918–8934 (2010)
Zhang, X., Shu, C.-W.: Positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations with source terms. J. Comput. Phys. 230(4), 1238–1248 (2011)
Zhang, X., Shu, C.-W.: Positivity-preserving high order finite difference WENO schemes for compressible Euler equations. J. Comput. Phys. 231(5), 2245–2258 (2012)
Zhang, Y., Zhang, X., Shu, C.-W.: Maximum-principle-satisfying second order discontinuous Galerkin schemes for convection–diffusion equations on triangular meshes. J. Comput. Phys. 234, 295–316 (2013)
Acknowledgements
We would like to thank Xiangxiong Zhang from Purdue University and Tao Xiong from Xiamen University for many fruitful discussions.
Author information
Authors and Affiliations
Corresponding author
Additional information
Research supported by ARO Grant W911NF-15-1-0226 and NSF Grant DMS-1719410.
Rights and permissions
About this article
Cite this article
Huang, J., Shu, CW. Positivity-Preserving Time Discretizations for Production–Destruction Equations with Applications to Non-equilibrium Flows. J Sci Comput 78, 1811–1839 (2019). https://doi.org/10.1007/s10915-018-0852-1
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10915-018-0852-1