Changes between Version 23 and Version 24 of u/massBElight
- Timestamp:
- 10/05/12 21:09:56 (12 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
u/massBElight
v23 v24 43 43 44 44 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. 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. 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..? 46 46 47 Now, it should also ,be noted that the crit Mbe is computed by the formula:47 Now, it should also be noted that the crit Mbe is computed by the formula: 48 48 49 49 … … 52 52 }}} 53 53 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. 54 where 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.