Graphing Densities

Data was taken from the "Standard Solar Model BS05(AGS, OP): this included solar density, solar pressure, mass and pressure for 800 points along a solar radius.

A code was then drafted which placed a 'bubble' of gas at some point on the radius, and the density was calculated via the Ideal Gas Law. This was done incorrectly for a few days, with the density of the hotter cell being calculated as much higher than the cooler surrounding sun. This occurred, as the data was not in Pascals or kg/m^3, but rather Dynes and g/cm^3. A factor of 0.6 was also added in, which represented the mean mass of the Sun.

The code used (bubble.py) asks the user for an initial mass in Mkg, inital radius and initial temperature in MK. From this, it calculates the density of the cell at every point along the radius until 0.95 solar radii. It shows the 'measured' density of the Sun above this on the same graph. It also calculates the volume of the cell at the initial point along the radius.

The output looks like this:

In [2]: % run -i bubble.py

Mass (Mkg): 500

Radial Fraction: 0.74

Temperature (MK): 2.7

For a cell of mass 500E6 kg, the volume is 4653384.02 cubic metres

http://www.sns.ias.edu/~jnb/SNdata/sndata.html#bs2005 - Data used.