Related papers: High order volume-preserving algorithms for relati…
We present a high-order conservative, positivity-preserving, and non-oscillatory scheme for solving the Vlasov equation. The scheme attains formal fifth-order accuracy through a convex combination of positive and non-oscillatory polynomials…
In this chapter, we aim at presenting the basic techniques necessary to go beyond the widely accepted paradigm of second-order numerics. We specifically focus on finite-volume schemes for hyperbolic conservation laws occuring in fluid…
In this paper, we introduce a Lagrange multiplier approach to construct linearly implicit energy-preserving schemes of arbitrary order for general Hamiltonian PDEs. Unlike the widely used auxiliary variable methods, this novel approach does…
The construction of high-order structure-preserving numerical schemes to solve hyperbolic conservation laws has attracted a lot of attention in the last decades and various different ansatzes exist. In this paper, we compare three…
We introduce a new electrostatic particle-in-cell algorithm capable of using large timesteps compared to particle gyro-period under a uniform external magnetic field. The algorithm extends earlier electrostatic fully implicit PIC…
In this paper we propose the first better than second order accurate method in space and time for the numerical solution of the resistive relativistic magnetohydrodynamics (RRMHD) equations on unstructured meshes in multiple space…
In this work, we propose a high-order multiscale method for an elliptic model problem with rough and possibly highly oscillatory coefficients. Convergence rates of higher order are obtained using the regularity of the right-hand side only.…
Computing the volume of a polytope in high dimensions is computationally challenging but has wide applications. Current state-of-the-art algorithms to compute such volumes rely on efficient sampling of a Gaussian distribution restricted to…
In this paper we design high-order positivity-preserving approximation schemes for an integro-differential model describing photochemical reactions. Specifically, we introduce and analyze three classes of dynamically consistent methods,…
The generalized-$\alpha$ time-marching method provides second-order accuracy in time and controls the numerical dissipation in the high-frequency region of the discrete spectrum. This method includes a wide range of time integrators. We…
We develop high-order numerical schemes to solve random hyperbolic conservation laws using linear programming. The proposed schemes are high-order extensions of the existing first-order scheme introduced in [{\sc S. Chu, M. Herty, M.…
In this paper, we define arbitrarily high-order energy-conserving methods for Hamiltonian systems with quadratic holonomic constraints. The derivation of the methods is made within the so-called line integral framework. Numerical tests to…
We consider high-order splitting schemes for large-scale differential Riccati equations. Such equations arise in many different areas and are especially important within the field of optimal control. In the large-scale case, it is critical…
Second-order Volume-preserving algorithms (VPAs) for simulating charged particle motion in electromagnetic fields have been generalized to a rotating angle formulation by using the matrix decomposition methods. Based on this method, the…
Structure-preserving algorithms and algorithms with uniform error bound have constituted two interesting classes of numerical methods. In this paper, we blend these two kinds of methods for solving nonlinear Hamiltonian systems with highly…
A basic leapfrog integrator and its energy-preserving and variational / symplectic variants are proposed and studied for the numerical integration of the equations of motion of relativistic charged particles in an electromagnetic field. The…
Symplectic integration methods based on operator splitting are well established in many branches of science. For Hamiltonian systems which split in more than two parts, symplectic methods of higher order have been studied in detail only for…
We consider the Vlasov-Poisson equation in a Hamiltonian framework and derive new time splitting methods based on the decomposition of the Hamiltonian functional between the kinetic and electric energy. Assuming smoothness of the solutions,…
Linear differential equations are ubiquitous in science and engineering. Quantum computers can simulate quantum systems, which are described by a restricted type of linear differential equations. Here we extend quantum simulation algorithms…
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…