Variable mass problems are the archetypal problems of `rocket science’ in which aspects of the motion of a rocket need to be calculated as the rocket burns up its fuel. Variable mass problems also arise in situations in which a body gains mass by accretion, e.g., the way a hailstone (pictured) grows by condensation as it falls. The calculations are based on the impulse-momentum principle which states that the change in linear momentum of a system in the time interval
is equal to the impulse of the external forces acting on the system in that interval. This principle derives from Newton’s second law for systems of particles, expressed in terms of momentum (intuitively,
gives us
).
In this note I want to explore a variable mass problem involving a falling hailstone gaining mass by accretion, focusing on the derivation and solution of the differential equation of motion. I find the problem interesting because it is illustrative of the fact that it is sometimes possible to solve for the required quantities by direct integration without going through the step of first finding a general solution of the differential equation. However, going through the longer process of first finding a general solution of the differential equation enables further insights to be obtained. In this note I will solve the problem in both ways, i.e., by direct integration and by finding a general solution of the differential equation. I will use the latter to delve a bit deeper into the situation.
The problem is as follows (the reader should attempt to solve it before reading on):
A spherical hailstone is falling under gravity in still air. At time the hailstone has speed
. The radius
increases by condensation at the rate
where
is a constant. Ignoring air resistance, derive the equation of motion of the hailstone and find the time taken for the speed of the hailstone to increase from
to
.
It helps to draw a sketch of the situation, such as the following:
Using the impulse-momentum principle () and noting that the external force here is equal to mass times acceleration due to gravity, we can write
which (treating products of differentials as zero) simplifies to
Note that the sign of the acceleration due to gravity on the right hand side is positive because it is in the same direction as the velocity of the falling hailstone. Dividing through by gives
At this point we can use the fact that mass is proportional to volume in this case, so
and therefore
so
which simplifies to
This is the required equation of motion for the hailstone. As stated earlier, we can now proceed to find the time taken for the speed of the hailstone to increase from to
either by direct integration, or by first solving the differential equation of motion to find an expression for
. The latter approach has the advantage that we can use the expression for
to obtain other insights. Using the direct integration approach, we would separate the variables and write
and therefore the time taken for the speed of the hailstone to increase from to
is
To get the same answer by first solving the differential equation to obtain an expression for , we would observe that the general solution of a first-order differential equation of the form
is
where
Here we have the differential equation
so
and therefore
So the general solution is
where is an arbitrary constant.
Now, when , we have
and solving for gives
Similarly, when , we have
and solving for gives
Then the time taken is
as before. Note that the arbitrary constant disappears in the course of the calculation.
Although this latter approach seems more long-winded, it has the advantage that we can now use the general solution of the differential equation we have obtained to explore other features of the problem. For example, although it is not specifically required in the original problem, we could work out how far the hailstone travelled in the time it took for its speed to increase from to
. To do this, we would rewrite the general solution of the differential equation obtained above as
Separating the variables and integrating on both sides then gives the following expression for the position of the hailstone as a function of time:
where is a second arbitrary constant. At the time
calculated above we have
and at the time calculated above we have
Therefore the distance travelled by the hailstone in the time it took for its speed to increase from to
is
Note that both the arbitrary constants and
have disappeared in the course of the calculation.