A structure-preserving, operator splitting scheme for reaction-diffusion equations with detailed balance

Abstract

In this paper, we propose and analyze a positivity-preserving, energy stable numerical scheme for certain type reaction-diffusion systems involving the Law of Mass Action with the detailed balance condition. The numerical scheme is constructed based on a recently developed energetic variational formulation, in which the reaction part is reformulated in terms of reaction trajectories. The fact that both the reaction and the diffusion parts dissipate the same free energy opens a path of an energy stable, operator splitting scheme for these systems. At the reaction stage, we solve equations of reaction trajectories by treating all the logarithmic terms in the reformulated form implicitly due to their convex nature. The positivity-preserving property and unique solvability can be theoretically proved, based on the singular behavior of the logarithmic function around the limiting value. Moreover, the energy stability of this scheme at the reaction stage can be proved by a careful convexity analysis. Similar techniques are used to establish the positivity-preserving property and energy stability for the standard semi-implicit solver at the diffusion stage. As a result, a combination of these two stages leads to a positivity-preserving and energy-stable numerical scheme for the original reaction-diffusion system. To our best knowledge, it is the first time to report an energy-dissipation-law-based operator splitting scheme to a nonlinear PDE with variational structures. Several numerical examples are presented to demonstrate the robustness of the proposed operator splitting scheme.

References in corpus (20)

• Complex Patterns in a Simple System
• The Hele-Shaw asymptotics for mechanical models of tumor growth
• Energy Stable and Efficient Finite-Difference Nonlinear Multigrid Schemes for the Modified Phase Field Crystal Equation
• An Energetic Variational Approach for the Cahn--Hilliard Equation with Dynamic Boundary Condition: Model Derivation and Mathematical Analysis
• A positivity-preserving, energy stable scheme for a Ternary Cahn-Hilliard system with the singular interfacial parameters
• Phase-field dynamics with transfer of materials: The Cahn--Hilliard equation with reaction rate dependent dynamic boundary conditions
• A positivity-preserving second-order BDF scheme for the Cahn-Hilliard equation with variable interfacial parameters
• Unconditionally positivity preserving and energy dissipative schemes for Poisson--Nernst--Planck equations
• Convergence analysis of Strang splitting for Vlasov-type equations
• Mathematical Formalism of Nonequilibrium Thermodynamics for Nonlinear Chemical Reaction Systems with General Rate Law
• Field Theory of Reaction-Diffusion: Mass Action with an Energetic Variational Approach
• Thermodynamics in the Limit of Irreversible Reactions
• On Lagrangian schemes for porous medium type generalized diffusion equations: a discrete energetic variational approach
• A Lagrangian scheme for the solution of nonlinear diffusion equations using moving simplex meshes
• An accurate front capturing scheme for tumor growth models with a free boundary limit
• Numerical methods for Porous Medium Equation by an Energetic Variational Approach
• Decay to equilibrium for energy-reaction-diffusion systems
• An almost symmetric Strang splitting scheme for nonlinear evolution equations
• Operator splitting for nonautonomous delay equations
• A two species micro-macro model of wormlike micellar solutions and its maximum entropy closure approximations: An energetic variational approach

Cited by in corpus (9)

• A positivity-preserving, energy stable scheme for a Ternary Cahn-Hilliard system with the singular interfacial parameters
• High order spatial discretization for variational time implicit schemes: Wasserstein gradient flows and reaction-diffusion systems
• A Uniquely Solvable, Positivity-Preserving and Unconditionally Energy Stable Numerical Scheme for the Functionalized Cahn-Hilliard Equation with Logarithmic Potential
• A two species micro-macro model of wormlike micellar solutions and its maximum entropy closure approximations: An energetic variational approach
• An energy stable finite element scheme for the three-component Cahn-Hilliard-type model for macromolecular microsphere composite hydrogels
• On a Reversible Gray-Scott Type System from Energetic Variational Approach and Its Irreversible Limit
• A Thermodynamically-Consistent Phase Field Crystal Model of Solidification with Heat Flux
• A Structure-Preserving Rational Integrator for the Replicator Dynamics on the Probability Simplex
• A brief overview of existence results and decay time estimates for a mathematical modeling of scintillating crystals