appletmagic.com Thayer Watkins Silicon Valley & Tornado Alley USA 


The purpose of this analysis is to derive the governing equations for three simple models of a star:
In the absence of rotation the body of the star is spherically symmetric so the only significant spatial variable is the distance r from the center. The gravitational attraction at r is balanced by the pressure gradient force; i.e.,
where G is the gravitational constant, M is the mass contained within a radius r, ρ is the mass density and p is pressure.
The mass contained within a radius r, M(r), is given by:
It is assumed that the star material obeys the ideal gas law; i.e.,
where R is the gas constant for the star material and T is the temperature in degrees Kelvin. For a steady state the temperature is uniform throughout the star. This means that
The governing equation for the internal structure of the star is:
The usual procedure for dealing with integraldifferential equations is to differentiate in order to eliminate the integration operation. But first a change of variable is in order. Let ln(ρ)=W so ρ=exp(W) and thus the governing equation is:
Now differentiation with respect to r gives:
Upon rearrangement this becomes:
In this model there is the creation of heat within the star matter. There is a balance between forces due to gravitation and the pressure gradient and the ideal gas law is assumed to apply as in Model I but the temperature is not uniform. Steady state conditions require that the heat created within the star be transferred to the surface and this requires a radial temperature gradient.
The heat energy passing through a surface is proportion to the area times the temperature gradient. The proportionality factor is negative because heat is transferred in the direction that temperature is declining. The net flow outflow from an infinitesimal volume is therefore proportional to the divergence of the temperature gradient. But for steady state conditions this net outflow must be equal to the heat generated within the infinitesimal volume. Thus
where c is the rate of production of heat per unit mass and D is the coefficient of heat conduction. As in Model I, ρ and T represent the mass density and temperature of the star material. This equation is of the form of a Poisson equation. The Laplacian of T, ∇^{2}T, for spherical coordinates when there is spherical symmetry is:
The full version of Model II is:
This last equation can be put into the form
which upon integration with respect to the radius variable gives
This means that
Since from the ideal gas equation
where γ = Gc/(4πDR).
This last equation above implies that
where ρ_{0} and T_{0} represents a standardized density and temperature.
Because from the ideal gas equation
The temperature profile is determined thus from the Poisson equation
When radiation is produced within the star radiation pressure as well as gas pressure support the gravitational attraction of the star matter. The radiation does not simply radiate away. The star material is opaque and a quantum of radiation is absorbed and reradiated randomly so the radiation must make its way to the surface by diffusion. A crucial parameter in this diffusion process is the average path length.
(To be continued.)
