اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدRadiative transfer and the energy equation in SPH simulations of star formation
60
0
0.0
(
0
)
تأليف
D. Stamatellos
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
We introduce and test a new and highly efficient method for treating the thermal and radiative effects influencing the energy equation in SPH simulations of star formation. The method uses the density, temperature and gravitational potential of each particle to estimate a mean optical depth, which then regulates the particles heating and cooling. The method captures -- at minimal computational cost -- the effects of (i) the rotational and vibrational degrees of freedom of H2, H2 dissociation, H0 ionisation, (ii) opacity changes due to ice mantle melting, sublimation of dust, molecular lines, H-, bound-free and free-free processes and electron scattering; (iv) external irradiation; and (v) thermal inertia. The new algorithm reproduces the results of previous authors and/or known analytic solutions. The computational cost is comparable to a standard SPH simulation with a simple barotropic equation of state. The method is easy to implement, can be applied to both particle- and grid-based codes, and handles optical depths 0<tau<10^{11}.
قيم البحث
اقرأ أيضاً
Abridged: We simulate a massive galaxy cluster in a LCDM Universe using three different approaches to solving the equations of non-radiative hydrodynamics: `classic Smoothed Particle Hydrodynamics (SPH); a novel SPH with a higher order dissipation sw
itch (SPHS); and adaptive mesh refinement (AMR). We find that SPHS and AMR are in excellent agreement, with both forming a well-defined entropy core that rapidly converges with increasing mass and force resolution. By contrast, SPH exhibits rather different behaviour. At low redshift, entropy decreases systematically with decreasing cluster-centric radius, converging on ever lower central values with increasing resolution. At higher redshift, SPH is in better agreement with SPHS and AMR but shows much poorer numerical convergence. We trace these discrepancies to artificial surface tension in SPH at phase boundaries. At early times, the passage of massive substructures close to the cluster centre stirs and shocks gas to build an entropy core. At later times, artificial surface tension causes low entropy gas to sink artificially to the centre of the cluster. We use SPHS to study the contribution of numerical versus physical dissipation on the entropy core, and argue that numerical dissipation is required to ensure single-valued fluid quantities in converging flows. However, provided this dissipation occurs only at the resolution limit, and provided that it does not propagate errors to larger scales, its effect is benign. There is no requirement to build `sub-grid models of unresolved turbulence for galaxy cluster simulations. We conclude that entropy cores in non-radiative simulations of galaxy clusters are physical, resulting from entropy generation in shocked gas during cluster assembly, putting to rest the long-standing puzzle of cluster entropy cores in AMR simulations versus their apparent absence in classic SPH simulations.
We discuss differences in simulation results that arise between the use of either the thermal energy or the entropy as an independent variable in smoothed particle hydrodynamics (SPH). In this context, we derive a new version of SPH that manifestly c
onserves both energy and entropy if smoothing lengths are allowed to adapt freely to the local mass resolution. To test various formu- lations of SPH, we consider point-like energy injection and find that powerful explosions are well represented by SPH even when the energy is deposited into a single particle, provided that the entropy equation is integrated. If the thermal energy is instead used as an independent variable, unphysical solutions can be obtained for this problem. We also examine the radiative cooling of gas spheres that collapse and virialize in isolation and of halos that form in cosmological simulations of structure formation. When applied to these problems, the thermal energy version of SPH leads to substantial overcooling in halos that are resolved with up to a few thousand particles, while the entropy formulation is biased only moderately low for these halos. For objects resolved with much larger particle numbers, the two approaches yield consistent results. We trace the origin of the differences to systematic resolution effects in the outer parts of cooling flows. The cumulative effect of this overcooling can be significant. In cosmological simulations of moderate size, we find that the fraction of baryons which cool and condense can be reduced by up to a factor ~2 if the entropy equation is employed rather than the thermal energy equation. We also demonstrate that the entropy method leads to a greatly reduced scatter in the density-temperature relation of the low-density Ly-alpha forest relative to the thermal energy approach, in accord with theoretical expectations.(abridged)
We study the effect of local stellar radiation and UVB on the physical properties of DLAs and LLSs at z=3 using cosmological SPH simulations. We post-process our simulations with the ART code for radiative transfer of local stellar radiation and UVB.
We find that the DLA and LLS cross sections are significantly reduced by the UVB, whereas the local stellar radiation does not affect them very much except in the low-mass halos. This is because clumpy high-density clouds near young star clusters effectively absorb most of the ionizing photons from young stars. We also find that the UVB model with a simple density threshold for self-shielding effect can reproduce the observed column density distribution function of DLAs and LLSs very well, and we validate this model by direct radiative transfer calculations of stellar radiation and UVB with high angular resolution. We show that, with a self-shielding treatment, the DLAs have an extended distribution around star-forming regions typically on ~ 10-30 kpc scales, and LLSs are surrounding DLAs on ~ 30-60 kpc scales. Our simulations suggest that the median properties of DLA host haloes are: Mh = 2.4*10^10 Msun, SFR = 0.3 Msun/yr, M* = 2.4*10^8 Msun, and Z/Zsun = 0.1. About 30 per cent of DLAs are hosted by haloes having SFR = 1 - 20 Msun/yr, which is the typical SFR range for LBGs. More than half of DLAs are hosted by the LBGs that are fainter than the current observational limit. Our results suggest that fractional contribution to LLSs from lower mass haloes is greater than for DLAs. Therefore the median values of LLS host haloes are somewhat lower with Mh = 9.6*10^9 Msun, SFR = 0.06 Msun/yr, M* = 6.5*10^7 Msun and Z/Zsun = 0.08. About 80 per cent of total LLS cross section are hosted by haloes with SFR < 1 Msun/yr, hence most LLSs are associated with low-mass halos with faint LBGs below the current detection limit.
We combine a Monte Carlo radiative transfer code with an SPH code, so that -- assuming thermal equilibrium -- we can calculate dust-temperature fields, spectral energy distributions, and isophotal maps, for the individual time-frames generated by an
SPH simulation. On large scales, the radiative transfer cells (RT cells) are borrowed from the tree structure built by the SPH code, and are chosen so that their size -- and hence the resolution of the calculated temperature field -- is comparable with the resolution of the density field. We refer collectively to these cubic RT cells as the global grid. The code is tested and found to treat externally illuminated dust configurations very well. However, when there are embedded discrete sources, i.e. stars, these produce very steep local temperature gradients which can only be modelled properly if -- in the immediate vicinity of, and centred on, each embedded star -- we supplement the global grid with a star grid of closely spaced concentric RT cells.
We present results of hydrodynamic simulations of star formation triggered by cloud-cloud collisions. During the early stages of star formation, low-mass objects form by gravitational instabilities in protostellar discs. A number of these low-mass ob
jects are in the sub-stellar mass range, including a few objects of planetary mass. The disc instabilities that lead to the formation of low-mass objects in our simulations are the product of disc-disc interactions and/or interactions between the discs and their surrounding gas.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
