Computational Fluid Dynamics Study of Swimmer's Hand Velocity, Orientation, and Shape: Contributions to Hydrodynamics

The aim of this paper is to determine the hydrodynamic characteristics of swimmer's scanned hand models for various combinations of both the angle of attack and the sweepback angle and shape and velocity of swimmer's hand, simulating separate underwater arm stroke phases of freestyle (front crawl) swimming. Four realistic 3D models of swimmer's hand corresponding to different combinations of separated/closed fingers positions were used to simulate different underwater front crawl phases. The fluid flow was simulated using FLUENT (ANSYS, PA, USA). Drag force and drag coefficient were calculated using (computational fluid dynamics) CFD in steady state. Results showed that the drag force and coefficient varied at the different flow velocities on all shapes of the hand and variation was observed for different hand positions corresponding to different stroke phases. The models of the hand with thumb adducted and abducted generated the highest drag forces and drag coefficients. The current study suggests that the realistic variation of both the orientation angles influenced higher values of drag, lift, and resultant coefficients and forces. To augment resultant force, which affects swimmer's propulsion, the swimmer should concentrate in effectively optimising achievable hand areas during crucial propulsive phases.


Introduction
The swimming propulsion is the result of the interaction of applied forces with water and is predominantly attributed to muscular force applied by hands and forearms. It was emphasised that the major part of about 85% to 90% of propulsion generation in water is created by the application of force by an arm [1,2]. It was suggested that, the shoulder has the least propulsive potential, the forearm and hand approximately equal according to linear and angular velocities which are least by the shoulder and greater by the forearm and hand [3]. Miller [3] used the mathematical model of the front crawl arm pull with bending elbow and showed that the ratio of hydrodynamic forces of hand was about 2.5 bigger when compared with forearm. Toussaint and Truijens [4] studied the visualization of flow tufts around arm and hand. They showed that a strong pressure gradient along the arm occurred that induced axial flow directed from elbow to the hand. Apart from that, the biggest influence of pressure (relative to atmospheric pressure) was identified on the palm of hand when swimming at sprint speed. The authors noted that the pressure was not corrected for differences in hydrostatic pressure due to differences in the depth of the sensors.
In reality, swimmers can change the depth, orientation, shape, and velocity of hand throughout underwater front crawl cycle. All these parameters have direct influence on propulsion force. The prevailing theory of propulsion generation relates to Newton's second and third laws of motion where propulsion is the vector sum of drag ( : force opposite to the direction of fluid flow) and lift ( : forces that are perpendicular to the fluid flow) forces [4,5]. These components depend on the density of the fluid, the velocity of the limb relative to the fluid, the projected surface area of the limb, and the coefficient of drag ( : a dimensionless constant used to show the resistance of the object in a fluid environment) and lift ( ) which vary according to the shape of the limb and its orientation (i.e., the angle of attack (or pitch angle) and the sweepback angle) (Figure 1). To uncover an influence of various parameters influencing propulsive force or its components, experimental studies [6][7][8] as well as numerical studies were carried out in the past [5,[9][10][11][12].
One of the approachs to calculate and assess forces acting on swimmer or on separated body segments is by an application of computational fluid dynamics (CFD) method. This methodology allows the analysis of the water flow with a reduced amount of complexity and is economical alternative to an experimental method. The validity studies of CFD results are usually carried out with comparative experimental studies. Gardano and Dabnichki [7] used replica of the entire human arm to assess drag and lift coefficients, measured in a low speed wind tunnel based experiment. The computational hand model was created from human body similar to experimental arm model. Authors showed that lift and drag coefficients obtained in wind tunnel experiments were in good correlation with CFD obtained from CFD simulations carried out using FLUENT (ANSYS, PA, USA). In similar lines, Bixler and Riewald [5] simulated the steady flow around a swimmer's hand and arm at various angles of attack in steady state by using CFD method. They improved validity of calculation and demonstrate that the force coefficients computed for the hand and arm compared well with steadystate coefficients determined experimentally by the previous studies.
In consideration of significant improvement in the methodology and validity of the CFD studies, this method is applied more often in swimming investigations. The calculations through CFD method were carried out to analyse the propulsive forces produced by the propelling segments and the drag force resisting forward motion. Rouboa et al. [11] calculated drag and lift coefficients and drag force for steady flow around a 2D swimmer's hand/forearm model placed at different angles of attack with application CFD technique. These results were compared with previous CFD studies and experimental results, which illustrate similar values under the steady-state flow conditions. It was explained that the increase in flow velocity did not have much influence on the variation of drag coefficient when the model of hand/forearm simulated is positioned at the same pitch angle. This was followed by the steady-state CFD analysis of the hydrodynamic characteristics of a true swimmer's hand model with the thumb in different positions, while the other fingers are kept close together [9]. Drag and lift coefficients were calculated for different angles of attack (the sweepback angle was equal to zero). The combination of drag and lift coefficients (resultant force coefficient) showed that the hand model positioned with the thumb fully abducted presented higher values than the other positions with the thumb partially abducted and adducted at angles of attack of 0 ∘ and 45 ∘ . However, at an angle of attack of 90 ∘ , the position with the thumb adducted presented the highest value of resultant force coefficient. Also, the lift and drag coefficients were steady despite the increase in the flow velocity, and they varied according to changes in pitch angle. Minetti et al. [12] used CFD calculation and in their short communication hypothesized that an intermediate finger spacing in the 3D hand model could increase a higher coefficient of drag providing swimmers with additional thrust. It was indicated that an optimal finger spacing (12 ∘ ), roughly corresponding to the resting hand posture increases, the drag coefficient (+8.8%), which is functionally equivalent to a greater hand palm area. In the following year, Marinho et al. [10] confirmed these findings using CFD method in their study of the effect of finger spread on the propulsive force production in swimming. The steady-state CFD and 3D scanned models of the hand at the different angles of attack (the sweepback angle equal zero) were used. These results confirmed that the model with a small spread between fingers presented higher values of drag coefficient than did the models with fingers closed and fingers with a large spread. And it was concluded that the optimum finger spread could allow the hand to create more propulsive force during swimming.
Researchers proved that the hand orientation as well as hand shape contributes to the values of hydrodynamics parameters. Rouboa et al. [11] agree that both the propulsion and drag forces will fluctuate during individual phases of stroke cycle along with the respective variation of angle of attack and sweepback angle. Schleihauf (as cited in [13]) summarised that the ideal pitch angle during the underwater motion of the hand will produce an optimal combination of lift and drag forces, which in turn will generate a resultant force that could be predominantly directed in forward direction. Gardano and Dabnichki [7] noticed that, in addition to the pitch angle having the influence on drag force, at the high Reynolds numbers (which swimmers usually undergo in the competitions), two main factors affect the drag force: shape and arm orientation. In previous CFD studies of free style arm strokes, the effect of swimmer's hand orientation and its effect on propulsion force were evaluated while varying exclusively the angle of attack and full consideration of influence of both the angles were seldom considered.
This quasi-steady state (with initial data is taken from real kinematic experimental research) CFD study is based on flow simulation around swimmer's scanned hand models by using FLUENT (ANSYS, PA, USA). It is hypothesize that the changes of hand shape, velocity, and both angles of hand orientation (pitch and sweepback) complement drag force and drag coefficient in underwater freestyle stroke. Also, it is hypothesised that the realistic variation of both the orientation angles contributes to higher values of lift, drag, and resultant coefficients and forces when compared with values from previous work based on CFD method. Therefore, the aim of this study was to evaluate drag force and drag coefficient corresponding to the variation in kinematic parameters (orientation, shape, and velocity) of swimmer's hand during complete underwater front crawl stroke cycle.

Methods
In this study, the hydrodynamic components of propulsion force ( , ) using realistic models of human hand were calculated. The FLUENT (ANSYS, PA, USA) software was used to simulate the fluid flow, allowing the analysis of distribution of pressure and flow around the hand model of the swimmer.

Hand's Geometries and Initial
Parameters. Artec L 3D scanner (Artec Group Inc., Luxembourg, Lux) was used to scan the left hand of international level swimmer. Working parameters of scanner are video frame rate of up to 15 fps and data acquisition speed up to 500000 points/s. The scanner was directly connected with computer, and all the data was transferred to image-processing program Artec Studio (Artec Group Inc., Luxembourg, Lux). Finally, using this software, cloud point nodal 3D hand models were created from scanned frames: adducted -with thumb adducted; abducted -with thumb fully abducted; spread -with spread fingers; adducted, spread -with adducted thumb and spread fingers. To optimize the models and to get good resolution by the selection of relevant parts of model and to transform the point cloud into a complete 3D polygon model, all the data from Artec Studio (Artec Group Inc, Luxembourg, Lux) was transferred to Leios 2 (EGS srl, Bologna, Italy) software. Typical models of hands with 1182, 1204, 1218, and 1142 surfaces corresponding to adducted , abducted , spread and adducted, spread were created ( Figure 2). The Leios2 (EGS srl, Bologna, Italy) software program can export the datasets into IGES (Initial Graphics Exchange Specification) format which allows file import by CFD software FLUENT (ANSYS, PA, USA).
Initial kinematic parameters corresponding to real swimmer's single front crawl underwater stroke cycle were obtained from previous experimental study [13]. The single front crawl stroke cycle is usually divided into four phases: glide-from the entry of the hand into the water to its maximal forward displacement in the longitudinal displacement; pull-from the maximal forward displacement of the hand in the longitudinal displacement (i.e. direction) corresponding to the stroke time when the hand is located exactly under the shoulder; push-from the end of the pull phase until the exit of the hand out of the water; recovery-from the hand's exit to its reentry into the water [14]. Each model of hand was positioned in nine different positions with the combinations of various pitch (defines angle between hand velocity vector and the plane of the hand) and sweepback angles (which define the inclination of the leading edge of hand) using SolidWorks (Dassault Systémes SolidWorks Corporation, MA, USA) software. Different values of average velocity correspond to different times (Table 1), resulting in three parts of underwater stroke time from glide phase, three from pull and three from push phases. Gourgoulis et al. [13] determined the velocity of the hand as the mean of the resultant velocities of the 2nd and the 5th metacarpophalangeal joints and transformed from the external (O; X, Y, Z) to the local reference system (O; x, y, z) of the swimmer's hand. In CFD calculations, this velocity was used as the resultant  average velocity of water flow while the hand model was kept stationary. To ensure that the model would provide accurate results, the grid was refined in presumed regions of high velocity and pressure gradients. A series of tests were carried out to estimate the independence of the results in relation to the grid resolution. adducted model was used at hand speed of 1.79 m⋅s −1 . The pressure force was calculated for each grid (Figure 3(b)). The number of grid elements equal to 844740 was chosen in further calculations to ensure reasonably fast and accurate computations, the independence of grid used, and repeatability of the results.

Computational Fluid Dynamics (CFD)
Governing Equations. CFD method was used to simulate drag force and resultant drag coefficient of the hand models. Steady velocity is implemented by keeping hand static with the fluid flowing at constant velocity. CFD method is based on the Navier-Stokes equations which fully govern fluid flow. These equations arise from applying Newton's second law to fluid motion, together with the assumption that the fluid stress is the sum of a diffusing viscous term (proportional to the gradient of velocity), plus a pressure term. The solution of the Navier-Stokes equations is a velocity field or flow field, which is a description of the velocity of the fluid at a given point in space and time. For the incompressible fluids, the continuity equation is only function of velocity and not a function of pressure. Only the momentum equations contain pressure term. A direct method is to discretize the equations of continuity and momentum and solve them simultaneously to obtain results of pressure.
CFD methodology consists of a mathematical model that replaces the Navier-Stokes equations with discretized algebraic expressions that can be solved by algorithms on  (1) Navier-Stokes (momentum) equations where ( ) ≡ + is the component of instantaneous velocity in -direction (m⋅s −1 ), is the component of time averaged mean velocity in -direction (m⋅s −1 ), is the component of fluctuating velocity in -direction (m⋅s −1 ), , are the coordinate direction vectors, is average fluid density ) is the turbulent kinetic energy per unit mass (m 2 ⋅s −2 ), and is the Kronecker delta with the condition that, = 1 if = and = 0 if ̸ = .

Boundary Conditions.
For the steady state fluid flow simulations, appropriate boundary conditions were considered. On the left side vertical surface of the domain (inlet velocity, Figure 3(a)), the horizontal component of the initial velocity was applied for all hand positions, respectively, (Table 1) and the vertical component of the velocity was assumed to be equal to zero. The pressure was set equal to zero Pascal on the right side of vertical surface (outlet pressure, Figure 3(a)).
The remaining side surfaces and bottom of the domain were considered as symmetry. Incompressible flow was assumed with turbulence intensity of 1.0% and turbulence scale of 0.10 m. The water temperature was 28 ∘ C with a density of 998.2 kg⋅m −3 and viscosity of 0.001 kg −1 ⋅(m⋅s) −1 with consideration of the gravity of 9.81 m⋅s −2 .

Numerical Scheme.
More accurate solution was considered with the choice of second-order numerical computational schemes. The simulations are based on finite volume method of discretization. In generic terms, the convergence of the calculation is checked by the value of the residuals of the various flow parameters. The convergence criteria in FLUENT (ANSYS, PA, USA) were set at 10 −6 . This criterion is assumed sufficient to ensure the convergence of the solution for the present study. The appropriate number of tetrahedral grids cells in the simulation model was arrived, which was an outcome of grid independence test carried out at the beginning of actual simulations. It was found that the difference in solutions for the drag coefficients for subsequent refinement in tetrahedral grid was less than 1%. In order to limit numerical dissipation, particularly when the geometry is complex consisting of an unstructured grid, as seen in Figure 3(a), the choice of second-order upwind discretization scheme for the convection terms in the solution equations and Pressure-Implicit with Splitting of Operators (PISO) pressure-velocity coupling scheme for the double precision pressure-based solver was chosen. The PISO pressurevelocity coupling scheme, part of the SIMPLE family of algorithms, is based on the higher degree of the approximate relation between the corrections for pressure and velocity.

Results
The calculated values of drag forces, drag coefficients, and areas of hands of all the studied phases at three different velocities (differences from initial velocity ± 0.5 m⋅s −1 ) are presented in Tables 2 and 3, along with the corresponding data of all hand models used. The largest mean of drag force  (40.85 N, = 18.02 N) by adducted was observed during pull phase, which also presented the largest projected hand area, whereas hand models abducted , spread , adducted, spread (appropriately: 43.95 N, = 4.56 N; 51.08 N, = 2.04 N; 41.51 N, = 2.80 N) presented largest mean of drag force during the push phase, when the initial mean velocity was greater than before. The lowest drag forces of all hand models were during the glide phase when initial flow velocity was decreased 0.5 m⋅s −1 and projected hand area was the smallest.
The maximum values of were calculated during the pull phase of all hand models (mean values: 1.96, = 0.36, 1.84, = 0.28, 1.65, = 0.07, and 1.75, = 0.19, accordingly adducted , abducted , spread , and adducted, spread ), when flow velocity was the lowest. Minimum mean values of drag coefficients were reached during the glide phase, when initial flow velocity was increased 0.5 m⋅s −1 and projected hand area was the smallest.
A variation of drag force and drag coefficient is observed due to different hand orientation, which depends on both hand angles.
slightly varied during the glide phase ( Figure 4). All drag forces rose in the beginning of pull phase and reached maximum values in the middle of this phase when projected hand area was the biggest (except spread max during the first point of push phase). The values of the drag forces fell after peaks due to decrease in projected hand area. varied more when initial velocity was increased, and there was less variation when initial velocity was decreased. Drag coefficient slightly varied in glide phase ( Figure 5). The values of rapidly increased from the 3rd position of the glide phase to the 2nd position of the pull phase due to significant changes in hand angles which in turn affected projected hand area. Drag coefficients decreased in push phase of all hand models.
There is a marked improvement in drag force and drag coefficient changes due to different shape of hand (Figures 4 and 5).
was not significant during the glide phase; however, the biggest means were observed in spread and adducted, spread . adducted and abducted reached maximum peaks in the pull phase though the values of drag force were not the biggest during push phase. All drag forces of adducted , abducted , and adducted, spread slightly decreased in push phase initially and increased initial velocities and were almost stable during decreased initial velocity. Meanwhile, of spread was negligibly increasing in push phase during decreased velocity and slightly varied during the initial and increased initial velocities. Drag coefficient varied insignificantly during the glide phase. Significant variation was observed in push phase of all hand models. adducted reached the biggest values of in the pull phase; however, values of were the biggest of abducted in push phase for every variation in velocity.    Pressure forces contribution to hand varied according to water flow velocity, hand orientation, and shape. Pressure visualisation considering different water flow velocity (1.82 ± 0.5 m⋅s −1 ) was proposed on adducted during 2nd position of pull phase ( Figure 6). Pressure contour fields depending on different hand orientation of spread during separated underwater freestyle phases were shown (Figure 7). The shape of hand caused different pressure force for the same water flow velocity (1.82 m⋅s −1 ) during 2nd position of pull phase (Figure 8).
The vectors of flow velocity distribution on the palm of the hand are presented (Figure 9). The velocity of flow decreased when in contact with hand surface. Flow dispersed to the corners of the palm with increasing in velocity. The

Discussion
The purpose of this study was to evaluate drag force and drag coefficient correspondence to variation in kinematic parameters (orientation, shape, and velocity) of swimmer's hand during complete underwater front crawl stroke cycle. Four scanned hand models with fingers in different positions were generated. The orientation of hand was set into nine different positions with accordance to average velocity. In the present study, the analysis of a particular case of front crawl underwater movement (Table 1) and initial kinematic data were taken from experimental kinematic research work [13]. Calculations were performed through FLUENT (ANSYS, PA, USA) software, which is based on computational fluid dynamics method applied to accurately solve fluid flow problems through numerical simulation. This method has been proven to generate accurate results with repeatability of results with identical values, when performed with similar initial conditions and settings [17]. This saves time, and also the results can be accessed in detail and analysed anytime, unlike repetitive experimental tests. The outcome of this study is that the modification of three variables: hand orientation, shape, and the average water flow velocity strongly affects drag coefficient and drag force of hand and with this in some measures affects the hand propulsion.
There is clear indication of progress in the understanding of the contribution of hand shape, orientation, and velocity on drag coefficient ( ) and drag force ( ) in underwater hand stroke. Obvious drag force and drag coefficient dependence on velocity are shown in Figures 4 and 5. The values of are seen to increase and decrease with respect to the initial velocity, which was accordingly increased and decreased in 0.5 m⋅s −1 intervals. A similar tendency of increase of with velocity was observed in a previous study performed under steady-state flow conditions of the 2D and 3D swimmer's hand model [11,18]. Rouboa   orientations where water flow velocity was increased from 0.5 m⋅s −1 to 4 m⋅s −1 . Sato and Hino [18] used hand model with spread fingers in a fixed direction with water flow, and velocity was set from 0.5 m⋅s −1 to 2 m⋅s −1 . However, the variation of current drag coefficient ( Figure 5) and previous findings disagree. It was indicated that the drag coefficient did not vary, while there was rise in water flow velocity [9,11,18]. Berger et al. [6] towed model of human hand in a towing tank at the similar orientation and evaluated hydrodynamic parameters. It was shown that the slightly decreased within the velocity range from 0.7 to 3.0 m⋅s −1 , and at the velocity lower than 0.7 m⋅s −1 strongly depends on velocity. Current findings indicated that the rise in velocity generated lower values of drag coefficients for all hand models during all phases. According to current findings, it can be concluded that the velocity affects values of drag force and drag coefficient during all front crawl underwater phases.
To our knowledge, there is no known work based on CFD method which compares different phases of underwater front crawl hand motion with appropriate changes of hand orientation. However, Gourgoulis et al. [19] calculated drag force of human hand by formulas while input data was taken from kinematic experiments. The drag force dependence on hand orientation appears with the variation of in both current ( Figure 4) and previous research. There is clear indication that, between different phases of underwater hand movement, when swimmer changes his hand orientation (i.e., pitch and sweepback angle) (Table 1), the changes directly affect the projected frontal area and hereupon the drag force. The variation of calculated from experimental data was similar with current computational hand model abducted with similar initial flow velocity. Therefore, we could imply that the swimmer maintained his hand with thumb abducted during underwater hand stroke experimental tests.
The calculation of drag coefficient during different underwater phases of hand model allowed evaluating the influence of hand orientation on . This study showed that small changes in pitch and sweepback angles during glide phase caused negligible variation of for all hand models ( Figure 5). The sharp variation of hand angles during the pull phase caused major changes in values of , and the dwindling of angles to smaller values shaped the reduction of drag coefficient during push phase. The consideration of the angle of attack and sweepback angle influenced contributions to drag coefficient during complete underwater cycle by different hand models.
The importance of hand shape variation and its contributory influence to drag force are presented (Figure 4). Different individual curve lines in the pictorial graph represent  Figure 10: Comparison of drag coefficient versus angle of attack acting on swimmer's hand: (a) current CFD results (with appropriate sweepback angle and flow velocity in each phase), Bixler and Riewald [5] CFD (velocity 2 m⋅s −1 , sweepback angle 0 ∘ ), Berger et al. [6] experimental results (velocity 1 m⋅s −1 , sweepback angle 0 ∘ ) and Schleihauf [15] experimental results (velocity 2 m⋅s −1 , sweepback angle 0 ∘ ); (b) current CFD study (with appropriate sweepback angle and flow velocity in each phase), Sato and Hino [12] CFD, and Takagi [16] (cited in [12]) experimental results where flow velocity is 1 m⋅s −1 and sweepback angle is 0 ∘ . separate models of the swimmer's hand. The highest peak value of drag force was obtained from the hand model with thumb adducted and the smallest peak from the model with all fingers spread in the middle of the pull phase regardless of variation in flow velocity. Although the hand area, that directly affects the drag force, was larger in the model with fingers spread than in that with the thumb adducted, the pitch angle was the largest in this position. However, in the end of pull phase and throughout the push phase, when the hand altered its leading edge, the hand with spread fingers generated the biggest projected area of total plane area of hand and values of , whereas the hand with the thumb adducted-the least. It means that the shape of the hand had different contributions on the drag force during separate underwater phases of front crawl stroke cycle. Likewise, it was summarised that higher drag force values could be reasonable by a hand model with optimal spacing between fingers [10,12]. According to the current results, it is clear that, in order to increase the drag force contribution throughout the underwater stroke cycle, the shape of the hand should be altered in separate individual phases. The drag coefficient values ( Figure 5) corresponding to different hand models were evaluated in recent CFD studies under steady state with spread fingers [10,12,18], with thumb in different positions [9], while current study used different values of velocity, pitch angle, and also different shapes of hand models. In previous studies, the drag coefficient was almost constant despite the changes in velocity. We observed that the drag coefficient was not constant and it varied throughout underwater front crawl cycle regardless of differences in hand model.
In the current study, we compared drag coefficients from the hand model with thumb adducted with previous studies obtained experimentally and from CFD studies, which were achieved by varying only the angle of attack in appropriate range while sweepback angle was equal to zero [5,6]. The comparative study of current and previous values versus angle of attack (Figure 10(a)) showed the similar tendency of variation in (the hydrodynamic characteristics between all research work were similar). However, the means of current results were higher when compared with previous ones. These perceived differences could be due to differences in the usage of distinct input conditions, and, according to Berger et al. [6], material, hand size, and shape of the model might also account for the observed differences in . Moreover, the distinct tendency of variation can also be caused by different hand orientation, which considered variation of both pitch and sweepback angles at the same time corresponding to that observed under real swimming conditions during underwater hand stroke path.
The drag coefficients of hand model with fingers spread obtained in this study were compared (Figure 10(b)) with those in Takagi's (as cited in [18]) study, obtained from experimental research, and with those in Sato and Hino's work [18] calculated through CFD method. The kinematic data from experimental work [13] and scanned human hand model were used for this study to evaluate drag force and drag coefficient; therefore, values and its variation present practical and further pragmatic results. Obvious differences are observed in values of the drag coefficients between present and previous CFD and experimental studies; whereas the tendency of variation was similar. Higher values in current work could be influenced by the choice of fluid flow model applied in the study: turbulent fluid flow was considered in this study, and laminar fluid flow was considered by Sato and Hino [18]. The values of were significantly higher in this study as the was dependent on not only the pitch angle, as considered in previous studies, but on both angles of hand orientation (pitch and sweepback). Furthermore, the consideration of experimental kinematic data of swimmer's hand motion (initial velocity, pitch, and sweepback angles) and realistic fluid flow model can affect realistic estimation of drag coefficient.
It was shown that the pressure force contribution on swimmer's hand acted differently vis-à-vis water flow velocity, hand orientation, and shape. The results of this study confirmed that the velocity of water flow directly affected pressure force on swimmer's hand ( Figure 6). The contribution of pressure force to hand increased when velocity was increased. There is strong evidence that the pressure force affected different parts of hand due to altered hand sweepback angle and bigger part of hand frontal cross area came up with maximum pressure force, when pitch angle was increased (Figure 7). This maximal pressure force shifted towards hand's leading edge; this is when the swimmer altered orientation of the hand, that is, the angle of attack and sweepback angle. Different hand shapes arrived at different pressure force distributions (Figure 8). In this case, the biggest maximal value of pressure (6250 Pa) acted on hand with spread fingers (Figure 8(c)) and the least (5720 Pa)-on hand with thumb adducted (Figure 8(a)). The present study corroborates with Marinho et al. [10] and Minetti et al. [12], with reference to the maximum values of the pressure force, suggesting that the hand with spread fingers goes through higher pressure force and in turn influences higher drag force than other forms of hand model.

Conclusions
This study, based on the CFD method, indicated that drag force and drag coefficient were significantly affected by velocity, hand frontal cross-section area, and both angles of hand orientation (pitch and sweepback) during complete underwater front crawl single cycle. Higher values of hydrodynamic parameters in this study are attributed to the consideration of real swimming conditions during simulations with the consideration of combination of both hand orientation angles observed during individual phases. To increase drag force contribution throughout the underwater stroke, the shape of the hand should be altered in separate phases.