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.