High order spatial discretization for variational time implicit schemes: Wasserstein gradient flows and reaction-diffusion systems
arXiv:2303.08950 · doi:10.1016/j.jcp.2023.112375
Abstract
We design and compute first-order implicit-in-time variational schemes with high-order spatial discretization for initial value gradient flows in generalized optimal transport metric spaces. We first review some examples of gradient flows in generalized optimal transport spaces from the Onsager principle. We then use a one-step time relaxation optimization problem for time-implicit schemes, namely generalized Jordan-Kinderlehrer-Otto schemes. Their minimizing systems satisfy implicit-in-time schemes for initial value gradient flows with first-order time accuracy. We adopt the first-order optimization scheme ALG2 (Augmented Lagrangian method) and high-order finite element methods in spatial discretization to compute the one-step optimization problem. This allows us to derive the implicit-in-time update of initial value gradient flows iteratively. We remark that the iteration in ALG2 has a simple-to-implement point-wise update based on optimal transport and Onsager's activation functions. The proposed method is unconditionally stable for convex cases. Numerical examples are presented to demonstrate the effectiveness of the methods in two-dimensional PDEs, including Wasserstein gradient flows, Fisher--Kolmogorov-Petrovskii-Piskunov equation, and two and four species reversible reaction-diffusion systems.
References in corpus (15)
- Complex Patterns in a Simple System
- Optimal Transport with Proximal Splitting
- The Hele-Shaw asymptotics for mechanical models of tumor growth
- Optimal transport in competition with reaction: the Hellinger-Kantorovich distance and geodesic curves
- A structure-preserving, operator splitting scheme for reaction-diffusion equations with detailed balance
- Unconditionally positivity preserving and energy dissipative schemes for Poisson--Nernst--Planck equations
- Fisher information regularization schemes for Wasserstein gradient flows
- Numerical simulation of nonlinear continuity equations by evolving diffeomorphisms
- 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
- Computational Mean-field information dynamics associated with Reaction diffusion equations
- Generalized Unnormalized Optimal Transport and its fast algorithms
- A dynamic mass transport method for Poisson-Nernst-Planck equations
- Controlling conservation laws II: compressible Navier-Stokes equations