Related papers: SEREN - A new SPH code for star and planet formati…
Stellar physics and evolution calculations enable a broad range of research in astrophysics. Modules for Experiments in Stellar Astrophysics (MESA) is a suite of open source libraries for a wide range of applications in computational…
We describe a newly developed cosmological hydrodynamics code based on the weighted essentially non-oscillatory (WENO) schemes for hyperbolic conservation laws. High order finite difference WENO schemes are designed for problems with…
An improved implementation of an N-body code for simulating collisionless cosmological dynamics is presented. TPM (Tree-Particle-Mesh) combines the PM method on large scales with a tree code to handle particle-particle interactions at small…
We compare two state-of-the-art numerical codes to study the overall accuracy in modeling the intergalactic medium and reproducing Lyman-$\alpha$ forest observables for DESI and high-resolution data sets. The codes employ different…
We substantially update the capabilities of the open-source software instrument Modules for Experiments in Stellar Astrophysics (MESA). MESA can now simultaneously evolve an interacting pair of differentially rotating stars undergoing…
We describe numerical methods for incorporating gas dynamics into cosmological simulations and present illustrative applications to the cold dark matter (CDM) scenario. Our evolution code, a version of TreeSPH (Hernquist \& Katz 1989)…
We describe an open source GPU implementation of a hybrid symplectic N-body integrator, GENGA (Gravitational ENcounters with Gpu Acceleration), designed to integrate planet and planetesimal dynamics in the late stage of planet formation and…
A fast $N$-body code has been developed for simulating a stellar disk embedded in a live dark matter halo. In generating its Poisson solver, a self-consistent field (SCF) code which inherently possesses perfect scalability is incorporated…
We present a 'two-fluid' implementation of dust in smoothed particle hydrodynamics (SPH) in the test particle limit. The scheme is able to handle both short and long stopping times and reproduces the short friction time limit, which is not…
We describe a new algorithm for including the dynamical effects of ionizing radiation in SPH simulations, and we present several examples of how the algorithm can be applied to problems in star formation. We use the HEALPix software to…
A unified Smoothed Particle Hydrodynamics (SPH) framework is proposed to simulate interaction dynamics involving thin shells modeled by a reduced-dimensional, single-layer particle discretization, as opposed to full-dimensional SPH solids.…
We describe and demonstrate a method for increasing the resolution locally in a Smoothed Particle Hydrodynamic (SPH) simulation, by splitting particles. We show that in simulations of self-gravitating collapse (of the sort which are…
We develop a general method for the self consistent calculation of the hydrodynamics of an astrophysical object irradiated by a radiation field with an arbitrary strength and spectral energy distribution (SED). Using the XSTAR…
We present and test a general-purpose code, called PPASPH, for evolving self-gravitating fluids in astrophysics, both with and without a collisionless component. In PPASPH, hydrodynamical properties are computed by using the SPH (Smoothed…
We present a new radiative transfer method (SPH-M1RT) that is coupled dynamically with smoothed particle hydrodynamics (SPH). We implement it in the (task-based parallel) SWIFT galaxy simulation code but it can be straightforwardly…
High-accuracy numerical simulations of merging neutron stars play an important role in testing and calibrating the waveform models used by gravitational wave observatories. Obtaining high-accuracy waveforms at a reasonable computational…
We report on a new open-source, user-friendly numerical relativity code package called SENR/NRPy+. Our code extends previous implementations of the BSSN reference-metric formulation to a much broader class of curvilinear coordinate systems,…
We describe a strategy for code modernisation of Gadget, a widely used community code for computational astrophysics. The focus of this work is on node-level performance optimisation, targeting current multi/many-core IntelR architectures.…
In this paper the problem of consistency of smoothed particle hydrodynamics (SPH) is solved. A novel error analysis is developed in $n$-dimensional space using the Poisson summation formula, which enables the treatment of the kernel and…
The fragmentation of shocked flows in a thermally bistable medium provides a natural mechanism to form turbulent cold clouds as precursors to molecular clouds. Yet because of the large density and temperature differences and the range of…