Related papers: The Onsager principle and structure preserving num…
We propose a novel Skew Gradient Embedding (SGE) framework for systematically reformulating thermodynamically consistent partial differential equation (PDE) models-capturing both reversible and irreversible processes-as generalized gradient…
A Lagrangian numerical scheme for solving nonlinear degenerate Fokker-Planck equations in space dimensions $d\ge2$ is presented. It applies to a large class of nonlinear diffusion equations, whose dynamics are driven by internal energies…
For gradient flows, the existing structure-preserving schemes are difficult to achieve arbitrary high-order accuracy in time while preserving maximum-principle (MBP) and energy dissipating simultaneously. In this paper, we develop a new…
Onsager's variational principle (OVP) provides us with a systematic way to derive dynamical equations for various soft matter and active matter. By reformulating the Onsager-Machlup variational principle (OMVP), which is a time-global…
We present an adaptation of the so-called structural method \cite{CMM23} for Hamiltonian systems, and redesign the method for this specific context, which involves two coupled differential systems. Structural schemes decompose the problem…
The numerical integration of phase-field equations is a delicate task which needs to recover at the discrete level intrinsic properties of the solution such as energy dissipation and maximum principle. Although the theory of energy…
We evaluate the Onsager matrix for a system under time-periodic driving by considering all its Fourier components. By application of the second law, we prove that all the fluxes converge to zero in the limit of zero dissipation. Reversible…
The Cahn-Hilliard equation and extensions, notably the Cahn-Hilliard-Darcy and Cahn-Hilliard-Navier-Stokes systems, provide widely used frameworks for coupling interfacial thermodynamics with flow. This review surveys the thermodynamic…
The Navier--Stokes equations are commonly used to model and to simulate flow phenomena. We introduce the basic equations and discuss the standard methods for the spatial and temporal discretization. We analyse the semi-discrete equations --…
We study the systematic numerical approximation of a class of Allen-Cahn type problems modeling the motion of phase interfaces. The common feature of these models is an underlying gradient flow structure which gives rise to a decay of an…
In this paper, we construct a semi-implicit finite difference method for the time dependent Poisson-Nernst-Planck system. Although the Poisson-Nernst-Planck system is a nonlinear system, the numerical method presented in this paper only…
We propose a machine learning framework to accelerate numerical computations of time-dependent ODEs and PDEs. Our method is based on recasting (generalizations of) existing numerical methods as artificial neural networks, with a set of…
The aim of this paper is the derivation of structure preserving schemes for the solution of the EPDiff equation, with particular emphasis on the two dimensional case. We develop three different schemes based on the Discrete Variational…
Numerical schemes for the solution of the Euler equations have recently been developed, which involve the discretisation of the internal energy equation, with corrective terms to ensure the correct capture of shocks, and, more generally,…
How to develop efficient numerical schemes while preserving the energy stability at the discrete level is a challenging issue for the three component Cahn-Hilliard phase-field model. In this paper, we develop first and second order temporal…
A simple variational Lagrangian is proposed for the time development of an arbitrary density matrix, employing the "factorization" of the density. Only the "kinetic energy" appears in the Lagrangian. The formalism applies to pure and mixed…
In this paper, we design a novel class of arbitrarily high-order structure-preserving numerical schemes for the time-dependent Gross-Pitaevskii equation with angular momentum rotation in three dimensions. Based on the idea of the scalar…
There is a qualitative difference between one-dimensional and multi-dimensional solutions to the Euler equations: new features that arise are vorticity and a nontrivial incompressible (low Mach number) limit. They present challenges to…
This paper detailedly discusses the locally one-dimensional numerical methods for efficiently solving the three-dimensional fractional partial differential equations, including fractional advection diffusion equation and Riesz fractional…
Employing a phase space which includes the (Riemann-Liouville) fractional derivative of curves evolving on real space, we develop a restricted variational principle for Lagrangian systems yielding the so-called restricted fractional…