Related papers: Exactly Conservative Integrators
A numerical scheme is developed for systems of conservation laws on manifolds which arise in high speed aerodynamics and magneto-aerodynamics. The systems are presented in an arbitrary coordinate system on the manifold and involve source…
This work aims to construct an efficient and highly accurate numerical method to address the time singularity at $t=0$ involved in a class of time-fractional parabolic integro-partial differential equations in one and two dimensions. The…
A Kepler solver is an analytical method used to solve a two-body problem. In this paper, we propose a new correction method by slightly modifying the Kepler solver. The only change to the analytical solutions is that the obtainment of the…
In this paper, we consider a fully-discrete approximation of an abstract evolution equation deploying a non-conforming spatial approximation and finite differences in time (Rothe-Galerkin method). The main result is the convergence of the…
Numerical methods that preserve geometric invariants of the system, such as energy, momentum or the symplectic form, are called geometric integrators. In this paper we present a method to construct symplectic-momentum integrators for…
We present a remarkable discretization of the classical Kepler problem which preserves its trajectories and all integrals of motion. The points of any discrete orbit belong to an appropriate continuous trajectory.
This paper presents a new strategy to deal with the excessive diffusion that standard finite volume methods for compressible Euler equations display in the limit of low Mach number. The strategy can be understood as using centered…
The explicit formulation of the general inverse problem on conservation laws is presented for the first time. In this problem one aims to derive the general form of systems of differential equations that admit a prescribed set of…
Incompressible two-dimensional flows such as the advection (Liouville) equation and the Euler equations have a large family of conservation laws related to conservation of area. We present two Eulerian numerical methods which preserve a…
The quasilinear theory describes the resonant interaction between particles and waves with two coupled equations: one for the evolution of the particle probability density function(\textit{pdf}), the other for the wave spectral energy…
We introduce a new methodology for adding localized, space-time smooth, artificial viscosity to nonlinear systems of conservation laws which propagate shock waves, rarefactions, and contact discontinuities, which we call the $C$-method. We…
The Active Flux method is a finite volume method for hyperbolic conservation laws that uses both cell averages and point values as degrees of freedom. Several versions of such methods are currently under development. We focus on third order…
Exponential integrability properties of numerical approximations are a key tool for establishing positive rates of strong and numerically weak convergence for a large class of nonlinear stochastic differential equations. It turns out that…
In this paper we present a formally fourth-order accurate hybrid-variable method for the Euler equations in the context of method of lines. The hybrid-variable (HV) method seeks numerical approximations to both cell-averages and nodal…
A novel structure-preserving numerical method to solve random hyperbolic systems of conservation laws is presented. The method uses a concept of generalized, measure-valued solutions to random conservation laws. This yields a linear partial…
In this paper we discuss energy conservation issues related to the numerical solution of the nonlinear wave equation. As is well known, this problem can be cast as a Hamiltonian system that may be autonomous or not, depending on the…
We present a systematic approach for the semiclassical treatment of many-body dynamics of interacting, open spin systems. Our approach overcomes some of the shortcomings of the recently developed discrete truncated Wigner approximation…
We present strongly convergent explicit and semi-implicit adaptive numerical schemes for systems of stiff stochastic differential equations (SDEs) where both the drift and diffusion are non-globally Lipschitz continuous. This stiffness may…
We introduce an algorithm to find possible constants of motion for a given replicator equation. The algorithm is inspired by Dirac geometry and a Hamiltonian description for the replicator equations with such constants of motion, up to a…
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…