Related papers: Structure-preserving numerical method for Maxwell-…
An asymptotic-preserving (AP) implicit-explicit PN numerical scheme is proposed for the gray model of the radiative transfer equation, where the first- and second-order numerical schemes are discussed for both the linear and nonlinear…
In this paper, we propose energy-conserving Eulerian solvers for the two-species Vlasov-Amp\`{e}re (VA) system and apply the methods to simulate current-driven ion-acoustic instability. The algorithm is generalized from our previous work…
The paper focuses on the development of numerical methods for the compressible Euler equations. It is well-known that if the Mach number is small, the system becomes stiff and hence explicit schemes suffer from severe time-step…
We study the non-canonical symplectic structure, or K-symplectic structure inherited by the charged particle dynamics. Based on the splitting technique, we construct non-canonical symplectic methods which is explicit and stable for the…
We introduce exponential numerical integration methods for stiff stochastic dynamical systems of the form $d\mathbf{z}_t = L(t)\mathbf{z}_tdt + \mathbf{f}(t)dt + Q(t)d\mathbf{W}_t$. We consider the setting of time-varying operators $L(t),…
In this paper, for the first time we propose two linear, decoupled, energy-stable numerical schemes for multi-component two-phase compressible flow with a realistic equation of state (e.g. Peng-Robinson equation of state). The methods are…
The problem of optimal mass transport arises in numerous applications including image registration, mesh generation, reflector design, and astrophysics. One approach to solving this problem is via the Monge-Amp\`ere equation. While recent…
We show quantitative stability results for the geometric "cells" arising in semi-discrete optimal transport problems. Our results show two types of stability, the first is stability of the associated Laguerre cells in measure, without any…
The most rigorous physical description of non-equilibrium gas dynamics is rooted in the numerical solution of the Boltzmann equation. Yet, the large number of degrees of freedom and the wide range of both spatial and temporal scales render…
The maximum entropy principle (MEP) apparently allows us to derive, or justify, fundamental results of equilibrium statistical mechanics. Because of this, a school of thought considers the MEP as a powerful and elegant way to make…
Lorentz invariant structure-preserving algorithms possess reference-independent secular stability, which is vital for simulating relativistic multi-scale dynamical processes. The splitting method has been widely used to construct…
In this paper, exponential energy-preserving methods are formulated and analysed for solving charged-particle dynamics in a strong and constant magnetic field. The resulting method can exactly preserve the energy of the dynamics. Moreover,…
We study the convergence and the stability of fictitious dynamical methods for electrons. First, we show that a particular damped second-order dynamics has a much faster rate of convergence to the ground-state than first-order steepest…
The single particle model (SPM) is a reduced electrochemical model that holds promise for applications in battery management systems due to its ability to accurately capture battery dynamics; however, the numerical discretization of the SPM…
We derive an implicit numerical scheme for the solution of advection equation where the roles of space and time variables are exchanged using the inverse Lax-Wendroff procedure. The scheme contains a linear weight for which it is always…
It is well known that open dynamical systems can admit an uncountable number of (absolutely continuous) conditionally invariant measures (ACCIMs) for each prescribed escape rate. We propose and illustrate a convex optimisation based…
We present a new multi-dimensional, robust, and cell-centered finite-volume scheme for the ideal MHD equations. This scheme relies on relaxation and splitting techniques and can be easily used at high order. A fully conservative version is…
In this paper, we consider numerical approximation of an electrically conductive ferrofluid model, which consists of Navier-Stokes equations, magnetization equation, and magnetic induction equation. To solve this highly coupled, nonlinear,…
This paper concerns the construction and analysis of a numerical scheme for a mixed discrete-continuous fragmentation equation. A finite volume scheme is developed, based on a conservative formulation of a truncated version of the…
Brownian ratchet like stochastic theory for the electrochemical membrane system of Hodgkin-Huxley (HH) is developed. The system is characterized by a continuous variable $Q_m(t)$, representing mobile membrane charge density, and a discrete…