Difference between revisions of "User:Tohline/Apps/GoldreichWeber80"
(Introduce equations (1)-(4) of Goldreich & Weber (1980)) |
(→Introduction: Identify constraints) |
||
Line 2: | Line 2: | ||
{{LSU_HBook_header}} | {{LSU_HBook_header}} | ||
==Introduction== | ==Introduction== | ||
This is principally a review of the dynamical model that Peter Goldreich & Stephen Weber [http://adsabs.harvard.edu/abs/1980ApJ...238..991G (1980, ApJ, 238, 991)] developed to describe the near-homologous collapse of stellar cores. Goldreich & Weber begin with the identical set of [[User:Tohline/PGE#Principal_Governing_Equations|principal governing equations]] that serves as the foundation for all of the discussions throughout this H_Book. In particular, as is documented by their equation (1), their adopted equation of state is adiabatic/polytropic, | This is principally a review of the dynamical model that Peter Goldreich & Stephen Weber [http://adsabs.harvard.edu/abs/1980ApJ...238..991G (1980, ApJ, 238, 991)] developed to describe the near-homologous collapse of stellar cores. | ||
===Governing Equations=== | |||
Goldreich & Weber begin with the identical set of [[User:Tohline/PGE#Principal_Governing_Equations|principal governing equations]] that serves as the foundation for all of the discussions throughout this H_Book. In particular, as is documented by their equation (1), their adopted equation of state is adiabatic/polytropic, | |||
<div align="center"> | <div align="center"> | ||
<math>~P = \kappa \rho^\gamma \, ,</math> | <math>~P = \kappa \rho^\gamma \, ,</math> | ||
Line 42: | Line 45: | ||
{{User:Tohline/Math/EQ_Poisson01}} | {{User:Tohline/Math/EQ_Poisson01}} | ||
</div> | </div> | ||
Tweaking the set of principal governing equations, as we have written them, to even more precisely match equations (1) - (4) in Goldreich & Weber (1980), we should replace the state variable <math>~P</math> (pressure) with <math>~H</math> ( | Tweaking the set of principal governing equations, as we have written them, to even more precisely match equations (1) - (4) in Goldreich & Weber (1980), we should replace the state variable <math>~P</math> (pressure) with <math>~H</math> (enthalpy), keeping in mind that, <math>~\gamma = 1 + 1/n</math>, and, [[User:Tohline/SR#Barotropic_Structure|as presented in our introductory discussion of barotropic supplemental relations]], | ||
<div align="center"> | <div align="center"> | ||
<math>~H = \biggl( \frac{\gamma}{\gamma-1} \biggr) \kappa \rho^{\gamma-1} \, ,</math> | <math>~H = \biggl( \frac{\gamma}{\gamma-1} \biggr) \kappa \rho^{\gamma-1} \, ,</math> | ||
Line 50: | Line 53: | ||
<math>~\nabla H = \frac{\nabla P}{\rho} \, .</math> | <math>~\nabla H = \frac{\nabla P}{\rho} \, .</math> | ||
</div> | </div> | ||
===Imposed Constraints=== | |||
[http://adsabs.harvard.edu/abs/1980ApJ...238..991G Goldreich & Weber (1980)] specifically choose to examine the spherically symmetric collapse of a <math>~\gamma = 4/3</math> fluid. With this choice of adiabatic index, the equation of state becomes, | |||
<div align="center"> | |||
<math>~H = 4 \kappa \rho^{1/3} \, .</math> | |||
</div> | |||
And because a strictly radial flow-field exhibits no vorticity (i.e., <math>\vec\zeta = 0</math>), the Euler equation can be rewritten as, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\frac{\partial v_r}{\partial t} </math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~-~ \nabla_r \biggl[ H + \Phi + \frac{1}{2}v^2 \biggr] \, .</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
Goldreich & Weber also realize that, because the flow is vorticity free, the velocity can be obtained from a stream function, <math>~\Psi</math>, via the relation, | |||
<div align="center"> | |||
<math>~\vec{v} = \nabla\Psi \, .</math> | |||
</div> | |||
Hence, the Euler equation becomes, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\frac{\partial \Psi}{\partial t} </math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~-~ \biggl[ H + \Phi + \frac{1}{2}\biggl( \nabla \Psi \biggr)^2 \biggr] \, ,</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
where, <math>~H</math>, <math>~\Phi</math>, and <math>~\Psi</math> are each functions only of the radial coordinate. | |||
=Related Discussions= | =Related Discussions= |
Revision as of 01:50, 1 September 2014
Homologously Collapsing Stellar Cores
| Tiled Menu | Tables of Content | Banner Video | Tohline Home Page | |
Introduction
This is principally a review of the dynamical model that Peter Goldreich & Stephen Weber (1980, ApJ, 238, 991) developed to describe the near-homologous collapse of stellar cores.
Governing Equations
Goldreich & Weber begin with the identical set of principal governing equations that serves as the foundation for all of the discussions throughout this H_Book. In particular, as is documented by their equation (1), their adopted equation of state is adiabatic/polytropic,
<math>~P = \kappa \rho^\gamma \, ,</math>
— where both <math>~\kappa</math> and <math>~\gamma</math> are constants — and therefore satisfies what we have referred to as the
Adiabatic Form of the
First Law of Thermodynamics
(Specific Entropy Conservation)
<math>~\frac{d\epsilon}{dt} + P \frac{d}{dt} \biggl(\frac{1}{\rho}\biggr) = 0</math> .
their equation (2) is what we have referred to as the
Eulerian Representation
or
Conservative Form
of the Continuity Equation,
<math>~\frac{\partial\rho}{\partial t} + \nabla \cdot (\rho \vec{v}) = 0</math>
their equation (3) is what we have referred to as the
Euler Equation
in terms of the Vorticity,
<math>~\frac{\partial\vec{v}}{\partial t} + \vec\zeta \times \vec{v}= - \frac{1}{\rho} \nabla P - \nabla \biggl[\Phi + \frac{1}{2}v^2 \biggr] </math>
where, <math>~\vec\zeta \equiv \nabla\times \vec{v}</math> is the fluid vorticity; and their equation (4) is the
Tweaking the set of principal governing equations, as we have written them, to even more precisely match equations (1) - (4) in Goldreich & Weber (1980), we should replace the state variable <math>~P</math> (pressure) with <math>~H</math> (enthalpy), keeping in mind that, <math>~\gamma = 1 + 1/n</math>, and, as presented in our introductory discussion of barotropic supplemental relations,
<math>~H = \biggl( \frac{\gamma}{\gamma-1} \biggr) \kappa \rho^{\gamma-1} \, ,</math>
and,
<math>~\nabla H = \frac{\nabla P}{\rho} \, .</math>
Imposed Constraints
Goldreich & Weber (1980) specifically choose to examine the spherically symmetric collapse of a <math>~\gamma = 4/3</math> fluid. With this choice of adiabatic index, the equation of state becomes,
<math>~H = 4 \kappa \rho^{1/3} \, .</math>
And because a strictly radial flow-field exhibits no vorticity (i.e., <math>\vec\zeta = 0</math>), the Euler equation can be rewritten as,
<math>~\frac{\partial v_r}{\partial t} </math> |
<math>~=</math> |
<math>~-~ \nabla_r \biggl[ H + \Phi + \frac{1}{2}v^2 \biggr] \, .</math> |
Goldreich & Weber also realize that, because the flow is vorticity free, the velocity can be obtained from a stream function, <math>~\Psi</math>, via the relation,
<math>~\vec{v} = \nabla\Psi \, .</math>
Hence, the Euler equation becomes,
<math>~\frac{\partial \Psi}{\partial t} </math> |
<math>~=</math> |
<math>~-~ \biggl[ H + \Phi + \frac{1}{2}\biggl( \nabla \Psi \biggr)^2 \biggr] \, ,</math> |
where, <math>~H</math>, <math>~\Phi</math>, and <math>~\Psi</math> are each functions only of the radial coordinate.
Related Discussions
© 2014 - 2021 by Joel E. Tohline |