The Structure, Stability, and Dynamics
of Self-Gravitating Systems

Joel E. Tohline
tohline@rouge.phys.lsu.edu

[ Download PDF file dated ]

Integrals of Motion in Time-Invariant Potentials

The time-evolution of dynamical systems is usually much easier to understand, or at least the future behavior of such systems is much easier to predict, once we realize that certain constants or "integrals" of the motion are preserved in time. The most familiar integrals of motion are energy and angular momentum. It is important to realize, however, that when particles (or fluid elements) move under the influence of an external potential, these familiar integrals may not be constant in time, but other integrals of the motion may exist. We present here a general discussion of integrals of motion for systems of noninteracting particles (i.e., effects of the gas pressure will be ignored) that move in time-invarient, external potentials.


Referring back to the principal governing equations, if the effects of pressure are ignored Euler's equation takes the form

Dv = - ÑF.

[Equation VI.I.1]

Our objective is to find one or more "integrals" of the motion Ii for which the mathematical statement,

D(Ii) = 0,

[Equation VI.I.2]

is both true and derivable from the equation of motion.


We'll focus, first, on the right-hand-side of the equation of motion. From the general (curvilinear form) of the gradient operator, we can write

- ÑF = - e1 ( 1/h1 ) x1F - e2 ( 1/h2 ) x2F - e3 ( 1/h3 ) x3F .

[Equation VI.I.3]


Note, also, that from the definition of the Lagrangian time derivative of any scalar function, in general we can write,

DF = tF + v×ÑF.

[Equation VI.I.4]

Here, physically, the Lagrangian (total) time-derivative DF identifies the time-variation of the gravitational potential as seen by a moving particle (fluid element) while the Eulerian (partial) time-derivative identifies the actual time-variation of the potential as viewed by a stationary observer. Our discussions will be confined to time-invarient potentials, as desired, by simply demanding that

tF = 0.

[Equation VI.I.5]

Hence, for our present discussions we may adopt the following relationship between the time- and spatial-variation of the gravitational potential, as seen by a moving particle:

DF = v×ÑF.

[Equation VI.I.6]


Now we turn our attention to the left-hand-side of the equation of motion. Realizing that the general (curvilinear coordinate) form of the velocity is

v = e1 (h1Dx1) + e2 (h2Dx2) + e3 (h3Dx3),

[Equation VI.M.18]

we can, quite generally, write the time-derivative of v as,

Dv = e1 { d(h1Dx1)/dt - (h2Dx2) A - (h3Dx3) B } +
e2 { d(h2Dx2)/dt + (h1Dx1) A - (h3Dx3) C } +
e3 { d(h3Dx3)/dt + (h1Dx1) B + (h2Dx2) C }.

[Equation VI.I.7]

where, as before [VI.M.17]:
A º f1Dx2 - f3Dx1 = [Dx2] (1/h1) x1h2 - [Dx1] (1/h2) x2h1
B º f2Dx3 - f5Dx1 = [Dx3] (1/h1) x1h3 - [Dx1] (1/h3) x3h1
C º f4Dx3 - f6Dx2 = [Dx3] (1/h2) x2h3 - [Dx2] (1/h3) x3h2

Utilizing the general expression for the time-derivative of any unit vector, derive the above expression [VI.I.7] for the acceleration (Dv) from the definition of v.


As a shorthand, it frequently is useful to adopt the variables,

v1 º h1Dx1,
v2 º h2Dx2,
v3 º h3Dx3,

[Equation VI.I.8]

to denote the three components of the particle velocity. In terms of these variables, the above expressions for the velocity vector [VI.M.18] and acceleration vector [VI.I.7] take the following, simpler forms:

v = e1 v1 + e2 v2 + e3 v3,

[Equation VI.I.9]

a º Dv = e1 { dv1 /dt - v2A - v3B } +
e2 { dv2 /dt + v1A - v3C } +
e3 { dv3 /dt + v1B + v2C }.

[Equation VI.I.10]

Combining the above expressions for the left-hand-side [VI.I.10] and the right-hand-side [VI.I.3] of the equation of motion, we can write in the following very general form the

Three Components of the Equation of Motion
[Equation VI.I.11]
e1 dv1 /dt - v2A - v3B = - ( 1/h1 ) x1F
e2 dv2 /dt + v1A - v3C = - ( 1/h2 ) x2F
e3 dv3 /dt + v1B + v2C = - ( 1/h3 ) x3F


Building upon this presentation of the equation of motion written in curvilinear coordinates, the accompanying pages present derivations and discussions of:


Other Integrals


Situations in which F = F( x1 )


Home Page | Preface | Context | Applications | Appendices | Search Index