Related papers: Efficient, positive, and energy stable schemes for…
We study a system of Maxwell's equations that describes the time evolution of electromagnetic fields with an additional electric scalar variable to make the system amenable to a mixed finite element spatial discretization. We demonstrate…
In this paper, we extend the positivity-preserving, entropy stable first-order finite volume-type scheme developed for the one-dimensional compressible Navier-Stokes equations in [1] to three spatial dimensions. The new first-order scheme…
In this paper, we develop a novel staggered mesh (SM) approach for general nonlinear dissipative systems with arbitrary energy distributions (including cases with known or unknown energy lower bounds). Based on this framework, we propose…
We design and analyse an energy stable, structure preserving and well-balanced scheme for the Ripa system of shallow water equations. The energy stability of the numerical solutions is achieved by introducing appropriate stabilisation terms…
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…
In this paper, a linear second order numerical scheme is developed and investigated for the Allen-Cahn equation with a general positive mobility. In particular, our fully discrete scheme is mainly constructed based on the Crank-Nicolson…
We present a set of linear, second order, unconditionally energy stable schemes for the Allen-Cahn equation with nonlocal constraints that preserves the total volume of each phase in a binary material system. The energy quadratization…
We have developed efficient numerical algorithms for solving 3D steady-state Poisson-Nernst-Planck (PNP) equations with excess chemical potentials described by the classical density functional theory (cDFT). The coupled PNP equations are…
In this paper, we consider numerical approximations for the viscous Cahn-Hilliard equation with hyperbolic relaxation. This type of equations processes energy-dissipative structure. The main challenge in solving such a diffusive system…
In this paper, we investigate the correlated diffusion of two ion species governed by a Poisson-Nernst-Planck (PNP) system. Here we further validate the numerical scheme recently proposed in \cite{astuto2025asymptotic}, where a time…
We present and analyze two stabilized finite element methods for solving numerically the Poisson--Nernst--Planck equations. The stabilization we consider is carried out by using a shock detector and a discrete graph Laplacian operator for…
Partial differential equations (PDEs) describing thermodynamically isolated systems typically possess conserved quantities (like mass, momentum, and energy) and dissipated quantities (like entropy). Preserving these conservation and…
In earlier work [H. Liu and Z. Wang, J. Comput. Phys., 328(2017)], an arbitrary high-order conservative and energy-dissipative direct discontinuous Galerkin (DDG) scheme was developed. Although this scheme enforced solution positivity using…
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…
The Poisson-Nernst-Planck equations with generalized Frumkin-Butler-Volmer boundary conditions (PNP-FBV) describe ion transport with Faradaic reactions and have applications in a wide variety of fields. Using an adaptive time-stepper based…
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,…
The analysis of structure-preserving numerical methods for the Poisson--Nernst--Planck (PNP) system has attracted growing interests in recent years. In this work, we provide an optimal rate convergence analysis and error estimate for finite…
We discuss structure-preserving numerical discretizations for repulsive and attractive Euler-Poisson equations that find applications in fluid-plasma and self-gravitation modeling. The scheme is fully discrete and structure preserving in…
Poisson-Nernst-Planck equations are widely used to describe the electrodiffusion of ions in a solvated biomolecular system. Two kinds of two-grid finite element algorithms are proposed to decouple the steady-state Poisson-Nernst-Planck…
The Poisson--Nernst--Planck (PNP) equations have been widely applied to describe ionic transport in ion channels, nanofluidic devices, and many electrochemical systems. Despite their wide applications, the PNP equations fail in predicting…