Related papers: Revisiting high-order Taylor methods for astrodyna…
We present IAS15, a 15th-order integrator to simulate gravitational dynamics. The integrator is based on a Gau\ss-Radau quadrature and can handle conservative as well as non-conservative forces. We develop a step-size control that can…
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…
We have designed and developed, from scratch, a global circulation model named THOR that solves the three-dimensional non-hydrostatic Euler equations. Our general approach lifts the commonly used assumptions of a shallow atmosphere and…
We consider the Hirota equation (the discrete analog of the generalized Toda system) over a finite field. We present the general algebro-geometric method of construction of solutions of the equation. As an example we construct analogs of…
The formation and evolution of protoplanetary systems, the breeding grounds of planet formation, is a complex dynamical problem that involves many orders of magnitudes. To serve this purpose, we present a new hybrid algorithm that combines…
We introduce a hybrid method to couple continuous Galerkin finite element methods and high-order finite difference methods in a nonconforming multiblock fashion. The aim is to optimize computational efficiency when complex geometries are…
We study the numerical approximation by space-time finite element methods of a multi-physics system coupling hyperbolic elastodynamics with parabolic transport and modeling poro- and thermoelasticity. The equations are rewritten as a…
We introduce a general framework for the construction of well-balanced finite volume methods for hyperbolic balance laws. We use the phrase well-balancing in a broader sense, since our proposed method can be applied to exactly follow any…
We present new splitting methods designed for the numerical integration of near-integrable Hamiltonian systems, and in particular for planetary N-body problems, when one is interested in very accurate results over a large time span. We…
We present WHFast, a fast and accurate implementation of a Wisdom-Holman symplectic integrator for long-term orbit integrations of planetary systems. WHFast is significantly faster and conserves energy better than all other Wisdom-Holman…
We present new almost time-reversible integrators for solution of planetary systems consisting of "planets" and a dominant mass ("star"). The algorithms can be considered adaptive generalizations of the Wisdom--Holman method, in which all…
OCTO-TIGER is an astrophysics code to simulate the evolution of self-gravitating and rotat-ing systems of arbitrary geometry based on the fast multipole method, using adaptive mesh refinement. OCTO-TIGER is currently optimised to simulate…
Our focus is on simulating the dynamics of non-interacting particles including the effects of an external potential, which, under certain assumptions, can be formally described by the Dean-Kawasaki equation. The Dean-Kawasaki equation can…
The iterative rational Krylov algorithm (\textsf{IRKA}) is a popular approach for producing locally optimal reduced-order $\mathcal{H}_2$-approximations to linear time-invariant (LTI) dynamical systems. Overall, \textsf{IRKA} has seen…
A novel numerical approach to solving the shallow-water equations on the sphere using high-order numerical discretizations in both space and time is proposed. A space-time tensor formalism is used to express the equations of motion…
Calculating the long term solution of ordinary differential equations, such as those of the $N$-body problem, is central to understanding a wide range of dynamics in astrophysics, from galaxy formation to planetary chaos. Because generally…
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…
Context. Integrating the motion of stars in a smoothed potential is necessary in many stellar and galactic studies. Previous works have often used numerical integrators that alternate between linear drifts and velocity kicks (such as the…
Gravitational wave signals from extreme mass ratio inspirals are a key target for space-based gravitational wave detectors. These systems are typically modeled as a distributionally-forced Teukolsky equation, where the smaller black hole is…
A space discretization of an integrable long wave-short wave interaction model, called the Yajima-Oikawa system, was proposed in the recent paper arXiv:1509.06996 using the Hirota bilinear method (see also…