Related papers: Multi-Point Hermite Methods for the N-Body Problem
Elegant integration schemes of second and fourth order for simulations of rigid body systems are presented which treat translational and rotational motion on the same footing. This is made possible by a recent implementation of the exact…
Hamiltonian systems are known to conserve the Hamiltonian function, which describes the energy evolution over time. Obtaining a numerical spatio-temporal scheme that accurately preserves the discretized Hamiltonian function is often a…
We design an accurate orbital integration scheme for the general N-body problem preserving all the conserved quantities but the angular momentum.This scheme is based on the chain concept (Mikkola & Aarseth 1993) and is regarded as an…
This article is concerned with Gaussian process quadratures, which are numerical integration methods based on Gaussian process regression methods, and sigma-point methods, which are used in advanced non-linear Kalman filtering and smoothing…
This paper is concerned with the theory, construction and application of variable-stepsize implicit Peer two-step methods that are super-convergent for variable stepsizes, i.e., preserve their classical order achieved for uniform stepsizes…
This work investigates finite differences and the use of interpolation models to obtain approximations to the first and second derivatives of a function. Here, it is shown that if a particular set of points is used in the interpolation…
Energy-conserving Hermite methods for solving Maxwell's equations in dielectric and dispersive media are described and analyzed. In three space dimensions methods of order $2m$ to $2m+2$ require $(m+1)^3$ degrees-of-freedom per node for…
Quantum master equations are common tools to describe the dynamics of many-body systems open to an environment. Due to the interaction with the latter, even for the case of non-interacting electrons, the computational cost to solve these…
In this paper, we propose a new set of midpoint-based high-order discretization schemes for computing straight and mixed nonlinear second derivative terms that appear in the compressible Navier-Stokes equations. Firstly, we detail a set of…
We propose a time-adaptive predictor/multi-corrector method to solve hyperbolic partial differential equations, based on the generalized-$\alpha$ scheme that provides user-control on the numerical dissipation and second-order accuracy in…
Quantum many-body theory has witnessed tremendous progress in various fields, ranging from atomic and solid-state physics to quantum chemistry and nuclear structure. Due to the inherent computational burden linked to the ab initio treatment…
First-order variational equations are widely used in N-body simulations to study how nearby trajectories diverge from one another. These allow for efficient and reliable determinations of chaos indicators such as the Maximal Lyapunov…
Rapid satellite-to-site visibility determination is of great significance to coverage analysis of satellite constellations as well as onboard mission planning of autonomous spacecraft. This paper presents a novel self-adaptive Hermite…
In this note we consider splitting methods based on linear multistep methods and stabilizing corrections. To enhance the stability of the methods, we employ an idea of Bruno & Cubillos (2016) who combine a high-order extrapolation formula…
By a high-order numerical homogenization method, a heterogeneous multiscale scheme was developed in Jin & Li (2022) for evolving differential equations containing two time scales. In this paper, we further explore the technique to propose…
Simulations of many rigid bodies colliding with each other sometimes yield particularly interesting results if the colliding objects differ significantly in size and are non-spherical. The most expensive part within such a simulation code…
We developed a Keplerian-based Hamiltonian splitting for solving the gravitational $N$-body problem. This splitting allows us to approximate the solution of a general $N$-body problem by a composition of multiple, independently evolved…
We present weak approximations schemes of any order for the Heston model that are obtained by using the method developed by Alfonsi and Bally (2021). This method consists in combining approximation schemes calculated on different random…
The neon isotopic chain displays a rich phenomenology, ranging from clustering in the ground-state of the self-conjugate doubly open-shell stable $^{20}$Ne isotope to the physics of the island of inversion around the neutron-rich $^{30}$Ne…
The Fast Multipole Method (FMM) offers an acceleration for pairwise interaction calculation, known as $N$-body problems, from $\mathcal{O}(N^2)$ to $\mathcal{O}(N)$ with $N$ particles. This has brought dramatic increase in the capability of…