*A Tale of two cheeses - *Adrian Oldknow a_oldknow@compuserve.com

*Chapter 5: Limits to growth: the falling Cheddar
*Chapter 1 used a model of acceleration which might be appropriate in a vacuum. However Cheddar is not very aerodynamic (drag coefficient = ??), and the nice warm summer's air would resist the cheese's motion, giving a deceleration which reduces the acceleration due to gravity. This resistance increases with the speed of the cheese and acts in the opposite direction to its velocity. At some critical speed the deceleration due to air resistance will cancel out the acceleration due to gravitational attraction entirely and the Cheddar will then have a constant velocity (called the terminal velocity). The refined model for the acceleration

Assuming that this acceleration stays roughly constant in each small time interval you can now make up a modified table of velocity

How long, approximately, will it take for *v* to reach 50 m/s? How much does this change if you compute the velocities every 0.1 s instead of every 1 s? (You may need a computer!) What does the graph of *v* against *t* look like now? Can you extend the table to include the distance fallen *d? *What should/does the graph of *d *against *t* look like? How far will the Cheddar have fallen when its speed reaches 50 m/s? A recent newspaper article stated that a Cheddar cheese will reach 99% of its terminal velocity of 53 m/s after about 14 seconds in which time it will have fallen about 570 m. See if you can find a value for *p* (and the corresponding value of *k*) for this data. If the balloon was at a height of 4000 m above the ground when the accidents occurred can you estimate the number of Listeria bacteria in the Camembert at the time the Cheddar hit the ground, and the speed at which the Cheddar was going?

*----------------------------------------------------------------------------------------------------------------------*

*
*The differential_{ }equation is now _{ }with *v* = 0 when *t* = 0. The difference equation_{ }is again non-linear (unless *p* = 1) and there is no analytic solution. The acceleration is zero when and the corresponding value of *v* is called the terminal velocity.

We'll start with the numerical solution:

The analytic solution is given by:

which we'll derive soon.

Again we just quote the analytic solution - read on!

Analytically we need to evaluate the definite integral:

With *p* = 2 this becomes, writing *tv* for terminal velocity:

With *p* = 1 this is:

And with *p* = 3 it is:

Ouch!!