Related papers: HydroSyMBA: a 1D hydrocode coupled with an N-body …
The central region of a circumstellar disk is difficult to resolve in global numerical simulations of collapsing cloud cores, but its effect on the evolution of the entire disk can be significant. We use numerical hydrodynamics simulations…
Symplectic integrators are a foundation to the study of dynamical $N$-body phenomena, at scales ranging from from planetary to cosmological. These integrators preserve the Poincar\'e invariants of Hamiltonian dynamics. The $N$-body…
As the number of discovered extrasolar planets has been increasing, diversity of planetary systems requires studies of new formation scenarios. It is important to study satellite formation in circumplanetary disks, which is often viewed as…
Cosmological hydrodynamical simulations, while the current state-of-the art methodology for generating theoretical predictions for the large scale structures of the Universe, are among the most expensive simulation tools, requiring upwards…
We present the first 2D hydrodynamical finite volume simulations in which dust is fully coupled with the gas, including its back-reaction onto it, and at the same time the dust size is evolving according to coagulation and fragmentation…
We show that the gravothermal collapse of self-interacting dark matter (SIDM) halos can deviate from local thermodynamic equilibrium. As a consequence, the self-similar evolution predicted by the commonly adopted conducting fluid model can…
We present a series of 2-dimensional hydrodynamic simulations of massive disks around protostars. We simulate the same physical problem using both a `Piecewise Parabolic Method' (PPM) code and a `Smoothed Particle Hydrodynamic' (SPH) code,…
It is now generally agreed that multidimensional, multigroup, radiation hydrodynamics is an indispensable element of any realistic model of stellar-core collapse, core-collapse supernovae, and protoneutron star instabilities. We have…
Many of the most intriguing features, including spirals and cavities, in the current disc observations are found in binary systems like GG Tau, HD 142527 or HD 100453. Such features are evidence of the dynamic interaction between binary…
For a rigid model satellite, Chandrasekhar's dynamical friction formula describes the orbital evolution quite accurately, when the Coulomb logarithm is chosen appropriately. However, it is not known if the orbital evolution of a real…
Hydrodynamical cosmological simulations are a powerful tool for accurately predicting the properties of the intergalactic medium (IGM) and for producing mock skies that can be compared against observational data. However, the need to…
In this paper, the second in a series, we document the algorithms and solvers for compressible nonrelativistic hydrodynamics implemented in GenASiS (General Astrophysical Simulation System)---a new code being developed initially and…
We analyze the dynamical and thermal evolution of dark matter and ICM in hydrodynamical Tree-SPH simulations of galaxy clusters. Starting from a sample of 17 high-resolution objects, with virial mass ranging from 3 10^{14} to 1.7 10^{15}…
In my thesis I present a study of the dynamics and observational characteristics of massive circumstellar disks in two dimensions ($r$, $\phi$) using two complimentary hydrodynamic codes: a `Smoothed Particle Hydrodynamic' (SPH) code and a…
Observations of protoplanetary disks show that some characteristics seem recurrent, even in star formation regions that are physically distant such as surface mass density profiles varying as $r^{-1}$, or aspect ratios about 0.03 to 0.23.…
We describe an implicit general relativistic hydrodynamics code. The evolution equations are formulated in comoving coordinates. A conservative finite differencing of the Einstein equations is outlined, and artificial viscosity and…
The reproducibility of experiments is one of the main principles of the scientific method. However, numerical N-body experiments, especially those of planetary systems, are currently not reproducible. In the most optimistic scenario, they…
State-of-the-art one-dimensional (1D) stellar evolution codes rely on simplifying assumptions, such as mixing length theory, in order to describe superadiabatic convection. As a result, 1D stellar structure models do not correctly recover…
The recently developed code for N-body/hydrodynamics simulations in Modified Newtonian Dynamics (MOND), known as RAyMOND, is used to investigate the consequences of MOND on structure formation in a cosmological context, with a particular…
To study the post-main sequence evolution of the Solar system and exoplanetary systems more accurately and efficiently, we introduce two new features to REBOUNDx, an extended library for the N-body integrator REBOUND. The first is a…