Role of Shape and Kinematics in the Hydrodynamics of a Fish-like Oscillating Hydrofoil

: In the present two-dimensional numerical study, we investigate the roles of geometrical parameters of a hydrofoil (shape/curvature of the leading and trailing edges and thickness) and kinematic parameters (phase difference between heave and pitch ( φ )) on the propulsive performance of different-shaped hydrofoils oscillating at maximum angles of attack up to α max = 30 ◦ . The study was carried out at a ﬁxed non-dimensional maximum heave to chord ratio h ◦ / C = 0.75, Strouhal number St = 0.25, and Reynolds number Re = 5000. Our ﬁndings reveal that hydrofoil performance and stability improve with leading and trailing edge curvatures but decline as thickness increases. By analyzing the near-wake structure, we establish that even minimal ﬂow separation increases power consumption while moderate ﬂow separation enhances thrust. Over the range of different-shaped hydrofoils at different α max and φ , maximum propulsion efﬁciency occurs for those parameters for which there is a small degree of ﬂow separation but with no roll-up of a separating vortex. In comparison, maximum thrust generation occurs when there is a moderately strong ﬂow separation but without induction of a signiﬁcant amount of ﬂuid around the trailing edge. These insights offer valuable knowledge for understanding ﬁsh propulsion efﬁciency and have applications in designing autonomous underwater vehicles (AUVs) and micro-air vehicles (MAVs).


Introduction
The problems pertaining to fluid-structure interaction (FSI) can be classified into two primary categories: one-way or weakly coupled, and two-way or strongly coupled problems.In one-way or weakly coupled problems, only the fluid flow field is influenced by the structure/body motion.On the other hand, in two-way or strongly coupled problems, both the structure/body and the fluid flow field exert mutual influence upon each other [1].The present work is on a one-way coupled FSI problem of an oscillating hydrofoil.
Especially over the past decade, the hydrodynamics of an oscillating hydrofoil has drawn the attention of many engineers, mathematicians, and biologists due to its association with understanding the biology of fish swimming and designing new innovative solutions for modern engineering applications [2].Some of these applications include micro-air and water propulsion systems [3][4][5], and energy-harvesting systems [6][7][8][9].The motivation of the present study is to understand the role of hydrofoil-shape on its propulsive performance.This knowledge is fundamental and holds significant potential for enhancing the precision and appropriateness of engineering designs for autonomous underwater vehicles (AUVs) and micro-air vehicles (MAVs).
For the present problem of a thrust-generating oscillating hydrofoil, the non-dimensional governing parameters can be categorized as two geometric parameters: the curvature (of leading and trailing edges) and thickness of the hydrofoil; four kinematic parameters: heave to chord ratio h • /C, Strouhal number St = 2h • f /u ∞ ( f is the frequency of oscillation), maximum angle of attack α max (Equation ( 2)), and the phase difference between heave and pitch φ (Equation ( 1)); and Reynolds number Re = ρu ∞ C/µ as the flow parameter based on free stream velocity u ∞ , chord length C, and dynamic viscosity µ.Studies on oscillating hydrofoils in the literature use a range of different NACA hydrofoils.To provide context, we first sequentially review previous experimental and numerical studies in this area.
In terms of studies focusing on experiments, Triantafyllou et al. [10,11] undertook a stability analysis of the wake of an oscillating hydrofoil and found that its stability is governed by a non-dimensional parameter, the Strouhal number, St, defined above.The wake was found to be stable for the frequency range of 0.25 < St < 0.4.Further, they determined the operational St for different fish and found that most propel themselves within this narrow St range.Later, through a parameter study, Anderson et al. [12] determined the optimal parameters for efficient oscillation of a NACA0012 hydrofoil, which replicates the caudal fin motion for carangiform fish.Their results showed that the optimal efficiency with high thrust also lies within this narrow St range.In addition to oscillation frequency, they also examined optimal values for three more non-dimensional governing parameters-heave amplitude h • /C, maximum angle of attack α max , and the phase difference between heave and pitch φ.For efficient propulsion, they proposed that the values of these parameters should be 0.25 < St < 0.4, h • /C ≈ 1, 15 • < α max < 25 • , and φ = 75 • .Following this, Read et al. [13] also undertook a parameter study for the same set of governing parameters as Anderson et al. [12].However, for almost all cases, they observed that the oscillating hydrofoil achieved high efficiency together with high thrust for 90 • < φ < 100 • .These two contradictory findings led to numerous further experimental and numerical studies with the prime objective of exploring the effect of φ on propulsive performance.Motivated by these results, Hover et al. [14] compared the effect of different angle of attack profiles on the hydrodynamic performance of a flapping hydrofoil.They observed that a cosine angle of attack profile provides an optimal combination of thrust and propulsive efficiency, as compared to a harmonic angle of attack profile.Recently, Van Buren et al. [15] also performed an experimental investigation on a pitching and heaving teardrop-shaped hydrofoil.They observed that a phase difference of 90 • between heave and pitch is required to achieve maximum propulsive efficiency.Additionally, they proposed scaling laws based on lift and added mass forces that align well with the experimental measurements.
Many numerical studies have also been undertaken on the problem of flow across oscillating hydrofoils.Isogai et al. [16] were among the first who performed two-dimensional simulations to understand the effect of φ on the propulsive performance of a NACA0012 oscillating hydrofoil.They concluded that the highest propulsive efficiency occurs at φ = 90 • .A similar conclusion was reached in other studies [17,18].Numerical studies have also been conducted to understand the effect of near-and far-wake vortices on the propulsive performance.For near-wake vortices, Tuncer and Platzer [19] concluded that the formation of even a small leading-edge vortex (LEV) should be prevented for efficient propulsion.Similar conclusions have resulted from other studies, including Tuncer and Kaya [20] and Young and Lai [21].In direct contrast, Guglielmini and Blondeaux [22] concluded that a small to moderately strong magnitude LEV is beneficial for propulsive efficiency.Thus, the role of LEV in efficient propulsion is contradictory in the literature.For the far-wake structures, Bose and Sarkar [23] investigated the transition in the flow pattern from periodic to chaotic through a quasi-periodic route by increasing the heave amplitude at Re = 1000.Numerical studies have also been conducted to understand the effect of various geometrical parameters on propulsive performance.Zhang et al. [24] studied the effect of chord and thickness ratio of an elliptical foil, while Ashraf et al. [25] studied the effect of thickness and camber using a variety of NACA hydrofoils.
The above brief literature survey shows that most previous studies focused on understanding the role of kinematic parameters on the performance of an oscillating hydrofoil, and few published works discuss the role of the hydrofoil shape on performance parameters.However, the effect of the hydrofoil curvature at the leading and trailing edges, along with its thickness, remains to be explored.Thus, one of the primary aims of the present study is to elucidate the effect of these geometrical parameters on the propulsive performance of an oscillating hydrofoil at various angles of attack, α max .Through this analysis, we are also able to comment on the role of LEV on the performance of a thrust-producing oscillating hydrofoil.Next, the phase difference between the heave and pitch φ is a crucial parameter that governs the propulsive performance.Hence, its role in the propulsive performance of a NACA00XX hydrofoil has been studied extensively, as discussed in the above paragraphs.Those studies reach similar conclusions that a phase difference of 90 • provides the maximum propulsive efficiency with an optimal magnitude of thrust.A change in the phase difference may enhance the thrust generation; however, it significantly reduces propulsion efficiency.Note that the motion of an oscillating hydrofoil at φ = 90 • seems to be the most fish-like, cleanly slicing through the water with the lowest angles of attack.Further, fish are also known for maintaining a phase difference of φ = 90 • between pitch and heave during swimming [26].Thus, the next objective is to understand whether φ = 90 • is a universal phase difference for attaining maximum efficiency or whether it depends on hydrofoil shape.Being the primary mechanism used by various fish for swimming and other engineering applications discussed above, the present study contributes from both fundamental and engineering points of view.
The paper structure is as follows.The methodology is detailed in Section 2, which provides a description of different hydrofoils considered for the study, oscillation kinematics, governing equations, validation and grid independence studies, and definitions of the input and output governing flow parameters.This is followed by a discussion of the results in Section 3. The main conclusions are summarized in Section 4.

Mapping the Hydrofoil Shape
For the study on the role of curvature, specifically the front and rear shape, and thickness of the hydrofoil, Figure 1 depicts various shapes that are considered in the present study.The different shapes are chosen in such a way that they connect the two extreme shapes (a NACA0012 airfoil and a thin flat plate) by considering the importance of the front and rear geometry, and thickness.Here, shape 1 corresponds to a NACA0012 hydrofoil with a bluff trailing edge; shape 2 corresponds to a NACA0012 hydrofoil with a bluff leading edge; the rectangular plate (effectively) corresponds to a NACA0012 hydrofoil with bluff leading as well as trailing edges; and the thin plate (effectively) corresponds to a NACA0012 hydrofoil or a rectangular plate with negligible thickness.Thus, comparing the NACA0012 hydrofoil with shape 1 explores the role of a tapering tail, a NACA0012 hydrofoil with shape 2 helps us understand the role of a rounded leading edge, a NACA0012 hydrofoil transformed into a rectangular plate helps to understand the combined role of a rounded edge and tapered tail, and finally a NACA0012 airfoil transformed into a thin plate investigates the role of thickness on the propulsive performance.

Oscillation Kinematics and Other Parameters of Interest
The oscillatory motion of the caudal (i.e., rear) fin of a fish can be modeled as a coupled pitching and heaving symmetric hydrofoil with a high aspect ratio [11,27].This high aspect ratio allows the derivation of the three-dimensional flow characteristics around the foil from the two-dimensional flow considerations [12].Figure 2 shows the schematic figure of a carangiform/thunniform fish about the medio-frontal plane.The hatched region shows the modeled caudal fin performing coupled pitching-and-heaving motions about the pivot point P in a freestream flow of velocity u ∞ .The corresponding motion can be replicated by sinusoidal motion for both the heave and pitch [13], given as where h 0 is the maximum heave amplitude, f is the frequency, θ 0 is the maximum pitching amplitude, and φ is the phase-angle difference between heave and pitch.Notice from the figure that due to heaving motions, the hydrofoil moves at an additional velocity, the heave velocity u heave = dh/dt, which results in a change in the resultant velocity (u res ) of the foil.Thus, the instantaneous angle of attack α(t) is no longer equal to pitching angle θ(t) [12] and can be given as: For analyzing the output performance, the non-dimensional parameters considered in the study are the thrust coefficient C T , lateral force coefficient C L , input power coefficient C P , and propulsive efficiency η P , defined as, respectively, Here, F T and F L are the net forces acting in the streamwise and lateral direction, respectively, and are calculated from the integration of pressure and viscous forces over the hydrofoil.Further, ρ is the density of the surrounding fluid, c L is the local lateral force coefficient per unit surface area of the hydrofoil, and V body is the lateral velocity of an oscillating hydrofoil.
The propulsive efficiency η P is defined as the ratio of the net output power to the input power [28,29], where an overbar represents time-averaged values obtained by using at least twenty oscillation cycles following statistical convergence of the simulations.

Governing Equations and Numerical Details
For an oscillating hydrofoil, Zurman et al. [30] found that the flow and forces are effectively two-dimensional for a NACA0016 hydrofoil at St = 0.15-0.45and Re = 5300.Thus, for the present cases at St = 0.25 and Re = 5000, 2D numerical simulations are performed using an in-house code based on a level-set function-based immersed interface method (LS-IIM).The LS-IIM-based code was developed in the CFD lab at IIT Bombay by Thekkethil [31], in which the temporal variation in the solid-fluid interface is tracked by using a normal signed distance-based level set function, Φ.At each time step, the value of Φ is calculated using the minimum distance and winding algorithms, governed by the predefined kinematics of oscillation (discussed in Section 2.2).The positive Φ values represent fluid cells, whereas negative values represent solid cells.The fluid dynamics of the present problem are governed by the incompressible Navier-Stokes equation.In non-dimensional form, these are given by Momentum : where − → U (≡ − → u /u ∞ , u ∞ is the freestream velocity) is the non-dimensional velocity vector and P (≡ p/(ρu 2 ∞ ), ρ is the density of the surrounding fluid) is the non-dimensional pressure.The LS-IIM uses the finite-volume method to discretize the above equations in space.Further, it is based on a fully implicit pressure projection method on a collocated grid.The non-linear advection term of the Navier-Stokes equation is integrated using the QUICK (Quadratic Upstream Interpolation for Convective Kinematics) scheme, while the CDS (Central Difference Scheme) is used for the diffusion term.More details of LS-IIM can be found in Thekkethil and Sharma [1].

Computational Details
Figure 3 shows the non-dimensionalized computational domain and the boundary conditions for the present problem on flow past an oscillating hydrofoil.The chord length C and the freestream velocity u ∞ are taken as the length and the velocity scale for the problem.The x-direction is aligned with the streamwise direction, while the y-direction is aligned with the lateral direction of the flow.The domain measures 16C in the streamwise direction, with the leading edge of the hydrofoil placed 5C from the inlet.The lateral dimension of the domain is 10C, with the hydrofoil located at the center.This provides a maximum lateral blockage (0.12C/10C) of 1.2%, ensuring a minimal effect of the lateral boundaries on the hydrodynamics of oscillation.In Figure 3, U and V correspond to the non-dimensional velocities in the streamwise and lateral direction, respectively.The boundary condition at the inlet corresponds to a constant velocity (equal to the freestream velocity, u ∞ ).A convective boundary condition, with a convective velocity U c equal to one, has been used at the outlet.A free-slip or symmetry boundary condition is used for the lateral boundaries.The hydrofoil is assumed to be non-deforming, and a no-slip (for the velocity) and a Neumann boundary condition (for pressure) have been employed at its surface.
For the grid independence study, we considered three non-homogeneous structured grids of sizes: 382 × 288, 768 × 552, and 1192 × 954.For each grid, a uniform size fine cell of size ∆x = ∆y = 0.01, 0.05, and 0.0025 is used in the region closer to the hydrofoil, and a uniform coarse cell size of ∆x = ∆y = 0.5, 0.25 and 0.125 is used in the far field.The finest and coarsest grid zones are connected using a hyperbolic stretching function, with a maximum expansion ratio of less than 1.01.For α max = 30 • , h • /C = 0.75, St = 0.25, and Re = 5000, Figure 4a shows no significant visual difference in the temporal variation in the lateral force coefficient C L between the two finest grids.Therefore, all further simulations for the present problem are carried out using the intermediate 768 × 552 grid size, noting that it has 200 points across the hydrofoil chord.The simulations are carried out in series on the local system with 32 GB RAM for a total of 48 simulation hours per case.
The numerical method used for the present study has already been applied successfully to various similar fluid-structure interaction (FSI) problems [1,[32][33][34][35][36][37].However, a further validation study has been carried out by comparing our results with that documented by Boiron et al. [38] and is presented in Figure 4b.The figure shows an excellent agreement between the present and published results for the temporal variation in the thrust and lift coefficients at α max = 20 • , h • /C = 0.75, St = 0.3 and Re = 4000.The percentage of difference is within 2% and 5% for the lift and thrust coefficients, respectively.

Results
The results section is divided into three subsections.Section 3.1 discusses the effect of varying the maximum angle of attack, α max , on the propulsive performance of different shaped hydrofoils.Section 3.2 discusses the correlation between the time-averaged engineering parameters and instantaneous flow parameters using a cause-and-effect-based study and instantaneous vorticity in the vicinity of the hydrofoils at the extreme rightward lateral position.Section 3.3 examines hydrofoil cases at maximum efficiency to study the role of the phase difference φ between heave and pitch.In both Sections 3.1 and 3.3, the time-averaged performance parameters are calculated over twenty oscillation cycles following the statistical convergence of the simulations.

Effect of Curvature Variations on Hydrodynamic Performance
Considering the five different shapes of the hydrofoil (Figure 1), Figure 5 presents the effect of curvature and thickness on the time-averaged propulsive performance of a hydrofoil for various angles of attack, α max = 5 • -30 • , at St = 0.25, h 0 /C = 0.75, φ = 90 • and Re = 5000.This leads to the pitching amplitude θ 0 ∼ 34 • -8 • (Equation ( 1)).For the mean thrust coefficient C T , Figure 5a shows that the magnitudes of C T for the NACA0012 and those with a bluff leading/trailing edge (shape 1 and shape 2) are bounded within the magnitude of C T for the thin plate and thick plate.With an increase in α max , the figure shows an increasing-decreasing trend for C T with its maximum close to α max = 22.5 • for the different shapes; the exception is shape 2, which achieves its maximum C T at α max = 17.5 • .Further, at lower values of α max , C T for bluff frontal edge-based shape 2 approaches that of a thin plate, while C T for bluff trailing edge-based shape 1 approaches that of a thick plate.However, at higher values of α max , the reverse applies; that is, C T for shape 1 approaches C T of a thin plate, while C T for shape 2 approaches that of a thick plate.The figure also shows the drag-to-thrust transition for different shapes with increasing α max , except for the thin plate with a net thrust even at a very low α max = 5 • .For other hydrofoils, the figure shows that the transition occurs first for hydrofoils with a tapered trailing edge (NACA0012 and shape 2), compared with hydrofoils with a bluff trailing edge (shape 1 and the thick plate).
For the mean input power coefficient C P , Figure 5b shows an apparent increasing C P with increasing α max , with C P for the different shapes bounded within the maximum C P for shape 2 and minimum C P for the NACA0012 airfoil.Further, the figure shows that C P is smaller for shapes with a rounded frontal edge (NACA0012 airfoil, shape 1) compared with bluff frontal edges (the flat plate, shape 2 and the thick plate), and the trend is independent of α max for the hydrofoil.Moreover, for the bluff frontal-edge hydrofoils, C P increases with frontal surface area, i.e., C P is larger for shape 2 and the thick plate compared to that for the thin plate at any α max .
Similar to the variation with C T , Figure 5c shows that η P for the different shapes is again bounded between the values for the thin and thick plate for any α max .Note that the efficiency curve contains only those cases that resulted in a net mean thrust.Furthermore, the figure shows that the nature of the η P curve is also of increasing-decreasing type.However, for three of the foil shapes, the maximum η P occurs at α max = 15 • -17.5 • ; this is not the case for the thin plate and shape 2, which have maximum η P at α max = 12.5 • .At larger angles of attack, α max ≥ 20.0 • , the figure shows that η P for the NACA0012 matches that of a thin plate, whereas the η P for shape 2 matches that of a thick plate.Further, by comparing shape 1 with the NACA0012 airfoil and shape 2 with the thick plate, it can be seen that η P is greater for shapes that possess a tapered end.However, the benefits of possessing a tapered end decrease with increasing α max .
For the rms variation in the lift coefficient, C L,rms , Figure 5d shows that the variation is similar to the variation in C P with α max , as discussed above.For instance, the magnitude of C L,rms increases with α max ; for any α max , the magnitude of C L,rms for different shapes is bounded within the maximum C L,rms for shape 2 and the minimum C L,rms for the NACA0012 airfoil, showing that shape 2 is the least stable and the NACA0012 airfoil is the most stable hydrofoil for propulsion.Further, it is important to notice that the C L,rms of shapes with a rounded frontal shape is smaller compared to shapes with a bluff frontal surface, which leads to a more stable propulsion for shapes with a rounded frontal area.
The above-discussed increasing-decreasing trend for C T and the increasing trends for C P and C L,rms with increasing α max are at least partially understandable as these parameters are directly proportional to the surrounding fluid displaced by the hydrofoil.Thus, these parameters are directly proportional to the projected lateral surface area (C sin(α)) of the hydrofoil.To discuss this, Figure 6 shows instantaneous schematics of a hydrofoil at three values of α max (α 1 < α 2 < α 3 ).For simplicity, we have assumed that the hydrofoil is oscillating in a still fluid, resulting in the vertical relative velocity u res .With increasing α max , the figure shows that the projected lateral surface area of the hydrofoil increases, resulting in an increased volume of displaced fluid.Thus, a larger C P is required for an oscillating hydrofoil to perform its movement at a larger α max .In contrast, for C T , there is a trade-off between increasing surface area and decreasing axial component, resulting in an increasing-decreasing variation.However, it is important to notice from Figure 5 that the magnitude of engineering parameters is different for different shapes.Further, the α max at which maximum C T and η p occur is different for different shapes.Thus, a more detailed study is needed; this is attempted in the next section (Section 3.2) by correlating the engineering parameters with the flow structures near the vicinity of the hydrofoil.

Flow Characteristics
The correlation between the above-discussed time-averaged engineering parameters with the instantaneous flow parameters is accomplished in two steps.In the first step, the instantaneous pressure contours in the vicinity of the hydrofoil are correlated with the corresponding instantaneous vortex structure and engineering parameters.For ease of understanding, the study is conducted only for an oscillating NACA0012 hydrofoil for a range of α max .The corresponding study is known as a unified cause-and-effect-based CFD analysis, and has been proposed by Sharma [39].The study helps to correlate the time-averaged engineering parameters with the instantaneous vortex structure when the NACA0012 hydrofoil is at the rightmost lateral position.In the second step, the correlation developed in the first step is applied to differently shaped hydrofoils.
Over one oscillation period, Figures 7 and 8 show a unified cause-and-effect-based analysis by plotting the time variation in 2D plots for flow properties and 1D line-plots for propulsive performance parameters for α max in the range 5 • ≤ α max ≤ 30 • .The 2D contour plots Figures 7a-i and 8a-i are for the vortex structure and pressure contours in the vicinity of the hydrofoil, respectively, while the 1D plots, Figure 7j,k, are for the thrust C T and power C P coefficients.
At τ/T = 0, Figure 7 shows the onset of the leftward movement of the hydrofoil from its extreme right lateral position.At this instant of time, a net drag acts on all the hydrofoils independent of α max ; however, the magnitude of drag is larger for hydrofoils with a higher α max , as shown in Figure 7j.The reason behind this increase in drag with α max can be seen in the instantaneous vortex contours of Figure 7(a1-a6).For α max = 5 • , the figure shows a fully attached flow across the hydrofoil, while the flow separates for hydrofoils with α max > 5 • -in this case, the flow starts to separate from the posterior region of the foil, and the point of separation moves towards the leading edge and becomes stronger with increasing α max .The stronger flow separation causes a stronger pressure gradient along the chord (Figure 8(a1-a6)), resulting in higher drag.
Figure 7(b1-b6) show that the hydrofoil rotates in the counter-clockwise (CCW) direction with further leftward movement.With this movement to τ/T = 0.1, Figure 7j shows the crossover from instantaneous drag to thrust for all α max , except for α max = 5 • for which crossover occurs later at τ/T = 0.25.Further, the figure shows that the magnitude of instantaneous C T at τ/T = 0.1 increases with increasing α max , reaches its maximum at α max = 20 • , and then decreases.In contrast, Figure 7k shows that the instantaneous power coefficient C P increases with α max .The corresponding vortex structure in Figure 7(b1-b6) shows that the increase in instantaneous C T and C P with α max is due to an increase in the strength of flow separation.However, above α max = 20 • , the strength of the flow separation is too large such that it draws the flow around the trailing edge of the hydrofoil, resulting in a reduction in net instantaneous thrust.With the further leftward movement to τ/T = 0.25, Figure 7j,k show that the hydrofoils reach their maximum instantaneous C T and C P , except for α max = 5 • .Further, the figure shows that the largest instantaneous C T occurs for α max = 25    With further leftward movement, Figure 7(d1-d5) show that the CCW rotation of the hydrofoil switches to clockwise (CW) rotation, and the hydrofoils decelerate, reaching the extreme leftmost position.This deceleration of the hydrofoil reduces the instantaneous positive leftward pressure for all the cases, except for α max = 5 • -compare Figure 8(c1-c5) with Figure 8(d1-d5).The reduction in pressure reduces the pressure gradient, resulting in a reduction in the instantaneous thrust and power for all cases, except for α max = 5 • , as shown in Figure 7j,k, respectively.Further, the switch in rotation from CCW to CW enhances the flow separation from the hydrofoil's rightward side, resulting in the continuation of the trend of increasing instantaneous C T with α max up to α max = 25 • and then decreasing, as observed for τ/T = 0.25.The retardation continues as the hydrofoil approaches the extreme leftmost position (completing half a cycle), resulting in a further decrease in instantaneous thrust and power, as shown in Figure 7j,k.
For the subsequent half cycle, where the hydrofoil moves from the extreme leftmost to the rightmost position, Figures 7 and 8 show an opposite nature of periodic variation in vortex structures and pressure on the two sides of the foil, resulting in similar instantaneous thrust and power curve.P:-0.20 -0.18 -0.16 -0.14 -0.12 -0.09 -0.07 -0.05 -0.03 -0.01 0.03 0.05 0.07 0.10 0.12 0.14 0.01 0.16 0.18 Therefore, by correlating the instantaneous vortex structure with the thrust and power coefficients for a cycle of oscillation of the NACA0012 hydrofoil at different α max , it can be concluded that any amount of flow separation enhances the power consumed by the hydrofoil for performing its oscillation motion.In contrast, a moderately strong degree of flow separation (until it becomes so strong that it starts pulling a significant amount of fluid around the trailing edge) enhances the thrust.For further investigation, Figure 9 compares the instantaneous vortex structure in the vicinity of different shaped hydrofoils at their extreme rightward lateral position for various α max .It can be seen from the highlighted rectangular regions in the figure that the η max for any shape occurs at that α max , for which there is a small degree of flow separation with almost no roll-up into a separating vortex.On the other hand, the highlighted cases for C T,max are those at which moderately strong flow separation occurs with no induction of a significant amount of fluid transferring around the trailing edge.Figure 9 also shows the effect of α max and geometrical parameters on the flow separation.For any shape, an increasing α max increases the magnitude of flow separation and brings the point of separation close to the leading edge.Geometrical features such as a bluff leading edge, tapered tail, and high thickness enhance the flow separation, while those such as a rounded leading edge and a flattened tail reduce the magnitude of flow separation.

Effect of Phase Difference between Heave and Pitch φ on Hydrodynamic Performance
In this section, the effect of the phase difference between heave and pitch, φ, on the propulsive performance of the hydrofoils is presented for cases at η p,max , as observed in Section 3.1.This corresponds to α max = 15 • for the NACA0012 airfoil and thick plate, α max = 17.5 • for shape 1, and α max = 12.5 • for shape 2 and the thin plate for phase differences φ = 75 • -105 • , all at constant St = 0.25, h o /C = 0.75 and Re = 5000.
Figure 10 shows the variation in the normalized C T and η p with phase difference φ for different shaped hydrofoils at η p,max , i.e., α max = 15 • for the NACA0012 airfoil and thick plate, α max = 17.5 • for shape 1, and α max = 12.5 • for shape 2 and the thin plate.Here, the normalized values are obtained by normalizing with corresponding values at φ = 90 • .For all cases, the figure shows that both the maximum thrust and efficiency occur in a narrow range of phase difference 90 • ≤ φ ≤ 100 • ; this indicates that the optimal phase difference for propulsion is independent of the shape and shows the hydrodynamic (and evolutionary) reason behind its adoption by different body-caudal fin (BCF) fish in nature.The occurrence of maximum efficiency and thrust coefficient at phase angles 90 • ≤ φ ≤ 100 • is explained here using the analogy developed in Section 3.1 between the instantaneous wake structure in the vicinity of hydrofoil (at maximum heave position) and the engineering parameters.Similar to Figure 9, Figure 11 also shows that the maximum efficiency occurs at those φ for which there is a small degree of flow separation with no roll-up of a separating vortex.Similarly, the maximum C T,max occurs at those φ for which moderately strong flow separation occurs with no significant induction of fluid around the trailing edge.
Figure 11 also shows the effect of the phase difference φ and geometrical parameters on the instantaneous vorticity in the vicinity of the hydrofoils.It is seen that the flow separation for all the shapes reduces by decreasing the phase difference below φ = 90 • .However, the flow separation enhances, and the point of separation travels posterior (towards the tail) by increasing the phase difference above φ = 90 • .Further, by increasing the phase difference above φ = 90 • , secondary flow separation occurs from the rightward leading edge for shapes with a bluff leading edge (shape 2, the thick plate and the thin plate).ω: -10.00 -8.24 -6.48 -4.72 -2.96 -1.21 0.55 2.31 4.07 5.83 7.59 9.35

Conclusions
The present study documents and contributes towards understanding the role of an oscillating hydrofoil's curvature/shape at the leading and trailing edges, along with its thickness, in the propulsive performance at various maximum angles of attack in the range α max = 5 • -30 • , and for phase differences between heave and pitch of φ = 75 • -105 • , at constant St = 0.25, h o /C = 0.75, and Re = 5000.The conclusions drawn are as follows.
In terms of the role of leading-and trailing-edge curvatures, a tapered tail provides (compare NACA0012 with shape 1) a more significant thrust, consumes slightly lower power, and generates a lower C L,rms .On the other hand, a rounded leading edge (compare the NACA0012 airfoil with shape 2) provides a smaller and larger thrust for α max < 15 • and α max > 15 • , respectively.However, the power consumption and C L,rms reduce significantly.Thus, these leading and trailing curvatures enhance the oscillating hydrofoil's overall propulsive performance and stability.Therefore, it is evident that the combination of the rounded leading edge and tapered tail (compare the NACA0012 airfoil and thick plate) provides better thrust and lower power consumption and C L,rms , providing overall enhanced propulsion.For the role of hydrofoil thickness (compare the thick and thin plates), a greater thickness reduces the net thrust and enhances the power consumption and C L,rms , resulting in degradation of the overall propulsive performance and stability of an oscillating hydrofoil.Finally, for the role of the phase difference between heave and pitch, it was found that both the maximum thrust and efficiency occur over a narrow range 90 • ≤ φ ≤ 100 • independent of the hydrofoil geometry.
The hydrodynamic reason behind the above variation is explained by correlating both the time-varying and time-instantaneous vortex structures in the vicinity of the hydrofoil with the corresponding engineering parameters.It is found that any amount of flow separation enhances the power consumed by the hydrofoil for performing its oscillation motion.In contrast, a moderate flow separation enhances the thrust generated.For all the different-shaped hydrofoils examined, the maximum efficiency occurs at those α max at which there is a small degree of flow separation with no roll-up of separating vortex.In comparison, maximum thrust generation occurs when there is a moderately strong flow separation with no induction of a significant quantity of fluid around the trailing edge.For any α max , flow separation can be controlled by the geometrical parameters of the hydrofoil.A rounded leading edge and flattened tail reduce the magnitude of flow separation, while a bluff leading edge, tapered tail, and increased thickness enhance it.
We anticipate that by contributing to an understanding of the correlation between different hydrofoil geometrical parameters with engineering parameters and flow separation, our results can aid in designing and developing more efficient modern autonomous underwater and micro air vehicles.Our results also have implications for understanding fish swimming.A fish possesses a symmetrical caudal fin with a somewhat rounded leading edge and tapered tail, and maintains a phase difference of 90 • while swimming.The results from the present study show that their NACA-like caudal fin generates high thrust, as high as that for a very thin plate, with minimum power requirement, resulting in higher propulsive efficiency.Moreover, their propulsion is very stable.A greater thickness without losing propulsive performance can give them a caudal fin with more structural strength, which is essential for directing water around it.The present study also shows that a phase difference of 90 • for any shape provides high thrust with high efficiency and stability of oscillation.Thus, an oscillating caudal fin of NACA-like shape with a 90 • phase difference between heave and pitch provides fish with an ideal combination for achieving high thrust, efficiency, and stability for propulsion.
In the current preliminary investigation, we have compared the hydrodynamics of hydrofoils with smooth, rounded leading edges and tapered tails to those with sharp leading and trailing edges.In the future, we plan to conduct a more comprehensive comparison using a non-dimensional geometric parameter that can assess the roundness and tapering of any hydrofoil shape to develop a flow law for optimization of thrust and propulsive efficiency.

Data Availability Statement:
The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest:
The authors declare no conflict of interest.

Figure 1 .
Figure 1.The different body shapes considered to understand the effects on propulsion and efficiency of curvature and thickness of hydrofoil-like geometries.

Figure 3 .
Figure 3. Non-dimensional computational setup along with boundary conditions for a coupled pitching-and-heaving NACA0012 hydrofoil in freestream flow.Here, U b and V b are the local velocity components of the hydrofoil-shaped body.

Figure 4 .
Figure 4. Validation of the computer code for an oscillating NACA0012 hydrofoil.(a) Grid independence study.Temporal variation of the lift coefficient C L at α max = 30 • , h • /C = 0.75, St = 0.25 and Re = 5000.(b) Code verification study comparing predictions against those of Boiron et al. [38].Temporal variation in thrust coefficient C T and lift coefficient C L at α max = 20 • , h • /C = 0.75, St = 0.3 and Re = 4000.The intermediate size grid, 768 × 552, is used.

Figure 5 .
Figure 5. Variation of the time-averaged (a) thrust coefficient C T , (b) power coefficient C P , (c) propulsive efficiency η P , and (d) rms of lift coefficient C L,rms , as a function of maximum angle of attack α max for the five different shapes considered.

Figure 6
Figure 6 also shows the normal direction in which the hydrodynamic force acts from the fluid to the hydrofoil.The axial component of this force provides the instantaneous • , in contrast to τ/T = 0.1.The change in trend is due to the flow remaining attached up to α max = 20 • (Figure 7(c1-c4)), a small separated flow at α max = 25 • (Figure 7(c5)), and a larger separation for α max = 30 • (Figure 7(c6)).

Figure 10 .
Figure 10.Variation of normalized (a) thrust coefficient C T , and (b) efficiency η p , with the phase difference between heave and pitch φ for cases at η p,max .

Figure 11 .
Figure 11.Time-instantaneous vortex structure (at maximum heave position) in the vicinity of the different hydrofoils at various φ for cases at η p,max .