Related papers: Structure-preserving numerical method for Maxwell-…
Modeling of kinetic plasmas using electromagnetic particle in cell methods (EM-PIC) is a problem that is well worn, in that methods developed have been used extensively both understanding physics and exploiting them for device design.…
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…
The Poisson-Nernst-Planck (PNP) equations are fundamental for modeling ion transport in electrochemical systems, capturing the intricate interplay of concentration gradients, electric fields, and ion fluxes essential for applications such…
The molecular mechanism of ion channel gating and substrate modulation is elusive for many voltage gated ion channels, such as eukaryotic sodium ones. The understanding of channel functions is a pressing issue in molecular biophysics and…
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…
As an effective emulator of ill-conditioned power flow, continuous Newton methods (CNMs) have been extensively investigated using explicit and implicit numerical integration algorithms. Explicit CNMs are prone to non-convergence issues due…
Understanding the system-level spectral immittance response of capacitive energy storage devices with analytically tractable physics-based models is not only important for the progress of the technology, but also allows to develop new…
A well known procedure to get quantum hydrodynamical models for charge transport is to resort to the Wigner equations and deduce the hierarchy of the moment equations as in the semiclassical approach. If one truncates the moment hierarchy…
In this paper we consider the development of numerical schemes for mean-field equations describing the collective behavior of a large group of interacting agents. The schemes are based on a generalization of the classical Chang-Cooper…
The PAMPA (Point-Average-Moment PolynomiAl-interpreted) method, proposed in [R. Abgrall, Commun. Appl. Math. Comput., 5: 370-402, 2023], combines conservative and non-conservative formulations of hyperbolic conservation laws to evolve cell…
This paper presents a novel structure-preserving scheme for Euler equations, focusing on the numerical conservation of entropy and kinetic energy. Explicit flux functions engineered to conserve entropy are introduced within the…
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 introduce and analyse numerical schemes for the homogeneous and the kinetic L\'evy-Fokker-Planck equation. The discretizations are designed to preserve the main features of the continuous model such as conservation of…
This paper presents a highly robust third-order accurate finite volume weighted essentially non-oscillatory (WENO) method for special relativistic hydrodynamics on unstructured triangular meshes. We rigorously prove that the proposed method…
Higher order finite volume schemes for magnetohydrodynamics (MHD) and relativistic magnetohydrodynamics (RMHD) are very valuable because they allow us to carry out astrophysical simulations with very high accuracy. However, astrophysical…
We present a novel hydrostatic and non-hydrostatic equilibria preserving Point-Average-Moment PolynomiAl-interpreted (PAMPA) method for solving the one-dimensional hyperbolic balance laws, with applications to the shallow water models…
Maxwell's equations describe the evolution of electromagnetic fields, together with constraints on the divergence of the magnetic and electric flux densities. These constraints correspond to fundamental physical laws: the nonexistence of…
We develop a structure-preserving computational framework for optimal mixing control in incompressible flows. Our approach exactly conserves the continuous system's key invariants (mass and $L^2$-energy), while also maintaining discrete…
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…
Accurate simulations of a condensed system of ions or polar molecules are concerned with proper handlings of the involved electrostatics. For such a Coulomb system at a charged planar interface, the Coulomb interaction averaged over the…