|
|
| Line 38: |
Line 38: |
| ===Structure:=== | | ===Structure:=== |
|
| |
|
| Here we show how the set of principal governing equations (PGEs) can be solved to determine the equilibrium structure of spherically symmetric fluid configurations — such as individual, nonrotating stars or protostellar gas clouds. After supplementing the PGEs by specifying an equation of state of the fluid, the system of equations is usually solved by employing one of three techniques to obtain a "detailed force-balanced" model that provides the radius, <math>~R_\mathrm{eq}</math>, of the equilibrium configuration — given its mass, <math>~M</math>, and central pressure, <math>~P_c</math>, for example — as well as details regarding the internal radial profiles of the fluid density and fluid pressure. As our various discussions illustrate, simply varying the power-law index, <math>~n</math>, in a ''polytropic'' equation of state, | | Here we show how the set of principal governing equations (PGEs) can be solved to determine the equilibrium structure of spherically symmetric fluid configurations — such as individual, nonrotating stars or protostellar gas clouds. After supplementing the PGEs by specifying an equation of state of the fluid, the system of equations is usually solved by employing one of three techniques to obtain a "detailed force-balanced" model that provides the radius, <math>~R_\mathrm{eq}</math>, of the equilibrium configuration — given its mass, <math>~M</math>, and central pressure, <math>~P_c</math>, for example — as well as details regarding the internal radial profiles of the fluid density and fluid pressure. As our various discussions illustrate, simply varying the power-law index in a ''polytropic'' equation of state gives rise to equilibrium configurations that have a wide variety of internal structural profiles. |
| <table border="0" align="right" cellpadding="5">
| |
| <tr>
| |
| <td align="center" bgcolor="white">
| |
| [[File:FreeNRGpressureRadiusIsothermal.png|250px|right|border|Whitworth's (1981) Isothermal Free-Energy Surface]]
| |
| </td>
| |
| </tr>
| |
| <tr>
| |
| <td align="center">
| |
| [[User:Tohline/SSC/Virial/PolytropesEmbeddedOutline#3DIsothermalSurface|Free-Energy Surface]]
| |
| </td>
| |
| </tr>
| |
| </table>
| |
| <div align="center">
| |
| <table border="0" cellpadding="5" align="center">
| |
|
| |
|
| <tr>
| | If one is not particularly concerned about details regarding the distribution of matter ''within'' the equilibrium configuration, a good estimate of the size of the equilibrium system can be determined by assuming a uniform-density structure then identifying local extrema in the system's global free energy … [[User:Tohline/SphericallySymmetricConfigurations/IntroductorySummary#Applications|<more>]] |
| <td align="right">
| |
| <math>~P</math>
| |
| </td>
| |
| <td align="center">
| |
| <math>~=</math>
| |
| </td>
| |
| <td align="left">
| |
| <math>~K \rho^{(n+1)/n} \, ,</math>
| |
| </td>
| |
| </tr>
| |
| </table>
| |
| </div> | |
| gives rise to equilibrium configurations that have a wide variety of internal structural profiles.
| |
|
| |
|
| If one is not particularly concerned about details regarding the distribution of matter ''within'' the equilibrium configuration, a good estimate of the size of the equilibrium system can be determined by assuming a uniform-density structure then identifying local extrema in the system's global free energy. An illustrative, undulating free-energy surface is displayed here, on the right; blue dots identify equilibria associated with a "valley" of the free-energy surface while white dots identify equilibria that lie along a "ridge" in the free-energy surface.
| |
|
| |
| In the astrophysics community, the mathematical relation that serves to define the properties of configurations that are associated with such free-energy extrema is often referred to as the ''scalar virial theorem''. Specifically, for ''isolated'' systems in virial equilibrium, the following relation between configuration parameters holds:
| |
| <div align="center">
| |
| <table border="0" cellpadding="5" align="center">
| |
|
| |
| <tr>
| |
| <td align="right">
| |
| <math>~\frac{GM^2}{P_c R_\mathrm{eq}^4}</math>
| |
| </td>
| |
| <td align="center">
| |
| <math>~=</math>
| |
| </td>
| |
| <td align="left">
| |
| <math>~\biggl(\frac{2^2\cdot 5 \pi}{3} \biggr) \frac{\mathfrak{f}_A \cdot \mathfrak{f}_M^2}{\mathfrak{f}_W}
| |
| \, ,</math>
| |
| </td>
| |
| </tr>
| |
| </table>
| |
| </div>
| |
| where all three of the dimensionless ''structural form factors'', <math>~\mathfrak{f}_M</math>, <math>~\mathfrak{f}_W</math>, and <math>~\mathfrak{f}_A</math>, are unity, under the assumption that the equilibrium configuration has uniform density and uniform pressure throughout, and are otherwise generically ''of order unity'' for detailed force-balanced models having a wide range of internal structures. Alternatively, if the specific entropy of fluid elements (set by the value of <math>~K</math>) throughout the system, rather than the central pressure, is held fixed while searching for extrema in the free-energy, the virial equilibrium relation for isolated polytropes is,
| |
| <div align="center">
| |
| <table border="0" cellpadding="5" align="center">
| |
|
| |
| <tr>
| |
| <td align="right">
| |
| <math>~M^{(n-1)/n} R_\mathrm{eq}^{(3-n)/n} \biggl( \frac{G}{K} \biggr) </math>
| |
| </td>
| |
| <td align="center">
| |
| <math>~=</math>
| |
| </td>
| |
| <td align="left">
| |
| <math>~\frac{5\mathfrak{f}_A \mathfrak{f}_M}{\mathfrak{f}_W} \biggl(\frac{3}{4\pi \mathfrak{f}_M}\biggr)^{1/n} \, .
| |
| </math>
| |
| </td>
| |
| </tr>
| |
| </table>
| |
| </div>
| |
|
| |
| If the physical system under consideration — such as a protostellar gas cloud — is not isolated but is, instead, embedded in a hot, tenuous medium that exerts on the system a confining external pressure, <math>~P_e</math>, the configuration's equilibrium parameters will be related via the expression,
| |
|
| |
| <div align="center">
| |
| <table border="0" cellpadding="5" align="center">
| |
|
| |
| <tr>
| |
| <td align="right">
| |
| <math>~\frac{GM^2}{P_c R_\mathrm{eq}^4}</math>
| |
| </td>
| |
| <td align="center">
| |
| <math>~=</math>
| |
| </td>
| |
| <td align="left">
| |
| <math>~\biggl(\frac{2^2\cdot 5 \pi}{3} \biggr) \frac{\mathfrak{f}_M^2}{\mathfrak{f}_W} \biggl[ \mathfrak{f}_A - \frac{P_e}{P_c} \biggr]
| |
| \, ;</math>
| |
| </td>
| |
| </tr>
| |
| </table>
| |
| </div>
| |
| or, fixing <math>~K</math> instead of <math>~P_c</math>, the relevant virial equilibrium expression is,
| |
| <div align="center">
| |
| <table border="0" cellpadding="5" align="center">
| |
|
| |
| <tr>
| |
| <td align="right">
| |
| <math>~P_e </math>
| |
| </td>
| |
| <td align="center">
| |
| <math>~=</math>
| |
| </td>
| |
| <td align="left">
| |
| <math>~K \mathfrak{f}_A \biggl( \frac{3M}{4\pi R_\mathrm{eq}^3} \cdot \frac{1}{\mathfrak{f}_M} \biggr)^{1 + 1/n}
| |
| - \frac{\mathfrak{f}_W}{5} \biggl(\frac{3GM^2}{4\pi R_\mathrm{eq}^4} \cdot \frac{1}{\mathfrak{f}_M^2} \biggr) \, .</math>
| |
| </td>
| |
| </tr>
| |
| </table>
| |
| </div>
| |
| It is a virial expression specifically of this form <math>~(</math>with <math>n = \infty</math> and <math>\mathfrak{f}_M = \mathfrak{f}_W = \mathfrak{f}_A = 1)</math> that identifies extrema (e.g., valleys or ridges) in the rainbow-colored free-energy surface, <math>\mathfrak{G}^*(R_\mathrm{eq}, P_e)</math>, displayed above. As can be determined from this algebraic expression and as the figure illustrates, for any specified mass no equilibrium states exist if <math>~P_e</math> is greater than some limiting value, <math>~P_\mathrm{crit}</math>; the equilibrium configuration associated with the limiting condition, <math>~P_e = P_\mathrm{crit}</math>, is marked by a red dot on the displayed free-energy surface. The astrophysical significance of this critical state was first discussed in the mid 1950s in the context of star formation and, specifically, [[User:Tohline/SSC/Structure/BonnorEbert#Pressure-Bounded_Isothermal_Sphere|Bonnor-Ebert spheres]].
| |
|
| |
| After rearranging terms, for any specified values of the parameters <math>~P_e</math> and <math>~K</math>, this virial equilibrium expression can also be viewed as a mass-radius relation of the form,
| |
| <table border="0" align="right" cellpadding="5">
| |
| <tr>
| |
| <td align="center" bgcolor="white">
| |
| [[File:MassRadiusVirialLabeled.png|200px|right|border|Virial Mass-Radius Relation]]
| |
| </td>
| |
| </tr>
| |
| </table>
| |
| <div align="center">
| |
| <table border="0" cellpadding="5" align="center">
| |
|
| |
| <tr>
| |
| <td align="right">
| |
| <math>~a R_\mathrm{eq}^4 - b M^{(n+1)/n} R_\mathrm{eq}^{(n-3)/n} + c M^2</math>
| |
| </td>
| |
| <td align="center">
| |
| <math>~=</math>
| |
| </td>
| |
| <td align="left">
| |
| <math>~ 0\, ,</math>
| |
| </td>
| |
| </tr>
| |
| </table>
| |
| </div>
| |
| where,
| |
| <div align="center">
| |
| <table border="0" cellpadding="5" align="center">
| |
|
| |
| <tr>
| |
| <td align="right">
| |
| <math>~a</math>
| |
| </td>
| |
| <td align="center">
| |
| <math>~\equiv</math>
| |
| </td>
| |
| <td align="left">
| |
| <math>~ \frac{4\pi}{3} \cdot P_e \, ,</math>
| |
| </td>
| |
| </tr>
| |
|
| |
| <tr>
| |
| <td align="right">
| |
| <math>~b</math>
| |
| </td>
| |
| <td align="center">
| |
| <math>~\equiv</math>
| |
| </td>
| |
| <td align="left">
| |
| <math>~ \biggl( \frac{3}{4\pi} \biggr)^{1/n} \cdot K \mathfrak{f}_A \mathfrak{f}_M^{-(n+1)/n} \, ,</math>
| |
| </td>
| |
| </tr>
| |
|
| |
| <tr>
| |
| <td align="right">
| |
| <math>~c</math>
| |
| </td>
| |
| <td align="center">
| |
| <math>~\equiv</math>
| |
| </td>
| |
| <td align="left">
| |
| <math>~ \frac{G\mathfrak{f}_W}{5\mathfrak{f}_M^2} \, .</math>
| |
| </td>
| |
| </tr>
| |
| </table>
| |
| </div>
| |
| Using this virial equilibrium relation (and assuming a = b = c = 1), the curves drawn in the figure that is displayed here, on the right, show how the equilibrium radius of an embedded, pressure-truncated polytropic sphere varies with mass for seven different adopted polytropic indexes. In direct analogy with the critical pressure that is associated with Bonnor-Ebert spheres, for systems having <math>~n \ge 3</math>, there is a mass, <math>~M_\mathrm{max}</math>, above which equilibrium configurations do not exist; and, when <math>~n > 3</math>, two equilibrium configurations having different radii can be constructed for any system having a mass, <math>~M < M_\mathrm{max}</math>.
| |
|
| |
|
| {{LSU_HBook_SSC_Structure}} | | {{LSU_HBook_SSC_Structure}} |
Preface from the original version of this HyperText Book (H_Book):
November 18, 1994
Much of our present, basic understanding of the structure, stability, and dynamical evolution of individual stars, short-period binary star systems, and the gaseous disks that are associated with numerous types of stellar systems (including galaxies) is derived from an examination of the behavior of a specific set of coupled, partial differential equations. These equations — most of which also are heavily utilized in studies of continuum flows in terrestrial environments — are thought to govern the underlying physics of all macroscopic "fluid" systems in astronomy. Although relatively simple in form, they prove to be very rich in nature... <more>
Pictorial Table of Contents
Context
Applications
Spherically Symmetric Configurations
Introduction
(Alternate Introduction)
Structure:
Here we show how the set of principal governing equations (PGEs) can be solved to determine the equilibrium structure of spherically symmetric fluid configurations — such as individual, nonrotating stars or protostellar gas clouds. After supplementing the PGEs by specifying an equation of state of the fluid, the system of equations is usually solved by employing one of three techniques to obtain a "detailed force-balanced" model that provides the radius, <math>~R_\mathrm{eq}</math>, of the equilibrium configuration — given its mass, <math>~M</math>, and central pressure, <math>~P_c</math>, for example — as well as details regarding the internal radial profiles of the fluid density and fluid pressure. As our various discussions illustrate, simply varying the power-law index in a polytropic equation of state gives rise to equilibrium configurations that have a wide variety of internal structural profiles.
If one is not particularly concerned about details regarding the distribution of matter within the equilibrium configuration, a good estimate of the size of the equilibrium system can be determined by assuming a uniform-density structure then identifying local extrema in the system's global free energy … <more>
|
Solution Strategies:
|
Detailed Force-Balance
(Introduction)
|
Virial Equilibrium
(Introduction)
|
| Example Solutions: |
|
|
| | Isolated … |
✓ | ✓ |
| | Embedded in an External Medium … |
✓ | ✓ |
|
|
| | Isolated … |
✓ | ✓ |
| | Embedded in an External Medium … |
✓ | ✓ 🎦 |
|
|
| | Isolated … |
✓ | ✓ |
| | Embedded in an External Medium (Bonnor-Ebert Sphere) … |
✓ | ✓ |
- Zero-temperature White Dwarf — Overview
- Power-law density distribution — Overview
|
- Other Analytically Definable Models
|
✓ |
… |
- BiPolytropes (also referred to as Composite Polytropes)
|
Overview |
Overview |
| | Core-Envelope Structure with <math>~(n_c,n_e) = (0,0)</math> … |
✓ |
✓ |
| | Core-Envelope Structure with <math>~(n_c,n_e) = (5,1)</math> … |
✓ |
✓ |
| | Core-Envelope Structure with <math>~(n_c,n_e) = (1,5)</math> … |
✓ |
… |
| | Core-Envelope Structure with <math>~(n_c,n_e) = (\tfrac{3}{2},3)</math> … |
✓ |
… |
|
|
|
Stability:
Two-Dimensional Configurations
Structure:
Stability:
Dynamics:
Three-Dimensional Configurations
Structure:
|
Solution Strategies
|
|
Example Solutions:
|
Stability:
Dynamics:
Related Projects Underway
Appendices
See Also
