# 8.7: Variable-mass systems

- Page ID
- 14101

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\id}{\mathrm{id}}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\kernel}{\mathrm{null}\,}\)

\( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\)

\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\)

\( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)

\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)

\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vectorC}[1]{\textbf{#1}} \)

\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)Lagrangian and Hamiltonian mechanics assume that the total mass and energy of the system are conserved. Variable-mass systems involve transferring mass and energy between donor and receptor bodies. However, such systems still can be conservative if the Lagrangian or Hamiltonian include all the active degrees of freedom for the combined donor-receptor system. The following examples of variable mass systems illustrate subtle complications that occur handling such problems using algebraic mechanics.

## Rocket propulsion:

Newtonian mechanics was used to solve the rocket problem in chapter \(2.12\). The equation of motion (\(2.12.23\)) relating the rocket thrust \(F_{ex}\) to the rate of change of the momentum separated into two terms,

\[F_{ex}=\dot{p}_{y}=m\ddot{y}+\dot{m}\dot{y}\]

The first term is the usual mass times acceleration, while the second term arises from the rate of change of mass times the velocity. The equation of motion for rocket motion is easily derived using either Lagrangian or Hamiltonian mechanics by relating the rocket thrust to the generalized force \(Q_{j}^{EXC}.\)

## Moving chains:

The motion of a flexible, frictionless, heavy chain that is falling in a gravitational field, often can be split into two coupled variable-mass partitions that have different chain-link velocities. These partitions are coupled at the moving intersection between the chain partitions. That is, these partitions share time-dependent fractions of the total chain mass. Moving chains were discussed first by Caley in \(1857\) and since then the moving chain problem has had a controversial history due to the frequent erroneous assumption that, in the gravitational field, the chain partitions fall with acceleration \(g\) rather than applying the correct energy conservation assumption for this conservative system. The following two examples of conservative falling-chain systems illustrate solutions obtained using variational principles applied to a single chain that is partitioned into two variable length sections.^{1}

Consider the following two possible scenarios for motion of a flexible, heavy, frictionless, chain located in a uniform gravitational field \(g\). The first scenario is the "folded chain" system which assumes that one end of the chain is held fixed, while the adjacent free end is released at the same altitude as the top of the fixed arm, and this free end is allowed to fall in the constant gravitational field \(g\). The second "falling chain", scenario assumes that one end of the chain is hanging down through a hole in a frictionless, smooth, rigid, horizontal table, with the stationary partition of the chain sitting on the table surrounding the hole. The falling section of this chain is being pulled out of the stationary pile by the hanging partition. Both of these systems are conservative since it is assumed that the total mass of the chain is fixed, and no dissipative forces are acting. The chains are assumed to be inextensible, flexible, and frictionless, and subject to a uniform gravitational field \(g\) in the vertical \(y\) direction. In both examples, the chain, with mass \(M\) and length \(L,\) is partitioned into a stationary segment, plus a moving segment, where the mass per unit length of the chain is \(\mu =\frac{M}{L}\). These partitions are strongly coupled at their intersection which propagates downward with time for the "folded chain" and propagates upward, relative to the lower end of the falling chain, for the "falling chain". For the "folded chain", the chain links are transferred from the moving segment to the stationary segment as the moving section falls. By contrast, for the "falling system", the chain links are transferred from the stationary upper section to the moving lower segment of the chain.

Example \(\PageIndex{1}\): Folded chain

The folded chain of length \(L\) and mass-per-unit-length \(\mu =\frac{M}{L}\) hangs vertically downwards in a gravitational field \(g\) with both ends held initially at the same height. The fixed end is attached to a fixed support while the free end of the chain is dropped at time \(t=0\) with the free end at the same height and adjacent to the fixed end. Let \(y\) be the distance the falling free end is below the fixed end. Using an idealized one-dimensional assumption, the Lagrangian \(\mathcal{L}\) is given by

\[\mathcal{L}(y,\dot{y})=\frac{M}{4L}(L-y)\dot{y}^{2}+Mg\frac{1}{4L} (L^{2}+2Ly-y^{2})\]

where the bracket in the second term is the height of the center of mass of the folded chain with respect to the fixed upper end of the chain.

The Hamiltonian is given by

\[H(y,p_{R})=p_{R}\dot{y}-\mathcal{L}(y,\dot{y})=\frac{p_{_{R}}}{\mu \left( L-y\right) }-Mg\frac{(L^{2}+2Ly-y^{2})}{4L}\]

where \(p_{R}\) is the linear momentum of the right-hand arm of the folded chain.

As shown in the discussion of the Generalized Energy Theorem, (chapters \(7.8\) and \(7.9\)), when all the active forces are included in the Lagrangian and the Hamiltonian, then the total mechanical energy \(E\) is given by \(E=H.\) Moreover, both the Lagrangian and the Hamiltonian are time independent, since

\[\frac{dE}{dt}=\frac{dH}{dt}=-\frac{\partial \mathcal{L}}{\partial t}=0\]

Therefore the "folded chain" Hamiltonian equals the total energy, which is a constant of motion. Energy conservation for this system can be used to give

\[\frac{\mu }{4}\left( L-y\right) \dot{y}^{2}-\frac{1}{4}\mu g(L^{2}+2Ly-y^{2})=-\frac{1}{4}\mu gL^{2}\] Solve for \(\dot{y}^{2}\) gives

\[\dot{y}^{2}=g\frac{(2Ly-y^{2})}{L-y} \label{8.74}\]

The acceleration of the falling arm, \(\ddot{y},\) is given by taking the time derivative of Equation \ref{8.74}

\[\ddot{y}=g+\frac{g\left( 2Ly-y^{2}\right) }{2\left( L-y\right) }\]

The rate of change in linear momentum for the moving right side of the chain, \(\dot{p}_{R}\), is given by

\[\dot{p}_{R}=m_{R}\ddot{y}+\dot{m}_{R}\dot{y}=m_{R}g+m_{R}g\frac{(2Ly-y^{2})}{ 2\left( L-y\right) } \label{8.76}\]

For this energy-conserving chain, the tension in the chain \(T_{0}\) at the fixed end of the chain is given by

\[T_{0}=\frac{\mu g}{2}\left( L+y\right) +\frac{1}{4}\mu \dot{y}^{2} \label{8.77}\]

Equations \ref{8.74} and \ref{8.76}, imply that the tension \(T_{o}\) diverges to infinity when \(y\rightarrow L\). Calkin and March measured the \(y\) dependence of the chain tension at the support for the folded chain and observed the predicted \(y\) dependence. The maximum tension was \(\simeq\) \(25Mg,\) which is consistent with that predicted using Equation \ref{8.77} after taking into account the finite size and mass of individual links in the chain. This result is very different from that obtained using the erroneous assumption that the right arm falls with the free-fall acceleration \(g\), which implies a maximum tension \(T_{0}=\) \(2Mg\). Thus the free-fall assumption disagrees with the experimental results, in addition to violating energy conservation and the tenets of Lagrangian and Hamiltonian mechanics. That is, the experimental result demonstrates unambiguously that the energy conservation predictions apply in contradiction with the erroneous free-fall assumption.

The unusual feature of variable mass problems, such as the folded chain problem, is that the rate of change of momentum in Equation \ref{8.76} includes two contributions to the force and rate of change of momentum, that is, it includes both the acceleration term \(m_{R}\ddot{y}\) plus the variable mass term \(\dot{m}_{R}\dot{y}\) that accounts for the transfer of matter at the intersection of the moving and stationary partitions of the chain. At the transition point of the chain, moving links are transferred from the moving section and are added to the stationary subsection. Since this moving section is falling downwards, and the stationary section is stationary, then the transferred momentum is in a downward direction corresponding to an increased effective downward force. Thus the measured acceleration of the moving arm actually is faster than \(g\). A related phenomenon is the loud cracking sound heard when cracking a whip.

Example \(\PageIndex{2}\): Falling chain

The "falling chain", scenario assumes that one end of the chain is hanging down through a hole in a frictionless, smooth, rigid, horizontal table, with the stationary partition of the chain lying on the frictionless table surrounding the hole. The falling section of this chain is being pulled out of the stationary pile by the hanging partition. The analysis for the problem of the falling chain behaves differently from the folded chain. For the "falling- chain" let \(y\) be the falling distance of the lower end of the chain measured with respect to the table top. The Lagrangian and Hamiltonian are given by \[\begin{aligned} \mathcal{L}(y,\dot{y}) &=&\frac{\mu }{2}y\dot{y}^{2}+\mu g\frac{y^{2}}{2} \\ p_{y} &=&\frac{\partial \mathcal{L}}{\partial \dot{y}}=\mu y\dot{y} \\ H &=&\frac{p_{y}^{2}}{2\mu y}-\frac{\mu gy^{2}}{2}=E\end{aligned}\]

The Lagrangian and Hamiltonian are not explicitly time dependent, and the Hamiltonian equals the initial total energy, \(E_{0}\). Thus energy conservation can be used to give that

\[E=\frac{1}{2}\mu y(\dot{y}^{2}-gy)=E_{0}\]

Lagrange’s equation of motion gives \[\dot{p}_{y}=m_{y}\ddot{y}+\dot{m}_{y}\dot{y}=m_{y}g+\frac{1}{2}\mu \dot{y} ^{2}=Mg-T_{0}\]

The important difference between the folded chain and falling chain is that the moving component of the falling chain is gaining mass with time rather than losing mass. Also the tension in the chain \(T_{0}\) reduces the acceleration of the falling chain making it less than the free-fall value \(g\). This is in contrast to that for the folded chain system where the acceleration exceeds \(g\).

The above discussion shows that Lagrangian and Hamiltonian can be applied to variable-mass systems if both the donor and receptor degrees of freedom are included to ensure that the total mass is conserved.

^{1}Discussions with Professor Frank Wolfs stimulated inclusion of these two examples of moving chains.