Extended Data for "Aerodynamic Performance of a Baja SAE Vehicle Using Hybrid RANS–LES Approach"

applied sciences Article Flow Visualization of Spinning and Nonspinning Soccer Balls Using Computational Fluid Dynamics Takeshi Asai * , Yasumi Nakanishi , Nakaba Akiyama and Sungchan Hong Faculty of Health and Sports Sciences, University of Tsukuba, Tsukuba 305-8574, Japan; [email protected] (Y.N.); [email protected] (N.A.); [email protected] (S.H.) *Correspondence: [email protected]; Tel./Fax: +81-29-853-2711 Received: 5 June 2020; Accepted: 28 June 2020; Published: 30 June 2020   Featured Application: This study indicates that the trajectory of a spinning soccer ball is regular and stable. This result is expected to help strikers to accurately aim a curled shot or pass at the outset of the ball’s course and help goalkeepers to easily judge and quickly react to the ball’s course. Abstract: Various studies have been conducted on the aerodynamic characteristics of nonspinning and spinning soccer balls. However, the vortex structures in the wake of the balls are almost unknown. One of the main computational fluid dynamics methods used for the analysis of vortex structures is the lattice Boltzmann method as it facilitates high-precision analysis. Studies to elucidate the dominant vortex structure are important because curled shots and passes involving spinning balls are frequently used in actual soccer games. In this study, we identify the large-scale dominant vortex structure of a soccer ball and investigate the stability of the structure using the lattice Boltzmann method, wind tunnel tests, and free-flight experiments. One of the dominant vortex structures in the wake of both nonspinning and spinning balls is a large-scale counter-rotating vortex pair. The side force acting on a spinning ball stabilizes when the fluctuation of the separation points of the ball is suppressed by the rotation of the ball. Thus, although a spinning soccer ball is deflected by the Magnus effect, its trajectory is regular and stable, suggesting that a spinning ball can be aimed accurately at the outset of its course. Keywords: aerodynamics; soccer ball; vortex; drag; lift 1. Introduction Thebehavior(e.g.,thetrajectory)ofsportsballsisknowntobesignificantlyaffectedbyaerodynamic characteristics such as drag and lift [ 1 , 2 ], and various studies (including wind tunnel experiments) have been conducted to investigate the effect of these characteristics [ 3 – 9 ]. However, relatively few studies have examined the air flow around a ball and the vortex structure and dynamics that generate these aerodynamic characteristics [ 10 , 11 ], although analysis of the vortex structure and dynamics is very important to elucidate the fundamental aerodynamics of sports balls. Many studies have previously been conducted on the vortex structure of a smooth sphere [ 12 – 14 ], and the results indicate that a horseshoe vortex, an alternating vortex, or a helical vortex is formed depending on the Reynolds number [ 15 – 18 ]. However, numerous aspects regarding the formation and transition mechanisms of vortex structures are still unknown. Various studies have been conducted on the aerodynamics of soccer balls. A drag crisis occurs under certain conditions during the drag and lift measurements of nonspinning balls for which wind tunnel experiments are relatively easy to set up [ 19 ]. Wind tunnel experiments have also been Appl. Sci. 2020,10, 4543; doi:10.3390/app10134543 www.mdpi.com/journal/applsci Appl. Sci. 2020,10, 4543 2 of 14 conducted on spinning balls to measure their Magnus force and aerodynamic characteristics [ 20 ]. Furthermore, attempts have been made to measure certain aerodynamic characteristics through free-flight experiments (without using ball-supporting devices), which provide conditions similar to the actual conditions encountered by a ball [ 21 ]. However, most of these studies measured the drag and lift acting on the ball through wind tunnel experiments, and very few studies have visualized and analyzed the air flow around the ball and the dominant vortex structure that generates the drag and lift forces [ 22 ]. In particular, the air flow around nonspinning (or low-spinning) and spinning balls in free flight and the associated dominant vortex structure have not yet been examined [ 23 ]. Moreover, curled shots and passes involving spinning balls are frequently used in actual soccer ball games, and hence, studies to elucidate the dominant vortex structure are important. Computational fluid dynamics (CFD) has been conventionally used to visualize and analyze the flow around bluffbodies; in general, the vortex structure is simulated and steady-state analysis iscarried out using turbulence models such as the realizable k– ε model [ 24 , 25 ]. Recently, the lattice Boltzmann method has become one of the main CFD methods for unsteady analysis [ 26 – 28 ]. The simulation results obtained from this method are closer to the actual phenomenon than those obtained from steady-state analysis and, hence, the method facilitates the realization of high-resolution visualization and high-precision analysis. In this study, we investigated the aerodynamic characteristics and large-scale dominant vortex structure of a soccer ball using a combination of the lattice Boltzmann method, wind tunnel tests, and free-flight experiments. This combination enabled us to analyze and visualize the drag force, lift force, and side force (Magnus force) acting on nonspinning and spinning balls in free flight as well as the air flow around the balls. Our results indicate that the dominant vortex structure in the wake of a spinning soccer ball is a large-scale counter-rotating vortex pair [ 29 ], which is similar to the vortex formed around a wing tip [ 30 ], and that the side force becomes stable if the fluctuation of the separation points of the ball is suppressed by the rotation of the ball. 2. Materials and Methods 2.1. CFD Analysis Using the Lattice Boltzmann Method A three-dimensional soccer ball model was constructed (Figure 1a) from data obtained by scanning a real soccer ball (Brazuca, Adidas) using a three-dimensional laser scanner (AICON 3D, Breuckmann GmbH, Germany). For a spinning ball, the flow speed at the velocity inlet was set to 28 m/s (Reynolds number (Re) =4.25 × 10 5 ), and the spinning rates of the ball were defined as 25.1, 50.2, and 75.4 rad/s (4, 8, and 12 rps); the cases corresponding to these spinning rates were denoted as Spin A, Spin B, and Spin C, respectively in this study. For a nonspinning ball, the flow speeds at the velocity inlet were set to 8.2 m/s (Re =1.25 × 10 5 ), 19.0 m/s (Re =2.8 × 10 5 ), and 27.0 m/s (Re =4.0 × 10 5 ); the cases corresponding to these flow speeds were denoted as Nonspin A, Nonspin B, and Nonspin C, respectively. The flow speeds (ball speeds) in this experiment were based on those observed in actual soccer games [ 23 ]. A Cartesian grid was adopted to generate a spatial grid with dimensions of 20 m × 20 m × 40 m (W × H × L) comprising nearly 500 million cells (Figure 1b). A sectional grid scale technique was used in this study with a minimum scale of 1 mm and a maximum scale of 4 mm for the cases involving a spinning ball. For the cases involving a nonspinning ball, we defined minimum scales of 1.63 × 10 -4 mm at Re =1.25 × 10 5 , 7.13 × 10 -5 mm at Re =2.80 × 10 5 , and 5.02 × 10 -5 mm at Re =4.00 × 10 5 , and a maximum scale of 4 mm. This grid structure could not represent detailed vortex formations perfectly, but it was used because of computational resource constraints [ 31 ]. The outlet pressure was defined as 1013.25 hPa (i.e., atmospheric pressure). The boundary wall of the soccer ball was assumed to obey a no-slip condition, and the outer walls (including the ground surface) were defined as slip walls. The time step for the calculation was 1.018 × 10 -5 s for the cases involving a nonspinning ball and 2.037 × 10 -4 s for the cases involving a spinning ball. In this study, aerodynamic simulations were performed using the incompressible flow model of commercial CFD software (PowerFLOW 5.1, Appl. Sci. 2020,10, 4543 3 of 14 Exa Corp., USA) based on the lattice Boltzmann method [ 26 ]. The behavior of many-particle kinetic systems can be expressed in terms of the basic mechanical laws governing single-particle motions at the molecular scale. The Boltzmann equation formulates the problem in terms of a distribution function f(x,v,t), which is the number density of molecules at position xand speed vat time t[ 32 ]. The equation (in the absence of external forces) can be written as Df Dt =∂f(x,v,t) ∂t+v∇f(x,v,t)=C(x,v,t)(1) Appl. Sci. 2020, 10, x FOR PEER REVIEW 4 of 15 Here, ρ is the density of air (1.2 kg/m3), U is the flow velocity (m/s), and A is the projected area of the soccer ball (given by πR2, where R is the radius of the soccer ball). Cd, Cl, and Cs are measured in the +X, +Z, and +Y directions, respectively. The ratio of the peripheral velocity to the velocity through the air, Sp, was calculated as  =   (8) where ω is the angular velocity of the soccer ball (rad/s) and R is the radius of the soccer ball (0.11 m). One-way analysis of variance was used to statistically test the average of the Cs values for the spinning and nonspinning balls. Fast Fourier transform (FFT) analysis was employed to compare the frequency characteristics of Cs. . Figure 1. (a) Three-dimensional soccer ball model obtained by laser scanning and (b) the Cartesian grid used to generate a spatial grid (W 20 m × H 20 m × L 40 m) for computational fluid dynamics (CFD) analysis. The flow is from the left to the right. 2.2. Wind Tunnel Test A low-speed circulating wind tunnel with a six-component balance (maximum wind velocity = 55 m/s; measuring section = 1.5 m × 1.5 m; turbulence level = 0.1%) was employed (Figure 2) to verify the drag force coefficient of the nonspinning ball using CFD. The tests were conducted on a thermally bonded soccer ball with six panels (Brazuca (size 5), Adidas; officially approved for use in international games). The ball was supported at the rear using a sting (fixture) fitted to the sixcomponent wind tunnel balance. The ball was fixed to the sting using an adhesive such that it could not rotate. The ball panel orientation in the wind tunnel experiment was the same as that in the CFD analysis. The flow speeds used in the wind tunnel experiment were in the range of 7–35 m/s (Re = 1.06 × 105 to 5.31 × 105) for the nonspinning ball [19]. Figure 1. ( a ) Three-dimensional soccer ball model obtained by laser scanning and ( b ) the Cartesian grid used to generate a spatial grid (W 20 m × H 20 m × L 40 m) for computational fluid dynamics (CFD) analysis. The flow is from the left to the right. Here, the total derivative on the left-hand side represents the convective motion of particles, whereas the right-hand side expresses complex intermolecular interactions (collisions). Integration of the distribution function makes it possible to obtain macroscopic variables such as fluid density, speed, and pressure. The collision operator’s main purpose is to drive the velocity distribution function towards its equilibrium distribution. The Bhatnagar, Gross and Krook (BGK) collision operator [ 32 ] can then be defined as C(x,v,t)=−1 τ[f(x,v,t)−feq(x,v,t) ] (2) where τis the relaxation time of the fluid, and feq(x,v,t)is the equilibrium distribution function. To solve these equations efficiently, we discretized them on a three-dimensional cubic lattice using the D3Q19 model [ 32 ]. This model discretized the velocity space into 19 discrete speeds. The discrete LB equation, using a specific finite-differencing of time (∆t=1), is written as fi(x+ci∆t,t+∆t)−fi(x,t)=Ci(x,t)(3) C(x,t)=−1 τ[fi(x,t)−fieq(x,t) ] (4) Appl. Sci. 2020,10, 4543 4 of 14 A volumetric boundary scheme was chosen as the fluid-structure interaction method. Here, the particle boundary condition was conducted at the surface itself (i.e., on the facets that make up the geometry description). Each of these facets had a set of extruded parallelograms corresponding to the discrete velocity directions. For the cases involving a spinning ball, the boundary layer was simulated using a sliding mesh model [ 33 ]. Turbulence was modeled according to the very large Eddy simulation (VLES) principle [ 32 ], which directly simulates resolvable flow scales. Unresolved scales were modeled using the renormalization group form of k– ε equations with proprietary extensions to achieve VLES time-accurate physics. The lattice in this solver was composed of voxels, which are three-dimensional cubic cells. The lattice also included surfels, which are surface elements that occur in areas where the surface of a body intersects with a fluid. For the cases involving a nonspinning ball, direct numerical simulation (DNS) was employed without a turbulent model. The average drag, lift, and side forces acting on the soccer ball model were calculated from the unsteady drag, lift, and side forces over a period of 1.0 s (the calculation ran from 0.2 s to 1.2 s). The following parameters were further calculated from the CFD and experimental data collected over a range of conditions: wind velocity (U); force acting in the opposite direction of the wind (i.e., drag D); force acting perpendicular to the wind direction (i.e., lift L); and force acting sideways (S) (i.e., Magnus force) with respect to the frontal view. The aerodynamic forces determined from CFD and through the experiments were converted into the drag force coefficient (Cd), lift force coefficient (Cl), and side force coefficient (Cs) as follows: Cd =D 1 2ρU2A(5) Cl =L 1 2ρU2A(6) Cs =S 1 2ρU2A(7) Here, ρ is the density of air (1.2 kg/m 3 ), Uis the flow velocity (m/s), and Ais the projected area of the soccer ball (given by π R 2 , where Ris the radius of the soccer ball). Cd,Cl, and Cs are measured in the +X,+Z, and +Ydirections, respectively. The ratio of the peripheral velocity to the velocity through the air, Sp, was calculated as Sp =ωR U(8) where ω is the angular velocity of the soccer ball (rad/s) and Ris the radius of the soccer ball (0.11 m). One-way analysis of variance was used to statistically test the average of the Cs values for the spinning and nonspinning balls. Fast Fourier transform (FFT) analysis was employed to compare the frequency characteristics of Cs. 2.2. Wind Tunnel Test A low-speed circulating wind tunnel with a six-component balance (maximum wind velocity = 55 m/s; measuring section =1.5 m × 1.5 m; turbulence level =0.1%) was employed (Figure 2) to verify the drag force coefficient of the nonspinning ball using CFD. The tests were conducted on a thermally bonded soccer ball with six panels (Brazuca (size 5), Adidas; officially approved for use in international games). The ball was supported at the rear using a sting (fixture) fitted to the six-component wind tunnel balance. The ball was fixed to the sting using an adhesive such that it could not rotate. The ball panel orientation in the wind tunnel experiment was the same as that in the CFD analysis. The flow speeds used in the wind tunnel experiment were in the range of 7–35 m/s (Re =1.06 × 10 5 to 5.31 × 10 5 ) for the nonspinning ball [19]. Appl. Sci. 2020,10, 4543 5 of 14 Appl. Sci. 2020, 10, x FOR PEER REVIEW 5 of 15 Figure 2. Set up of wind tunnel experiment. 2.3. Free-Flight Test The flow (in the frontal cross-sectional plane) behind a soccer ball during flight was visualized using titanium tetrachloride. The soccer ball was placed directly in front of a soccer goal at a distance of 25 m, and a research participant performed a nonspinning kick with slight rotation (knuckle ball) as well as a side-spinning curled kick aimed at the goal. Both placement kicks were delivered at similar velocities to mimic the conditions in a soccer game. Two high-speed video cameras (Photron SA2, Photron Limited) were set up on one side of and behind the ball trajectory, and photographs were taken at 1000 fps. The ball speed and spin rate were measured using high-speed video images. We defined a ball with a rotation of less than 1 rps as a nonspinning ball and a ball with a rotation of more than 4 rps as a spinning ball [34]. The experimental procedure was as follows. Each soccer ball was brush painted with titanium tetrachloride, placed on a designated spot, and kicked toward the goal. As the ball moved toward the goal, the air flow around it was revealed by white smoke produced by the titanium tetrachloride. Photographs were taken using a high-speed video camera. Finally, the ball was collected and cleaned. 3. Results and Discussion 3.1. Steady Cd and Cs Obtained from CFD The steady drag force coefficient (average Cd) of a nonspinning soccer ball obtained from CFD was ~0.39 in the subcritical regime of Re = 1.25 × 105, ~0.19 in the supercritical regime of Re = 2.80 × 105, and ~0.14 at Re = 4.00 × 105 (Figure 3). The average Cd of the nonspinning ball obtained in wind tunnel tests was ~0.45 in the subcritical regime of Re = 1.24 × 105 and ~0.18 in the supercritical regime of Re = 2.82 × 105 (the critical Reynolds number was 2.67 × 105 (average Cd = 0.15)), which were close to the average Cd values obtained from CFD. We observed a peculiar spike in the subcritical regime owing to the effect of the irregular surface panel shape on the ball, which can affect turbulent reattachment. The maximum measurement error of the average Cd for verifying the measurement repeatability at the same ball orientation was less than 5% in the subcritical regime, whereas those in the critical and supercritical regimes were lower. Moreover, the maximum measurement error of Cd at different ball orientations was also less than 5% in the subcritical regime, whereas those in the critical and supercritical regimes were lower. It is considered that the entire surface roughness of the ball in this wind tunnel test was not significantly different from that of a conventional 32-panel ball (Vantaggio; the critical Reynolds number is 2.20–2.50 × 105) as the total bond length of the six-panel ball used in this wind tunnel test is similar to that of the conventional ball [35,36]. Our experimental setup could not measure the average Cd and Cs of the spinning ball. Therefore, the average Cd and Cs of the spinning ball obtained from CFD were verified using the experimental results reported in a previous study in which a soccer ball with 14 panels was used (Teamgeist (size 5); the critical Reynolds number of the ball was 2.6 × 105) [19]. The average Cd of a spinning ball obtained from CFD was ~0.28 at spin parameter (Sp) of 0.1, ~0.34 at Sp of 0.2, and ~0.36 at Sp of 0.3 (Figure 4a). The average Cd tended to increase as Sp increased. These values were slightly larger than Figure 2. Set up of wind tunnel experiment. 2.3. Free-Flight Test The flow (in the frontal cross-sectional plane) behind a soccer ball during flight was visualized using titanium tetrachloride. The soccer ball was placed directly in front of a soccer goal at a distance of 25 m, and a research participant performed a nonspinning kick with slight rotation (knuckle ball) as well as a side-spinning curled kick aimed at the goal. Both placement kicks were delivered at similar velocities to mimic the conditions in a soccer game. Two high-speed video cameras (Photron SA2, Photron Limited) were set up on one side of and behind the ball trajectory, and photographs were taken at 1000 fps. The ball speed and spin rate were measured using high-speed video images. We defined a ball with a rotation of less than 1 rps as a nonspinning ball and a ball with a rotation of more than 4 rps as a spinning ball [ 34 ]. The experimental procedure was as follows. Each soccer ball was brush painted with titanium tetrachloride, placed on a designated spot, and kicked toward the goal. As the ball moved toward the goal, the air flow around it was revealed by white smoke produced by the titanium tetrachloride. Photographs were taken using a high-speed video camera. Finally, the ball was collected and cleaned. 3. Results and Discussion 3.1. Steady Cd and Cs Obtained from CFD The steady drag force coefficient (average Cd) of a nonspinning soccer ball obtained from CFD was ~0.39 in the subcritical regime of Re =1.25 × 10 5 , ~0.19 in the supercritical regime of Re =2.80 ×105 , and ~0.14 at Re =4.00 × 10 5 (Figure 3). The average Cd of the nonspinning ball obtained in wind tunnel tests was ~0.45 in the subcritical regime of Re =1.24 ×105 and ~0.18 in the supercritical regime of Re =2.82 ×105 (the critical Reynolds number was 2.67 × 10 5 (average Cd =0.15)), which were close to the average Cd values obtained from CFD. We observed a peculiar spike in the subcritical regime owing to the effect of the irregular surface panel shape on the ball, which can affect turbulent reattachment. The maximum measurement error of the average Cd for verifying the measurement repeatability at the same ball orientation was less than 5% in the subcritical regime, whereas those in the critical and supercritical regimes were lower. Moreover, the maximum measurement error of Cd at different ball orientations was also less than 5% in the subcritical regime, whereas those in the critical and supercritical regimes were lower. It is considered that the entire surface roughness of the ball in this wind tunnel test was not significantly different from that of a conventional 32-panel ball (Vantaggio; the critical Reynolds number is 2.20–2.50 × 10 5 ) as the total bond length of the six-panel ball used in this wind tunnel test is similar to that of the conventional ball [35,36]. Our experimental setup could not measure the average Cd and Cs of the spinning ball. Therefore, the average Cd and Cs of the spinning ball obtained from CFD were verified using the experimental results reported in a previous study in which a soccer ball with 14 panels was used (Teamgeist (size 5); the critical Reynolds number of the ball was 2.6 × 10 5 ) [ 19 ]. The average Cd of a spinning ball obtained Appl. Sci. 2020,10, 4543 6 of 14 from CFD was ~0.28 at spin parameter (Sp) of 0.1, ~0.34 at Sp of 0.2, and ~0.36 at Sp of 0.3 (Figure 4a). The average Cd tended to increase as Sp increased. These values were slightly larger than the average Cd of the spinning ball obtained in wind tunnel tests [ 19 ]. The average side force coefficient (average Cs) of the spinning ball obtained from CFD was ~0.22 at Sp of 0.1, ~0.28 at Sp of 0.2, and ~0.30 at Sp of 0.3 (Figure 4b). The average Cs also tended to increase as Sp increased, and the values were similar to the average Cs of the spinning ball obtained in wind tunnel tests [19]. Appl. Sci. 2020, 10, x FOR PEER REVIEW 6 of 15 the average Cd of the spinning ball obtained in wind tunnel tests [19]. The average side force coefficient (average Cs) of the spinning ball obtained from CFD was ~0.22 at Sp of 0.1, ~0.28 at Sp of 0.2, and ~0.30 at Sp of 0.3 (Figure 4b). The average Cs also tended to increase as Sp increased, and the values were similar to the average Cs of the spinning ball obtained in wind tunnel tests [19]. Figure 3. Comparison of the drag force coefficients of a nonspinning ball determined from CFD analyses and wind tunnel tests (the error bars indicate the standard deviation of Cd). Figure 4. Comparison of (a) drag force coefficients and (b) side force coefficients of a spinning ball determined from CFD analyses and wind tunnel tests. The average Cd of the nonspinning ball obtained from CFD in this study is similar to that obtained from wind tunnel tests in the present study as well as in previous studies [19,20] and, hence, it is reasonable to conclude that the CFD calculations are valid. Further, the average Cs of the spinning Figure 3. Comparison of the drag force coefficients of a nonspinning ball determined from CFD analyses and wind tunnel tests (the error bars indicate the standard deviation of Cd). Appl. Sci. 2020, 10, x FOR PEER REVIEW 6 of 15 the average Cd of the spinning ball obtained in wind tunnel tests [19]. The average side force coefficient (average Cs) of the spinning ball obtained from CFD was ~0.22 at Sp of 0.1, ~0.28 at Sp of 0.2, and ~0.30 at Sp of 0.3 (Figure 4b). The average Cs also tended to increase as Sp increased, and the values were similar to the average Cs of the spinning ball obtained in wind tunnel tests [19]. Figure 3. Comparison of the drag force coefficients of a nonspinning ball determined from CFD analyses and wind tunnel tests (the error bars indicate the standard deviation of Cd). Figure 4. Comparison of (a) drag force coefficients and (b) side force coefficients of a spinning ball determined from CFD analyses and wind tunnel tests. The average Cd of the nonspinning ball obtained from CFD in this study is similar to that obtained from wind tunnel tests in the present study as well as in previous studies [19,20] and, hence, it is reasonable to conclude that the CFD calculations are valid. Further, the average Cs of the spinning Figure 4. Comparison of ( a ) drag force coefficients and ( b ) side force coefficients of a spinning ball determined from CFD analyses and wind tunnel tests. Appl. Sci. 2020,10, 4543 7 of 14 The average Cd of the nonspinning ball obtained from CFD in this study is similar to that obtained from wind tunnel tests in the present study as well as in previous studies [ 19 , 20 ] and, hence, it is reasonable to conclude that the CFD calculations are valid. Further, the average Cs of the spinning ball obtained from CFD is slightly larger than the value obtained in wind tunnel tests in this study, but the trends are similar in that the average Cs increases with an increase in Sp. Some reports [ 34 , 37 ] have claimed that the average Cs of the spinning ball increases as Re increases, and considering that the average Cs depends on the ball speed, these calculations may be deemed to be within the permissible range. Thus, the CFD results for the nonspinning and spinning balls in this study are in broad agreement with the results of wind tunnel tests and can be considered reliable. 3.2. Flow Visualization from CFD The boundary layer separation points (line) and pressure distribution in the case of Nonspin B of a nonspinning ball ( ω =0 rad/s, Re =2.80 × 10 5 ) fluctuated irregularly along the time axis, and the accompanying ball wake streamlines also exhibited a deflecting tendency (Figure 5a). In addition, a large-scale counter-rotating vortex pair was very often formed in the wake of the ball (Figure 5b,c), and the position of the vortex pair changed as the separation points and pressure distribution fluctuated. The other cases—Nonspin A and Nonspin C—showed similar trends. It is considered that the simulations did not accurately model the separation and reattachment phenomena, although the separation point in the case of the subcritical regime (Nonspin A) shifted to the stagnation point of the ball. Appl. Sci. 2020, 10, x FOR PEER REVIEW 7 of 15 ball obtained from CFD is slightly larger than the value obtained in wind tunnel tests in this study, but the trends are similar in that the average Cs increases with an increase in Sp. Some reports [34,37] have claimed that the average Cs of the spinning ball increases as Re increases, and considering that the average Cs depends on the ball speed, these calculations may be deemed to be within the permissible range. Thus, the CFD results for the nonspinning and spinning balls in this study are in broad agreement with the results of wind tunnel tests and can be considered reliable. 3.2. Flow Visualization from CFD The boundary layer separation points (line) and pressure distribution in the case of Nonspin B of a nonspinning ball (ω = 0 rad/s, Re = 2.80 × 105) fluctuated irregularly along the time axis, and the accompanying ball wake streamlines also exhibited a deflecting tendency (Figure 5a). In addition, a large-scale counter-rotating vortex pair was very often formed in the wake of the ball (Figure 5b,c), and the position of the vortex pair changed as the separation points and pressure distribution fluctuated. The other cases—Nonspin A and Nonspin C—showed similar trends. It is considered that the simulations did not accurately model the separation and reattachment phenomena, although the separation point in the case of the subcritical regime (Nonspin A) shifted to the stagnation point of the ball. Figure 5. Example of static pressure and x-vorticity seen on streamlines around a nonspinning soccer ball obtained using CFD from (a) side view, (b) top view, and (c) back view. On the other hand, the boundary layer separation points and pressure distribution in the case of Spin B of a spinning ball (ω = 50.2 rad/s, Re = 4.25 × 105) deflected under the influence of rotation, and Figure 5. Example of static pressure and x-vorticity seen on streamlines around a nonspinning soccer ball obtained using CFD from (a) side view, (b) top view, and (c) back view. Appl. Sci. 2020,10, 4543 8 of 14 On the other hand, the boundary layer separation points and pressure distribution in the case of Spin B of a spinning ball ( ω =50.2 rad/s, Re =4.25 × 10 5 ) deflected under the influence of rotation, and the accompanying ball wake streamlines also exhibited a deflecting tendency (Figure 6a,b). Similar to the case of a nonspinning ball, a large-scale counter-rotating vortex pair was observed in the wake, but no large fluctuations were seen, and the pair largely remained at the same position (Figure 6c). The other cases—Spin A and Spin C—exhibited similar trends. Appl. Sci. 2020, 10, x FOR PEER REVIEW 8 of 15 the accompanying ball wake streamlines also exhibited a deflecting tendency (Figure 6a,b). Similar to the case of a nonspinning ball, a large-scale counter-rotating vortex pair was observed in the wake, but no large fluctuations were seen, and the pair largely remained at the same position (Figure 6c). The other cases—Spin A and Spin C—exhibited similar trends. Figure 6. Example of static pressure and x-vorticity observed on streamlines around a spinning soccer ball using CFD from (a) side view, (b) top view, and (c) back view. The asymmetry in the flow direction of the separation points (line) in the abovementioned boundary layers of the ball affected the pressure distribution (Figure 7a,b) and vortex structure, resulting in the side and lift forces acting on the ball [23]. In addition, it is considered that the largescale counter-rotating vortex pair produced traces of deflection in the separation points and pressure distribution, and the vortex structure may be similar to that of a wing tip [38,39]. Figure 6. Example of static pressure and x-vorticity observed on streamlines around a spinning soccer ball using CFD from (a) side view, (b) top view, and (c) back view. The asymmetry in the flow direction of the separation points (line) in the abovementioned boundary layers of the ball affected the pressure distribution (Figure 7a,b) and vortex structure, resulting in the side and lift forces acting on the ball [ 23 ]. In addition, it is considered that the large-scale counter-rotating vortex pair produced traces of deflection in the separation points and pressure distribution, and the vortex structure may be similar to that of a wing tip [38,39]. In the near wake structure of the ball obtained from the back view using CFD, a large-scale counter-rotating vortex pair was formed in the wake of the nonspinning ball (Nonspin B; ω =0 rad/s, Re =2.80 × 10 5 ), and unstable movements such as rotation around the axis of travel, breakdown [ 40 , 41 ], and re-formation [ 42 ] could be observed (Figure 8a–c). Similar trends were observed in the other cases involving a nonspinning ball. In the free-flight test, a large-scale counter-rotating vortex pair was observed in the wake of the nonspinning ball ( ω =5 rad/s, Re =4.10 × 10 5 ) (Figure 8d–f)). The counter-rotating vortex pair exhibited fluctuations with highly unstable movements such as rotation around the axis of travel or breakdown, similar to the cases in the CFD analysis. Appl. Sci. 2020,10, 4543 9 of 14 Appl. Sci. 2020, 10, x FOR PEER REVIEW 9 of 15 Figure 7. Example of static pressure observed on (a) nonspinning and (b) spinning soccer balls using CFD (back view). The spin direction of the spinning soccer ball in (b) is anticlockwise from the top view. In the near wake structure of the ball obtained from the back view using CFD, a large-scale counter-rotating vortex pair was formed in the wake of the nonspinning ball (Nonspin B; ω = 0 rad/s, Re = 2.80 × 105), and unstable movements such as rotation around the axis of travel, breakdown [40,41], and re-formation [42] could be observed (Figure 8a–c). Similar trends were observed in the other cases involving a nonspinning ball. In the free-flight test, a large-scale counter-rotating vortex pair was observed in the wake of the nonspinning ball (ω = 5 rad/s, Re = 4.10 × 105) (Figure 8d–f)). The counter-rotating vortex pair exhibited fluctuations with highly unstable movements such as rotation around the axis of travel or breakdown, similar to the cases in the CFD analysis. Figure 7. Example of static pressure observed on ( a ) nonspinning and ( b ) spinning soccer balls using CFD (back view). The spin direction of the spinning soccer ball in ( b ) is anticlockwise from the top view. Appl. Sci. 2020, 10, x FOR PEER REVIEW 10 of 15 Figure 8. Flow visualization of nonspinning soccer balls using (a–c) CFD and (d–f) free-flight tests, and that of spinning soccer balls using (g–i) CFD and (j–l) free-flight tests. The ball spins with a positive rotation about the Z-axis. The Re values for the nonspinning and spinning balls are 2.80 × 105 and 4.25 × 105 (spin rate = 8 rps), respectively. In the CFD results, the color of the surface of the soccer ball varies according to the pressure coefficient, whereas the color of the streamlines varies according to the vorticity. The view is from behind the ball (back view). In the case of a spinning ball, the CFD results (Spin B; ω = 50.2 rad/s, Re = 4.25 × 105) indicated a large-scale counter-rotating vortex pair in the wake downstream of the rotating direction (right side of the back view in Figure 8), but the vortex pair remained relatively stable and exhibited only minor fluctuations (Figure 8g–i). Similar trends were observed in the other cases involving a spinning ball. In the free-flight test, a large-scale counter-rotating vortex pair, similar to that observed in the CFD analysis, was formed in the wake of the spinning ball (ω = 37.7 rad/s, Re = 3.94 × 105) on the downstream side of the rotating direction (counterclockwise when viewed from the top). Again, this vortex pair structure was relatively stable and had relatively minor fluctuations (Figure 8j–l). Various vortex structures containing a large-scale counter-rotating vortex have been observed in the wake of a smooth sphere, such as a horseshoe vortex, an alternating vortex, or a hairpin vortex Figure 8. Flow visualization of nonspinning soccer balls using ( a–c ) CFD and ( d–f ) free-flight tests, and that of spinning soccer balls using ( g–i ) CFD and ( j–l ) free-flight tests. The ball spins with a positive rotation about the Z-axis. The Re values for the nonspinning and spinning balls are 2.80 × 10 5 and 4.25 × 10 5 (spin rate =8 rps), respectively. In the CFD results, the color of the surface of the soccer ball varies according to the pressure coefficient, whereas the color of the streamlines varies according to the vorticity. The view is from behind the ball (back view).