What's New
About the Project
NIST
22 Jacobian Elliptic FunctionsApplications

§22.19 Physical Applications

Contents

§22.19(i) Classical Dynamics: The Pendulum

With appropriate scalings, Newton’s equation of motion for a pendulum with a mass in a gravitational field constrained to move in a vertical plane at a fixed distance from a fulcrum is

22.19.1 d2θ(t)dt2=-sinθ(t),

θ being the angular displacement from the point of stable equilibrium, θ=0. The bounded (-πθπ) oscillatory solution of (22.19.1) is traditionally written

22.19.2 sin(12θ(t))=sin(12α)sn(t+K,sin(12α)),

for an initial angular displacement α, with dθ/dt=0 at time 0; see Lawden (1989, pp. 114–117). The period is 4K(sin(12α)). The periodicity and symmetry of the pendulum imply that the motion in each four intervals θ(0,±α) and θ(±α,0) have the same “quarter periods” K=K(sin(12α)). Thus the offset t+K in 22.19.2 as the motion starts θ(0)=α, rather than θ(0)=0 as in 22.19.3, which follows. The angle α=π is a separatrix, separating oscillatory and unbounded motion. With the same initial conditions, if the sign of gravity is reversed then the new period is 4K(sin(12α)); see Whittaker (1964, §44).

Alternatively, Sala (1989) writes:

22.19.3 θ(t)=2am(tE/2,2/E),

for the initial conditions θ(0)=0, the point of stable equilibrium for E=0, and dθ(t)/dt=2E. Here E=12(dθ(t)/dt)2+1-cosθ(t) is the energy, which is a first integral of the motion. This formulation gives the bounded and unbounded solutions from the same formula (22.19.3), for k1 and k1, respectively. Also, θ(t) is not restricted to the principal range -πθπ. Figure 22.19.1 shows the nature of the solutions θ(t) of (22.19.3) by graphing am(x,k) for both 0k1, as in Figure 22.16.1, and k1, where it is periodic.

See accompanying text
Figure 22.19.1: Jacobi’s amplitude function am(x,k) for 0x10π and k=0.5,0.9999,1.0001,2. When k<1, am(x,k) increases monotonically indicating that the motion of the pendulum is unbounded in θ, corresponding to free rotation about the fulcrum; compare Figure 22.16.1. As k1-, plateaus are seen as the motion approaches the separatrix where θ=nπ, n=±1,±2,, at which points the motion is time independent for k=1. This corresponds to the pendulum being “upside down” at a point of unstable equilibrium. For k>1, the motion is periodic in x, corresponding to bounded oscillatory motion. Magnify

§22.19(ii) Classical Dynamics: The Quartic Oscillator

Classical motion in one dimension is described by Newton’s equation

22.19.4 d2x(t)dt2=-dV(x)dx,

where V(x) is the potential energy, and x(t) is the coordinate as a function of time t. The potential

22.19.5 V(x)=±12x2±14βx4

plays a prototypal role in classical mechanics (Lawden (1989, §5.2)), quantum mechanics (Schulman (1981, Chapter 29)), and quantum field theory (Pokorski (1987, p. 203), Parisi (1988, §14.6)). Its dynamics for purely imaginary time is connected to the theory of instantons (Itzykson and Zuber (1980, p. 572), Schäfer and Shuryak (1998)), to WKB theory, and to large-order perturbation theory (Bender and Wu (1973), Simon (1982)).

For β real and positive, three of the four possible combinations of signs give rise to bounded oscillatory motions. We consider the case of a particle of mass 1, initially held at rest at displacement a from the origin and then released at time t=0. The subsequent position as a function of time, x(t), for the three cases is given with results expressed in terms of a and the dimensionless parameter η=12βa2.

Case I: V(x)=12x2+14βx4

This is an example of Duffing’s equation; see Ablowitz and Clarkson (1991, pp. 150–152) and Lawden (1989, pp. 117–119). The subsequent time evolution is always oscillatory with period 4K(k)/1+2η and modulus k=1/2+η-1:

22.19.6 x(t)=acn(t1+2η,k).

Case II: V(x)=12x2-14βx4

There is bounded oscillatory motion near x=0, with period 4K(k)/1-η, and modulus k=1/η-1-1, for initial displacements with |a|1/β.

22.19.7 x(t)=asn(t1-η,k).

As a1/β from below the period diverges since a=±1/β are points of unstable equilibrium.

Case III: V(x)=-12x2+14βx4

Two types of oscillatory motion are possible. For an initial displacement with 1/β|a|<2/β, bounded oscillations take place near one of the two points of stable equilibrium x=±1/β. Such oscillations, of period 2K(k)/η, with modulus k=1/2-η-1 are given by:

22.19.8 x(t)=adn(tη,k).

As a2/β from below the period diverges since x=0 is a point of unstable equlilibrium. For initial displacement with |a|2/β the motion extends over the full range -axa:

22.19.9 x(t)=acn(t2η-1,k),

with period 4K(k)/2η-1 and modulus k=1/2-η-1. As |a|1/β from above the period again diverges. Both the dn and cn solutions approach asecht as a2/β from the appropriate directions.

§22.19(iii) Nonlinear ODEs and PDEs

Many nonlinear ordinary and partial differential equations have solutions that may be expressed in terms of Jacobian elliptic functions. These include the time dependent, and time independent, nonlinear Schrödinger equations (NLSE) (Drazin and Johnson (1993, Chapter 2), Ablowitz and Clarkson (1991, pp. 42, 99)), the Korteweg–de Vries (KdV) equation (Kruskal (1974), Li and Olver (2000)), the sine-Gordon equation, and others; see Drazin and Johnson (1993, Chapter 2) for an overview. Such solutions include standing or stationary waves, periodic cnoidal waves, and single and multi-solitons occurring in diverse physical situations such as water waves, optical pulses, quantum fluids, and electrical impulses (Hasegawa (1989), Carr et al. (2000), Kivshar and Luther-Davies (1998), and Boyd (1998, Appendix D2.2)).

§22.19(iv) Tops

The classical rotation of rigid bodies in free space or about a fixed point may be described in terms of elliptic, or hyperelliptic, functions if the motion is integrable (Audin (1999, Chapter 1)). Hyperelliptic functions u(z) are solutions of the equation z=0u(f(x))-1/2dx, where f(x) is a polynomial of degree higher than 4. Elementary discussions of this topic appear in Lawden (1989, §5.7), Greenhill (1959, pp. 101–103), and Whittaker (1964, Chapter VI). A more abstract overview is Audin (1999, Chapters III and IV), and a complete discussion of analytical solutions in the elliptic and hyperelliptic cases appears in Golubev (1960, Chapters V and VII), the original hyperelliptic investigation being due to Kowalevski (1889).

§22.19(v) Other Applications

Numerous other physical or engineering applications involving Jacobian elliptic functions, and their inverses, to problems of classical dynamics, electrostatics, and hydrodynamics appear in Bowman (1953, Chapters VII and VIII) and Lawden (1989, Chapter 5). Whittaker (1964, Chapter IV) enumerates the complete class of one-body classical mechanical problems that are solvable this way.