Related papers: Accelerating self-gravitating hydrodynamics simula…
This paper extends the gas-kinetic scheme for one-dimensional inviscid shallow water equations (J. Comput. Phys. 178 (2002), pp. 533-562) to multidimensional gas dynamic equations under gravitational fields. Four important issues in the…
In gravitationally stratified fluids, length scales are normally much greater in the horizontal direction than in the vertical one. When modelling these fluids it can be advantageous to use the hydrostatic approximation, which filters out…
The equations of relativistic hydrodynamics are transformed so that steps forward in time preserves local simultaneity. In these variables, the space-time coordinates of neighboring points on the mesh are simultaneous according to co-moving…
We present a data-driven method for reconstructing the galactic acceleration field from phase-space measurements of stellar streams. Our approach is based on a flexible and differentiable fit to the stream in phase-space, enabling a direct…
The transport of energy through radiation is very important in many astrophysical phenomena. In dynamical problems the time-dependent equations of radiation hydrodynamics have to be solved. We present a newly developed…
The numerical study of relativistic magnetohydrodynamics (MHD) plays a crucial role in high-energy astrophysics, but unfortunately is computationally demanding, given the complex physics involved (high Lorentz factor flows, extreme…
Numerical simulations in modified gravity have commonly been performed under the quasi-static approximation -- that is, by neglecting the effect of time derivatives in the equation of motion of the scalar field that governs the fifth force…
We propose IMEX HDG-DG schemes for planar and spherical shallow water systems. Of interest is subcritical flow, where the speed of the gravity wave is faster than that of nonlinear advection. In order to simulate these flows efficiently, we…
The tree method is a widely implemented algorithm for collisionless $N$-body simulations in astrophysics well suited for GPU(s). Adopting hierarchical time stepping can accelerate $N$-body simulations; however, it is infrequently…
There is an increasing need to develop a two-dimensional (2D) water entry model including the slamming and transition stages for the 2.5-dimensional (2.5D) method being used on the take-off and water landing of seaplanes, and for the strip…
We present Cholla (Computational Hydrodynamics On ParaLLel Architectures), a new three-dimensional hydrodynamics code that harnesses the power of graphics processing units (GPUs) to accelerate astrophysical simulations. Cholla models the…
Shallow free surface flows are often characterized by both subdomains that require high modeling complexity and subdomains that can be sufficiently accurately modeled with low modeling complexity. Moreover, these subdomains may change in…
The migration of planets within their nascent protoplanetary disks is a fundamental process that shapes the final architecture of planetary systems. However, studying this phenomenon through direct hydrodynamical simulations is…
A method of calculating a new class of symmetries is presented for partial differential equations. The method give a new dynamical solution for an isothermal and cylindrically symmetric hydrodynamics equations under self-gravity. The…
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…
Simulating and predicting dynamics of quantum many-body systems is extremely challenging, even for state-of-the-art computational methods, due to the spread of entanglement across the system. However, in the long-wavelength limit, quantum…
The ideal MHD equations are a central model in astrophysics, and their solution relies upon stable numerical schemes. We present an implementation of a new method, which possesses excellent stability properties. Numerical tests demonstrate…
OCTO-TIGER is an astrophysics code to simulate the evolution of self-gravitating and rotat-ing systems of arbitrary geometry based on the fast multipole method, using adaptive mesh refinement. OCTO-TIGER is currently optimised to simulate…
Hydrodynamic simulations have been used to study accretion disks consisting of counterrotating components with an intervening shear layer(s). Configurations of this type can arise from the accretion of newly supplied counterrotating matter…
We present results based on an implementation of the Godunov Smoothed Particle Hydrodynamics (GSPH), originally developed by Inutsuka (2002), in the GADGET-3 hydrodynamic code. We first review the derivation of the GSPH discretization of…