Smoothed Particle Hydrodynamics coupled with Radiation Transfer

We have developed an accurate, one-dimensional, spherically symmetric, Lagrangian hydrodynamics/gravity code, designed to study the effects of radiative cooling and photo-ionization on the formation of protogalaxies. We examine the ability of collapsing perturbations to cool within the age of the universe. In contrast to some studies based on order-of-magnitude estimates, we find that cooling arguments alone cannot explain the sharp upper cutoff observed in the galaxy luminos...

The chemical and photometric evolution of star forming disk galaxies is investigated. Numerical simulations of the complex gasdynamical flows are based on our own coding of the Chemo - Dynamical Smoothed Particle Hydrodynamical (CD - SPH) approach, which incorporates the effects of star formation. As a first application, the model is used to describe the chemical and photometric evolution of a disk galaxy like the Milky Way.

We propose an explicit-implicit scheme for numerically solving Special Relativistic Radiation Hydrodynamic (RRHD) equations, which ensures a conservation of total energy and momentum (matter and radiation). In our scheme, 0th and 1st moment equations of the radiation transfer equation are numerically solved without employing a flux-limited diffusion (FLD) approximation. For an hyperbolic term, of which the time scale is the light crossing time when the flow velocity is compar...

We describe the structure and implementation of a radiation hydrodynamic solver for MANGA, the moving-mesh hydrodynamics module of the large-scale parallel code, Charm N-body GrAvity solver (ChaNGa). We solve the equations of time dependent radiative transfer using a reduced speed of light approximation following the algorithm of Jiang et al (2014). By writing the radiative transfer equations as a generalized conservation equation, we solve the transport part of these equatio...

We describe an extension of the Enzo code to enable fully-coupled radiation hydrodynamical simulation of inhomogeneous reionization in large $\sim (100 Mpc)^3$ cosmological volumes with thousands to millions of point sources. We solve all dynamical, radiative transfer, thermal, and ionization processes self-consistently on the same mesh, as opposed to a postprocessing approach which coarse-grains the radiative transfer. We do, however, employ a simple subgrid model for star f...

Radiation transport plays a crucial role in star formation models, as certain questions within this field cannot be accurately addressed without taking it into account. Given the high complexity of the interstellar medium from which stars form, numerical simulations are frequently employed to model the star formation process. This study reviews recent methods for incorporating radiation transport into star formation simulations, discussing them in terms of the used algorithms...

We report first results of an implementation of a chemical model in a cosmological code, based on the Smoothed Particle Hydrodynamics (SPH) technique. We show that chemical SPH simulations are a promising tool to provide clues for the understanding of the chemical properties of galaxies in relation to their formation and evolution in a cosmological framework.

We develop a new numerical scheme for solving the radiative transfer equation in a spherically symmetric system. This scheme does not rely on any kind of diffusion approximation and it is accurate for optically thin, thick, and intermediate systems. In the limit of a homogeneously distributed extinction coefficient, our method is very accurate and exceptionally fast. We combine this fast method with a slower, but more generally applicable method to describe realistic problems...

The advantages and disadvantages of two approaches to astrophysical hydrodynamics, Smoothed Particle Hydrodynamics and Adaptive Mesh Refinement, are briefly discussed together with some current problems of computational hydrodynamics.

We describe a novel adaptive ray tracing scheme to solve the equation of radiative transfer around point sources in hydrodynamical simulations. The angular resolution adapts to the local hydrodynamical resolution and hence is of use for adaptive meshes as well as adaptive smooth particle hydrodynamical simulations. Recursive creation of rays ensures ease of implementation. The multiple radial integrations needed to solve the time dependent radiative transfer are sped up signi...