اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدProtostellar collapse: A comparison between SPH and AMR calculations
40
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
The development of parallel supercomputers allows today the detailed study of the collapse and the fragmentation of prestellar cores with increasingly accurate numerical simulations. Thanks to the advances in sub-millimeter observations, a wide range of observed initial conditions enable us to study the different modes of low-mass star formation. The challenge for the simulations is to reproduce the observational results. Two main numerical methods, namely AMR and SPH, are widely used to simulate the collapse and the fragmentation of prestellar cores. We compare thoroughly these two methods within their standard framework. We use the AMR code RAMSES and the SPH code DRAGON. Our physical model is as simple as possible and consists of an isothermal sphere rotating around the z-axis. We first study the conservation of angular momentum as a function of the resolution. Then, we explore a wide range of simulation parameters to study the fragmentation of prestellar cores. There seems to be a convergence between the two methods, provided resolution in each case is sufficient. Resolution criteria adapted to our physical cases, in terms of resolution per Jeans mass, for an accurate description of the formation of protostellar cores are deduced from the present study. This convergence is encouraging for future work in simulations of low-mass star formation, providing the aforementioned criteria are fulfilled. Higher resolution figures can be downloaded at http://www-dapnia.cea.fr/Projets/COAST/paper_amrvssph.pdf
قيم البحث
اقرأ أيضاً
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 propose a new analytic model for the initial conditions of protostellar collapse in relatively isolated regions of star formation. The model is non-magnetic, and is based on a Plummer-like radial density profile as its initial condition. It fits:
the observed density profiles of pre-stellar cores and Class 0 protostars; recent observations in pre-stellar cores of roughly constant contraction velocities over a wide range of radii; and the lifetimes and accretion rates derived for Class 0 and Class I protostars. However, the model is very simple, having in effect only 2 free parameters, and so should provide a useful framework for interpreting observations of pre-stellar cores and protostars, and for calculations of radiation transport and time-dependent chemistry. As an example, we model the pre-stellar core L1544.
The core collapse supernova (CCSN) rate provides a strong lower limit for the star formation rate (SFR). Progress in using it as a cosmic SFR tracer requires some confidence that it is consistent with more conventional SFR diagnostics in the nearby U
niverse. This paper compares standard SFR measurements based on Halpha, FUV and TIR galaxy luminosities with the observed CCSN rate in the same galaxy sample. The comparison can be viewed from two perspectives. Firstly, by adopting an estimate of the minimum stellar mass to produce a CCSN one can determine a SFR from SN numbers. Secondly, the radiative SFRs can be assumed to be robust and then the SN statistics provide a constrain on the minimum stellar mass for CCSN progenitors. The novel aspect of this study is that Halpha, FUV and TIR luminosities are now available for a complete galaxy sample within the local 11Mpc volume and the number of discovered SNe in this sample within the last 13 years is high enough to perform a meaningful statistical comparison. We exploit the multi-wavelength dataset from 11HUGS, a volume-limited survey designed to provide a census of SFR in the local Volume. Assuming a lower limit for CCSN progenitors of 8 Msun, the CCSN rate matches the SFR from the FUV luminosity. However the SFR based on Halpha luminosity is lower than these two estimates by a factor of nearly 2. If we assume that the FUV or Halpha based luminosities are a true reflection of the SFR, we find that the minimum mass for CCSN progenitors is 8 +/- 1 Msun, and 6 +/- 1 Msun, respectively. The estimate of the minimum mass for CCSN progenitors obtained exploiting FUV data is in good agreement with that from the direct detection of CCSN progenitors. The concordant results by these independent methods point toward a constraint of 8 +/- 1 Msun on the lower mass limit for progenitor stars of CCSNe.
We investigate protostellar collapse of molecular cloud cores by numerical simulations, taking into account turbulence and magnetic fields. By using the adaptive mesh refinement technique, the collapse is followed over a wide dynamic range from the s
cale of a turbulent cloud core to that of the first core. The cloud core is lumpy in the low density region owing to the turbulence, while it has a smooth density distribution in the dense region produced by the collapse. The shape of the dense region depends mainly on the mass of the cloud core; a massive cloud core tends to be prolate while a less massive cloud core tends to be oblate. In both cases, anisotropy of the dense region increases during the isothermal collapse. The minor axis of the dense region is always oriented parallel to the local magnetic field. All the models eventually yield spherical first cores supported mainly by the thermal pressure. Most of turbulent cloud cores exhibit protostellar outflows around the first cores. These outflows are classified into two types, bipolar and spiral flows, according to the morphology of the associated magnetic field. Bipolar flow often appears in the less massive cloud core. The rotation axis of the first core is oriented parallel to the local magnetic field for bipolar flow, while the orientation of the rotation axis from the global magnetic field depends on the magnetic field strength. In spiral flow, the rotation axis is not aligned with the local magnetic field.
The formation and collapse of a protostar involves the simultaneous infall and outflow of material in the presence of magnetic fields, self-gravity, and rotation. We use self-similar techniques to self-consistently model the anisotropic collapse and
outflow by a set of angle-separated self-similar equations. The outflow is quite strong in our model, with the velocity increasing in proportion to radius, and material formally escaping to infinity in the finite time required for the central singularity to develop. Analytically tractable collapse models have been limited mainly to spherically symmetric collapse, with neither magnetic field nor rotation. Other analyses usually employ extensive numerical simulations, or either perturbative or quasistatic techniques. Our model is unique as an exact solution to the non-stationary equations of self-gravitating MHD, which features co-existing regions of infall and outflow. The velocity and magnetic topology of our model is quadrupolar, although dipolar solutions may also exist. We provide a qualitative model for the origin and subsequent evolution of such a state. However, a central singularity forms at late times, and we expect the late time behaviour to be dominated by the singularity rather than to depend on the details of its initial state. Our solution may, therefore, have the character of an attractor among a much more general class of self-similarity.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
