January 27, 2006
We have been constructed a brand-new radiation hydrodynamics solver based upon Smoothed Particle Hydrodynamics (SPH), which works on parallel computer system. The code is designed to investigate the formation and evolution of the first generation objects at $z \gtrsim 10$, where the radiative feedback from various sources play important roles. The code can compute the fraction of chemical species e, H$^+$, H, H$^-$, H$_2$, and H$_2^+$ by fully implicit time integration. It also can deal with multiple sources of ionizing radiation, as well as the radiation at Lyman-Werner band. We compare the results for a few test calculations with the results of one dimensional simulations, in which we find good agreements with each other. We also evaluate the speedup by parallelization, that is found to be almost ideal, as far as the number of sources is comparable to the number of processors.
Similar papers 1
May 28, 2010
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 algorithm that is hitherto used to speed up the gravitational force calculation in an N-body system. It is demonstrated that the radiative transfer calculations can be dramatically accelerated, where the computational time is scaled as Np lo...
February 26, 2008
We introduce SPHRAY, a Smoothed Particle Hydrodynamics (SPH) ray tracer designed to solve the 3D, time dependent, radiative transfer (RT) equations for arbitrary density fields. The SPH nature of SPHRAY makes the incorporation of separate hydrodynamics and gravity solvers very natural. SPHRAY relies on a Monte Carlo (MC) ray tracing scheme that does not interpolate the SPH particles onto a grid but instead integrates directly through the SPH kernels. Given initial conditions ...
February 13, 2008
We present TRAPHIC, a novel radiative transfer scheme for Smoothed Particle Hydrodynamics (SPH) simulations. TRAPHIC is designed for use in simulations exhibiting a wide dynamic range in physical length scales and containing a large number of light sources. It is adaptive both in space and in angle and can be employed for application on distributed memory machines. The commonly encountered computationally expensive scaling with the number of light sources in the simulation is...
February 13, 2009
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 tessellate the sky and to solve the equation of ionization equilibrium along a ray towards each of the resulting tesserae. We exploit the hierarchical nature of HEALPix to make the algorithm adaptive, so that fine angular resolution is invoked ...
February 11, 2000
A new method for the inclusion of ionizing radiation from uniform radiation fields into 3D Smoothed Particle Hydrodynamics (SPHI) simulations is presented. We calculate the optical depth for the Lyman continuum radiation from the source towards the SPHI particles by ray-tracing integration. The time-dependent ionization rate equation is then solved locally for the particles within the ionizing radiation field. Using test calculations, we explore the numerical behaviour of the...
June 17, 2004
We describe the implementation and testing of a smoothed particle hydrodynamics (SPH) code that solves the equations of radiation hydrodynamics in the flux-limited diffusion (FLD) approximation. The SPH equations of radiation hydrodynamics for an explicit integration scheme are derived and tested. We also discuss the implementation of an implicit numerical scheme for solving the radiation equations that allows the system to be evolved using timesteps much longer than the expl...
February 16, 2021
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 implemented in other SPH codes. Our moment-based method simultaneously solves the radiation energy and flux equations in SPH, making it adaptive in space and time. We modify the M1 closure relation to stabilize radiation fronts in the optically thi...
November 7, 2017
In this paper we introduce a new radiative transfer code SPRAI (Simplex Photon Radiation in the Arepo Implementation) based on the SimpleX radiation transfer method. This method, originally used only for post-processing, is now directly integrated into the Arepo code and takes advantage of its adaptive unstructured mesh. Radiated photons are transferred from the sources through the series of Voronoi gas cells within a specific solid angle. From the photon attenuation we deriv...
January 30, 2020
The radiation hydrodynamics equations for smoothed particle hydrodynamics are derived by operator splitting the radiation and hydrodynamics terms, including necessary terms for material motion, and discretizing each of the sets of equations separately in time and space. The implicit radiative transfer discussed in the first paper of this series is coupled to explicit smoothed particle hydrodynamics. The result is a multi-material meshless radiation hydrodynamics code with arb...
April 26, 2013
We present a new implementation of radiation hydrodynamics (RHD) in the adaptive mesh refinement (AMR) code RAMSES. The multi-group radiative transfer (RT) is performed on the AMR grid with a first-order Godunov method using the M1 closure for the Eddington tensor, and is coupled to the hydrodynamics via non-equilibrium thermochemistry of hydrogen and helium. This moment-based approach has the large advantage that the computational cost is independent of the number of radiati...