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, α=tan1[q2tani0], 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=(1q2)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

=

ycosi0+(zz0)sini0,

z

=

(zz0)cosi0ysini0.

x

=

x,

y

=

ycosi0zsini0,

zz0

=

zcosi0+ysini0.

x˙

=

x˙,

y˙

=

y˙cosi0+z˙sini0,

z˙0

=

z˙cosi0y˙sini0.

x˙

=

x˙,

y˙

=

y˙cosi0z˙0sini0,

z˙

=

z˙0cosi0+y˙sini0.

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

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

=

A1[(x)2+(ycosi0zsini0)2]+A3[z0+zcosi0+ysini0]2

 

=

A1[(x)2+(y)2cos2i0+(z)2sin2i02(yz)sini0cosi0]

 

 

+A3[z02+2zz0cosi0+2z0ysini0+(z)2cos2i0+2yzsini0cosi0+(y)2sin2i0]

 

=

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

 

 

+z0A3[z0+2zcosi0+2ysini0]+2(A3A1)yzsini0cosi0.

See Also[edit]

Tiled Menu

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