User:Tohline/SSC/UniformDensity

Uniform-Density Sphere (structure)

Here we're going to derive the interior structural properties of a uniform-density sphere using all three solution strategies. While deriving essentially the same solution three different ways might seem like a bit of overkill, this approach proves to be instructive because (a) it forces us to examine the structural behavior of a number of different physical parameters, and (b) it illustrates how to work through the different solution strategies for one model whose structure can in fact be derived analytically using any of the techniques. As we shall see when studying other self-gravitating configurations, the three strategies are not always equally fruitful.

Solution Technique 1

Adopting solution technique #1, we need to solve the integro-differential equation,

<math>~\frac{dP}{dr} = - \frac{GM_r \rho}{r^2}</math>

appreciating that,

<math> M_r \equiv \int_0^r 4\pi r^2 \rho dr </math> .

For a uniform-density configuration, <math>~\rho</math> = <math>\rho_c</math> = constant, so the density can be pulled outside the mass integral and the integral can be completed immediately to give,

<math> M_r = \frac{4\pi}{3}\rho_c r^3 </math> .

Hence, the differential equation describing hydrostatic balance becomes,

<math> \frac{dP}{dr} = - \frac{4\pi G}{3} \rho_c^2 r </math> .

Integrating this from the center of the configuration — where <math>r=0</math> and <math>P = P_c</math> — out to an arbitrary radius <math>r</math> that is still inside the configuration, we obtain,

<math> \int_{P_c}^P dP = - \frac{4\pi G}{3} \rho_c^2 \int_0^r r dr </math>
<math>\Rightarrow ~~~~ P = P_c - \frac{2\pi G}{3} \rho_c^2 r^2 </math>

We expect the pressure to drop to zero at the surface of our spherical configuration — that is, at <math>r=R</math> — so the central pressure must be,

<math>P_c = \frac{2\pi G}{3} \rho_c^2 R^2 = \frac{3G}{8\pi}\biggl( \frac{M^2}{R^4} \biggr)</math> ,

where <math>M</math> is the total mass of the configuration. Finally, then, we have,

<math>P(r) = P_c\biggl[1 - \biggl(\frac{r}{R}\biggr)^2 \biggr] </math> .

Solution Technique 3

Adopting solution technique #3, we need to solve the algebraic expression,

<math>H + \Phi = C_\mathrm{B}</math> .

in conjunction with the Poisson equation,

<math>\frac{1}{r^2} \frac{d }{dr} \biggl( r^2 \frac{d \Phi}{dr} \biggr) = 4\pi G \rho </math> .

Appreciating that, as shown above, for a uniform density (<math>~\rho</math> = <math>\rho_c</math> = constant) configuration,

<math> M_r = \int_0^r 4\pi r^2 \rho dr = \frac{4\pi}{3}\rho_c r^3 </math> ,

we can integrate the Poisson equation once to give,

<math> \frac{d\Phi}{dr} = \frac{4\pi G}{3} \rho_c r </math> ,

everywhere inside the configuration. Integrating this expression from any point inside the configuration to the surface, we find that,

<math> \int_{\Phi(r)}^{\Phi_\mathrm{surf}} d\Phi = \frac{4\pi G}{3} \rho_c \int_r^R r dr </math>
<math>\Rightarrow ~~~~~ \Phi_\mathrm{surf} - \Phi(r) = \frac{2\pi G}{3} \rho_c R^2 \biggl[ 1- \biggl(\frac{r}{R} \biggr)^2 \biggr] </math>

Turning to the above algebraic condition, we will adopt the convention that <math>~H</math> is set to zero at the surface of a barotropic configuration, in which case the constant, <math>C_\mathrm{B}</math>, must be,

<math>C_\mathrm{B} = (H + \Phi)_\mathrm{surf} = \Phi_\mathrm{surf}</math> .

Therefore, everywhere inside the configuration <math>~H</math> must be given by the expression,

<math>H(r) = \Phi_\mathrm{surf} - \Phi(r)</math> .

Matching this with our solution of the Poisson equation, we conclude that, throughout the configuration,

<math> H(r) = \frac{2\pi G}{3} \rho_c R^2 \biggl[ 1- \biggl(\frac{r}{R} \biggr)^2 \biggr]</math> .

Comparing this result with the result we obtained using solution technique #1, it is clear that throughout a uniform-density, self-gravitating sphere,

<math>\frac{P}{H} = \rho</math> .

Solution Technique 2

Adopting solution technique #2, we need to solve the following single, <math>2^\mathrm{nd}</math>-order ODE:

<math>\frac{1}{r^2} \frac{d }{dr} \biggl( r^2 \frac{d H}{dr} \biggr) = - 4\pi G \rho </math> .

Appreciating again that, for a uniform density (<math>~\rho</math> = <math>\rho_c</math> = constant) configuration,

<math> M_r = \int_0^r 4\pi r^2 \rho dr = \frac{4\pi}{3}\rho_c r^3 </math> ,

we can integrate the <math>2^\mathrm{nd}</math>-order ODE once to give,

<math> \frac{dH}{dr} = -\frac{4\pi G}{3} \rho_c r </math> ,

everywhere inside the configuration. Integrating this expression from any point inside the configuration to the surface — where, again, we adopt the convention that <math>~H</math> = 0 — we find that,

<math> \int_{H(r)}^{0} dH = - \frac{4\pi G}{3} \rho_c \int_r^R r dr </math>
<math>\Rightarrow ~~~~~ H(r) = \frac{2\pi G}{3} \rho_c R^2 \biggl[ 1- \biggl(\frac{r}{R} \biggr)^2 \biggr] </math> .

Summary

From the above derivations, we can describe the properties of a uniform-density, self-gravitating sphere as follows:

• Mass:

Given the density, <math>\rho_c</math>, and the radius, <math>R</math>, of the configuration, the total mass is,

<math>M = \frac{4\pi}{3} \rho_c R^3 </math> ;

and, expressed as a function of <math>M</math>, the mass that lies interior to radius <math>r</math> is,

<math>\frac{M_r}{M} = \biggl(\frac{r}{R} \biggr)^3</math> .

• Pressure:

Given values for the pair of model parameters <math>( \rho_c , R )</math>, or <math>( M , R )</math>, or <math>( \rho_c , M )</math>, the central pressure of the configuration is,

<math>P_c = \frac{2\pi G}{3} \rho_c^2 R^2 = \frac{3G}{8\pi}\biggl( \frac{M^2}{R^4} \biggr) = \biggl[ \frac{\pi}{6} G^3 \rho_c^4 M^2 \biggr]^{1/3}</math> ;

and, expressed in terms of the central pressure <math>P_c</math>, the variation with radius of the pressure is,

<math>P(r) = P_c \biggl[ 1 -\biggl(\frac{r}{R} \biggr)^2 \biggr]</math> .

• Enthalpy:

Throughout the configuration, the enthalpy is given by the relation,

<math>H(r) = \frac{P(r)}{ \rho_c} = \frac{GM}{2R} \biggl[ 1 -\biggl(\frac{r}{R} \biggr)^2 \biggr]</math> .

• Gravitational potential:

Throughout the configuration — that is, for all <math>r \leq R</math> — the gravitational potential is given by the relation,

<math>\Phi_\mathrm{surf} - \Phi(r) = H(r) = \frac{G M}{2R} \biggl[ 1- \biggl(\frac{r}{R} \biggr)^2 \biggr] </math> .

Outside of this spherical configuration— that is, for all <math>r \geq R</math> — the potential should behave like a point mass potential, that is,

<math>\Phi(r) = - \frac{GM}{r} </math> .

Matching these two expressions at the surface of the configuration, that is, setting <math>\Phi_\mathrm{surf} = - GM/R</math>, we have what is generally considered the properly normalized prescription for the gravitational potential inside a uniform-density, spherically symmetric configuration:

<math>\Phi(r) = - \frac{G M}{R} \biggl\{ 1 + \frac{1}{2}\biggl[ 1- \biggl(\frac{r}{R} \biggr)^2 \biggr] \biggr\} = - \frac{3G M}{2R} \biggl[ 1 - \frac{1}{3} \biggl(\frac{r}{R} \biggr)^2 \biggr] </math> .

• Mass-Radius relationship:

We see that, for a given value of <math>\rho_c</math>, the relationship between the configuration's total mass and radius is,

<math>M \propto R^3 ~~~~~\mathrm{or}~~~~~R \propto M^{1/3} </math> .

• Central- to Mean-Density Ratio:

Because this is a uniform-density structure, the ratio of its central density to its mean density is unity, that is,

<math>\frac{\rho_c}{\bar{\rho}} = 1 </math> .

© 2014 - 2021 by Joel E. Tohline |   H_Book Home   |   YouTube   | Appendices: | Equations | Variables | References | Ramblings | Images | myphys.lsu | ADS | Recommended citation:   Tohline, Joel E. (2021), The Structure, Stability, & Dynamics of Self-Gravitating Fluids, a (MediaWiki-based) Vistrails.org publication, https://www.vistrails.org/index.php/User:Tohline/citation