Related papers: Crash: A Block-Adaptive-Mesh Code for Radiative Sh…
In recent years, dynamical relativistic jet simulation techniques have progressed to a point where it is becoming possible to fully numerically resolve gamma-ray burst (GRB) blast-wave evolution across scales. However, the modeling of…
We compare two different codes for simulations of cosmological structure formation to investigate the sensitivity of hydrodynamical instabilities to numerics, in particular, the hydro solver and the application of adaptive mesh refinement…
This article condenses current endeavors and improvements in the expansion of applications of the DualSPHysics code to analyze heat transfer in a nuclear reactor core. This includes the essential conservation equations and certain physical…
We implemented sink particles in the adaptive mesh refinement (AMR) hydrodynamics code FLASH. Sink particles are created in regions of local gravitational collapse, and their trajectories and accretion can be followed over many dynamical…
Many astrophysical explosions, such as type Ia supernovae, classical novae, and X-ray bursts, are dominated by thermonuclear runaway. To model these processes accurately, one must evolve nuclear reactions concurrently with hydrodynamics. We…
Exascale supercomputing unleashes the potential for simulations of astrophysical systems with unprecedented resolution. Taking full advantage of this computing power requires the development of new algorithms and numerical methods that are…
We present a novel radiation hydrodynamics code, START, which is a smoothed particle hydrodynamics (SPH) scheme coupled with accelerated radiative transfer. The basic idea for the acceleration of radiative transfer is parallel to the tree…
Solar flares involve complex processes that are coupled and span a wide range of temporal, spatial, and energy scales. Modeling such processes self-consistently has been a challenge in the past. Here we present results from simulations that…
Recent three-dimensional radiative hydrodynamics simulations of protoplanetary disks report disparate disk behaviors, and these differences involve the importance of convection to disk cooling, the dependence of disk cooling on metallicity,…
There is great interest in numerical relativity simulations involving matter due to the likelihood that binary compact objects involving neutron stars will be detected by gravitational wave observatories in the coming years, as well as to…
We present a study of cooling in radiative shocks simulated with smoothed particle hydrodynamics (SPH) and adaptive mesh refinement codes. We obtain a similarity solution for a shock-tube problem in the presence of radiative cooling, and…
A module for the ZEUS-2D code is described which may be used to solve the equations of radiation hydrodynamics to order unity in v/c, in the flux-limited diffusion (FLD) approximation. In this approximation, the tensor Eddington factor f…
The relativistic hydrodynamics (RHD) equations can give rise to solutions which have shocks, contact discontinuities, and other sharp structures, which interact and evolve over time. Capturing these sharp waves effectively requires a mesh…
A new combined sub-filter scale turbulence and shock-capturing model is developed for high-order finite volume numerics, extending previous work to unstructured solvers. Block Spectral Stresses (BSS) method relies on the spectra of the…
We present an implementation of an adaptive ray tracing (ART) module in the Athena hydrodynamics code that accurately and efficiently handles the radiative transfer involving multiple point sources on a three-dimensional Cartesian grid. We…
Modeling the dynamics of most astrophysical structures requires an adequate description of the radiation-matter interaction. Several numerical (magneto)hydrodynamics codes were upgraded with a radiation module to fulfill this request.…
The smoothed particle hydrodynamics (SPH) method has been increasingly used to study fluid problems in recent years; but its computational cost can be high if high resolution is required. In this study, an adaptive resolution method based…
High-order implicit shock tracking (fitting) is a class of high-order, optimization-based numerical methods to approximate solutions of conservation laws with non-smooth features by aligning elements of the computational mesh with…
One of the key questions to understanding the efficiency of diffusive shock acceleration of the cosmic rays (CRs) is the injection process from thermal particles. A self-consistent injection model based on the interactions of the…
We develop a new relativistic radiation hydrodynamics code based on the Monte-Carlo algorithm. In this code, we implement a new scheme to achieve the second-order accuracy in time in the limit of a large packet number for solving the…