A discontinuous Galerkin method for nonlinear parabolic equations and gradient flow problems with interaction potentials

Abstract

We consider a class of time dependent second order partial differential equations governed by a decaying entropy. The solution usually corresponds to a density distribution, hence positivity (non-negativity) is expected. This class of problems covers important cases such as Fokker-Planck type equations and aggregation models, which have been studied intensively in the past decades. In this paper, we design a high order discontinuous Galerkin method for such problems. If the interaction potential is not involved, or the interaction is defined by a smooth kernel, our semi-discrete scheme admits an entropy inequality on the discrete level. Furthermore, by applying the positivity-preserving limiter, our fully discretized scheme produces non-negative solutions for all cases under a time step constraint. Our method also applies to two dimensional problems on Cartesian meshes. Numerical examples are given to confirm the high order accuracy for smooth test cases and to demonstrate the effectiveness for preserving long time asymptotics.

Cited by in corpus (22)

• Solving Fokker-Planck equation using deep learning
• Fisher information regularization schemes for Wasserstein gradient flows
• Positivity-preserving third order DG schemes for Poisson--Nernst--Planck equations
• Primal dual methods for Wasserstein gradient flows
• Is the Classic Convex Decomposition Optimal for Bound-Preserving Schemes in Multiple Dimensions?
• Functional flows for complex effective actions
• Bound-Preserving Finite-Volume Schemes for Systems of Continuity Equations with Saturation
• Positive and free energy satisfying schemes for diffusion with interaction potentials
• A Finite Volume Method for Continuum Limit Equations of Nonlocally Interacting Active Chiral Particles
• Well-balanced finite volume schemes for hydrodynamic equations with general free energy
• A scalable DG solver for the electroneutral Nernst-Planck equations
• On structure-preserving discontinuous Galerkin methods for Hamiltonian partial differential equations: Energy conservation and multi-symplecticity
• Proximal Galerkin: A structure-preserving finite element method for pointwise bound constraints
• Enforcing strong stability of explicit Runge--Kutta methods with superviscosity
• THIRD ORDER MAXIMUM-PRINCIPLE-SATISFYING DG SCHEMES Third Order Maximum-Principle-Satisfying DG schemes for Convection-Diffusion problems with Anisotropic Diffusivity DIFFUSIVITY
• Positivity-preserving and energy-dissipative finite difference schemes for the Fokker-Planck and Keller-Segel equations
• An entropy stable high-order discontinuous Galerkin method for cross-diffusion gradient flow systems
• Structure-Preserving Numerical Methods for Nonlinear Fokker--Planck Equations with Nonlocal Interactions by an Energetic Variational Approach
• Strong stability of explicit Runge-Kutta time discretizations
• A structure-preserving discontinuous Galerkin scheme for the Fischer-KPP equation
• A structure-preserving LDG discretization of the Fisher-Kolmogorov equation for modeling neurodegenerative diseases
• Drift-diffusion equations with saturation