Changes between Version 23 and Version 24 of u/massBElight


Ignore:
Timestamp:
10/05/12 21:09:56 (12 years ago)
Author:
Erica Kaminski
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • u/massBElight

    v23 v24  
    4343
    4444
    45 I see that Mbe crit is below the minima of M(t). This is no doubt due to some SMALL round-off error somewhere in computing the total masses in the grid. If for instance, I used simply the formula above to find the mass of the BE sphere, it does not give me exactly what my code gives. The Visit calculation gives me ~0.78 when I subtract the mass of the ambient from the mass of the sphere, whereas my code says it is ~0.72 -- in line with the Crit Mbe as expected.
     45I see that Mbe crit is below the minima of M(t). This is no doubt due to some SMALL round-off error somewhere in computing the total masses in the grid. If for instance, I used simply the formula above to find the mass of the BE sphere, it does not give me exactly what my code gives. The Visit calculation gives me ~0.78 when I subtract the mass of the ambient from the mass of the sphere, whereas my code says it is ~0.72 -- in line with the Crit Mbe as expected. This seems to indicate that these curves may be closer to each other than the are here. It would be interesting to see where these curves intersected. At least for now, it is safe to take away from this plot that M(t) does not fluctuate wildly, and never goes more than 100% over the crit Mbe..?
    4646
    47 Now, it should also, be noted that the crit Mbe is computed by the formula:
     47Now, it should also be noted that the crit Mbe is computed by the formula:
    4848
    4949
     
    5252}}}
    5353
    54 This describes __ blah blah, whereas my sphere is _blah Blah.
    55 
    56 
    57 I guess what can be taken home safely from this curve is that M(t) does not fluctuate wildly, and never goes more than 10% over the crit Mbe.       
     54where Cs = isothermal sound speed, Po = external pressure, and G is gravitational constant. This Mbe is a measure of the maximum mass a BE sphere can hold up against its own gravity. For the values used to compute the crit Mbe, I used Cs for the BE sphere (r=Rbe), whereas the Cs would actually increase 10-fold outside of the sphere, since the boundary at Rbe is held at pressure equilibrium.