Related papers: Fast and reliable symplectic integration for plane…
We present a statistical approximate solution of the bound, non-hierarchical three-body problem, and extend it to a general analysis of encounters between hard binary systems and single stars. Any such encounter terminates when one of the…
Astrophysical Challenges which demand the solution of the one million (or more) gravitating body problem are briefly discussed for the fields of cosmology, galactic nuclei and globular star clusters. Results from the classical three-body…
We propose a natural strategy to deal with compatible and incompatible binary questions, and with their time evolution. The strategy is based on the simplest, non-commutative, Hilbert space $\mathcal{H}=\mathbb{C}^2$, and on the (commuting…
Symplectic integrators separate a problem into parts that can be solved in isolation, alternately advancing these sub-problems to approximate the evolution of the complete system. Problems with a single, dominant mass can use mixed-variable…
The existence of explicit symplectic integrators for general nonseparable Hamiltonian systems is an open and important problem in both numerical analysis and computing in science and engineering, as explicit integrators are usually more…
Simulations of isolated binary mergers of galaxy clusters are a useful tool to study the evolution of these objects. For exceptionally massive systems they even represent the only viable way of simulation, because these are rare in typical…
The H(div)-conforming approach for the Brinkman equation is studied numerically, verifying the theoretical a priori and a posteriori analysis in previous work of the authors. Furthermore, the results are extended to cover a non-constant…
We derive a new criterion for estimating characteristic dynamical timescales in N-body simulations. The criterion uses the second, third, and fourth derivatives of particle positions: acceleration, jerk, and snap. It can be used for…
Numerical integrations of the Solar System have been carried out for decades. Their results have been used, for example, to determine whether the Solar System is chaotic, whether Mercury's orbit is stable, or to help discern Earth's climate…
We construct an advanced model for interacting multiple stellar systems in which we compute all trajectories with a numerical N-body integrator, namely the Bulirsch--Stoer from the SWIFT package. We can then derive various observables:…
A small deformation controlled by four free parameters to the Schwarzschild metric could be referred to a nonspinning black hole solution in alternative theories of gravity. Because such a non-Schwarzschild metric can be changed into a…
An efficient geometric integrator is proposed for solving the perturbed Kepler motion. This method is stable and accurate over long integration time, which makes it appropriate for treating problems in astrophysics, like solar system…
It is attempted to obtain the masses of the celestial bodies, the initial conditions of their motion, and the constant of gravitation, by a global parameter optimization. First, a numerical solution of the N-bodies problem for mass points…
In this work, we forecast the number of, and requirements on, N-body simulations needed to train hybrid effective field theory (HEFT) emulators for a range of use cases, using a hybrid of HMcode and perturbation theory as a surrogate model.…
Quasi-periodic responses composed of multiple base frequencies widely exist in science and engineering problems. The multiple harmonic balance (MHB) method is one of the most commonly used approaches for such problems. However, it is…
This study addresses optimal impulsive trajectory design within the Circular Restricted Three-Body Problem (CR3BP), presenting a global optimization-based approach to identify minimum $\Delta V$ transfers between periodic orbits, including…
Symplectic integrators are widely used for long-term integration of conservative astrophysical problems due to their ability to preserve the constants of motion; however, they cannot in general be applied in the presence of nonconservative…
We describe a new algorithm for the integration of self-gravitating fluid systems using SPH method. We split the Hamiltonian of a self-gravitating fluid system to the gravitational potential and others (kinetic and internal energies) and…
The motion of three bodies can be solved perturbatively when a tightly bound inner binary is orbited by a distant perturber, giving rise for example to the well-known Kozai-Lidov oscillations. We propose to study the relativistic…
Given the compact binary evolution problem of numerical relativity, in the finite-difference, block-based, adaptive mesh refinement context, choices must be made on how evolved fields are to be discretized. In GR-Athena++, the space-time…