Related papers: Construction of second-order six-dimensional Hamil…
We show in this paper that third- and fourth-order low storage Runge-Kutta algorithms can be built specifically for quadratic nonlinear operators, at the expense of roughly doubling the time needed for evaluating the temporal derivatives.…
Direct observations of compact objects, in the form of radiation spectra, gravitational waves from VIRGO/LIGO, and forthcoming direct imaging, are currently one of the primary source of information on the physics of plasmas in extreme…
In this paper we study arbitrarily high-order energy-conserving methods for simulating the dynamics of a charged particle. They are derived and studied within the framework of Line Integral Methods (LIMs), previously used for defining…
An explicit high-order noncanonical symplectic algorithm for ideal two-fluid systems is developed. The fluid is discretized as particles in the Lagrangian description, while the electromagnetic fields and internal energy are treated as…
Time-reversible symplectic methods, which are precisely compatible with Liouville's phase-volume-conservation theorem, are often recommended for computational simulations of Hamiltonian mechanics. Lack of energy drift is an apparent…
Energy and decay estimates for the wave equation on the exterior region of slowly rotating Kerr spacetimes are proved. The method used is a generalization of the vector-field method, which allows the use of higher-order symmetry operators.…
We take into account the dynamics of a complete third post-Newtonian conservative Hamiltonian of two spinning black holes, where the orbital part arrives at the third post-Newtonian precision level and the spin-spin part with the spin-orbit…
Many interesting applications of hyperbolic systems of equations are stiff, and require the time step to satisfy restrictive stability conditions. One way to avoid small time steps is to use implicit time integration. Implicit integration…
Black hole perturbation theory beyond second order is not well understood because typically one defines the meaning of gauge invariance order by order which is ambiguous. In this series of works we therefore developed a new approach which…
We present a numerical scheme for solving a sixth-order Cahn-Hilliard type equation that captures the dynamics of phase transitions in a ternary mixture consisting of two immiscible fluids and a surface active molecule that is amphiphilic.…
First-order fully implicit as well as implicit--explicit schemes for coupled elliptic-parabolic systems are discussed in [Ern and Meunier, ESAIM: M2AN, 2009] and [Altmann et al., Math.\ Comp., 2021], respectively. The extension of the…
When dealing with stiff conservation laws, explicit time integration forces to employ very small time steps, due to the restrictive CFL stability condition. Implicit methods offer an alternative, yielding the possibility to choose the time…
We give a possible splitting method to a Hamiltonian for the description of charged particles moving around the Reissner-Nordstrom-(anti)-de Sitter black hole with an external magnetic field. This Hamiltonian can be separated into six…
A new methodology is developed to integrate numerically the equations of motion for classical many-body systems in molecular dynamics simulations. Its distinguishable feature is the possibility to preserve, independently on the size of the…
A high-order well-balanced scheme for the Euler equations with gravitation is presented. The scheme is able to preserve a spatially high-order accurate discrete representation of a large class of hydrostatic equilibria. It is based on a…
This paper deals with the numerical solution of conservation laws in the two dimensional case using a novel compact implicit time discretization that enables applications of fast algebraic solvers. We present details for the second order…
We investigate linear dynamical systems of second order. Uncertainty quantification is applied, where physical parameters are substituted by random variables. A stochastic Galerkin method yields a linear dynamical system of second order…
This study proposes an extension of the high-order compact gas-kinetic scheme (CGKS) to compressible flow simulation in an arbitrary Lagrangian-Eulerian (ALE) formulation in unstructured mesh. The ALE method is achieved by subdividing…
We present a new numerical code, ECHO, based on an Eulerian Conservative High Order scheme for time dependent three-dimensional general relativistic magnetohydrodynamics (GRMHD) and magnetodynamics (GRMD). ECHO is aimed at providing a…
The Kerr-MOG black hole has recently attracted significant research attention and has been extensively applied in various fields. To accurately characterize the long-term dynamical evolution of charged particles around Kerr-MOG black hole,…