Related papers: Provably Positivity-Preserving Constrained Transpo…
Ion transport, often described by the Poisson--Nernst--Planck (PNP) equations, is ubiquitous in electrochemical devices and many biological processes of significance. In this work, we develop conservative, positivity-preserving, energy…
This paper aims at developing exactly energy-conservative and structure-preserving finite volume schemes for the discretisation of first-order symmetric-hyperbolic and thermodynamically compatible (SHTC) systems of partial differential…
We construct a structure-preserving finite element method and time-stepping scheme for compressible barotropic magnetohydrodynamics (MHD) both in the ideal and resistive cases, and in the presence of viscosity. The method is deduced from…
We present a novel structure-preserving semi-implicit finite volume method on vertex-based staggered meshes for the compatible discretization of first order systems of time-dependent partial differential equations (PDEs). The method…
A third order shock-capturing numerical scheme for three-dimensional special relativistic magnetohydrodynamics (3-D RMHD) is presented and validated against several numerical tests. The simple and efficient central scheme described in Paper…
Maintaining conservation laws in the fully discrete setting is critical for accurate long-time behavior of numerical simulations and requires accounting for discrete conservation properties in both space and time. This paper derives…
First-order energy dissipative schemes in time are available in literature for the Poisson-Nernst-Planck (PNP) equations, but second-order ones are still in lack. This work proposes novel second-order discretization in time and finite…
We present an adaptive-order positivity-preserving conservative finite-difference scheme that allows a high-order solution away from shocks and discontinuities while guaranteeing positivity and robustness at discontinuities. This is…
We present a new numerical code, ECHO, based on an Eulerian Conservative High Order scheme for time dependent three-dimensional general relativistic magnetohydrodynamics (GRMHD) and magnetodynamics (GRMD). ECHO is aimed at providing a…
We present the implementation of a three-dimensional, second order accurate Godunov-type algorithm for magneto-hydrodynamic (MHD), in the adaptive-mesh-refinement (AMR) cosmological code {\tt CHARM}. The algorithm is based on the full…
A positivity-preserving conservative semi-Lagrangian transport model by multi-moment finite volume method has been developed on the cubed-sphere grid. In this paper, two kinds of moments, i.e. point values (PV moment) at cell boundaries and…
Numerical simulation of the spherically symmetric Einstein--Euler (EE) system faces severe challenges due to the stringent physical admissibility constraints of relativistic fluids and the geometric singularities inherent in metric…
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…
A spatial second-order scheme for the nonlinear radiative transfer equations is introduced in this paper. The discretization scheme is based on the filtered spherical harmonics ($FP_N$) method for the angular variable and the unified gas…
In this paper, we propose and analyze a temporally second-order accurate, fully discrete finite element method for the magnetohydrodynamic (MHD) equations. A modified Crank--Nicolson method is used to discretize the model and appropriate…
We present compatible finite element space discretizations for the ideal compressible magnetohydrodynamic equations. The magnetic field is considered both in div- and curl-conforming spaces, leading to a strongly or weakly preserved…
In this paper, we propose and analyze a second order accurate (in both time and space) numerical scheme for the Poisson-Nernst-Planck-Navier-Stokes system, which describes the ion electro-diffusion in fluids. In particular, the…
In this paper, we develop parametrized positivity satisfying flux limiters for the high order finite difference Runge-Kutta weighted essentially non-oscillatory (WENO) scheme solving compressible Euler equations to maintain positive density…
In this manuscript we present a novel and efficient numerical method for the compressible viscous and resistive MHD equations for all Mach number regimes. The time-integration strategy is a semi-implicit splitting, combined with a hybrid…
In this paper, we construct an explicit, second-order, and maximum-principle-preserving Crouzeix-Raviart (CR) finite element method for two-dimensional time-dependent transport equation. The key observation is that the mass matrix of the CR…