Siacci's theorem

In kinematics, the acceleration of a particle moving along a curve in space is the time derivative of its velocity. In most applications, the acceleration vector is expressed as the sum of its normal and tangential components, which are orthogonal to each other. Siacci's theorem, formulated by the Italian mathematician Francesco Siacci (1839–1907), is the kinematical decomposition of the acceleration vector into its radial and tangential components. In general, the radial and tangential components are not orthogonal to each other. Siacci's theorem is particularly useful in motions where the angular momentum is constant.

Siacci's theorem in the plane

Motion of a particle P in a plane.

Let a particle P of mass m  move in a two-dimensional Euclidean space (planar motion). Suppose that C is the curve traced out by P and s is the arc length of C corresponding to time t. Let O be an arbitrary origin in the plane and {i,j} be a fixed orthonormal basis. The position vector of the particle is

$\mathbf {r} =r\mathbf {e} _{r}.$

The unit vector e_r is the radial basis vector of a polar coordinate system in the plane. The velocity vector of the particle is

$\mathbf {v} ={\frac {d\mathbf {r} }{dt}}={\dot {s}}\mathbf {e} _{t}=v\mathbf {e} _{t},$

where e_t is the unit tangent vector to C. Define the angular momentum of P as

$\mathbf {h} =\mathbf {r} \times m\mathbf {v} =h\mathbf {k} ,$

where k = i x j. Assume that h ≠ 0. The position vector r may then be expressed as

$\mathbf {r} =q\mathbf {e} _{t}-p\mathbf {e} _{n}$

in the Serret-Frenet Basis {e_t, e_n, e_b}. The magnitude of the angular momentum is h = mpv, where p is the perpendicular from the origin to the tangent line ZP. According to Siacci's theorem, the acceleration a of P can be expressed as

$\mathbf {a} =-{\frac {\kappa v^{2}r}{p}}\mathbf {e} _{r}+{\frac {(h^{2})'}{2p^{2}}}\mathbf {e} _{t}=S_{r}\mathbf {e} _{r}+S_{t}\mathbf {e} _{t}.$

where the prime denotes differentiation with respect to the arc length s, and κ is the curvature function of the curve C. In general, S_r and S_t are not equal to the orthogonal projections of a onto e_r and e_t.

Example: Central forces

Suppose that the angular momentum of the particle P is a nonzero constant and that S_r is a function of r. Then

$S_{r}=-f(r)=-{\frac {\kappa rh^{2}}{p^{3}}},\qquad S_{t}=0.$

Because the curvature at a point in an orbit is given by

$\kappa ={\frac {1}{r}}{\frac {dp}{dr}},$

the function f can be conveniently written as a first order ODE

$f(r)={\frac {h^{2}}{p^{3}}}{\frac {dp}{dr}}.$

The energy conservation equation for the particle is then obtained if f(r) is integrable.

${\frac {1}{2}}{\frac {h^{2}}{p^{2}}}+\int f(r)=\mathrm {constant} .$

Siacci's theorem in space

Siacci's theorem can be extended to three-dimensional motions. Thus, let C be a space curve traced out by P and s is the arc length of C corresponding to time t. Also, suppose that the binormal component of the angular momentum does not vanish. Then the acceleration vector of P can be expressed as

$\mathbf {a} =-{\frac {\kappa v^{2}r}{p}}\mathbf {e} _{r}+\left(v{\frac {dv}{ds}}+{\frac {\kappa v^{2}q}{p}}\right)\mathbf {e} _{t}.$

The tangential component is tangent to the curve C. The radial component is directed from the point P to the point where the perpendicular from an arbitrary fixed origin meets the osculating plane. Other expressions for a can be found in, where a new proof of Siacci's theorem is given.

See also

• Acceleration
• Areal velocity
• Central force
• Serret-Frenet equations

References

• F. Siacci. Moto per una linea plana. Atti della Reale Accademia della Scienze di Torino, XIV, 750–760, 1879.
• F. Siacci. Moto per una linea gobba. Atti della Reale Accademia della Scienze di Torino, XIV, 946–951, 1879.
• E. T. Whittaker. A Treatise on the Analytical Dynamics of Particles and Rigid Bodies. 4th edition, Cambridge University Press, Cambridge. Reprinted by Dover Publications, Inc., New York (1944).
• Nathaniel Grossman. The sheer joy of celestial mechanics. Birkhäuser, Basel, 1996.