Version 65 (modified by 10 years ago) ( diff ) | ,
---|
Numerics
AMR
-Self-Gravitational hydrodynamics with 3d AMR: methodolgy and applications to molecular cloud collapse and fragmentation. Truelove, Klein, Mc Kee, et al. (‘98)
Sink Particles
-Modeling Collapse and Accretion in Turbulent Gas Clouds: Implementation and Comparison of Sink Particles in AMR and SPHFFedderath, et al. (‘10)
Self Similarity Solutions for Collapse
-Self-similar collapse of isothermal spheres and star formation Shu (‘77)
-The collapse of unstable isothermal spheres Hunter, (‘77)
Larson & Penston et al. (‘04)
[http://adsabs.harvard.edu/abs/1969MNRAS.144..425P Dynamics of self-gravitating gaseous spheres-III. Analytical results in the free-fall of isothermal case] [
Bonnor Ebert Spheres
Banerjee & Pudritz et al. (‘04)
Summary
3d hydro, AMR simulations of rotating Bonnor Ebert spheres with varying degrees of initial rotation. Final collapse structures depend on the degree of rotation, with the general trend being slower rotation leads to bars/filaments, faster rotation leads to rings. Radiative cooling is incorporated into the model and a more complicated EOS.
Methods
Both high and low mass BE spheres were studied. All runs had an m=2 density perturbation. The EOS was more sophisticated than earlier works in that it was a function of density and time.
Results
No fragmentation/non-central collapse in the isothermal regime.
Bars or rings form after the 'non-isothermal core' accretes surrounding gas and forms a disc.
Temperature shocks form after the isothermal phase; as the core begins to hit up and halt the collapse, cold material from the envelope falls in and forms a shock front of temperature discontinuity. The inner shock front sets the conditions that control subsequent disc structure and dynamics.
Foster and Chevalier (‘93)
Gravitational Collapse of an Isothermal Sphere
Motivation
While the collapse problem of an isothermal sphere is well studied, there remains some conflict on the correct analytic solution for the collapse. Larson and Penston (1969) produced a solution for the collapse problem that was later criticized by Shu (1977) as being ad hoc in some of the mathematical formulations. Later numerical studies by Hunter (1977), however, showed that the collapse does follow the LP solution. Despite this, the Shu solution remains the prominent treatment of the problem. The authors here set out to check which solution indeed describes the collapse problem.
Summary
The Shu solution begins from core formation (born out of a singular isothermal sphere) and initiates collapse into the medium through an 'expansion wave'. In contrast, the LP solution follows the collapse of a uniformly dense, static sphere prior to and up to core formation. Hunter, however, found that the LP solution can in fact be extended past core formation. The Shu solution evolves to a solution similar to homologous inflow like the LP solution.
Foster and Chevalier ran simulations of a Bonner Ebert sphere with an outer density of 1/1000 the density at the outermost edge of the clump (my simulations used 1/100). The ambient therefore had a large sound speed throughout the constant pressure medium, that remained nearly constant throughout. Their results are robustly similar to the LP solution.
Stability Tests
The 'marginally stable sphere' had a xi=xiCrit=6.451. A density enhancement of 1% above equilibrium value was needed to bring out expected stability behavior:
- xi = 3, flow oscillated about v=0 for ~ 15 free fall times
- xi = 6.451, collapse happened in 5 ff times
- xi = 20, collapse happened in 6 ff times
Collapse Results
- Inflow is initially marked by velocity proportional to r for small r, peaking at xi = 2, and turning back to 0 velocity out toward edge of sphere.
- Simulation approaches LP solution in center region.
- Velocity smoothly transitions from homologous inflow to an inflow rate of ~ -3Cs after core formation.
- Inner region evolves fast, outer evolves slow.
- Core formation is not a subsonic adjustment of density with radius as Shu said. The inner regions collapse to a core supersonically.
- Varying the boundary conditions/initial conditions did not change the behavior from tending to the LP solution. Likewise, neither did changing the perturbation conditions. The mass accretion rate, however, if it became constant, was proportional to the magnitude of the perturbation as Shu found.
- Should compare anything to mass accretion rate?
Hennebelle Whitworth Gladwin Andre (‘2002)
Protostellar collapse by compression
The paper begins by providing some background information on different phases of the pre-stellar and protostellar phases. Definitions of class 0, 1, 2 are provided as well as their relative properties (flat inner core regions, envelope properties that differ in high density star forming regions compared to more sparse regions, infall velocity maps, etc.) and lifetimes (class 0 < class I < class II)- as inferred by statistical arguments. The observations seem to tell a coherent story of star formation that is at odds with the "standard theory of star formation" of Shu, Adams, and Lizano (1987), which describes collapse as occurring by way of the similarity solution of the SIS, predicting a prestellar cloud that is centrally peaked and static (vrad ~ 0 for all r), that collapses from the inside out in an expansion wave, and forms a central protostar that accretes surrounding material at a CONSTANT rate. However, these predictions of the SIS similarity solution are not supported by observation. While, the standard theory may be an okay 0th order approximation in quiescent, sparse regions of star formation such as Taurus, certainly more dynamic models in densely packed areas like rho-Ophiuchus are needed.
The paper then begins to set up its motivation - namely that research on flat topped isothermal spheres have begun to be pursued, such as F&C. They say F&C gets more accurate results, but that the collapse is triggered artificially. They plan to take the studies further on non-singular spheres, but by initiating collapse in a more physically relevant way (increasing Pext as suggested by Myers et al). They claim that most of the observational constraints (such as decreasing accretion rate, velocity fields, initial condition) are recovered by this model (although in my opinion - they also seem to have been recovered in F&C model too… so maybe they find MORE accurate lifetimes for their collapsing BE spheres). I find it interesting that they do not call the sphere a Bonnor Ebert sphere, but rather an isothermal sphere… is there any significance to this?
In their results section they present very weak quantitative results/explanations/insight of their plots (interestingly they spend very little time developing equations of BE sphere and also do not mention the critical values of the BE sphere such as Pcrit, Mcrit). Their different runs are broken down into the rate of compression, in how much time (measured in units of sound crossing time — ostensibly in terms of the BE sphere) does the Pext double. They initiate the collapse of the critical BE sphere by increasing Pext (note - any increase in Pext would exceed Pcrit as their initial sphere is already critical), but then they continue to increase Pext throughout the course of the simulation (this potentially is similar to what we did in practice by just allowing the simulation to evolve "naturally" with no forced increase in Pext, simply Pext at the sphere's outer edge increased naturally by the infall of material onto that outer edge). In the first couple of cases they show, compression is slow. In the first 3 panels of Fig. 1, a sink has not yet formed, and so they identify this phase as the pre-stellar (aka the pre-protostellar) phase. It appears to me that they get a nice outside in flow beginning to develop in these plots, although the radius of the BE sphere is unclear. No compression wave seems evident in these plots, but instead, it seems to be a re-equilibration of matter in terms of the language I use in my paper. This could be due to early subsonic adjustment of the material into a modified BE profile that exceeds the critical mass, and hence collapses in the canonical fashion. They use no such language, and make no such identification which is curious. They describe this, as well as all other cases, as being a compression wave solution. This might make sense if you follow Whitworth language that all collapse problems are due to compression waves, but some of the waves have 0 amplitude..? This case looks most like some of our lighter ambient runs, where the Pext increases very slowly due to a slow accumulation of matter.. In the next 2 panes of the plot, a central sink has formed that accretes matter in a free-fall manner (vrad~r½) that they say moves outward in time. I can't tell from this plot that this type of flow is moving outward in time. They talk about the mean cruising velocity (Whitworth language) for the different runs, but it is unclear if this is an average over the entire sphere at a given time, or if it is at a given radius over time, or what.. They say that the mean cruising velocity is set up by the compression wave, but how do they know when it penetrated the sphere if they can't even see it in the density plot? They do say that in this case the compression wave converges on the center, at which time the sink forms, so maybe in this limiting case of the Whitworth solution, the compression wave looks like the classic result?
Skipping to figure 3, for which the Pext doubles in a crossing time (aka strong compression). I have a hard time understanding their description in the text, p 874, right column, 2nd paragraph. They say a small compression wave moves in (I see this in plot), but they go on to say that before the wave converges on the center to form the sink (and class 0 phase - last 2 panes in figure), the central density has hardly changed "for the inner regions were unaware of Pext increasing". First of all, in the plot the central density is clearing increasing in panes 1-3. Second of all, they would be unaware only for a supersonic compression wave — for which they did not establish exists. They say again that around the sink, a free fall v field is set up, but in the outer regions a uniform sonic field is established (v~0.12 km/s). Again, I do not see a uniform v field in the outer regions.. I assume what they mean is uniform in r… not for a fixed r over time.. This plot does however, seem to show a change in the velocity field from the early "compression wave phase" to the later "classic phase", in the language of my paper. Although I think these plots overall do a poor job showing the evolution of the flow since they only sample 3 time states before sink formation, they don't capture the dynamics very well.
In Figure 4 for even stronger compression, the compression wave is now evident in the plots. I see in their velocity plots a turn-over similar to wee see for the different phases of evolution. The last plot they show before the sink forms, seems to be indicating a compression wave is about to "converge on the center". However, we also saw this type of evolution, but between this time and sink formation in our runs, the density profile evened out to a collapsing BE profile. I wonder if Hennebelle et al found similar results, and if so, why they did not choose to report it… This case seems most like our Matched case.
In Figure 5, the compression wave seems to completely rip through the BE sphere, although no mention of the speed of the incoming compression wave is given (which would have been interesting to compare to the triggered star formation stuff Shule has been looking at, which for mach > 20 has been shown to disrupt the BE sphere from being able to form a sink). It seems like the compression wave keeps building mass because of the abrupt increasing change in Pext that adds density to the wave in an isothermal fashion.
The results section concludes with the authors saying that when Pext is increased but then kept steady the results are different, citing figure 6. They do not ellaborate, so the reader is left to guess at what features they may be referring to. Looking at these plots, I am curious why the turn-up in vrad at large radii… We also saw a turn up in vrad in some cases…
The concluding remarks go like this — in each of the cases, key observational features are reproduced, namely the velocity fields. If collapse was triggered by a compression wave for the observational cases cited in the intro, then the Pext increased on such and such timescales. They go on to talk about the significance of the results in terms of compression wave solutions of Whitworth in the 3rd paragraph on page 881. They say, the generation of a uniform inward velocity field appears to be a general feature of compression wave solutions (again - I did not see this generated in their simulations from their plots), and that ALL OF THE INFINITY OF WHITWORTH SOLUTIONS involve a compression wave being driven from outside - in, leaving behind them vrad~constant and a 1/r2 density envelope. This needs to be clarified — is this vrad (R,t) = constant or is this vrad(r,T) = constant? Moreover, they say that the convergence of this wave on the center signifies the formation of the protostar - and that it accretes at a constant rate, developing an outward traveling expansion wave that leaves in its wake a v~r-½ and rho~r-3/2! He says Shu's solution is one limiting case of this behavior - of a 0 amplitude compression wave leaving behind a uniform v=0 velocity field. I am unsure though if the constant accretion rate is in fact a general feature of the Whitworth solution, as I have always heard that referred to in conjunction with the Shu solution. It is hard to connect the constant accretion rate to the Hennebelle results and others, given the accretion rates decline over time. Why they didn't clarify this, I am unsure…
The rest of the paper then just says an external increase in P is likely to happen from such and such astro phenomenon, and that the pressure increase is unlikely to be isotropic, so next steps should be adding some asymmetric forcing..
Molecular Clouds
Turbulence
Stellar Feedback
Stellar feedback by radiation pressure and photoionization
Sales, Marinacci, Springel, Petkova, 2013
Whereas radiation pressure has widely been regarded as the dominant source of feedback, this group suggests that photoionization may be a more efficient source of regulation, operating on shorter timescales than radiation pressure. Also read a little about their code, Arepo, which is a so-called 'moving mesh' code that has some adaptive refinement capabilities, uses a Godunov method + exact Riemann solver to calculate fluxes across the cells, and has self-gravity and radiation transfer in the 'single-scattering'/optically thin limit.
Attachments (4)
- Hennebelle1.png (385.0 KB ) - added by 12 years ago.
- Hennebell2.jpg (371.7 KB ) - added by 12 years ago.
- Hennebelle1.jpg (383.1 KB ) - added by 12 years ago.
- Intro2ISM.pdf (1.4 MB ) - added by 10 years ago.