Appendix/Ramblings/ForPaulFisher

From JETohlineWiki
Jump to navigation Jump to search


For Paul Fisher[edit]

Overview of Dissertation[edit]

Paul Fisher's (1999) doctoral dissertation (accessible via the LSU Digital Commons) is titled, Nonaxisymmetric Equilibrium Models for Gaseous Galaxy Disks. Its abstract reads, in part:

Three-dimensional hydrodynamic simulations show that, in the absence of self-gravity, an axisymmetric, gaseous galaxy disk whose angular momentum vector is initially tipped at an angle, i0, to the symmetry axis of a fixed spheroidal dark matter halo potential does not settle to the equatorial plane of the halo. Instead, the disk settles to a plane that is tipped at an angle, α=tan−1[q2tan⁡i0], to the equatorial plane of the halo, where q is the axis ratio of the halo equipotential surfaces. The equilibrium configuration to which the disk settles appears to be flat but it exhibits distinct nonaxisymmetric features. .

All three-dimensional hydrodynamic simulations employ Richstone's (1980) time-independent "axisymmetric logarithmic potential" that is prescribed by the expression,

Φ(x,y,z)

=

v022ln⁡[x2+y2+z2q2].

Thoughts Moving Forward[edit]

Let's continue to examine a collection of Lagrangian fluid elements that are orbiting in an (axisymmetric) oblate-spheroidal potential with flattening "q." But rather than adopting the Richstone potential, we will consider the potential generated inside an homogeneous (i.e., Maclaurin) spheroid whose eccentricity is, e=(1−q2)1/2, namely,

Φ(ϖ,z)=−πGρ[IBTa12−(A1ϖ2+A3z2)],

[ST83], §7.3, p. 169, Eq. (7.3.1)

where, the coefficients A1, A3, and IBT are functions only of the spheroid's eccentricity. What does the potential field look like from the perspective of a particle/fluid-element whose orbital angular momentum vector is tipped at an angle, i0, to the symmetry axis of the oblate-spheroidal potential? Presumably the potential is "observed" to vary with position around the orbit as though the underlying potential is non-axisymmetric. Does it appear to be the potential inside a Riemann S-Type ellipsoid? If so, what values of (b/a,c/a) correspond to the chosen parameter pair, (q,i0)?

Well, let's define a primed (Cartesian) coordinate system whose z'-axis is tipped at this angle, i0, with respect to the symmetry axis of the oblate-spheroidal potential. Drawing from a discussion in which we have presented a closely analogous methodical derivation of orbital parameters, we have,

x′

=

x,

y′

=

ycos⁡i0+(z−z0)sin⁡i0,

z′

=

(z−z0)cos⁡i0−ysin⁡i0.

x

=

x′,

y

=

y′cos⁡i0−z′sin⁡i0,

z−z0

=

z′cos⁡i0+y′sin⁡i0.

x˙′

=

x˙,

y˙′

=

y˙cos⁡i0+z˙sin⁡i0,

z˙′0

=

z˙cos⁡i0−y˙sin⁡i0.

x˙

=

x˙′,

y˙

=

y˙′cos⁡i0−z˙′0sin⁡i0,

z˙

=

z˙′0cos⁡i0+y˙′sin⁡i0.

When viewed from this primed frame, the potential associated with a Maclaurin spheroid becomes,

(πGρ)−1Φ(x′,y′,z′)+IBTa12

=

A1[(x′)2+(y′cos⁡i0−z′sin⁡i0)2]+A3[z0+z′cos⁡i0+y′sin⁡i0]2

 

=

A1[(x′)2+(y′)2cos2i0+(z′)2sin2i0−2(y′z′)sin⁡i0cos⁡i0]

 

 

+A3[z02+2z′z0cos⁡i0+2z0y′sin⁡i0+(z′)2cos⁡2i0+2y′z′sin⁡i0cos⁡i0+(y′)2sin⁡2i0]

 

=

A1(x′)2+(y′)2[A1cos2i0+A3sin2i0]+(z′)2[A1sin2i0+A3cos2i0]

 

 

+z0A3[z0+2z′cos⁡i0+2y′sin⁡i0]+2(A3−A1)y′z′sin⁡i0cos⁡i0.

See Also[edit]

Tiled Menu

Appendices: | VisTrailsEquations | VisTrailsVariables | References | Ramblings | VisTrailsImages | myphys.lsu | ADS |