On-site cycling drag analysis with the Ring of Fire

The Ring of Fire (RoF) measurement concept, introduced by Terra et al. (Exp Fluids 58:83. https://doi.org/10.1007/s00348-017-2331-0, 2017; Experiments in Fluids 59:120, 2018), is applied to real cyclists to enable the aerodynamic drag determination during sport action. This principle is based on large-scale stereoscopic particle image velocimetry (PIV) measurements over a plane crossed by the athlete during cycling. The momentum before and after the passage of the athlete poses the basis for the control volume analysis in the athlete’s frame of reference, which returns the aerodynamic drag. This approach extrapolates aerodynamic studies towards more realistic conditions, compared to experiments performed in wind tunnels with scaled or stationary athletes. The measurement concept is termed Ring of Fire as the rider crosses a region of intense light. Two experiments are conducted, indoor and outdoor, with attention placed on the effects of the environmental conditions and the confinement of the measurement region. Stereo-PIV measurements feature a plane of approximately 2 × 2 m2, using neutrally buoyant sub-millimeter helium-filled soap bubbles (HFSB) as flow tracers. The drag measurement is obtained examining the wake produced by the athlete. It is observed that the drag value becomes independent of time after about 5 torso lengths from the passage. A statistical estimate of the drag is produced combining the results of several passages. Fluctuations of the drag value during a single passage are associated with the unsteady wake flow. Overall fluctuations among different transits are ascribed to the varying conditions of the airflow prior to the passage of the athlete. The experiments conducted outdoor exhibit significantly larger dispersion of the drag value, compared to the quieter conditions indoor. Repetition of the transit 10–30 times yields a basis for statistical convergence of the average drag value. The flow topology past the cyclist compares satisfactorily between both experiments and with wind tunnel experiments reported in literature. The current measurements clearly separate drag values from upright and time–trial athlete’s positions, indicating the suitability of this principle for aerodynamic analysis and optimization studies. .


Introduction
Most experimental research in sport aerodynamics is performed in wind tunnels, despite the fact that the dynamical situation to be simulated (e.g., a cycling or running athlete) poses challenges related to the athlete motion and its control. As a result, the problem is often simplified reverting to a stationary scaled model to match the constraints posed by the wind tunnel size and the measurement techniques used for the aerodynamic analysis. The aerodynamic force is directly measured by connecting the model to a force balance. Alternatively, the drag force can be derived from velocity measurements in the wake of the object, carried out either via Pitot rakes (e.g., Jones 1936) or particle image velocimetry more recently (Kurtulus et al. 2007;van Oudheusden et al. 2007;David et al. 2009). The latter principle invokes conservation of momentum in a control volume that encloses the object. The deficit of momentum flux past the object corresponds to the aerodynamic drag acting on it (van Oudheusden et al. 2006).
In some cases, wind tunnel experiments cannot reproduce the aerodynamic conditions with sufficient accuracy and experiments whereby the model moves in quiescent air are considered. Typical limitations encountered in wind tunnel tests of sport aerodynamics range from inaccurate scaling of shape and roughness, model blockage, interference of the support (Barlow et al. 1999) or, for instance, when the flow behind accelerating objects is to be dealt with (Coutanceau and Bouard 1977). Additionally, the wake development far downstream of the model can be performed with more advantages when the model is towed (Scarano et al. 2002) or catapulted in quiescent fluid (Von Carmer et al. 2008). The transiting model approach has been successfully adopted to investigate also ground vehicle aerodynamics (Jönsson and Loose 2016) or animals in free flight (Hedenström and Johansson 2015;Ben-Gida et al. 2013).
The recent works of Terra et al. (2017Terra et al. ( , 2018 with tomographic PIV in combination with helium-filled soap bubbles to determine the drag of a towed sphere can be seen as preliminary to the current study. The use of HFSB offers the potential to upscale the measurement region up to several square meters as demonstrated by Bosbach et al. (2009).
In the present study, a measurement apparatus is realized that quantifies the aerodynamic drag of a full-scale cyclist during sport action. The experimental procedure to achieve drag measurements follows the same principles discussed by Terra et al. (2017). The measurements are performed by large-scale stereoscopic PIV over a field of view of about 4 m 2 . For such approach, where the rider crosses the illuminated measurement plane, the experimental method is referred with the name "Ring of Fire" (RoF).
Here, the RoF concept is applied for the study of the aerodynamic drag in cycling. The latter dominates the forces opposing the athlete's motion, which justifies the attention devoted to aerodynamic drag in several studies (Kyle and Burke 1984;Wilson 2004;Lukes et al. 2005;Crouch et al. 2017 among others). The flow field around a pedaling cyclist features a complex system of vortices, in turn depending also upon the cyclist's torso and legs position along the crank cycle. These latter variations in aerodynamic drag cannot be solely ascribed to changes in frontal area, but result from complex aerodynamic interactions leading to different flow and vortex topology. Crouch et al. (2014) have produced a detailed aerodynamic survey in the wake of a cyclist by wind tunnel experiments. The work resulted in the identification of the most prominent streamwise vortices emanating from the athlete at different positions during pedaling.
The detailed velocity field around a full-scale cyclist model has been recently measured with robotic particle image velocimetry (Jux et al. 2018), where also the near field flow topology has been characterized.
The RoF experiments aim at determining the aerodynamic drag of the cyclist during sport action, such to obtain an estimate close to on-site conditions. The results are to be compared not only to the above mentioned wind tunnel studies and computational fluid dynamics (CFD) simulations, but also to other techniques currently practiced for on-site measurements (coast down, Petrushov 1998; torque power output, Grappe et al. 1997). Moreover, they support the correlation between the aerodynamic drag and the flow field by quantitative visualizations of the velocity field in the cyclist's wake.
The present work describes the realization of the RoF concept for full-scale sport aerodynamics and discusses the experimental procedures for indoor and outdoor experiments, mimicking, respectively, track and road cycling. The aerodynamic drag estimation from cyclists during sport action is compared to literature data from wind tunnel experiments and other techniques. Furthermore, the experiments cover different postures of the cyclist (time trial and upright) with the aim to directly measure the effect of posture on aerodynamic drag and its detectability with the RoF.

Working principle
The aerodynamic drag of an object moving in a fluid can be evaluated invoking the conservation of momentum expressed in a control volume. A recent review of the problem has been provided by Rival and Van Oudheusden (2017). The formulation of the problem is simplified from its unsteady form to the steady condition when applying a Galilean transformation (Arnold 1989), whereby the frame Page 3 of 16 90 of reference moves with the object and the fluid is considered in uniform motion upstream of the object (Terra et al. 2018). If the control surfaces S 1 and S 2 are sufficiently far from the object surface, it can be shown that the viscous stress is negligible (Kurtulus et al. 2007). The drag force can, thus, be expressed as follows: In the expression above is the air density. The cyclist is assumed to move at constant speed u C with respect to the laboratory frame of reference. Let us now consider a control surface normal to the direction of motion of the cyclist. Prior to the passage of the cyclist, the air motions feature a chaotic velocity u env , resulting from the environmental effects, as depicted in Fig. 1-top. Assuming uniform and quiescent conditions prior to the passage would largely simplify the problem formulation. However, even in scaled experiments, the disturbances in the air motion induced by the environment and the seeding generation are reported not to be negligible (Terra et al. 2018). (1) After the passage of the cyclist, the flow velocity features a coherent wake with a velocity profile u wake that follows the moving cyclist. Making use of a Galilean transformation, the representation of velocity and momentum changes from the laboratory to the cyclist frame of reference moving at speed u C . As a result, the air flow velocity ahead of the cyclist can be written as U ∞ = u env − u C , while that at its back as ( Fig. 1-bottom): This expression is valid at the condition that the mass flow is conserved across S 1 and S 2 . This is ensured by shrinking the inlet plane (S 1 ) from the outer edges, starting from the equal size as that of the outlet plane (S 2 ).
Equation (2) yields the instantaneous aerodynamic drag from the surface integral of momentum and pressure over a fixed plane before and after passage of the cyclist. Ensemble averaging (Eq. 3) of the drag among multiple passages is performed to achieve a higher degree of statistical convergence.
where N is the number of model passages. The aerodynamic drag exhibits temporal fluctuations associated with the unsteady nature of the flow around the cyclist. However, these unsteady fluctuations are little relevant to the evaluation of the cyclist's drag, given their short time scale. Time averaging is, therefore, performed within the ensemble average (Eq. 4) with the objective of reducing the effect of the unsteady fluctuations on the evaluation of the time-average drag.
where T is the total number of time steps and D(t i ) is the ensemble average drag at each time step in the wake.

Experimental setup and procedure
Experiments were conducted with a cyclist riding a timetrial (TT) bike. For the indoor case, the cyclist was male, 1.89 m tall, with weight of 68 kg. He wore a short-sleeve  time-trial suit from Team Sunweb and a Giant Rivet TT helmet. The cyclist riding through the outdoor setup was 1.84 m tall and weighed 83 kg at the moment of testing. He was equipped with a long sleeved time-trial suit from Team Blanco and a Lazer Wasp TT helmet. Moreover, for safety reasons, both cyclists wore a pair of laser goggles. The approximate torso chord length, c, for both athletes, is 600 mm. In the indoor experiment, a Giant Trinity TT Advanced Pro bike with 2 × 11 gears was used, while a Ridley Cheetah TT bike with 2 × 9 gears was used for the outdoor experiment.

Experimental facilities and cycling conditions
The experimental facilities and test conditions are presented in Table 1. The top view of the sport hall and of the outdoor site is shown in Fig. 2. The flow tracers are generated and confined within a tunnel of 4 m × 3 m and 3 m × 2 m (width × height) for the indoor and outdoor experiments, respectively. Curtains are used to maintain a high concentration of tracers within the duct. The entrance and the exit in the outdoor experiment are closed during accumulation and opened prior to the transit of the cyclist. For the indoor experiment, a curtain at the exit was sufficient. The measurement plane is near the half of the duct. Considering the small blockage ratios of 3.5 and 7% for the indoor and outdoor experiments, respectively, a non-confined environment is assumed for the control volume approach. The floor is covered with a thin carpet (polypropylene, 3 mm) to avoid ground slipperiness due to the PIV seeding. A photograph of the setup during experiments is shown in Figs. 3 and 4. During the indoor experiment, the cameras were positioned 6 m upstream of the duct entrance.
Although the two experiments have similar acceleration length before the measurement plane (Fig. 2), the limited available braking length in the indoor experiments requires conducting the tests at lower velocity (5.3 m/s). The crank angle is defined as the angular position of the right foot (forward) with respect to horizontal crank position .
In both the cases, measurements are conducted with the cyclist in upright and time-trial position (see Fig. 5). Following Crouch et al. (2014), the pedaling frequency (cadence) is normalized with the advancing speed, k = 2r f u C , where r is the bike crank length, f the cadence and u C the cyclist velocity, as reported in Table 1. The reduced frequency is k = 0.12 indoor and k = 0.23 outdoor, respectively.

PIV instrumentation, imaging and data processing
Velocity measurements are performed with a large-scale stereoscopic-PIV system. The experimental parameters are presented in Table 2. Neutrally buoyant helium-filled soap bubbles (HFSB) are used with an average diameter between 0.3 and 0.4 mm, providing sufficient light scattering to visualize a field of view (FOV) of the order of 4 m 2 . The tracers In the indoor experiment, a low-repetition rate PIV system is used, whereas the outdoor experiment features high-speed PIV (see Table 2 for specifications). The results are not affected by the selection of the hardware, which is different only due to the availability at the time of the  experiment. The low-speed system benefits from the higher pulse energy and sensor resolution with well-resolved particle images (diffraction disk covered with approximately 2 pixels). On the other hand, the high-speed system offers three orders of magnitude higher temporal resolution, enabling more advanced data processing, at the cost, however, of a lower imaging resolution (diffraction disk imaged over 0.5 pixels). The pulse separation with the low-speed system is chosen considering the out-of-plane loss-of-correlation factor (Keane and Adrian 1992). A cross-correlation analysis with multigrid image deformation (Scarano and Riethmuller 2000) is employed. A typical recording of particle images is shown in Fig. 6 for both experiments. The more controlled environment conditions in the indoor experiment result in a more uniform dispersion of the tracers and PIV images with homogeneous concentration. Achieving uniform seeding distribution in the outdoor experiment is hampered by the effect of wind gusts. From the raw PIV images, the cyclist's crank angle at the moment of the passage through the laser sheet is determined with an accuracy of ± 10°.

Data processing
The recorded images are analyzed with the LaVision DaVis 8 software. The pre-processing removes background light by subtracting the minimum intensity over time at each pixel. The recordings from the indoor experiment are analyzed with dual-frame cross-correlation. The time separation between frames is set to 3 ms. A sliding sum-of-correlation algorithm (Sciacchitano et al. 2012) is employed for the outdoor experiment. For the latter, the analysis performs an average of the correlation maps from seven pairs of frames sliding a time interval of 3.5 ms. The time separation between frames is set to 2 ms. To quantify the range of resolvable velocity scales, the dynamic velocity range (DVR) is determined as the ratio between the maximum velocity in the near wake of the cyclist and the standard deviation of the velocity distribution in the quiescent flow prior to the cyclist's passage. Details of the image processing parameters and estimates of the measurement dynamic range are summarized in Table 3. The drag force evaluation after one passage of the cyclist is obtained via Eq. (2). The velocity field prior to the passage of the cyclist is significantly weaker than in his wake. Averaging the measurements before passage over a short time interval (1.25 s and 0.1 s for the indoor and outdoor experiment, respectively) reduces the effect of measurement noise in the determination of u env . To further reduce the measurement noise in the drag estimate, a wake contouring approach is applied which isolates the cyclist's wake from the outer flow region. The wake is defined as the flow region whose velocity is below a certain fixed percentage (5% in the present case) of the minimum velocity in the flow field. Such region is then dilated by two adjacent vectors to include also the shear layers, thus obtaining the outlet surface S 2 of Eq.
(2). The inlet surface S 1 is obtained by shrinking S 1 in all directions up to the point that the conservation of mass is satisfied.
The cyclist's speed is monitored measuring the bicycle transit time across the light sheet. In the indoor experiment, a magnetic sensor provides the cyclist speed in real time, additionally.
The wake past the cyclist exhibits unsteady behavior. Consequently, also the evaluation of the drag force yields temporal variations. A statistically significant estimate of the cyclist's average drag is produced by ensemble averaging (Eq. 3) the velocity field obtained from 10 and 28 repeated measurements for the outdoor and indoor conditions, respectively.
Two main repeatability issues are identified that require a specific treatment of the instantaneous data to retrieve ensemble average flow fields: (1) since the cyclist crosses the measurement plane at a different Y coordinate for every passage, the measured velocity field is relocated in the Y-direction to compensate for such shift; (2) the relative distance between cyclist and measured wake planes is not exactly the same among different passages; the exact streamwise relocation is obtained examining the position of the cyclist when he crosses the measurement plane. For the latter problem, the high-speed PIV system resolves the motion of the cyclist within few millimeters in the streamwise direction; therefore, any error associated with variations of the relative distance between cyclist and wake planes can be neglected.
The results are presented in the coordinate system as shown in Figs. 3 and 4 with t = 0 defined as when the rearmost point of the saddle crosses the laser sheet. To make the comparison between results from both experiments possible, the flow field variables and time are made dimensionless in the following way: The dimensionless streamwise velocity u * x is written in the frame of reference of the cyclist, meaning that when u * x = 0 , the velocity deficit is equal to the cyclist velocity, and when u * x = 1 there is no velocity deficit (equivalent to freestream conditions).
The uncertainty of the estimated C d A values is analyzed a posteriori, based on the standard deviation of the instantaneous drag area estimates, and the number of independent samples (considering both the number of passages of the cyclist, and the number of independent flow measurements in the wake of a cyclist during one passage). A detailed analysis of the measurement uncertainty and drag resolution of the Ring of Fire system for small-scale applications is reported in a recent work of Terra et al. (2018), where the effect of simplifications in the conservation of momentum equation is considered.

Measurement procedure
Before the passage of the athlete, the duct curtains are closed and the HFSB accumulate for approximately 2 min. Atmospheric wind conditions require continuous operation of the seeding generator for the outdoor experiment. Instead, in the indoor experiment the bubbles production is paused prior to the passage of the cyclist and the momentum disturbance introduced by the seeding rake micro jets decays.
Dimensionless streamwise vorticity * The cyclist starts from the same predefined distance and crank angle for each passage, to have well-matching athlete posture (leg position) in the measured area between passages. In the indoor experiment, the image acquisition is triggered by a photoelectric sensor, while the user manually triggers the image acquisition in the outdoor experiment. Transferring the acquired images to mass storage requires 5 min with the high-speed PIV system, whereas typically 40 image pairs are recorded with the low-speed PIV system permitting to repeat the experiment within 1 min.

Air flow conditions before cyclist transit
The conditions before the passage of the cyclist rarely exhibit fully quiescent air. The environmental flow motions feature a velocity u env , which is in general non-zero, nonuniform and non-stationary, mostly due to external conditions and the seeding injection. An instantaneous flow field before the cyclist's passage is illustrated in Fig. 7 for both the indoor (left) and outdoor (right) experiments. To reduce the noise in the data, the velocity is averaged in time during 1.25 s (indoor) and 0.1 s (outdoor) before the passage of the athlete.
The indoor experiment was performed in a closed, thus, quieter environment; whereas during the outdoor experiment, the presence of moderate wind (0.5-1 m/s) could only be partly attenuated by the walls of the tunnel. This is clearly visible in Fig. 7, where the environment velocity is of the order of a 5 cm/s in the indoor experiment and attains 30 cm/s outdoor.
The velocity distribution prior to the passage is taken into account for the drag computation via Eq. (2) as it contributes to the overall momentum budget, as also discussed by Terra et al. (2018). Furthermore, unsteady effect may influence the interaction of the wake with the initial velocity field, resulting in variations of the measured drag. The latter effects, however, are neglected and cannot be directly observed with the current experimental apparatus.

Velocity field in the cyclist wake
The flow fields in the wake of the cyclist are discussed for the indoor upright and time-trial configuration as well as for the outdoor time-trial configuration. Figure 8 shows a comparison of the instantaneous streamwise velocity u * x at t* = 3. Note that the cyclist contours in Fig. 8 are meant to indicate the general cross section of the athlete and do not reproduce the exact position of the legs. The development of both an indoor as well as an outdoor instantaneous wake over time is available online as supplementary material. First, a comparison between time-trial ( Fig. 8-left) and upright positions (Fig. 8-middle) for the indoor experiment is given. The magnitude and location of the peak momentum deficit are similar in both cases. The out-of-plane velocity contour of the wake ( u * x = 0.95), however, is clearly wider for the upright case. Interestingly, it has the same height for the time-trial position as it has for the upright position, despite a higher height of the cyclist in upright position.
Next, the time-trial position is compared between the indoor (Fig. 8-middle) and outdoor experiments (Fig. 8right). The wake observed in the outdoor experiment is wider and shows a slightly higher peak momentum deficit. Despite that the heights of both cyclists in time-trial position were very similar, in the indoor experiment, the u * x = 0.95 contour is consistently higher (see also Fig. 9). A reason for this can be the different inclination angle of the torso of both cyclists, generating a different amount of downwash over the back.
The temporal development of the ensemble average streamwise velocity field (u * x ) past the cyclist is shown in Fig. 9. The ensemble average is obtained from 28 and 10 individual runs from, respectively, the indoor and outdoor experiments. The maximum deficit in the wake (~ 45%) is observed at the shortest time delay after the passage. The deficit is not uniformly distributed and attains its maximum behind the legs. The turbulent diffusion causes a rapid redevelopment of the flow in the wake, as it is seen for the individual runs as well. Considering its boundary by the contour where the streamwise velocity attains 95% of the undisturbed value, one observes that the flow entrainment smoothens the fine details of the streamwise velocity distribution and internally to the wake, the peak velocity deficit reduces. The diffusion process causes the wake to exceed the measurement region, with consequences on the uncertainty of the drag estimate. This occurs earlier for the outdoor experiment (t* ~ 9) than for the indoor experiment (t* ~ 13), which is ascribed to the higher intensity of velocity fluctuations in the surrounding environment. The higher acquisition frequency of the outdoor experiment provides a more detailed look into the temporal development of the wake, however, at the cost of a lower accuracy and higher amount of erroneous vectors.
Next to the out-of-plane velocity, the similarity between the flow fields is also assessed by looking at the in-plane streamlines. It is apparent that the primary features are consistent throughout Fig. 9, in that, close to the cyclist, a strong downwash exists near the vertical centreline. It can be reasoned that this characteristic is responsible for the downward movement over time of the wake structure. Furthermore, a strong inwash between 0.8 m and 1.2 m from the floor is induced by the main hip vortices in both experiments, which is further increased by the head vortices as seen in Fig. 10. Over time, the hip/thigh vortex structure seems to outlast the smaller vortex structures, which in turn means that the former will dominate the wake behavior in the far wake. There, the induced inwash causes a narrowing of the upper wake, while the broadening of the lower wake structure can be assigned to the induced outwash by the vortex pair, as well as the present ground, which constrains the downwash.
The analysis of the wake in terms of vorticity elucidates some of the characteristic aspects of the flow developing around and past the cyclist. Figure 10 illustrates and compares the distribution of streamwise vortices as measured indoor (upright and time trial) and outdoor (time trial). Positive vorticity relates to counter-clockwise rotating vortices, while negative vorticity to clockwise ones. The flow structures characterizing the upright and timetrial wakes from the indoor experiment are compared in Fig. 10-left and -middle. There is substantial equivalence in the vortex structure strength and position, with the exception of the hip-thighs and the head vortices. In the former, the upright position shows higher vorticity on both sides. In the latter, the upright helmet vortices are negligible structures.
Moreover, the upright posture shows new large-scale structures, namely the shoulder vortex and the arm vortex couple. It is hypothesized that for each shoulder, one outer vortex is shed. Its generation mechanism is proper of what has been called as a 3D separation. In fact, on both sides, they are co-rotating with the hip vortices. This structure arises as a consequence of the very low pressure in the upper back of the cyclist. The arm vortex couples consist of a counterrotating vortex with respect to the shoulder vortex on the outside of the arm and a co-rotating one on the inside of the arm. They are assumed to be originated from the forward extended arms towards the brake hoods.
The vorticity field of the indoor and outdoor time-trial position exhibits an overall agreement, although some details are not exactly reproduced. This may be ascribed to the torso angle, not fully repeated during indoor and outdoor experiments. The vorticity structure presented in Fig. 10 also shows a good similarity with that reported in the studies of Crouch et al. (2014Crouch et al. ( , 2016.

Ensemble average drag area
Following the authoritative review article from Crouch et al. (2017), the drag results are presented as drag area (C d A). In fact, the overall aerodynamic efficiency of the cyclist is governed by both frontal area of the cyclists and the bike and the drag coefficient (shape of the cyclist and bike). Based on Eq. (2) and on the procedures described in Sect. 3.3, the instantaneous drag area is computed for each passage as a function of the dimensionless time. In Fig. 11, the drag area evaluation is given for five passages with the cyclist in upright posture. In the outdoor experiment, half a crank cycle is spanned along ∆t * ≈ 4, while along ∆t * ≈ 7.5 for the indoor experiment.
For t * ≤ 5 , the drag area computed via Eq.
(2) is underestimated as the contribution of the static pressure in the measurement plane is wrongly estimated too close to the cyclist due to the larger velocity gradients (the in-plane gradients are modulated due to the limited spatial resolution; the outof-plane ones are completely neglected because stereo PIV is used). In case of the outdoor experiment, a C d A plateau persists until approximately t * ≤ 10 , when a sudden drop in the drag area occurs. This can be related to part of the wake moving out of the measurement domain in several runs. In the outdoor experiment, the external atmospheric conditions and a narrower field of view cause the problem. Moreover, the outdoor experiment generally exhibits larger fluctuations, especially in the near wake, which indicate a poorer control and repeatability of experimental conditions.
The comparison between two distinct postures of the athlete is shown in Fig. 12 to illustrate the overall sensitivity of the Ring of Fire system to macroscopic variations of the drag area. Together with the time-average C d A, a shaded band wide 2 C d A represents the experimental uncertainty at 95% confidence level. Interestingly, although both the experiments were designed to obtain phase-locked average data, no clear cyclic trend depending on the crank angle is visible. This result differs from the findings of Crouch et al. (2014), who highlighted a 20% drag area variation with the crank angle, for a fixed t * . This outcome shows that the wake diffusion and turbulent mixing is the main phenomenon affecting the streamwise wake trend.

Time ensemble average drag area
The measurements of drag area value for several configurations are summarized in Fig. 13. The interval 6 ≤ t * ≤ 9 is considered, where systematic errors due to the pressure term and wake exit from the measurement region can be neglected.
The drag area of the cyclist in the outdoor experiment is higher for both time-trial and upright positions. These results are in agreement with the wake contours in Figs. 8 and 9, where a wider contour with higher peak momentum deficit is observed for the outdoor cyclist. A relative difference between 20 and 35% is measured between time-trial and upright positions, which is in agreement with literature. The   Fig. 11 Instantaneous drag area measurements with the cyclist in upright posture bigger difference between the two experiments is observed when comparing the mean drag area in upright position, with the outdoor experiment returning a higher C d A value. It is hypothesized that this is due to a bigger difference in frontal area between the upright postures compared to the difference between the time-trial postures.
Finally, the current results are compared to the data collected from literature. The results of aerodynamic research in cycling exhibit a large scatter due to differences in riders, bicycle models, postures and garment and general experimental conditions. Figure 14 (top) and (bottom) compare drag areas versus velocity, measured in time-trial position and upright position, respectively, for different experiments and some computer simulations. Measurements obtained during races are obtained at velocity between 12 and 16 m/s. In our experiments, the limited space for accelerating and braking led to a lower velocity for the tests (5-8 m/s). Conversely, Grappe (2009) showed that in the range of 5-20 m/s, the drag area of a cyclist remains approximately constant. The results from the current experiments fall within this large cloud of data and correlate favorably with wind tunnel and on-site experiments. In contrast, results from CFD simulations yield systematically lower values of drag area.

Conclusions
Large-scale stereo-PIV measurements are conducted to determine the aerodynamic drag of a moving cyclist in indoor and outdoor on-site conditions using the control volume approach. The flow is measured in the wake of a cyclist moving at 5 m/s and 8 m/s for, respectively the indoor and outdoor experiments. Instantaneous as well as ensemble average streamwise velocity fields have been obtained. Despite the differences between the two experiments in the cyclist geometry, bike model and the cycling speed, the flow fields in the near wake of the riders compare well between both experiments and literature. The instantaneous and ensemble average aerodynamic drag is evaluated via a control volume approach along the wake