Introducing a Hybrid Method of Radiative Transfer for Smoothed Particle Hydrodynamics
arXiv:0812.0304 · doi:10.1111/j.1365-2966.2008.14373.x
Abstract
A new means of incorporating radiative transfer into smoothed particle hydrodynamics (SPH) is introduced, which builds on the success of two previous methods - the polytropic cooling approximation as devised by Stamatellos et al (2007), and flux limited diffusion (e.g. Mayer et al 2007). This hybrid method preserves the strengths of its individual components, while removing the need for atmosphere matching or other boundary conditions to marry optically thick and optically thin regions. The code uses a non-trivial equation of state to calculate temperatures and opacities of SPH particles, which captures the effects of molecular hydrogen dissociation, atomic hydrogen ionisation, He0 and He+ ionisation, ice evaporation, dust sublimation, molecular absorption, bound-free and free-free transitions and electron scattering. The method is tested in several scenarios, including: (1) the evolution of a 0.07 solar mass protoplanetary disc surrounding a 0.5 solar mass star; (2) the collapse of a 1 solar mass protostellar cloud, and (3) the thermal relaxation of temperature fluctuations in a static homogeneous sphere.
11 pages, 22 figures, accepted for publication in MNRAS
References in corpus (8)
- SPLASH: An interactive visualisation tool for Smoothed Particle Hydrodynamics simulations
- Fragmentation of gravitationally unstable gaseous protoplanetary disks with radiative transfer
- The Thermal Regulation of Gravitational Instabilities in Protoplanetary Disks III. Simulations with Radiative Cooling and Realistic Opacities
- Radiative transfer and the energy equation in SPH simulations of star formation
- Numerical requirements for simulations of self gravitating and non-self gravitating disks
- The Thermal Regulation of Gravitational Instabilities in Protoplanetary Disks. IV. Simulations with Envelope Irradiation
- 3D Radiative Hydrodynamics for Disk Stability Simulations: A Proposed Testing Standard and New Results
- Flux-Limited Diffusion Approximation Models of Giant Planet Formation by Disk Instability