Difference between revisions of "User:Tohline/ThreeDimensionalConfigurations/FerrersPotential"
Line 129: | Line 129: | ||
== | ==Derivation== | ||
<table border="1" cellpadding="8" width="80%" align="center"> | |||
<tr><td align="left"> | |||
Other references to Ferrers Potential: | |||
<ul> | |||
<li>[https://babel.hathitrust.org/cgi/pt?id=uc1.$b417536&view=1up&seq=15 Ferrers, N. M. (1877, Quarterly Journal of Pure and Applied Mathematics, 14, 1)] … from [[User:Tohline/Appendix/References#BT87|BT87]] References (p. 711)</li> | |||
<li>[https://docs.galpy.org/en/v1.6.0/reference/potentialferrers.html Galpy methods]</li> | |||
<li>[https://www.mdpi.com/2075-4434/5/4/101 Lucas Antonio Caritá, Irapuan Rodriguez, & I. Puerari (2017)] <i>Explicit Second Partial Derivatives of the Ferrers Potential</i></li> | |||
<li>[https://www.researchgate.net/publication/225744594_Anisotropic_and_inhomogeneous_S-type_Riemann_ellipsoids_inside_spheroidal_halos_II Martin G. Abrahamyan (2006, Astrophysics, 49(3), 306-319)], ''Anisotropic and inhomogeneous S-type Riemann ellipsoids inside spheroidal halos. II''</li> | |||
</ul> | |||
</td> | |||
</tr> | |||
</table> | |||
Following §2.3.2 (beginning on p. 60) of [[User:Tohline/Appendix/References#BT87|BT87]], let's examine ''inhomogeneous'' configurations whose isodensity surfaces (including the surface, itself) are defined by triaxial ellipsoids on which the Cartesian coordinates <math>~(x_1, x_2, x_3)</math> satisfy the condition that, | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~m^2</math> | |||
</td> | |||
<td align="center"> | |||
<math>~\equiv</math> | |||
</td> | |||
<td align="left"> | |||
<math>~a_1^2 \sum_{i=1}^{3} \frac{x_i^2}{a_i^2} \, ,</math> | |||
</td> | |||
</tr> | |||
<tr><td align="center" colspan="3"> | |||
[ [[User:Tohline/Appendix/References#EFE|EFE]], <font color="#00CC00">Chapter 3, §20, p. 50, Eq. (75)</font> ]<br /> | |||
[ [[User:Tohline/Appendix/References#BT87|BT87]], <font color="#00CC00">§2.3.2, p. 61, Eq. (2-97)</font> ] | |||
</td></tr> | |||
</table> | |||
be constant. <span id="DensitySpecification">More specifically,</span> let's consider the case (related to the so-called ''Ferrers'' potentials) in which the configuration's density distribution is given by the expression, | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\rho(m^2)</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\rho_c \biggl[1 - \frac{m^2}{a_1^2}\biggr]^n </math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\rho_c \biggl[1 - | |||
\sum_{i=1}^{3} \frac{x_i^2}{a_i^2} \biggr]^n | |||
</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\rho_c \biggl[1 - \biggl( \frac{x^2}{a^2} + \frac{y^2}{b^2} + \frac{z^2}{c^2}\biggr) \biggr]^n | |||
\, .</math> | |||
</td> | |||
</tr> | |||
</table> | |||
<table border="1" align="center" width="80%" cellpadding="10"><tr><td align="left"> | |||
<font color="red">'''NOTE:'''</font> In our [[User:Tohline/ThreeDimensionalConfigurations/Challenges#Trial_.232|accompanying discussion]] of compressible analogues of Riemann S-type ellipsoids, we have discovered that — at least in the context of infinitesimally thin, nonaxisymmetric disks — this heterogeneous density profile can be nicely paired with an analytically expressible stream function, at least for the case where the integer exponent is, n = 1. | |||
</td></tr></table> | |||
According to Theorem 13 of [[User:Tohline/Appendix/References#EFE|EFE]] — see his Chapter 3, §20 (p. 53) — the potential at any point inside a triaxial ellipsoid with this specific density distribution is given by the expression, | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\Phi_\mathrm{grav}(\bold{x})</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
- \frac{\pi G \rho_c a_1 a_2 a_3}{(n+1)} \int_0^\infty \frac{ du}{\Delta } Q^{n+1} \, , | |||
</math> | |||
</td> | |||
</tr> | |||
<tr><td align="center" colspan="3"> | |||
[ [[User:Tohline/Appendix/References#EFE|EFE]], <font color="#00CC00">Chapter 3, §20, p. 53, Eq. (101)</font> ] | |||
</td></tr> | |||
</table> | |||
where, <math>~\Delta</math> has the same definition as [[User:Tohline/ThreeDimensionalConfigurations/HomogeneousEllipsoids#Derivation_of_Expressions_for_Ai|above]], and, | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~Q</math> | |||
</td> | |||
<td align="center"> | |||
<math>~\equiv</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
1 - \sum_{\ell = 1}^3 \frac{x_\ell^2}{ a_\ell^2 + u } \, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
For purposes of illustration, in what follows we will assume that, <math>~a_1 > a_2 > a_3</math>. | |||
===The Case Where n = 0=== | |||
When <math>~n = 0</math>, we have a uniform-density configuration, and the "interior" potential will be given by the expression, | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\Phi_\mathrm{grav}(\bold{x})</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
- \pi G \rho_c a_1 a_2 a_3 \int_0^\infty \frac{ du}{\Delta } \biggl[ 1 - \sum_{\ell = 1}^3 \frac{x_\ell^2}{ a_\ell^2 + u } \biggr] | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
- \pi G \rho_c a_1 a_2 a_3 \biggl\{ | |||
\int_0^\infty \frac{ du}{\Delta } | |||
- \int_0^\infty \frac{ du}{\Delta } \biggl( \frac{x^2}{ a_1^2 + u } \biggr) | |||
- \int_0^\infty \frac{ du}{\Delta } \biggl( \frac{y^2}{ a_2^2 + u } \biggr) | |||
- \int_0^\infty \frac{ du}{\Delta } \biggl( \frac{z^2}{ a_3^2 + u } \biggr) | |||
\biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
- \pi G \rho_c a_1 a_2 a_3 \biggl\{ | |||
\int_0^\infty \frac{ du}{\Delta } | |||
~ - ~x^2 \int_0^\infty \frac{ du}{\Delta (a_1^2 + u) } | |||
~ - ~y^2 \int_0^\infty \frac{ du}{\Delta (a_2^2 + u) } | |||
~ - ~ \int_0^\infty \frac{ du}{\Delta (a_3^2 + u) } | |||
\biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
-\pi G \rho_c \biggl[ I_\mathrm{BT} a_1^2 - \biggl(A_1 x^2 + A_2 y^2 +A_3 z^2 \biggr) \biggr] \, .</math> | |||
</td> | |||
</tr> | |||
</table> | |||
As a check, let's see if this scalar potential satisfies the differential form of the | |||
<div align="center"> | |||
<span id="PGE:Poisson"><font color="#770000">'''Poisson Equation'''</font></span><br /> | |||
{{User:Tohline/Math/EQ_Poisson01}} | |||
</div> | |||
<span id="SumTo2">Given that,</span> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\sum_{\ell = 1}^3 A_\ell</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~2 \, ,</math> | |||
</td> | |||
</tr> | |||
<tr><td align="center" colspan="3">[ [[User:Tohline/Appendix/References#EFE|EFE]], <font color="#00CC00">§21, Eq. (108)</font> ]</td></tr> | |||
</table> | |||
we find, | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\nabla^2\Phi_\mathrm{grav} = \biggl[\frac{\partial^2}{\partial x^2} + \frac{\partial^2}{\partial y^2} + \frac{\partial^2}{\partial z^2}\biggr]\Phi_\mathrm{grav}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
+ 2\pi G \rho_c (A_1 + A_2 + A_3) = 4\pi G\rho_c \, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
Q.E.D. | |||
===The Case Where n = 1=== | |||
When <math>~n = 1</math>, we have a specific heterogeneous density configuration, and the "interior" potential will be given by the expression, | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\frac{ \Phi_\mathrm{grav}(\bold{x})}{(-\pi G\rho_c)} </math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{1}{2} a_1 a_2 a_3 \int_0^\infty \frac{ du}{\Delta } \biggl[ 1 - \sum_{\ell = 1}^3 \frac{x_\ell^2}{ a_\ell^2 + u } \biggr]^2 | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{1}{2} a_1 a_2 a_3 \biggl\{ | |||
\int_0^\infty \frac{ du}{\Delta } \biggl[ 1 - \sum_{\ell = 1}^3 \frac{x_\ell^2}{ a_\ell^2 + u } \biggr] | |||
~- ~ x^2 \int_0^\infty \frac{ du}{\Delta (a_1^2 + u)} \biggl[ 1 - \sum_{\ell = 1}^3 \frac{x_\ell^2}{ a_\ell^2 + u } \biggr] | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
| |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
~- ~y^2 \int_0^\infty \frac{ du}{\Delta (a_2^2 + u)} \biggl[ 1 - \sum_{\ell = 1}^3 \frac{x_\ell^2}{ a_\ell^2 + u } \biggr] | |||
~ - ~z^2 \int_0^\infty \frac{ du}{\Delta (a_3^2 + u)} \biggl[ 1 - \sum_{\ell = 1}^3 \frac{x_\ell^2}{ a_\ell^2 + u } \biggr] | |||
\biggr\} \, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
The first definite-integral expression inside the curly braces is, to within a leading factor of <math>~\tfrac{1}{2}</math>, identical to the entire expression for the normalized potential that was derived in the case where n = 0. That is, we can write, | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\frac{ \Phi_\mathrm{grav}(\bold{x})}{(-\pi G\rho_c)} </math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{1}{2} \biggl[ I_\mathrm{BT} a_1^2 - \biggl(A_1 x^2 + A_2 y^2 +A_3 z^2 \biggr) \biggr] | |||
~- \frac{1}{2} a_1 a_2 a_3 \biggl\{ ~ x^2 \int_0^\infty \frac{ du}{\Delta (a_1^2 + u)} \biggl[ 1 - \sum_{\ell = 1}^3 \frac{x_\ell^2}{ a_\ell^2 + u } \biggr] | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
| |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
~+~y^2 \int_0^\infty \frac{ du}{\Delta (a_2^2 + u)} \biggl[ 1 - \sum_{\ell = 1}^3 \frac{x_\ell^2}{ a_\ell^2 + u } \biggr] | |||
~ + ~z^2 \int_0^\infty \frac{ du}{\Delta (a_3^2 + u)} \biggl[ 1 - \sum_{\ell = 1}^3 \frac{x_\ell^2}{ a_\ell^2 + u } \biggr] | |||
\biggr\} \, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
Then, from §22, p. 56 of [[User:Tohline/Appendix/References#EFE|EFE]], we see that, | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~a_1 a_2 a_3 \int_0^\infty \frac{ du}{\Delta (a_i^2 + u)} \biggl[ 1 - \sum_{\ell = 1}^3 \frac{x_\ell^2}{ a_\ell^2 + u } \biggr]</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\biggl( A_i - \sum_{\ell=1}^3 A_{i\ell} x_\ell^2 \biggr) \, . | |||
</math> | |||
</td> | |||
</tr> | |||
<tr><td align="center" colspan="3"> | |||
[ [[User:Tohline/Appendix/References#EFE|EFE]], <font color="#00CC00">Chapter 3, §22, p. 53, Eq. (125)</font> ] | |||
</td></tr> | |||
</table> | |||
<span id="GravFor1">Applying this result to each of the other three definite integrals gives us,</span> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\frac{ \Phi_\mathrm{grav}(\bold{x})}{(-\pi G\rho_c)} </math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{1}{2} \biggl[ I_\mathrm{BT} a_1^2 - \biggl(A_1 x^2 + A_2 y^2 +A_3 z^2 \biggr) \biggr] | |||
~- \frac{x^2}{2} \biggl( A_1 - \sum_{\ell=1}^3 A_{1\ell} x_\ell^2 \biggr) | |||
~- \frac{y^2}{2} \biggl( A_2 - \sum_{\ell=1}^3 A_{2\ell} x_\ell^2 \biggr) | |||
~- \frac{z^2}{2} \biggl( A_3 - \sum_{\ell=1}^3 A_{3\ell} x_\ell^2 \biggr) \, . | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{1}{2} \biggl[ I_\mathrm{BT} a_1^2 - \biggl(A_1 x^2 + A_2 y^2 +A_3 z^2 \biggr) \biggr] | |||
~- \frac{x^2}{2} \biggl[ A_1 - \biggl( A_{11}x^2 + A_{12}y^2 + A_{13}z^2 \biggr) \biggr] | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
| |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
~- \frac{y^2}{2} \biggl[ A_2 - \biggl( A_{21}x^2 + A_{22}y^2 + A_{23}z^2 \biggr) \biggr] | |||
~- \frac{z^2}{2} \biggl[ A_3 - \biggl( A_{31}x^2 + A_{32}y^2 + A_{33}z^2 \biggr) \biggr] | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{1}{2} I_\mathrm{BT} a_1^2 | |||
- \biggl(A_1 x^2 + A_2 y^2 +A_3 z^2 \biggr) | |||
~+ \biggl( A_{12} x^2y^2 + A_{13} x^2z^2 + A_{23} y^2z^2\biggr) | |||
~+ \frac{1}{2} \biggl(A_{11}x^4 + A_{22}y^4 + A_{33}z^4 \biggr) | |||
\, , | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
where, | |||
<table border="1" align="center" cellpadding="10" width="80%"> | |||
<tr> | |||
<td align="center" width="50%"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr><td align="center" colspan="3">for <math>~i \ne j</math></td></tr> | |||
<tr> | |||
<td align="right"> | |||
<math>~A_{ij}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~\equiv</math> | |||
</td> | |||
<td align="left"> | |||
<math>~-\frac{A_i-A_j}{(a_i^2 - a_j^2)} </math> | |||
</td> | |||
</tr> | |||
<tr><td align="center" colspan="3">[ [[User:Tohline/Appendix/References#EFE|EFE]], <font color="#00CC00">§21, Eq. (107)</font> ]</td></tr> | |||
</table> | |||
</td> | |||
<td align="center" width="50%"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr><td align="center" colspan="3">for <math>~i = j</math></td></tr> | |||
<tr> | |||
<td align="right"> | |||
<math>~2A_{ii} + \sum_{\ell = 1}^3 A_{i\ell}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\frac{2}{a_i} </math> | |||
</td> | |||
</tr> | |||
<tr><td align="center" colspan="3">[ [[User:Tohline/Appendix/References#EFE|EFE]], <font color="#00CC00">§21, Eq. (109)</font> ]</td></tr> | |||
</table> | |||
</td> | |||
</tr> | |||
</table> | |||
and we have made use of the symmetry relation, <math>~A_{ij} = A_{ji}</math>. Again, as a check, let's see if this scalar potential satisfies the differential form of the | |||
<div align="center"> | |||
<span id="PGE:Poisson"><font color="#770000">'''Poisson Equation'''</font></span><br /> | |||
{{User:Tohline/Math/EQ_Poisson01}} | |||
</div> | |||
We find, | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\nabla^2 \biggl[ \frac{\Phi_\mathrm{grav}}{-2\pi G \rho_c} \biggr] </math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{1}{2}\biggl[\frac{\partial^2}{\partial x^2} + \frac{\partial^2}{\partial y^2} + \frac{\partial^2}{\partial z^2}\biggr] | |||
\biggl[- \biggl(A_1 x^2 + A_2 y^2 +A_3 z^2 \biggr) | |||
~+ \biggl( A_{12} x^2y^2 + A_{13} x^2z^2 + A_{23} y^2z^2\biggr) | |||
~+ \frac{1}{2} \biggl(A_{11}x^4 + A_{22}y^4 + A_{33}z^4 \biggr) | |||
\biggr] | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{\partial}{\partial x} \biggl[- A_1 x ~+ A_{12} x y^2 + A_{13} x z^2 ~+ A_{11}x^3 \biggr] | |||
+\frac{\partial}{\partial y} \biggl[- A_2 y ~+ A_{12} x^2y + A_{23} y z^2~+ A_{22}y^3 \biggr] | |||
+ \frac{\partial}{\partial z}\biggl[- A_3 z ~+ A_{13} x^2z + A_{23} y^2z~+ A_{33}z^3 \biggr] | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\biggl[- A_1 + A_{12} y^2 + A_{13} z^2 ~+ 3A_{11}x^2 \biggr] | |||
+ \biggl[- A_2 + A_{12} x^2 + A_{23} z^2~+ 3A_{22}y^2 \biggr] | |||
+ \biggl[- A_3 + A_{13} x^2 + A_{23} y^2~+ 3A_{33}z^2 \biggr] | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ - (A_1 + A_2 + A_3) + x^2(3A_{11} + A_{12} + A_{13}) + y^2( 3A_{22} + A_{12} + A_{23}) + z^2( 3A_{33} + A_{13} + A_{23})\, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
In addition to recognizing, as [[#SumTo2|stated above]], that <math>~(A_1 + A_2 + A_3) = 2</math>, and making explicit use of the relation, | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~2A_{ii} + \sum_{\ell = 1}^3 A_{i\ell}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\frac{2}{a_i} \, ,</math> | |||
</td> | |||
</tr> | |||
</table> | |||
this last expression can be simplified to discover that, | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\nabla^2 \biggl[ \frac{\Phi_\mathrm{grav}}{-2\pi G \rho_c} \biggr] </math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ - (2) + \frac{2x^2}{a_1^2} + \frac{2y^2}{a_2^2} + \frac{2z^2}{a_3^2} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
<math>~\Rightarrow ~~~ \nabla^2 \Phi_\mathrm{grav}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ 4\pi G \rho_c \biggl[ 1 - \biggl( \frac{x^2}{a_1^2} + \frac{y^2}{a_2^2} + \frac{z^2}{a_3^2}\biggr) \biggr] \, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
This does indeed demonstrate that the derived gravitational potential is consistent with [[#DensitySpecification|our selected mass distribution in the case where n = 1]], namely, | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\rho</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ \rho_c \biggl[ 1 - \biggl( \frac{x^2}{a_1^2} + \frac{y^2}{a_2^2} + \frac{z^2}{a_3^2}\biggr) \biggr] \, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
Q.E.D. | |||
=See Also= | =See Also= |
Revision as of 03:39, 11 December 2020
Ferrers (1877) Gravitational Potential for Inhomogeneous Ellipsoids
| Tiled Menu | Tables of Content | Banner Video | Tohline Home Page | |
In an accompanying chapter titled, Properties of Homogeneous Ellipsoids (1), we have shown how analytic expressions may be derived for the gravitational potential inside of uniform-density ellipsoids. In that discussion, we largely followed the derivations of EFE. In the latter part of the nineteenth-century, N. M. Ferrers, (1877, Quarterly Journal of Pure and Applied Mathematics, 14, 1 - 22) showed that very similar analytic expressions can be derived for ellipsoids that have certain, specific inhomogeneous mass distributions. Here we specifically discuss the case of configurations that exhibit concentric ellipsoidal iso-density surfaces of the form,
SUMMARY — copied from accompanying, Trial #2 Discussion | ||||||||||||||||||||||||||||||||
After studying the relevant sections of both EFE and BT87 — this is an example of a heterogeneous density distribution whose gravitational potential has an analytic prescription. As is discussed in a separate chapter, the potential that it generates is sometimes referred to as a Ferrers potential, for the exponent, n = 1. In our accompanying discussion we find that,
where,
More specifically, in the three cases where the indices, <math>~i=j</math>,
|
Derivation
Other references to Ferrers Potential:
|
Following §2.3.2 (beginning on p. 60) of BT87, let's examine inhomogeneous configurations whose isodensity surfaces (including the surface, itself) are defined by triaxial ellipsoids on which the Cartesian coordinates <math>~(x_1, x_2, x_3)</math> satisfy the condition that,
<math>~m^2</math> |
<math>~\equiv</math> |
<math>~a_1^2 \sum_{i=1}^{3} \frac{x_i^2}{a_i^2} \, ,</math> |
[ EFE, Chapter 3, §20, p. 50, Eq. (75) ] |
be constant. More specifically, let's consider the case (related to the so-called Ferrers potentials) in which the configuration's density distribution is given by the expression,
<math>~\rho(m^2)</math> |
<math>~=</math> |
<math>~\rho_c \biggl[1 - \frac{m^2}{a_1^2}\biggr]^n </math> |
<math>~=</math> |
<math>~ \rho_c \biggl[1 - \sum_{i=1}^{3} \frac{x_i^2}{a_i^2} \biggr]^n </math> |
<math>~=</math> |
<math>~ \rho_c \biggl[1 - \biggl( \frac{x^2}{a^2} + \frac{y^2}{b^2} + \frac{z^2}{c^2}\biggr) \biggr]^n \, .</math> |
NOTE: In our accompanying discussion of compressible analogues of Riemann S-type ellipsoids, we have discovered that — at least in the context of infinitesimally thin, nonaxisymmetric disks — this heterogeneous density profile can be nicely paired with an analytically expressible stream function, at least for the case where the integer exponent is, n = 1. |
According to Theorem 13 of EFE — see his Chapter 3, §20 (p. 53) — the potential at any point inside a triaxial ellipsoid with this specific density distribution is given by the expression,
<math>~\Phi_\mathrm{grav}(\bold{x})</math> |
<math>~=</math> |
<math>~ - \frac{\pi G \rho_c a_1 a_2 a_3}{(n+1)} \int_0^\infty \frac{ du}{\Delta } Q^{n+1} \, , </math> |
[ EFE, Chapter 3, §20, p. 53, Eq. (101) ] |
where, <math>~\Delta</math> has the same definition as above, and,
<math>~Q</math> |
<math>~\equiv</math> |
<math>~ 1 - \sum_{\ell = 1}^3 \frac{x_\ell^2}{ a_\ell^2 + u } \, . </math> |
For purposes of illustration, in what follows we will assume that, <math>~a_1 > a_2 > a_3</math>.
The Case Where n = 0
When <math>~n = 0</math>, we have a uniform-density configuration, and the "interior" potential will be given by the expression,
<math>~\Phi_\mathrm{grav}(\bold{x})</math> |
<math>~=</math> |
<math>~ - \pi G \rho_c a_1 a_2 a_3 \int_0^\infty \frac{ du}{\Delta } \biggl[ 1 - \sum_{\ell = 1}^3 \frac{x_\ell^2}{ a_\ell^2 + u } \biggr] </math> |
|
<math>~=</math> |
<math>~ - \pi G \rho_c a_1 a_2 a_3 \biggl\{ \int_0^\infty \frac{ du}{\Delta } - \int_0^\infty \frac{ du}{\Delta } \biggl( \frac{x^2}{ a_1^2 + u } \biggr) - \int_0^\infty \frac{ du}{\Delta } \biggl( \frac{y^2}{ a_2^2 + u } \biggr) - \int_0^\infty \frac{ du}{\Delta } \biggl( \frac{z^2}{ a_3^2 + u } \biggr) \biggr\} </math> |
|
<math>~=</math> |
<math>~ - \pi G \rho_c a_1 a_2 a_3 \biggl\{ \int_0^\infty \frac{ du}{\Delta } ~ - ~x^2 \int_0^\infty \frac{ du}{\Delta (a_1^2 + u) } ~ - ~y^2 \int_0^\infty \frac{ du}{\Delta (a_2^2 + u) } ~ - ~ \int_0^\infty \frac{ du}{\Delta (a_3^2 + u) } \biggr\} </math> |
|
<math>~=</math> |
<math>~ -\pi G \rho_c \biggl[ I_\mathrm{BT} a_1^2 - \biggl(A_1 x^2 + A_2 y^2 +A_3 z^2 \biggr) \biggr] \, .</math> |
As a check, let's see if this scalar potential satisfies the differential form of the
Given that,
<math>~\sum_{\ell = 1}^3 A_\ell</math> |
<math>~=</math> |
<math>~2 \, ,</math> |
[ EFE, §21, Eq. (108) ] |
we find,
<math>~\nabla^2\Phi_\mathrm{grav} = \biggl[\frac{\partial^2}{\partial x^2} + \frac{\partial^2}{\partial y^2} + \frac{\partial^2}{\partial z^2}\biggr]\Phi_\mathrm{grav}</math> |
<math>~=</math> |
<math>~ + 2\pi G \rho_c (A_1 + A_2 + A_3) = 4\pi G\rho_c \, . </math> |
Q.E.D.
The Case Where n = 1
When <math>~n = 1</math>, we have a specific heterogeneous density configuration, and the "interior" potential will be given by the expression,
<math>~\frac{ \Phi_\mathrm{grav}(\bold{x})}{(-\pi G\rho_c)} </math> |
<math>~=</math> |
<math>~ \frac{1}{2} a_1 a_2 a_3 \int_0^\infty \frac{ du}{\Delta } \biggl[ 1 - \sum_{\ell = 1}^3 \frac{x_\ell^2}{ a_\ell^2 + u } \biggr]^2 </math> |
|
<math>~=</math> |
<math>~ \frac{1}{2} a_1 a_2 a_3 \biggl\{ \int_0^\infty \frac{ du}{\Delta } \biggl[ 1 - \sum_{\ell = 1}^3 \frac{x_\ell^2}{ a_\ell^2 + u } \biggr] ~- ~ x^2 \int_0^\infty \frac{ du}{\Delta (a_1^2 + u)} \biggl[ 1 - \sum_{\ell = 1}^3 \frac{x_\ell^2}{ a_\ell^2 + u } \biggr] </math> |
|
|
<math>~ ~- ~y^2 \int_0^\infty \frac{ du}{\Delta (a_2^2 + u)} \biggl[ 1 - \sum_{\ell = 1}^3 \frac{x_\ell^2}{ a_\ell^2 + u } \biggr] ~ - ~z^2 \int_0^\infty \frac{ du}{\Delta (a_3^2 + u)} \biggl[ 1 - \sum_{\ell = 1}^3 \frac{x_\ell^2}{ a_\ell^2 + u } \biggr] \biggr\} \, . </math> |
The first definite-integral expression inside the curly braces is, to within a leading factor of <math>~\tfrac{1}{2}</math>, identical to the entire expression for the normalized potential that was derived in the case where n = 0. That is, we can write,
<math>~\frac{ \Phi_\mathrm{grav}(\bold{x})}{(-\pi G\rho_c)} </math> |
<math>~=</math> |
<math>~ \frac{1}{2} \biggl[ I_\mathrm{BT} a_1^2 - \biggl(A_1 x^2 + A_2 y^2 +A_3 z^2 \biggr) \biggr] ~- \frac{1}{2} a_1 a_2 a_3 \biggl\{ ~ x^2 \int_0^\infty \frac{ du}{\Delta (a_1^2 + u)} \biggl[ 1 - \sum_{\ell = 1}^3 \frac{x_\ell^2}{ a_\ell^2 + u } \biggr] </math> |
|
|
<math>~ ~+~y^2 \int_0^\infty \frac{ du}{\Delta (a_2^2 + u)} \biggl[ 1 - \sum_{\ell = 1}^3 \frac{x_\ell^2}{ a_\ell^2 + u } \biggr] ~ + ~z^2 \int_0^\infty \frac{ du}{\Delta (a_3^2 + u)} \biggl[ 1 - \sum_{\ell = 1}^3 \frac{x_\ell^2}{ a_\ell^2 + u } \biggr] \biggr\} \, . </math> |
Then, from §22, p. 56 of EFE, we see that,
<math>~a_1 a_2 a_3 \int_0^\infty \frac{ du}{\Delta (a_i^2 + u)} \biggl[ 1 - \sum_{\ell = 1}^3 \frac{x_\ell^2}{ a_\ell^2 + u } \biggr]</math> |
<math>~=</math> |
<math>~ \biggl( A_i - \sum_{\ell=1}^3 A_{i\ell} x_\ell^2 \biggr) \, . </math> |
[ EFE, Chapter 3, §22, p. 53, Eq. (125) ] |
Applying this result to each of the other three definite integrals gives us,
<math>~\frac{ \Phi_\mathrm{grav}(\bold{x})}{(-\pi G\rho_c)} </math> |
<math>~=</math> |
<math>~ \frac{1}{2} \biggl[ I_\mathrm{BT} a_1^2 - \biggl(A_1 x^2 + A_2 y^2 +A_3 z^2 \biggr) \biggr] ~- \frac{x^2}{2} \biggl( A_1 - \sum_{\ell=1}^3 A_{1\ell} x_\ell^2 \biggr) ~- \frac{y^2}{2} \biggl( A_2 - \sum_{\ell=1}^3 A_{2\ell} x_\ell^2 \biggr) ~- \frac{z^2}{2} \biggl( A_3 - \sum_{\ell=1}^3 A_{3\ell} x_\ell^2 \biggr) \, . </math> |
|
<math>~=</math> |
<math>~ \frac{1}{2} \biggl[ I_\mathrm{BT} a_1^2 - \biggl(A_1 x^2 + A_2 y^2 +A_3 z^2 \biggr) \biggr] ~- \frac{x^2}{2} \biggl[ A_1 - \biggl( A_{11}x^2 + A_{12}y^2 + A_{13}z^2 \biggr) \biggr] </math> |
|
|
<math>~ ~- \frac{y^2}{2} \biggl[ A_2 - \biggl( A_{21}x^2 + A_{22}y^2 + A_{23}z^2 \biggr) \biggr] ~- \frac{z^2}{2} \biggl[ A_3 - \biggl( A_{31}x^2 + A_{32}y^2 + A_{33}z^2 \biggr) \biggr] </math> |
|
<math>~=</math> |
<math>~ \frac{1}{2} I_\mathrm{BT} a_1^2 - \biggl(A_1 x^2 + A_2 y^2 +A_3 z^2 \biggr) ~+ \biggl( A_{12} x^2y^2 + A_{13} x^2z^2 + A_{23} y^2z^2\biggr) ~+ \frac{1}{2} \biggl(A_{11}x^4 + A_{22}y^4 + A_{33}z^4 \biggr) \, , </math> |
where,
|
|
and we have made use of the symmetry relation, <math>~A_{ij} = A_{ji}</math>. Again, as a check, let's see if this scalar potential satisfies the differential form of the
We find,
<math>~\nabla^2 \biggl[ \frac{\Phi_\mathrm{grav}}{-2\pi G \rho_c} \biggr] </math> |
<math>~=</math> |
<math>~ \frac{1}{2}\biggl[\frac{\partial^2}{\partial x^2} + \frac{\partial^2}{\partial y^2} + \frac{\partial^2}{\partial z^2}\biggr] \biggl[- \biggl(A_1 x^2 + A_2 y^2 +A_3 z^2 \biggr) ~+ \biggl( A_{12} x^2y^2 + A_{13} x^2z^2 + A_{23} y^2z^2\biggr) ~+ \frac{1}{2} \biggl(A_{11}x^4 + A_{22}y^4 + A_{33}z^4 \biggr) \biggr] </math> |
|
<math>~=</math> |
<math>~ \frac{\partial}{\partial x} \biggl[- A_1 x ~+ A_{12} x y^2 + A_{13} x z^2 ~+ A_{11}x^3 \biggr] +\frac{\partial}{\partial y} \biggl[- A_2 y ~+ A_{12} x^2y + A_{23} y z^2~+ A_{22}y^3 \biggr] + \frac{\partial}{\partial z}\biggl[- A_3 z ~+ A_{13} x^2z + A_{23} y^2z~+ A_{33}z^3 \biggr] </math> |
|
<math>~=</math> |
<math>~ \biggl[- A_1 + A_{12} y^2 + A_{13} z^2 ~+ 3A_{11}x^2 \biggr] + \biggl[- A_2 + A_{12} x^2 + A_{23} z^2~+ 3A_{22}y^2 \biggr] + \biggl[- A_3 + A_{13} x^2 + A_{23} y^2~+ 3A_{33}z^2 \biggr] </math> |
|
<math>~=</math> |
<math>~ - (A_1 + A_2 + A_3) + x^2(3A_{11} + A_{12} + A_{13}) + y^2( 3A_{22} + A_{12} + A_{23}) + z^2( 3A_{33} + A_{13} + A_{23})\, . </math> |
In addition to recognizing, as stated above, that <math>~(A_1 + A_2 + A_3) = 2</math>, and making explicit use of the relation,
<math>~2A_{ii} + \sum_{\ell = 1}^3 A_{i\ell}</math> |
<math>~=</math> |
<math>~\frac{2}{a_i} \, ,</math> |
this last expression can be simplified to discover that,
<math>~\nabla^2 \biggl[ \frac{\Phi_\mathrm{grav}}{-2\pi G \rho_c} \biggr] </math> |
<math>~=</math> |
<math>~ - (2) + \frac{2x^2}{a_1^2} + \frac{2y^2}{a_2^2} + \frac{2z^2}{a_3^2} </math> |
<math>~\Rightarrow ~~~ \nabla^2 \Phi_\mathrm{grav}</math> |
<math>~=</math> |
<math>~ 4\pi G \rho_c \biggl[ 1 - \biggl( \frac{x^2}{a_1^2} + \frac{y^2}{a_2^2} + \frac{z^2}{a_3^2}\biggr) \biggr] \, . </math> |
This does indeed demonstrate that the derived gravitational potential is consistent with our selected mass distribution in the case where n = 1, namely,
<math>~\rho</math> |
<math>~=</math> |
<math>~ \rho_c \biggl[ 1 - \biggl( \frac{x^2}{a_1^2} + \frac{y^2}{a_2^2} + \frac{z^2}{a_3^2}\biggr) \biggr] \, . </math> |
Q.E.D.
See Also
- Our Speculation6 effort to develop a "Concentric Ellipsoidal (T6) Coordinate System."
- Challenges Constructing Ellipsoidal-Like Configurations
© 2014 - 2021 by Joel E. Tohline |