**ABSTRACT **

**Key words:**Training, human body, drag, tandem, finite element modeling.

**INTRODUCTION **

The few experimental studies conducted in swimmers and triathlon athletes showed that the distance between swimmers significantly influences the energy cost of the swimmer submitted to the suction effect (Bassett et al., 1991; Chatard and Wilson, 2003; Hausswirth et al., 1999; 2001) and it helps proper technique maintenance when fatigue appears (Chollet et al., 2000). On the other hand, it also helps sparing energetic resources that can be used as an advantage in later phases of the competition, as it happens in triathlon where the athletes pass from the swimming to the cycling course (Delextrat et al., 2003) and afterwards to the race course (Hausswirth et al., 2001). In competitive swimmers, the wake generated by the leading swimmer induces significant reductions in energy cost (from a mean value of 3.12 ± 0.66 to 2.85 ± 0.63 l.O2.min-1), in blood lactate concentration (from a mean value of 5.0 ± 0.5 to 3.4 ± 0.6 mmol.l-1) and in perceived exertion ratings (from a mean value of 14.9 ± 0.5 to 11.7 ± 0.4) in the back swimmer (Bassett et al. , 1991).

The typical approach to study the dynamical interaction between bodies moving in a queue in a fluid is to experimentally investigate the forces generated as a function of the distance between two or more bodies. An alternative option to deal with this problem is to apply numerical simulation techniques to determine the forces exerted by the fluid on moving bodies. The Computational Fluid Dynamics (CFD) models started to be used in the middle of the 90s in the study of insects and birds during flight (Liu et al., 1995; 1997), as well as in the computation of the aerodynamic and hydrodynamic forces involved in the propulsion of animals moving through body undulation (Cheng and Chahine, 2001; Liu et al., 1996). Recently, medical applications were also described by this method, analyzing the fluid flow inside the human body (Berthier et al., 2002; Marshall et al., 2004).

In swimming research, the numerical simulation started to solve problems that were expensive and time-spending. The pioneer study of Bixler and Schloder, 1996 with CFD tools was performed in conditions of steady flow and unsteady conditions (with acceleration) in an external surface of a disc with an area similar to a typical swimmer hand. The same authors revisited this topic (Bixler and Riewald, 2002) and forwarded to the calculation of the drag and lift coefficients through the application of a two dimensional CFD hand/forearm analysis. These simulations confirmed the importance of the CFD in swimming research. Nevertheless few publications have addressed this issue, either studying the production of propulsive forces (Alves et al., 2007; Rouboa et al., 2006; Sato and Hino, 2003; Silva et al., 2005) or the drag suffered by the swimmers' body (Lyttle and Keys, 2004).

Therefore, the aims of the present study were: (i) to describe the pressure profile around the swimmers in drafting through CFD analysis; (ii) to assess the drafting distances effect in the drag coefficient; (iii) to assess the distance in which the drag coefficient of the back swimmer is equal to the drag coefficient of the leading swimmer, considering various flow velocities.

**METHODS **

Cd = FD / (0.5.?.S.V2) (1)

Where FD is the drag force, ? is the fluid density, S is the projection surface of the model and V is the flow velocity.

Computational Fluid Dynamics methodology consists in a mathematical model applied to the fluid flow. In a given computational domain, it implements the Navier-Stokes equations with discretized algebraic expressions and solves them by iterative calculations. This computational domain consists in a two dimensional grid or mesh of cells that simulate the fluid flow. The fluid mechanical properties, the flow characteristics along the outside grid boundaries and the mathematical relationship to account for the turbulence were considered.

**Mathematical model**

It is usually stated that a system of equations representing turbulent models can not predict the details of an unsteady flow, such as the flow around structures, even when including a low Reynolds number (Re). To overcome this difficulty, the incompressible Reynolds averaged Navier-Stokes equations with the standard k-epsilon (k-?) model were considered and implemented in the commercial code Fluent® (Moreira et al., 2006). Accordingly, the continuity equation, the momentum conservation equations and the turbulent and dissipated energy (k-?) conservation equations (for an incompressible fluid in Cartesian coordinates) were written in a conservative form as:

divV = 0 (2)

(3)

(4)

(5)

Where k is the turbulent kinetic energy and ? is the turbulent kinetic energy dissipation ratio. x and y represent the horizontal and vertical coordinates, respectively, Vx and Vy represent the x and y components of the velocity V. ? is the turbulent viscosity, ? represents the fluid density, v is the kinematic viscosity and Ф is the pressure strain. C2, C?, ?? and ?K are model constants, 1.92, 0.09, 1.30 and 1.00, respectively.

**Resolution method**

The numerical method used by Fluent® was based on the finite elements approach. The steady solutions of the governing system equations were given in each square element of the discretized full domain. To solve the linear system, Fluent® code adopted an AMG (Algebraic Multi-Grid) solver. Horizontal and vertical velocities, pressure (p), turbulent kinetic energy (k) and turbulent kinetic energy dissipation ratio (?) were considered as degrees of freedom (DOF) for each grid element.

The convergence criterion of AMG was set at 10-3 for the relative error between successive iterations and for all the variables (Vx, Vy, p, k and ?).

The numerical simulation was carried out in two dimensions (2D) for the computational domain in steady flow conditions.

**Application**

The software Fluent® needs a first approach for treatment of the geometry and conformation of the mesh for the model, requiring the use of the Fluent® pre-processor Gambit. The accomplishment of the model in two dimensions in the Gambit allowed to generate meshes and to define the necessary boundary conditions for the application in Fluent®.

**Boundary conditions**

The boundary conditions are those presented in Figure 1. In the upper and lower walls the velocity value in the two coordinates was considered null. The lower wall represents the bottom of the swimming pool and the upper wall represents the top of the pool. In the interface between air and water the fluid velocity is considered equal to the swimmer's absolute velocity of displacement. In the entry (left), the velocity assumed, on the y coordinate (perpendicular to the model), a null value and, on the x coordinate (in the model direction), the value ranged from 1.60 m·s-1 to 2.00 m·s-1, with increments of 0.10 m·s-1. In the exit (right) the relative pressure (difference between the absolute pressure or the pressure of the fluid and the environmental pressure, P-Pa) was zero to assure that there was not any flow conflict by reflection. The horizontal distance between swimmers changed from 0.50 m to 8.00 m. This distance was measured from the toe tips of the leading swimmer to the head vertex of the back swimmer. The position of the leading swimmer remained constant during all simulations. The model representing the swimmers was 1.80 m height and 0.5 m width, with the arms aside the trunk. These anatomical characteristics were chosen to be within the anthropometrical data of a group of elite national level swimmers.

Under these conditions, the Cd values of both swimmers were calculated for different flow velocities (1.60 m·s-1 < V < 2.00 m·s-1) and different distances between swimmers (0.50 m < x < 8. 00 m). The Cd of the leading swimmer was considered constant for all distances between swimmers.

The Cd values were calculated by integration of the relative pressure through the body surface using equation 6:

Cd = [?s(P-Pa)dS] / (0.5.?.S.V2) (6)

In the CFD application the mesh with a spacing of 0.05 generated 12812 cells (Figure 2). The computation time was about 3 min for each simulation, changing the number of iterations from 100 to 130.

**Statistics**

To estimate the Cd of the back swimmer and compare it with the assumed Cd of the leading swimmer, a curve fitting based on the application of various polynomials was performed. The evaluation of the quality of the adjustment and the selection of the polynomial degree were done aiming the lower prediction error (dispersion of the residues around the adjustment).

**RESULTS **

The Cd of the leading swimmer (on left, in Figure 1) was considered constant for all the distances between swimmers and for the same flow velocity. However, as the flow velocity increased the Cd of the leading swimmer decreased (see Table 1).

The CFD approach also showed that, for a 0.50 m distance, the back swimmer's Cd corresponds to a mean value of 55.97 ± 1.03% of that of the leading swimmer. This value progressively increased until the distance between swimmers reaches 6.0 m, where the Cd of the back swimmer corresponded to a mean value of 84.38 ± 0.71% of that by the leading swimmer (see Figure 4).

However, it was possible to verify that for all the velocities and distances higher than 6 m (until 8 m), the back swimmer's Cd remained almost constant and did not tend to the value presented by the leading swimmer. To calculate the distances in which the back swimmer's Cd equals the value of the leading swimmer, a fitting of the Cd curves of the back swimmer was carried out (according to a polynomial function of the values found until the 6 m distance). Table 2 shows the polynomials retained by this model.

The results indicated that the Cd of the back swimmer (considering the various flow velocities) was equal to that by the leading swimmer in distances ranging between 6.45 m and 8.90 m (see Table 1).

**DISCUSSION **

The major finding of this work was that the Cd of the back swimmer progressively increases concomitantly with the distance between swimmers, until a distance between 6.45 m and 8.90 m is reached. Afterwards the Cd of both swimmers remains equal.

In the present study we only analyzed a passive drag situation. Therefore, we attempted to compare our numerical results with experimental data found in passive drag analyses.

Comparing the values of the Cd obtained by the CFD analysis with experimental data from other studies, we found values that agree with the previous research (Clarys, 1979; Mollendorf et al., 2004). Clarys, 1979 reported Cd values between 0.58 and 1.04 in the human body in passive towing with a Reynolds number varying between 6.6 x 105 and 3.9 x 106. Recently, Mollendorf et al., 2004 has found Cd mean values in passive towing between 0.83 ± 0.04 and 0.90 ± 0.08 in a group of male adult swimmers using various types of competition suits.

The comparison between our data and data from active drag measurements is not simple, since there are different approaches to evaluate active drag, which can lead to different results (Toussaint et al., 2004). On the other hand, it is difficult to experimentally evaluate active drag with the swimmer totally submerged, which is what happens in our study. In the future, we aim to improve numerical simulations, allowing, for instance, to evaluate underwater active drag while the swimmer is kicking. With these data, we will be able to compare passive and active drag using the same methodology.

The Cd behavior verified in the present study with different flow velocities showed that when velocity increases, Cd decreases, which is consistent with previous observations with non-human models. In lifeless objects with regular geometry (not biological) the Cd varies with the inverse rate of Re (Wu, 1971). The same author presented Cd data of dead fishes (equivalent to passive towing) in a moving water tank. It was verified that a gradual reduction of this variable occurred when the Re increased. In that study the slope of the regression line was about -0.40, changing the Cd values between 0.1 and 0.5. Jiskoot and Clarys, 1975 described, in the human body when completely submerged at a 0.60 m depth and in passive towing (velocity range between 1.50 m.s-1 and 1.90 m·s-1), an inverse relationship between the Cd and the velocity values calculated by drag (D) and the mean body dimensions of the sample reported by Clarys, 1978. The Cd values changed between the unit and 0.95, with a slope in the regression with the velocity of -0.17. This value is close to the value (-0.16) found by Lyttle et al., 1999, when they calculated the Cd by the D values reported for swimmers in passive towing with an equal depth of 0.6 m. In that study, the range of velocities varied between 1.60 m·s-1 and 3.10 m·s-1. Therefore, the inverse association between the Cd and the velocity of the flow that we have found seems to correspond to what occurs in an experimental situation. Though, few studies have addressed the human body completely submerged, in passive towing. However, the values presented in the present study showed a gradual reduction of Cd with velocity increase and presented, also, a higher slope (-3.96). This slope is similar to experimental values found in lifeless objects with a regular geometry or in fishes with bodies of higher hydrodynamic profile. Therefore, it is possible that the human body representation (model) that was used is still much stylized. On the other hand, the use of a three dimensional model with the arms extended at the front, in the streamlined position, could lead to improve CFD data. This is another issue that should be addressed in future studies.

Comparing our Cd values observed for different distances between swimmers with the experimental studies, a general agreement can be found. Chatard and Wilson, 2003 studied the effect of the distance between swimmers displacing in a queue on the metabolic and hydrodynamic responses of the swimmer displacing immediately behind. The distances that were evaluated were 0.0, 0.5, 1.0 and 1.5 m. The drag was explored by a passive towing situation. The authors found a decrease of 20% to 21% in passive drag in distances between 0 m and 0.5 m and concluded that these distances were the most advantageous to the swimmer submitted to the suction effect. In the present study, as the distance between swimmers increased, the Cd value of the back swimmer approaches that of the leading swimmer. To a minimal distance between swimmers (0.50 m), the Cd mean value of the back swimmer is only 55.97 ± 1.03% of the value of the leading swimmer. Thus, the increase of this distance may induce an increase in the effort of the back swimmer. When the distance reached a value of 6.00 m, the Cd percentage difference between both swimmers increased to values from 83.77% up to 85.45%, depending on the velocity, as it can be verified in Figure 4.

Considering only the calculations of the drafting effect based on the application of numerical simulations we concluded that Cd values of both swimmers do not equalize. This may be because it is a free flow and/or due to the fact that we did not consider other factors that interact with the flow (i.e. the swimmers roughness and the distance to the wall - right side of the domain). This latter situation, particularly, could have an important effect on Cd. Indeed, the approach of the back swimmer to the wall of the right side could have prevented the velocity convergence and the pressure reestablishment. This fact leaded us to remove the values higher than 6 m from the model, as those values could influence the natural convergence observed until the 6.00 m distance between swimmers. Therefore, in future studies, this issue should be considered and a domain with larger dimensions may be used to enable pressure reestablishment and velocity convergence. Nevertheless, it was interesting to observe that our data revealed that the distance in which the swimmers were in the same hydrodynamics conditions ranged from 6.45 m to 8.90 m.

This study is also limited to the analysis of entirely submerged bodies. In fact, this is a different situation from that observed at the water surface. To our knowledge, up to the present, no work was able to analyze the human body in the interface water-air. We chose to analyze the body totally submerged in the water as a first approach to study the human movement in swimming, as this approach simplifies the simulation conditions. Moreover, we simulated a domain with a high depth to minimize the wave drag and to include the water flow under and above the model. Future models will be applied for bodies on/at the water surface, taking into account the above and under water body volumes and fluid characteristics.

**CONCLUSION **

With the present study we confirmed the importance of the CFD methodology of numerical simulation to enable the analysis of the fluid forces around objects in water, as it happens in swimming.

**ACKNOWLEDGMENTS **

**KEY POINTS **

The relative drag coefficient of the back swimmer was least (about 56% of the leading swimmer) for the smallest inter-swimmer distance (0.5 m).

The drag coefficient values of both swimmers in drafting were equal to distances ranging between 6.45 m and 8.90 m, considering the different flow velocities.

The numerical simulation techniques could be a good approach to enable the analysis of the fluid forces around objects in water, as it happens in swimming.

**AUTHORS BIOGRAPHY **

Antonio Jose SILVA

Employment: Associated and Aggregated Professor at the Sport Sciences Department of the University of Tras-os-Montes and Alto Douro (UTAD - Portugal).

Degree: PhD.

Research interests: Physiological and biomechanical indicators of energy cost during physical activities, namely in swimming.

E-mail: ajsilva@utad.pt

Abel ROUBOA

Employment: Associated Professor at the Engineering Department of the University of Tras-os-Montes and Alto Douro, Portugal.

Degree: PhD.

Research interests: Applied Mathematics, Parallel Computation and Biomechanics, namely, applied on Sport and Biomedical domains.

E-mail: rouboa@utad.pt

Antonio MOREIRA

Employment: Assistant Professor at the Sports Science Institute of Rio Maior (Portugal).

Degree: PhD.

Research interests: The application of numerical simulation techniques on sports domain.

E-mail: antmoreira@esdrm.pt

Victor Machado REIS

Employment: Assistant Professor at the Sport Sciences Department of the University of Tras-os-Montes and Alto Douro, Portugal.

Degree: PhD.

Research interests: The physiological and biomechanical indicators of energy cost during physical activities.

E-mail: vreis@utad.pt

Francisco ALVES

Employment: Associated Professor at the Faculty of Human Kinetics of the Technical University of Lisbon (FMH-UTL, Portugal).

Degree: PhD.

Research interests: The biomechanical and physiological determinants of swimming economy and the evaluation of the training process.

E-mail: falves@fmh.utl.pt

Joao Paulo VILAS-BOAS

Employment: Full Professor at the Faculty of Sport of the University of Porto (Portugal).

Degree: PhD.

Research interests: Swimming Biomechanics and Bioenergetics.

E-mail: jpvb@fade.up.pt

Daniel Almeida MARINHO

Employment: PhD student of the University of Tras-os-Montes and Alto Douro, Portugal.

Degree: MSc.

Research interests: The biomechanical and physiological determinant factors of the sports performance, specially the swimming performance.

E-mail: dmarinho@utad.pt