Difference between revisions of "User:Tohline/2DStructure/AxisymmetricInstabilities"

From VistrailsWiki
Jump to navigation Jump to search
 
(34 intermediate revisions by the same user not shown)
Line 2: Line 2:
<!-- __NOTOC__ will force TOC off -->
<!-- __NOTOC__ will force TOC off -->
=Axisymmetric Instabilities to Avoid=
=Axisymmetric Instabilities to Avoid=
Here we draw heavily from the extensive discussion of instabilities that appears in [<b>[[User:Tohline/Appendix/References#T78|<font color="red">T78</font>]]</b>]




{{LSU_HBook_header}}
{{LSU_HBook_header}}


When constructing rotating equilibrium configurations that obey a [[User:Tohline/SR#Time-Independent_Problems|barotropic equation of state]], keep in mind that certain parameter ranges should be avoided because they will lead to structures that are unstable toward the dynamical development of local, ''convective''-type motions.  Here are a few well-known examples.
When constructing rotating equilibrium configurations that obey a [[User:Tohline/SR#Time-Independent_Problems|barotropic equation of state]], keep in mind that certain physical variable profiles should be avoided because they will lead to structures that are unstable toward the dynamical development of shape-distorting or ''convective''-type motions.  Here are a few well-known examples.


==Rayleigh-Taylor Instability==
==Rayleigh-Taylor Instability==
This bouyancy-driven instability arises when the condition,
Referencing both p. 101 of [ <b>[[User:Tohline/Appendix/References#Shu92|<font color="red">Shu92</font>]] </b>], and volume I, p. 410 of [ <b>[[User:Tohline/Appendix/References#P00|<font color="red">P00</font>]] </b>], a [https://en.wikipedia.org/wiki/Rayleigh%E2%80%93Taylor_instability Rayleigh-Taylor] instability is a bouyancy-driven instability that arises when a heavy fluid rests on top of a light fluid in an effective gravitational field, <math>~\vec{g}</math>.  In the simplest case of spherically symmetric, self-gravitating configurations, the condition for stability against a Rayleigh-Taylor instability may be written as,
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 28: Line 27:
</table>
</table>
</div>
</div>
is violated, where, <math>~\vec{g}</math> is the local gravitational acceleration vectorIn the simplest case of spherically symmetric configurations, this means that the mass density ''must decrease outward.''
that is to say, the mass density ''must decrease outward.''  In an expanded discussion, [ <b>[[User:Tohline/Appendix/References#P00|<font color="red">P00</font>]] </b>] &#8212; see pp. 410 - 413 &#8212; derives a dispersion relation that describes the development of the Rayleigh-Taylor instability in a plane-parallel fluid layer that initially contains a discontinuous jump/drop in the density.
 
In addition, [ <b>[[User:Tohline/Appendix/References#P00|<font color="red">P00</font>]] </b>] &#8212; see pp. 413 - 416 &#8212; and [ <b>[[User:Tohline/Appendix/References#Shu92|<font color="red">Shu92</font>]] </b>] &#8212; see pp. 101 - 105 &#8212; both discuss circumstances that can give rise to the so-called ''Kelvin-Helmholtz'' instability.  It is another common two-fluid instability, but one that depends on the existence of transverse velocity flows and that is independent of the gravitational field.
 
==Poincar&eacute;-Wavre Theorem==
As [ <b>[[User:Tohline/Appendix/References#T78|<font color="red">T78</font>]] </b>] points out &#8212; see his pp. 78 - 81 &#8212; Poincar&eacute; and Wavre were the first to, effectively, prove the following theorem:
<table border="0" align="center" width="75%" cellpadding="5">
<tr><td align="left">For rotating, self-gravitating configurations<font color="darkgreen"> "any of the following statements implies the three others: &nbsp; (i) the angular velocity is a constant over cylinders centered about the axis of rotation, (ii) the effective gravity can be derived from a potential, (iii) the effective gravity is normal to the isopycnic surfaces, (iv) the isobaric- and isopycnic-surfaces coincide."
</font></td></tr>
</table>
Among other things, this implies that for rotating barotropic configurations not only is the equation of state given by a function of the form, <math>~P = P(\rho)</math>, but it must also be true that,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{\partial \dot\varphi}{\partial z}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~0 \, .</math>
  </td>
</tr>
<tr>
  <td align="center" colspan="3">
[ <b>[[User:Tohline/Appendix/References#T78|<font color="red">T78</font>]] </b>], &sect;4.3, Eq. (30)
  </td>
</tr>
</table>
</div>
 
NOTE: We should investigate how this theorem comes into play in the context of our [[User:Tohline/ThreeDimensionalConfigurations/RiemannTypeI#Riemann_Type_1_Ellipsoids|accompanying discussion of Type 1 Riemann ellipsoids]].  These are equilibrium triaxial, uniform-density configurations in which the system's internal vorticity vector does not align with the tumble-axis of the ellipsoid and, therefore apparently <math>~\dot\varphi</math> is not independent of <math>~z</math>.


==H&oslash;iland Criterion==
==H&oslash;iland Criterion==


According to [ <b>[[User:Tohline/Appendix/References#T78|<font color="red">T78</font>]] </b>] &#8212; see p. 168 &#8212; this pair of mathematically expressed conditions has the following meaning:
<table border="0" cellpadding="3" align="center" width="60%">
<table border="0" align="center" width="75%" cellpadding="8">
<tr><td align="left">
<font color="indianred">
"For an incompressible liquid contained between concentric cylinders and rotating with angular velocity <math>~\Omega(\varpi)</math>, Rayleigh's criterion &#8212; that <math>~\varpi^4 \Omega^2</math> increase outwards &#8212; is necessary and sufficient for stability to axisymmetric disturbances.  In a star rotating with angular velocity <math>~\Omega(\varpi)</math> if we continue to restrict attention to axisymmetric disturbances, this criterion must be modified by buoyancy effects; that is, some combination of the Rayleigh and Schwarzschild criteria should obtain.  Such a combination has been found, for example, by [https://en.wikipedia.org/wiki/Einar_Høiland H&oslash;iland (1941)] (see also [https://archive.org/details/AllerStellarStructure Ledoux's Chapter 10, pp. 499-574 of ''Stellar Structure'' (1965)], and indicates, as one would expect, that a stable stratification of angular velocity exerts a stabilizing influence on an unstable distribution of temperature, and vice versa.  The combined criterion has not been placed on the solid analytical foundation of its two component criteria, however."</font>
</td></tr>
<tr><td align="right">
&#8212; Drawn from p. 475 of [https://ui.adsabs.harvard.edu/abs/1967ARA%26A...5..465L/abstract N. R. Lebovitz (1967)], ARAA, 5, 465
</td></tr></table>
 
As is stated on p. 166 of [ <b>[[User:Tohline/Appendix/References#T78|<font color="red">T78</font>]] </b>], in rotating barotropic configurations, axisymmetric stability requires the simultaneous satisfaction of the following pair of conditions:
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\biggl(\frac{1}{\varpi^3} \biggr) \frac{\partial j^2}{\partial \varpi} + \frac{1}{c_P} \biggl( \frac{\gamma - 1}{\Gamma_3 - 1}\biggr)
(- \vec{g} ) \cdot \nabla s</math>
  </td>
  <td align="center">
<math>~></math>
  </td>
  <td align="left">
<math>~0 </math>
  </td>
  <td align="center">&nbsp; &nbsp;&nbsp;&nbsp;[stable] ;</td>
</tr>
<tr>
  <td align="center" colspan="4">
[ <b>[[User:Tohline/Appendix/References#T78|<font color="red">T78</font>]] </b>], &sect;7.3, Eq. (41)
<br />see also<br />[ <b>[[User:Tohline/Appendix/References#KW94|<font color="red">KW94</font>]] </b>], &sect;43.2, Eq. (43.22)
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~-g_z \biggl[ \frac{\partial j^2}{\partial \varpi} \biggl(\frac{\partial s}{\partial z} \biggr)
- \frac{\partial j^2}{\partial z} \biggl(\frac{\partial s}{\partial \varpi} \biggr)\biggr]</math>
  </td>
  <td align="center">
<math>~></math>
  </td>
  <td align="left">
<math>~0 </math>
  </td>
  <td align="center">&nbsp; &nbsp;&nbsp;&nbsp;[stable] .</td>
</tr>
<tr>
  <td align="center" colspan="4">
[ <b>[[User:Tohline/Appendix/References#T78|<font color="red">T78</font>]] </b>], &sect;7.3, Eq. (42)
<br />see also<br />[ <b>[[User:Tohline/Appendix/References#KW94|<font color="red">KW94</font>]] </b>], &sect;43.2, Eq. (43.23)
  </td>
</tr>
</table>
</div>
 
where, <math>~s</math>, is the local specific entropy, and <math>~j \equiv \dot\varphi \varpi^2</math>, is the local specific angular momentum of the fluid.  According to [ <b>[[User:Tohline/Appendix/References#T78|<font color="red">T78</font>]] </b>] &#8212; see p. 168 &#8212; this pair of mathematically expressed conditions has the following meaning:
<table border="0" align="center" width="75%" cellpadding="5">
<tr><td align="left"><font color="darkgreen">
<tr><td align="left"><font color="darkgreen">
"A baroclinic star in permanent rotation is dynamically stable with respect to axisymmetric motions if and only if the two following conditions are satisfied:  &nbsp; (i) the entropy per unit mass, <math>~s</math>, never decreases outward, and (ii) on each surface <math>~s</math> = constant, the angular momentum per unit mass, <math>~\Omega \varpi^2</math>, increases as we move from the poles to the equator."
"A baroclinic star in permanent rotation is dynamically stable with respect to axisymmetric motions if and only if the two following conditions are satisfied:  &nbsp; (i) the entropy per unit mass, <math>~s</math>, never decreases outward, and (ii) on each surface <math>~s</math> = constant, the angular momentum per unit mass, <math>~j</math>, increases as we move from the poles to the equator."
</font></td></tr>
</font></td></tr>
</table>
</table>
===Schwarzschild Criterion===
===Schwarzschild Criterion===
In the case of nonrotating equilibrium configurations, the H&oslash;iland Criterion reduces to the Schwarzschild criterion.  That is, thermal convection arises when the condition,
In the case of nonrotating equilibrium configurations, the H&oslash;iland Criterion reduces to the Schwarzschild criterion.  That is, thermal convection arises when the condition,
Line 47: Line 135:
<math>~(- \vec{g} ) \cdot \nabla s</math>
<math>~(- \vec{g} ) \cdot \nabla s</math>
   </td>
   </td>
   <td align="center">
   <td align="center" width="20px">
<math>~> </math>
  </td>
  <td align="left">
<math>~0 </math>
  </td>
  <td align="left">&nbsp; &nbsp;&nbsp;&nbsp;[stable] ,</td>
</tr>
<tr>
  <td align="center" colspan="4">
[ <b>[[User:Tohline/Appendix/References#T78|<font color="red">T78</font>]] </b>], &sect;7.3, Eq. (43)
<br />see also<br />[ <b>[[User:Tohline/Appendix/References#KW94|<font color="red">KW94</font>]] </b>], &sect;6.1, Eq. (6.13) &hellip; or &hellip; pp. 93 - 98 of [ <b>[[User:Tohline/Appendix/References#Shu92|<font color="red">Shu92</font>]] </b>]
  </td>
</tr>
</table>
</div>
is violated.  This means that, in order for a spherical system to be stable against dynamical convective motions, the specific entropy ''must increase outward.''
 
===Solberg/Rayleigh Criterion===
In the case of an homentropic equilibrium configuration, the H&oslash;iland Criterion reduces to the Solberg criterion.  That is, an axisymmetric exchange of fluid "rings" will occur on a dynamical time scale if the condition,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{dj^2}{d\varpi} </math>
  </td>
  <td align="center" width="20px">
<math>~> </math>
<math>~> </math>
   </td>
   </td>
Line 53: Line 168:
<math>~0 </math>
<math>~0 </math>
   </td>
   </td>
   <td align="center">&nbsp; &nbsp;&nbsp;&nbsp;[stable] ,</td>
   <td align="left">&nbsp; &nbsp;&nbsp;&nbsp;[stable] ,</td>
</tr>
</tr>
<tr>
<tr>
   <td align="center" colspan="4">
   <td align="center" colspan="4">
[ <b>[[User:Tohline/Appendix/References#T78|<font color="red">T78</font>]] </b>], p. 167, Eq. (43)
[ <b>[[User:Tohline/Appendix/References#T78|<font color="red">T78</font>]] </b>], &sect;7.3, Eq. (44)
<br />see also<br />[ <b>[[User:Tohline/Appendix/References#KW94|<font color="red">KW94</font>]] </b>], &sect;43.2, Eq. (43.18) &hellip; or &hellip; pp. 98 - 101 of [ <b>[[User:Tohline/Appendix/References#Shu92|<font color="red">Shu92</font>]] </b>]
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>
</div>
is violated, where, <math>~s</math> is the local specific entropy of the fluid.  This means that, in order for a spherical system to be stable against dynamical convective motions, the specific entropy ''must increase outward.''
is violated.  This means that, for stability, the specific angular momentum ''must necessarily increase outward.'' As [ <b>[[User:Tohline/Appendix/References#T78|<font color="red">T78</font>]] </b>] points out, this <font color="darkgreen">"Solberg criterion generalizes to homentropic bodies the well-known [http://adsabs.harvard.edu/abs/1917RSPSA..93..148R Rayleigh (1917)] criterion for an inviscid, incompressible fluid."</font>


==Poincar&eacute;-Wavre Theorem==
=Modeling Implications and Advice=
Here are some recommendations to keep in mind as you attempt to construct equilibrium models of self-gravitating astrophysical fluids.
 
<font color="purple">Rayleigh-Taylor instability:</font>&nbsp; In order to avoid constructing configurations that are subject to the Rayleigh-Taylor instability, be sure that lower density material is never placed "beneath" higher density material.  For example &hellip;
<ul><li>When building a bipolytropic configuration, a value for the (imposed) discontinuous jump in the mean-molecular weight will need to be specified at the interface between the envelope and the core.  If you choose a ratio, <math>~{\bar\mu}_e/{\bar\mu}_c</math>, that is greater than unity, the resulting equilibrium model will exhibit a discontinuous density jump that makes the density higher at the base of the envelope than it is at the surface of the core.  The core/envelope interface of this configuration will be unstable to the Rayleigh-Taylor instability, but you won't know that until and unless you examine the hydrodynamic stability of the configuration.</li></ul>
 
<font color="purple">Schwarzschild criterion:</font>&nbsp; In order to avoid constructing configurations that violate the Schwarzschild criterion, be sure that the specific entropy of the fluid is uniform (marginally stable) or increases outward throughout the equilibrium structure, where the word "outward" is only meaningful when referenced against the direction that the effective gravity points.  For example &hellip;
<ul><li>Suppose that you build a spherically symmetric polytropic configuration whose ''structural'' index is, <math>~n</math>, but then you want to test the stability of the configuration assuming that compressions/expansions of individual fluid elements occur along adiabats for which the adiabatic index is, <math>~\gamma \ne (n+1)/n</math>.  Although we generally think of polytropes as being homentropic configurations, if <math>~\gamma \ne (n+1)/n</math>, then different fluid elements will, in practice, evolve along adiabats that are characterized by ''different'' values of the specific entropy; that is, throughout the equilibrium model, <math>~\nabla s</math> will not be zero.  Whether the specific entropy increase (stable) or decreases (unstable) outward will depend on whether you select a value for the ''evolutionary'' <math>~\gamma</math> that is greater than (stable) or less than (unstable) <math>~(n+1)/n</math>.</li>
 
<li>According to [http://www.ucolick.org/~woosley/ay112-14/lectures/lecture8.4x.pdf Woosley's class lecture notes],
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{d\ln \rho}{d\ln P}\biggr|_\mathrm{structure} = \frac{n}{n+1}</math>
  </td>
  <td align="center">
<math>~></math>
  </td>
  <td align="left">
<math>~\frac{1}{\gamma_g}</math>
  </td>
  <td align="center">
<math>~\Rightarrow</math>&nbsp; &nbsp; &nbsp; stable
  </td>
</tr>


==Lord Rayleigh Instability==
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~<</math>
  </td>
  <td align="left">
<math>~\frac{1}{\gamma_g}</math>
  </td>
  <td align="center">
<math>~\Rightarrow</math>&nbsp; &nbsp; &nbsp; unstable
  </td>
</tr>
</table>
This is another way of expressing the same stability criterion for polytopes.
</li>
<li>
Examples:  &nbsp; An n = 1 polytope is ''unstable toward convection'' if expansions (or contractions) occur along an adiabat with  <math>~\gamma_g < 2</math>.  Alternatively, an n = 5 polytope is  ''unstable toward convection'' if expansions (or contractions) occur along an adiabat with <math>~\gamma_g < \tfrac{6}{5}</math>.
</li></ul>


=See Also=
=See Also=
 
* [http://adsabs.harvard.edu/abs/1917RSPSA..93..148R Lord Rayleigh (1917, Proc. Royal Society of London. Series A, 93, 148-154)] &#8212; ''On the Dynamics of Revolving Fluids''
* [http://adsabs.harvard.edu/abs/2016AJ....152...35F T. Fukushima (2016, AJ, 152, article id. 35, 31 pp.)] &#8212; ''Zonal Toroidal Harmonic Expansions of External Gravitational Fields for Ring-like Objects''
* [http://adsabs.harvard.edu/abs/2016ApJ...829...45K W.-T. Kim &amp; S. Moon (2016, ApJ, 829, article id. 45, 22 pp.)] &#8212; ''Equilibrium Sequences and Gravitational Instability of Rotating Isothermal Rings''
* [http://adsabs.harvard.edu/abs/2008MNRAS.389..156P D. Petroff &amp; S. Horatschek (2008, MNRAS, 389,156 - 172)] &#8212; ''Uniformly Rotating Homogeneous and Polytropic Rings in Newtonian Gravity''
* [http://adsabs.harvard.edu/abs/2006IJMPB..20.3113C P. H. Chavanis (2006, International Journal of Modern Physics B, 20, 3113 - 3198)] &#8212; ''Phase Transitions in Self-Gravitating Systems''
* [http://adsabs.harvard.edu/abs/2001A%26A...375.1091L M. Lombardi &amp; G. Bertin (2001, Astronomy &amp; Astrophysics, 375, 1091 - 1099)] &#8212; ''Boyle's Law and Gravitational Instability''
* [http://adsabs.harvard.edu/abs/1994ApJ...420..247W J. W. Woodward, J. E. Tohline, &amp; I. Hachisu (1994, ApJ, 420, 247 - 267)] &#8212; ''The Stability of Thick, Self-Gravitating Disks in Protostellar Systems''
* [http://adsabs.harvard.edu/abs/1991ApJ...374..610B I. Bonnell &amp; P. Bastien (1991, ApJ, 374, 610 - 622)] &#8212; ''The Collapse of Cylindrical Isothermal and Polytropic Clouds with Rotation''
* [http://adsabs.harvard.edu/abs/1990ApJ...361..394T J. E. Tohline &amp; I. Hachisu (1990, ApJ, 361, 394 - 407)] &#8212; ''The Breakup of Self-Gravitating Rings, Tori, and Thick Accretion Disks''
* [http://adsabs.harvard.edu/abs/1988A%26A...200..127S F. Schmitz (1988, Astronomy &amp; Astrophysics, 200, 127 - 134)] &#8212; ''Equilibrium Structures of Differentially Rotating Self-Gravitating Gases''
* [http://adsabs.harvard.edu/abs/1985Ap%26SS.109...45V P. Veugelen (1985, Astrophysics &amp; Space Science, 109, 45 - 55)] &#8212; ''Equilibrium Models of Differentially Rotating Polytropic Cylinders''
* [http://adsabs.harvard.edu/abs/1984MNRAS.208..279A M. A. Abramowicz, A. Curir, A. Schwarzenberg-Czerny, &amp; R. E. Wilson (1984, MNRAS, 208, 279 - 291)] &#8212; ''Self-Gravity and the Global Structure of Accretion Discs''
* [http://adsabs.harvard.edu/abs/1983A%26A...119..109B P. Bastien (1983, Astronomy &amp; Astrophysics, 119, 109 - 116)] &#8212; ''Gravitational Collapse and Fragmentation of Isothermal, Non-Rotating, Cylindrical Clouds''
* [http://adsabs.harvard.edu/abs/1981PThPh..65.1870E Y. Eriguchi &amp; D. Sugimoto (1981, Progress of Theoretical Physics, 65, 1870 - 1875)] &#8212; ''Another Equilibrium Sequence of Self-Gravitating and Rotating Incompressible Fluid''
* [http://adsabs.harvard.edu/abs/1980ApJ...236..160T J. E. Tohline (1980, ApJ, 236, 160 - 171)] &#8212; ''Ring Formation in Rotating Protostellar Clouds''
* [http://adsabs.harvard.edu/abs/1980PThPh..63.1957F T. Fukushima, Y. Eriguchi, D. Sugimoto, &amp; G. S. Bisnovatyi-Kogan (1980, Progress of Theoretical Physics, 63, 1957 - 1970)] &#8212; ''Concave Hamburger Equilibrium of Rotating Bodies''
* [http://adsabs.harvard.edu/abs/1978MNRAS.184..709K J. Katz &amp; D. Lynden-Bell (1978, MNRAS, 184, 709 - 712)] &#8212; ''The Gravothermal Instability in Two Dimensions''
* [http://adsabs.harvard.edu/abs/1977ApJ...214..584M P. S. Marcus, W. H. Press, &amp; S. A. Teukolsky (1977, ApJ, 214, 584- 597)] &#8212; ''Stablest Shapes for an Axisymmetric Body of Gravitating, Incompressible Fluid'' (includes torus with non-uniform rotation)
* [http://adsabs.harvard.edu/abs/1977ApJ...214..584M P. S. Marcus, W. H. Press, &amp; S. A. Teukolsky (1977, ApJ, 214, 584- 597)] &#8212; ''Stablest Shapes for an Axisymmetric Body of Gravitating, Incompressible Fluid'' (includes torus with non-uniform rotation)
<ol type="a"><li>Shortly after their equation (3.2), Marcus, Press &amp; Teukolsky make the following statement: "&hellip; we know that an ''equilibrium'' incompressible configuration must rotate uniformly on cylinders (the famous "Poincar&eacute;-Wavre" theorem, cf. Tassoul 1977, &Sect;4.3) &hellip;"</li></ol>
<ol type="a"><li>Shortly after their equation (3.2), Marcus, Press &amp; Teukolsky make the following statement: "&hellip; we know that an ''equilibrium'' incompressible configuration must rotate uniformly on cylinders (the famous "Poincar&eacute;-Wavre" theorem, cf. Tassoul 1977, &Sect;4.3) &hellip;"</li></ol>
* [http://adsabs.harvard.edu/abs/1976ApJ...207..736H C. J. Hansen, M. L. Aizenman, &amp; R. L. Ross (1976, ApJ, 207, 736 - 744)] &#8212; ''The Equilibrium and Stability of Uniformly Rotating, Isothermal Gas Cylinders''
* Referring to our [[User:Tohline/ThreeDimensionalConfigurations/RiemannTypeI#Riemann_Type_1_Ellipsoids|accompanying discussion of Type 1 Riemann ellipsoids]], it seems that uniform rotation on cylinders is not required. What's going on?
* [http://adsabs.harvard.edu/abs/1974ApJ...190..675W C.-Y. Wong (1974, ApJ, 190, 675 - 694)] &#8212; ''Toroidal Figures of Equilibrium''
* [http://adsabs.harvard.edu/abs/1973AnPhy..77..279W C.-Y. Wong (1973, Annals of Physics, 77, 279 - 353)] &#8212; ''Toroidal and Spherical Bubble Nuclei''
* [http://adsabs.harvard.edu/abs/1964ApJ...140.1056O J. Ostriker (1964, ApJ, 140, 1056)] &#8212; ''The Equilibrium of Polytropic and Isothermal Cylinders''
* [http://adsabs.harvard.edu/abs/1964ApJ...140.1067O J. Ostriker (1964, ApJ, 140, 1067)] &#8212; ''The Equilibrium of Self-Gravitating Rings''
* [http://adsabs.harvard.edu/abs/1964ApJ...140.1529O J. Ostriker (1964, ApJ, 140, 1529)] &#8212; ''On the Oscillations and the Stability of a Homogeneous Compressible Cylinder''
* [http://adsabs.harvard.edu/abs/1965ApJS...11..167O J. Ostriker (1965, ApJ Supplements, 11, 167)] &#8212; ''Cylindrical Emden and Associated Functions''
* [http://adsabs.harvard.edu/abs/1942ApJ....95...88R Gunnar Randers (1942, ApJ, 95, 88)] &#8212; ''The Equilibrium and Stability of Ring-Shaped 'barred SPIRALS'.''
* [http://adsabs.harvard.edu/abs/1893RSPTA.184.1041D F. W. Dyson (1893, Philosophical Transaction of the Royal Society London. A., 184, 1041 - 1106)] &#8212; ''The Potential of an Anchor Ring. Part II.'' <ol type="a"><li>In this paper, Dyson derives the gravitational potential ''inside'' the ring mass distribution</li></ol>
* [http://adsabs.harvard.edu/abs/1893RSPTA.184...43D F. W. Dyson (1893, Philosophical Transaction of the Royal Society London. A., 184, 43 - 95)] &#8212; ''The Potential of an Anchor Ring. Part I.''<ol type="a"><li>In this paper, Dyson derives the gravitational potential ''exterior to'' the ring mass distribution</li></ol>
* [http://adsabs.harvard.edu/abs/1885AN....111...37K S. Kowalewsky (1885, Astronomische Nachrichten, 111, 37)] &#8212; ''Zus&auml;tze und Bemerkungen zu Laplace's Untersuchung &uuml;ber die Gestalt der Saturnsringe''
* Poincar&eacute; (1885a, C. R. Acad. Sci., 100, 346), (1885b, Bull. Astr., 2, 109), (1885c, Bull. Astr. 2, 405). &#8212; references copied from paper by [http://adsabs.harvard.edu/abs/1974ApJ...190..675W Wong (1974)]




&nbsp;<br />
&nbsp;<br />
{{LSU_HBook_footer}}
{{LSU_HBook_footer}}

Latest revision as of 22:07, 25 June 2020

Axisymmetric Instabilities to Avoid

Whitworth's (1981) Isothermal Free-Energy Surface
|   Tiled Menu   |   Tables of Content   |  Banner Video   |  Tohline Home Page   |

When constructing rotating equilibrium configurations that obey a barotropic equation of state, keep in mind that certain physical variable profiles should be avoided because they will lead to structures that are unstable toward the dynamical development of shape-distorting or convective-type motions. Here are a few well-known examples.

Rayleigh-Taylor Instability

Referencing both p. 101 of [ Shu92 ], and volume I, p. 410 of [ P00 ], a Rayleigh-Taylor instability is a bouyancy-driven instability that arises when a heavy fluid rests on top of a light fluid in an effective gravitational field, <math>~\vec{g}</math>. In the simplest case of spherically symmetric, self-gravitating configurations, the condition for stability against a Rayleigh-Taylor instability may be written as,

<math>~(- \vec{g} ) \cdot \nabla\rho</math>

<math>~< </math>

<math>~0 </math>

     [stable] ,

that is to say, the mass density must decrease outward. In an expanded discussion, [ P00 ] — see pp. 410 - 413 — derives a dispersion relation that describes the development of the Rayleigh-Taylor instability in a plane-parallel fluid layer that initially contains a discontinuous jump/drop in the density.

In addition, [ P00 ] — see pp. 413 - 416 — and [ Shu92 ] — see pp. 101 - 105 — both discuss circumstances that can give rise to the so-called Kelvin-Helmholtz instability. It is another common two-fluid instability, but one that depends on the existence of transverse velocity flows and that is independent of the gravitational field.

Poincaré-Wavre Theorem

As [ T78 ] points out — see his pp. 78 - 81 — Poincaré and Wavre were the first to, effectively, prove the following theorem:

For rotating, self-gravitating configurations "any of the following statements implies the three others:   (i) the angular velocity is a constant over cylinders centered about the axis of rotation, (ii) the effective gravity can be derived from a potential, (iii) the effective gravity is normal to the isopycnic surfaces, (iv) the isobaric- and isopycnic-surfaces coincide."

Among other things, this implies that for rotating barotropic configurations not only is the equation of state given by a function of the form, <math>~P = P(\rho)</math>, but it must also be true that,

<math>~\frac{\partial \dot\varphi}{\partial z}</math>

<math>~=</math>

<math>~0 \, .</math>

[ T78 ], §4.3, Eq. (30)

NOTE: We should investigate how this theorem comes into play in the context of our accompanying discussion of Type 1 Riemann ellipsoids. These are equilibrium triaxial, uniform-density configurations in which the system's internal vorticity vector does not align with the tumble-axis of the ellipsoid and, therefore apparently <math>~\dot\varphi</math> is not independent of <math>~z</math>.

Høiland Criterion

"For an incompressible liquid contained between concentric cylinders and rotating with angular velocity <math>~\Omega(\varpi)</math>, Rayleigh's criterion — that <math>~\varpi^4 \Omega^2</math> increase outwards — is necessary and sufficient for stability to axisymmetric disturbances. In a star rotating with angular velocity <math>~\Omega(\varpi)</math> if we continue to restrict attention to axisymmetric disturbances, this criterion must be modified by buoyancy effects; that is, some combination of the Rayleigh and Schwarzschild criteria should obtain. Such a combination has been found, for example, by Høiland (1941) (see also Ledoux's Chapter 10, pp. 499-574 of Stellar Structure (1965), and indicates, as one would expect, that a stable stratification of angular velocity exerts a stabilizing influence on an unstable distribution of temperature, and vice versa. The combined criterion has not been placed on the solid analytical foundation of its two component criteria, however."

— Drawn from p. 475 of N. R. Lebovitz (1967), ARAA, 5, 465

As is stated on p. 166 of [ T78 ], in rotating barotropic configurations, axisymmetric stability requires the simultaneous satisfaction of the following pair of conditions:

<math>~\biggl(\frac{1}{\varpi^3} \biggr) \frac{\partial j^2}{\partial \varpi} + \frac{1}{c_P} \biggl( \frac{\gamma - 1}{\Gamma_3 - 1}\biggr) (- \vec{g} ) \cdot \nabla s</math>

<math>~></math>

<math>~0 </math>

     [stable] ;

[ T78 ], §7.3, Eq. (41)
see also
[ KW94 ], §43.2, Eq. (43.22)

<math>~-g_z \biggl[ \frac{\partial j^2}{\partial \varpi} \biggl(\frac{\partial s}{\partial z} \biggr) - \frac{\partial j^2}{\partial z} \biggl(\frac{\partial s}{\partial \varpi} \biggr)\biggr]</math>

<math>~></math>

<math>~0 </math>

     [stable] .

[ T78 ], §7.3, Eq. (42)
see also
[ KW94 ], §43.2, Eq. (43.23)

where, <math>~s</math>, is the local specific entropy, and <math>~j \equiv \dot\varphi \varpi^2</math>, is the local specific angular momentum of the fluid. According to [ T78 ] — see p. 168 — this pair of mathematically expressed conditions has the following meaning:

"A baroclinic star in permanent rotation is dynamically stable with respect to axisymmetric motions if and only if the two following conditions are satisfied:   (i) the entropy per unit mass, <math>~s</math>, never decreases outward, and (ii) on each surface <math>~s</math> = constant, the angular momentum per unit mass, <math>~j</math>, increases as we move from the poles to the equator."

Schwarzschild Criterion

In the case of nonrotating equilibrium configurations, the Høiland Criterion reduces to the Schwarzschild criterion. That is, thermal convection arises when the condition,

<math>~(- \vec{g} ) \cdot \nabla s</math>

<math>~> </math>

<math>~0 </math>

     [stable] ,

[ T78 ], §7.3, Eq. (43)
see also
[ KW94 ], §6.1, Eq. (6.13) … or … pp. 93 - 98 of [ Shu92 ]

is violated. This means that, in order for a spherical system to be stable against dynamical convective motions, the specific entropy must increase outward.

Solberg/Rayleigh Criterion

In the case of an homentropic equilibrium configuration, the Høiland Criterion reduces to the Solberg criterion. That is, an axisymmetric exchange of fluid "rings" will occur on a dynamical time scale if the condition,

<math>~\frac{dj^2}{d\varpi} </math>

<math>~> </math>

<math>~0 </math>

     [stable] ,

[ T78 ], §7.3, Eq. (44)
see also
[ KW94 ], §43.2, Eq. (43.18) … or … pp. 98 - 101 of [ Shu92 ]

is violated. This means that, for stability, the specific angular momentum must necessarily increase outward. As [ T78 ] points out, this "Solberg criterion generalizes to homentropic bodies the well-known Rayleigh (1917) criterion for an inviscid, incompressible fluid."

Modeling Implications and Advice

Here are some recommendations to keep in mind as you attempt to construct equilibrium models of self-gravitating astrophysical fluids.

Rayleigh-Taylor instability:  In order to avoid constructing configurations that are subject to the Rayleigh-Taylor instability, be sure that lower density material is never placed "beneath" higher density material. For example …

  • When building a bipolytropic configuration, a value for the (imposed) discontinuous jump in the mean-molecular weight will need to be specified at the interface between the envelope and the core. If you choose a ratio, <math>~{\bar\mu}_e/{\bar\mu}_c</math>, that is greater than unity, the resulting equilibrium model will exhibit a discontinuous density jump that makes the density higher at the base of the envelope than it is at the surface of the core. The core/envelope interface of this configuration will be unstable to the Rayleigh-Taylor instability, but you won't know that until and unless you examine the hydrodynamic stability of the configuration.

Schwarzschild criterion:  In order to avoid constructing configurations that violate the Schwarzschild criterion, be sure that the specific entropy of the fluid is uniform (marginally stable) or increases outward throughout the equilibrium structure, where the word "outward" is only meaningful when referenced against the direction that the effective gravity points. For example …

  • Suppose that you build a spherically symmetric polytropic configuration whose structural index is, <math>~n</math>, but then you want to test the stability of the configuration assuming that compressions/expansions of individual fluid elements occur along adiabats for which the adiabatic index is, <math>~\gamma \ne (n+1)/n</math>. Although we generally think of polytropes as being homentropic configurations, if <math>~\gamma \ne (n+1)/n</math>, then different fluid elements will, in practice, evolve along adiabats that are characterized by different values of the specific entropy; that is, throughout the equilibrium model, <math>~\nabla s</math> will not be zero. Whether the specific entropy increase (stable) or decreases (unstable) outward will depend on whether you select a value for the evolutionary <math>~\gamma</math> that is greater than (stable) or less than (unstable) <math>~(n+1)/n</math>.
  • According to Woosley's class lecture notes,

    <math>~\frac{d\ln \rho}{d\ln P}\biggr|_\mathrm{structure} = \frac{n}{n+1}</math>

    <math>~></math>

    <math>~\frac{1}{\gamma_g}</math>

    <math>~\Rightarrow</math>      stable

     

    <math>~<</math>

    <math>~\frac{1}{\gamma_g}</math>

    <math>~\Rightarrow</math>      unstable

    This is another way of expressing the same stability criterion for polytopes.

  • Examples:   An n = 1 polytope is unstable toward convection if expansions (or contractions) occur along an adiabat with <math>~\gamma_g < 2</math>. Alternatively, an n = 5 polytope is unstable toward convection if expansions (or contractions) occur along an adiabat with <math>~\gamma_g < \tfrac{6}{5}</math>.

See Also

  1. Shortly after their equation (3.2), Marcus, Press & Teukolsky make the following statement: "… we know that an equilibrium incompressible configuration must rotate uniformly on cylinders (the famous "Poincaré-Wavre" theorem, cf. Tassoul 1977, &Sect;4.3) …"


 

Whitworth's (1981) Isothermal Free-Energy Surface

© 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