Related papers: Explicit volume-preserving numerical schemes for r…
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…
Many systems of interest in general relativistic astrophysics, including neutron stars, accreting compact objects in X-ray binaries and active galactic nuclei, core collapse, and collapsars, are assumed to be approximately spherically…
We present a fully relativistic numerical method for the study of cosmological problems using the Baumgarte-Shapiro-Shibata-Nakamura formalism on a dynamical Friedmann-Lema\^itre-Robertson-Walker background. This has many potential…
In recent years it has been recognized that the hyperbolic numbers (an extension of complex numbers, defined as z=x+h*y with h*h=1 and x,y real numbers) can be associated to space-time geometry as stated by the Lorentz transformations of…
We present an energy-conserving numerical scheme to solve the Vlasov-Maxwell (VM) system based on the regularized moment method proposed in [Z. Cai, Y. Fan, and R. Li. CPAM, 2014]. The globally hyperbolic moment system is deduced for the…
We address a class of schemes for the Euler equations with the following features: the space discretization is staggered, possible upwinding is performed with respect to the material velocity only and the internal energy balance is solved,…
Kinetic or Boltzmann schemes are interesting alternatives to the macroscopic numerical methods for solving the hyperbolic conservation laws of gas dynamics. They utilize the particle-based description instead of the wave propagation models.…
This paper is part of a program to combine a staggered time and staggered spatial discretization of continuum wave equations so that important properties of the continuum that are proved using vector calculus can be proven in an analogous…
We present a dynamical framework for modeling the motion of point-like charged particles, with or without mass, in general external electromagnetic fields. A key feature of this formulation is the treatment of time coordinate as a dynamical…
A second order accurate (in time) numerical scheme is proposed and analyzed for the Poisson-Nernst-Planck equation (PNP) system, reformulated as a non-constant mobility $H^{-1}$ gradient flow in the Energetic Variational Approach (EnVarA).…
The present work concerns the derivation of a numerical scheme to approximate weak solutions of the Euler equations with a gravitational source term. The designed scheme is proved to be fully well-balanced since it is able to exactly…
This paper describes a numerical scheme for multi-fluid hydrodynamics in the limit of small mass densities of the charged particles. The inertia of the charged particles can then be neglected, which makes it possible to write an evolution…
Numerical simulation of compressible fluid flows is performed using the Euler equations. They include the scalar advection equation for the density, the vector advection equation for the velocity and a given pressure dependence on the…
A new scheme for numerical integration of motion for classical systems composed of rigid polyatomic molecules is proposed. The scheme is based on a matrix representation of the rotational degrees of freedom. The equations of motion are…
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…
A Lagrangian formalism is used to study the motion of a spinning massive particle in Friedmann--Robertson--Walker and G\"odel spacetimes, as well as in a general Schwarzschild-like spacetime and in static spherically symmetric conformally…
Classical Finite Volume methods for multi-dimensional problems include stabilization (e.g.\ via a Riemann solver), that is derived by considering several one-dimensional problems in different directions. Such methods therefore ignore a…
We discuss structure-preserving numerical discretizations for repulsive and attractive Euler-Poisson equations that find applications in fluid-plasma and self-gravitation modeling. The scheme is fully discrete and structure preserving in…
We consider the dynamical behavior of Martin-L\"of random points in dynamical systems over metric spaces with a computable dynamics and a computable invariant measure. We use computable partitions to define a sort of effective symbolic…
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…