Aerodynamics of the newly approved football for the English Premier League 2020–21 season

Footballs are typically constructed with 32 panels. Recently, the number of panels has been successively reduced to 14, 8, and 6 panels, and official balls have been adopted with complex panel shapes and aerodynamics that differ from those of 32-panel balls. The official ball for the 2020–21 season of the English Premier League comprises just four panels with a complex panel shape and surface groove design; however, its aerodynamics have not yet been clarified. This study aims to clarify the aerodynamic characteristics (drag, side force, lift force, their deviations, and critical Reynolds number) of the new 4-panel ball (Flight 2020, Nike) in comparison to a 6-panel ball (Tsubasa 2020, Adidas) and conventional 32-panel ball (Pelada 2020, Molten) using a wind tunnel test, surface design measurement, and a simple 2D flight simulation. The results showed that Flight 2020 has greater surface roughness and smaller critical Reynolds number than Pelada 2020 and Tsubasa 2020, resulting to its marginally greater drag force in the supercritical region, and slightly smaller fluctuations of the side and lift forces. Furthermore, Flight with a symmetrical orientation exhibits a significantly higher drag coefficient in the supercritical region, suggesting its greater air resistance during flight under this condition.


Scientific Reports
| (2021) 11:9578 | https://doi.org/10.1038/s41598-021-89162-y www.nature.com/scientificreports/ transition velocity 28,29 , the turbulence intensity (< 0.1%) of this experiment is presumed to be slightly smaller than that of an actual football stadium 30,31 . However, to obtain the basic aerodynamic characteristics of the ball, the experiment was conducted using a wind tunnel setup. The wind speed was incremented from 7 m/s (Re ≈ 1.0 × 10 5 ) to 35 m/s (Re ≈ 5.0 × 10 5 ) in intervals of 1 m/s every 10 s (100 Hz). The main flow speed was continuously controlled using a computer, and the aerodynamic force applied to the ball was measured for 10 s after the main flow reached the set flow speed and became sufficiently stable. The forces acting on the footballs were measured with a sting-type six-component force measurement balance (LMC-61256, Nissho Electric Works). The measured aerodynamic drag force (D), side force (S), and lift force (L) were converted to the drag force coefficient (C d ), side force coefficient (Cs), and lift force coefficient (Cl), respectively: where ρ is the air density (ρ = 1.2 kg/m 3 ), U is the flow speed, and A is the projected area of the football (A = π × 0.11 2 = 0.038 m 2 ).
The back of the ball was supported by a sting with a length of 0.8 m and width of 0.02 m in the wind tunnel. Since placing the six-component force measurement balance behind the football reduced the effect of the sting vibrations to a negligible value, the magnitude of the sting forces was neglected.
Three types of balls were used in this experiment: the new 4-panel ball (Flight 2020, Nike), a 6-panel ball (Tsubasa 2020, Adidas), and a conventional 32-panel ball (Pelada 2020, Molten). Tsubasa 2020 was the official ball of the 2020 FIFA U-20 Women's World Cup and Pelada 2020 is a conventional ball approved by FIFA. We defined two faces for the ball orientation: (A) symmetric (about the horizontal axis) and (B) asymmetric (rotated 30° horizontally, clockwise from the top view). We then measured the aerodynamic forces acting on each ball face. The symmetric orientation cases were designated as Pelada A, Tsubasa A, and Flight A, and the asymmetric orientation cases were designated as Pelada B, Tsubasa B, and Flight B. In order to examine the accuracy of the measurement balance and validity of the experiment, we also measured the drag coefficient of a smooth sphere made of plastic (radius is 0.11 m).
Tested footballs. For the statistical analysis, we employed a multiple comparison test (Tukey-Kramer method) to compare the average drag coefficients and average standard deviation of the side and lift coefficients in the supercritical region.
Panel shape and surface groove design measurement. To study the characteristics of the football surface shape in more detail, the surface shape parameters that affect the surface roughness were measured. The parameters obtained in this study were the length of the panel seams (bonds and grooves) and the width and depth of the grooves. The length of the panel joints was measured using a curvimeter (Concurve 10; Koizumi Sokki Mfg. Co., Ltd., Japan). The width and depth of the grooves were measured using a high-speed 2D laser scanner (LJ-V7000; Keyence Corp., Japan). To measure these parameters, all the grooves of the football were Ball trajectory simulation. We conducted a simple 2D flight simulation to compare the effects of the drag coefficients and orientations on the flight distance and trajectory of the different balls in the subcritical regime (U = 15 m/s) and supercritical regime (U = 30 m/s). For all cases, the initial flight trajectory was set to an angle of 30°1 8 . The following equations were used for the simulation: where m is the mass of the football, a h is the horizontal acceleration of the ball, a v is the vertical acceleration of the ball, g is the gravitational acceleration, and γ is the initial attack angle of the ball flight trajectory. The 2D ball velocity and displacement were obtained by the explicit Euler method by interpolating the drag coefficient at each time step (0.001 s) using a sixth-degree polynomial approximation 23,25 . The lift and side forces acting on the ball were omitted in this ball trajectory simulation.
We employed a multiple comparison test (Tukey-Kramer method) to compare the average standard deviations of the side and lift forces at the initial ball (flow) speeds of 15 and 30 m/s (Re of ~ 2.1 × 10 5 and ~ 4.2 × 10 5 , respectively). Extended total length of the panel seams (bonds and grooves) on the ball surface. The extended total length of the panel seams (bonds and grooves) on the ball surface are ~ 10,910 mm for Flight 2020, which is longer than that of the other balls (Table 1). Flight 2020 has a groove depth of 0.7 mm, which is similar to that of the other balls, and a width of 5.7 mm, which is substantially wider than that of the other balls.

Results
The correlation coefficients of the critical Reynolds number to the extended total length, groove width, and groove depth are 0.81 (p = 0.052), 0.87 (p = 0.023), and 0.20 (p = 0.702) respectively in the symmetric and asymmetric orientation cases of the three types of balls.

Ball trajectory simulation. For the symmetric orientation cases at an initial ball speed of 15 m/s, Tsubasa
A has the longest flight distance at 17.6 m, followed by Flight A at 17.1 m and Pelada A at 16.0 m (Fig. 4). Meanwhile for the asymmetric orientation cases at an initial ball speed of 15 m/s, Tsubasa B has the longest flight distance at 17.  (Fig. 5). Similar tendencies are demonstrated at an initial ball speed of 30 m/s (Fig. 6). The significance of the difference (P < 0.05) was determined through a multiple comparison test. The standard deviations of the side force (S) and lift force (L) tend to increase with speed for all balls and orientations (Fig. 7). The standard deviations of the side and lift forces are slightly smaller for Flight 2020 than for Pelada 2020. Moreover, the standard deviations of the side and lift forces are considerably smaller for Flight B than for Pelada B and Tsubasa B. In

Discussion
In the wind tunnel test, the smooth sphere has a critical Reynolds number of approximately 4.0 × 10 5 , which is in good agreement with the drag coefficient curve of Achenbach's smooth sphere 16 (Fig. 3).    www.nature.com/scientificreports/ In both the symmetric and asymmetric orientation cases, Flight 2020 has the lowest critical Reynolds number, followed by Tsubasa 2020 and then Pelada 2020. The critical Reynolds number of a spherical bluff body tends to decrease as the surface roughness increases 16 . The total length of the panel bonds and grooves, groove width, and groove depth on the surface has been suggested as an index of the surface roughness of a ball 23,24 . Notably, Flight 2020 has the widest groove width and longest total length, followed by Tsubasa 2020 and then Pelada 2020. These results indicate that Flight 2020 has the highest surface roughness, followed by Tsubasa 2020 and then Pelada 2020.
In both the symmetric and asymmetric orientation cases at the initial ball speed of 15 m/s (Re ≈ 2.10 × 10 5 ), Flight 2020 has the longest flight distance than Pelada 2020. This is can be due to the initial ball speed, which is in the supercritical region for Flight 2020, where C d is relatively low, whereas it is in the critical region for Pelada 2020, where C d is increasing. In contrast, at an initial ball speed of 30 m/s (Re ≈ 4.20 × 10 5 ), Flight 2020 has a shorter flight distance than Pelada 2020 in both the symmetric and asymmetric orientation cases. This may be attributed to the initial ball speed that is in the supercritical region for both balls. Hence, Flight 2020 has a shorter flight distance because it has a larger C d than Pelada 2020.
The drag coefficients extracted by tracking a football in flight using high-speed video cameras and the subsequent trajectory analysis 18,33 approximately confirmed the present drag coefficient curves obtained in the study. Flight A, which has a symmetrical ball face orientation, exhibits a substantially larger drag coefficient in the supercritical region compared to the other balls. Even in the ball trajectory simulation, the flight distance of Flight A is approximately 10% less than that of the other balls. Further, in the supercritical region, the average drag coefficient of Flight 2020 tends to be higher than that of the other balls (p < 0.001). Therefore, during its actual flight, where there is no rotation and the fluctuation of the stagnation point is extremely small, unlike in the wind tunnel experiment, it is suggested that Flight A could fly with higher air resistance than the other balls. However, this condition rarely occurs.
Based on the results, we concluded that Flight 2020 has a smaller critical Reynolds number and greater surface roughness than Pelada 2020, resulting in its slightly higher air resistance in the supercritical region. This may be attributed to the higher excitation occurring at the separation of the boundary layer on the surface of Flight 2020 owing to the presence of the panel bonds and grooves. Further, Flight 2020 exhibits smaller fluctuations (i.e. standard deviations) of its C s and C l values at initial ball speeds of 15 and 30 m/s than those of Pelada 2020 for both orientations. In addition, Pelada 2020 tends to have larger deviations of S and L than those of Flight 2020 at the different ball speeds. The average deviations of the S and L of Flight 2020 are lower than that of the other balls in supercritical region (p < 0.001), which can be partially attributed to its higher surface roughness 23 . The deviations of S and L are inferred to be related to the instability of the flight trajectory without spinning or minimal spinning 21,22 . Therefore, Flight 2020 tends to suppress the instability of flight trajectory more than Pelada 2020 owing to its smaller deviations of S and L. In volleyball, Wei et al. 9 found that the deviations of S and L are greater with the asymmetric orientation than the symmetric orientation. However, our results showed the opposite; the deviations tend to be smaller with the asymmetric orientation than with the symmetric orientation. This difference could be due to measurement region limited to the vicinity of the critical Reynolds number in the study of Wei et al. 9 , whereas a wider region from the subcritical to supercritical Reynolds numbers is considered in our study. Further, the different ball shapes are also another factor. In summary, Flight 2020 has a greater surface roughness and smaller critical Reynolds number than Pelada 2020, resulting to its marginally greater drag force in the supercritical region, and slightly smaller fluctuations of the side and lift forces.
Owing to the smaller critical Reynolds number of Flight 2020 than those of Tsubasa 2020 and Pelada 2020, it is expected to be travel slightly more for short passes (represented by the initial ball speed of approximately 15 m/s). However, because Flight 2020 has a large air resistance in the supercritical region, it should carry slightly less for strong shots, long free kicks, and other situations that required a longer distance (represented by the ball speed of approximately 30 m/s). Since Flight 2020 has smaller fluctuations in the side and lift forces than Tsubasa 2020 and Pelada 2020, it is expected to experience smaller irregular changes during the ball trajectory, thereby leading to its stability during flight.
This study has several limitations. First, we measured the aerodynamics of a non-spinning football fixed in a wind tunnel and did not take measurements using a spinning ball. However, in actual applications, the ball inevitably spins. Hence, the aerodynamics of a spinning ball needs to be considered. Second, the aerodynamic forces acting on the ball were measured, but the vortex structure around the ball and dynamics of the boundary layer that cause these forces have not been visualised or studied. These should also be examined through visualisation techniques, such as particle image velocimetry.

Conclusions
This study aimed to clarify the aerodynamic characteristics (drag, side force, lift force, their deviations, and critical Reynolds number) of the new 4-panel ball (Flight 2020, Nike) in comparison to a 6-panel ball (Tsubasa 2020, Adidas) and conventional 32-panel ball (Pelada 2020, Molten) using a wind tunnel test, surface design measurement, and a simple 2D flight simulation. Flight 2020 had a greater surface roughness and smaller critical Reynolds number than Pelada 2020 and Tsubasa 2020, which led to a marginally greater drag force in the supercritical region and slightly smaller fluctuations of the side and lift forces. In particular, the symmetrical orientation of Flight 2020 demonstrated a considerably higher drag coefficient in the supercritical region than the other studied balls, suggesting it could experience greater air resistance during flight under this condition. On the other hand, the smaller fluctuations in the side and lift forces of Flight 2020 than those of Pelada and Tsubasa 2020 indicates that it is more likely to experience smaller irregular changes in the ball trajectory, thereby possibly leading to its stability during flight.