Structure-preserving local discontinuous Galerkin discretization of conformational conversion systems
By: Paola F. Antonietti , Mattia Corti , Sergio Gómez and more
Potential Business Impact:
Makes computer models of nature more accurate.
We investigate a two-state conformational conversion system and introduce a novel structure-preserving numerical scheme that couples a local discontinuous Galerkin space discretization with the backward Euler time-integration method. The model is first reformulated in terms of auxiliary variables involving suitable nonlinear transformations, which allow us to enforce positivity and boundedness at the numerical level. Then, we prove a discrete entropy-stability inequality, which we use to show the existence of discrete solutions, as well as to establish the convergence of the scheme by means of some discrete compactness arguments. As a by-product of the theoretical analysis, we also prove the existence of global weak solutions satisfying the system's physical bounds. Numerical results validate the theoretical results and assess the capabilities of the proposed method in practice.
Similar Papers
Structure-Preserving High-Order Methods for the Compressible Euler Equations in Potential Temperature Formulation for Atmospheric Flows
Numerical Analysis
Makes weather forecasts more accurate by better predicting air movement.
Sequential symmetric interior penalty discontinuous Galerkin method for fully coupled quasi-static thermo-poroelasticity problems
Numerical Analysis
Solves hard math problems faster for science.
Structure-preserving nodal DG method for the Euler equations with gravity: well-balanced, entropy stable, and positivity preserving
Numerical Analysis
Simulates gravity and balance in fluid flow.