ParabolicDensity/Axisymmetric/Structure

Parabolic Density Distribution

Part I:   Gravitational Potential

Part II:   Spherical Structures

Part III:   Axisymmetric Equilibrium Structures  Old: 1st thru 7th tries  Old: 8th thru 10th tries

Part IV:   Triaxial Equilibrium Structures (Exploration)

Axisymmetric (Oblate) Equilibrium Structures

Tentative Summary

Known Relations

Density:	$\frac{\rho (\varpi ,z)}{{\rho }_c}$	$=$	$[1-{\chi }^2-{\zeta }^2(1-e^2{)}^{-1}],$
Gravitational Potential:	$\frac{{\mathrm{\Phi }}_{\mathrm{g}\mathrm{r}\mathrm{a}\mathrm{v}}(\varpi ,z)}{(-\pi G{\rho }_c{a}_{\ell }^2)}$	$=$	$\frac{1}{2}{I}_{\mathrm{B}\mathrm{T}}-A_{\ell }{\chi }^2-A_s{\zeta }^2+\frac{1}{2}[(A_{ss}{a}_{\ell }^2){\zeta }^4+2(A_{\ell s}{a}_{\ell }^2){\chi }^2{\zeta }^2+(A_{\ell \ell }{a}_{\ell }^2){\chi }^4].$
Vertical Pressure Gradient:	$\left[\frac{1}{(\pi G{\rho }_c^2{a}_{\ell }^2)}\right]\frac{\partial P}{\partial \zeta }$	$=$	$\frac{\rho }{{\rho }_c}\cdot [2A_{\ell s}{a}_{\ell }^2{\chi }^2\zeta -2A_s\zeta +2A_{ss}{a}_{\ell }^2{\zeta }^3]$

where, $\chi \equiv \varpi /a_{\ell }$ and $\zeta \equiv z/a_{\ell }$, and the relevant index symbol expressions are:

$I_{\mathrm{B}\mathrm{T}}$	$=$	$2A_{\ell }+A_s(1-e^2)=2(1-e^2{)}^{1/2}[\frac{{\sin }^{-1}e}{e}];$
$A_{\ell }$	$=$	$\frac{1}{e^2}[\frac{{\sin }^{-1}e}{e}-(1-e^2{)}^{1/2}](1-e^2{)}^{1/2};$
$A_s$	$=$	$\frac{2}{e^2}[(1-e^2{)}^{-1/2}-\frac{{\sin }^{-1}e}{e}](1-e^2{)}^{1/2};$
$a_{\ell }^2{A}_{\ell \ell }$	$=$	$\frac{1}{4e^4}\{-(3+2e^2)(1-e^2)+3(1-e^2{)}^{1/2}\left[\frac{{\sin }^{-1}e}{e}\right]\};$
$\frac{3}{2}{a}_{\ell }^2{A}_{ss}$	$=$	$\frac{(4e^2-3)}{e^4(1-e^2)}+\frac{3(1-e^2{)}^{1/2}}{e^4}\left[\frac{{\sin }^{-1}e}{e}\right];$
$a_{\ell }^2{A}_{\ell s}$	$=$	$\frac{1}{e^4}\{(3-e^2)-3(1-e^2{)}^{1/2}\left[\frac{{\sin }^{-1}e}{e}\right]\},$

where the eccentricity,

$e\equiv [1-(\frac{a_s}{a_{\ell }}{)}^2{]}^{1/2}.$

Drawing from our separate "6^th Try" discussion — and as has been highlighted here for example — for the axisymmetric configurations under consideration, the ${\hat{e}}_z$ and ${\hat{e}}_{\varpi }$ components of the Euler equation become, respectively,

${\hat{e}}_z$:     $0$ = $[\frac{1}{\rho }\frac{\partial P}{\partial z}+\frac{\partial \mathrm{\Phi }}{\partial z}]$ ${\hat{e}}_{\varpi }$:     $\frac{j^2}{{\varpi }^3}$ = $[\frac{1}{\rho }\frac{\partial P}{\partial \varpi }+\frac{\partial \mathrm{\Phi }}{\partial \varpi }]$

Multiplying through by length $(a_{\ell })$ and dividing through by the square of the velocity $(\pi G{\rho }_c{a}_{\ell }^2)$, we have,

${\hat{e}}_z$:	$0$	=	$[\frac{1}{\rho }\frac{\partial P}{\partial z}+\frac{\partial \mathrm{\Phi }}{\partial z}]\frac{a_{\ell }}{(\pi G{\rho }_c{a}_{\ell }^2)}$
		=	$\frac{{\rho }_c}{\rho }\cdot \frac{\partial }{\partial \zeta }\left[\frac{P}{(\pi G{\rho }_c^2{a}_{\ell }^2)}\right]-\frac{\partial }{\partial \zeta }\left[\frac{\mathrm{\Phi }}{(-\pi G{\rho }_c{a}_{\ell }^2)}\right]$
${\hat{e}}_{\varpi }$:	$\frac{j^2}{{\varpi }^3}\cdot \frac{a_{\ell }}{(\pi G{\rho }_c{a}_{\ell }^2)}$	=	$[\frac{1}{\rho }\frac{\partial P}{\partial \varpi }+\frac{\partial \mathrm{\Phi }}{\partial \varpi }]\frac{a_{\ell }}{(\pi G{\rho }_c{a}_{\ell }^2)}$
	$\Rightarrow \frac{1}{{\chi }^3}\cdot \frac{j^2}{(\pi G{\rho }_c{a}_{\ell }^4)}$	=	$\frac{{\rho }_c}{\rho }\cdot \frac{\partial }{\partial \chi }\left[\frac{P}{(\pi G{\rho }_c^2{a}_{\ell }^2)}\right]-\frac{\partial }{\partial \chi }\left[\frac{\mathrm{\Phi }}{(-\pi G{\rho }_c{a}_{\ell }^2)}\right]$

9^th Try

Starting Key Relations

Density:	$\frac{\rho (\varpi ,z)}{{\rho }_c}$	$=$	$[1-{\chi }^2-{\zeta }^2(1-e^2{)}^{-1}],$
Gravitational Potential:	$\frac{{\mathrm{\Phi }}_{\mathrm{g}\mathrm{r}\mathrm{a}\mathrm{v}}(\varpi ,z)}{(-\pi G{\rho }_c{a}_{\ell }^2)}$	$=$	$\frac{1}{2}{I}_{\mathrm{B}\mathrm{T}}-A_{\ell }{\chi }^2-A_s{\zeta }^2+\frac{1}{2}[(A_{ss}{a}_{\ell }^2){\zeta }^4+2(A_{\ell s}{a}_{\ell }^2){\chi }^2{\zeta }^2+(A_{\ell \ell }{a}_{\ell }^2){\chi }^4].$
Vertical Pressure Gradient:	$\left[\frac{1}{(\pi G{\rho }_c^2{a}_{\ell }^2)}\right]\frac{\partial P}{\partial \zeta }$	$=$	$\frac{\rho }{{\rho }_c}\cdot [2A_{\ell s}{a}_{\ell }^2{\chi }^2\zeta -2A_s\zeta +2A_{ss}{a}_{\ell }^2{\zeta }^3]$

Play With Vertical Pressure Gradient

$\left[\frac{1}{(\pi G{\rho }_c^2{a}_{\ell }^2)}\right]\frac{\partial P}{\partial \zeta }$	$=$	$[1-{\chi }^2-{\zeta }^2(1-e^2{)}^{-1}\left]\right[2A_{\ell s}{a}_{\ell }^2{\chi }^2\zeta -2A_s\zeta +2A_{ss}{a}_{\ell }^2{\zeta }^3]$
	$=$	$[(2A_{\ell s}{a}_{\ell }^2{\chi }^2-2A_s)\zeta +2A_{ss}{a}_{\ell }^2{\zeta }^3]-{\chi }^2[(2A_{\ell s}{a}_{\ell }^2{\chi }^2-2A_s)\zeta +2A_{ss}{a}_{\ell }^2{\zeta }^3]-{\zeta }^2(1-e^2{)}^{-1}[(2A_{\ell s}{a}_{\ell }^2{\chi }^2-2A_s)\zeta +2A_{ss}{a}_{\ell }^2{\zeta }^3]$
	$=$	$(2A_{\ell s}{a}_{\ell }^2{\chi }^2-2A_s)\zeta +2A_{ss}{a}_{\ell }^2{\zeta }^3-(2A_{\ell s}{a}_{\ell }^2{\chi }^4-2A_s{\chi }^2)\zeta -2A_{ss}{a}_{\ell }^2{\chi }^2{\zeta }^3-(1-e^2{)}^{-1}[(2A_{\ell s}{a}_{\ell }^2{\chi }^2-2A_s){\zeta }^3+2A_{ss}{a}_{\ell }^2{\zeta }^5]$
	$=$	$[(2A_{\ell s}{a}_{\ell }^2{\chi }^2-2A_s)-(2A_{\ell s}{a}_{\ell }^2{\chi }^4-2A_s{\chi }^2)]\zeta +[2A_{ss}{a}_{\ell }^2-2A_{ss}{a}_{\ell }^2{\chi }^2-(1-e^2{)}^{-1}(2A_{\ell s}{a}_{\ell }^2{\chi }^2-2A_s)]{\zeta }^3+[-(1-e^2{)}^{-1}2A_{ss}{a}_{\ell }^2]{\zeta }^5.$

Integrate over $\zeta$ gives …

$\left[\frac{1}{(\pi G{\rho }_c^2{a}_{\ell }^2)}\right]\int \left[\frac{\partial P}{\partial \zeta }\right]d\zeta$	$=$	$[(A_{\ell s}{a}_{\ell }^2{\chi }^2-A_s)-(A_{\ell s}{a}_{\ell }^2{\chi }^4-A_s{\chi }^2)]{\zeta }^2+\frac{1}{2}[A_{ss}{a}_{\ell }^2-A_{ss}{a}_{\ell }^2{\chi }^2-(1-e^2{)}^{-1}(A_{\ell s}{a}_{\ell }^2{\chi }^2-A_s)]{\zeta }^4+\frac{1}{3}[-(1-e^2{)}^{-1}{A}_{ss}{a}_{\ell }^2]{\zeta }^6+\mathrm{c}\mathrm{o}\mathrm{n}\mathrm{s}\mathrm{t}$
	$=$	$[-A_s{\zeta }^2+\frac{1}{2}{A}_{ss}{a}_{\ell }^2{\zeta }^4+\frac{1}{2}(1-e^2{)}^{-1}{A}_s{\zeta }^4-\frac{1}{3}(1-e^2{)}^{-1}{A}_{ss}{a}_{\ell }^2{\zeta }^6]{\chi }^0+[A_{\ell s}{a}_{\ell }^2{\zeta }^2+A_s{\zeta }^2-\frac{1}{2}{A}_{ss}{a}_{\ell }^2{\zeta }^4-\frac{1}{2}(1-e^2{)}^{-1}(A_{\ell s}{a}_{\ell }^2{\zeta }^4)]{\chi }^2+[-A_{\ell s}{a}_{\ell }^2{\zeta }^2]{\chi }^4+\mathrm{c}\mathrm{o}\mathrm{n}\mathrm{s}\mathrm{t}\mathrm{.}$

Now Play With Radial Pressure Gradient

$\left[\frac{1}{(-\pi G{\rho }_c{a}_{\ell }^2)}\right]\frac{\partial \mathrm{\Phi }}{\partial \chi }$	$=$	$\frac{\rho }{{\rho }_c}\cdot \{-2A_{\ell }\chi +\frac{1}{2}[4(A_{\ell s}{a}_{\ell }^2){\zeta }^2\chi +4(A_{\ell \ell }{a}_{\ell }^2){\chi }^3\left]\right\}$
	$=$	$2[1-{\chi }^2-{\zeta }^2(1-e^2{)}^{-1}\left]\right[(A_{\ell s}{a}_{\ell }^2{\zeta }^2-A_{\ell })\chi +A_{\ell \ell }{a}_{\ell }^2{\chi }^3]$
	$=$	$2[(A_{\ell s}{a}_{\ell }^2{\zeta }^2-A_{\ell })\chi +A_{\ell \ell }{a}_{\ell }^2{\chi }^3]-2{\chi }^2[(A_{\ell s}{a}_{\ell }^2{\zeta }^2-A_{\ell })\chi +A_{\ell \ell }{a}_{\ell }^2{\chi }^3]-2{\zeta }^2(1-e^2{)}^{-1}[(A_{\ell s}{a}_{\ell }^2{\zeta }^2-A_{\ell })\chi +A_{\ell \ell }{a}_{\ell }^2{\chi }^3]$
	$=$	$2(A_{\ell s}{a}_{\ell }^2{\zeta }^2-A_{\ell })\chi +2[A_{\ell \ell }{a}_{\ell }^2+(A_{\ell }-A_{\ell s}{a}_{\ell }^2{\zeta }^2)]{\chi }^3-2A_{\ell \ell }{a}_{\ell }^2{\chi }^5+2(1-e^2{)}^{-1}[(A_{\ell }{\zeta }^2-A_{\ell s}{a}_{\ell }^2{\zeta }^4)\chi -A_{\ell \ell }{a}_{\ell }^2{\zeta }^2{\chi }^3]$
	$=$	$2[(A_{\ell s}{a}_{\ell }^2{\zeta }^2-A_{\ell })+(1-e^2{)}^{-1}(A_{\ell }{\zeta }^2-A_{\ell s}{a}_{\ell }^2{\zeta }^4)]\chi +2[A_{\ell \ell }{a}_{\ell }^2+(A_{\ell }-A_{\ell s}{a}_{\ell }^2{\zeta }^2)-(1-e^2{)}^{-1}{A}_{\ell \ell }{a}_{\ell }^2{\zeta }^2]{\chi }^3-2A_{\ell \ell }{a}_{\ell }^2{\chi }^5$

10^th Try

Repeating Key Relations

Density:	$\frac{\rho (\varpi ,z)}{{\rho }_c}$	$=$	$[1-{\chi }^2-{\zeta }^2(1-e^2{)}^{-1}],$
Gravitational Potential:	$\frac{{\mathrm{\Phi }}_{\mathrm{g}\mathrm{r}\mathrm{a}\mathrm{v}}(\varpi ,z)}{(-\pi G{\rho }_c{a}_{\ell }^2)}$	$=$	$\frac{1}{2}{I}_{\mathrm{B}\mathrm{T}}-A_{\ell }{\chi }^2-A_s{\zeta }^2+\frac{1}{2}[(A_{ss}{a}_{\ell }^2){\zeta }^4+2(A_{\ell s}{a}_{\ell }^2){\chi }^2{\zeta }^2+(A_{\ell \ell }{a}_{\ell }^2){\chi }^4].$
Vertical Pressure Gradient:	$\left[\frac{1}{(\pi G{\rho }_c^2{a}_{\ell }^2)}\right]\frac{\partial P}{\partial \zeta }$	$=$	$\frac{\rho }{{\rho }_c}\cdot [2A_{\ell s}{a}_{\ell }^2{\chi }^2\zeta -2A_s\zeta +2A_{ss}{a}_{\ell }^2{\zeta }^3]$

From the above (9^th Try) examination of the vertical pressure gradient, we determined that a reasonably good approximation for the normalized pressure throughout the configuration is given by the expression,

$\left[\frac{1}{(\pi G{\rho }_c^2{a}_{\ell }^2)}\right]\int \left[\frac{\partial P}{\partial \zeta }\right]d\zeta$

$=$

$[-A_s{\zeta }^2+\frac{1}{2}{A}_{ss}{a}_{\ell }^2{\zeta }^4+\frac{1}{2}(1-e^2{)}^{-1}{A}_s{\zeta }^4-\frac{1}{3}(1-e^2{)}^{-1}{A}_{ss}{a}_{\ell }^2{\zeta }^6]{\chi }^0+[A_{\ell s}{a}_{\ell }^2{\zeta }^2+A_s{\zeta }^2-\frac{1}{2}{A}_{ss}{a}_{\ell }^2{\zeta }^4-\frac{1}{2}(1-e^2{)}^{-1}(A_{\ell s}{a}_{\ell }^2{\zeta }^4)]{\chi }^2+[-A_{\ell s}{a}_{\ell }^2{\zeta }^2]{\chi }^4+\mathrm{c}\mathrm{o}\mathrm{n}\mathrm{s}\mathrm{t}\mathrm{.}$

If we set $\chi =0$ — that is, if we look along the vertical axis — this approximation should be particularly good, resulting in the expression,

$P_z\equiv \left\{\right[\frac{1}{(\pi G{\rho }_c^2{a}_{\ell }^2)}]\int [\frac{\partial P}{\partial \zeta }]d\zeta {\}}_{\chi =0}$

$=$

$P_c^{*}-A_s{\zeta }^2+\frac{1}{2}{A}_{ss}{a}_{\ell }^2{\zeta }^4+\frac{1}{2}(1-e^2{)}^{-1}{A}_s{\zeta }^4-\frac{1}{3}(1-e^2{)}^{-1}{A}_{ss}{a}_{\ell }^2{\zeta }^6.$

Note that in the limit that $z\to a_s$ — that is, at the pole along the vertical (symmetry) axis where the $P_z$ should drop to zero — we should set $\zeta \to (1-e^2{)}^{1/2}$. This allows us to determine the central pressure. $P_c^{*}$ $=$ $A_s(1-e^2)-\frac{1}{2}{A}_{ss}{a}_{\ell }^2(1-e^2{)}^2-\frac{1}{2}(1-e^2{)}^{-1}{A}_s(1-e^2{)}^2+\frac{1}{3}(1-e^2{)}^{-1}{A}_{ss}{a}_{\ell }^2(1-e^2{)}^3$   $=$ $A_s(1-e^2)-\frac{1}{2}{A}_s(1-e^2)+\frac{1}{3}{A}_{ss}{a}_{\ell }^2(1-e^2{)}^2-\frac{1}{2}{A}_{ss}{a}_{\ell }^2(1-e^2{)}^2$   $=$ $\frac{1}{2}{A}_s(1-e^2)-\frac{1}{6}{A}_{ss}{a}_{\ell }^2(1-e^2{)}^2.$

This means that, along the vertical axis, the pressure gradient is,

$P_z\equiv \left\{\right[\frac{1}{(\pi G{\rho }_c^2{a}_{\ell }^2)}]\int [\frac{\partial P}{\partial \zeta }]d\zeta {\}}_{\chi =0}$

$=$

$P_c^{*}-A_s{\zeta }^2+\frac{1}{2}{A}_{ss}{a}_{\ell }^2{\zeta }^4+\frac{1}{2}(1-e^2{)}^{-1}{A}_s{\zeta }^4-\frac{1}{3}(1-e^2{)}^{-1}{A}_{ss}{a}_{\ell }^2{\zeta }^6.$

$\frac{\partial P_z}{\partial \zeta }$

$=$

$-2A_s\zeta +2A_{ss}{a}_{\ell }^2{\zeta }^3+2(1-e^2{)}^{-1}{A}_s{\zeta }^3-2(1-e^2{)}^{-1}{A}_{ss}{a}_{\ell }^2{\zeta }^5.$

This should match the more general "vertical pressure gradient" expression when we set, $\chi =0$, that is,

$\left\{\right[\frac{1}{(\pi G{\rho }_c^2{a}_{\ell }^2)}]\frac{\partial P}{\partial \zeta }{\}}_{\chi =0}$	$=$	$[1-{\xcancel{{\chi }^2}}^0-{\zeta }^2(1-e^2{)}^{-1}]\cdot [2A_{\ell s}{a}_{\ell }^2\zeta {\xcancel{{\chi }^2}}^0-2A_s\zeta +2A_{ss}{a}_{\ell }^2{\zeta }^3]$
	$=$	$[-2A_s\zeta +2A_{ss}{a}_{\ell }^2{\zeta }^3]+{\zeta }^2(1-e^2{)}^{-1}[2A_s\zeta -2A_{ss}{a}_{\ell }^2{\zeta }^3]$

Yes! The expressions match!

Test04

From above, we understand that, analytically,

$\left[\frac{1}{(\pi G{\rho }_c^2{a}_{\ell }^2)}\right]\frac{\partial P}{\partial \zeta }$	$=$	$[1-{\chi }^2-{\zeta }^2(1-e^2{)}^{-1}\left]\right[2A_{\ell s}{a}_{\ell }^2{\chi }^2\zeta -2A_s\zeta +2A_{ss}{a}_{\ell }^2{\zeta }^3]$
	$=$	$[(2A_{\ell s}{a}_{\ell }^2{\chi }^2-2A_s)-(2A_{\ell s}{a}_{\ell }^2{\chi }^4-2A_s{\chi }^2)]\zeta +[2A_{ss}{a}_{\ell }^2-2A_{ss}{a}_{\ell }^2{\chi }^2-(1-e^2{)}^{-1}(2A_{\ell s}{a}_{\ell }^2{\chi }^2-2A_s)]{\zeta }^3+[-(1-e^2{)}^{-1}2A_{ss}{a}_{\ell }^2]{\zeta }^5$
	$=$	$[2A_s({\chi }^2-1)+2A_{\ell s}{a}_{\ell }^2(1-{\chi }^2){\chi }^2]\zeta +[2A_{ss}{a}_{\ell }^2(1-{\chi }^2)-2A_{\ell s}{a}_{\ell }^2(1-e^2{)}^{-1}{\chi }^2+2(1-e^2{)}^{-1}{A}_s]{\zeta }^3+[-2A_{ss}{a}_{\ell }^2(1-e^2{)}^{-1}]{\zeta }^5.$

Also from above, we have shown that if,

$\frac{P_{\mathrm{t}\mathrm{e}\mathrm{s}\mathrm{t}\mathrm{0}\mathrm{2}}}{P_c}$

$=$

$p_2(\frac{\rho }{{\rho }_c}{)}^2+p_3(\frac{\rho }{{\rho }_c}{)}^3$

SUMMARY from test02: $p_2$ $=$ $3-e^2-(1+e^2)\mathrm{{\rm Y}}=e^4(A_{\ell s}{a}_{\ell }^2)-e^2\mathrm{{\rm Y}},$ $p_3$ $=$ $\frac{2}{3}[(4e^2-3)+\mathrm{{\rm Y}}]=e^4(A_{ss}{a}_{\ell }^2)+\frac{2}{3}{e}^2\mathrm{{\rm Y}}.$

$\Rightarrow \frac{\partial }{\partial \zeta }\left[\frac{P_{\mathrm{t}\mathrm{e}\mathrm{s}\mathrm{t}\mathrm{0}\mathrm{2}}}{P_c}\right]$	$=$	$\left(\frac{\rho }{{\rho }_c}\right)(1-e^2{)}^{-2}\{6p_3{\chi }^2\zeta (1-e^2)-2(2p_2+3p_3)(1-e^2)\zeta +6p_3{\zeta }^3\}$
	$=$	$\left(\frac{\rho }{{\rho }_c}\right)(1-e^2{)}^{-2}\{6[e^4(A_{ss}{a}_{\ell }^2)+\frac{2}{3}{e}^2\mathrm{{\rm Y}}]{\chi }^2\zeta (1-e^2)-2[2e^4(A_{\ell s}{a}_{\ell }^2)+3e^4(A_{ss}{a}_{\ell }^2)](1-e^2)\zeta +6[e^4(A_{ss}{a}_{\ell }^2)+\frac{2}{3}{e}^2\mathrm{{\rm Y}}]{\zeta }^3\}$

---

Here (test04), we add a term that is linear in the normalized density, which means,

$\frac{P_{\mathrm{t}\mathrm{e}\mathrm{s}\mathrm{t}\mathrm{0}\mathrm{4}}}{P_c}$	$=$	$\frac{P_{\mathrm{t}\mathrm{e}\mathrm{s}\mathrm{t}\mathrm{0}\mathrm{2}}}{P_c}+p_1\left(\frac{\rho }{{\rho }_c}\right)$
$\Rightarrow \frac{\partial }{\partial \zeta }\left[\frac{P_{\mathrm{t}\mathrm{e}\mathrm{s}\mathrm{t}\mathrm{0}\mathrm{4}}}{P_c}\right]$	$=$	$\frac{\partial }{\partial \zeta }\left[\frac{P_{\mathrm{t}\mathrm{e}\mathrm{s}\mathrm{t}\mathrm{0}\mathrm{2}}}{P_c}\right]+\frac{\partial }{\partial \zeta }\left[p_1\right(\frac{\rho }{{\rho }_c}\left)\right]=\frac{\partial }{\partial \zeta }\left[\frac{P_{\mathrm{t}\mathrm{e}\mathrm{s}\mathrm{t}\mathrm{0}\mathrm{2}}}{P_c}\right]+p_1\frac{\partial }{\partial \zeta }[1-{\chi }^2-{\zeta }^2(1-e^2{)}^{-1}]$

See Also