Optimizing running a race on a curved track

In order to determine the optimal strategy to run a race on a curved track according to the lane number, we introduce a model based on differential equations for the velocity, the propulsive force and the anaerobic energy which takes into account the centrifugal force. This allows us to analyze numerically the different strategies according to the types of track since different designs of tracks lead to straights of different lengths. In particular, we find that the tracks with shorter straights lead to better performances, while the double bend track with the longest straight leads to the worst performances and the biggest difference between lanes. Then for a race with two runners, we introduce a psychological interaction: there is an attraction to follow someone just ahead, but after being overtaken, there is a delay before any benefit from this interaction occurs. We provide numerical simulations in different cases. Overall, the results agree with the IAAF rules for lane draws in competition, where the highest ranked athletes get the center lanes, the next ones the outside lanes, while the lowest ranked athletes get the inside lanes.


Introduction
In athletics, inside lanes are considered a disadvantage due to curvature, while in outside lanes, there is no one to chase.The aim of this paper is to understand from a physical and mathematical point of view the effect of the curved part of a track and of the lane number on the running performance both for a single runner and for a two-runner race.
To our knowledge, no optimal control problem including these effects has been studied.There is a huge literature on the way of running on a curved track, see for instance [1][2][3][4][5][6][7][8].Nevertheless, it is never coupled with the psychological effect to have a neighbor on the next lane, which is mentioned as important.Furthermore, though the IAAF regulations [9] do not impose a fixed shape of track, but allow the straights to vary between 80m and 100m, we are not aware of any study discussing the effect of the the lane and the track coupled with the psychological effect.
In this paper, we will build on a model introduced by Keller [10] and extended by [11,12], to investigate how the shape of the track and the centrifugal force change the optimal strategy in a race: this leads to longer race times for higher curvatures, and therefore favors the outer lanes.Estimating the performance of champions based on the modeling of Keller [10] has been developed by various authors [6,[13][14][15][16], but never taking into account so many parameters as in this paper.We will also introduce a model taking into account the psychological effect between two runners.This is made up of two effects: on the one hand, the attraction by a runner close ahead, and on the other hand, the delay before any benefit from the interaction occurs again after being overtaken.This delay model is inspired by a paper on walking [17].Let us point out that the mathematical problem encompassing delay in the equations is quite involved.We model the attraction by a runner close ahead as a decreased friction, since the focus on chasing someone ahead improves the runner's economy.Due to the staggered start positions in the curved part of the track, this "rabbit effect" is less favorable on the outer lanes.
After introducing the model, we perform simulations using the optimal control toolbox Bocop [18].Since the IAAF regulations do not impose a single shape of track, we analyze the effect of the shape of the track on the optimal velocity profile, as well as the influence of the various parameters of the runner for a single runner.Then we perform simulations for two runners and our results show that the combination of the centrifugal force and the two runners interaction brings a numerical justification to the fact that the central lanes are the most favorable to win a race.

Race model
Model for a single runner race Single runner on a straight track.When a runner is running on a straight, as used by Keller [10], according to Newton's second law, the acceleration is equal to the sum of forces.We can list two forces, the propulsive force f (t) in the direction of motion, and the friction force, that we assume to be linear in velocity.This leads to the first equation of motion for the velocity v(t) written by unit of mass: where τ is the friction coefficient.This coefficient models the friction due to various effects such as joints friction, the runner's economy and the elasticity of the track.
Other friction effects can be encompassed in the equation such as wind or slopes, that is going upwards or downwards, or banking.Because the runner has a limited capacity, the propulsive force is bounded from above by a constant f M , that is Typical values for f M range from 13 for a sprint at the world level [19] to 5 for a marathon [20].
The power developed by the propulsive force is f (t) • v(t), which is to be taken into account in the energy balance.This energy balance leads to the definition of the oxygen uptake σ introduced in [11], which depends on the anaerobic energy e(t).Note that at the beginning of the race, e(0) = e 0 , the available energy at initial time, and σ depends on the accumulated oxygen deficit at time t, that is e 0 − e(t).The function σ depends on the length of the race [21]: for short races (up to 400m), σ is a linear increasing function of e 0 − e.When the race gets longer, σ reaches its maximal value σ max in the central part of the race, but is increasing at the beginning of the race, and decreasing at the end.Note that σ is the energetic equivalent of V O2, the volume of oxygen used by a unit of time and σ max is related to V O2 max .For the shorter races considered in this paper, we assume a linear function σ and note σ f the final value, thus This leads to the energy model A champion-level runner has a V O2 max about 75ml/mn/kg.Since one liter of oxygen produces an energy about 21.1kJ via aerobic cellular mechanisms, this provides an estimate of the available energy per kg per second σ max = 75/60 * 21.1 ≈ 26m 2 s −3 .Furthermore, on a 200m race, V O2 and σ reach only about 75% of their maximal values [21,22], so we set σ f = 20.We point out that this term is of lower order than ė: we will see below in our simulations that the anaerobic part is roughly 87% of the total energy, which is perfectly consistent with [23].
For a fixed value of the final distance d, the optimal strategy to run the race is obtained by solving the control problem ( 1)-( 2)-( 3)-( 4) under the constraint: minimize T, such that This problem has been studied in [11,12,24].The parameters are matched to reproduce champions' races.For a race less than 400m, when the function σ(e) is decreasing, the velocity is increasing and then decreasing.Indeed, the runner never has enough energy to maintain his maximal force for the whole duration of the race.Therefore, the optimal strategy is to start at maximal force, and then the force decreases, and so does the velocity.
Since in the optimal control problem, it is usually the distance which is prescribed, in this paper, we choose to take the distance s instead of the time t as variable.We define y(s) to be the time required to run the distance s so that, if x(t) is the distance run in time t, we have We call f (s) the propulsive force needed at distance s and e(s) the energy.This allows us to derive the equations for y(s), f (s), e(s), from (1)-( 4), which are ÿ(s) = −f (s) ẏ3 (s) + 1 τ ẏ2 (s), y(0) = 0, and ẏ(0 ė(s) = σ(e(s)) ẏ(s) − f (s), e(s) ≥ 0, and e(0) = e 0 .
This formulation requires an initial velocity v 0 which is not zero, but given the effect of the starting blocks where our dynamical model is not correct, assuming an initial velocity of 3 or 4m/s is quite consistent with the effect of the beginning of the race, 10m from the start [25].The constraint on the force is The optimal control problem as such would lead to variations of the force which are too strong.In order to take into account the impossibility for the runner to vary his propulsive force instantaneously, we instead take df /ds as a bounded control.We seek the optimal race strategy to minimize T = y(d).
Centrifugal force on a curved track.For races of 200m or more, the track is not straight but includes one or more bends.While on a bend, the runner has to move against the centrifugal force, which, by unit of mass, is f c = v 2 /R where v is the velocity of the runner and R the curvature radius.In order to produce a mathematical model for the dynamics in the curved part, we have to take into account the centrifugal force in Newton's law of motion and project this equation on the 3 directions of motion.
Even on straights, there is an equation to be written in the z direction: the reaction of the ground, N , is equal to the weight.By the principle of action/reaction, the reaction of the ground is equal to the runner's propulsive force in the z direction.Note that the runner does not have his feet on the ground all the time in the stride: he rather 0, 2019 3/21 pushes (propulsive force) only for some time in a stride [19].Some remarks in [26] can be found related to this issue.We point out that there is an interesting explanation of the effect of arms to counterbalance the torque, and that since there are two legs, the reaction on each leg is not exactly the same [1].In this paper, we do not include these effects as we believe them to be of minor importance.The specificity of our work is that although we consider a mean force and mean velocity in a stride, our model allows us to compute an instantaneous force and speed along the race.On a curve, the runner makes an angle α with the vertical axis to balance the centrifugal force.The runner is subject to gravity g, to the reaction of the ground N along the angle α, and to the centrifugal force f c = v 2 /R (see Fig 1).One has to consider the equations of motion in the centrifugal direction and the z direction, which lead to which provides the angle according to the velocity and the value of N : By the principle of action/reaction, the propulsive force in the transverse direction is the opposite of the reaction of the ground in the horizontal direction, hence is equal to N sin α.Moreover, the propulsive force in the vertical direction is N cos α.The total propulsive force F is therefore such that F 2 = f 2 + N 2 where we recall that f is in the direction of movement.From (10), we find Since F has to be bounded and g is constant, this leads to the new constraint We point out that eventually the effect of the centrifugal force is taken into account in the force constraint.It cannot have an energy effect directly since the centrifugal force does not produce any work.Study of the track shape.It is important to know the exact shape of the track since it influences the runner's optimal pacing strategy and performance.However, there is no fixed regulation to build an athletic track.Actually, as indicated in the IAAF manual [9], the length of the straight part can vary between 80 and 100m, while the curved part can be a half circle ('standard' track) or two different circular sections ('double bend' tracks).We choose to study a standard track with an 84.39m straight part, and then two double bend tracks with straight parts of 79.996m and 98.52m respectively.The shapes and dimensions of theses tracks are detailed in Fig 2 .Note that for races longer than 100m, runners start the race in the curved part.The starting positions are therefore adjusted in order to have the same total distance for all lanes ('staggered start').
Note that each runner is assumed to run at a distance of 30cm from the inner limit of the lane.This is how the radius for the circular parts is set in order to obtain a 400m distance for lane 1.Then the width of each lane is 1.22m.This leads to different radii of curvature R k (s) depending on the lane k and the distance s run on the lane since the start.On the straight part, 1/R k (s) = 0.For more details on the value of R k (s) according to the track, we refer to the Appendix "Track Shape Details".
We want to point out that at the junction between the circular and straight parts, the runner will experience a discontinuity in the centrifugal force.This force is 0 on the 0, 2019 4/21  straight part and can reach a value of the order of 2.5N per kilo on the circular part (since v ∼ 10m/s and R ∼ 40m), which is about one quarter of the weight.We will see on the numerical simulations that this may lead to an acceleration of the runner when reaching the straight part of the track.One could think that it would be better to build a track where the curvature goes smoothly from 0 to the value of the matched circle so that the runner experiences a continuous variation of his centrifugal force.This type of curve, known as a clothoid, is used for instance for railways and roads.The simulations in Section "Numerical simulations for a single runner" indicate that the final time is actually larger on a clothoid, because the smooth transition leads to a smaller radius for the circular part, therefore a larger centrifugal force.
One of the main results of our simulations is that the tracks with shorter straights lead to better performances (see Section "Effect of different track shapes").
Final model for a single runner race.The optimal problem is to minimize 0, 2019 5/21 T = y(d) with y(s), f (s), e(s) solving ( 6)-( 7), σ being given by ( 3), with the bounded control df ds ≤ 0.015 (13) and the force constraint coming from ( 12) where the curvature radius R k (s) is prescribed according to the lane k and the track shape, see the Appendix "Track shape Details".We use the convention R k (s) = +∞ on a straight.Finally, introducing a state variable for the inverse of speed z(s) = 1/v(s), the optimal control problem for a single runner is

Model for a two-runner race
When two runners are involved, we label them with i, i = 1, 2 and define y i (s), f i (s), e i (s) respectively the time to reach the distance s, the propulsive force at distance s and the anaerobic energy left at distance s.We also label by i the parameters of each runner: τ i the friction coefficient, f M,i the maximal force, e 0,i the initial energy, v 0 i the initial velocity.Finally, we call T i the final time to reach the distance d that is T i = y i (d).
Objective function.We want to solve the race problem where both runners try to obtain their minimum time and win the race.The issue is to define a good mathematical problem.Minimizing min(T 1 , T 2 ) is not enough since it could lead to a situation where one of the runner stops optimizing his race once he knows he will lose.Then, minimizing the sum of the times T 1 + T 2 could lead to some cooperative interaction where the faster runner would wait for the slower one to optimize the global time.This is why we choose to minimize a combination of these two objectives, namely minimize with k w being a small parameter such that the second term does not modify the value of the leading order min(T 1 , T 2 ), but yet does not let max(T 1 , T 2 ), which is the time of the slower runner, be too big.In our simulations, values of k w ranging from 10 −3 to 10 −4 provide this kind of behaviour.
We point out that some authors [27,28] have tried to settle a stochastic description in the framework of game theory but they are not able to handle as many parameters as this model.Also in a short race, we do not believe that there is time to think and adapt one's strategy on the course of the race.
Psychological interaction.When two runners race against each other, we introduce an interaction term which mollifies the friction term of each runner ẏ2 i /τ i .This term is equal to 1 in case of no interaction, and is lower than 1 in case of a beneficial interaction.It models the psychological benefit that comes from chasing someone just ahead.Note that this interaction is not an aerodynamic effect ('drafting') 0, 2019 6/21 as in bicycle or car racing, because the velocity is too small.Cognitive effects are known to reduce perceived exertion: shielding has a psychological basis for runners and the focus on chasing produces better running economy [29][30][31].This psychological effect is indeed acknowledged by runners (sometimes called "rabbit effect") and can allegedly have an effect as high as 1 second per 400m lap [31].
The differential equations for y 1 and y 2 are therefore where F (y 1 , y 2 ) is to be determined as a function of r(s) which is the distance between the two runners.The detailed expression of r(s) is presented in the Appendix "Distance gap between two runners".Basic interaction.We choose the function F of r to be equal to 0.04 when r is roughly between 0 and −2.5m and 0 outside this interval, which corresponds to the distance for which an effect can be felt.A lot of possible functions can match this goal.We choose for instance the interaction function illustrated in Fig 3 where γ = 0.04, H a smoothed Heaviside function defined by and with the values for the offsets and slopes Lateral attenuation.It seems reasonable to assume that the positive interaction only occurs when the two runners are close enough to each other.Thus we introduce a limitation of the interaction based on the lane gap between the runners.In practice, the interaction is multiplied by an attenuation function A(k 1 , k 2 ) with A defined by Inhibition and delay problem.
A refinement of the interaction model is that the benefit should not hold just after being overtaken, since there is a delay in reacting.In order to build a mathematical model for this, we use the behavior of following in pedestrian traffic introduced in [17].The delay is meaningful in terms of human perception: perception of speed actually comes from successive perceptions of distance over time, and this integration process introduces a delay, while the perception of distances can be considered instantaneous.
We therefore introduce an inhibition formula that suppresses the interaction for a short duration after being overtaken.Since our model is formulated with distance as the independent variable, the delay is also expressed in terms of a distance frame η.The basic idea is to multiply the interaction term F by an characteristic function I η defined by For the numerical simulations, I η is smoothed using the Heaviside approximation H defined by (17).The detection of an overtaking is performed by checking for sign changes of r over [s − η, s].
Note that this check relies on past values of the state variables used to compute r, thus leading to a delay optimal control problem.Delay problems are a quite involved class of optimal control problems, and we refer the interested readers to [32,33] for recent theory developments.A classical way to solve delay problems is to reformulate them as non-delayed problems, see [34], but the manipulation is rather cumbersome.In our case, we take advantage of the feature from the toolbox Bocop to handle delays automatically in the fixed final time case (which we have since we use distance instead of time as the independent variable).
Final model for a 2-runner race.We define T i = y i (d) and F from (18).The optimal control problem becomes: It is worth pointing out that this optimal control problem has several families of local solutions, typically with a different number of overtakings.In the numerical simulations, we overcome this difficulty by trying several initial points and picking the best solutions.Using a global optimization method would of course solve this problem, however in our case the dimension of the state variables is too high.

Numerical simulations for a single runner
In the numerical simulations, we chose to study the 200m race.For reference, in 2018, the world record for 200m is 19.19s (Usain Bolt, Berlin World Championships, 2009).In all the following, we will simulate races with fictitious runners whose parameters (see Table1) are chosen so that their race times are close to 20s.We have chosen the maximal force f M to range between 6.5 and 13 which is the range in the 0, 2019 8/21 literature [3,[6][7][8]25].Then τ is chosen so that the peak velocity which is close to f M τ is roughly 11.1, which is the world's level.The initial energy e 0 is such that the ratio of the anaerobic contribution to the total contribution is 87% for runner A 1 , which is consistent with [22,23], and the final time is around 20 seconds.The initial velocity is taken to be 1/0.43 to take into account the departure in the starting blocks [25].The bound on the maximal variation of f is taken of order σ f /τ e 0 , which is close to what we expect as a singular control [35].

Single runner on a straight track
We start with a simple straight 200m race to illustrate the effect of parameters f M , τ , and e 0 .We take as reference athlete A 1 of Table 1.The corresponding speed and force profiles are shown with black lines in Fig 4 .The velocity increases to its peak value v m ∼ f M τ and then decreases.The runner does not have enough energy to run the whole duration of the race at maximal force.
The propulsive force starts at its maximal value f M , then decreases at the constant rate |df /ds| max .The time at which the force begins to decrease depends on the values of f M and e 0 .Indeed, increasing e 0 does not change the beginning of the race but allows to run longer at f = f M .On the other hand, increasing f M increases the peak velocity but does not change much the second part of the race.Finally, increasing τ has a more uniform effect and increases the velocity for the whole duration of the race.

Single runner on a standard curved track
We simulate the same runner on the so-called standard track, i.e. 115.61m half circle of radius 36.80mfollowed by a 84.39m straight.Fig 4 shows the race profiles obtained for the inner and outer lanes (respectively 1 and 8), and the straight race.The time splits for 50 − 100m, 100 − 150m and 150 − 200m are indicated in the figure: we have chosen the parameters for A 1 so that they match the order of magnitude of time splits for athletes in World Championships.The velocity profiles of the curved track are quite different from the straight track: i) the runner starts slower because of the curvature: even though he puts his maximal propulsive force at the start, part of it is used to counterbalance the centrifugal force, resulting into a lower effective force and a lower velocity ii) in the middle part of the race, the maximal propulsive force is reached and we can derive from (1) and (12) the relation between f and v: with R k the curvature radius on lane k.We can compare this formula with our simulations: on the straight v s = f M τ , while from (20), the velocity in the middle of the race on lane k is The numerical simulations indicate an extremely good consistency with this expression: in this case equation ( 21 iii) after the curved part, there is no more centrifugal force so that the runner can increase both his propulsive force and velocity.iv) finally, at the end, the runner slows down again, because he does not have enough energy left to sustain his maximal force.
If we compare lane 1 and lane 8, on lane 1 the runner starts slower since the centrifugal force is stronger due to larger curvature.On the other hand, he puts a slightly larger force in the second part of the race, having more energy left, yet he is slower overall.Final times are: 20.43 for the straight track, 20.46 for lane 8, and 20.48 for lane 1.Let us point out that our simulations are consistent with the experiments in [3], where runners are asked to run 60m on a straight path and on a curved path.The authors observe the existence of two groups, one "good" group who manages to reach the same velocity in the curved path as in the straight path and the other "poor" group who is strongly affected by the curve.Our parameters values of runner A 1 corresponds to a runner of the "good" group.Now we simulate several runners (see Table 1 for parameters) in order to assess the influence of the maximal force f M .
Runner with large maximal force f M .We want to point out that due to the way the curvature is taken into account in the model, see (12), a runner with a greater maximal force f M will be less sensitive to the curvature of the track.We illustrate this with the runner A 2 defined in Table 1 whose f M = 13; final times are: 20.31s for straight track, 20.32s for lane 8 and 1.In this extreme case, the runner is basically unaffected by the curvature of the track, that is the curves of velocity and force versus distance are almost the same for straights, line 1 and 8.
Runner with small maximal force f M .With a low maximal force f M = 6.5, the runner A 3 of Table 1 can increase his force and velocity when he reaches the 0, 2019 10/21 straight part of the track, since he has not spent as much energy as the others at the beginning, yet he is slower overall.This is illustrated in Fig 5 .Final times are: 20.32s for straights, 20.54s for lane 8, and 20.72s for lane 1.This runner A 3 corresponds to a runner of the "poor" group of [3], with as much as 0.2s gap between extreme lanes.It is also consistent with performances for runners in [6][7][8]36].Force in N/kg vs distance on the right graph.Speed vs distance on the left graph, with the constant speeds given by eq. ( 21) in dashed lines.Force and velocity increase when the centrifugal force disappears.
Let us point out that for these three runners, we have computed the ratio of the anaerobic energy versus the total energy and find: 87% for A 1 , 90% for A 2 and 82% for A 3 , which is quite consistent with [23].

Effect of different track shapes
Now we study the effect of different track shapes defined in Fig 2 : standard with 84.39m straight (STD), double bend 1 with 80m straight (DB1), double bend 2 with 100m straight (DB2), and two modified standard tracks with smoothed curvature, including clothoid junctions of 10m (CL1) and 30m (CL2).For the clothoid tracks, we choose a straight of 84.39m as the standard track.As explained in the Appendix, the length of the junctions provides the radius of the circle and the angle, which are respectively 33.32m and 164°for (CL1) and 29.95m and 118°for (CL2).
For each track shape, we simulate the race on the inner and outer lanes (1 and 8).The results are summarized in Table 2, with the races for runner A 1 (on lane 5) shown in Fig 6 .Reference athlete A 1 has a difference of 0.17s between the best (DB1 track, lane 8) and worst (DB2, lane 1) case.As mentioned previously, runner A 2 with a very high force f M = 13 is almost unaffected by the curvature, with times varying only between 20.32s and 20.36s.Yet, the DB2 track is still worse than the others.Conversely, athlete A 3 , with a lower force f M = 6.5, is more affected, with 1.01s between the best and worst cases.The DB2 is his worst track and his best performance is on the standard track.
Our results show a time difference between inner and outer lanes ranging from 0.02s for the standard track to 0.15s for the worst double bend track.This is consistent with [6] who also finds the double bend track to be the worst, using a simplified model based on constant mean velocity and curvature.Let us point out that in the next section, we will study a two-runner race where the effect of the lane becomes more pronounced: we find a larger difference between the best and worst mean time per lane.
Focusing on runner A 1 in Fig 6, we analyze more closely the effect of the track shape and lane: 0, 2019 11/21 -DB1 is the quickest track for the outside lane, though it is very close to STD.
-The standard track has the smallest difference between lanes.
-DB2 is the slowest track, from 0.01s on the outside lane to 0.14s on the inner lane.When on the outer radius of curvature 24m, the velocity significantly decreases.
-CL1 is quite close to DB1 and STD, though a little slower.CL2 is slower than DB2 on the outer lanes, although not as bad in terms of difference between lanes.
It may seem surprising that the tracks with smoothed curvature do not perform better than the ones with a discontinuous curvature.This comes from the fact that the clothoid junction actually results in a smaller radius for the circular part, and thus a greater curvature.The longer the clothoid junction, the more pronounced the effect, and the slower the times.
To conclude the single runner races, it appears that the track with the shortest straight is the quickest track for strong athletes in outer lanes.The standard track shape is the one with the best race times overall, and also the smallest time gap between the inner and outer lanes.On the opposite, the double bend with the long 100m straight (DB2) yields the worst times overall, and the highest gap between the inner and outer lanes.These conclusions seem consistent with runners' feelings though there is no study yet of what the ideal shape of track would be for a specific runner.

Numerical simulations for two runners
We move to the simulations for two-runner races, combining the interaction effect with the curvature effect previously studied for the single runner case.Firstly, we study races with two runners competing in adjacent lanes, to see the effect of the interaction.Then we compute the mean times corresponding to all possible races of a runner versus himself and find that the best lanes are indeed the center ones.

Races on different lanes and illustration of the interaction effect
We perform simulations for the optimal control problem (OCP2) for two runners, combining the interaction effect with the centrifugal force.We first set A 1 to be the runner on each lane 1 and 2. We recall that if A 1 runs alone, his time on lane 1 is 20.485s and on lane 2 20.480s, so of course because of the centrifugal effect, lane 2 is quicker.Due to staggered starts, as soon as we set the interaction, the runner on lane 1 benefits from the interaction at the beginning of the race.First, we set the interaction term γ = 0.04 but with no inhibition η = 0.The results are illustrated in Fig 7, with the velocity profile in lane 1 on the left and the interaction for each runner and relative distance on the right.When the relative distance is negative, the runner in lane 1 is behind.So in this case, the runner in lane 2 wins the race and they overtake each other twice: lane 1 starts behind because of the staggered starts, benefits from interaction and overtakes at 50m; then lane 2 benefits from interaction right away and is able to overtake again at 150m.Then they are on the straight, very close to each other, lane 1 benefits from interaction and is ready to overtake again but loses in the end by 0.04s.0, 2019 12/21 Then in Fig 8, the interaction term is set at γ = 0.04, and the inhibition frame is η = 20m.This means that the positive interaction is disabled when a runner is overtaken in the previous 20m of race.Fig 8 shows the speed profile (left graph) and interaction / inhibition terms (right graph).Compared to the race without inhibition in Fig 9, we observe a different behaviour with only one overtaking and the runner on the inside winning by 0.27s.Note that since we optimize the whole race, there is no reason for the race with inhibition to coincide with the race without inhibition, even before any overtaking occurs.We observe that the inhibition (on the right graph) correctly detects the overtaking and suppresses the interaction accordingly.This prevents the overtaken runner at lane 2 to keep up (and eventually catch up) with the one at lane 1, as we see that the distance gap increases after the overtaking.In the race without inhibition, the overtaken runner was benefiting from the interaction right away, which allowed him to catch up and take the lead back.With inhibition, the runner on lane 1 manages to win the race, though he is on a disadvantageous lane.In the full race, of course, the runner in lane 2 has a neighbour on the other side which changes the total result.
There are cases where, though the inhibition η = 20, there are still two overtakings.We study for instance the 5-4 race, with the speed and force profile of the runner at lane 5 shown in Fig 9 .Without interaction (γ = 0), lane 5 wins without any overtakings, with final time 22.47s.With interaction (γ = 0.04 and η = 20), lane 5 still wins after 2 overtakings, with final time 22.23s.At the start, the runner on lane 4 benefits from the interaction due to the runner at lane 5 being ahead (staggered start).He catches up then overtakes the outer runner, who in turn gains the interaction, catches up and overtakes the inner runner again.At the end the inner runner, being behind, has the interaction again and is catching up with the outer runner, but too late.
We have also made simulations with runner A 1 vs runner A 2 , and though runner A 2 0, 2019 13/21 is stronger in force, on some lanes, runner A 1 can benefit from interaction to be able to win.
0, 2019 14/21 We point out that the interaction parameters can be runner dependent since some may be very sensitive to this effect and others much less.

Mean time per lane
In a real race, there are eight runners, however our model is only for two.Therefore, we simulate a set of races with two identical runners, the first on a fixed lane, the second on each possible other lane.We define T k1,k2 1 to be the time for the winner in the race between two identical runners A 1 in lanes k 1 and k 2 .We want to compute the average performance at lane i as the mean time Ti = 1 7 j=1..8,j =i .
First, we compute the times T i,j 1 : the best times in j for each i are indicated in Table 3.In Fig 10, we have plotted the times for i = 1, 5, and 8.The best times are obtained for the maximal interaction, namely with the second runner on an adjacent lane.For runner A 1 on lane 5, his best performance is obtained with a neighbor on lane 4 rather than 6.We recall that the model includes a lateral attenuation for the interaction, which is 0 when runners are more than 3 lanes apart.If we compare the best time for each case, it is decreasing with the lane.against himself, with an interaction weight γ = 0.04 and η = 20 when he runs on each lane i.If we look at the overall performance then lane 5 is the best, followed by lane 6, 4, 7, 3, 8, 2 and lane 1 is by far the worst.We compare with the solo case (γ = 0) where of course the outside lane is the quickest.
The results are nicely consistent with the IAAF rules for the lane drawn.Indeed, according to the IAAF rules [9], starting lanes are drawn in three lots: -a first draw is made for the four best runners in the center lanes 3,4,5 and 6.
0, 2019 15/21 -a second draw is made for the next two runners between the outer lanes 7 and 8.
-a last draw is made between the runners with the lowest performance to get the inside lanes 1 and 2.
Nevertheless, we find that the inside lanes 1,2 are a real disadvantage, the more so as if the runners are not as strong.
In [4] the authors recall some average time data for Olympics 1996 and 2000, and World Championship 2001: they indicate an average time gap of 0.16s between inside lanes 1 and 2 and outside lanes 7 and 8.We obtain a smaller gap of 0.047s, which may be due to the fact that we consider identical runners in our simulations, while in actual races the athletes in the outside lanes were supposedly stronger than those in the inside lanes.11.Mean times per lane for runner A 1 when in lane i vs himself in all other lanes.Without interaction (γ = 0), and with interaction γ = 0.04, η = 20.Lane performance (sorted by mean time): T5 < T6 < T4 < T7 < T3 < T8 < T2 < T1 .Gap T1 − T8 = 0.050859.

Conclusion
In this paper, we have studied how the geometry of the track and the psychological interaction between runners affect performances.We have introduced an optimal control model taking into account the centrifugal force as a limiting factor for the maximal propulsive force.We couple this with a new model describing the positive interaction exerted by a runner close ahead and the delay to benefit from it after being overtaken.We carry out numerical simulations for several runner profiles on different track shapes.The results indicate that the track with the shortest straights is the quickest for strong athletes in outside lanes.The so-called standard track (two straights and half circles) yields the best performances overall.The double bend tracks with longer straights (DB2) are significantly slower.In particular running on lane 1 on the DB2 track appears to be an overwhelming disadvantage.
Furthermore, the combination of the centrifugal and interaction effects leads to the center lanes being the most favorable, followed by the outside lanes, with the inside lanes being the worst.These results fit very well with the IAAF rules for lane draws, which follow this preference order.0, 2019 16/21 Denoting by θ 0 k = µ k π Ro the starting angular position, the angle after running s meters is

Clothoid track
Let us study the design of a modified standard track in which the circular part of radius R is bracketed by two smoother junctions of length l with continuous curvature.We choose to use a junction whose curvature is linear with respect to the distance, called a clothoid (also known as Euler curve or Cornu spiral).The angle with the tangent ϕ is ϕ(s) = sc(s) 2 where c(s) is the curvature at distance s.Then, since the total angle for the circular part and the two clothoid junctions is π, and using the same notations as before that is l s is the length of the straight part and l c the length on the circle, we find Since the total angle for one clothoid is 2ϕ( l) = l/R, this equation leads to l c + l = Rπ.Moreover l s + l c + 2 l = d where d is the distance of the race, that is 200m in our case, thus Rπ = d − l s − l.We find therefore that when there is a junction with a clothoid, the radius of the circular part gets smaller than in the case of a full half circle.On the clothoid, for s ∈ [s begin , s end ] , the expression of the curvature is linear: c(s) = c begin s end − s s end − s begin + c end s − s begin s end − s begin .
In our case the clothoids will join the straight part (curvature 0) and circular part (curvature 1/R).
Similarly to the double bend track, we denote by C 1 , C 2 , C 3 , S respectively the first clothoid, circular, second clothoid, and straight parts.We denote by l s the straight length, R, l c the radius and length of the circular part, and l1,2 the length of the two clothoid junctions.As before, we call R k = lc+ Taking these variable lengths into account, the curvature at lane k after running s meters is

Distance gap between two runners
For the interaction term, we need to define the relative distance between the runner on lane 1 and on lane 2, taken by convention at time y 1 (s): ρ(s) = x 1 (y 1 (s)) − x 2 (y 1 (s)) = s − x 2 (y 1 (s)).
On a curved track, this approximation is adjusted by projecting the two runners on a median circle, while also taking into account the staggered start on different lanes: , where c i (s) = 1/R i (s) is the curvature at distance s.

Figure 1 .
Figure 1.Illustration of the forces on the runner

Figure 4 .
Figure 4. Single runner A 1 on a standard track, lanes 1 and 8, and straight track.Force in N/kg vs distance on the right graph.Speed vs distance on the left graph, with the constant speeds given by eq.(21) in dashed lines.Time splits for 50 − 100m, 100 − 150m and 150 − 200m are indicated.

Figure 5 .
Figure 5. Single runner A 3 on a standard track, lanes 1 and 8, and straight track.Force in N/kg vs distance on the right graph.Speed vs distance on the left graph, with the constant speeds given by eq.(21) in dashed lines.Force and velocity increase when the centrifugal force disappears.

Figure 6 .
Figure 6.Effect of the track shape on the race time for Runner A 1 vs lane number (left graph).Speed profile for lane 5 (right graph) with zoom.The tracks are standard (STD), double bend 1 (DB1) with short straights, double bend2 with long straights (DB2), and curves with short and long clothoid junctions (CL1 and CL2).

Figure 7 .Figure 8 .Figure 9 .
Figure 7. Race A 1 vs himself at lanes 1-2, with interaction γ = 0.04 and a frame η = 0 that is no inhibition.Left graph: speed profile and time splits of the runner at lane 1.Right graph: distance gap and interaction term for both runners.The sign change of the distance gap corresponds to the overtaking.Lane 2 wins by 0.04s.

l π + 1 .
22(k − 1) and l c,k = l c (1 + 1.22(k − 1)/R) the radius and length of the circular part at lane k.On lane 1, both clothoids have same length l, while for k > 1 the first clothoid is shorter in order to keep the same total length.The second clothoid has full length lk,2 = lk = l(1 + 1.22(k − 1)/R).Thus the first clothoid has length lk,1 = d/2 − l s − l c − lk,2 .

Table 1 .
Athletes' parameters.A 1 is a good runner, A 2 has a very strong propulsive force and A 3 is a poor runner.

Table 2 .
Times for different runners and track shapes.

Table 3 .
Athlete A 1 at lane k running against himself at lane k − 1. Interaction γ = 0.04 with inhibition η = 20m.Race time and gain with respect to solo race time.