Next: Euler Method Up: Kinematics and Dynamics of Previous: Hamilton-Jacobi Equation

## Solving the Hamilton-Jacobi Equation

In this section we're going to use the Hamilton-Jacobi equation in order to derive analytical formulae for a motion of a material point in the central Newtonian M/r potential.

The derived formulae can then be used to compare an approximate numerical solution against an analytical, i.e., exact solution to the problem.

The problem is easiest to describe in spherical coordinates, r, and . In these coordinates the Hamiltonian assumes the following form:

 (4.41)

where we have assumed that the gravitational constant G = 1, and the mass of the material point m = 1 too.

The Hamilton-Jacobi equation that corresponds to that Hamiltonian is

 (4.42)

The method that is commonly use here is called the separation of variables. The Hamiltonian does not depend explicitly on
• time, t
• angle
Therefore a solution is sought in the following form:

 (4.43)

where C is a constant.

Substituting this into the Hamilton-Jacobi equation yields:

 (4.44)

Multiply this equation by 2 r2. Now we can rewrite this equation placing all terms that depend on on the left hand side and all terms that depend on r on the right hand side:

 (4.45)

Because expressions on both sides of this equation depend on different variables, the equality can hold only if they are equal to the same constant, L2:

 (4.46)

And this implies that:
 = (4.47) = (4.48)

These, in turn, are first order ordinary differential equations, which can be readily integrated.

But recall that there are additional conditions that function S must satisfy. In this case

 (4.49)

where Pr, and are constants and
 = Qr (4.50) = (4.51) = (4.52)

where Qr, and are also constants. In our case we have , , and , and we can set Qr, and to zero so that:

 = 0 (4.53) = 0 (4.54) = 0 (4.55)

The last equation (4.59) yields:

 (4.56)

It can be proven that this relation is satisfied by a flat motion, i.e., that the material point moves in a plane with vector L perpendicular to that plane. We can therefore change our system of coordinates so that, say, . Then

 (4.57)

Equation (4.58) yields the shape of the orbit:

 (4.58)

with the following solution:

 (4.59)

where e is the eccentricity of the orbit:

 (4.60)

Other parameters pertaining to the orbit are:

 (4.61)

Remember that for a trapped particle the energy is negative. a is the semimajor axis of the orbit (when elliptic).

 (4.62)

b is the semiminor axis of the orbit (when elliptic)

 (4.63)

is the distance of the closest approach.

Equation (4.57) yields the time dependence of r versus t:

 (4.64)

The solution to this equation is quite complicated and can be given in terms of Bessel functions and harmonic motion with the mean circular frequency of

 (4.65)

Next: Euler Method Up: Kinematics and Dynamics of Previous: Hamilton-Jacobi Equation
Zdzislaw Meglicki
2001-02-26