Difference between revisions of "User:Tohline/AxisymmetricConfigurations/HSCF"

From VistrailsWiki
Jump to navigation Jump to search
 
(22 intermediate revisions by the same user not shown)
Line 1: Line 1:
__FORCETOC__  
__FORCETOC__  
<!-- __NOTOC__ will force TOC off -->
<!-- __NOTOC__ will force TOC off -->
=Hachisu Self-Consistent-Field Technique=
=Introduction=
 
 
{{LSU_HBook_header}}
 
Text colored dark green has been extracted ''verbatim'' from a relatively early review by [http://adsabs.harvard.edu/abs/1969ARA%26A...7..665S P. A. Strittmatter (1969)]: &nbsp;<font color="darkgreen">The study of the internal structure of rotating stars is complicated by the necessity to introduce circulation currents and to solve Poisson's equation in two dimensions with an unknown boundary to the density distribution &hellip; With the advent of large-capacity, fast computers it has become possible to tackle the problem by direct numerical integration of the appropriate partial-differential equations.  This method has been adopted by [http://adsabs.harvard.edu/abs/1964ApJ...140..552J R. A. James (1964)] in considering uniformly rotating polytopes and by [http://adsabs.harvard.edu/abs/1965ApJ...142..208S R. Stoeckly (1965)] who treated the case when the rotation velocity <math>~\Omega</math> varies as a Gaussian function of the distance from the rotation axis &hellip; Ostriker and his collaborators (cf. [http://adsabs.harvard.edu/abs/1968ApJ...151.1075O J. P. Ostriker &amp; J. W.-K. Mark 1968]) have developed an approach (called, by them, the self-consistent field or SCF method) in which Poisson's equation is replaced by its formal integral solution and an iterative procedure established in which the potential is derived from a guessed density distribution.  A new density distribution is then obtained from the equation of hydrostatic support &hellip; While the present author has no experience in using the SCF method, it would appear to combine considerable power with reasonable speed of application and must probably be considered the most practical method of solution so far developed.</font>
 
 
<div align="center">
<table border="1" cellpadding="8" align="center" width="80%">
<tr><th align="center" colspan="2"><font size="+0">Table 1: &nbsp;Poisson Equation</font></th></tr>
<tr>
  <th align="center">Integral Relation</th>
  <th align="center">Differential Relation </th>
</tr>
<tr>
  <td align="center">
<table border="0" align="center">
<tr>
  <td align="right">
<math>~ \Phi(\vec{x})</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ -G \int \frac{\rho(\vec{x}^{~'})}{|\vec{x}^{~'} - \vec{x}|} d^3x^' \, .</math>
  </td>
</tr>
</table>
  </td>
  <td align="center">
{{User:Tohline/Math/EQ_Poisson01}}
  </td>
</tr>
</table>
</div>
 
 
In this paragraph, teal-colored text has been extracted ''verbatim'' from &sect;II.b of [http://adsabs.harvard.edu/abs/1968ApJ...151.1075O J. P. Ostriker and J. W.-K. Mark (1968, ApJ, 151, 1075 - 1088)]: &nbsp;<font color="#009999">&hellip;  we shall alternately solve each of the two problems as exactly as numerical techniques permit, and then iterate to obtain self-consistency. 
 
<div align="center">
i) &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; <math>~\rho(\vec{x})</math>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; <math>~\xrightarrow{~~~~~\mathrm{potential ~theory}~~~~~~} ~ \Phi(\vec{x})</math><br />
ii) &nbsp; &nbsp;<math>~\Phi(\vec{x})+\Psi(\varpi) ~\xrightarrow{~~~\mathrm{equilibrium ~condition}~~} ~ \rho(\vec{x})</math>
</div>


<table border="0" cellpadding="10" align="left"><tr><th>[[Image:Hai.gif|90px|center|Izumi Hachisu]]<br />I. Hachisu</th></tr></table> In 1986, Izumi Hachisu published two papers in '''The Astrophysical Journal Supplement Series''' ([http://adsabs.harvard.edu/abs/1986ApJS...61..479H vol. 61, pp. 479-507], and [http://adsabs.harvard.edu/abs/1986ApJS...62..461H vol. 62, pp. 461-499]) describing ''"A Versatile Method for Obtaining Structures of Rapidly Rotating Stars."''  (Henceforth, we will refer to this method as the Hachisu Self-Consistent-Field, or HSCF, technique.) We have found the HSCF technique to be an extremely powerful tool for constructing equilibrium configurations of self-gravitating fluid systems under a wide variety of different circumstances.  This chapter has been built upon [http://www.phys.lsu.edu/astro/H_Book.current/Applications/Structure/HSCF_Code/HSCF.outline.html an (''ca.'' 1999) outline of the HSCF technique] that appeared in our [http://www.phys.lsu.edu/astro/H_Book.current/H_Book.html original version of this HyperText Book] (H_Book).  The photo of Professor Izumi Hachisu shown here, on the left, dates from the mid-1980s &#8212; about the time he developed this remarkably useful numerical technique; a more recent photo can be found on the [http://gpes.c.u-tokyo.ac.jp/faculty-staff/materials-systems-and-dynamics/izumi-hachisu.html web page associated with Professor Hachisu's current faculty appointment at the University of Tokyo, Komaba].
The integral relation between potential and density</font> &#8212; as opposed to the ''differential relation'', see Table 1 &#8212; <font color="#009999">encourages us to think that a "bad" density will, in step (i), lead to a "good" potential, which, by step (ii), will yield a "good" density which when substituted in (i) produces a "very good" potential, etc.  The method works.</font> It should be noted that <font color="#009999">an SCF approach is commonly used for determining molecular structure under the name [https://en.wikipedia.org/wiki/Hartree–Fock_method Hartree-Fock]</font>.


=Hachisu Self-Consistent-Field Technique=


{{LSU_HBook_header}}
<table border="0" cellpadding="10" align="left"><tr><th>[[Image:Hai.gif|90px|center|Izumi Hachisu]]<br />I. Hachisu</th></tr></table>
In 1986, Izumi Hachisu published two papers in '''The Astrophysical Journal Supplement Series''' ([http://adsabs.harvard.edu/abs/1986ApJS...61..479H vol. 61, pp. 479-507], and [http://adsabs.harvard.edu/abs/1986ApJS...62..461H vol. 62, pp. 461-499]) describing ''"A Versatile Method for Obtaining Structures of Rapidly Rotating Stars."''  (Henceforth, we will refer to this method as the Hachisu Self-Consistent-Field, or HSCF, technique.)  The HSCF technique has been built upon &#8212; and has further improved &#8212; the SCF method introduced by [http://adsabs.harvard.edu/abs/1968ApJ...151.1075O J. P. Ostriker &amp; J. W.-K. Mark (1968)].  We have found the HSCF technique to be an extremely powerful tool for constructing equilibrium configurations of self-gravitating fluid systems under a wide variety of different circumstances. 
<!-- This chapter has been built upon [http://www.phys.lsu.edu/astro/H_Book.current/Applications/Structure/HSCF_Code/HSCF.outline.html an (''ca.'' 1999) outline of the HSCF technique] that appeared in our [http://www.phys.lsu.edu/astro/H_Book.current/H_Book.html original version of this HyperText Book] (H_Book).  -->
The photo of Professor Izumi Hachisu shown here, on the left, dates from the mid-1980s &#8212; about the time he developed this remarkably useful numerical technique; a more recent photo can be found on the [http://gpes.c.u-tokyo.ac.jp/faculty-staff/materials-systems-and-dynamics/izumi-hachisu.html web page associated with Professor Hachisu's current faculty appointment at the University of Tokyo, Komaba].


==Constructing Two-Dimensional, Axisymmetric Structures==
==Constructing Two-Dimensional, Axisymmetric Structures==
As has been explained in [[User:Tohline/AxisymmetricConfigurations/SolutionStrategies#Technique|an accompanying discussion]], our objective is to solve an ''algebraic'' expression for hydrostatic balance,
As has been explained in [[User:Tohline/AxisymmetricConfigurations/SolutionStrategies#Technique|an accompanying discussion]], our objective is to solve an ''algebraic'' expression for hydrostatic balance,
<div align="center">
<div align="center">
<math>~H + \Phi + \Psi = C_\mathrm{B}</math> ,
<math>~H + \Phi + \Psi = C_0</math> ,
</div>
</div>


Line 22: Line 70:


===Steps to Follow===
===Steps to Follow===
 
<table border="0" width="33%" align="right" cellpadding="5"><tr><td align="center">
<table border="1" cellpadding="5" align="center">
<tr><td align="center">
Figure 1 extracted from [http://adsabs.harvard.edu/abs/1986ApJS...61..479H I. Hachisu (1986)]<p></p>
"''A Versatile Method for Obtaining Structures of Rapidly Rotating Stars''"<p></p>
ApJS, vol. 61, pp. 479-507 &copy; American Astronomical Society
</td></tr>
<tr><td>
[[File:HSCF MeridionalPlane.png|center|250px|HSCF Meridional Plane]]</td></tr>
<tr><td align="left">Two points A and B on the surface specify a model uniquely: (a) spheroidal structure; (b) toroidal structure.</td></tr></table>
</td></tr></table>
<ol>
<ol>
   <li>'''Choose a particular [[User:Tohline/SR#Barotropic_Structure|barotropic equation of state]].''' &nbsp; More specifically, functionally define the density-enthalpy relationship, <math>~\rho(H)</math>, and identify what value, <math>~H_\mathrm{surface}</math>, the enthalpy will have at the surface of your configuration.  For example, if a ''polytropic'' equation of state is adopted, <math>~H_\mathrm{surface} = 0</math> is a physically reasonable prescription.</li>
   <li>'''Choose a particular [[User:Tohline/SR#Barotropic_Structure|barotropic equation of state]].''' &nbsp; More specifically, functionally define the density-enthalpy relationship, <math>~\rho(H)</math>, and identify what value, <math>~H_\mathrm{surface}</math>, the enthalpy will have at the surface of your configuration.  For example, if a ''polytropic'' equation of state is adopted, <math>~H_\mathrm{surface} = 0</math> is a physically reasonable prescription.</li>
   <li>Choosing from, for example, a list of astrophysically relevant ''[[User:Tohline/AxisymmetricConfigurations/SolutionStrategies#Simple_Rotation_Profile_and_Centrifugal_Potential|simple rotation profiles]],'' '''specify the corresponding functional form of the centrifugal potential''', <math>~\Psi(\varpi)</math>, that will define the radial distribution of specific angular momentum in your equilibrium configuration.</li>
   <li>Choosing from, for example, a list of astrophysically relevant ''[[User:Tohline/AxisymmetricConfigurations/SolutionStrategies#Simple_Rotation_Profile_and_Centrifugal_Potential|simple rotation profiles]],'' '''specify the corresponding functional form of the centrifugal potential''', <math>~\Psi(\varpi)</math>, that will define the radial distribution of specific angular momentum in your equilibrium configuration.  If the choice is uniform rotation, then <math>~\Psi = - \varpi^2 \omega_0^2/2 \, ,</math> where <math>~\omega_0</math> is a constant to be determined.</li>
   <li>On your chosen computational lattice &#8212; for example, on a cylindrical-coordinate mesh &#8212; '''identify two boundary points''', A and B, that will lie on the surface of your equilibrium configuration.  These two points should remain fixed in space during the HSCF iteration cycle and ultimately will confine the volume and define the geometry of the derived equilibrium object.  Note that, by definition, the enthalpy at these two points is, <math>~H_A = H_B = H_\mathrm{surface}</math>.</li>
   <li>On your chosen computational lattice &#8212; for example, on a cylindrical-coordinate mesh &#8212; '''identify two boundary points''', A and B, that will lie on the surface of your equilibrium configuration.  These two points should remain fixed in space during the HSCF iteration cycle and ultimately will confine the volume and define the geometry of the derived equilibrium object.  Note that, by definition, the enthalpy at these two points is, <math>~H_A = H_B = H_\mathrm{surface}</math>.</li>
   <li>Throughout the volume of your computational lattice, ''' ''guess'' a trial distribution of the mass density,''' <math>~\rho(\varpi,z)</math>, such that no material falls outside a volume defined by the two boundary points, A and B, that were identified in Step #3.  Usually an initially uniform density distribution will suffice to start the SCF iteration.</li>
   <li>Throughout the volume of your computational lattice, ''' ''guess'' a trial distribution of the mass density,''' <math>~\rho(\varpi,z)</math>, such that no material falls outside a volume defined by the two boundary points, A and B, that were identified in Step #3.  Usually an initially uniform density distribution will suffice to start the SCF iteration.</li>
   <li>Via ''some'' accurate numerical algorithm, '''solve the Poisson equation''' to determine the gravitational potential, <math>~\Phi(\varpi,z)</math>, throughout the computation lattice that corresponds to the trial mass-density distribution that was specified in Step #4 (or in Step #9).</li>
   <li>Via ''some'' accurate numerical algorithm, '''solve the Poisson equation''' to determine the gravitational potential, <math>~\Phi(\varpi,z)</math>, throughout the computational lattice corresponding to the trial mass-density distribution that was specified in Step #4 (or in Step #9).</li>
  <li>From the gravitational potential determined in Step #5, '''identify the values of <math>~\Phi_A</math> and <math>~\Phi_B</math>''' at the two boundary points that were selected in Step #3.</li>
  <li>From the "known" values of the enthalpy (Step #3) and the gravitational potential (Step #6) at the two selected surface boundary points A and B, '''determine the values of the constants, <math>~C_0</math> and <math>~\omega_0</math>,''' that appear in the algebraic equation that defines hydrostatic equilibrium.</li>
  <li>From the most recently determined values of the gravitational potential, <math>~\Phi(\varpi,z)</math> (Step #5), and the values of the two constants, <math>~C_0</math> and <math>~\omega_0</math> just determined (Step #7), '''determine the enthalpy distribution throughout the computational lattice.'''</li>
  <li>From <math>~H(\varpi,z)</math> and the selected barotropic equation of state (Step #1), '''calculate an "improved guess" of the density distribution,''' <math>~\rho(\varpi,z)</math>, throughout the computational lattice.</li>
  <li>'''Has the model converged to a satisfactory equilibrium solution?'''  (Usually a satisfactory solution has been achieved when the derived model parameters &#8212; for example, the values of <math>~C_0</math> and <math>~\omega_0</math> &#8212; change very little between successive iterations and the viral error is sufficiently small.)
  <ul>
<li>If the answer is, "NO": &nbsp; Repeat steps 5 through 10.</li>
<li>If the answer is, "YES": &nbsp; Stop iteration.</li>
  </ul>
  </li>
</ol>
</ol>


=Related Discussions=
=Related Discussions=
* [[User:Tohline/SSC/Structure/BiPolytropes#BiPolytropes|Constructing BiPolytropes]]
==Reviews==
* [[User:Tohline/SSC/Structure/BiPolytropes/Analytic5_1#BiPolytrope_with_nc_.3D_5_and_ne_.3D_1|Analytic description of BiPolytrope with <math>(n_c, n_e) = (5,1)</math>]]
 
* [[User:Tohline/SSC/Structure/BonnorEbert#Pressure-Bounded_Isothermal_Sphere|Bonnor-Ebert spheres]]
* [http://adsabs.harvard.edu/abs/1969ARA%26A...7..665S P. A. Strittmatter (1969, Annual Review of Astronomy and Astrophysics, 7, 665 - 684)] &#8212; ''Stellar Rotation''
** [http://en.wikipedia.org/wiki/Bonnor-Ebert_mass Bonnor-Ebert Mass] according to Wikipedia
* [http://adsabs.harvard.edu/abs/1967ARA%26A...5..465L N. R. Lebovitz (1967, Annual Review of Astronomy and Astrophysics, 5, 465 - 480)] &#8212; ''Rotating Fluid Masses''
** [http://www.astro.umd.edu/~cychen/MATLAB/ASTR320/matlabFrom320spring2011/Bonnor-EbertSphere/html/BonnorEbert.html A MATLAB script to determine the Bonnor-Ebert Mass coefficient] developed by [http://www.astro.umd.edu/people/cychen.html Che-Yu Chen] as a graduate student in the University of Maryland Department of Astronomy
 
* [[User:Tohline/SSC/Structure/LimitingMasses#Sch.C3.B6nberg-Chandrasekhar_Mass|Sch&ouml;nberg-Chandrasekhar limiting mass]]
==Solution Methods==
* [[User:Tohline/SSC/Structure/LimitingMasses#Relationship_Between_the_Bonnor-Ebert_and_Sch.C3.B6nberg-Chandrasekhar_Critical_Masses|Relationship between Bonnor-Ebert and Sch&ouml;nberg-Chandrasekhar limiting masses]]
 
* [http://adsabs.harvard.edu/abs/1985A%26A...146..260E Y. Eriguchi &amp; E. Mueller (1985, A&amp;A, 146, 260 - 268)] &#8212; ''A General Computational Method for Obtaining Equilibria of Self-Gravitating and Rotating Gases''
* [http://adsabs.harvard.edu/abs/1983ApJ...268..155S S. W. Stahler (1983, ApJ, 268, 155 - 184)] &#8212; ''The Equilibria of Rotating, Isothermal Clouds. I. - Method of Solution''
* [http://adsabs.harvard.edu/abs/1978PASJ...30..507E Y. Eriguchi (1978, PASJ, 30, 507 - 518)] &#8212; ''Hydrostatic Equilibria of Rotating Polytropes''
* [http://adsabs.harvard.edu/abs/1975SvA....19..151B S. I. Blinnikov (1975, Soviet Astronomy, 19, 151 - 156)] &#8212; ''Self-Consistent Field Method in the Theory of Rotating Stars''
* [http://adsabs.harvard.edu/abs/1974ApJ...194..709C M. J. Clement (1974, ApJ, 194, 709 - 714)] &#8212; ''On the Solution of Poisson's Equation for Rapidly Rotating Stars''
* [http://adsabs.harvard.edu/abs/1970ApJ...161..579J S. Jackson (1970, ApJ, 161, 579 - 585)] &#8212; ''Rapidly Rotating Stars. V. The Coupling of the Henyey and the Self-Consistent Methods''
* [http://adsabs.harvard.edu/abs/1968ApJ...151.1075O J. P. Ostriker &amp; J. W.-K. Mark (1968, ApJ, 151, 1075 - 1088)] &#8212; ''Rapidly Rotating Stars. I. The Self-Consistent-Field Method''
* [http://adsabs.harvard.edu/abs/1964ApJ...140..552J R. A. James (1964, ApJ, 140, 552 - 582)] &#8212; ''The Structure and Stability of Rotating Gas Masses''
 
==Early Eriguchi Applications==
 
* [http://adsabs.harvard.edu/abs/1985A%26A...147..161E Y. Eriguchi &amp; E. Mueller (1985, A&amp;A, 147, 161 - 168)] &#8212; ''Equilibrium Models of Differentially Rotating Polytropes and the Collapse of Rotating Stellar Cores''
* [http://adsabs.harvard.edu/abs/1984PASJ...36..497H I. Hachisu &amp; Y. Eriguchi (1984, PASJ, 36, 497 - 503)] &#8212; ''Bifurcation Points on the Maclaurin Sequence''
* [http://adsabs.harvard.edu/abs/1984PASJ...36..259H I. Hachisu &amp; Y. Eriguchi (1984, PASJ, 36, 259 - 276)] &#8212; ''Binary Fluid Star''
* [http://adsabs.harvard.edu/abs/1984PASJ...36..239H I. Hachisu &amp; Y. Eriguchi (1984, PASJ, 36, 239 - 257)] &#8212; ''Fission of Dumbbell Equilibrium and Binary State of Rapidly Rotating Polytropes''
* [http://adsabs.harvard.edu/abs/1983MNRAS.204..583H I. Hachisu &amp; Y. Eriguchi (1983, MNRAS, 204, 583 - 589)] &#8212; ''Bifurcations and Phase Transitions of Self-Gravitating and Uniformly Rotating Fluid''
* [http://adsabs.harvard.edu/abs/1982PThPh..68..206H I. Hachisu &amp; Y. Eriguchi (1982, Prog. Theor. Phys., 68, 206 - 221)] &#8212; ''Bifurcation and Fission of Three Dimensional, Rigidly Rotating and Self-Gravitating Polytropes''
* [http://adsabs.harvard.edu/abs/1982PThPh..68..191H I. Hachisu, Y. Eriguchi, &amp; D. Sugimoto (1982, Prog. Theor. Phys., 68, 191 - 205)] &#8212; ''Rapidly Rotating Polytropes and Concave Hamburger Equilibrium''
* [http://adsabs.harvard.edu/abs/1981PThPh..65.1870E Y. Eriguchi &amp; D. Sugimoto (1981, Prog. Theor. Phys., 65, 1870 - 1875)] &#8212; ''Another Equilibrium Sequence of Self-Gravitating and Rotating Incompressible Fluid''
* [http://adsabs.harvard.edu/abs/1980PThPh..63.1957F T. Fukushima, Y. Eriguchi, D. Sugimoto, G. S. Bisnovatyi-Kogan (1980, Prog. Theor. Phys., 63, 1957 - 1970)] &#8212; ''Concave Hamburger Equilibrium of Rotating Bodies''
 
==Other Example Applications==
 
* [http://adsabs.harvard.edu/abs/1985A%26A...152..325M E. Mueller &amp; Y. Eriguchi (1985, A&amp;A, 152, 325 - 335)] &#8212; ''Equilibrium Models of Differentially Rotating, Completely Catalyzed, Zero-Temperature Configurations With Central Densities Intermediate to White Dwarf and Neutron Star Densities''
* [http://adsabs.harvard.edu/abs/1981ApJ...250..362I J. R. Ipser &amp; R. A. Managan (1981, ApJ, 250, 362 - 372)] &#8212; ''On the Existence and Structure of Inhomogeneous Analogs of the Dedekind and Jacobi Ellipsoids''
* [http://adsabs.harvard.edu/abs/1977ApJ...211..568C C. T. Cunningham (1977, ApJ, 211, 568 - 578)] &#8212; ''Rapidly Rotating Spheroids of Polytropic Index n = 1''
* [http://adsabs.harvard.edu/abs/1975ApJ...199..179D R. H. Durisen (1975, ApJ, 199, 179 - 183)] &#8212; ''Upper Mass Limits for Stable Rotating White Dwarfs''
* [http://adsabs.harvard.edu/abs/1973ApJ...180..159B P. Bodenheimer &amp; J. P. Ostriker (1973, ApJ, 180, 159 - 170)] &#8212; ''Rapidly Rotating Stars.  VIII. Zero-Viscosity Polytropic Sequences'
* [http://adsabs.harvard.edu/abs/1971ApJ...167..153B P. Bodenheimer (1971, ApJ, 167, 153 - 163)] &#8212; ''Rapidly Rotating Stars.  VII. Effects of Angular Momentum on Upper-Main-Sequence Models''
* [http://adsabs.harvard.edu/abs/1970ApJ...161.1101B P. Bodenheimer &amp; J. P. Ostriker (1970, ApJ, 161, 1101 - 1113)] &#8212; ''Rapidly Rotating Stars.  VI. Pre-Main-Sequence Evolution of Massive Stars''
* R. Kippenhahn &amp; H.-C. Thomas (1970) in [http://adsabs.harvard.edu/abs/1970IAUCo...4.....S Proceedings of the 4th IAU Colloquium], held at the Ohio State University, Columbus, Ohio, September 8 - 11, 1969, Dordrecht:  Riedel Publishing Co., edited by A. Slettebak &#8212; ''Stellar Rotation'' ==> Purchase proceedings from [http://www.springer.com/us/book/9789027701565 Springer], from [https://trove.nla.gov.au/work/21340423?selectedversion=NBD179107 Australia], or from [https://books.google.com/books/about/Stellar_Rotation.html?id=4NnsCAAAQBAJ Google]
** In the introductory section of his paper, [http://adsabs.harvard.edu/abs/1970ApJ...161..579J S. Jackson (1970)] references this article by Kippenhahn &amp; Thomas in the context of uniformly rotating, and therefore only mildly distorted, structures.
* [http://adsabs.harvard.edu/abs/1969ApJ...155..987O J. P. Ostriker &amp; J. L. Tassoul (1969, ApJ, 155, 987 - )] &#8212; ''On the Oscillations and Stability of Rotating Stellar Models. II. Rapidly Rotating White Dwarfs''
* [http://adsabs.harvard.edu/abs/1969ApJ...156.1051C M. J. Clement (1969, ApJ, 156, 1051 - 1068)] &#8212; ''Differential Rotation in Stars on the Upper Main Sequence''
* [http://adsabs.harvard.edu/abs/1968ApJ...154..627M J. W.-K. Mark (1968, ApJ, 154, 627 - )] &#8212; ''Rapidly Rotating Stars. III. Massive Main-Sequence Stars''
* [http://adsabs.harvard.edu/abs/1968ApJ...151.1089O J. P. Ostriker &amp; P. Bodenheimer (1968, ApJ, 151, 1089 - )] &#8212; ''Rapidly Rotating Stars. II. Massive White Dwarfs''
* [http://adsabs.harvard.edu/abs/1968ApJ...151..203F J. Faulkner, I. W. Roxburgh, &amp; P. A. Strittmatter (1968, ApJ, 151, 203 - 216)] &#8212; ''Uniformly Rotating Main-Sequence Stars''
* [http://adsabs.harvard.edu/abs/1965ApJ...142..208S R. Stoeckly (1965, ApJ, 142, 208 - 228)] &#8212; ''Polytropic Models with Fast, Non-Uniform Rotation''
** In the introductory section of his paper, [http://adsabs.harvard.edu/abs/1970ApJ...161..579J S. Jackson (1970)] states that a differentially rotating polytropic structure with a rotationally induced ''extreme distortion'' was first illustrated in this article by Stoeckly.
 
==Henyey Technique for Nonrotating Stars==
* [http://adsabs.harvard.edu/abs/1959ApJ...129..628H L. G. Henyey, L. Wilets, K. H. B&ouml;hm, R. Lelevier, &amp; R. D. Levee (1959, ApJ, 129, 628 - )] &#8212; ''A Method for Automatic Computation of Stellar Evolution''
* [http://adsabs.harvard.edu/abs/1964ApJ...139..306H L. G. Henyey, J. E. Forbes, &amp; N. L. Gould (1964, ApJ, 139, 306 - )] &#8212; ''A New Method of Automatic Computation of Stellar Evolution''


* Wikipedia introduction to the [http://en.wikipedia.org/wiki/Lane-Emden_equation Lane-Emden equation]
* Wikipedia introduction to [http://en.wikipedia.org/wiki/Polytrope Polytropes]


{{LSU_HBook_footer}}
{{LSU_HBook_footer}}

Latest revision as of 22:13, 2 May 2018

Introduction

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

Text colored dark green has been extracted verbatim from a relatively early review by P. A. Strittmatter (1969):  The study of the internal structure of rotating stars is complicated by the necessity to introduce circulation currents and to solve Poisson's equation in two dimensions with an unknown boundary to the density distribution … With the advent of large-capacity, fast computers it has become possible to tackle the problem by direct numerical integration of the appropriate partial-differential equations. This method has been adopted by R. A. James (1964) in considering uniformly rotating polytopes and by R. Stoeckly (1965) who treated the case when the rotation velocity <math>~\Omega</math> varies as a Gaussian function of the distance from the rotation axis … Ostriker and his collaborators (cf. J. P. Ostriker & J. W.-K. Mark 1968) have developed an approach (called, by them, the self-consistent field or SCF method) in which Poisson's equation is replaced by its formal integral solution and an iterative procedure established in which the potential is derived from a guessed density distribution. A new density distribution is then obtained from the equation of hydrostatic support … While the present author has no experience in using the SCF method, it would appear to combine considerable power with reasonable speed of application and must probably be considered the most practical method of solution so far developed.


Table 1:  Poisson Equation
Integral Relation Differential Relation

<math>~ \Phi(\vec{x})</math>

<math>~=</math>

<math>~ -G \int \frac{\rho(\vec{x}^{~'})}{|\vec{x}^{~'} - \vec{x}|} d^3x^' \, .</math>

LSU Key.png

<math>\nabla^2 \Phi = 4\pi G \rho</math>


In this paragraph, teal-colored text has been extracted verbatim from §II.b of J. P. Ostriker and J. W.-K. Mark (1968, ApJ, 151, 1075 - 1088):  … we shall alternately solve each of the two problems as exactly as numerical techniques permit, and then iterate to obtain self-consistency.

i)             <math>~\rho(\vec{x})</math>          <math>~\xrightarrow{~~~~~\mathrm{potential ~theory}~~~~~~} ~ \Phi(\vec{x})</math>
ii)    <math>~\Phi(\vec{x})+\Psi(\varpi) ~\xrightarrow{~~~\mathrm{equilibrium ~condition}~~} ~ \rho(\vec{x})</math>

The integral relation between potential and density — as opposed to the differential relation, see Table 1 — encourages us to think that a "bad" density will, in step (i), lead to a "good" potential, which, by step (ii), will yield a "good" density which when substituted in (i) produces a "very good" potential, etc. The method works. It should be noted that an SCF approach is commonly used for determining molecular structure under the name Hartree-Fock.

Hachisu Self-Consistent-Field Technique

Izumi Hachisu

I. Hachisu

In 1986, Izumi Hachisu published two papers in The Astrophysical Journal Supplement Series (vol. 61, pp. 479-507, and vol. 62, pp. 461-499) describing "A Versatile Method for Obtaining Structures of Rapidly Rotating Stars." (Henceforth, we will refer to this method as the Hachisu Self-Consistent-Field, or HSCF, technique.) The HSCF technique has been built upon — and has further improved — the SCF method introduced by J. P. Ostriker & J. W.-K. Mark (1968). We have found the HSCF technique to be an extremely powerful tool for constructing equilibrium configurations of self-gravitating fluid systems under a wide variety of different circumstances. The photo of Professor Izumi Hachisu shown here, on the left, dates from the mid-1980s — about the time he developed this remarkably useful numerical technique; a more recent photo can be found on the web page associated with Professor Hachisu's current faculty appointment at the University of Tokyo, Komaba.

Constructing Two-Dimensional, Axisymmetric Structures

As has been explained in an accompanying discussion, our objective is to solve an algebraic expression for hydrostatic balance,

<math>~H + \Phi + \Psi = C_0</math> ,

in conjunction with the Poisson equation in a form that is appropriate for two-dimensional, axisymmetric systems, namely,

<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>

Steps to Follow

Figure 1 extracted from I. Hachisu (1986)

"A Versatile Method for Obtaining Structures of Rapidly Rotating Stars"

ApJS, vol. 61, pp. 479-507 © American Astronomical Society

HSCF Meridional Plane
Two points A and B on the surface specify a model uniquely: (a) spheroidal structure; (b) toroidal structure.
  1. Choose a particular barotropic equation of state.   More specifically, functionally define the density-enthalpy relationship, <math>~\rho(H)</math>, and identify what value, <math>~H_\mathrm{surface}</math>, the enthalpy will have at the surface of your configuration. For example, if a polytropic equation of state is adopted, <math>~H_\mathrm{surface} = 0</math> is a physically reasonable prescription.
  2. Choosing from, for example, a list of astrophysically relevant simple rotation profiles, specify the corresponding functional form of the centrifugal potential, <math>~\Psi(\varpi)</math>, that will define the radial distribution of specific angular momentum in your equilibrium configuration. If the choice is uniform rotation, then <math>~\Psi = - \varpi^2 \omega_0^2/2 \, ,</math> where <math>~\omega_0</math> is a constant to be determined.
  3. On your chosen computational lattice — for example, on a cylindrical-coordinate mesh — identify two boundary points, A and B, that will lie on the surface of your equilibrium configuration. These two points should remain fixed in space during the HSCF iteration cycle and ultimately will confine the volume and define the geometry of the derived equilibrium object. Note that, by definition, the enthalpy at these two points is, <math>~H_A = H_B = H_\mathrm{surface}</math>.
  4. Throughout the volume of your computational lattice, guess a trial distribution of the mass density, <math>~\rho(\varpi,z)</math>, such that no material falls outside a volume defined by the two boundary points, A and B, that were identified in Step #3. Usually an initially uniform density distribution will suffice to start the SCF iteration.
  5. Via some accurate numerical algorithm, solve the Poisson equation to determine the gravitational potential, <math>~\Phi(\varpi,z)</math>, throughout the computational lattice corresponding to the trial mass-density distribution that was specified in Step #4 (or in Step #9).
  6. From the gravitational potential determined in Step #5, identify the values of <math>~\Phi_A</math> and <math>~\Phi_B</math> at the two boundary points that were selected in Step #3.
  7. From the "known" values of the enthalpy (Step #3) and the gravitational potential (Step #6) at the two selected surface boundary points A and B, determine the values of the constants, <math>~C_0</math> and <math>~\omega_0</math>, that appear in the algebraic equation that defines hydrostatic equilibrium.
  8. From the most recently determined values of the gravitational potential, <math>~\Phi(\varpi,z)</math> (Step #5), and the values of the two constants, <math>~C_0</math> and <math>~\omega_0</math> just determined (Step #7), determine the enthalpy distribution throughout the computational lattice.
  9. From <math>~H(\varpi,z)</math> and the selected barotropic equation of state (Step #1), calculate an "improved guess" of the density distribution, <math>~\rho(\varpi,z)</math>, throughout the computational lattice.
  10. Has the model converged to a satisfactory equilibrium solution? (Usually a satisfactory solution has been achieved when the derived model parameters — for example, the values of <math>~C_0</math> and <math>~\omega_0</math> — change very little between successive iterations and the viral error is sufficiently small.)
    • If the answer is, "NO":   Repeat steps 5 through 10.
    • If the answer is, "YES":   Stop iteration.

Related Discussions

Reviews

Solution Methods

Early Eriguchi Applications

Other Example Applications

Henyey Technique for Nonrotating Stars


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