In the framework of the small perturbations theory, we study the incompressible inviscid flow of a uniform stream past an oscillatory/undulatory thin hydrofoil including floor effects. A Green function is used to deduce the integral equation for the jump of the pressure past the foil. The integral equation is numerically solved and the average drag coefficient is calculated. For some wings there appears a propulsive force and this force increases when the hydrofoil is close to the floor.
MSC: 76B10, 65R20, 45H99, 31A10.
Keywords:flexible hydrofoil; hypersingular integral equation; floor effect; thrust
In the present paper we study the small-amplitude oscillatory/undulatory motion of an incompressible fluid past a thin flexible plate which performs prescribed oscillations in the presence of a wall (floor) with which it is parallel in its undisturbed state, and relative to which it is moving with constant speed. We shall limit the analysis to bodies large enough so that the Reynolds number is large. As is stated by Eloy et al. in , when the flexible surface has a typical speed of several body lengths per second, the flow can be considered irrotational, meaning that the flow vorticity is concentrated in thin boundary layers adjacent to the body surface and in a thin wake (vortex sheet) behind the body. Since the effects of viscosity manifest inside the thin boundary layers, we may treat the fluid as inviscid in the rest of the flow domain. Recalling Lagrange-Cauchy’s theorem which states that if the flow is potential in a certain configuration, it remains potential in every configuration arising from the initial one, we deduce that the theory of the unsteady motion of lifting wings as well as the theory of potential flow can be successfully utilized (see the papers of Carabineanu [2-6], Dowell and Hall , Dragoş , Homentcovschi [9,10], Lighthill , Street , Taylor , Wu [14,15], Watkins et al.).
The periodic motion of a flexible foil is oscillatory if the foil or parts of it remain rigid during the motion. The undulatory motion involves a traveling wave down the foil (Street ). As we know from aerodynamics and hydrodynamics studies (Dragoş and Carabineanu [17,18], Dragoş et al.), the hydrodynamic coefficients of a hydrofoil are influenced by the presence of the floor. The aim of the paper is to predict the drag or the thrust enhancement generated by the presence of the floor. We employ like in [2-6,8,9] the linearized Euler equations for the incompressible flow. For taking into account the floor effect we use the Green function of the Laplacean for the Neumann problem in the half-space. We use the integral representation for the harmonic functions and the slipping condition to obtain an integral equation for the jump of the pressure over the hydrofoil. In order to discretize the integral equation, we split the kernel of the equation into several kernels for which we provide appropriate approximation formulas depending on the type of singularity of the kernel. Assuming that the hydrofoil is subjected to harmonic oscillations, we simplify the integral equation making it independent of time. By solving the discretized integral equation we calculate the jump of the pressure over the wing.
After obtaining the pressure field, the average drag is calculate by performing a numerical integration. We study an example of undulatory motion of the flexible thin delta wing. When the frequency surpasses a critical value, the drag becomes negative i.e. it appears a propulsive force. We notice that the distance between the wing and the floor influences the drag and the thrust.
The statement of the problem
We consider the continuity and Euler equations for incompressible flow in a fixed Cartesian frame of reference having the versors i, j, k. At the moment t, the hydrofoil (see Figure 1) has the equation:
Figure 1. The flexible hydrofoil and the vortex sheet.
and we obtain
We consider the floor
where the slipping condition is also imposed:
Let represent the thin vortex wake behind the hydrofoil . In the small-amplitude approximation theory the wake remains planar (see the demonstration of Homentcovschi in ). We linearize the equations of motion around the rest state (neglecting the products of the perturbation quantities) in the domain where and we have
where p is the pressure and is the constant density of the fluid. The aim of the present paper is to use the boundary conditions (7), (9), and the partial differential equations (10) for obtaining an integral equation for the jump of the pressure across . In order to ensure the uniqueness of the solution we shall impose a certain periodic in time behavior to the unknowns.
The Green function. The integral equation of the problem
Eliminating u, v, w, from (10) we get the Laplace equation
Obviously on the floor F
We assume that the harmonic function p vanishes at infinity and we have the integral representation formula
where stands for the outward normal derivative and . Here and are, respectively, the upper face and the lower face of and similarly for . As usual in the small perturbation theory, we replace with and we assume that is planar (as Homentcovschi stated in ). Hence
From the boundary condition (7) and from the fourth equation (10) we deduce that
From (9) it follows that
Taking into account that
we get from the fourth equation of (10) and from (18)
In order to obtain a domain of integration that does not vary in time, we shall introduce (like in  and ) and a new system of coordinates , related to the lifting wing. We consider the Galilean transformations
We also denote
In the new system of coordinates, the equation of the flexible hydrofoil is
and the perturbation velocity is
In (27) may have complex values. By convention means the real part of . For the sake of simplicity, we shall calculate (as is usual in the oscillatory hydrofoil theory) complex values for the jump of pressure and then we shall consider the real part.
We also have
From (7), (26), and (27) we get
We reuse the notations which must not be confounded with the notations for the dimensional variables corresponding to the fixed frame . Let . The velocity field (with respect to the fixed frame) is , where is the perturbation velocity of the fluid. Introducing the dimensionless functions and variables (reduced frequency), , , , and taking into account the linearized slip condition, the integral equation (30) becomes
We discretize the hypersingular integral equation (32) in order to solve it numerically. In Appendix A we split the kernel into several kernels and describe the type of singularity for each one. In Appendix B, depending on the kind of singularity, we deliver appropriate approximation formulas. In order to ensure the uniqueness of the solution, we impose a certain behavior of the pressure jump in the vicinity of the leading edge.
The propulsive force. Numerical results. Floor effects
We introduce the pressure jump coefficient
We consider the undulatory delta hydrofoil whose equation is
whence, using nondimensional coordinates,
We used the values , and the nondimensional distances to floor and . In Figure 2 we present (as three-dimensional surfaces with contour plots beneath the surface) the flexible hydrofoil and the pressure jump coefficient fields (divided by 2α) for the nondimensional moments where is the nondimensional period of oscillation.
Figure 2. Pressure jump past the flexible hydrofoil.
We are also interested in calculating the drag coefficient
We introduce and we calculate the average drag coefficient:
In Figure 3 we present the average drag coefficient (divided by ) against the reduced frequency for the delta hydrofoil in ground effects. We consider , . We notice that if the reduced frequency surpasses a certain value, the average drag coefficient is negative, i.e. it appears a propulsive force. We also notice that the propulsive force is bigger for (dash dot line) than for (continuous line) i.e. it is bigger when the oscillatory wing is closer to the floor.
Figure 3. Average drag coefficient for the oscillatory hydrofoil.
Appendix 1: The singularities of the kernel of the integral equation
For solving numerically the integral equation (32) we have to discretize the left hand member in order to obtain a linear algebraic system of equations. To this aim we split, like in  and , the kernel
into several kernels in order to show the kind of singularities we are dealing with and to find afterwards the most convenient approximation formulas.
We have step by step
The integrals from the right hand part of (43) represent the sine and cosine Fourier transforms of and in  one shows that
The integral equation becomes
The kernels and have strong singularities of order . The kernel has a polar singularity. The kernels and have integrable logarithmic singularities. Taking into account the series expansions of , , , , , we may easily prove that the kernels , , , , , have no singularity and they are continuous functions. , , are also continuous functions. We notice that for .
Appendix 2: The discretization of the integral equation
We consider the undulatory delta hydrofoil. The equations of the leading edge are
In order to ensure the uniqueness of the solution of the integral equation, some analytical results from  suggest to presume that there exists a continuous finite function g such that . We have therefore
where FP stands for the finite part of the hypersingular integral as it is introduced by Ch. Fox in . Since the inner integral vanishes for , we assume that
where is finite for . We consider on D a net consisting of the nodes (grid points, control points) , , . For the hypersingular integral occurring in (65) we may use the quadrature formula for equidistant control points given by Dumitrescu ,
whence it follows that
Finally we deduce
At last we find
where , . and are integrals which are evaluated numerically with the trapezoidal rule. For calculating the Bessel (MacDonald) functions and we may utilize the series expansions. We may also utilize the libraries offered by MATLAB. For calculating the kernels and we use the integral representations
The integrals are evaluated numerically with the trapezoidal rule. The approximation formulas for the kernels , , , , , , , , , , and were also given in  and . The new approximation formulas for the kernels , , are given for the first time herein. Denoting , we obtain, discretizing the two-dimensional integral equation (32):
The authors declare that they have no competing interests.
All authors read and approved the final manuscript.
Eloy, C, Doaré, O, Duchemin, L, Scouveiller, L: A unified introduction to fluid mechanics of flying or swimming at high Reynolds number. Exp. Mech. (2009). Publisher Full Text
Carabineanu, A: Incompressible flow past oscillatory wings of low aspect ratio by the integral equations method. Int. J. Numer. Methods Eng.. 45, 1187–1201 (1999). Publisher Full Text
Carabineanu, A: Self-propulsion of oscillating wings in incompressible flow. Int. J. Numer. Methods Fluids. 56, 1–21 (2008). Publisher Full Text
Carabineanu, A: Self-propulsion of aquatic animals by undulatory or oscillatory motion. In: Stelian I, Marinoschi G (eds.) Proc. 3rd Workshop on Math. Modeling of Environmental and Life Sciences, Constanta, May 2004. Series on Mathematical Modelling of Environmental and Life Sciences Problems. pp. 131–145. Ed. Acad. Romane, Bucharest (2004)
Carabineanu, A: Self-propulsion of an oscillatory wing including ground effects. In: Stelian I, Marinoschi G, Popa C (eds.) Proc. 4th Workshop on Math. Modeling of Environmental and Life Sciences, Constanta, September 2005. Series on Mathematical Modelling of Environmental and Life Sciences Problems. pp. 39–54. Ed. Acad. Romane, Bucharest (2004)
Dowell, EH, Hall, KC: Modeling of fluid-structure interaction. Annu. Rev. Fluid Mech.. 33, 445–490 (2001). Publisher Full Text
Dragoş, L: The theory of oscillating thick wings in subsonic flow. Lifting line theory. Acta Mech.. 54, 221–238 (1985). Publisher Full Text
Homentcovschi, D: Theory of the lifting surface in unsteady motion in an inviscid fluid. Acta Mech.. 27, 205–216 (1977). Publisher Full Text
Homentcovschi, D: Determination of the vortex associated with the motion of a body in an ideal fluid. Mech. Res. Commun.. 3(3), 191–195 (1976). Publisher Full Text
Street, A: Preliminary Finite Element Modeling of a Piezometric Actuated Marin Propulsion Fish. Thesis, Department of Mechanical Engineering, Rochester Institute of Technology, Rocester, New York (2006)
Wu, TY: Mathematical biofluiddynamics and mechanophysiology of fish locomotion. Math. Methods Appl. Sci.. 24, 1541–1564 (2001). Publisher Full Text
Dragoş, L, Carabineanu, A: A numerical solution for the equation of the lifting surface in ground effects. Commun. Numer. Methods Eng.. 18, 177–187 (2002). Publisher Full Text
Dragoş, L, Carabineanu, A: The supersonic flow past a thin profile including ground and tunnel effects. Z. Angew. Math. Mech.. 82, 649–652 (2002). Publisher Full Text
Dragoş, L, Carabineanu, A, Dumitrache, R: A numerical solution for the equation of the lifting line including ground and tunnel effects. Proc. Rom. Acad., Ser. A: Math. Phys. Tech. Sci. Inf. Sci.. 9(3), 1–7 (2008)
Fox, C: A generalisation of the Cauchy principal value. Can. J. Math.. 9, 110–115 (1957). Publisher Full Text