Difference between revisions of "User:Tohline/SphericallySymmetricConfigurations/SolutionStrategies"
(→Solution Strategies: copy edits from older version of "parent" page) |
|||
(23 intermediate revisions by the same user not shown) | |||
Line 1: | Line 1: | ||
<!-- __FORCETOC__ will force the creation of a Table of Contents --> | <!-- __FORCETOC__ will force the creation of a Table of Contents --> | ||
<!-- __NOTOC__ will force TOC off --> | <!-- __NOTOC__ will force TOC off --> | ||
=Spherically Symmetric Configurations (Part II)= | |||
{{LSU_HBook_header}} | {{LSU_HBook_header}} | ||
<!-- [[Image:LSU_Structure_still.gif|74px|left]] --> | |||
Equilibrium, spherically symmetric '''structures''' are obtained by searching for time-independent solutions to the [[User:Tohline/SphericallySymmetricConfigurations/PGE#Spherically_Symmetric_Configurations_.28Part_I.29|identified set of simplified governing equations]]. The steady-state flow field that must be adopted to satisfy both a spherically symmetric geometry and the time-independent constraint is, | |||
<div align="center"> | |||
<math>~\vec{v} = \hat{e}_r v_r = 0 \, .</math> | |||
</div> | |||
After setting the radial velocity, <math>~v_r</math>, and all time-derivatives to zero, we see that the 1<sup>st</sup> (continuity) and 3<sup>rd</sup> (first law of thermodynamics) equations are trivially satisfied while the 2<sup>nd</sup> (Euler) and 4<sup>th</sup> give, respectively, | |||
<div align="center"> | <div align="center"> | ||
<span id=" | <span id="HydrostaticBalance"><font color="#770000">'''Hydrostatic Balance'''</font></span><br /> | ||
<math>\frac{1}{\rho}\frac{dP}{dr} =- \frac{d\Phi}{dr} </math> ,<br /> | <math>\frac{1}{\rho}\frac{dP}{dr} =- \frac{d\Phi}{dr} </math> ,<br /> | ||
</div> | </div> | ||
and, | |||
<div align="center"> | <div align="center"> | ||
<span id="Poisson"><font color="#770000">'''Poisson Equation'''</font></span><br /> | <span id="Poisson"><font color="#770000">'''Poisson Equation'''</font></span><br /> | ||
<math>\frac{1}{r^2} | <math>\frac{1}{r^2} \frac{d }{dr} \biggl( r^2 \frac{d \Phi}{dr} \biggr) = 4\pi G \rho </math> .<br /> | ||
</div> | </div> | ||
(We recognize the first of these expressions as being the statement of [[User:Tohline/PGE/ConservingMomentum#Time-independent_Behavior|hydrostatic balance]] appropriate for spherically symmetric configurations.) | |||
We need one supplemental relation to close this set of equations because there are two equations, but three unknown functions — {{User:Tohline/Math/VAR_Pressure01}}<math>(r)~</math>, {{User:Tohline/Math/VAR_Density01}}<math>(r)~</math>, and {{User:Tohline/Math/VAR_NewtonianPotential01}}<math>(r)~</math>. As has been outlined in our discussion of [[User:Tohline/SR#Time-Independent_Problems|supplemental relations for time-independent problems]] — and as is discussed further, below — in the context of this H_Book we will close this set of equations by specifying a structural, barotropic relationship between {{User:Tohline/Math/VAR_Pressure01}} and {{User:Tohline/Math/VAR_Density01}}. | |||
==Solution Strategies== | ==Solution Strategies== | ||
When attempting to solve the identified pair of simplified governing differential equations, it will be useful to note that, in a spherically symmetric configuration (where {{User:Tohline/Math/VAR_Density01}} is not a function of <math>\theta</math> or <math>\varphi</math>), the differential mass <math>dm_r</math> that is enclosed within a spherical shell of thickness <math>dr</math> is, | When attempting to solve the identified pair of simplified governing differential equations, it will be useful to note that, in a spherically symmetric configuration (where {{User:Tohline/Math/VAR_Density01}} is not a function of <math>~\theta</math> or <math>~\varphi</math>), the differential mass <math>~dm_r</math> that is enclosed within a spherical shell of thickness <math>~dr</math> is, | ||
<div align="center"> | <div align="center"> | ||
<math>dm_r = \rho dr \oint dS = r^2 \rho dr \int_0^\pi \sin\theta d\theta \int_0^{2\pi} d\varphi = 4\pi r^2 \rho dr</math> , | <math>~dm_r = \rho dr \oint dS = r^2 \rho dr \int_0^\pi \sin\theta d\theta \int_0^{2\pi} d\varphi = 4\pi r^2 \rho dr</math> , | ||
</div> | </div> | ||
where we have pulled from the [http://en.wikipedia.org/wiki/Spherical_coordinate_system#Integration_and_differentiation_in_spherical_coordinates Wikipedia discussion of integration and differentiation in spherical coordinates] to define the spherical surface element <math>dS</math>. Integrating from the center of the spherical configuration | |||
where we have pulled from the [http://en.wikipedia.org/wiki/Spherical_coordinate_system#Integration_and_differentiation_in_spherical_coordinates Wikipedia discussion of integration and differentiation in spherical coordinates] to define the spherical surface element, <math>~dS</math>. Integrating from the center of the spherical configuration <math>~(r=0)</math> out to some finite radius, <math>~r</math>, that is still inside the configuration gives the mass enclosed within that radius, <math>~M_r</math>; specifically, | |||
<div align="center"> | <div align="center"> | ||
<math>M_r \equiv \int_0^r dm_r = \int_0^r 4\pi r^2 \rho dr</math> . | <math>~M_r \equiv \int_0^r dm_r = \int_0^r 4\pi r^2 \rho dr</math> . | ||
</div> | </div> | ||
We can also state that, | We can also state that, | ||
<div align="center"> | <div align="center"> | ||
{{User:Tohline/Math/EQ_SSmassConservation01}} | |||
</div> | </div> | ||
This differential relation is often identified as a statement of mass conservation that replaces the equation of continuity for spherically symmetric, static equilibrium structures. | This differential relation is often identified as a statement of mass conservation that replaces the equation of continuity for spherically symmetric, static equilibrium structures. | ||
===Technique | ===Technique 1=== | ||
Integrating the Poisson equation once, from the center of the configuration | Integrating the Poisson equation once, from the center of the configuration <math>~(r=0)</math> out to some finite radius, <math>~r</math>, that is still inside the configuration, gives, | ||
<div align="center"> | <div align="center"> | ||
<math> | <math> | ||
\int_0^r d\biggl( r^2 \frac{d \Phi}{dr} \biggr) = \int_0^r 4\pi G r^2 \rho dr | ~\int_0^r d\biggl( r^2 \frac{d \Phi}{dr} \biggr) = \int_0^r 4\pi G r^2 \rho dr | ||
</math><br /> | </math><br /> | ||
<math> | <math> | ||
\Rightarrow ~~~~~ r^2 \frac{d \Phi}{dr} \biggr|_0^r = GM_r | \Rightarrow ~~~~~ r^2 \frac{d \Phi}{dr} \biggr|_0^r = GM_r \, . | ||
</math> | </math> | ||
</div> | </div> | ||
Now, as long as <math>d\Phi/dr</math> increases less steeply than <math>r^{-2}</math> as we move toward the center of the configuration — indeed, we will find that | Now, as long as <math>~d\Phi/dr</math> increases less steeply than <math>~r^{-2}</math> as we move toward the center of the configuration — indeed, we will find that <math>~d\Phi/dr</math> usually goes smoothly to zero at the center — the term on the left-hand-side of this last expression will go to zero at <math>~r=0</math>. Hence, this first integration of the Poisson equation gives, | ||
<div align="center"> | <div align="center"> | ||
<math> | <math> | ||
\frac{d \Phi}{dr} = \frac{G M_r}{r^2} | ~\frac{d \Phi}{dr} = \frac{G M_r}{r^2} \, . | ||
</math> | </math> | ||
</div> | </div> | ||
Substituting this expression into the hydrostatic | |||
Substituting this expression into the hydrostatic balance equation gives, | |||
<div align="center"> | <div align="center"> | ||
{{User:Tohline/Math/EQ_SShydrostaticBalance01}} | |||
</div> | </div> | ||
that is, a single governing integro-differential equation which depends only on the two unknown functions, {{User:Tohline/Math/VAR_Pressure01}} and {{User:Tohline/Math/VAR_Density01}} . | that is, a single governing integro-differential equation which depends only on the two unknown functions, {{User:Tohline/Math/VAR_Pressure01}} and {{User:Tohline/Math/VAR_Density01}} . | ||
===Technique 2=== | |||
As long as we are examining only barotropic structures, we can replace <math>~dP/\rho</math> by <math>~dH</math> in the hydrostatic balance relation to obtain, | |||
<div align="center"> | |||
<math>~\frac{dH}{dr} =- \frac{d\Phi}{dr} \, .</math> | |||
</div> | |||
If we multiply this expression through by <math>~r^2</math> then differentiate it with respect to <math>~r</math>, we obtain, | |||
<div align="center"> | |||
<math>\frac{d}{dr}\biggl( r^2 \frac{dH}{dr} \biggr) =- \frac{d}{dr} \biggl( r^2 \frac{d\Phi}{dr} \biggr) \, ,</math> | |||
</div> | |||
which can be used to replace the left-hand-side of the Poisson equation and give, | |||
<div align="center"> | |||
<math>~\frac{1}{r^2} \frac{d}{dr}\biggl( r^2 \frac{dH}{dr} \biggr) =- 4\pi G \rho \, ,</math> | |||
</div> | |||
that is, a single second-order governing differential equation which depends only on the two unknown functions, {{User:Tohline/Math/VAR_Enthalpy01}} and {{User:Tohline/Math/VAR_Density01}}. | |||
<b>Numerical integration examples:</b> | |||
* Isothermal sphere — | |||
** [[User:Tohline/SSC/Structure/IsothermalSphere#Emden.27s_Numerical_Solution|Emden's (1907) tabulation]]. | |||
** Tabulation by [http://adsabs.harvard.edu/abs/1949ApJ...109..551C Chandrasekhar & Wares (1949, ApJ, 109, 551)]. | |||
** Outline of [[User:Tohline/SSC/Structure/IsothermalSphere#Our_Numerical_Integration|our numerical integration scheme]]. | |||
* Spherical polytropes — | |||
** [[User:Tohline/SSC/Structure/Polytropes#Tabulated_Properties|Some tabulated global properties]]. | |||
** Outline of [[User:Tohline/SSC/Structure/Polytropes#Straight_Numerical_Integration|our numerical integration scheme]]. | |||
===Technique 3=== | |||
As in Technique #2, we replace <math>~dP/\rho</math> by <math>~dH</math> in the hydrostatic balance relation, but this time we realize that the resulting expression can be written in the form, | |||
<div align="center"> | |||
<math>~\frac{d}{dr}(H+\Phi) = 0 \, .</math> | |||
</div> | |||
This means that, throughout our configuration, the functions {{User:Tohline/Math/VAR_Enthalpy01}}<math>(\rho)~</math> and {{User:Tohline/Math/VAR_NewtonianPotential01}}<math>(\rho)~</math> must sum to a constant value, call it <math>~C_\mathrm{B}</math>. That is to say, the statement of hydrostatic balance reduces to the ''algebraic'' expression, | |||
<div align="center"> | |||
<math>H + \Phi = C_\mathrm{B} \, .</math> | |||
</div> | |||
This relation must be solved in conjunction with the Poisson equation, | |||
<div align="center"> | |||
<math>~\frac{1}{r^2} \frac{d }{dr} \biggl( r^2 \frac{d \Phi}{dr} \biggr) = 4\pi G \rho \, ,</math> | |||
</div> | |||
giving us two equations (one algebraic and the other a <math>2^\mathrm{nd}</math>-order ODE) that relate the three unknown functions, {{User:Tohline/Math/VAR_Enthalpy01}}, {{User:Tohline/Math/VAR_Density01}}, and {{User:Tohline/Math/VAR_NewtonianPotential01}} to one another. | |||
<b>Self-Consistent Field (SCF) Technique:</b> | |||
* Spherical polytropes — | |||
** [[User:Tohline/SSC/Structure/Polytropes#Tabulated_Properties|Some tabulated global properties]]. | |||
** Outline of [[User:Tohline/SSC/Structure/Polytropes#HSCF_Technique|our implementation of the HSCF scheme]]. | |||
=See Also= | |||
* Part I of ''Spherically Symmetric Configurations'': [[User:Tohline/SphericallySymmetricConfigurations/PGE#Spherically_Symmetric_Configurations_.28Part_I.29|Simplified Governing Equations]] | |||
{{LSU_HBook_footer}} | {{LSU_HBook_footer}} |
Latest revision as of 21:58, 10 April 2017
Spherically Symmetric Configurations (Part II)
| Tiled Menu | Tables of Content | Banner Video | Tohline Home Page | |
Equilibrium, spherically symmetric structures are obtained by searching for time-independent solutions to the identified set of simplified governing equations. The steady-state flow field that must be adopted to satisfy both a spherically symmetric geometry and the time-independent constraint is,
<math>~\vec{v} = \hat{e}_r v_r = 0 \, .</math>
After setting the radial velocity, <math>~v_r</math>, and all time-derivatives to zero, we see that the 1st (continuity) and 3rd (first law of thermodynamics) equations are trivially satisfied while the 2nd (Euler) and 4th give, respectively,
Hydrostatic Balance
<math>\frac{1}{\rho}\frac{dP}{dr} =- \frac{d\Phi}{dr} </math> ,
and,
Poisson Equation
<math>\frac{1}{r^2} \frac{d }{dr} \biggl( r^2 \frac{d \Phi}{dr} \biggr) = 4\pi G \rho </math> .
(We recognize the first of these expressions as being the statement of hydrostatic balance appropriate for spherically symmetric configurations.)
We need one supplemental relation to close this set of equations because there are two equations, but three unknown functions — <math>~P</math><math>(r)~</math>, <math>~\rho</math><math>(r)~</math>, and <math>~\Phi</math><math>(r)~</math>. As has been outlined in our discussion of supplemental relations for time-independent problems — and as is discussed further, below — in the context of this H_Book we will close this set of equations by specifying a structural, barotropic relationship between <math>~P</math> and <math>~\rho</math>.
Solution Strategies
When attempting to solve the identified pair of simplified governing differential equations, it will be useful to note that, in a spherically symmetric configuration (where <math>~\rho</math> is not a function of <math>~\theta</math> or <math>~\varphi</math>), the differential mass <math>~dm_r</math> that is enclosed within a spherical shell of thickness <math>~dr</math> is,
<math>~dm_r = \rho dr \oint dS = r^2 \rho dr \int_0^\pi \sin\theta d\theta \int_0^{2\pi} d\varphi = 4\pi r^2 \rho dr</math> ,
where we have pulled from the Wikipedia discussion of integration and differentiation in spherical coordinates to define the spherical surface element, <math>~dS</math>. Integrating from the center of the spherical configuration <math>~(r=0)</math> out to some finite radius, <math>~r</math>, that is still inside the configuration gives the mass enclosed within that radius, <math>~M_r</math>; specifically,
<math>~M_r \equiv \int_0^r dm_r = \int_0^r 4\pi r^2 \rho dr</math> .
We can also state that,
This differential relation is often identified as a statement of mass conservation that replaces the equation of continuity for spherically symmetric, static equilibrium structures.
Technique 1
Integrating the Poisson equation once, from the center of the configuration <math>~(r=0)</math> out to some finite radius, <math>~r</math>, that is still inside the configuration, gives,
<math>
~\int_0^r d\biggl( r^2 \frac{d \Phi}{dr} \biggr) = \int_0^r 4\pi G r^2 \rho dr
</math>
<math> \Rightarrow ~~~~~ r^2 \frac{d \Phi}{dr} \biggr|_0^r = GM_r \, . </math>
Now, as long as <math>~d\Phi/dr</math> increases less steeply than <math>~r^{-2}</math> as we move toward the center of the configuration — indeed, we will find that <math>~d\Phi/dr</math> usually goes smoothly to zero at the center — the term on the left-hand-side of this last expression will go to zero at <math>~r=0</math>. Hence, this first integration of the Poisson equation gives,
<math> ~\frac{d \Phi}{dr} = \frac{G M_r}{r^2} \, . </math>
Substituting this expression into the hydrostatic balance equation gives,
that is, a single governing integro-differential equation which depends only on the two unknown functions, <math>~P</math> and <math>~\rho</math> .
Technique 2
As long as we are examining only barotropic structures, we can replace <math>~dP/\rho</math> by <math>~dH</math> in the hydrostatic balance relation to obtain,
<math>~\frac{dH}{dr} =- \frac{d\Phi}{dr} \, .</math>
If we multiply this expression through by <math>~r^2</math> then differentiate it with respect to <math>~r</math>, we obtain,
<math>\frac{d}{dr}\biggl( r^2 \frac{dH}{dr} \biggr) =- \frac{d}{dr} \biggl( r^2 \frac{d\Phi}{dr} \biggr) \, ,</math>
which can be used to replace the left-hand-side of the Poisson equation and give,
<math>~\frac{1}{r^2} \frac{d}{dr}\biggl( r^2 \frac{dH}{dr} \biggr) =- 4\pi G \rho \, ,</math>
that is, a single second-order governing differential equation which depends only on the two unknown functions, <math>~H</math> and <math>~\rho</math>.
Numerical integration examples:
- Isothermal sphere —
- Emden's (1907) tabulation.
- Tabulation by Chandrasekhar & Wares (1949, ApJ, 109, 551).
- Outline of our numerical integration scheme.
- Spherical polytropes —
Technique 3
As in Technique #2, we replace <math>~dP/\rho</math> by <math>~dH</math> in the hydrostatic balance relation, but this time we realize that the resulting expression can be written in the form,
<math>~\frac{d}{dr}(H+\Phi) = 0 \, .</math>
This means that, throughout our configuration, the functions <math>~H</math><math>(\rho)~</math> and <math>~\Phi</math><math>(\rho)~</math> must sum to a constant value, call it <math>~C_\mathrm{B}</math>. That is to say, the statement of hydrostatic balance reduces to the algebraic expression,
<math>H + \Phi = C_\mathrm{B} \, .</math>
This relation must be solved 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>
giving us two equations (one algebraic and the other a <math>2^\mathrm{nd}</math>-order ODE) that relate the three unknown functions, <math>~H</math>, <math>~\rho</math>, and <math>~\Phi</math> to one another.
Self-Consistent Field (SCF) Technique:
- Spherical polytropes —
See Also
- Part I of Spherically Symmetric Configurations: Simplified Governing Equations
© 2014 - 2021 by Joel E. Tohline |