Related papers: A Positive and Energy Stable Numerical Scheme for …
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,…
We present several first-order and second-order numerical schemes for the Cahn-Hilliard equation with discrete unconditional energy stability. These schemes stem from the generalized Positive Auxiliary Variable (gPAV) idea, and require only…
In this paper, we design and analyze second order positive and free energy satisfying schemes for solving diffusion equations with interaction potentials. The semi-discrete scheme is shown to conserve mass, preserve solution positivity, and…
We study numerical methods for porous media equation (PME). There are two important characteristics: the finite speed propagation of the free boundary and the potential waiting time, which make the problem not easy to handle. Based on…
We design an arbitrary-order free energy satisfying discontinuous Galerkin (DG) method for solving time-dependent Poisson-Nernst-Planck systems. Both the semi-discrete and fully discrete DG methods are shown to satisfy the corresponding…
This work presents new parallelizable numerical schemes for the integration of Dissipative Particle Dynamics with Energy conservation (DPDE). So far, no numerical scheme introduced in the literature is able to correctly preserve the energy…
We present new linear energy-stable numerical schemes for numerical simulation of complex polymer-solvent mixtures. The mathematical model proposed by Zhou, Zhang and E (Physical Review E 73, 2006) consists of the Cahn-Hilliard equation…
In this work, the MMC-TDGL equation, a stochastic Cahn-Hilliard equation is solved numerically by using the finite difference method in combination with a convex splitting technique of the energy functional. For the non-stochastic case, we…
An energy-conserving and an energy-and-enstrophy conserving numerical schemes are derived, by approximating the Hamiltonian formulation, based on the Poisson brackets and the vorticity-divergence variables, of the inviscid shallow water…
We develop a structure-preserving numerical discretization for the electrostatic Euler-Poisson equations with a constant magnetic field. The scheme preserves positivity of the density, positivity of the internal energy and a minimum…
We propose an edge averaged finite element(EAFE) discretization to solve the Heat-PNP (Poisson-Nernst-Planck) equations approximately. Our method enforces positivity of the computed charged density functions and temperature function. Also…
We present a linear, second order fully discrete numerical scheme on a staggered grid for a thermodynamically consistent hydrodynamic phase field model of binary compressible fluid flow mixtures derived from the generalized Onsager…
In this paper, we investigate numerically a diffuse interface model for the Navier-Stokes equation with fluid-fluid interface when the fluids have different densities \cite{Lowengrub1998}. Under minor reformulation of the system, we show…
Varieties of energy-stable numerical methods have been developed for incompressible two-phase flows based on the Navier-Stokes-Cahn-Hilliard (NSCH) model in the Eulerian framework, while few investigations have been made in the Lagrangian…
Small-scale plasticity problems are often characterised by different patterning behaviours ranging from macroscopic down to the atomistic scale. In successful models of such complex behaviour, its origin lies within non-convexity of the…
In this paper, we study a second-order accurate and linear numerical scheme for the nonlocal Cahn-Hilliard equation. The scheme is established by combining a modified Crank-Nicolson approximation and the Adams-Bashforth extrapolation for…
We propose a new fully-discretized finite difference scheme for a quantum diffusion equation, in both one and two dimensions. This is the first fully-discretized scheme with proven positivity-preserving and energy stable properties using…
Energy conserving particle-in-cell schemes are constructed for a class of reduced relativistic Vlasov--Maxwell equations of laser-plasma interaction. Discrete Poisson equation is also satisfied by the numerical solution. Specifically,…
In this article, we present and analyze a finite element numerical scheme for a three-component macromolecular microsphere composite (MMC) hydrogel model, which takes the form of a ternary Cahn-Hilliard-type equation with…
We study the numerical solution of a Cahn-Hilliard/Allen-Cahn system with strong coupling through state and gradient dependent non-diagonal mobility matrices. A fully discrete approximation scheme in space and time is proposed which…