Changes between Version 26 and Version 27 of u/EricasLibrary


Ignore:
Timestamp:
07/02/12 15:55:01 (12 years ago)
Author:
Erica Kaminski
Comment:

—

Legend:

Unmodified
Added
Removed
Modified

  • u/EricasLibrary

    v26 v27  
    1 [[CollapsibleStart(Banerjee & Pudritz, et al. (‘04))]]
    2 
    3 
    4 [http://adsabs.harvard.edu/abs/2004MNRAS.355..248B The formation and evolution of protostellar discs; three-dimensional adaptive mesh refinement hydrosimulations of collapsing, rotating Bonnor-Ebert spheres]
    5 
    6 === Summary ===
    7 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.
    8 
    9 === Methods ===
    10 
    11 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.
    12 
    13 === Results ===
    14 
    15 No fragmentation/non-central collapse in the isothermal regime.
    16 
    17  Bars or rings form after the 'non-isothermal core' accretes surrounding gas and forms a disc.
    18 
    19 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.
    20 
    21 
    22 
    23 
    24 
    25 
    26 
    27 [[CollapsibleEnd]]
    28 
    29 
    30 
    31 [[CollapsibleStart(Fedderath, et al. (‘10))]]
    32 
    33 [http://adsabs.harvard.edu/abs/2010ApJ...713..269F Modeling Collapse and Accretion in Turbulent Gas Clouds: Implementation and Comparison of Sink Particles in AMR and SPH]
    34 
    35 === Background ===
    36 
    37 The numerical difficulty with modeling the collapse of a clump, while keeping track of the entire cloud, is given by the fact that the free-fall time, Tff, where
    38 
    39 {{{#!latex
    40 $ Tff = (3\pi / 32G\rho)^{1/2} $
    41 }}}
    42 
    43 decreases with increasing density and so resolution of the subgrids is demanded over many dynamical time scales. The two methods thus far designed to deal with this matter are 'Jeans heating' and 'sink particles'. Sink particles are a more realistic methodology and was first developed for SPH by Bate et al. ('95). The algorithm was later adapted by Krumholz et al. ('04) for Eulerian AMR. This paper describes a more rigorous series of checks for sink formation.
    44 
    45 === Sink Implentation ===
    46 
    47 Sink particles enable the star formation rate/star formation efficiency, and mass distribution, to be addressed in a robust and quantitative way.
    48 
    49 Sink algorithms originated from the notion of a density threshold. In earlier work, such a threshold was defined, and once surpassed, a sink particle was placed in the grid.
    50 
    51 The present work, however, has added in addition to this criterion, a series of checks to insure that sinks are formed only in gravitationally bound and collapsing structures. These are listed as follows:
    52 
    53 * Converging flow
    54 * Bound system
    55 * Jean's unstable
    56 * Gravitational potential minima
    57 
    58 Else, under conditions such as shear, a sink may erroneously form.
    59 
    60 === Tests ===
    61 Tests included collapse of Bonner Ebert and singular isothermal spheres.
    62 
    63 [[CollapsibleEnd]]
    64 
    65 [[CollapsibleStart(Foster and Chevalier (‘93))]]
    66 
    67 [http://articles.adsabs.harvard.edu/cgi-bin/nph-iarticle_query?1993ApJ...416..303F&data_type=PDF_HIGH&whole_paper=YES&type=PRINTER&filetype=.pdf Gravitational Collapse of an Isothermal Sphere]
    68 
    69 === Motivation ===
    70 
    71 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.
    72 
    73 === Summary ===
    74 
    75 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.
    76 
    77 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.
    78 
    79 === Stability Tests ===
    80 
    81 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:
    82 
    83 * xi = 3, flow oscillated about v=0 for ~ 15 free fall times
    84 * xi = 6.451, collapse happened in 5 ff times
    85 * xi = 20, collapse happened in 6 ff times
    86 
    87 === Collapse Results ===
    88 
    89 - 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.
    90 
    91 - Simulation approaches LP solution in center region.
    92 
    93 - Velocity smoothly transitions from homologous inflow to an inflow rate of ~ -3Cs after core formation.
    94 
    95 - Inner region evolves fast, outer evolves slow.
    96 
    97 - Core formation is not a subsonic adjustment of density with radius as Shu said. The inner regions collapse to a core supersonically.
    98 
    99 - 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.
    100 
    101 - Should compare anything to mass accretion rate?
    102 
    103 
    104 
    105 [[CollapsibleEnd]]
    106 
    107 [[CollapsibleStart(Truelove, Klein, Mc Kee, et al. (‘98))]]
    108 
    109 [http://iopscience.iop.org/0004-637X/495/2/821/pdf/36882.pdf Self-Gravitational hydrodynamics with 3d AMR: methodolgy and applications to molecular cloud collapse and fragmentation.]
    110 
    111 
    112 === Summary ===
    113 Develops methods for AMR, dynamic grids that allot finer resolution over many length scales that is imperative for studying problems of gravitational collapse. The criteria for refinement is crucial, and it is the jeans condition. This sets the resolution smaller than the local Jean’s condition, allowing new benchmarks in the probing of the dynamics. They find uniformly rotating spherical clouds to collapse along the equatorial plane. When perturbed, these form ‘filamentary singularties’ that don’t fragment when isothermal.
    114 
    115 === Methodology ===
    116 This paper goes through in extensive detail on the 3 components of their code methodology: the hyperbolic solvers that employ the Gudonov method for solution of the hydro equations,  elliptic solvers that utilize AMR multigrid method to solve Poisson’s equation, and finally these two methods operation within an AMR framework. The use of stencils as different layers of cell-centered quantities used for averaging the node centered quantities (for self - gravity) is detailed. Paper refers to Almgren for discussion on AMR multigrid cycle procedures.
    117 
    118 === AMR and refinement criteria ===
    119 
    120 Paper refers to Pember et al. (‘98) for a summary of the procedure for advancing grids. Their methods are based on Berger and Oliger (‘84), Berger and Collela (‘89), and Bell et al. (‘94).
    121 
    122 To trigger refinement from level 0 to level 1, the code employs a density criterion to cells that contain gas of the original cloud. Such cells are identified if they contain a density greater than or equal to 1/2 the original cloud density at the outer edge. This results in the level 0 cells effectively removing the boundary of the computational volume from the surface of the cloud. Jean's condition is than used for refinement to 2nd and higher levels (discussed below).
    123 
    124 
    125 === The Truelove criteria ===
    126 
    127 Jeans' analysis of the linearized 1D gravitohydrodynamic (GHD) equations of a medium of infinite extent lead to the expression for the jeans length:
    128 
    129 {{{#!latex
    130 $ \lambda_j  = (\frac{\pi c_s^2}{G \rho})^{1/2}$
    131 }}}
    132 
    133 which showed that perturbations that are larger than this are gravitationally unstable and collapse. Truelove ('97) showed that these instabilities can be triggered from numerical errors of the GHD solvers. Errors introduced on the cell scale at coarser grids can be transmitted to finer levels and can lead to 'artificial fragmentation'. A way to avoid this is to maintain resolution of the local jeans length. By defining:
    134 
    135 {{{#!latex
    136 $ J =  \frac{\triangle x}{\lambda_j} $
    137 }}}
    138 
    139 Truelove ('97) showed that keeping
    140 
    141 {{{#!latex
    142 $ J \leqq 0.25$ 
    143 }}}
    144 
    145 prevented artificial fragmentation of an isothermal cloud spanning 7 decades of density. It is expected to be necessary albeit not necessarily sufficient for isothermal collapse. By not having the resolution of gradients can trigger artificial viscosity (used for numerical stability), which can lead to an incorrect formulation of the problem originally deemed inviscid.
    146  
    147 
    148 
    149 === Poisson's boundary conditions ===
    150 
    151 The boundary condition they employ for self-gravity is periodic. They describe this as mimicking an infinite region of repeating boxes. Thus, when they model a cloud in the box using these boundary conditions, the solver for the gravity within this box, takes into account the gravitational effects of an infinite number of other spheres in boxes. To minimize the effect of these other spheres, they suggest making the distance between the center-to-center images d=4R. By increasing the d=8R, results in less gravitational pull by the mirrored clouds that hastens the collapse by only 1%.
    152 
    153 [[CollapsibleEnd]]
    154 
    155 
    156