Related papers: Octo-Tiger: A New, 3D Hydrodynamic Code for Stella…
[Abridged] We present the results of a highly parallel Kepler equation solver using the Graphics Processing Unit (GPU) on a commercial nVidia GeForce 280GTX and the "Compute Unified Device Architecture" programming environment. We apply…
The state-of-the-art in optimal control from timed temporal logic specifications, including Metric Temporal Logic (MTL) and Signal Temporal Logic (STL), is based on Mixed-Integer Convex Programming (MICP). The standard MICP approach is…
Hydrodynamical simulations of stellar interactions require stable models of stars as initial conditions. Such initial models, however, are difficult to construct for giant stars because of the wide range in spatial scales of the hydrostatic…
Mixed-precision methods combine low and high precision arithmetics to exploit low precision computational speed and high precision accuracy. Large ODE systems that contain many heterogeneous interactions lead to a high computational cost…
We describe the newly written code GADGET which is suitable both for cosmological simulations of structure formation and for the simulation of interacting galaxies. GADGET evolves self-gravitating collisionless fluids with the traditional…
Merging haloes with similar masses (i.e., major mergers) pose significant challenges for halo finders. We compare five halo finding algorithms' (AHF, HBT, Rockstar, SubFind, and VELOCIraptor) recovery of halo properties for both isolated…
Neutron star mergers are very violent events involving extreme physical processes: dynamic, strong-field gravity, large magnetic field, very hot, dense matter, and the copious production of neutrinos. Accurate modeling of such a system and…
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…
Galaxy mergers have been investigated for decades using smoothed particle hydrodynamics (SPH), but recent work highlighting inaccuracies inherent in the traditional SPH technique calls into question the reliability of previous studies. We…
We explore the effect of using different kinematic tracers (HI and H$\alpha$) on reconstructing the encounter parameters of the Mice major galaxy merger (NGC 4676A/B). We observed the Mice using the SparsePak Integral Field Unit (IFU) on…
According to the modern paradigms of software engineering, standard tasks are best accomplished by reusable open source libraries. We describe OpenOrbitalOptimizer: a reusable open source C++ library for the iterative solution of coupled…
Predicting the merger timescale ($\tau_{\rm merge}$) of merging dark matter halos, based on their orbital parameters and the structural properties of their hosts, is a fundamental problem in gravitational dynamics that has important…
This paper presents an accurate density computation approach for large dark matter simulations, based on a recently introduced phase-space tessellation technique and designed for massively parallel, heterogeneous cluster architectures. We…
Future gravitational wave detections of merging binary neutron star systems have the possibility to tightly constrain the equation of state of dense nuclear matter. In order to extract such constraints, gravitational waveform models need to…
A six-dimensional parameter space based on high-resolution numerical simulations of isolated binary galaxy collisions has been constructed to investigate the dynamical friction timescales, $\tau_{\rm mer}$, for major mergers. Our…
Convective overshooting in super asymptotic giant branch stars has been suggested to lead to the formation of hybrid white dwarfs with carbon-oxygen cores and oxygen-neon mantles. As the white dwarf cools, this core-mantle configuration…
Recent studies indicate that mergers of a white dwarf (WD) with a neutron star (NS) or a stellar-mass black hole (BH) may be a potential progenitor channel for certain merger-kind, but long-duration $\gamma$-ray bursts (GRBs), e.g., GRBs…
The demand for substantial increases in the spatial resolution of global weather- and climate- prediction models makes it necessary to use numerically efficient and highly scalable algorithms to solve the equations of large scale…
A crucial ingredient for numerically solving the 3D radiative transfer problem is the choice of the grid that discretizes the transfer medium. Many modern radiative transfer codes, whether using Monte Carlo or ray tracing techniques, are…
We present a new parallel supercomputer implementation of the Monte-Carlo method for simulating the dynamical evolution of globular star clusters. Our method is based on a modified version of Henon's Monte-Carlo algorithm for solving the…