Related papers: Efficient, positive, and energy stable schemes for…
In this work we design and analyze a free energy satisfying finite difference method for solving Poisson-Nernst-Planck equations in a bounded domain. The algorithm is of second order in space, with numerical solutions satisfying all three…
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…
An asymptotic preserving and energy stable scheme for the Euler-Poisson system under the quasineutral scaling is designed and analysed. Correction terms are introduced in the convective fluxes and the electrostatic potential, which lead to…
In this paper, we present a first-order finite element scheme for the viscoelastic electrohydrodynamic model. The model incorporates the Poisson-Nernst-Planck equations to describe the transport of ions and the Oldroyd-B constitutive model…
In this paper we consider unconditionally energy stable numerical schemes for the nonstationary 3D magneto-micropolar equations that describes the microstructure of rigid microelements in electrically conducting fluid flow under some…
Efficient and energy stable high order time marching schemes are very important but not easy to construct for the study of nonlinear phase dynamics. In this paper, we propose and study two linearly stabilized second order semi-implicit…
We propose fully discrete, implicit-in-time finite-volume schemes for a general family of non-linear and non-local Fokker-Planck equations with a gradient-flow structure, usually known as aggregation-diffusion equations, in any dimension.…
We introduce a new structure preserving, second order in time relaxation-type scheme for approximating solutions of the Schr\"odinger-Poisson system. More specifically, we use the Crank-Nicolson scheme as a time stepping mechanism, whilst…
In this paper, we present a class of nonuniform time-stepping, high-order linear stabilized schemes that can preserve both the discrete energy stability and maximum-bound principle (MBP) for the time-fractional Allen-Cahn equation. To this…
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 introduce second order and fourth order space discretization via finite difference implementation of the finite element method for solving Fokker-Planck equations associated with irreversible processes. The proposed…
We propose and analyze a novel approach to construct structure preserving approximations for the Poisson-Nernst-Planck equations, focusing on the positivity preserving and mass conservation properties. The strategy consists of a standard…
We generalize the explicit high-order positivity-preserving entropy stable spectral collocation schemes developed in Upperman 2023 and Yamaleev 2023 for the three-dimensional (3D) compressible Navier Stokes equations to a time implicit…
We develop new numerical schemes for Vlasov--Poisson equations with high-order accuracy. Our methods are based on a spatially monotonicity-preserving (MP) scheme and are modified suitably so that positivity of the distribution function is…
This paper introduces a robust reformulation of the incompressible Navier-Stokes equations, establishing a foundational framework for designing efficient, structure-preserving algorithms that strictly conserve the original energy…
We study a two-grid strategy for decoupling the time-dependent Poisson-Nernst-Planck equations describing the mass concentration of ions and the electrostatic potential. The computational system is decoupled to smaller systems by using…
We propose and analyze a linearly stabilized semi-implicit diffusive Crank--Nicolson scheme for the Cahn--Hilliard gradient flow. In this scheme, the nonlinear bulk force is treated explicitly with two second-order stabilization terms. This…
With a small parameter $\epsilon$, Poisson-Nernst-Planck (PNP) systems over a finite one-dimensional (1D) spatial domain have steady state solutions, called 1D boundary layer solutions, which profiles form boundary layers near boundary…
We propose and analyse numerical schemes for a system of quasilinear, degenerate evolution equations modelling biofilm growth as well as other processes such as flow through porous media and the spreading of wildfires. The first equation in…
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…