A more accurate 2-person cycling simulation than my "falling cyclist" is given by a "constant power" approximation, so the differential equation is now m*dv/dt = P/v - (1/2)*D*v^2 where the terminal velocity is v_T=(2*P/D)^(1/3). Semi-reasonable numbers are m=75kg, v_T=9m/s, P=216watts, D=0.5926. While the above differential equation appears not to be solvable in any reasonable closed form, we can approximate solutions of small deviations from v_T. When this is done, we find that these small deviations decay exponentially to 0 with a time constant t_C = m*v_T^2/(3*P) ~ 9.375 seconds. Associated with this time constant is a distance constant t_C*v_T ~ 84.375 meters. Also, an energy constant t_C*P ~ 2.025 kJ. -- Suppose the shaded cyclist has a drag D=0.432 instead of D=0.5926. Then the terminal velocity for the shaded cyclist will be v_T=10m/s instead of 9m/s. The time constant for the shaded cyclist is 11.574 secs and the distance constant for the shaded cyclist is 115.74 meters. We can now estimate a speed for the two cyclists "working together" as a mini-peloton. (Actually, they aren't "working together" at all, because their behavior emerges from the equations.) If the leading cyclist is going faster than his/her terminal velocity (9m/s), then he/she will slow down with an exponential decay towards 9m/s. If the following cyclist is going slower than his/her terminal velocity (10m/s), then he/she will speed up with an exponential decay towards 10m/s. The average velocity of the two cyclists should end up somewhere in the range of 9-10m/s. Since there are only 2 cyclists, and they are both identical, symmetry requires that both spend 50% of the time in the lead and 50% of the time being shaded. But this is equivalent to facing a drag constant which is an arithmetic _average_ of the two drag constants, hence D~0.5123. So we can now calculate the average terminal velocity as (2*216/0.5123)^(1/3)=9.448m/s. Thus, the velocity of this mini-peloton is _not_ the arithmetic average of the two terminal velocities, but slightly less than that. In this simple model, the "dwell time" of each cyclist in the lead gets shorter & shorter, so the cyclists exchange places faster & faster until they are riding "together" -- the distance has converged to zero. Notice that there is no "coordination" between the two cyclists -- they simply produce power at a constant unvarying wattage. At 02:00 PM 7/17/2010, Henry Baker wrote:
Thanks very much, Hans! This paper is quite germane, although it is trying to be an accurate model, while I'm merely trying to capture the essential behavior.
---
Consider a one-dimensional system where the cyclist "molecules" are falling vertically, and can pass one another without any collisions or interactions with one another. However, the cyclist molecule at the front bears the full air resistance, while the cyclist molecules further behind are shaded to some extent. According to the paper Hans referenced, the shading is quite short-range; it falls to zero after about 3 meters.
We now consider 1, 2, 3, etc. numbers of cyclist molecules.
Start with a single cyclist molecule falling through air.
According to Wikipedia, an object of the density of water falling through the air reaches a terminal velocity of approx. 90*sqrt(d) meters/second ~= 201*sqrt(d) mph, where d is the approx diameter of the object. Wikipedia claims the terminal velocity of a human body is approx 70 m/s ~= 157 mph, terminal velocity of a cat approx 40 m/s ~= 89 mph, the terminal velocity of a bird approx 20 m/s ~= 45 mph.
http://en.wikipedia.org/wiki/Drag_%28physics%29
http://en.wikipedia.org/wiki/Free_fall
Consider the differential equation, where D=rho*Cd*A, rho=density of air, Cd=drag coefficient, A=cross sectional area.
m * dv/dt = 1/2 * D * v^2 - m * g
Solving, we get v as a function of t:
v(t) = - v_T * tanh(g*t/v_T), where v_T = sqrt(2*m*g/D)
Integrating, we get y as a function of t:
y(t) = y_0 - v_T^2/g * log(cosh(g*t/v_T))
Consider the following thought experiment with two identical cyclist molecules. One cyclist molecule is already falling at its terminal velocity v_T, and as it passes another cyclist molecule at t=0, the other cyclist starts accelerating at g (assume no shading). While the second cyclist will never catch up, it will asymptote falling at v_T at distance ln(2)/g*v_T^2 behind. Furthermore, this distance can be covered at speed v_T in time ln2/g*v_T; this time is a kind of "time constant": the second cyclist reaches 60% of terminal velocity within this time constant, and at 4 times this time constant, the second cyclist exceeds 99% of terminal velocity. For reference, the difference in energy between the two cyclists molecules will asymptote at ln2*m*v_T^2.
If we use v_T=70m/s, then the distance will be approx 346 meters and the time constant will be approx 4.95 secs. If we use v_T=13m/s (approx. 29mph), then the distance will be approx 11.94 meters and the time constant will be approx 0.919 seconds. Of course, no real cyclist is powerful enough to accelerate that fast, so real distances and time constants would be somewhat larger.
This thought experiment demonstrates a "worst-case" scenario in which a cyclist gets no shading and starts with a velocity differential of -v_T to the leading cyclist. (A true "worst-case" scenario is one in which a cyclist has to stop for 30 seconds or so to change a flat, and must then catch up to the peloton by applying substantially greater power than normal to get back to the peloton; during this time, the cyclist is unshaded, unless he/she is lucky enough to have a teammate come back to share the load.)
Now if the cyclist "A" in front is travelling at terminal velocity v_T and the cyclist "B" behind is shaded, then the following cyclist B will accelerate relative to the front cyclist A and eventually pass him/her A. Since A was already travelling at terminal velocity, B will be travelling at greater than terminal velocity when B passes A, so after the pass, B will start decelerating and (due to shading) A will start accelerating. Now, since A is accelerating from v_T and B is decelerating from a velocity greater than v_T, A will now pass B again, but this time they will both be going faster than v_T. Clearly, the average speed of the two cyclists will converge at some speed slightly faster than v_T, as the time between passes of A & B converges to zero.
The details of whether the "cycle" of A&B passings converge, become cyclic, or diverge, depends upon the details of the shading function and the effective "inertial mass" of the riders. If the passings diverge, then at some point the shading function falls to zero and the cyclists become independent & never pass again. If the passings converge, then the cyclists asymptotically ride side by side at some terminal velocity greater than v_T. If the passings become cyclic, then the cyclists exchange places every n seconds and achieve an average speed slightly larger than v_T.
It would be a lot of fun to find the conditions under which the passings become "cyclic", so that the riders "orbit" one another (in phase space), and the behavior of sharing the lead is an emergent behavior!
At 05:35 AM 7/12/2010, Hans Havermann wrote:
While not necessarily germane to your approach, perhaps this online paper is of use: