Difference between revisions of "User:Tohline/Apps/RotatingPolytropes/BarmodeLinearTimeDependent"

From VistrailsWiki
Jump to navigation Jump to search
 
(49 intermediate revisions by the same user not shown)
Line 1: Line 1:
__FORCETOC__  
__FORCETOC__  
<!-- __NOTOC__ will force TOC off -->
<!-- __NOTOC__ will force TOC off -->
=Simulating the Onset of a Barmode Instability=
=Onset of Bar-mode Instability in Initially Axisymmetric &amp; Differentially Rotating Polytropes=
<table border="0" cellpadding="3" align="center" width="60%">
<tr><td align="left">
<font color="darkgreen">
"&hellip; the onset of instability is not very sensitive to the compressibility or angular momentum distribution of the polytrope when the models are parameterized by T/|W|</font> &#8212; [in particular, the m = 2 barmode becomes unstable at T/|W| &sim; 0.26 - 0.28. ] <font color="darkgreen">The polytrope eigenfunctions are &hellip; qualitatively different from the Maclaurin eigenfunctions in one respect: they develop strong spiral arms. The spiral arms are stronger for more compressible polytropes and for polytropes whose angular momentum distributions deviate significantly from those of the Maclaurin spheroids."
</font>
</td></tr>
<tr><td align="right">
&#8212; Drawn from [https://ui.adsabs.harvard.edu/abs/1998ApJ...497..370T/abstract Toman, Imamura, Pickett &amp; Durisen (1998)], ApJ, 497, 370
</td></tr></table>


{{LSU_HBook_header}}
{{LSU_HBook_header}}
==Summary ''circa'' 2000==
<table border="0" align="right" cellpadding="10"><tr><td align="center">
<table border="1" align="right" cellpadding="5">
<tr><td align="center" bgcolor="black">[[File:Dissertation.fig3cropped.png|250px|Model A Composite Images]]</td></tr>
<tr><td align="center">Cazes' Model A</td></tr>
</table>
</td></tr></table>
This summary is drawn &#8212; largely verbatim &#8212; from a [http://www.phys.lsu.edu/astro/movie_captions/two_armed_spiral.html separate web page] that was originally intended to highlight animation sequences produced by LSU's astrophysicsl theory group.
In the mid-1980s, using three different 3D hydrodynamic simulation tools, [https://ui.adsabs.harvard.edu/abs/1986ApJ...305..281D/abstract R. H. Durisen, R. A. Gingold, J. E. Tohline &amp; A. P. Boss (1986)] &#8212; hereafter and in the references, below, DGTB86 &#8212; followed the development to nonlinear amplitudes of a two-armed, spiral-mode instability that arises naturally in rapidly rotating, self-gravitating fluids.  The particular initial model that was examined had an <math>~n = \tfrac{3}{2}</math> polytropic equation of state, an n' = 0 angular momentum distribution, and an initial ratio of rotational to gravitational potential energy T/|W| = 0.33.  The simulation tools that were used at the time were relatively crude, compared to tools that are available today.
Employing a significantly improved finite-volume simulation code and improved spatial resolution (128<sup>3</sup> grid zones), a very similar simulation has been carried out and reported by [https://ui.adsabs.harvard.edu/abs/2000ApJ...532.1051C/abstract J. E. Cazes &amp; J. E. Tohline (2000)]; see also the [https://digitalcommons.lsu.edu/gradschool_disstheses/6982/ Cazes (1999)] doctoral dissertation.  This was done, in part, to check the validity of the earlier work and, in part, to permit us to conduct a much more thorough analysis of the end-state configuration.  We will refer to this as the ''Model A'' simulation.  The initial model for this simulation was constructed using the [[User:Tohline/AxisymmetricConfigurations/HSCF#Hachisu_Self-Consistent-Field_Technique|Hachisu Self-Consistent-Field]] (HSCF) technique; it had an <math>~n = \tfrac{3}{2}</math> polytropic equation of state, an n' = 0 angular momentum distribution, and an initial ratio of rotational to gravitational potential energy T/|W| = 0.30.  As has been cataloged in [[#Table1|Table 1, below]], a YouTube animation carrying ID = [https://youtu.be/BhRUqZe0Ly4 BhRUqZe0Ly4] shows the nonlinear development of the bar-mode &#8212; actually, two-armed, spiral-mode &#8212; instability in Model A.  The evolution is shown in the inertial reference frame and covers 46 dynamical times (20 central initial rotation periods) as defined by the properties of the initial, equilibrium axisymmetric model.  Each frame of the ''Model A'' movie displays four nested isodensity contours at 80%, 40%, 4% and 0.4% of the maximum density.  Via the trailing spiral structure, gravitational torques are able to effectively redistribute angular momentum on a dynamical time scale; a relatively small amount of material is shed into an equatorial disk.  (This disk material is not visible in the ''Model A'' movie because the disk material all has a mass density less than 0.4% of the maximum density.)  Over time, the central object (containing most of the initial object's mass) settles down into a new quasi-equilibrium configuration.  The six images displayed in the upper portion of color figure shown here, on the right, have been extracted from this animation at six different points in time (measured in terms of the central initial rotation period), as labeled.
For comparison, we also have studied the development of the same type of bar-mode (two-armed, spiral) instability in a model which will be referred to here as ''Model B''.  It had an <math>~n = \tfrac{3}{2}</math> polytropic equation of state and an initial ratio of rotational to gravitational potential energy T/|W| = 0.28, but with an angular momentum distribution specified to produce a uniform ''vortensity'' profile in the initial model.  (Vortensity is defined as the ratio of vorticity to mass density.)  A YouTube animation carrying ID = [https://youtu.be/Qh-bMto2e3k Qh-bMto2e3k] shows this evolution from an inertial frame and covers approximately 32 dynamical times as defined by the mean density of the initial model.  The instability in this model has less of a pronounced spiral character, but ultimately results in the formation of a new triaxial, quasi-equilibrium configuration with properties that are similar to the end state of the ''Model A'' evolution.
As illustrated and discussed elsewhere, we are convinced that the ellipsoidal-like configurations that have been formed through both of these model simulations are compressible analogs of Riemann S-type ellipsoids.
<div id="Table1">
<table border="1" width="85%" align="center" cellpadding="8">
<tr>
<td align="center" colspan="4"><b>Table 1: &nbsp; Accompanying YouTube Animations</b> (click on ID)</td>
</tr>
<tr>
<td align="right" colspan="2" width="50%">Model:</td>
<td align="center" colspan="1">A</td>
<td align="center" colspan="1">B</td>
</tr>
<tr>
<td align="right" colspan="2">Nonlinear Development of Bar-mode Instability (isodensity surfaces):</td>
<td align="center" colspan="1"><b><font color="darkgreen">[https://youtu.be/BhRUqZe0Ly4 BhRUqZe0Ly4]</font></b></td>
<td align="center" colspan="1"><b><font color="darkgreen">[https://youtu.be/Qh-bMto2e3k Qh-bMto2e3k]</font></b></td>
</tr>
<tr>
<td align="right" colspan="1" rowspan="4">Quasi-Steady-State "CARE"</td>
<td align="right" colspan="1">Isodensity Surfaces:</td>
<td align="center" colspan="1"><b><font color="darkgreen">[https://youtu.be/z3RDthHnY3A z3RDthHnY3A]</font></b></td>
<td align="center" colspan="1"><b><font color="darkgreen">[https://youtu.be/Eyl81O9s62s Eyl81O9s62s]</font></b></td>
</tr>
<tr>
<td align="right" colspan="1">Equatorial-Plane Fluid Flow:</td>
<td align="center" colspan="1"><b><font color="darkgreen">[https://youtu.be/fgjjzu3Nmo8 fgjjzu3Nmo8]</font></b></td>
<td align="center" colspan="1"><b><font color="darkgreen">[https://youtu.be/yggbz9z8JQ4 yggbz9z8JQ4]</font></b></td>
</tr>
<tr>
<td align="right" colspan="1">Andalib's Analogous ''Prograde'' Flow-Field:</td>
<td align="center" colspan="2"><b><font color="darkgreen">[https://youtu.be/p550IDP8h-s p550IDP8h-s]</font></b></td>
</tr>
<tr>
<td align="right" colspan="1">Meridional Fluid Flow:</td>
<td align="center" colspan="1"><b><font color="darkgreen">[https://youtu.be/_nXKv0k0h-Q _nXKv0k0h-Q]</font></b></td>
<td align="center" colspan="1">n/a</td>
</tr>
<tr>
<td align="right" colspan="1" rowspan="4">''Cooling'' Evolution</td>
<td align="right" colspan="1">Isodensity Surfaces:</td>
<td align="center" colspan="1"><b><font color="darkgreen">[https://youtu.be/erYWHgaKYyc erYWHgaKYyc]</font></b></td>
<td align="center" colspan="1">n/a</td>
</tr>
<tr>
<td align="right" colspan="1">Equatorial-Plane Fluid Flow:</td>
<td align="center" colspan="1"><b><font color="darkgreen">[https://youtu.be/nC4adQzlXrI nC4adQzlXrI]</font></b></td>
<td align="center" colspan="1">n/a</td>
</tr>
<tr>
<td align="right" colspan="1">Andalib's Analogous ''Binary'' Flow-Field:</td>
<td align="center" colspan="1"><b><font color="darkgreen">[https://youtu.be/nkwRJh2tgos nkwRJh2tgos]</font></b></td>
<td align="center" colspan="1">n/a</td>
</tr>
<tr>
<td align="right" colspan="1">Meridional Fluid Flow:</td>
<td align="center" colspan="1"><b><font color="darkgreen">[https://youtu.be/S0tg2jURzGE S0tg2jURzGE]</font></b></td>
<td align="center" colspan="1">n/a</td>
</tr>
</table>
</div>


==Index of Relevant Publications==
==Index of Relevant Publications==
Here is a list of relevant research papers as enumerated by [https://ui.adsabs.harvard.edu/abs/2008PhRvD..78l4001K/abstract Y. Kojima &amp; M. Saijo (2008)], Phys. Rev. D, vol. 78, Issue 12, id. 124001:  ''Amplification of azimuthal modes with odd wave numbers during dynamical bar-mode growth in rotating stars''  
Here is a list of relevant research papers as largely enumerated by [https://ui.adsabs.harvard.edu/abs/2008PhRvD..78l4001K/abstract Y. Kojima &amp; M. Saijo (2008)] &hellip;
 
* [TDM85] [https://ui.adsabs.harvard.edu/abs/1985ApJ...298..220T/abstract J. E. Tohline, R. H. Durisen &amp; M. McCollough (1985)], ApJ, 298, 220:  ''The linear and nonlinear dynamic stability of rotating N = 3/2 polytropes''
<table border="0" align="center" width="100%" cellpadding="1"><tr>
<td align="center" width="5%">&nbsp;</td><td align="left">
HYDROCODE:  Newtonian, 3D Eulerian, 1<sup>st</sup>-order donor-cell on a cylindrical grid; &pi;-symmetry plus reflection symmetry through equatorial plane; <math>~\Gamma=5/3</math>; Poisson solved with FFT + Buneman cyclic reduction<br />
 
MODEL(s): Constructed using Ostriker-Mark SCF method; axisymmetric, n = 3/2 polytrope; n' = 0 rotation law; four different equilibrium configurations having T/|W| = 0.28, 0.30, 0.33, 0.35.
</td></tr></table>
 
* [DGTB86] [https://ui.adsabs.harvard.edu/abs/1986ApJ...305..281D/abstract R. H. Durisen, R. A. Gingold, J. E. Tohline &amp; A. P. Boss (1986)], ApJ, 305, 281: ''Dynamic Fission Instabilities in Rapidly Rotating N = 3/2 Polytropes: A Comparison of Results from Finite-Difference and Smoothed Particle Hydrodynamics Codes''
* [WT87] [https://ui.adsabs.harvard.edu/abs/1987ApJ...315..594W/abstract H. A. Williams &amp; J. E. Tohline (1987)], ApJ, 315, 594:  ''Linear and Nonlinear Dynamic Instability of Rotating Polytropes''
<table border="0" align="center" width="100%" cellpadding="1"><tr>
<td align="center" width="5%">&nbsp;</td><td align="left">
HYDROCODE:  Same as in [https://ui.adsabs.harvard.edu/abs/1985ApJ...298..220T/abstract Tohline, Durisen &amp; McCollough (1985)], above.<br />
 
MODEL(s): Constructed using Ostriker-Mark SCF method; axisymmetric, n' = 0 rotation law; five different equilibrium configurations having (see column 1 of their Table 1) n = 0.8, 1.0, 1.3, 1.5, 1.8, all having T/|W| = 0.310.
</td></tr></table>
 
* [WT88] [https://ui.adsabs.harvard.edu/abs/1988ApJ...334..449W/abstract H. A. Williams &amp; J. E. Tohline (1988)], ApJ, 334, 449:  ''Circumstellar Ring Formation in Rapidly Rotating Protostars''
<table border="0" align="center" width="100%" cellpadding="1"><tr>
<td align="center" width="5%">&nbsp;</td><td align="left">
Two of the models that were studied in Williams &amp; Tohline (1987) &#8212; specifically, the models having n = 0.8 and 1.8 &#8212; <font color="green">
&hellip; are shown as they evolve to extremely nonlinear amplitudes: the end result in both cases is &hellip; The models shed a fraction of their mass and angular momentum, producing a ring which surrounds a more centrally condensed object &hellip; The central object is a triaxial figure that is rotating about its shortest axis.
</font><br />
</td></tr></table>
* [HCS94] [https://ui.adsabs.harvard.edu/abs/1994PhRvL..72.1314H/abstract J. L. Houser, J. M. Centrella &amp; S. C. Smith (1994)], PRL, 72, 1314: ''Gravitational radiation from nonaxisymmetric instability in a rotating star''
<table border="0" align="center" width="100%" cellpadding="1"><tr>
<td align="center" width="5%">&nbsp;</td><td align="left">
Using Newtonian dynamics and Newtonian gravity <font color="green">&hellip; we have carried out computer simulations of a differentially rotating compact star with a polytropic equation of state undergoing the dynamical bar instability.  This instability has previously been modeled numerically by Tohline and collaborators in the context of star formation &hellip; Our work is the first to calculate</font> [using a post-Newtonian approximation] <font color="green"> the gravitational radiation produced by this instability, including wave forms and luminosities.  It is also a significant advance over the earlier studies because, in addition to using better numerical techniques, we model the fluid correctly using an energy equation.  This is essential due to the generation of entropy by shocks during the later stages of the evolution.
</font><br />
 
HYDROCODE:  Newtonian, 3D Lagrangian-Cartesian (SPH as implemented by [https://ui.adsabs.harvard.edu/abs/1989ApJS...70..419H/abstract Hernquist &amp; Katz 1989]), Includes energy equation and does ''not'' impose &pi;-symmetry; <math>~\Gamma=5/3</math>; Poisson solved with TREESPH<br />
 
MODEL(s): Constructed using the SCF method of Smith &amp; Centrella (1992; see article #15 in [https://ui.adsabs.harvard.edu/abs/2005anr..book.....D/abstract R. d'Inverno]), which is based on earlier work of [https://ui.adsabs.harvard.edu/abs/1968ApJ...151.1075O/abstract Ostriker-Mark 1968] and [https://ui.adsabs.harvard.edu/abs/1986ApJS...61..479H/abstract Hachisu 1986]; axisymmetric, n = 3/2 polytrope; n' = 0 rotation law; only one equilibrium configuration, with T/|W| = 0.30.
</td></tr></table>
 
* [SHC95] [https://ui.adsabs.harvard.edu/abs/1996ApJ...458..236S/abstract S. C. Smith, J. L. Houser &amp; J. M. Centrella (1995)], ApJ, 458, 236: ''Simulations of Nonaxisymmetric Instability in a Rotating Star: A Comparison between Eulerian and Smooth Particle Hydrodynamics''
* [HC96] [https://ui.adsabs.harvard.edu/abs/1996PhRvD..54.7278H/abstract J. L. Houser &amp; J. M. Centrella (1996)], Phys. Rev. D, 54, 7278:  ''Gravitational radiation from rotational instabilites in compact stellar cores with stiff equations of state''
<table border="0" align="center" width="100%" cellpadding="1"><tr>
<td align="center" width="5%">&nbsp;</td><td align="left">
Extending the work of Houser, Centrella &amp; Smith (1994) just mentioned &#8212; and closely paralleling the work of Williams &amp; Tohline (1987, 1988) &#8212; they examine the development of the dynamical bar instability in configurations that have a range of different equations of state.
 
HYDROCODE:  Same as Houser, Centrella &amp; Smith (1994) <br />
 
MODEL(s): Axisymmetric, n' = 0 rotation law; a total of 3 models examined with, respectively,  (n, T/|W|) = (0.5, 0.31), (1.0, 0.32), (1.5, 0.32) &#8212; see their Table II.  The SCF-based method used to construct these models and the (rather contrived) method used to map the continuum distribution of mass into the Lagrangian-based SPH code is discussed in detail in &sect;III of the paper.
</td></tr></table>
 
====Break 1====
* [TIPD98] [https://ui.adsabs.harvard.edu/abs/1998ApJ...497..370T/abstract J. Toman, J. N. Imamura, B. J. Pickett &amp; R. H. Durisen (1998)], ApJ, 497, 370:  ''Nonaxisymmetric Dynamic Instabilities of Rotating Polytropes. I. The Kelvin Modes''
<table border="1" align="center" cellpadding="5">
<tr>
  <th align="center" colspan="5">Equilibrium Models constucted and examined by<br />
[https://ui.adsabs.harvard.edu/abs/1998ApJ...497..370T/abstract TIPD98]
</th>
</tr>
<tr>
  <td align="center"><math>~n</math></td>
  <td align="center"><math>~n'</math></td>
  <td align="center"><math>~T/|W|</math></td>
  <td align="center"><math>~\Omega_p/\Omega_\mathrm{eq}</math></td>
  <td align="center"><math>~R_p/R_\mathrm{eq}</math></td>
</tr>
<tr>
  <td align="center" rowspan=6><math>~\tfrac{5}{4}</math></td>
  <td align="center" rowspan="15">0</td>
  <td align="right">0.274</td>
  <td align="right">2.44</td>
  <td align="left">0.258</td>
</tr>
<tr>
  <td align="right">0.279</td>
  <td align="right">2.45</td>
  <td align="left">0.250</td>
</tr>
<tr>
  <td align="right">0.288</td>
  <td align="right">2.49</td>
  <td align="left">0.234</td>
</tr>
<tr>
  <td align="right">0.319</td>
  <td align="right">2.65</td>
  <td align="left">0.187</td>
</tr>
<tr>
  <td align="right">0.346</td>
  <td align="right">2.88</td>
  <td align="left">0.147</td>
</tr>
<tr>
  <td align="right">0.363</td>
  <td align="right">3.11</td>
  <td align="left">0.123</td>
</tr>
<tr>
  <td align="center" rowspan="16"><math>~\tfrac{3}{2}</math></td>
  <td align="right">0.266</td>
  <td align="right">2.97</td>
  <td align="left">0.258</td>
</tr>
<tr>
  <td align="right">0.280</td>
  <td align="right">3.03</td>
  <td align="left">0.234</td>
</tr>
<tr>
  <td align="right">0.290</td>
  <td align="right">3.09</td>
  <td align="left">0.218</td>
</tr>
<tr>
  <td align="right">0.300</td>
  <td align="right">3.16</td>
  <td align="left">0.202</td>
</tr>
<tr>
  <td align="right">0.316</td>
  <td align="right">3.28</td>
  <td align="left">0.179</td>
</tr>
<tr>
  <td align="right">0.327</td>
  <td align="right">3.35</td>
  <td align="left">0.163</td>
</tr>
<tr>
  <td align="right">0.344</td>
  <td align="right">3.61</td>
  <td align="left">0.139</td>
</tr>
<tr>
  <td align="right">0.368</td>
  <td align="right">4.04</td>
  <td align="left">0.107</td>
</tr>
<tr>
  <td align="right">0.396</td>
  <td align="right">4.79</td>
  <td align="left">0.0754</td>
</tr>
<tr>
  <td align="center" rowspan="4">1</td>
  <td align="right">0.262</td>
  <td align="right">6.85</td>
  <td align="left">0.139</td>
</tr>
<tr>
  <td align="right">0.283</td>
  <td align="right">8.90</td>
  <td align="left">0.107</td>
</tr>
<tr>
  <td align="right">0.292</td>
  <td align="right">9.81</td>
  <td align="left">0.0958</td>
</tr>
<tr>
  <td align="right">0.311</td>
  <td align="right">12.4 &nbsp;</td>
  <td align="left">0.0754</td>
</tr>
<tr>
  <td align="center" rowspan="3"><math>~\tfrac{3}{2}</math></td>
  <td align="right">0.272</td>
  <td align="right">15.1 &nbsp;</td>
  <td align="left">0.0754</td>
</tr>
<tr>
  <td align="right">0.280</td>
  <td align="right">16.9 &nbsp;</td>
  <td align="left">0.0675</td>
</tr>
<tr>
  <td align="right">0.290</td>
  <td align="right">19.0 &nbsp;</td>
  <td align="left">0.0595</td>
</tr>
<tr>
  <td align="center" rowspan="10"><math>~\tfrac{5}{2}</math></td>
  <td align="center" rowspan="10">0</td>
  <td align="right">0.262</td>
  <td align="right">7.19</td>
  <td align="left">0.202</td>
</tr>
<tr>
  <td align="right">0.268</td>
  <td align="right">7.26</td>
  <td align="left">0.194</td>
</tr>
<tr>
  <td align="right">0.273</td>
  <td align="right">7.34</td>
  <td align="left">0.187</td>
</tr>
<tr>
  <td align="right">0.285</td>
  <td align="right">7.52</td>
  <td align="left">0.171</td>
</tr>
<tr>
  <td align="right">0.297</td>
  <td align="right">7.73</td>
  <td align="left">0.155</td>
</tr>
<tr>
  <td align="right">0.304</td>
  <td align="right">7.85</td>
  <td align="left">0.147</td>
</tr>
<tr>
  <td align="right">0.324</td>
  <td align="right">8.29</td>
  <td align="left">0.123</td>
</tr>
<tr>
  <td align="right">0.338</td>
  <td align="right">8.68</td>
  <td align="left">0.107</td>
</tr>
<tr>
  <td align="right">0.371</td>
  <td align="right">9.85</td>
  <td align="left">0.0754</td>
</tr>
<tr>
  <td align="right">0.389</td>
  <td align="right">10.8 &nbsp;</td>
  <td align="left">0.0595</td>
</tr>
</table>
 
====Break 2====
* [https://ui.adsabs.harvard.edu/abs/2000PhRvD..62f4019N/abstract K. C. B. New, J. M. Centrella &amp; J. E. Tohline (2000)], Phys. Rev. D, 62, 064019:  ''Gravitational waves from long-duration simulations of the dynamical bar instability''
* [https://ui.adsabs.harvard.edu/abs/2000ApJ...542..453S/abstract M. Shibata, T. W. Baumgarte &amp; S. L. Shapiro (2000)], ApJ, 542, 453:  ''The Bar-Mode Instability in Differentially Rotating Neutron Stars: Simulations in Full General Relativity''
* [https://ui.adsabs.harvard.edu/abs/2001MNRAS.324.1063L/abstract Y.-T. Liu &amp; L. Lindblom (2001)], MNRAS, 324, 1063:  ''Models of rapidly rotating neutron stars:  remnants of accretion-induced collapse''
* [https://ui.adsabs.harvard.edu/abs/2001ApJ...548..919S/abstract M. Saijo, M. Shibata, T. W. Baumgarte &amp; S. L. Shapiro (2001)], ApJ, 548, 919: ''Dynamical Bar Instability in Rotating Stars: Effect of General Relativity''
* [https://ui.adsabs.harvard.edu/abs/2001ApJ...550L.193C/abstract J. M. Centrella, K. C. B. New, L. L. Lowe &amp; J. D. Brown (2001)], ApJ, 550, L193: ''Dynamical Rotational Instability at Low T/W''
* [https://ui.adsabs.harvard.edu/abs/2002PhRvD..65l4003L/abstract Y.-T. Liu (2002)], Phys. Rev. D, 65, 124003: ''Dynamical instability of new-born neutron stars as sources of gravitational radiation''
* [https://ui.adsabs.harvard.edu/abs/2003ApJ...595..352S/abstract M. Saijo, T. W. Baumgarte &amp; S. L. Shapiro (2003)], ApJ, 595, 352:  ''One-armed Spiral Instability in Differentially Rotating Stars''
* [https://ui.adsabs.harvard.edu/abs/2006ApJ...651.1068O/abstract S. Ou &amp; J. E. Tohline (2006)], ApJ, 651, 1068:  ''Unexpected Dynamical Instabilities in Differentially Rotating Neutron Stars''
* [https://ui.adsabs.harvard.edu/abs/2007PhRvD..75d4023B/abstract L. Baiotti, R. De Pietri, G. M. Manco &amp; L. Rezzolla (2007)], Phys. Rev. D, 75, 044023: ''Accurate simulations of the dynamical bar-mode instability in full general relativity''
* [https://ui.adsabs.harvard.edu/abs/2007CoPhC.177..288C/abstract P. Cerda-Duran, V. Quilos &amp; J. A. Font (2007)], Comp. Phys. Comm., 177, 288: ''AMR simulations of the low T/|W| bar-mode instability of neutron stars''
* [https://ui.adsabs.harvard.edu/abs/2007ApJ...665.1074O/abstract S. Ou, J. E. Tohline &amp; P. M. Motl (2007)], ApJ, 665, 1074:  ''Further Evidence for an Elliptical Instability in Rotating Fluid Bars and Ellipsoidal Stars''
* [https://ui.adsabs.harvard.edu/abs/2008PhRvD..77f3002S/abstract M. Saijo &amp; Y. Kojima (2008)], Phys. Rev. D, 77, 063002:  ''Faraday resonance in dynamical bar instability of differentially rotating stars''
* [https://ui.adsabs.harvard.edu/abs/2008PhRvD..78l4001K/abstract Y. Kojima &amp; M. Saijo (2008)], Phys. Rev. D, vol. 78, Issue 12, id. 124001:  ''Amplification of azimuthal modes with odd wave numbers during dynamical bar-mode growth in rotating stars''  
<table border="0" align="center" width="100%" cellpadding="1"><tr>
<td align="center" width="5%">&nbsp;</td><td align="left">
<font color="green">
Nonlinear growth of the bar-mode deformation is studied for a differentially rotating star with supercritical rotational energy. In particular, the growth mechanism of some azimuthal modes with odd wave numbers is examined &hellip; Mode coupling to even modes, i.e., the bar mode and higher harmonics, significantly enhances the amplitudes of odd modes &hellip;
</font><br />
 
HYDROCODE:  Newtonian, 3D Eulerian, Cartesian, with entropy tracer; reflection symmetry through equatorial plane; <math>~\Gamma=2</math>; Poisson solved with preconditioned conjugate gradient (PCG) method<br />
 
MODEL(s):  axisymmetric, n = 1 polytrope; [[User:Tohline/AxisymmetricConfigurations/SolutionStrategies#SRPtable|j-constant rotation law]] with A = 1; their Table I lists four different equilibrium configurations having T/|W| = 0.256, 0.268, 0.277, 0.281.
</td></tr></table>


* J. E. Tohline, R. H. Durisen &amp; M. McCollough (1985)], ApJ, 298, 220: 
Additional references identified through the above set of references:
* R. H. Durisen, R. A. Gingold, J. E. Tohline &amp; A. P. Boss (1986)], ApJ, 305, 281:
* [https://ui.adsabs.harvard.edu/abs/2018PhRvD..98b4003S/abstract M. Saijo (2018)], Phys. Rev. D, 98, 024003: ''Determining the stiffness of the equation of state using low T/W dynamical instabilities in differentially rotating stars''
* H. A. Williams &amp; J. E. Tohline (1988)], ApJ, 334, 449:
<table border="0" align="center" width="100%" cellpadding="1"><tr>
* J. L. Houser, J. M. Centrella &amp; S. C. Smith (1994)], PRL, 72, 1314:
<td align="center" width="5%">&nbsp;</td><td align="left">
* S. C. Smith, J. L. Houser &amp; J. M. Centrella (1995)], ApJ, 458, 236:
<font color="green">
* J. L. Houser &amp; J. M. Centrella (1996)], Phys. Rev. D, 54, 7278:
We investigate the nature of low T/W dynamical instabilities in various ranges of the stiffness of the equation of state in differentially rotating stars &hellip; We analyze these instabilities in both a linear perturbation analysis and a three-dimensional hydrodynamical simulation &hellip; the nature of the eigenfunction that oscillates between corotation and the surface for an unstable star requires reinterpretation of pulsation modes in differentially rotating stars.
* J. Toman, J. N. Imamura, B. J. Pickett &amp; R. H. Durisen (1998)], ApJ, 497, 370:
</font>
* K. C. B. New, J. M. Centrella &amp; J. E. Tohline (2000)], Phys. Rev. D, 62, 064019:
</td></tr></table>
* Y.-T. Liu &amp; L. Lindblom (2001)], MNRAS, 324, 1063:
* Y.-T. Liu (2002)], Phys. Rev. D, 65, 124003:
* M. Saijo &amp; Y. Kojima (2008)], Phys. Rev. D, 77, 063002:


=See Also=
=See Also=

Latest revision as of 23:06, 20 February 2020

Onset of Bar-mode Instability in Initially Axisymmetric & Differentially Rotating Polytropes

"… the onset of instability is not very sensitive to the compressibility or angular momentum distribution of the polytrope when the models are parameterized by T/|W| — [in particular, the m = 2 barmode becomes unstable at T/|W| ∼ 0.26 - 0.28. ] The polytrope eigenfunctions are … qualitatively different from the Maclaurin eigenfunctions in one respect: they develop strong spiral arms. The spiral arms are stronger for more compressible polytropes and for polytropes whose angular momentum distributions deviate significantly from those of the Maclaurin spheroids."

— Drawn from Toman, Imamura, Pickett & Durisen (1998), ApJ, 497, 370

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

Summary circa 2000

Model A Composite Images
Cazes' Model A

This summary is drawn — largely verbatim — from a separate web page that was originally intended to highlight animation sequences produced by LSU's astrophysicsl theory group.

In the mid-1980s, using three different 3D hydrodynamic simulation tools, R. H. Durisen, R. A. Gingold, J. E. Tohline & A. P. Boss (1986) — hereafter and in the references, below, DGTB86 — followed the development to nonlinear amplitudes of a two-armed, spiral-mode instability that arises naturally in rapidly rotating, self-gravitating fluids. The particular initial model that was examined had an <math>~n = \tfrac{3}{2}</math> polytropic equation of state, an n' = 0 angular momentum distribution, and an initial ratio of rotational to gravitational potential energy T/|W| = 0.33. The simulation tools that were used at the time were relatively crude, compared to tools that are available today.

Employing a significantly improved finite-volume simulation code and improved spatial resolution (1283 grid zones), a very similar simulation has been carried out and reported by J. E. Cazes & J. E. Tohline (2000); see also the Cazes (1999) doctoral dissertation. This was done, in part, to check the validity of the earlier work and, in part, to permit us to conduct a much more thorough analysis of the end-state configuration. We will refer to this as the Model A simulation. The initial model for this simulation was constructed using the Hachisu Self-Consistent-Field (HSCF) technique; it had an <math>~n = \tfrac{3}{2}</math> polytropic equation of state, an n' = 0 angular momentum distribution, and an initial ratio of rotational to gravitational potential energy T/|W| = 0.30. As has been cataloged in Table 1, below, a YouTube animation carrying ID = BhRUqZe0Ly4 shows the nonlinear development of the bar-mode — actually, two-armed, spiral-mode — instability in Model A. The evolution is shown in the inertial reference frame and covers 46 dynamical times (20 central initial rotation periods) as defined by the properties of the initial, equilibrium axisymmetric model. Each frame of the Model A movie displays four nested isodensity contours at 80%, 40%, 4% and 0.4% of the maximum density. Via the trailing spiral structure, gravitational torques are able to effectively redistribute angular momentum on a dynamical time scale; a relatively small amount of material is shed into an equatorial disk. (This disk material is not visible in the Model A movie because the disk material all has a mass density less than 0.4% of the maximum density.) Over time, the central object (containing most of the initial object's mass) settles down into a new quasi-equilibrium configuration. The six images displayed in the upper portion of color figure shown here, on the right, have been extracted from this animation at six different points in time (measured in terms of the central initial rotation period), as labeled.

For comparison, we also have studied the development of the same type of bar-mode (two-armed, spiral) instability in a model which will be referred to here as Model B. It had an <math>~n = \tfrac{3}{2}</math> polytropic equation of state and an initial ratio of rotational to gravitational potential energy T/|W| = 0.28, but with an angular momentum distribution specified to produce a uniform vortensity profile in the initial model. (Vortensity is defined as the ratio of vorticity to mass density.) A YouTube animation carrying ID = Qh-bMto2e3k shows this evolution from an inertial frame and covers approximately 32 dynamical times as defined by the mean density of the initial model. The instability in this model has less of a pronounced spiral character, but ultimately results in the formation of a new triaxial, quasi-equilibrium configuration with properties that are similar to the end state of the Model A evolution.

As illustrated and discussed elsewhere, we are convinced that the ellipsoidal-like configurations that have been formed through both of these model simulations are compressible analogs of Riemann S-type ellipsoids.


Table 1:   Accompanying YouTube Animations (click on ID)
Model: A B
Nonlinear Development of Bar-mode Instability (isodensity surfaces): BhRUqZe0Ly4 Qh-bMto2e3k
Quasi-Steady-State "CARE" Isodensity Surfaces: z3RDthHnY3A Eyl81O9s62s
Equatorial-Plane Fluid Flow: fgjjzu3Nmo8 yggbz9z8JQ4
Andalib's Analogous Prograde Flow-Field: p550IDP8h-s
Meridional Fluid Flow: _nXKv0k0h-Q n/a
Cooling Evolution Isodensity Surfaces: erYWHgaKYyc n/a
Equatorial-Plane Fluid Flow: nC4adQzlXrI n/a
Andalib's Analogous Binary Flow-Field: nkwRJh2tgos n/a
Meridional Fluid Flow: S0tg2jURzGE n/a

Index of Relevant Publications

Here is a list of relevant research papers as largely enumerated by Y. Kojima & M. Saijo (2008)

 

HYDROCODE: Newtonian, 3D Eulerian, 1st-order donor-cell on a cylindrical grid; π-symmetry plus reflection symmetry through equatorial plane; <math>~\Gamma=5/3</math>; Poisson solved with FFT + Buneman cyclic reduction

MODEL(s): Constructed using Ostriker-Mark SCF method; axisymmetric, n = 3/2 polytrope; n' = 0 rotation law; four different equilibrium configurations having T/|W| = 0.28, 0.30, 0.33, 0.35.

 

HYDROCODE: Same as in Tohline, Durisen & McCollough (1985), above.

MODEL(s): Constructed using Ostriker-Mark SCF method; axisymmetric, n' = 0 rotation law; five different equilibrium configurations having (see column 1 of their Table 1) n = 0.8, 1.0, 1.3, 1.5, 1.8, all having T/|W| = 0.310.

 

Two of the models that were studied in Williams & Tohline (1987) — specifically, the models having n = 0.8 and 1.8 — … are shown as they evolve to extremely nonlinear amplitudes: the end result in both cases is … The models shed a fraction of their mass and angular momentum, producing a ring which surrounds a more centrally condensed object … The central object is a triaxial figure that is rotating about its shortest axis.

 

Using Newtonian dynamics and Newtonian gravity … we have carried out computer simulations of a differentially rotating compact star with a polytropic equation of state undergoing the dynamical bar instability. This instability has previously been modeled numerically by Tohline and collaborators in the context of star formation … Our work is the first to calculate [using a post-Newtonian approximation] the gravitational radiation produced by this instability, including wave forms and luminosities. It is also a significant advance over the earlier studies because, in addition to using better numerical techniques, we model the fluid correctly using an energy equation. This is essential due to the generation of entropy by shocks during the later stages of the evolution.

HYDROCODE: Newtonian, 3D Lagrangian-Cartesian (SPH as implemented by Hernquist & Katz 1989), Includes energy equation and does not impose π-symmetry; <math>~\Gamma=5/3</math>; Poisson solved with TREESPH

MODEL(s): Constructed using the SCF method of Smith & Centrella (1992; see article #15 in R. d'Inverno), which is based on earlier work of Ostriker-Mark 1968 and Hachisu 1986; axisymmetric, n = 3/2 polytrope; n' = 0 rotation law; only one equilibrium configuration, with T/|W| = 0.30.

 

Extending the work of Houser, Centrella & Smith (1994) just mentioned — and closely paralleling the work of Williams & Tohline (1987, 1988) — they examine the development of the dynamical bar instability in configurations that have a range of different equations of state.

HYDROCODE: Same as Houser, Centrella & Smith (1994)

MODEL(s): Axisymmetric, n' = 0 rotation law; a total of 3 models examined with, respectively, (n, T/|W|) = (0.5, 0.31), (1.0, 0.32), (1.5, 0.32) — see their Table II. The SCF-based method used to construct these models and the (rather contrived) method used to map the continuum distribution of mass into the Lagrangian-based SPH code is discussed in detail in §III of the paper.

Break 1

Equilibrium Models constucted and examined by

TIPD98

<math>~n</math> <math>~n'</math> <math>~T/|W|</math> <math>~\Omega_p/\Omega_\mathrm{eq}</math> <math>~R_p/R_\mathrm{eq}</math>
<math>~\tfrac{5}{4}</math> 0 0.274 2.44 0.258
0.279 2.45 0.250
0.288 2.49 0.234
0.319 2.65 0.187
0.346 2.88 0.147
0.363 3.11 0.123
<math>~\tfrac{3}{2}</math> 0.266 2.97 0.258
0.280 3.03 0.234
0.290 3.09 0.218
0.300 3.16 0.202
0.316 3.28 0.179
0.327 3.35 0.163
0.344 3.61 0.139
0.368 4.04 0.107
0.396 4.79 0.0754
1 0.262 6.85 0.139
0.283 8.90 0.107
0.292 9.81 0.0958
0.311 12.4   0.0754
<math>~\tfrac{3}{2}</math> 0.272 15.1   0.0754
0.280 16.9   0.0675
0.290 19.0   0.0595
<math>~\tfrac{5}{2}</math> 0 0.262 7.19 0.202
0.268 7.26 0.194
0.273 7.34 0.187
0.285 7.52 0.171
0.297 7.73 0.155
0.304 7.85 0.147
0.324 8.29 0.123
0.338 8.68 0.107
0.371 9.85 0.0754
0.389 10.8   0.0595

Break 2

 

Nonlinear growth of the bar-mode deformation is studied for a differentially rotating star with supercritical rotational energy. In particular, the growth mechanism of some azimuthal modes with odd wave numbers is examined … Mode coupling to even modes, i.e., the bar mode and higher harmonics, significantly enhances the amplitudes of odd modes …

HYDROCODE: Newtonian, 3D Eulerian, Cartesian, with entropy tracer; reflection symmetry through equatorial plane; <math>~\Gamma=2</math>; Poisson solved with preconditioned conjugate gradient (PCG) method

MODEL(s): axisymmetric, n = 1 polytrope; j-constant rotation law with A = 1; their Table I lists four different equilibrium configurations having T/|W| = 0.256, 0.268, 0.277, 0.281.

Additional references identified through the above set of references:

  • M. Saijo (2018), Phys. Rev. D, 98, 024003: Determining the stiffness of the equation of state using low T/W dynamical instabilities in differentially rotating stars
 

We investigate the nature of low T/W dynamical instabilities in various ranges of the stiffness of the equation of state in differentially rotating stars … We analyze these instabilities in both a linear perturbation analysis and a three-dimensional hydrodynamical simulation … the nature of the eigenfunction that oscillates between corotation and the surface for an unstable star requires reinterpretation of pulsation modes in differentially rotating stars.

See Also


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