User:Tohline/AxisymmetricConfigurations/PoissonEq
Solving the (Multi-dimensional) Poisson Equation Numerically
|
|---|
| | Tiled Menu | Tables of Content | Banner Video | Tohline Home Page | |
Overview
The set of Principal Governing Equations that serves as the foundation of our study of the structure, stability, and dynamical evolution of self-gravitating fluids contains an equation of motion (the Euler equation) that includes an acceleration due to local gradients in the (Newtonian) gravitational potential, <math>~\Phi</math>. As has been pointed out in an accompanying chapter that discusses the origin of the Poisson equation, the mathematical definition of this acceleration is fundamentally drawn from Isaac Newton's inverse-square law of gravitation, but takes into account that our fluid systems are not point-mass sources but, rather, are represented by a continuous distribution of mass via the function, <math>~\rho(\vec{x},t)</math>. As indicated, in our study, <math>~\rho</math> may depend on time as well as space. The acceleration felt at any point in space is obtained by integrating over the accelerations exerted by each differential mass element. As has been explicitly demonstrated in, respectively, Step 1 and Step 3 of the same accompanying chapter, at any point in time the spatial variation of the gravitational potential, <math>~\Phi(\vec{x})</math>, may be determined from <math>~\rho(\vec{x})</math> via either an integral or a differential equation as follows:
| Table 1: Poisson Equation | ||||||
|---|---|---|---|---|---|---|
| Integral Representation | Differential Representation | |||||
|
|
|||||
While it is possible in some restricted situations to determine analytic expressions for the matched pair of functions, <math>~\Phi </math> and <math>~\rho</math>, that satisfy the Poisson equation, modeling the vast majority of interesting astrophysical problems requires the develop of a numerical scheme to solve the Poisson equation. In what follows, our aim is twofold: (a) To recount — in a reasonable amount of detail — the steps that we have taken over the past, approximately forty years to develop more and more accurate and efficient ways to solve the Poisson equation in full three-dimensional generality; and (b) To list, if not summarize, alternative techniques that have been successfully employed by other research groups over the years.
|
|
Dimensionality: | 2D or 3D | ||
|
|
Green's Function for Integral Representation: | [Sph]erical [Tor]oidal Gridless ([Lag]rangian) | ||
|
|
Computational Mesh Used for Differential Representation: | [Car]tesian [Cyl]indrical [Sph]erical Gridless ([Lag]rangian) | ||
Our Approach — Initially Heavily Influenced by Black & Bodenheimer (1975)
|
D. C. Black & P. Bodenheimer (1975, ApJ, 199, 619 - 632) published a detailed description of the numerical techniques that they implemented in order to model, in two dimensions, the dynamical collapse of axisymmetric interstellar gas clouds. Among the techniques was their method of solving, in cylindrical coordinates, the two-dimensional (2D) Poisson equation. As is acknowledged in the last paragraph of their paper, Black & Bodenheimer were introduced "to the pitfalls of two-dimensional hydrodynamics" by Drs. J. LeBlanc and J. Wilson who, at the time, were both staff scientists at the Lawrence Livermore Laboratory. In 1976, having completed my formal graduate-level course work in the astronomy program at UC, Santa Cruz, Peter Bodenheimer (on the faculty of Lick Observatory and UC, Santa Cruz) asked me if I would be interested in working with him and David Black (a planetary scientist at NASA, Ames Research Center) on the development of a fully three-dimensional (3D) hydrodynamics code to study, not only the collapse, but also the fragmentation of self-gravitating gas clouds. I jumped at the opportunity. As a result, an effort to model the process of spontaneous cloud fragmentation became the focus of my doctoral dissertation research.
|
Borrowing from the Black & Bodenheimer (1975) work, I decided to solve the set of 3D principal governing equations on a cylindrical coordinate <math>~(\varpi, \theta, z)</math> mesh:
- whose outermost radial and vertical boundaries were placed entirely outside of the mass distribution (by way of illustration, see the green-dashed lines in Figure 1);
- that exhibited reflection symmetry through the <math>~(z = 0)</math> equatorial plane ;
- that allowed for non-uniform (logarithmically stretched) grid spacing in both the radial and vertical directions; and,
- (extending the Black & Bodenheimer work from 2D to 3D) with strictly uniform spacing in the azimuthal-coordinate direction.
Mesh Choice: In retrospect, I am still comfortable with the choice that was made, at the time, of a cylindrical coordinate mesh with uniform spacing in the azimuthal-coordinate direction. When expressed in terms of cylindrical coordinates — as opposed to cartesian coordinates, for example — the azimuthal component of the Euler equation offers a natural means by which the conservation of angular momentum can be monitored, if not enforced. And, by enabling the implementation of FFTs — hence, providing the ability to rapidly transform functions, like <math>~\Phi(\theta)</math> and <math>~\rho(\theta)</math>, back and forth between real space and Fourier space — a uniform descritization of the azimuthal grid facilitates the development of an efficient 3D Poisson solver as well as tools to straightforwardly analyze the behavior — e.g., the exponential growth — of individual nonaxisymmetric modes. |
Also, following the lead of Black & Bodenheimer (1975), a hybrid scheme was developed to solve the Poisson equation. As defined in Table 1 above, the integral expression was evaluated at grid locations along the outermost radial and vertical boundaries — illustrated by the dashed green lines displayed in Figure 1 — by integrating over the "pink" mass distribution lying inside of these grid boundaries. Separately, the differential expression was used to evaluate the potential at all interior grid locations; this differential expression was supplemented by the implementation of Neumann boundary conditions (reflection symmetry) along the equatorial plane, and by using the values of the potential just determined along the outermost grid boundaries to provide Dirichlet boundary conditions along those grid boundaries.
Determining Values of the Potential on the Mesh Boundary
Let's determine the potential, <math>~\Phi_B</math>, at all points along the boundary of the cylindrical coordinate mesh by evaluating Table 1's integral expression for the Poisson equation.
Using a Spherical Harmonic Expansion
Following the lead of Black & Bodenheimer (1975), we will insert into this integral relation the Green's function expression for <math>~|\vec{x} - \vec{x}^{~'}|^{-1} </math> as given in terms of Spherical Harmonics, <math>~Y_{\ell m}</math>, which in turn can be written in terms of Associated Legendre Functions. Table 2, below, provides the primary details.
Written in the context of a spherical coordinate system we have,
|
<math>~ \Phi_B(r,\theta,\phi)</math> |
<math>~=</math> |
<math>~ -G \int \frac{1}{|\vec{x}^{~'} - \vec{x}|} ~\rho(\vec{x}^{~'}) d^3x^' </math> |
|
|
<math>~=</math> |
<math>~ -G \int \sum_{\ell=0}^{\infty} \sum_{m=-\ell}^{+\ell} \frac{4\pi}{2\ell+1} \biggl[ \frac{r_<^\ell}{r_>^{\ell+1}} \biggr] Y_{\ell m}^*(\theta^', \phi^') Y_{\ell m}(\theta,\phi) ~\rho(r^', \theta^', \phi^') d^3x^' </math> |
|
|
<math>~=</math> |
<math>~ -4\pi G \sum_{\ell=0}^{\infty} \sum_{m=-\ell}^{+\ell} \frac{Y_{\ell m}(\theta,\phi)}{(2\ell+1)} \biggl[ \frac{1}{r^{\ell+1}}\int_0^r (r^')^\ell Y_{\ell m}^*(\theta^', \phi^') ~\rho(r^', \theta^', \phi^') d^3x^' </math> |
|
|
|
<math>~ + r^\ell \int_r^\infty (r^')^{-(\ell+1)} Y_{\ell m}^*(\theta^', \phi^') ~\rho(r^', \theta^', \phi^') d^3x^' \biggr] \, . </math> |
|
[BT87], p. 66, Eq. (2-122) |
||
If the distance from the origin, <math>~r</math>, of a boundary point (i.e., any point lying along the dashed green lines in Figure 1) is greater than the distance from the origin, <math>~r^'</math>, of all of the (pink) mass elements, the second integral can be completely ignored. We have, then, an expression that we will henceforth refer to as,
|
Form A of the Boundary Potential |
||
|
<math>~ \Phi_B(r,\theta,\phi)</math> |
<math>~=</math> |
<math>~ -4\pi G \sum_{\ell=0}^{\infty} \sum_{m=-\ell}^{+\ell} \frac{Y_{\ell m}(\theta,\phi)}{(2\ell+1)} \biggl[ \frac{1}{r^{\ell+1}}\int_0^r (r^')^\ell Y_{\ell m}^*(\theta^', \phi^') ~\rho(r^', \theta^', \phi^') d^3x^' \biggr] \, . </math> |
|
Jackson (1975), p. 137, Eq. (4.2) |
||
Rewriting this expression for <math>~\Phi_B</math> in terms of cylindrical coordinates — which aligns with our chosen grid coordinate system — and admitting that in practice our summation over the index, <math>~\ell</math>, cannot extend to infinity, we have,
|
<math>~ \Phi_B(\varpi, \phi, z)</math> |
<math>~=</math> |
<math>~ -4\pi G \sum_{\ell=0}^{\ell_\mathrm{max}} \sum_{m=-\ell}^{+\ell} \frac{Y_{\ell m}}{(2\ell+1)} \biggl[ \varpi^2 + z^2 \biggr]^{-(\ell+1)/2} \int Y_{\ell m}^* \biggl[ (\varpi^')^2 + (z^')^2 \biggr]^{\ell/2} ~\rho(\varpi^', \phi^', z^') d^3x^' \, . </math> |
Note that, as a consequence of assuming that our configurations have equatorial-plane symmetry, the weighted integral over the mass distribution necessarily goes to zero anytime the sum of the two indexes, <math>~(\ell + m)</math>, is an odd number. This is because, in each of these situations, the <math>~Y_{\ell m}</math> includes an overall factor of <math>~\cos\theta</math>, which necessarily switches signs between the two hemispheres. After setting <math>~\ell_\mathrm{max}=4</math> — and dropping all terms in the summation for which the index sum, <math>~(\ell + m)</math>, is odd — this expression becomes precisely the relation that was used to determine the boundary values of the gravitational potential in our earliest set of simulations; see, for example, Tohline (1978), Tohline (1980), and Bodenheimer, Tohline, & Black (1980).
Simplification for 2D, Axisymmetric Systems
It is easy to show that this last expression for <math>~\Phi_B</math> — which has been used in our 3D simulations — is a generalization of the expression for <math>~\Phi_B</math> that was employed by Black & Bodenheimer (1975) for 2D, axisymmetric simulations. In axisymmetric systems, by definition, physical variables exhibit no variation in the azimuthal coordinate direction. Hence, in the expression for <math>~\Phi_B</math>:
- the azimuthal coordinate, <math>~\phi</math>, need not appear explicitly as an independent variable;
- the index, <math>~m</math>, must be set to zero, so there is no summation over this index; and,
- every surviving spherical harmonic can be written more simply in terms of a Legendre function, namely,
|
<math>~Y_{\ell m} \rightarrow Y_{\ell 0}</math> |
<math>~=</math> |
<math>~\sqrt{\frac{(2\ell+1 )}{4\pi}} P_\ell(\chi) \, ,</math> |
where, |
<math>~\chi \equiv \frac{z}{(\varpi^2 + z^2)^{1 / 2}} \, .</math> |
Note that the argument, <math>~\chi</math>, is still the spherical-coordinate expression, <math>~\cos\theta</math>, but here it has been written in terms of cylindrical coordinates. We have, therefore,
|
<math>~ \Phi_B(\varpi, z)\biggr|_{2D}</math> |
<math>~=</math> |
<math>~ - G \sum_{\ell=0}^{\ell_\mathrm{max}} P_\ell(\chi) \biggl[ \varpi^2 + z^2 \biggr]^{-(\ell+1)/2} \int P_\ell(\chi^') \biggl[ (\varpi^')^2 + (z^')^2 \biggr]^{\ell/2} ~\rho(\varpi^', z^') d^3x^' \, , </math> |
where, <math>~d^3x^' = 2\pi \varpi^' d\varpi^' dz^'</math>. This is precisely the same as equation (5) from Black & Bodenheimer (1975).
Using Toroidal Functions
|
Written in terms of toroidal functions we have,
|
<math>~ \Phi_B(\varpi,\theta,z)</math> |
<math>~=</math> |
<math>~ -G \int \frac{1}{|\vec{x}^{~'} - \vec{x}|} ~\rho(\vec{x}^{~'}) d^3x^' </math> |
|
|
<math>~=</math> |
<math>~ -G \int \biggl\{ \frac{1}{\pi \sqrt{\varpi \varpi^'}} \sum_{m=-\infty}^{\infty} e^{im(\phi - \phi^')}Q_{m- 1 / 2}(\chi) \biggr\}~ \rho(\vec{x}^{~'}) d^3x^' </math> |
|
H. S. Cohl & J. E. Tohline (1999), p. 88, Eqs. (15) & (16) |
||
where:
<math>~\chi \equiv \frac{\varpi^2 + (\varpi^')^2 + (z - z^')^2}{2\varpi \varpi^'} \, .</math>
Extension from 2D to 3D
An extension from 2D to 3D was accomplished by using separate, finite Fourier series expansion to represent the azimuthal dependence of <math>~\rho</math> and <math>~\Phi</math>. This allowed us to decouple the differential expression for the Poisson equation into a finite number of independent Fourier modes, and to straightforwardly solve a finite set of 2D Helmholtz equations, instead of solving a single 2D Poisson equation.
Intermediate
|
Much of the material enclosed in this text box has been drawn from §3.1 of Jackson (1975). In spherical coordinates, <math>~(r,\theta,\phi)</math>, the Poisson equation takes the form,
As it turns out, the potential — which, generally, is a function of all three coordinates — can be written as the product of three functions, each of which is a function of only one coordinate. To demonstrate this, specifically suppose that,
Then the differential form of the Poisson equation can be rewritten as,
Setting, <math>~Q = e^{\pm im\phi}</math>, gives,
Now, if for each selected pair of the integer indexes, <math>~(\ell, m)</math>, the function <math>~P</math> is identified as the associated Legendre function, <math>~P_\ell^m</math>, as defined above in Table 2, then the set of terms inside the square brackets on the right-hand-side of this last expression can be set equal to <math>~[-\ell(\ell + 1)]</math> because every associated Legendre function satisfies the differential equation,
Rewriting our adopted separable expression for the potential as the sum over all possible <math>~(\ell,m)</math>-component solutions, that is,
a solution to the Poisson equation then imposes the requirement that, for each <math>~(\ell,m)</math> pair,
|
Other Approaches
Deupree (1974)
|
|
"The fully nonlinear, nonradial, nonadiabatic calculation of stellar oscillations has not as yet been attempted by anyone, to the best of the author's knowledge. Only Deupree (1974b, 1975) so far seems to have taken some steps in this direction. He has carried out numerical calculations of the nonlinear axisymmetric oscillations in the adiabatic (Deupree 1974b) and nonadiabatic (Deupree 1975) approximations, with apparently encouraging results." |
|
— Drawn from §3.5 of the review article by J. P. Cox (1976, ARAA, 14, 247 - 273) |
Apparently, R. G. Deupree (1974, ApJ, 194, 393 - 402) was the first astrophysicist to employ self-gravitating, numerical hydrodynamic techniques to model nonlinear, nonradial stellar pulsations. His earliest simulations were restricted to the examination of axisymmetric (2D) configurations and dynamical motions.
Norman and Wilson (1978)
|
Following the discussion presented in §4.1 of Jackson (1975), it is sometimes useful to rewrite Form A of the boundary potential as,
|
<math>~ \Phi_B(r,\theta,\phi)</math> |
<math>~=</math> |
<math>~ -4\pi G \sum_{\ell=0}^{\infty} \sum_{m=-\ell}^{+\ell} \frac{q_{\ell m}}{r^{\ell+1}} \frac{Y_{\ell m}(\theta,\phi)}{(2\ell+1)} \, , </math> |
where we have introduced what is commonly referred to as the,
|
Multipole Moments of the Mass Distribution |
||
|
<math>~q_{\ell m}</math> |
<math>~\equiv</math> |
<math>~ \int (r^')^\ell Y_{\ell m}^*(\theta^', \phi^') ~\rho(r^', \theta^', \phi^') d^3x^' \, . </math> |
|
Jackson (1975), p. 137, Eq. (4.3) |
||
When written explicitly in terms of cartesian coordinates — see Table 2, below, for each of the relevant <math>~Y_{\ell m}</math> expressions — the first few of these moments have the functional representations derived below. For the cases that correspond to positive values of the index, <math>~m</math>, the set of multipole moment expressions that have been included in our Table 3 summary exactly matches the set of expressions presented as equations (4.4), (4.5), and (4.6) in Jackson (1975). (Jackson's equation 4.7 explains how to map these expressions to the cases corresponding to negative values of the index, <math>~m</math>.)
Let's now look at various terms in the summed expression for the boundary potential with each term expressing the contribution for a separate value of the index, <math>~\ell</math>. Isolating the first three terms from all the rest, for example, we have,
|
<math>~ \Phi_B(r,\theta,\phi)</math> |
<math>~=</math> |
<math>~ \Phi_B(r,\theta,\phi) \biggr|_{\ell=0} +~ \Phi_B(r,\theta,\phi)\biggr|_{\ell=1} + ~\Phi_B(r,\theta,\phi)\biggr|_{\ell=2} -~4\pi G \sum_{\ell=3}^{\infty} \sum_{m=-\ell}^{+\ell} \frac{q_{\ell m}}{r^{\ell+1}} \frac{Y_{\ell m}(\theta,\phi)}{(2\ell+1)} \, . </math> |
We recognize, first, that as we consider boundary points that lie farther and farther away from the mass distribution, the magnitude of the first <math>(\ell = 0)</math> term drops off as <math>~r^{-1}</math> — the expected behavior of the potential outside of a point mass; the second <math>~(\ell=1)</math> term drops off as <math>~r^{-3}</math>; and the third <math>~(\ell=2)</math> term drops off as <math>~r^{-5}</math>. Below, we have evaluated in more detail the behavior of these first three terms. This evaluation gives us the,
|
Boundary Potential Written in Terms of Multipole Moments |
||
|
<math>~ \Phi_B(r,\theta,\phi)</math> |
<math>~=</math> |
<math>~ ~-~~\frac{GM}{r} ~~- ~~\frac{GM}{r^3} \biggl[ \vec{x} \cdot \vec{x}_\mathrm{com} \biggr] ~~-~~\frac{G}{2r^5} \sum_{i=1}^3 \sum_{j=1}^3 Q_{i,j} \biggl[ x_i x_j \biggr] ~~-~~4\pi G \sum_{\ell=3}^{\infty} \sum_{m=-\ell}^{+\ell} \frac{q_{\ell m}}{r^{\ell+1}} \frac{Y_{\ell m}(\theta,\phi)}{(2\ell+1)} \, , </math> |
where, as defined below, <math>~\vec{x}_\mathrm{com}</math> is the center-of-mass location, and <math>~Q_{i,j} </math> is the traceless quadrupole moment tensor.
If a modeled mass-density distribution, <math>~\rho( \vec{x}^{~'})</math>, has been configured such that the center-of-mass of the system coincides with the origin of the coordinate system — that is, if it has been configured such that <math>~\vec{x}_\mathrm{com} = 0</math> — then the second term in this series can be set to zero. If, in addition, all terms having <math>~\ell \ge 3</math> are ignored because their values drop off rapidly with distance — specifically, inverse distance to the <math>~(2\ell + 1)</math> power — then a reasonably good approximation for the potential on the boundary of the modeled system is given by the expression,
|
<math>~ \Phi_B</math> |
<math>~\approx</math> |
<math>~ -\frac{GM}{r} ~-\frac{G}{2r^5} \sum_{i=1}^3 \sum_{j=1}^3 Q_{i,j} \biggl[ x_i x_j \biggr] \, . </math> |
This is precisely the relation that was adopted by M. L. Norman & J. R. Wilson (1978, ApJ, 224, 497 - 511) when, in the context of star formation, they modeled the Fragmentation of Isothermal Rings — see specifically their equation (18).
Boss (1980)
|
Building upon equations (5) and (7) from A. P. Boss (1980, ApJ, 236, 619 - 627) — see also Jackson's equation (3.58) — we can write,
|
<math>~\Phi(r,\theta,\phi)</math> |
<math>~=</math> |
<math>~ \sum_{\ell=0}^\infty \sum_{m = -\ell}^{\ell} \Phi_{\ell m} ( r ) Y_{\ell m}(\theta,\phi) \, , </math> |
and,
|
<math>~\rho(r,\theta,\phi)</math> |
<math>~=</math> |
<math>~ \sum_{\ell=0}^\infty \sum_{m = -\ell}^{\ell} \rho_{\ell m} ( r ) Y_{\ell m}(\theta,\phi) \, , </math> |
where the coefficients in the series expansion of <math>~\rho</math> can each be obtained from the known spatial density distribution via the integral expression,
|
<math>~\rho_{\ell m}(r)</math> |
<math>~=</math> |
<math>~ \int d\Omega Y_{\ell m}^*(\theta,\phi) \rho(r,\theta,\phi) \, , </math> |
and the differential solid angle is,
<math>~d\Omega \equiv \sin\theta d\theta d\phi \, .</math>
Stahler (1983)
|
In the chapter of this H_Book that focuses on a discussion of Dyson-Wong tori, we have included the expression for the gravitational potential of a thin ring of mass, <math>~M</math>, that passes through the meridional plane at coordinate location, <math>~(\varpi^', z^') = (a, 0)</math>, as derived, for example, by O. D. Kellogg (1929) and by W. D. MacMillan (1958; originally, 1930), namely,
|
<math>~\Phi(\varpi, z)</math> |
<math>~=</math> |
<math>~ - \biggl[ \frac{2GMc}{\pi\rho_1}\biggr] K(k) </math> |
|
|
<math>~=</math> |
<math>~ - \biggl[ \frac{2GM}{\pi } \biggr]\frac{1}{\sqrt{(\varpi+a)^2 + z^2}} \times K(k) \, , </math> |
where,
|
<math>~k</math> |
<math>~=</math> |
<math>~ \biggl[ \frac{4a\varpi }{(\varpi+a)^2 + z^2} \biggr]^{1 / 2} \, . </math> |
Stahler (1983a) has argued that a reasonably good approximation to the gravitational potential due to any extended axisymmetric mass distribution can be obtained by adding up the contributions due to many thin rings — <math>~\delta M(\varpi^', z^')</math> being the appropriate differential mass contributed by each ring element — positioned at various meridional coordinate locations throughout the mass distribution. According to his independent derivation, the differential contribution to the potential, <math>~\delta\Phi_g(\varpi, z)</math>, due to each differential mass element is (see his equation 11 and the explanatory text that follows it):
|
<math>~\delta\Phi_g(\varpi,z)</math> |
<math>~=</math> |
<math>~ - \biggl[\frac{2G}{\pi \varpi^'}\biggr] \frac{\delta M}{[(\alpha + 1)^2 + \beta^2]^{1 / 2}} \times K\biggl\{ \biggl[ \frac{4\alpha}{(\alpha+1)^2 + \beta^2} \biggr]^{1 / 2} \biggr\} </math> |
|
|
<math>~=</math> |
<math>~ - \biggl[\frac{2G}{\pi }\biggr] \frac{\delta M}{[(\varpi + \varpi^')^2 + (z^' - z)^2]^{1 / 2}} \times K\biggl\{ \biggl[ \frac{4\varpi^' \varpi}{(\varpi +\varpi^')^2 + (z^' - z)^2} \biggr]^{1 / 2} \biggr\} \, . </math> |
Stahler's expression for each thin ring contribution is clearly the same as the expressions presented by Kellogg (1929) and MacMillan (1958) for a ring that cuts through the meridional plane at <math>~(\varpi^', z^') = (a, 0)</math>.
From our broad analysis of the integral Poisson equation expressed in cylindrical coordinates, we can independently state — see the discussion of our original derivation, or separate summary — that the exact integral expression for the gravitational potential due to any axisymmetric mass-density distribution, <math>~\rho(\varpi^', z^')</math>, is,
|
<math>~\Phi(\varpi, z)</math> |
<math>~=</math> |
<math>~- \frac{2G}{\varpi^{1 / 2}} \int\int (\varpi^')^{1 / 2} \mu K(\mu) \rho(\varpi^', z^') d\varpi^' dz^' \, ,</math> |
where,
|
<math>~\mu^2</math> |
<math>~=</math> |
<math>~ \biggl[\frac{4\varpi^' \varpi}{(\varpi + \varpi^')^2 + (z^' - z)^2} \biggr] \, . </math> |
Recognizing that, for axisymmetric structures, the differential mass element is, <math>~dM^' = 2\pi \rho(\varpi^', z^') \varpi^' d\varpi^' dz^'</math>, this integral expression may be rewritten as,
|
<math>~\Phi(\varpi, z)</math> |
<math>~=</math> |
<math>~- \frac{2G}{\varpi^{1 / 2}} \int\int (\varpi^')^{1 / 2} \mu K(\mu) \biggl[ \frac{dM^'}{2\pi \varpi^'} \biggr] </math> |
|
|
<math>~=</math> |
<math>~- \frac{G}{\pi} \int\int \biggl[ \frac{1}{\varpi^'\varpi}\biggr]^{1 / 2} \biggl[\frac{4\varpi^' \varpi}{(\varpi + \varpi^')^2 + (z^' - z)^2} \biggr]^{1 / 2} K(\mu) dM^' </math> |
|
|
<math>~=</math> |
<math>~- \frac{2G}{\pi} \int\int \biggl[\frac{1}{(\varpi + \varpi^')^2 + (z^' - z)^2} \biggr]^{1 / 2} K(\mu) dM^' \, .</math> |
We see, therefore, that our differential contribution to the potential exactly matches Stahler's (1983a).
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
|
- 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.
- 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.
- 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>.
- 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.
- 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).
- 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.
- 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.
- 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.
- 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.
- 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.
Special Functions & Other Broadly Used Representations
Spherical Harmonics and Associated Legendre Functions
| Table 2: Green's Function in Terms of Associated Legendre Functions, <math>~P_\ell^m(\cos\theta)</math>, and the Spherical Harmonics, <math>~Y_{\ell m}(\theta,\phi)</math> |
|||||||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
J. D. Jackson (1975, 2nd Edition), Eq. (3.70) |
|||||||||||||||||||||||||||||||
|
Note: <math>~Y_{\ell m}(\theta,\phi) = \biggl[ \frac{(2\ell + 1 )(\ell - m)!}{4\pi(\ell + m)!} \biggr]^{1 / 2} P_\ell^m(\cos\theta)e^{im\phi}</math>
Jackson (1975), Eq. (3.53)
<math>~P_{\ell}^m(x) = (-1)^m (1-x^2)^{m/2} ~ \frac{d^m}{dx^m} P_\ell(x)</math>
Jackson (1975), Eq. (3.49)
<math>~P_{\ell}(x) = \frac{1}{2^\ell \ell !} ~ \frac{d^\ell}{dx^\ell} (x^2 - 1)^\ell</math>
Jackson (1975), Eq. (3.16) |
Leading Legendre Functions
Jackson (1975), Eq. (3.15) |
||||||||||||||||||||||||||||||
|
|||||||||||||||||||||||||||||||
| Leading Spherical Harmonics
Jackson (1975), §3.5 |
|||||||||||||||||||||||||||||||
|
<math>~Y_{\ell 0}(\theta,\phi) = \sqrt{\frac{(2\ell + 1 )}{4\pi} } ~P_\ell(\cos\theta)</math>
Jackson (1975), Eq. (3.57) |
||||||||||||||||||||||||||||||
Multipole Expansions
Mass Multipole Moments
As an extension of the above discussion of
|
Multipole Moments of the Mass Distribution |
||
|
<math>~q_{\ell m}</math> |
<math>~\equiv</math> |
<math>~ \int (r^')^\ell Y_{\ell m}^*(\theta^', \phi^') ~\rho(r^', \theta^', \phi^') d^3x^' \, , </math> |
|
Jackson (1975), p. 137, Eq. (4.3) |
||
here we evaluate a set of the leading order mass-multipole moments in terms of their cartesian-coordinate representations.
|
<math>~q_{00}</math> |
<math>~=</math> |
<math>~\frac{1}{\sqrt{4\pi}} \int \rho(\vec{x}^{~'}) d^3x^' \, ; </math> |
|
<math>~q_{11}</math> |
<math>~=</math> |
<math>~ \int r^' \biggl[ Y_{11}^* \biggr]~\rho(\vec{x}^{~'}) d^3x^' </math> |
|
|
<math>~=</math> |
<math>~ - ~\int r^' \biggl[ \sqrt{\frac{3}{8\pi}}\sin\theta ~\biggl(\cos\phi^' - i\sin\phi^' \biggr) \biggr]~\rho(\vec{x}^{~'}) d^3x^' </math> |
|
|
<math>~=</math> |
<math>~ - \sqrt{\frac{3}{8\pi}}~\int r^' \biggl[ \frac{\sqrt{(x^')^2+(y^')^2}}{r^'} ~\biggl(\frac{x^'}{\sqrt{(x^')^2+(y^')^2}} - \frac{iy^'}{\sqrt{(x^')^2+(y^')^2}} \biggr)\biggr] ~\rho(\vec{x}^{~'}) d^3x^' </math> |
|
|
<math>~=</math> |
<math>~ - \sqrt{\frac{3}{8\pi}}~\int (x^' - iy^') ~\rho(\vec{x}^{~'}) d^3x^' \, ; </math> |
|
<math>~q_{1,- 1}</math> |
<math>~=</math> |
<math>~ + \sqrt{\frac{3}{8\pi}}~\int (x^' + iy^') ~\rho(\vec{x}^{~'}) d^3x^' \, ; </math> |
|
<math>~q_{10}</math> |
<math>~=</math> |
<math>~ \int r^' \biggl[ Y_{10}^* \biggr]~\rho(\vec{x}^{~'}) d^3x^' </math> |
|
|
<math>~=</math> |
<math>~ \int r^' \biggl[ \sqrt{\frac{3}{4\pi}}~\cos\theta^' \biggr]~\rho(\vec{x}^{~'}) d^3x^' </math> |
|
|
<math>~=</math> |
<math>~ \sqrt{\frac{3}{4\pi}} \int z^' \rho(\vec{x}^{~'}) d^3x^' \, ; </math> |
|
<math>~q_{22}</math> |
<math>~=</math> |
<math>~ \int (r^')^2 \biggl[ Y_{22}^*(\theta^', \phi^') \biggr]~\rho(\vec{x}^{~'}) d^3x^' </math> |
|
|
<math>~=</math> |
<math>~ \int (r^')^2 \biggl\{ \frac{1}{4}\sqrt{\frac{15}{2\pi}} ~\sin^2\theta^' ~\biggl[\cos(2\phi^') - i\sin(2\phi^')\biggr] \biggr\}~\rho(\vec{x}^{~'}) d^3x^' </math> |
|
|
<math>~=</math> |
<math>~\frac{1}{4}\sqrt{\frac{15}{2\pi}} \int (r^')^2 \biggl\{ \frac{(x^')^2+(y^')^2}{(r^')^2} ~\biggl[\frac{ (x^')^2 - (y^')^2 - 2i ( x^' y^' )}{(x^')^2+(y^')^2} \biggr] \biggr\}~\rho(\vec{x}^{~'}) d^3x^' </math> |
|
|
<math>~=</math> |
<math>~\frac{1}{4}\sqrt{\frac{15}{2\pi}} \int (x^' -iy^')^2~\rho(\vec{x}^{~'}) d^3x^' \, ; </math> |
|
<math>~q_{2,-2}</math> |
<math>~=</math> |
<math>~\frac{1}{4}\sqrt{\frac{15}{2\pi}} \int (x^' + iy^')^2~\rho(\vec{x}^{~'}) d^3x^' \, ; </math> |
|
<math>~q_{21}</math> |
<math>~=</math> |
<math>~ \int (r^')^2 \biggl[ Y_{21}^*(\theta^', \phi^') \biggr]~\rho(\vec{x}^{~'}) d^3x^' </math> |
|
|
<math>~=</math> |
<math>~ - \sqrt{\frac{15}{8\pi}} \int (r^')^2 \biggl[ \sin\theta^' \cos\theta^' (\cos\phi^' -i\sin\phi^') \biggr]~\rho(\vec{x}^{~'}) d^3x^' </math> |
|
|
<math>~=</math> |
<math>~ - \sqrt{\frac{15}{8\pi}} \int (r^')^2 \biggl[ \frac{z^' \sqrt{(x^')^2+(y^')^2}}{(r^')^2} ~\biggl(\frac{x^'}{\sqrt{(x^')^2+(y^')^2}} - \frac{iy^'}{\sqrt{(x^')^2+(y^')^2}} \biggr)\biggr] ~\rho(\vec{x}^{~'}) d^3x^' </math> |
|
|
<math>~=</math> |
<math>~ - \sqrt{\frac{15}{8\pi}} \int z^'(x^' - iy^') ~\rho(\vec{x}^{~'}) d^3x^' \, ; </math> |
|
<math>~q_{2,-1}</math> |
<math>~=</math> |
<math>~ + \sqrt{\frac{15}{8\pi}} \int z^'(x^' + iy^') ~\rho(\vec{x}^{~'}) d^3x^' \, ; </math> |
|
<math>~q_{20}</math> |
<math>~=</math> |
<math>~ \int (r^')^2 \biggl[ Y_{20}^*(\theta^', \phi^') \biggr]~\rho(\vec{x}^{~'}) d^3x^' </math> |
|
|
<math>~=</math> |
<math>~ \int (r^')^2 \biggl[ \sqrt{\frac{5}{16\pi}} (3\cos^2\theta^' - 1) \biggr]~\rho(\vec{x}^{~'}) d^3x^' </math> |
|
|
<math>~=</math> |
<math>~\sqrt{\frac{5}{16\pi}} \int \biggl[ 3(z^')^2 - (r^')^2 \biggr]~\rho(\vec{x}^{~'}) d^3x^' \, . </math> |
| Table 3: Summary of Cartesian Expressions for Leading Mass-Multipole Moments | ||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
|
Contributions to the Gravitational Potential in Terms of Multipole Moments
Expanding on the above discussion of the gravitational potential, written as a series expansion in terms of mass-multipole moments and spherical harmonics,
|
<math>~ \Phi_B(r,\theta,\phi)</math> |
<math>~=</math> |
<math>~ \Phi_B(r,\theta,\phi) \biggr|_{\ell=0} +~ \Phi_B(r,\theta,\phi)\biggr|_{\ell=1} + ~\Phi_B(r,\theta,\phi)\biggr|_{\ell=2} -~4\pi G \sum_{\ell=3}^{\infty} \sum_{m=-\ell}^{+\ell} \frac{q_{\ell m}}{r^{\ell+1}} \frac{Y_{\ell m}(\theta,\phi)}{(2\ell+1)} \, , </math> |
here we evaluate the first three, leading order terms.
First Term
First (term with <math>~\ell = 0</math>):
|
<math>~ \Phi_B(r,\theta,\phi)\biggr|_{\ell = 0}</math> |
<math>~=</math> |
<math>~ -4\pi G \biggl[ q_{00} \biggr] \frac{ Y_{00}(\theta,\phi)}{r} </math> |
|
|
<math>~=</math> |
<math>~ -4\pi G \biggl[ \frac{1}{\sqrt{4\pi}} \int \rho(\vec{x}^{~'}) d^3x^' \biggr] \frac{1}{\sqrt{4\pi}~r} </math> |
|
|
<math>~=</math> |
<math>~ -\frac{GM}{r} \, . </math> |
Second Term
Second (terms with <math>~\ell = 1</math>):
|
<math>~ \Phi_B(r,\theta,\phi)\biggr|_{\ell = 1}</math> |
<math>~=</math> |
<math>~ -\frac{4\pi G}{3} \sum_{m=-1}^{+1} \biggl[ q_{1 m}\biggr] \frac{Y_{1 m}(\theta,\phi)}{r^2} </math> |
|
|
<math>~=</math> |
<math>~ -\frac{4\pi G}{3r^2} \biggl\{ \biggl[ q_{1,- 1}\biggr] Y_{1, -1}(\theta,\phi) + \biggl[ q_{1 0}\biggr] Y_{1 0}(\theta,\phi) + \biggl[ q_{1 1}\biggr] Y_{1 1}(\theta,\phi) \biggr\} </math> |
|
|
<math>~=</math> |
<math>~ -\frac{4\pi G}{3r^2} \biggl\{ \biggl[ + \frac{3}{8\pi}~\int (x^' + iy^')~\rho(\vec{x}^{~'}) d^3x^' \biggr] \frac{(x - iy) }{r} </math> |
|
|
|
<math>~ + \biggl[ \frac{3}{4\pi} \int z^' \rho(\vec{x}^{~'}) d^3x^' \biggr] \frac{z}{r} + \biggl[ \frac{3}{8\pi}~\int (x^' - iy^')~\rho(\vec{x}^{~'}) d^3x^' \biggr] \frac{(x + iy) }{r} \biggr\} </math> |
|
|
<math>~=</math> |
<math>~ -\frac{G}{r^3} \biggl\{ \biggl[ \frac{1}{2}~\int (x^' + iy^')~\rho(\vec{x}^{~'}) d^3x^' \biggr] (x - iy) </math> |
|
|
|
<math>~ + \biggl[ \int z^' \rho(\vec{x}^{~'}) d^3x^' \biggr] z + \biggl[ \frac{1}{2}~\int (x^' - iy^')~\rho(\vec{x}^{~'}) d^3x^' \biggr] (x + iy) \biggr\} </math> |
|
|
<math>~=</math> |
<math>~ -\frac{G}{r^3} \biggl\{ \biggl[ \int (x^')~\rho(\vec{x}^{~'}) d^3x^' \biggr] (x ) + \biggl[ \int (y^')~\rho(\vec{x}^{~'}) d^3x^' \biggr] (y ) + \biggl[ \int z^' \rho(\vec{x}^{~'}) d^3x^' \biggr] z \biggr\} </math> |
|
<math>~</math> |
<math>~=</math> |
<math>~- \frac{GM}{r^3} \biggl[ \vec{x} \cdot \vec{x}_\mathrm{com} \biggr]</math> |
|
|
<math>~ \vec{x}_\mathrm{com} \equiv \frac{1}{M} \int ~\vec{x}^{~'} \rho(\vec{x}^{~'}) d^3x^'\, . </math> |
Third Term
Third (terms with <math>~\ell = 2</math>):
|
<math>~ \Phi_B(r,\theta,\phi)\biggr|_{\ell = 2}</math> |
<math>~=</math> |
<math>~ -\frac{4\pi G}{5r^3} \biggl\{ \biggl[ q_{2, -2} \biggr] Y_{2,-2} + \biggl[ q_{2, 2} \biggr] Y_{2,2} + \biggl[ q_{2,-1} \biggr] Y_{2, -1} + \biggl[ q_{2,1} \biggr] Y_{2, 1} + \biggl[ q_{20} \biggr] Y_{20} \biggr\} </math> |
|
|
<math>~=</math> |
<math>~ -\frac{4\pi G}{5r^3} \biggl\{ \biggl[ \frac{1}{4}\sqrt{\frac{15}{2\pi}} \int (x^' + iy^')^2~\rho(\vec{x}^{~'}) d^3x^' \biggr] \frac{1}{4}\sqrt{\frac{15}{2\pi}} \sin^2\theta e^{-2i\phi} + \biggl[\frac{1}{4}\sqrt{\frac{15}{2\pi}} \int (x^' -iy^')^2~\rho(\vec{x}^{~'}) d^3x^' \biggr] \frac{1}{4}\sqrt{\frac{15}{2\pi}} \sin^2\theta e^{2i\phi} </math> |
|
|
|
<math>~ + \biggl[ \sqrt{\frac{15}{8\pi}} \int z^'(x^' + iy^') ~\rho(\vec{x}^{~'}) d^3x^' \biggr] \sqrt{\frac{15}{8\pi}} \sin\theta \cos\theta e^{-i\phi} + \biggl[ \sqrt{\frac{15}{8\pi}} \int z^'(x^' - iy^') ~\rho(\vec{x}^{~'}) d^3x^' \biggr] \sqrt{\frac{15}{8\pi}} \sin\theta \cos\theta e^{i\phi} </math> |
|
|
|
<math>~ + \biggl[ \sqrt{\frac{5}{16\pi}} \int \biggl[ 3(z^')^2 - (r^')^2 \biggr]~\rho(\vec{x}^{~'}) d^3x^' \biggr] \sqrt{\frac{5}{16\pi}} (3\cos^2\theta - 1) \biggr\} </math> |
|
|
<math>~=</math> |
<math>~ -\frac{4\pi G}{5r^3} \biggl\{ \frac{3\cdot 5}{2^5 \pi} \biggl(\frac{\varpi}{r}\biggr)^2 \biggl[ \int (x^' + iy^')^2~\rho(\vec{x}^{~'}) d^3x^' \biggr] \biggl[ \cos(2\phi) - i\sin(2\phi) \biggr] + \frac{3\cdot 5}{2^5 \pi} \biggl(\frac{\varpi}{r}\biggr)^2 \biggl[\int (x^' -iy^')^2~\rho(\vec{x}^{~'}) d^3x^' \biggr] \biggl[ \cos(2\phi) + i\sin(2\phi) \biggr] </math> |
|
|
|
<math>~ + \frac{3\cdot 5}{2^3\pi} \biggl( \frac{\varpi z}{r^2} \biggr) \biggl[\int z^'(x^' + iy^') ~\rho(\vec{x}^{~'}) d^3x^' \biggr] \biggl[ \cos\phi - i\sin\phi \biggr] + \frac{3\cdot 5}{2^3\pi} \biggl( \frac{\varpi z}{r^2} \biggr) \biggl[ \int z^'(x^' - iy^') ~\rho(\vec{x}^{~'}) d^3x^' \biggr] \biggl[ \cos\phi + i\sin\phi \biggr] </math> |
|
|
|
<math>~ + \frac{5}{2^4\pi} \biggl[ \int \biggl[ 3(z^')^2 - (r^')^2 \biggr]~\rho(\vec{x}^{~'}) d^3x^' \biggr] \biggl( \frac{3z^2 - r^2}{r^2} \biggr) \biggr\} </math> |
|
|
<math>~=</math> |
<math>~ -\frac{4G}{r^5} \biggl\{ \frac{3}{2^5 } \biggl[ \int (x^' + iy^')^2~\rho(\vec{x}^{~'}) d^3x^' \biggr] \biggl[ x^2 - y^2 - 2i xy\biggr] + \frac{3}{2^5 } \biggl[\int (x^' -iy^')^2~\rho(\vec{x}^{~'}) d^3x^' \biggr] \biggl[ x^2 - y^2 + 2i xy \biggr] </math> |
|
|
|
<math>~ + \frac{3z}{2^3} \biggl[\int z^'(x^' + iy^') ~\rho(\vec{x}^{~'}) d^3x^' \biggr] \biggl[x - i y \biggr] + \frac{3z}{2^3} \biggl[ \int z^'(x^' - iy^') ~\rho(\vec{x}^{~'}) d^3x^' \biggr] \biggl[x + i y \biggr] </math> |
|
|
|
<math>~ + \biggl( \frac{3z^2 - r^2}{2^4} \biggr) \biggl[ \int \biggl[ 3(z^')^2 - (r^')^2 \biggr]~\rho(\vec{x}^{~'}) d^3x^' \biggr] \biggr\} </math> |
|
|
<math>~=</math> |
<math>~ -\frac{G}{2r^5} \biggl\{ \frac{3( x^2 - y^2 )}{2^2 } \biggl[ \int (x^' + iy^')^2~\rho(\vec{x}^{~'}) d^3x^' + \int (x^' -iy^')^2~\rho(\vec{x}^{~'}) d^3x^' \biggr] + \frac{3i xy}{2 } \biggl[\int (x^' -iy^')^2~\rho(\vec{x}^{~'}) d^3x^' - \int (x^' + iy^')^2~\rho(\vec{x}^{~'}) d^3x^' \biggr] </math> |
|
|
|
<math>~ + 3xz \biggl[\int z^'(x^' + iy^') ~\rho(\vec{x}^{~'}) d^3x^' + \int z^'(x^' - iy^') ~\rho(\vec{x}^{~'}) d^3x^'\biggr] + 3i yz\biggl[ \int z^'(x^' - iy^') ~\rho(\vec{x}^{~'}) d^3x^' - \int z^'(x^' + iy^') ~\rho(\vec{x}^{~'}) d^3x^' \biggr] </math> |
|
|
|
<math>~ + \biggl( \frac{3z^2 - r^2}{2} \biggr) \biggl[ \int \biggl[ 3(z^')^2 - (r^')^2 \biggr]~\rho(\vec{x}^{~'}) d^3x^' \biggr] \biggr\} </math> |
|
|
<math>~=</math> |
<math>~ -\frac{G}{2r^5} \biggl\{ \frac{3( x^2 - y^2 )}{2^2 } \biggl[ \int [2(x^')^2 - 2(y^')^2]~\rho(\vec{x}^{~'}) d^3x^' \biggr] + \frac{3i xy}{2 } \biggl[\int -4i x^' y^'~\rho(\vec{x}^{~'}) d^3x^' \biggr] </math> |
|
|
|
<math>~ + 3xz \biggl[\int 2x^' z^' ~\rho(\vec{x}^{~'}) d^3x^' \biggr] + 3i yz\biggl[ \int -2iy^' z^' ~\rho(\vec{x}^{~'}) d^3x^' \biggr] + \biggl( \frac{3z^2 - r^2}{2} \biggr) \biggl[ \int \biggl[ 3(z^')^2 - (r^')^2 \biggr]~\rho(\vec{x}^{~'}) d^3x^' \biggr] \biggr\} </math> |
|
|
<math>~=</math> |
<math>~ -\frac{G}{2r^5} \biggl\{ \frac{3( x^2 - y^2 )}{2} \biggl[ \int [(x^')^2 - (y^')^2]~\rho(\vec{x}^{~'}) d^3x^' \biggr] + 6 xy \biggl[\int x^' y^'~\rho(\vec{x}^{~'}) d^3x^' \biggr] </math> |
|
|
|
<math>~ + 6xz \biggl[\int x^' z^' ~\rho(\vec{x}^{~'}) d^3x^' \biggr] + 6yz\biggl[ \int y^' z^' ~\rho(\vec{x}^{~'}) d^3x^' \biggr] + \biggl( \frac{3z^2 - r^2}{2} \biggr) \biggl[ \int \biggl[ 3(z^')^2 - (r^')^2 \biggr]~\rho(\vec{x}^{~'}) d^3x^' \biggr] \biggr\} \, . </math> |
Rearranging terms, we have,
|
<math>~ \Phi_B(r,\theta,\phi)\biggr|_{\ell = 2}</math> |
<math>~=</math> |
<math>~ -\frac{G}{2r^5} \biggl\{ 2 xy \biggl[\int 3x^' y^'~\rho(\vec{x}^{~'}) d^3x^' \biggr] + 2xz \biggl[\int 3x^' z^' ~\rho(\vec{x}^{~'}) d^3x^' \biggr] + 2yz\biggl[ \int 3y^' z^' ~\rho(\vec{x}^{~'}) d^3x^' \biggr] </math> |
|
|
|
<math>~ +\frac{3( x^2 - y^2 )}{2} \biggl[ \int [(x^')^2 - (y^')^2]~\rho(\vec{x}^{~'}) d^3x^' \biggr] + \biggl( \frac{3z^2 - r^2}{2} \biggr) \biggl[ \int \biggl[ 3(z^')^2 - (r^')^2 \biggr]~\rho(\vec{x}^{~'}) d^3x^' \biggr] \biggr\} </math> |
|
|
<math>~=</math> |
<math>~ -\frac{G}{2r^5} \biggl\{ 2 xy \biggl[\int 3x^' y^'~\rho(\vec{x}^{~'}) d^3x^' \biggr] + 2xz \biggl[\int 3x^' z^' ~\rho(\vec{x}^{~'}) d^3x^' \biggr] + 2yz\biggl[ \int 3y^' z^' ~\rho(\vec{x}^{~'}) d^3x^' \biggr] </math> |
|
|
|
<math>~ +\frac{3( x^2 )}{2} \biggl[ \int [(x^')^2 - (y^')^2]~\rho(\vec{x}^{~'}) d^3x^' \biggr] - \biggl( \frac{x^2}{2} \biggr) \biggl[ \int \biggl[ 3(z^')^2 - (r^')^2 \biggr]~\rho(\vec{x}^{~'}) d^3x^' \biggr] </math> |
|
|
|
<math>~ -\frac{3( y^2 )}{2} \biggl[ \int [(x^')^2 - (y^')^2]~\rho(\vec{x}^{~'}) d^3x^' \biggr] - \biggl( \frac{y^2}{2} \biggr) \biggl[ \int \biggl[ 3(z^')^2 - (r^')^2 \biggr]~\rho(\vec{x}^{~'}) d^3x^' \biggr] </math> |
|
|
|
<math>~ + z^2 \biggl[ \int \biggl[ 3(z^')^2 - (r^')^2 \biggr]~\rho(\vec{x}^{~'}) d^3x^' \biggr] \biggr\} </math> |
|
|
<math>~=</math> |
<math>~ -\frac{G}{2r^5} \biggl\{ 2 xy \biggl[\int 3x^' y^'~\rho(\vec{x}^{~'}) d^3x^' \biggr] + 2xz \biggl[\int 3x^' z^' ~\rho(\vec{x}^{~'}) d^3x^' \biggr] + 2yz\biggl[ \int 3y^' z^' ~\rho(\vec{x}^{~'}) d^3x^' \biggr] </math> |
|
|
|
<math>~ + z^2 \biggl[ \int \biggl[ 3(z^')^2 - (r^')^2 \biggr]~\rho(\vec{x}^{~'}) d^3x^' \biggr] </math> |
|
|
|
<math>~ +\frac{( x^2 )}{2} \biggl[ \int [3(x^')^2 - 3(y^')^2-3(z^')^2 + (r^')^2 ]~\rho(\vec{x}^{~'}) d^3x^' \biggr] </math> |
|
|
|
<math>~ +\frac{( y^2 )}{2} \biggl[ \int [3(x^')^2 - 3(y^')^2 -3(z^')^2 + (r^')^2 ]~\rho(\vec{x}^{~'}) d^3x^' \biggr] \biggr\} </math> |
|
|
<math>~=</math> |
<math>~ -\frac{G}{2r^5} \biggl\{ 2 xy \biggl[\int 3x^' y^'~\rho(\vec{x}^{~'}) d^3x^' \biggr] + 2xz \biggl[\int 3x^' z^' ~\rho(\vec{x}^{~'}) d^3x^' \biggr] + 2yz\biggl[ \int 3y^' z^' ~\rho(\vec{x}^{~'}) d^3x^' \biggr] </math> |
|
|
|
<math>~ + z^2 \biggl[ \int \biggl[ 3(z^')^2 - (r^')^2 \biggr]~\rho(\vec{x}^{~'}) d^3x^' \biggr] </math> |
|
|
|
<math>~ +\frac{( x^2 )}{2} \biggl[ \int [3(x^')^2 - 3(y^')^2-3(z^')^2 + 3(r^')^2 - 3(r^')^2 + (r^')^2 ]~\rho(\vec{x}^{~'}) d^3x^' \biggr] </math> |
|
|
|
<math>~ +\frac{( y^2 )}{2} \biggl[ \int [-3(x^')^2 + 3(y^')^2 -3(z^')^2 + 3(r^')^2 - 3(r^')^2 + (r^')^2 ]~\rho(\vec{x}^{~'}) d^3x^' \biggr] \biggr\} </math> |
|
|
<math>~=</math> |
<math>~ -\frac{G}{2r^5} \biggl\{ 2 xy \biggl[\int 3x^' y^'~\rho(\vec{x}^{~'}) d^3x^' \biggr] + 2xz \biggl[\int 3x^' z^' ~\rho(\vec{x}^{~'}) d^3x^' \biggr] + 2yz\biggl[ \int 3y^' z^' ~\rho(\vec{x}^{~'}) d^3x^' \biggr] </math> |
|
|
|
<math>~ + z^2 \biggl[ \int \biggl[ 3(z^')^2 - (r^')^2 \biggr]~\rho(\vec{x}^{~'}) d^3x^' \biggr] +\frac{( x^2 )}{2} \biggl[ \int [6(x^')^2 - 2(r^')^2 ]~\rho(\vec{x}^{~'}) d^3x^' \biggr] </math> |
|
|
|
<math>~ +\frac{( y^2 )}{2} \biggl[ \int [6(y^')^2 - 2(r^')^2 ]~\rho(\vec{x}^{~'}) d^3x^' \biggr] \biggr\} </math> |
|
|
<math>~=</math> |
<math>~ -\frac{G}{2r^5} \biggl\{ 2 xy \biggl[\int 3x^' y^'~\rho(\vec{x}^{~'}) d^3x^' \biggr] + 2xz \biggl[\int 3x^' z^' ~\rho(\vec{x}^{~'}) d^3x^' \biggr] + 2yz\biggl[ \int 3y^' z^' ~\rho(\vec{x}^{~'}) d^3x^' \biggr] </math> |
|
|
|
<math>~ + x^2 \biggl[ \int [3(x^')^2 - (r^')^2 ]~\rho(\vec{x}^{~'}) d^3x^' \biggr] + y^2 \biggl[ \int [3(y^')^2 - 2(r^')^2 ]~\rho(\vec{x}^{~'}) d^3x^' \biggr] + z^2 \biggl[ \int \biggl[ 3(z^')^2 - (r^')^2 \biggr]~\rho(\vec{x}^{~'}) d^3x^' \biggr] \biggr\} </math> |
|
|
<math>~=</math> |
<math>~ -\frac{G}{2r^5} \sum_{i=1}^3 \sum_{j=1}^3 Q_{i,j} \biggl[ x_i x_j \biggr] \, , </math> |
where, following Jackson (1975), we have introduced the,
|
Traceless Quadrupole Moment Tensor |
||
|
<math>~Q_{i,j}</math> |
<math>~\equiv</math> |
<math>~ \int [ 3(x_i^') (x_j^') - (r^')^2 \delta_{ij} ] \rho(\vec{x}^{~'}) d^3x^' \, . </math> |
|
Jackson (1975), p. 138, Eq. (4.9) |
||
Toroidal Functions
| Table 4: Green's Function in Terms of Half-Integer Degree Legendre Functions of the 2nd Kind, <math>~Q_{m-1 / 2}(\chi)</math> (also referred to as Toroidal Functions) |
|||
|---|---|---|---|
where: <math>~\chi \equiv \frac{\varpi^2 + (\varpi^')^2 + (z - z^')^2}{2\varpi \varpi^'}</math>
H. S. Cohl & J. E. Tohline (1999), p. 88, Eqs. (15) & (16) |
|||
Related Discussions
Reviews
- P. A. Strittmatter (1969, Annual Review of Astronomy and Astrophysics, 7, 665 - 684) — Stellar Rotation
- N. R. Lebovitz (1967, Annual Review of Astronomy and Astrophysics, 5, 465 - 480) — Rotating Fluid Masses
Solution Methods
- Y. Eriguchi & E. Mueller (1985, A&A, 146, 260 - 268) — A General Computational Method for Obtaining Equilibria of Self-Gravitating and Rotating Gases
- S. W. Stahler (1983, ApJ, 268, 155 - 184) — The Equilibria of Rotating, Isothermal Clouds. I. - Method of Solution
- Y. Eriguchi (1978, PASJ, 30, 507 - 518) — Hydrostatic Equilibria of Rotating Polytropes
- S. I. Blinnikov (1975, Soviet Astronomy, 19, 151 - 156) — Self-Consistent Field Method in the Theory of Rotating Stars
- M. J. Clement (1974, ApJ, 194, 709 - 714) — On the Solution of Poisson's Equation for Rapidly Rotating Stars
- S. Jackson (1970, ApJ, 161, 579 - 585) — Rapidly Rotating Stars. V. The Coupling of the Henyey and the Self-Consistent Methods
- J. P. Ostriker & J. W.-K. Mark (1968, ApJ, 151, 1075 - 1088) — Rapidly Rotating Stars. I. The Self-Consistent-Field Method
- R. A. James (1964, ApJ, 140, 552 - 582) — The Structure and Stability of Rotating Gas Masses
Early Eriguchi Applications
- Y. Eriguchi & E. Mueller (1985, A&A, 147, 161 - 168) — Equilibrium Models of Differentially Rotating Polytropes and the Collapse of Rotating Stellar Cores
- I. Hachisu & Y. Eriguchi (1984, PASJ, 36, 497 - 503) — Bifurcation Points on the Maclaurin Sequence
- I. Hachisu & Y. Eriguchi (1984, PASJ, 36, 259 - 276) — Binary Fluid Star
- I. Hachisu & Y. Eriguchi (1984, PASJ, 36, 239 - 257) — Fission of Dumbbell Equilibrium and Binary State of Rapidly Rotating Polytropes
- I. Hachisu & Y. Eriguchi (1983, MNRAS, 204, 583 - 589) — Bifurcations and Phase Transitions of Self-Gravitating and Uniformly Rotating Fluid
- I. Hachisu & Y. Eriguchi (1982, Prog. Theor. Phys., 68, 206 - 221) — Bifurcation and Fission of Three Dimensional, Rigidly Rotating and Self-Gravitating Polytropes
- I. Hachisu, Y. Eriguchi, & D. Sugimoto (1982, Prog. Theor. Phys., 68, 191 - 205) — Rapidly Rotating Polytropes and Concave Hamburger Equilibrium
- Y. Eriguchi & D. Sugimoto (1981, Prog. Theor. Phys., 65, 1870 - 1875) — Another Equilibrium Sequence of Self-Gravitating and Rotating Incompressible Fluid
- T. Fukushima, Y. Eriguchi, D. Sugimoto, G. S. Bisnovatyi-Kogan (1980, Prog. Theor. Phys., 63, 1957 - 1970) — Concave Hamburger Equilibrium of Rotating Bodies
Other Example Applications
- E. Mueller & Y. Eriguchi (1985, A&A, 152, 325 - 335) — Equilibrium Models of Differentially Rotating, Completely Catalyzed, Zero-Temperature Configurations With Central Densities Intermediate to White Dwarf and Neutron Star Densities
- J. R. Ipser & R. A. Managan (1981, ApJ, 250, 362 - 372) — On the Existence and Structure of Inhomogeneous Analogs of the Dedekind and Jacobi Ellipsoids
- C. T. Cunningham (1977, ApJ, 211, 568 - 578) — Rapidly Rotating Spheroids of Polytropic Index n = 1
- R. H. Durisen (1975, ApJ, 199, 179 - 183) — Upper Mass Limits for Stable Rotating White Dwarfs
- P. Bodenheimer & J. P. Ostriker (1973, ApJ, 180, 159 - 170) — Rapidly Rotating Stars. VIII. Zero-Viscosity Polytropic Sequences'
- P. Bodenheimer (1971, ApJ, 167, 153 - 163) — Rapidly Rotating Stars. VII. Effects of Angular Momentum on Upper-Main-Sequence Models
- P. Bodenheimer & J. P. Ostriker (1970, ApJ, 161, 1101 - 1113) — Rapidly Rotating Stars. VI. Pre-Main-Sequence Evolution of Massive Stars
- R. Kippenhahn & H.-C. Thomas (1970) in 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 — Stellar Rotation ==> Purchase proceedings from Springer, from Australia, or from Google
- In the introductory section of his paper, S. Jackson (1970) references this article by Kippenhahn & Thomas in the context of uniformly rotating, and therefore only mildly distorted, structures.
- J. P. Ostriker & J. L. Tassoul (1969, ApJ, 155, 987 - ) — On the Oscillations and Stability of Rotating Stellar Models. II. Rapidly Rotating White Dwarfs
- M. J. Clement (1969, ApJ, 156, 1051 - 1068) — Differential Rotation in Stars on the Upper Main Sequence
- J. W.-K. Mark (1968, ApJ, 154, 627 - ) — Rapidly Rotating Stars. III. Massive Main-Sequence Stars
- J. P. Ostriker & P. Bodenheimer (1968, ApJ, 151, 1089 - ) — Rapidly Rotating Stars. II. Massive White Dwarfs
- J. Faulkner, I. W. Roxburgh, & P. A. Strittmatter (1968, ApJ, 151, 203 - 216) — Uniformly Rotating Main-Sequence Stars
- R. Stoeckly (1965, ApJ, 142, 208 - 228) — Polytropic Models with Fast, Non-Uniform Rotation
- In the introductory section of his paper, 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
- L. G. Henyey, L. Wilets, K. H. Böhm, R. Lelevier, & R. D. Levee (1959, ApJ, 129, 628 - ) — A Method for Automatic Computation of Stellar Evolution
- L. G. Henyey, J. E. Forbes, & N. L. Gould (1964, ApJ, 139, 306 - ) — A New Method of Automatic Computation of Stellar Evolution
|
|
|---|
|
© 2014 - 2021 by Joel E. Tohline |