The High-Order Variable-Coefficient Explicit-Implicit-Null Method for Diffusion and Dispersion Equations

Abstract

For the high-order diffusion and dispersion equations, the general practice of the explicit-implicit-null (EIN) method is to add and subtract an appropriately large linear highest derivative term with a constant coefficient at one side of the equation, and then apply the standard implicit-explicit method to the equivalent equation. We call this approach the constant-coefficient EIN method in this paper and hereafter denote it by “CC-EIN”. To reduce the error in the CC-EIN method, the variable-coefficient explicit-implicit-null (VC-EIN) method, which is obtained by adding and subtracting a linear highest derivative term with a variable coefficient, is proposed and studied in this paper. Coupled with the local discontinuous Galerkin (LDG) spatial discretization, the VC-EIN method is shown to be unconditionally stable and can achieve high order of accuracy for both one-dimensional and two-dimensional quasi-linear and nonlinear equations. In addition, although the computational cost slightly increases, the VC-EIN method can obtain more accurate results than the CC-EIN method, if the diffusion coefficient or the dispersion coefficient has a few high and narrow bumps and the bumps only account for a small part of the whole computational domain.

Keywords

Diffusion equation / Dispersion equation / Stability / Explicit-implicit-null (EIN) time discretization / Local discontinuous Galerkin (LDG) method / 65M60 / 65M12

References

[1] Ascher UM, Ruuth SJ, Spiteri RJ. Implicit-explicit Runge-Kutta methods for time-dependent partial differential equations. Appl. Numer. Math.. 1997, 25: 151-167

[2] Bassi F, Botti L, Colombo A, Ghidoni A, Massa F. Linearly implicit Rosenbrock-type Runge-Kutta schemes applied to the Discontinuous Galerkin solution of compressible and incompressible unsteady flows. Comput. Fluids. 2015, 118: 305-320

[3] Benjamin TB, Bona JL, Mahony JJ. Model equations for long waves in nonlinear dispersive systems. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci.. 1972, 272: 47-78

[4] Boscarino S, LeFloch P-G, Russo G. On a class of uniformly accurate IMEX Runge-Kutta schemes and application to hyperbolic systems with relaxation. SIAM J. Sci. Comput.. 2014, 36: 377-395

[5] Cavaliere P, Zavarise G, Perillo M. Modeling of the carburizing and nitriding processes. Comput. Mater. Sci.. 2009, 46: 26-35

[6] Cockburn B, Shu C-W. The local discontinuous Galerkin method for time-dependent convection-diffusion systems. SIAM J. Numer. Anal.. 1998, 35: 2440-2463

[7] Cockburn B, Shu C-W. Runge-Kutta discontinuous Galerkin methods for convection-dominated problems. J. Sci. Comput.. 2001, 16: 173-261

[8] Douglas JJr, Dupont T. Alternating-direction Galerkin methods on rectangles, Numerical Solution of Partial Differential Equations-II. 1971, New York, Academic Press: 133214

[9] Duchemin L, Eggers J. The explicit-implicit-null method: removing the numerical instability of PDEs. J. Comput. Phys.. 2014, 263: 37-52

[10] Filbet F, Jin S. A class of asymptotic-preserving schemes for kinetic equations and related problems with stiff sources. J. Comput. Phys.. 2010, 229: 7625-7648

[11] Giraldo FX, Kelly JF, Constantinescu EM. Implicit-explicit formulations of a three-dimensional nonhydrostatic unified model of the atmosphere (NUMA). SIAM J. Sci. Comput.. 2013, 35: B1162-B1194

[12] Shi H, Li Y. Local discontinuous Galerkin methods with implicit-explicit multistep time-marching for solving the nonlinear Cahn-Hilliard equation. J. Comput. Phys.. 2019, 394: 719-731

[13] Smereka P. Semi-implicit level set methods for curvature and surface diffusion motion. J. Sci. Comput.. 2003, 19: 439-456

[14] Tan M, Cheng J, Shu C-W. Stability of high order finite difference and local discontinuous Galerkin schemes with explicit-implicit-null time-marching for high order dissipative and dispersive equations. J. Comput. Phys.. 2022, 464 111314

[15] Tan M, Cheng J, Shu C-W. Stability of spectral collocation schemes with explicit-implicit-null time-marching for convection-diffusion and convection-dispersion equations.. East Asian J. Appl. Math.. 2023, 13: 464-498

[16] Wang H, Zhang Q, Wang S, Shu C-W. Local discontinuous Galerkin methods with explicit-implicit-null time discretizations for solving nonlinear diffusion problems. Sci. China Math.. 2020, 63: 183-204

[17] Xu Y, Shu C-W. Local discontinuous Galerkin methods for three classes of nonlinear wave equations. J. Comput. Math.. 2004, 22: 250-274

[18] Xu Y, Shu C-W. Optimal error estimates of the semidiscrete local discontinuous Galerkin methods for high order wave equations. SIAM J. Numer. Anal.. 2012, 50: 79-104

[19] Yan J, Shu C-W. A local discontinuous Galerkin method for KdV type equations. SIAM J. Numer. Anal.. 2002, 40: 769-791

[20] Yan J, Shu C-W. A local discontinuous Galerkin methods for partial differential equations with higher order derivatives. J. Sci. Comput.. 2002, 17: 27-47

[21] Zhang Q, Wu Z. Numerical simulation for porous medium equation by local discontinuous Galerkin finite element method. J. Sci. Comput.. 2009, 38: 127-148

Funding

NSFC(12031001)

National Key R &D Program of China(2022YFA1004500)

NSF(DMS-2010107)

RIGHTS & PERMISSIONS

Shanghai University