AxisymmetricConfigurations/PGE - Revision history

Jet53man: /* Axisymmetric Configurations (Governing Equations) */

2021-11-01T21:37:55Z

Axisymmetric Configurations (Governing Equations)

← Older revision		Revision as of 17:37, 1 November 2021
Line 4:		Line 4:
	{\| class="AxisymmetricConfigurationsPGE" style="float:left; margin-right: 20px; border-style: solid; border-width: 3px border-color: black"		{\| class="AxisymmetricConfigurationsPGE" style="float:left; margin-right: 20px; border-style: solid; border-width: 3px border-color: black"
	\|-		\|-
	! style="height: ~~150px~~; width: ~~150px~~; background-color:lightgreen;" \|[[H_BookTiledMenu#Two-Dimensional_Configurations_.28Axisymmetric.29\|<b>PGEs<br />for<br />Axisymmetric<br />Systems</b>]]		! style="height: 125px; width: 125px; background-color:lightgreen;" \|[[H_BookTiledMenu#Two-Dimensional_Configurations_.28Axisymmetric.29\|<b>PGEs<br />for<br />Axisymmetric<br />Systems</b>]]
	\|}		\|}
	If the self-gravitating configuration that we wish to construct is axisymmetric, then the coupled set of multidimensional, partial differential equations that serve as our [[PGE\|principal governing equations]] can be simplified, resulting in a coupled set of two-dimensional PDEs.		If the self-gravitating configuration that we wish to construct is axisymmetric, then the coupled set of multidimensional, partial differential equations that serve as our [[PGE\|principal governing equations]] can be simplified, resulting in a coupled set of two-dimensional PDEs.

Jet53man at 21:07, 9 September 2021

2021-09-09T21:07:16Z

← Older revision		Revision as of 17:07, 9 September 2021
Line 6:		Line 6:
	! style="height: 150px; width: 150px; background-color:lightgreen;" \|[[H_BookTiledMenu#Two-Dimensional_Configurations_.28Axisymmetric.29\|<b>PGEs<br />for<br />Axisymmetric<br />Systems</b>]]		! style="height: 150px; width: 150px; background-color:lightgreen;" \|[[H_BookTiledMenu#Two-Dimensional_Configurations_.28Axisymmetric.29\|<b>PGEs<br />for<br />Axisymmetric<br />Systems</b>]]
	\|}		\|}
	If the self-gravitating configuration that we wish to construct is axisymmetric, then the coupled set of multidimensional, partial differential equations that serve as our [[~~User:Tohline/~~PGE\|principal governing equations]] can be simplified, resulting in a coupled set of two-dimensional PDEs.		If the self-gravitating configuration that we wish to construct is axisymmetric, then the coupled set of multidimensional, partial differential equations that serve as our [[PGE\|principal governing equations]] can be simplified, resulting in a coupled set of two-dimensional PDEs.
	<br />		<br />
	<br />		<br />
Line 44:		Line 44:
	<tr>		<tr>
	<td align="center" colspan="3">		<td align="center" colspan="3">
	[<b>[[~~User:Tohline/~~Appendix/References#BT87\|<font color="red">BT87</font>]]</b>], p. 649, Eq. (1B-37)		[<b>[[Appendix/References#BT87\|<font color="red">BT87</font>]]</b>], p. 649, Eq. (1B-37)
	</td>		</td>
	</tr>		</tr>
Line 65:		Line 65:
	<tr>		<tr>
	<td align="center" colspan="3">		<td align="center" colspan="3">
	[<b>[[~~User:Tohline/~~Appendix/References#BT87\|<font color="red">BT87</font>]]</b>], p. 650, Eq. (1B-50)		[<b>[[Appendix/References#BT87\|<font color="red">BT87</font>]]</b>], p. 650, Eq. (1B-50)
	</td>		</td>
	</tr>		</tr>
Line 104:		Line 104:
	<tr>		<tr>
	<td align="center" colspan="3">		<td align="center" colspan="3">
	[<b>[[~~User:Tohline/~~Appendix/References#BT87\|<font color="red">BT87</font>]]</b>], p. 650, Eq. (1B-45)		[<b>[[Appendix/References#BT87\|<font color="red">BT87</font>]]</b>], p. 650, Eq. (1B-45)
	</td>		</td>
	</tr>		</tr>
Line 130:		Line 130:
	<tr>		<tr>
	<td align="center" colspan="3">		<td align="center" colspan="3">
	[<b>[[~~User:Tohline/~~Appendix/References#BT87\|<font color="red">BT87</font>]]</b>], p. 651, Eq. (1B-54)		[<b>[[Appendix/References#BT87\|<font color="red">BT87</font>]]</b>], p. 651, Eq. (1B-54)
	</td>		</td>
	</tr>		</tr>
Line 216:		Line 216:
	<tr>		<tr>
	<td align="center" colspan="3">		<td align="center" colspan="3">
	[<b>[[~~User:Tohline/~~Appendix/References#BT87\|<font color="red">BT87</font>]]</b>], p. 647, Eq. (1B-23)		[<b>[[Appendix/References#BT87\|<font color="red">BT87</font>]]</b>], p. 647, Eq. (1B-23)
	</td>		</td>
	</tr>		</tr>
Line 226:		Line 226:
	===Governing Equations (CYL.)===		===Governing Equations (CYL.)===

	Introducing the above expressions into the [[~~User:Tohline/~~PGE\|principal governing equations]] gives,		Introducing the above expressions into the [[PGE\|principal governing equations]] gives,

	<div align="center">		<div align="center">
Line 249:		Line 249:
	<font color="#770000">'''First Law of Thermodynamics'''</font></span><br />		<font color="#770000">'''First Law of Thermodynamics'''</font></span><br />

	{{~~User:Tohline/~~Math/EQ_FirstLaw02}}		{{Math/EQ_FirstLaw02}}


Line 428:		Line 428:
	<tr>		<tr>
	<td align="center" colspan="3">		<td align="center" colspan="3">
	[<b>[[~~User:Tohline/~~Appendix/References#BT87\|<font color="red">BT87</font>]]</b>], p. 649, Eq. (1B-38)		[<b>[[Appendix/References#BT87\|<font color="red">BT87</font>]]</b>], p. 649, Eq. (1B-38)
	</td>		</td>
	</tr>		</tr>
Line 451:		Line 451:
	<tr>		<tr>
	<td align="center" colspan="3">		<td align="center" colspan="3">
	[<b>[[~~User:Tohline/~~Appendix/References#BT87\|<font color="red">BT87</font>]]</b>], p. 650, Eq. (1B-51)		[<b>[[Appendix/References#BT87\|<font color="red">BT87</font>]]</b>], p. 650, Eq. (1B-51)
	</td>		</td>
	</tr>		</tr>
Line 492:		Line 492:
	<tr>		<tr>
	<td align="center" colspan="3">		<td align="center" colspan="3">
	[<b>[[~~User:Tohline/~~Appendix/References#BT87\|<font color="red">BT87</font>]]</b>], p. 650, Eq. (1B-46)		[<b>[[Appendix/References#BT87\|<font color="red">BT87</font>]]</b>], p. 650, Eq. (1B-46)
	</td>		</td>
	</tr>		</tr>
Line 546:		Line 546:
	<tr>		<tr>
	<td align="center" colspan="3">		<td align="center" colspan="3">
	[<b>[[~~User:Tohline/~~Appendix/References#BT87\|<font color="red">BT87</font>]]</b>], p. 651, Eq. (1B-55)		[<b>[[Appendix/References#BT87\|<font color="red">BT87</font>]]</b>], p. 651, Eq. (1B-55)
	</td>		</td>
	</tr>		</tr>
Line 652:		Line 652:
	<tr>		<tr>
	<td align="center" colspan="3">		<td align="center" colspan="3">
	[<b>[[~~User:Tohline/~~Appendix/References#BT87\|<font color="red">BT87</font>]]</b>], p. 648, Eq. (1B-30)		[<b>[[Appendix/References#BT87\|<font color="red">BT87</font>]]</b>], p. 648, Eq. (1B-30)
	</td>		</td>
	</tr>		</tr>
Line 660:		Line 660:
	===Governing Equations (SPH.)===		===Governing Equations (SPH.)===

	Introducing the above expressions into the [[~~User:Tohline/~~PGE\|principal governing equations]] gives,		Introducing the above expressions into the [[PGE\|principal governing equations]] gives,

	<div align="center">		<div align="center">
Line 743:		Line 743:
	<font color="#770000">'''First Law of Thermodynamics'''</font></span><br />		<font color="#770000">'''First Law of Thermodynamics'''</font></span><br />

	{{~~User:Tohline/~~Math/EQ_FirstLaw02}}		{{Math/EQ_FirstLaw02}}

Jet53man: Created page with "FORCETOC =Axisymmetric Configurations (Governing Equations)= {| class="AxisymmetricConfigurationsPGE" style="float:left; margin-righ..."

2021-09-09T21:03:41Z

Created page with "__FORCETOC__  =Axisymmetric Configurations (Governing Equations)= {| class="AxisymmetricConfigurationsPGE" style="float:left; margin-righ..."

New page

__FORCETOC__

=Axisymmetric Configurations (Governing Equations)=
{| class="AxisymmetricConfigurationsPGE" style="float:left; margin-right: 20px; border-style: solid; border-width: 3px border-color: black"
|-
! style="height: 150px; width: 150px; background-color:lightgreen;" |[[H_BookTiledMenu#Two-Dimensional_Configurations_.28Axisymmetric.29|<b>PGEs<br />for<br />Axisymmetric<br />Systems</b>]]
|}
If the self-gravitating configuration that we wish to construct is axisymmetric, then the coupled set of multidimensional, partial differential equations that serve as our [[User:Tohline/PGE|principal governing equations]] can be simplified, resulting in a coupled set of two-dimensional PDEs.
  <br />
  <br />
  <br />
  <br />
  <br />
  <br />

==Cylindrical Coordinate Base==
Here we choose to …

<ol>
<li>Express each of the multidimensional spatial operators in cylindrical coordinates (<math>\varpi, \varphi, z</math>) (see, for example, the [http://en.wikipedia.org/wiki/Del_in_cylindrical_and_spherical_coordinates Wikipedia discussion of vector calculus formulae in cylindrical coordinates]) and set to zero all spatial derivatives that are taken with respect to the angular coordinate <math>\varphi</math>:

<table align="center" border="0" cellpadding="5">
<tr>
<td colspan="3" align="center">
<font color="#770000"><b>Spatial Operators in Cylindrical Coordinates</b></font>
</td>
</tr>

<tr>
<td align="right">
<math>
\nabla f
</math>
</td>
<td align="center">
=
</td>
<td align="left">
<math>
{\hat{e}}_\varpi \biggl[ \frac{\partial f}{\partial\varpi} \biggr] + {\hat{e}}_\varphi \cancel{\biggl[ \frac{1}{\varpi} \frac{\partial f}{\partial\varphi} \biggr]} + {\hat{e}}_z \biggl[ \frac{\partial f}{\partial z} \biggr] ;
</math>
</td>
</tr>
<tr>
<td align="center" colspan="3">
[<b>[[User:Tohline/Appendix/References#BT87|<font color="red">BT87</font>]]</b>], p. 649, Eq. (1B-37)
</td>
</tr>

<tr>
<td align="right">
<math>
\nabla^2 f
</math>
</td>
<td align="center">
=
</td>
<td align="left">
<math>
\frac{1}{\varpi} \frac{\partial }{\partial\varpi} \biggl[ \varpi \frac{\partial f}{\partial\varpi} \biggr] + \cancel{\frac{1}{\varpi^2} \frac{\partial^2 f}{\partial\varphi^2}} + \frac{\partial^2 f}{\partial z^2} ;
</math>
</td>
</tr>
<tr>
<td align="center" colspan="3">
[<b>[[User:Tohline/Appendix/References#BT87|<font color="red">BT87</font>]]</b>], p. 650, Eq. (1B-50)
</td>
</tr>

<tr>
<td align="right">
<math>
(\vec{v}\cdot\nabla)f
</math>
</td>
<td align="center">
=
</td>
<td align="left">
<math>
\biggl[ v_\varpi \frac{\partial f}{\partial\varpi} \biggr] +
\cancel{\biggl[ \frac{v_\varphi}{\varpi} \frac{\partial f}{\partial\varphi} \biggr]} +
\biggl[ v_z \frac{\partial f}{\partial z} \biggr] ;
</math>
</td>
</tr>

<tr>
<td align="right">
<math>
\nabla \cdot \vec{F}
</math>
</td>
<td align="center">
=
</td>
<td align="left">
<math>
\frac{1}{\varpi} \frac{\partial (\varpi F_\varpi)}{\partial\varpi} + \cancel{\frac{1}{\varpi} \frac{\partial F_\varphi}{\partial\varphi}} + \frac{\partial F_z}{\partial z} ;
</math>
</td>
</tr>
<tr>
<td align="center" colspan="3">
[<b>[[User:Tohline/Appendix/References#BT87|<font color="red">BT87</font>]]</b>], p. 650, Eq. (1B-45)
</td>
</tr>
</table>

<table border="0" cellpadding="5" align="center">

<tr>
<td align="right">
<math>
(\vec{F} \cdot \nabla )\vec{B}
</math>
</td>
<td align="center">
=
</td>
<td align="left">
<math>
\hat{e}_\varpi \biggl[ F_\varpi \frac{\partial B_\varpi}{\partial\varpi} + \cancel{\frac{F_\varphi}{\varpi} \frac{\partial B_\varpi}{\partial\varphi}} + F_z \frac{\partial B_\varpi}{\partial z} - \frac{F_\varphi B_\varphi}{\varpi} \biggr]
+ \hat{e}_\varphi \biggl[ F_\varpi \frac{\partial B_\varphi}{\partial \varpi} + \cancel{ \frac{F_\varphi}{\varpi} \frac{\partial B_\varphi}{\partial\varphi}} + F_z \frac{\partial B_\varphi}{\partial z} + \frac{F_\varphi B_\varpi}{\varpi} \biggr]
+ \hat{e}_z \biggl[ F_\varpi \frac{\partial B_z}{\partial\varpi} +\cancel{ \frac{F_\varphi}{\varpi} \frac{\partial B_z}{\partial \varphi}} + F_z \frac{\partial B_z}{\partial z} \biggr] \, .
</math>
</td>
</tr>
<tr>
<td align="center" colspan="3">
[<b>[[User:Tohline/Appendix/References#BT87|<font color="red">BT87</font>]]</b>], p. 651, Eq. (1B-54)
</td>
</tr>
</table>
<span id="CYLconvectiveOperator">From this last expression — the so-called ''convective operator'' — we conclude as well that, for axisymmetric systems,</span>

<table border="0" cellpadding="5" align="center">

<tr>
<td align="right">
<math>
(\vec{v} \cdot \nabla )\vec{v}
</math>
</td>
<td align="center">
=
</td>
<td align="left">
<math>
\hat{e}_\varpi \biggl[ v_\varpi \frac{\partial v_\varpi}{\partial\varpi} + v_z \frac{\partial v_\varpi}{\partial z} - \frac{v_\varphi v_\varphi}{\varpi} \biggr]
+ \hat{e}_\varphi \biggl[ v_\varpi \frac{\partial v_\varphi}{\partial \varpi} + v_z \frac{\partial v_\varphi}{\partial z} + \frac{v_\varphi v_\varpi}{\varpi} \biggr]
+ \hat{e}_z \biggl[ v_\varpi \frac{\partial v_z}{\partial\varpi} + v_z \frac{\partial v_z}{\partial z} \biggr] \, .
</math>
</td>
</tr>
</table>

<li>Express all vector time-derivatives in cylindrical coordinates:

<table align="center" border="0" cellpadding="5">
<tr>
<td colspan="3" align="center">
<font color="#770000"><b>Vector Time-Derivatives in Cylindrical Coordinates</b></font>
</td>
</tr>

<tr>
<td align="right">
<math>
\frac{d}{dt}\vec{F}
</math>
</td>
<td align="center">
=
</td>
<td align="left">
<math>
{\hat{e}}_\varpi \frac{dF_\varpi}{dt} + F_\varpi \frac{d{\hat{e}}_\varpi}{dt} + {\hat{e}}_\varphi \frac{dF_\varphi}{dt} + F_\varphi \frac{d{\hat{e}}_\varphi}{dt} + {\hat{e}}_z \frac{dF_z}{dt} + F_z \frac{d{\hat{e}}_z}{dt}
</math>
</td>
</tr>

<tr>
<td align="right">
 
</td>
<td align="center">
=
</td>
<td align="left">
<math>
{\hat{e}}_\varpi \biggl[ \frac{dF_\varpi}{dt} - F_\varphi \dot\varphi \biggr] + {\hat{e}}_\varphi \biggl[ \frac{dF_\varphi}{dt} + F_\varpi \dot\varphi \biggr] + {\hat{e}}_z \frac{dF_z}{dt} ;
</math>
</td>
</tr>

<tr>
<td align="right">
<math>
\vec{v} = \frac{d\vec{x}}{dt} = \frac{d}{dt}\biggl[ \hat{e}_\varpi \varpi + \hat{e}_z z \biggr]
</math>
</td>
<td align="center">
=
</td>
<td align="left">
<math>
{\hat{e}}_\varpi \biggl[ \dot\varpi \biggr] +
{\hat{e}}_\varphi \biggl[ \varpi \dot\varphi \biggr] +
{\hat{e}}_z \biggl[ \dot{z} \biggr] .
</math>
</td>
</tr>
<tr>
<td align="center" colspan="3">
[<b>[[User:Tohline/Appendix/References#BT87|<font color="red">BT87</font>]]</b>], p. 647, Eq. (1B-23)
</td>
</tr>
</table>

</ol>

===Governing Equations (CYL.)===

Introducing the above expressions into the [[User:Tohline/PGE|principal governing equations]] gives,

<div align="center">
<span id="Continuity"><font color="#770000">'''Equation of Continuity'''</font></span><br />

<math>\frac{d\rho}{dt} + \frac{\rho}{\varpi} \frac{\partial}{\partial\varpi} \biggl[ \varpi \dot\varpi \biggr]
+ \rho \frac{\partial}{\partial z} \biggl[ \rho \dot{z} \biggr] = 0 </math><br />

<span id="PGE:Euler">
<font color="#770000">'''Euler Equation'''</font>
</span><br />

<math>
{\hat{e}}_\varpi \biggl[ \frac{d \dot\varpi}{dt} - \varpi {\dot\varphi}^2 \biggr] + {\hat{e}}_\varphi \biggl[ \frac{d(\varpi\dot\varphi)}{dt} + \dot\varpi \dot\varphi \biggr] + {\hat{e}}_z \biggl[ \frac{d \dot{z}}{dt} \biggr] = -
{\hat{e}}_\varpi \biggl[ \frac{1}{\rho}\frac{\partial P}{\partial\varpi} + \frac{\partial \Phi}{\partial\varpi}\biggr] - {\hat{e}}_z \biggl[ \frac{1}{\rho}\frac{\partial P}{\partial z} + \frac{\partial \Phi}{\partial z} \biggr]
</math><br />

<span id="PGE:AdiabaticFirstLaw">Adiabatic Form of the<br />
<font color="#770000">'''First Law of Thermodynamics'''</font></span><br />

{{User:Tohline/Math/EQ_FirstLaw02}}

<span id="PGE:Poisson"><font color="#770000">'''Poisson Equation'''</font></span><br />

<math>
\frac{1}{\varpi} \frac{\partial }{\partial\varpi} \biggl[ \varpi \frac{\partial \Phi}{\partial\varpi} \biggr] + \frac{\partial^2 \Phi}{\partial z^2} = 4\pi G \rho .
</math><br />
</div>

===Conservation of Specific Angular Momentum (CYL.)===

The <math>\hat{e}_\varphi</math> component of the Euler equation leads to a statement of conservation of specific angular momentum, <math>j</math>, as follows.

<div align="center">
<math>
\frac{d(\varpi\dot\varphi)}{dt} + \dot\varpi \dot\varphi = \frac{1}{\varpi}\biggl[ \varpi \frac{d(\varpi\dot\varphi)}{dt} + \varpi \dot\varpi \dot\varphi \biggr] =0
</math><br />

<math>
\Rightarrow ~~~~~ \frac{d(\varpi^2 \dot\varphi)}{dt} = 0
</math><br />

<math>
\Rightarrow ~~~~~ j(\varpi,z) \equiv \varpi^2 \dot\varphi = \mathrm{constant} ~(\mathrm{i.e.,}~\mathrm{independent~of~time})
</math><br />
</div>

<span id="RelevantCylindricalComponents">So, for axisymmetric configurations, the <math>\hat{e}_\varpi</math> and <math>\hat{e}_z</math> components of the Euler equation become, respectively,</span>
<table border="1" align="center" cellpadding="10"><tr><td align="center">
<table border="0" cellpadding="5" align="center">
<tr>
<td align="right"><math>~{\hat{e}}_\varpi</math>:    </td>
<td align="right">
<math>
\frac{d \dot\varpi}{dt} - \frac{j^2}{\varpi^3}
</math>
</td>
<td align="center">
=
</td>
<td align="left">
<math>
- \biggl[ \frac{1}{\rho}\frac{\partial P}{\partial\varpi} + \frac{\partial \Phi}{\partial\varpi}\biggr]
</math>
</td>
</tr>
<tr>
<td align="right"><math>~{\hat{e}}_z</math>:    </td>
<td align="right">
<math>
\frac{d \dot{z}}{dt}
</math>
</td>
<td align="center">
=
</td>
<td align="left">
<math>
- \biggl[ \frac{1}{\rho}\frac{\partial P}{\partial z} + \frac{\partial \Phi}{\partial z} \biggr]
</math>
</td>
</tr>
</table>
</td></tr></table>

===Eulerian Formulation (CYL.)===

Each of the above simplified governing equations has been written in terms of Lagrangian time derivatives. An Eulerian formulation of each equation can be obtained by replacing each Lagrangian time derivative by its Eulerian counterpart. Specifically, for any scalar function, <math>f</math>,

<div align="center">
<math>
\frac{df}{dt} \rightarrow \frac{\partial f}{\partial t} + (\vec{v}\cdot \nabla)f =
\frac{\partial f}{\partial t} + \biggl[ \dot\varpi \frac{\partial f}{\partial\varpi} \biggr] +
\biggl[ \dot{z} \frac{\partial f}{\partial z} \biggr] .
</math>
</div>

Making this substitution throughout the set of governing relations gives:

<div align="center">
<span id="Continuity"><font color="#770000">'''Equation of Continuity'''</font></span>

<math>\frac{\partial\rho}{\partial t} + \frac{1}{\varpi} \frac{\partial}{\partial\varpi} \biggl[ \rho \varpi \dot\varpi \biggr]
+ \frac{\partial}{\partial z} \biggl[ \rho \dot{z} \biggr] = 0 </math><br />

<span id="PGE:Euler">The Two Relevant Components of the<br />
<font color="#770000">'''Euler Equation'''</font>
</span>
<table border="0" cellpadding="5">
<tr>
<td align="right"><math>~{\hat{e}}_\varpi</math>:    </td>
<td align="right">
<math>~
\frac{\partial \dot\varpi}{\partial t} + \biggl[ \dot\varpi \frac{\partial \dot\varpi}{\partial\varpi} \biggr] +
\biggl[ \dot{z} \frac{\partial \dot\varpi}{\partial z} \biggr]
</math>
</td>
<td align="center">
<math>~=</math>
</td>
<td align="left">
<math>~
- \biggl[ \frac{1}{\rho}\frac{\partial P}{\partial\varpi} + \frac{\partial \Phi}{\partial\varpi}\biggr] + \frac{j^2}{\varpi^3}
</math>
</td>
</tr>
<tr>
<td align="right"><math>~{\hat{e}}_z</math>:    </td>
<td align="right">
<math>~
\frac{\partial \dot{z}}{\partial t} + \biggl[ \dot\varpi \frac{\partial \dot{z}}{\partial\varpi} \biggr] +
\biggl[ \dot{z} \frac{\partial \dot{z}}{\partial z} \biggr]
</math>
</td>
<td align="center">
<math>~=</math>
</td>
<td align="left">
<math>~
- \biggl[ \frac{1}{\rho}\frac{\partial P}{\partial z} + \frac{\partial \Phi}{\partial z} \biggr]
</math>
</td>
</tr>
</table>

<span id="PGE:AdiabaticFirstLaw">Adiabatic Form of the<br />
<font color="#770000">'''First Law of Thermodynamics'''</font></span><br />

<math>~
\biggl\{\frac{\partial \epsilon}{\partial t} + \biggl[ \dot\varpi \frac{\partial \epsilon}{\partial\varpi} \biggr] + \biggl[ \dot{z} \frac{\partial \epsilon}{\partial z} \biggr]\biggr\} +
P \biggl\{\frac{\partial }{\partial t}\biggl(\frac{1}{\rho}\biggr) +
\biggl[ \dot\varpi \frac{\partial }{\partial\varpi}\biggl(\frac{1}{\rho}\biggr) \biggr] +
\biggl[ \dot{z} \frac{\partial }{\partial z}\biggl(\frac{1}{\rho}\biggr) \biggr] \biggr\} = 0
</math>

<span id="PGE:Poisson"><font color="#770000">'''Poisson Equation'''</font></span><br />

<math>
\frac{1}{\varpi} \frac{\partial }{\partial\varpi} \biggl[ \varpi \frac{\partial \Phi}{\partial\varpi} \biggr] + \frac{\partial^2 \Phi}{\partial z^2} = 4\pi G \rho .
</math><br />
</div>

==Spherical Coordinate Base==
Here we choose to …

<ol>
<li>Express each of the multidimensional spatial operators in spherical coordinates (<math>r, \theta, \varphi</math>) (see, for example, the [http://en.wikipedia.org/wiki/Del_in_cylindrical_and_spherical_coordinates Wikipedia discussion of vector calculus formulae in spherical coordinates]) and set to zero all spatial derivatives that are taken with respect to the angular coordinate <math>\varphi</math>:

<table align="center" border="0" cellpadding="5">
<tr>
<td colspan="3" align="center">
<font color="#770000"><b>Spatial Operators in Spherical Coordinates</b></font>
</td>
</tr>

<tr>
<td align="right">
<math>
\nabla f
</math>
</td>
<td align="center">
=
</td>
<td align="left">
<math>
{\hat{e}}_r \biggl[ \frac{\partial f}{\partial r} \biggr]
+ {\hat{e}}_\theta \biggl[ \frac{1}{r} \frac{\partial f}{\partial\theta} \biggr]
+ {\hat{e}}_\varphi \cancel{\biggl[\frac{1}{r\sin\theta}~ \frac{\partial f}{\partial \varphi} \biggr]} ;
</math>
</td>
</tr>
<tr>
<td align="center" colspan="3">
[<b>[[User:Tohline/Appendix/References#BT87|<font color="red">BT87</font>]]</b>], p. 649, Eq. (1B-38)
</td>
</tr>

<tr>
<td align="right">
<math>
\nabla^2 f
</math>
</td>
<td align="center">
=
</td>
<td align="left">
<math>
\frac{1}{r^2} \frac{\partial }{\partial r} \biggl[ r^2 \frac{\partial f}{\partial r} \biggr]
+ \frac{1}{r^2 \sin\theta} \frac{\partial }{\partial \theta}\biggl(\sin\theta \frac{\partial f}{\partial\theta}\biggr)
+ \cancel{ \biggl[\frac{1}{r^2 \sin^2\theta} \frac{\partial^2 f}{\partial \varphi^2} \biggr]} ;
</math>
</td>
</tr>
<tr>
<td align="center" colspan="3">
[<b>[[User:Tohline/Appendix/References#BT87|<font color="red">BT87</font>]]</b>], p. 650, Eq. (1B-51)
</td>
</tr>

<tr>
<td align="right">
<math>
(\vec{v}\cdot\nabla)f
</math>
</td>
<td align="center">
=
</td>
<td align="left">
<math>
\biggl[ v_r \frac{\partial f}{\partial r} \biggr]
+ \biggl[ \frac{v_\theta}{r} \frac{\partial f}{\partial\theta} \biggr]
+ \cancel{\biggl[\frac{v_\varphi}{r\sin\theta}~ \frac{\partial f}{\partial \varphi} \biggr]} ;
</math>
</td>
</tr>

<tr>
<td align="right">
<math>
\nabla \cdot \vec{F}
</math>
</td>
<td align="center">
=
</td>
<td align="left">
<math>
\frac{1}{r^2} \frac{\partial (r^2 F_r)}{\partial r}
+ \frac{1}{r\sin\theta} \frac{\partial }{\partial\theta} \biggl( F_\theta \sin\theta \biggr)
+ \cancel{ \biggl[ \frac{1}{r\sin\theta}~\frac{\partial F_\varphi}{\partial \varphi} \biggr]} ;
</math>
</td>
</tr>
<tr>
<td align="center" colspan="3">
[<b>[[User:Tohline/Appendix/References#BT87|<font color="red">BT87</font>]]</b>], p. 650, Eq. (1B-46)
</td>
</tr>
</table>

<table border="0" cellpadding="5" align="center">

<tr>
<td align="right">
<math>
(\vec{F} \cdot \nabla )\vec{B}
</math>
</td>
<td align="center">
=
</td>
<td align="left">
<math>
\hat{e}_r \biggl[ F_r \frac{\partial B_r}{\partial r} + \frac{F_\theta}{r} \frac{\partial B_r}{\partial \theta} + \cancel{ \frac{F_\varphi}{r\sin\theta} \frac{\partial B_r}{\partial \varphi} } - \frac{(F_\theta B_\theta + F_\varphi B_\varphi)}{r}\biggr]
</math>
</td>
</tr>

<tr>
<td align="right">
 
</td>
<td align="center">
 
</td>
<td align="left">
<math>
+ \hat{e}_\theta \biggl[ F_r \frac{\partial B_\theta}{\partial r} + \frac{F_\theta}{r} \frac{\partial B_\theta}{\partial \theta } + \cancel{ \frac{F_\varphi}{r\sin\theta} \frac{\partial B_\theta}{\partial \varphi} }
+ \frac{F_\theta B_r}{r} - \frac{F_\varphi B_\varphi \cot\theta}{r} \biggr]
</math>
</td>
</tr>

<tr>
<td align="right">
 
</td>
<td align="center">
 
</td>
<td align="left">
<math>
+ \hat{e}_\varphi \biggl[ F_r \frac{\partial B_\varphi}{\partial r} + \frac{F_\theta}{r} \frac{\partial B_\varphi}{\partial \theta} + \cancel{ \frac{F_\varphi}{r\sin\theta} \frac{\partial B_\varphi}{\partial \varphi} }
+ \frac{F_\varphi B_r}{r} + \frac{F_\varphi B_\theta \cot\theta}{r} \biggr] \, .
</math>
</td>
</tr>
<tr>
<td align="center" colspan="3">
[<b>[[User:Tohline/Appendix/References#BT87|<font color="red">BT87</font>]]</b>], p. 651, Eq. (1B-55)
</td>
</tr>
</table>
From this last expression — the so-called ''convective operator'' — we conclude as well that, for axisymmetric systems,

<table border="0" cellpadding="5" align="center">

<tr>
<td align="right">
<math>
(\vec{v} \cdot \nabla )\vec{v}
</math>
</td>
<td align="center">
=
</td>
<td align="left">
<math>
\hat{e}_r \biggl[ v_r \frac{\partial v_r}{\partial r} + \frac{v_\theta}{r} \frac{\partial v_r}{\partial \theta} - \frac{(v_\theta^2 + v_\varphi^2 )}{r}\biggr]
+ \hat{e}_\theta \biggl[ v_r \frac{\partial v_\theta}{\partial r} + \frac{v_\theta}{r} \frac{\partial v_\theta}{\partial \theta }
+ \frac{v_\theta v_r}{r} - \frac{v_\varphi^2 \cot\theta}{r} \biggr]
+ \hat{e}_\varphi \biggl[ v_r \frac{\partial v_\varphi}{\partial r} + \frac{v_\theta}{r} \frac{\partial v_\varphi}{\partial \theta}
+ \frac{v_\varphi v_r}{r} + \frac{v_\varphi v_\theta \cot\theta}{r} \biggr] \, .
</math>
</td>
</tr>
</table>

<li>Express all vector time-derivatives in spherical coordinates:

<table align="center" border="0" cellpadding="5">
<tr>
<td colspan="3" align="center">
<font color="#770000"><b>Vector Time-Derivatives in Spherical Coordinates</b></font>
</td>
</tr>

<tr>
<td align="right">
<math>
\frac{d}{dt}\vec{F}
</math>
</td>
<td align="center">
=
</td>
<td align="left">
<math>
{\hat{e}}_r \frac{dF_r}{dt} + F_r \frac{d{\hat{e}}_r}{dt} + {\hat{e}}_\theta \frac{dF_\theta}{dt} + F_\theta \frac{d{\hat{e}}_\theta}{dt} + {\hat{e}}_\varphi \frac{dF_\varphi}{dt} + F_\varphi \frac{d{\hat{e}}_\varphi}{dt}
</math>
</td>
</tr>

<tr>
<td align="right">
 
</td>
<td align="center">
=
</td>
<td align="left">
<math>
{\hat{e}}_r \frac{dF_r}{dt} + F_r \biggl[ {\hat{e}}_\theta \dot\theta + {\hat{e}}_\varphi \dot\varphi \sin\theta \biggr]
+ {\hat{e}}_\theta \frac{dF_\theta}{dt} + F_\theta \biggl[ - {\hat{e}}_r \dot\theta + {\hat{e}}_\varphi \dot\varphi \cos\theta \biggr]
+ {\hat{e}}_\varphi \frac{dF_\varphi}{dt} + F_\varphi \biggl[ - {\hat{e}}_r \dot\varphi \sin\theta - {\hat{e}}_\theta \dot\varphi \cos\theta \biggr]
</math>
</td>
</tr>

<tr>
<td align="right">
 
</td>
<td align="center">
=
</td>
<td align="left">
<math>
{\hat{e}}_r \biggl[ \frac{dF_r}{dt} - F_\theta \dot\theta - F_\varphi \dot\varphi \sin\theta \biggr]
+ {\hat{e}}_\theta \biggl[ \frac{dF_\theta}{dt} + F_r \dot\theta - F_\varphi \dot\varphi \cos\theta \biggr]
+ {\hat{e}}_\varphi \biggl[ \frac{dF_\varphi}{dt} + F_r \dot\varphi \sin\theta + F_\theta \dot\varphi \cos\theta \biggr] ;
</math>
</td>
</tr>

<tr>
<td align="right">
<math>
\vec{v} = \frac{d\vec{x}}{dt}
</math>
</td>
<td align="center">
=
</td>
<td align="left">
<math>
\frac{d}{dt}\biggl[ \hat{e}_r r \biggr]
= {\hat{e}}_r \dot{r}
+ {\hat{e}}_\theta~ r \dot\theta
+ {\hat{e}}_\varphi ~r \sin\theta ~ \dot\varphi .
</math>
</td>
</tr>
<tr>
<td align="center" colspan="3">
[<b>[[User:Tohline/Appendix/References#BT87|<font color="red">BT87</font>]]</b>], p. 648, Eq. (1B-30)
</td>
</tr>
</table>

</ol>
===Governing Equations (SPH.)===

Introducing the above expressions into the [[User:Tohline/PGE|principal governing equations]] gives,

<div align="center">
<span id="Continuity"><font color="#770000">'''Equation of Continuity'''</font></span><br />

<table border="0" cellpadding="5" align="center">

<tr>
<td align="right">
<math>~\frac{d\rho}{dt} + \rho \biggl[ \frac{1}{r^2} \frac{\partial (r^2 \dot{r})}{\partial r}
+ \frac{1}{r\sin\theta} \frac{\partial }{\partial\theta} \biggl( \dot\theta r \sin\theta \biggr)
\biggr]</math>
</td>
<td align="center">
<math>~=</math>
</td>
<td align="left">
<math>~0</math>
</td>
</tr>
</table>

<span id="PGE:Euler">
<font color="#770000">'''Euler Equation'''</font>
</span><br />

<table border="0" cellpadding="5" align="center">


<tr>
<td align="right">
<math>~
{\hat{e}}_r \biggl[ \frac{d\dot{r}}{dt} - r {\dot\theta}^2 - r {\dot\varphi}^2 \sin^2\theta \biggr]
+ {\hat{e}}_\theta \biggl[ \frac{d(r\dot\theta)}{dt} + \dot{r} \dot\theta - r { \dot\varphi }^2 \sin\theta \cos\theta \biggr]
+ {\hat{e}}_\varphi \biggl[ \frac{d(r \sin\theta \dot\varphi)}{dt} + \dot{r} \dot\varphi \sin\theta + r \dot\theta \dot\varphi \cos\theta \biggr]
</math>
</td>
<td align="center">
<math>~=</math>
</td>
<td align="left">
<math>~- {\hat{e}}_r \biggl[ \frac{1}{\rho} \frac{\partial P}{\partial r}+ \frac{\partial \Phi }{\partial r} \biggr]
- {\hat{e}}_\theta \biggl[ \frac{1}{\rho r} \frac{\partial P}{\partial\theta} + \frac{1}{r} \frac{\partial \Phi}{\partial\theta} \biggr]
</math>
</td>
</tr>
</table>

<span id="PGE:AdiabaticFirstLaw">Adiabatic Form of the<br />
<font color="#770000">'''First Law of Thermodynamics'''</font></span><br />

{{User:Tohline/Math/EQ_FirstLaw02}}

<span id="PGE:Poisson"><font color="#770000">'''Poisson Equation'''</font></span><br />

<table border="0" cellpadding="5" align="center">

<tr>
<td align="right">
<math>~
\frac{1}{r^2} \frac{\partial }{\partial r} \biggl[ r^2 \frac{\partial \Phi }{\partial r} \biggr]
+ \frac{1}{r^2 \sin\theta} \frac{\partial }{\partial \theta}\biggl(\sin\theta ~ \frac{\partial \Phi}{\partial\theta}\biggr)
</math>
</td>
<td align="center">
<math>~=</math>
</td>
<td align="left">
<math>~4\pi G\rho</math>
</td>
</tr>
</table>

</div>

===Conservation of Specific Angular Momentum (SPH.)===

The <math>\hat{e}_\varphi</math> component of the Euler equation leads to a statement of conservation of specific angular momentum, <math>~j</math>, as follows.

<div align="center">

<table border="0" cellpadding="5" align="center">

<tr>
<td align="right">
<math>~0</math>
</td>
<td align="center">
<math>~=</math>
</td>
<td align="left">
<math>~
\frac{d(r \sin\theta \dot\varphi)}{dt} + \dot{r} \dot\varphi \sin\theta + r \dot\theta \dot\varphi \cos\theta
</math>
</td>
</tr>

<tr>
<td align="right">
 
</td>
<td align="center">
<math>~=</math>
</td>
<td align="left">
<math>~
\frac{1}{r\sin\theta} \biggl[ r\sin\theta\frac{d(r \sin\theta \dot\varphi)}{dt} + r\sin\theta \dot\varphi ( \dot{r}\sin\theta + r\dot\theta \cos\theta) \biggr]
</math>
</td>
</tr>

<tr>
<td align="right">
 
</td>
<td align="center">
<math>~=</math>
</td>
<td align="left">
<math>~
\frac{1}{r\sin\theta} \biggl[\frac{d(r^2 \sin^2\theta \dot\varphi )}{dt} \biggr] \, .
</math>
</td>
</tr>
</table>

<math>
\Rightarrow ~~~~~ j(r,\theta) \equiv (r\sin\theta)^2 \dot\varphi = \mathrm{constant} ~(\mathrm{i.e.,}~\mathrm{independent~of~time})
</math><br />

</div>

<span id="RelevantSphericalComponents">So, for axisymmetric configurations, the <math>\hat{e}_r</math> and <math>\hat{e}_\theta</math> components of the Euler equation become, respectively,</span>

<table border="1" align="center" cellpadding="10"><tr><td align="center">
<table border="0" cellpadding="5" align="center">
<tr>
<td align="right"><math>~{\hat{e}}_r</math>:    </td>
<td align="right">
<math>
\frac{d\dot{r}}{dt} - r {\dot\theta}^2 - \biggl[ \frac{j^2}{r^3 \sin^3\theta} \biggr]\sin\theta
</math>
</td>
<td align="center">
=
</td>
<td align="left">
<math>
- \biggl[ \frac{1}{\rho} \frac{\partial P}{\partial r}+ \frac{\partial \Phi }{\partial r} \biggr] \, ,
</math>
</td>
</tr>

<tr>
<td align="right"><math>~{\hat{e}}_\theta</math>:    </td>
<td align="right">
<math>
\frac{d(r\dot\theta)}{dt} + \dot{r} \dot\theta - \biggl[ \frac{j^2}{r^3 \sin^3\theta} \biggr] \cos\theta
</math>
</td>
<td align="center">
=
</td>
<td align="left">
<math>
- \biggl[ \frac{1}{\rho r} \frac{\partial P}{\partial\theta} + \frac{1}{r} \frac{\partial \Phi}{\partial\theta} \biggr] \, .
</math>
</td>
</tr>
</table>
</td></tr></table>

===Eulerian Formulation (SPH.)===

Each of the above simplified governing equations has been written in terms of Lagrangian time derivatives. An Eulerian formulation of each equation can be obtained by replacing each Lagrangian time derivative by its Eulerian counterpart. Specifically, for any scalar function, <math>f</math>,

<div align="center">
<math>
\frac{df}{dt} \rightarrow \frac{\partial f}{\partial t} + (\vec{v}\cdot \nabla)f =
\frac{\partial f}{\partial t} +
\biggl[ v_r \frac{\partial f}{\partial r} \biggr]
+ \biggl[ \frac{v_\theta}{r} \frac{\partial f}{\partial\theta} \biggr] =
\frac{\partial f}{\partial t} +
\biggl[ \dot{r} \frac{\partial f}{\partial r} \biggr]
+ \biggl[ \dot\theta \frac{\partial f}{\partial\theta} \biggr]
\, .
</math>
</div>

Making this substitution throughout the set of governing relations gives:

<div align="center">
<span id="Continuity"><font color="#770000">'''Equation of Continuity'''</font></span><br />

<table border="0" cellpadding="5" align="center">

<tr>
<td align="right">
<math>~
\frac{\partial \rho}{\partial t}
+ \biggl[ \frac{1}{r^2} \frac{\partial (\rho r^2 \dot{r})}{\partial r}
+ \frac{1}{r\sin\theta} \frac{\partial }{\partial\theta} \biggl( \rho \dot\theta r \sin\theta \biggr)
\biggr]</math>
</td>
<td align="center">
<math>~=</math>
</td>
<td align="left">
<math>~0</math>
</td>
</tr>
</table>

<span id="PGE:Euler">The Two Relevant Components of the<br />
<font color="#770000">'''Euler Equation'''</font>
</span><br />

<table border="0" cellpadding="5" align="center">
<tr>
<td align="right"><math>~{\hat{e}}_r</math>:    </td>
<td align="right">
<math>
\biggl\{ \frac{\partial \dot{r}}{\partial t} + \biggl[ \dot{r} \frac{\partial \dot{r}}{\partial r} \biggr] + \biggl[ \dot\theta \frac{\partial \dot{r}}{\partial\theta} \biggr] \biggr\}
- r {\dot\theta}^2
</math>
</td>
<td align="center">
=
</td>
<td align="left">
<math>
- \biggl[ \frac{1}{\rho} \frac{\partial P}{\partial r}+ \frac{\partial \Phi }{\partial r} \biggr] + \biggl[ \frac{j^2}{r^3 \sin^2\theta} \biggr]
</math>
</td>
</tr>

<tr>
<td align="right"><math>~{\hat{e}}_\theta</math>:    </td>
<td align="right">
<math>
r \biggl\{ \frac{\partial \dot{\theta}}{\partial t} + \biggl[ \dot{r} \frac{\partial \dot{\theta}}{\partial r} \biggr] + \biggl[ \dot\theta \frac{\partial \dot{\theta}}{\partial\theta} \biggr] \biggr\} + 2\dot{r} \dot\theta
</math>
</td>
<td align="center">
=
</td>
<td align="left">
<math>
- \biggl[ \frac{1}{\rho r} \frac{\partial P}{\partial\theta} + \frac{1}{r} \frac{\partial \Phi}{\partial\theta} \biggr] + \biggl[ \frac{j^2}{r^3 \sin^3\theta} \biggr] \cos\theta
</math>
</td>
</tr>
</table>

<span id="PGE:AdiabaticFirstLaw">Adiabatic Form of the<br />
<font color="#770000">'''First Law of Thermodynamics'''</font></span><br />

<table border="0" cellpadding="5" align="center">

<tr>
<td align="right">
<math>~
\biggl\{ \frac{\partial \epsilon}{\partial t} + \biggl[ \dot{r} \frac{\partial \epsilon}{\partial r} \biggr] + \biggl[ \dot\theta \frac{\partial \epsilon}{\partial\theta} \biggr] \biggr\}
+ P\biggl\{ \frac{\partial }{\partial t} \biggl( \frac{1}{\rho}\biggr)
+ \biggl[ \dot{r} \frac{\partial }{\partial r} \biggl( \frac{1}{\rho}\biggr) \biggr]
+ \biggl[ \dot\theta \frac{\partial }{\partial\theta} \biggl( \frac{1}{\rho}\biggr) \biggr] \biggr\}
</math>
</td>
<td align="center">
<math>~=</math>
</td>
<td align="left">
<math>~0</math>
</td>
</tr>
</table>

<span id="PGE:Poisson"><font color="#770000">'''Poisson Equation'''</font></span><br />

<table border="0" cellpadding="5" align="center">

<tr>
<td align="right">
<math>~
\frac{1}{r^2} \frac{\partial }{\partial r} \biggl[ r^2 \frac{\partial \Phi }{\partial r} \biggr]
+ \frac{1}{r^2 \sin\theta} \frac{\partial }{\partial \theta}\biggl(\sin\theta ~ \frac{\partial \Phi}{\partial\theta}\biggr)
</math>
</td>
<td align="center">
<math>~=</math>
</td>
<td align="left">
<math>~4\pi G\rho</math>
</td>
</tr>
</table>

</div>

=See Also=

{{ SGFfooter }}

AxisymmetricConfigurations/PGE - Revision history

Jet53man: /* Axisymmetric Configurations (Governing Equations) */

Jet53man at 21:07, 9 September 2021

Jet53man: Created page with "__FORCETOC__ =Axisymmetric Configurations (Governing Equations)= {| class="AxisymmetricConfigurationsPGE" style="float:left; margin-righ..."

Jet53man: Created page with "FORCETOC =Axisymmetric Configurations (Governing Equations)= {| class="AxisymmetricConfigurationsPGE" style="float:left; margin-righ..."