Related papers: A structure-preserving implicit exponential time d…
A dynamic mass-transport method is proposed for approximately solving the Poisson-Nernst-Planck(PNP) equations. The semi-discrete scheme based on the JKO type variational formulation naturally enforces solution positivity and the energy law…
In this paper, we consider a nonlinear and nonlocal parabolic model for multi-species ionic fluids and introduce a semi-implicit finite volume scheme, which is second order accurate in space, first order in time and satisfies the following…
It is well-known that the Allen-Cahn equation not only satisfies the energy dissipation law but also possesses the maximum bound principle (MBP) in the sense that the absolute value of its solution is pointwise bounded for all time by some…
This paper aims to develop a linearly implicit structure-preserving numerical scheme for the space fractional sine-Gordon equation, which is based on the newly developed invariant energy quadratization method. First, we reformulate the…
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 our latest studies, by introducing the novel order-preserving (OP) criterion, we have successfully addressed the widely concerned issue of the previously published mapped weighted essentially non-oscillatory (WENO) schemes that it is…
This paper proposes a charge-conserving, variational, spatio-temporal discretization for the drift-kinetic Vlasov-Maxwell system, utilizing finite-elements for the electromagnetic fields and the particle-in-cell approach for the Vlasov…
We present a divergence-free semi-implicit finite volume scheme for the simulation of the ideal magnetohydrodynamics (MHD) equations which is stable for large time steps controlled by the local transport speed at all Mach and Alfv\'en…
In this work, we aim at constructing numerical schemes, that are as efficient as possible in terms of cost and conservation of invariants, for the Vlasov--Fokker--Planck system coupled with Poisson or Amp\`ere equation. Splitting methods…
We present a novel class of high-order space-time finite element schemes for the Poisson-Nernst-Planck (PNP) equations. We prove that our schemes are mass conservative, positivity preserving, and unconditionally energy stable for any order…
We develop a general framework for designing conservative numerical methods based on summation by parts operators and split forms in space, combined with relaxation Runge-Kutta methods in time. We apply this framework to create new classes…
Conventional explicit electromagnetic particle-in-cell (PIC) algorithms do not conserve discrete energy exactly. Time-centered fully implicit PIC algorithms can conserve discrete energy exactly, but may introduce large dispersion errors in…
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 consider the development of exponential methods for the robust time discretization of space inhomogeneous Boltzmann equations in stiff regimes. Compared to the space homogeneous case, or more in general to the case of splitting based…
We present a conservative/dissipative time integration scheme for nonlinear mechanical systems. Starting from a weak form, we derive algorithmic forces and velocities that guarantee the desired conservation/dissipation properties. Our…
The main objective of this paper is to present an efficient structure-preserving scheme, which is based on the idea of the scalar auxiliary variable approach, for solving the space fractional nonlinear Schr\"{o}dinger equation. First, we…
In this paper, we consider numerical approximation of constrained gradient flows of planar closed curves, including the Willmore and the Helfrich flows. These equations have energy dissipation and the latter has conservation properties due…
The parabolic-elliptic Keller-Segel equation with sensitivity saturation, because of its pattern formation ability, is a challenge for numerical simulations. We provide two finite-volume schemes whose goals are to preserve, at the discrete…
In this paper, we propose two efficient fully-discrete schemes for Q-tensor flow of liquid crystals by using the first- and second-order stabilized exponential scalar auxiliary variable (sESAV) approach in time and the finite difference…
This paper proposes a new class of mass or energy conservative numerical schemes for the generalized Benjamin-Ono (BO) equation on the whole real line with arbitrarily high-order accuracy in time. The spatial discretization is achieved by…