Related papers: The accuracy of mutual potential approximations in…
A new technique that utilizes surface integrals to find the force, torque and potential energy between two non-spherical, rigid bodies is presented. The method is relatively fast, and allows us to solve the full rigid two-body problem for…
We review various approaches to approximating general relativistic effects in hydrodynamic simulations of stellar core collapse and post-bounce evolution. Different formulations of a modified Newtonian gravitational potential are presented.…
Asteroid binaries found amongst the Near-Earth objects are believed to have formed from rotational fission. In this paper, we aim to study the dynamical evolution of asteroid systems the moment after fission. The initial condition is…
In asteroid rendezvous missions, the dynamical environment near the asteroid's surface should be made clear prior to the mission launch. However, most of the asteroids have irregular shapes, which lower the efficiency of calculating their…
Spin-orbit coupling is widespread in binary asteroid systems and it has been widely studied for the case of ellipsoidal secondary. Due to angular momentum exchange, dynamical coupling is stronger when the orbital and rotational angular…
The determination of the gravitational potential by the polyhedral method is revisited in the case where the surface of a body is composed of triangular facets. Based upon six test-shapes of astrophysical interest (sphere, spheroid,…
Accurate gravity field calculations are necessary for landing on planets, moons, asteroids, minimoons, or other irregularly shaped bodies, but current methods become increasingly inaccurate and slow near the surface. We present high…
Simulation of materials at the atomistic level is an important tool in studying microscopic structure and processes. The atomic interactions necessary for the simulation are correctly described by Quantum Mechanics. However, the…
We develop a fourth-order numerical integrator to simulate the coupled spin and orbital motions of two rigid bodies having arbitrary mass distributions under the influence of their mutual gravitational potential. We simulate the dynamics of…
A number of problems arise when long-range forces, such as those governed by Bessel functions, are used in particle-particle simulations. If a simple cut-off for the interaction is used, the system may find an equilibrium configuration at…
This paper presents a robust linear method for impact probability estimation of near-Earth asteroids with the Earth. This method is a significantly modified and improved method, which uses a special curvilinear coordinate system associated…
The paper presents a numerical implementation of the gravitational N-body problem with contact interactions between non-spherically shaped bodies. The work builds up on a previous implementation of the code and extends its capabilities. The…
The minimum orbital intersection distance is used as a measure to assess potential close approaches and collision risks between astronomical objects. Methods to calculate this quantity have been proposed in several previous publications.…
We provide a generalized discussion of tidal evolution to arbitrary order in the expansion of the gravitational potential between two spherical bodies of any mass ratio. To accurately reproduce the tidal evolution of a system at separations…
We construct time-evolving gravitational potential models for a Milky Way-mass galaxy from the FIRE-2 suite of cosmological-baryonic simulations using basis function expansions. These models capture the angular variation with spherical…
Direct gravitational simulations of n-body systems have a time complexity O(n^2), which gets computationally expensive as the number of bodies increases. Distributing this workload to multiple cores significantly speeds up the computation…
Astrophysical black holes could be nearly extremal (that is, rotating nearly as fast as possible); therefore, nearly extremal black holes could be among the binaries that current and future gravitational-wave observatories will detect.…
We simulate the collapse of a turbulent gas cloud with the particle-based code Gadget2. We choose two sub-clouds, formed by those particles located around the centers $\vec{r}_L$ and $\vec{r}_R$ and within the radii $r_L$ and $r_R$,…
There are two popular ways to speed up simulations of planet formation via increasing the collision probability: ({\it i}) confine motion to 2D, ({\it ii}) artificially enhance the physical radii of the bodies by an expansion factor. In…
Bayesian adaptive designs enable flexible clinical trials by adapting features based on accumulating data. Among these, Bayesian Response-Adaptive Randomization (BRAR) skews patient allocation towards more promising treatments based on…