Entropy dissipative higher order accurate positivity preserving time-implicit discretizations for nonlinear degenerate parabolic equations

DOI10.1016/J.CAM.2023.115674MaRDI QIDQ6145196zbMATH OpenFDO

Authors Fengna Yan, J. J. W. van der Vegt, Yinhua Xia, Yan Xu

Publication date 30 January 2024

Published in Journal of Computational and Applied Mathematics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/2301.01427

zbMATH Keywords

DIRK methodsentropy dissipationKKT limiterlocal discontinuous Galerkin discretizationsnonlinear degenerate parabolic equations

Mathematics Subject Classification ID

Positive solutions to PDEs (35B09) Nonlinear initial, boundary and initial-boundary value problems for nonlinear parabolic equations (35K61) Degenerate parabolic equations (35K65) Numerical mathematical programming methods (65K05) Multistep, Runge-Kutta and extrapolation methods for ordinary differential equations (65L06) Finite difference methods for initial value and initial-boundary value problems involving PDEs (65M06) Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs (65M60) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30) Nonlinear programming (90C30)

Abstract: We develop entropy dissipative higher order accurate local discontinuous Galerkin (LDG) discretizations coupled with Diagonally Implicit Runge-Kutta (DIRK) methods for nonlinear degenerate parabolic equations with a gradient flow structure. Using the simple alternating numerical flux, we construct DIRK-LDG discretizations that combine the advantages of higher order accuracy, entropy dissipation and proper long-time behavior. The implicit time-discrete methods greatly alleviate the time-step restrictions needed for the stability of the numerical discretizations. Also, the larger time step significantly improves computational efficiency. We theoretically prove the unconditional entropy dissipation of the implicit Euler-LDG discretization. Next, in order to ensure the positivity of the numerical solution, we use the Karush-Kuhn-Tucker (KKT) limiter, which couples the positivity inequality constraint with higher order accurate DIRK-LDG discretizations using Lagrange multipliers. In addition, mass conservation of the positivity-limited solution is ensured by imposing a mass conservation equality constraint to the KKT equations. The unique solvability and unconditional entropy dissipation for an implicit first order accurate in time, but higher order accurate in space, KKT-LDG discretizations are proved, which provides a first theoretical analysis of the KKT limiter. Finally, numerical results demonstrate the higher order accuracy and entropy dissipation of the KKT-DIRK-LDG discretizations for problems requiring a positivity limiter.

Recommendations

• Positivity preserving limiters for time-implicit higher order accurate discontinuous Galerkin discretizations
• A discontinuous Galerkin method for nonlinear parabolic equations and gradient flow problems with interaction potentials
• An entropy satisfying discontinuous Galerkin method for nonlinear Fokker-Planck equations
• Entropy-dissipating semi-discrete Runge-Kutta schemes for nonlinear diffusion equations
• Positivity preserving high-order local discontinuous Galerkin method for parabolic equations with blow-up solutions

Cites work

• A finite volume scheme for nonlinear degenerate parabolic equations
• A finite-volume method for nonlinear nonlocal equations with a gradient flow structure
• A high order adaptive time-stepping strategy and local discontinuous Galerkin method for the modified phase field crystal equation
• A mixed finite element method for nonlinear diffusion equations
• A new Lagrange multiplier approach for constructing structure preserving schemes. I: Positivity preserving
• An \(h\)-adaptive local discontinuous Galerkin method for the Navier-Stokes-Korteweg equations
• An entropy satisfying discontinuous Galerkin method for nonlinear Fokker-Planck equations
• Bound/positivity preserving and energy stable scalar auxiliary variable schemes for dissipative systems: applications to Keller-Segel and Poisson-Nernst-Planck equations
• Diagonally Implicit Runge–Kutta Methods for Stiff O.D.E.’s
• Discontinuous Galerkin method for Krause's consensus models and pressureless Euler equations
• Fermi-Dirac-Fokker-Planck equation: well-posedness and long-time asymptotics
• Finite time blow up in Kaniadakis-Quarati model of Bose-Einstein particles
• Finite-Dimensional Variational Inequalities and Complementarity Problems
• Implicit positivity-preserving high-order discontinuous Galerkin methods for conservation laws
• Maximum-principle-satisfying third order discontinuous Galerkin schemes for Fokker-Planck equations
• Numerical simulation for porous medium equation by local discontinuous Galerkin finite element method
• On maximum-principle-satisfying high order schemes for scalar conservation laws
• On positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations on rectangular meshes
• On the minimization problem of sub-linear convex functionals
• Positivity preserving limiters for time-implicit higher order accurate discontinuous Galerkin discretizations
• Solving ordinary differential equations. II: Stiff and differential-algebraic problems.
• Stability analysis and error estimates of semi-implicit spectral deferred correction coupled with local discontinuous Galerkin method for linear convection-diffusion equations
• Superconvergence of the local discontinuous Galerkin method for elliptic problems on Cartesian grids
• The entropy satisfying discontinuous Galerkin method for Fokker-Planck equations
• The Local Discontinuous Galerkin Method for Time-Dependent Convection-Diffusion Systems
• The porous medium equation. Mathematical theory

Cited in

(4)

• An energy stable local discontinuous Galerkin method for the isothermal Navier-Stokes-Korteweg equations
• An energy stable local discontinuous Galerkin method for a binary compressible flow
• Analysis of average bound preserving time-implicit discretizations for convection-diffusion-reaction equation
• A second-order bound-preserving exponential scheme for degenerate parabolic equations

This page was built for publication: Entropy dissipative higher order accurate positivity preserving time-implicit discretizations for nonlinear degenerate parabolic equations

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6145196)