Skip to main content

A Green function approach for the investigation of the incompressible flow past an oscillatory thin hydrofoil including floor effects

Abstract

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: Primary 76B10; secondary 65R20; 45H99; 31A10.

Introduction

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 [1], 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 [26], Dowell and Hall [7], Dragoş [8], Homentcovschi [9, 10], Lighthill [11], Street [12], Taylor [13], Wu [14, 15], Watkins et al. [16]).

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 [12]). As we know from aerodynamics and hydrodynamics studies (Dragoş and Carabineanu [17, 18], Dragoş et al. [19]), 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 [26, 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 Oxyz having the versors i, j, k. At the moment t, the hydrofoil S(t) (see Figure 1) has the equation:

zh(x+ V 0 t,y,t)=0,(x,y,0)D(t).
(1)

D(t) is the projection of the hydrofoil onto the Oxy plane. We assume that

| h x |1,| h y |1,| h t |1
(2)

and that the hydrofoil moves into the Ox direction, producing small perturbations v=ui+vj+wk of the vanishing velocity of the surrounding fluid.

Figure 1
figure 1

The flexible hydrofoil and the vortex sheet.

Let V 0 i be the average velocity of the hydrofoil S(t) i.e. the velocity of D(t). In a point of D(t) we have

x(t)= x V 0 t,y(t)= y , ( x , y ) D =D(0).
(3)

The velocity of an arbitrary point of (x,y,z)S(t) is

V= V 0 i+ d z d t k= V 0 i+ d h ( x ( t ) , y , t ) d t k
(4)

and the normal vector on S(t) is

n= h x i h y j+k.
(5)

We linearize the slipping condition on the two sides of S(t)

( v V )n=0,
(6)

and we obtain

w ( x ( t ) , y , ± 0 , t ) = V 0 h x ( x ( t ) , y , t ) + d h d t ( x ( t ) , y , t ) , ( x ( t ) , y , 0 ) D(t).
(7)

We consider the floor

F= { ( x , y , χ ) ; ( x , y ) R 2 , χ > 0 } ,
(8)

where the slipping condition is also imposed:

w(x,y,χ,t)=0.
(9)

Let Σ(t) represent the thin vortex wake behind the hydrofoil S(t). In the small-amplitude approximation theory the wake remains planar (see the demonstration of Homentcovschi in [10]). We linearize the equations of motion around the rest state (neglecting the products of the perturbation quantities) in the domain Ξ={(x,y,z,t);tR,(x,y,z)Ω(t)} where Ω(t)={(x,y,z) R 3 ;zχ}(S(t)Σ(t)) and we have

u x + v y + w z = 0 , u t + 1 ρ 0 p x = 0 , v t + 1 ρ 0 p y = 0 , w t + 1 ρ 0 p z = 0 ,
(10)

where p is the pressure and ρ 0 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 S(t). 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

Δp=0.
(11)

Let us consider the Green function of the Laplacean for the Neumann problem in the half-space z>χ:

G ( x , y , z ; x , y , z ) = 1 4 π | ( x x ) i + ( y y ) j + ( z z ) k | 1 4 π | ( x x ) i + ( y y ) j + ( z + z + 2 χ ) k | .
(12)

Obviously on the floor F

G ( x , y , χ ; x , y , z ) z =0.
(13)

We assume that the harmonic function p vanishes at infinity and we have the integral representation formula

p ( x , y , z , t ) = Ω ( t ) [ p ( x , y , z , t ) G ( x , y , z ; x , y , z ) n p ( x , y , z , t ) n G ( x , y , z ; x , y , z ) ] d σ ,
(14)

where n stands for the outward normal derivative and Ω(t)=F S + (t) S (t) Σ + (t) Σ (t). Here S + (t) and S (t) are, respectively, the upper face and the lower face of S(t) and similarly for Σ(t). As usual in the small perturbation theory, we replace S(t) with D(t) and we assume that Σ(t) is planar (as Homentcovschi stated in [10]). Hence

n | D + ( t ) Σ + ( t ) = z , n | D ( t ) Σ ( t ) = z , n | F = z .
(15)

From the boundary condition (7) and from the fourth equation (10) we deduce that

p n | D + ( t ) Σ + ( t ) D ( t ) Σ ( t ) =0.
(16)

From (9) it follows that

p n | F =0.
(17)

Since the pressure is continuous over Σ(t), from (9)-(14) we obtain

p(x,y,z,t)= D ( t ) [ p ( x , y , + 0 , t ) p ( x , y , 0 , t ) ] G ( x , y , z ; x , y , z ) z d x d y .
(18)

Taking into account that

z 1 | ( x x ) i + ( y y ) j + ( z z ) k | = z 1 | ( x x ) i + ( y y ) j + ( z z ) k | ,
(19)
z 1 | ( x x ) i + ( y y ) j + ( z + z + 2 χ ) k | = z 1 | ( x x ) i + ( y y ) j + ( z + z + 2 χ ) k | ,
(20)

we get from the fourth equation of (10) and from (18)

w ( x , y , z , t ) = 1 4 π t D ( t ) f ( x , y , t ) 2 z 2 [ 1 | ( x x ) i + ( y y ) j + ( z z ) k | 1 | ( x x ) i + ( y y ) j + ( z + z + 2 χ ) k | ] d x d y d t ,
(21)

where f(x,y,t)= p ( x , y , + 0 , t ) p ( x , y , 0 , t ) ρ 0 .

In order to obtain a domain of integration that does not vary in time, we shall introduce (like in [5] and [6]) and a new system of coordinates O x y z , related to the lifting wing. We consider the Galilean transformations

x =x+ V 0 t, y =y, z =z, t =t.
(22)

We also denote

χ =χ, x = x + V 0 t , y = y , λ = x x V 0 ( t t ) .
(23)

In the new system of coordinates, the equation of the flexible hydrofoil is

z = h ( x , y , t ) =h(x+ V 0 t,y,t)
(24)

and the perturbation velocity is

v ( x , y , t ) =v(x+ V 0 t,y,t), v = u i+ v j+ w k.
(25)

Obviously,

h t = d h d t .
(26)

Let D be the projection of the lifting surface onto the O x y plane. Considering the lifting surface subject to harmonic oscillations, we impose

h ( x , y , t ) = h ( x , y ) exp ( i ω t ) .
(27)

In (27) h ( x , y ) may have complex values. By convention h ( x , y )exp(iω t ) means the real part of h ( x , y )exp(iω t ). 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

f ( x , y , t ) =f ( x , y ) exp ( i ω t ) .
(28)

From (7), (26), and (27) we get

w ( x , y , ± 0 , t ) = V 0 h x exp ( i ω t ) +iω h ( x , y ) exp ( i ω t ) .
(29)

Performing the change of variables ( x , y , t )( x , y , λ ) and considering z 0, from (28) and (29) it follows that

V 0 h x + i ω h ( x , y ) = 1 4 π V 0 D f ( x , y ) exp ( i ω x 0 V 0 ) × [ x 0 exp ( i ω λ V 0 ) ( 1 ( λ 2 + y 0 2 ) 3 / 2 + 8 χ 2 λ 2 y 0 2 ( λ 2 + y 0 2 + 4 χ 2 ) 5 / 2 ) d λ ] d x d y
(30)

with x 0 = x x , y 0 = y y . The sign indicates the finite part in the Hadamard sense of the integral. Denoting by a the half-span of the wing, we introduce the dimensionless variables

(x,y,z,λ,ξ,η,χ,h,t)= ( x a , y a , z a , λ a , x a , y a , χ a , h a , V 0 t a ) .
(31)

We reuse the notations (x,y,z,λ,χ,h,t) which must not be confounded with the notations for the dimensional variables corresponding to the fixed frame Oxyz. Let D={(x,y,0);(ax,ay,0) D }. The velocity field (with respect to the fixed frame) is V= V 0 i+v, where v=ui+vj+wk is the perturbation velocity of the fluid. Introducing the dimensionless functions and variables ω ˜ = ω a V 0 (reduced frequency), f ˜ (x,y)= f ( a x , b y ) V 0 2 , x 0 =xξ, y 0 =yη, and taking into account the linearized slip condition, the integral equation (30) becomes

h ( x , y ) x + i ω ˜ h ( x , y ) = 1 4 π D f ˜ ( ξ , η ) exp ( i ω ˜ x 0 ) × [ x 0 exp ( i ω ˜ λ ) ( 1 ( λ 2 + y 0 2 ) 3 / 2 + 8 χ 2 λ 2 y 0 2 ( λ 2 + y 0 2 + 4 χ 2 ) 5 / 2 ) d λ ] d ξ d η .
(32)

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

C p (x,y,t)=2Re [ f ˜ ( x , y ) exp ( i ω t ) ] .
(33)

We consider the undulatory delta hydrofoil whose equation is

0= z α x exp ( i ω 1 x + i ω t ) ; ( x , y ) D
(34)

whence, using nondimensional coordinates,

h(x,y)=αxexp(i ω ˜ 1 x), ω ˜ 1 =a ω 1 ;(x,y)D.
(35)

We used the values ω ˜ =π/6, ω ˜ 1 =2π/3 and the nondimensional distances to floor χ=1 and χ=10. 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 t{0, T 3 , 2 T 3 } where T= 2 π ω ˜ = 2 π V 0 a ω is the nondimensional period of oscillation.

Figure 2
figure 2

Pressure jump past the flexible hydrofoil.

We are also interested in calculating the drag coefficient

C D (t)= D n x C p (x,y,t)dxdy.
(36)

Since

n x =αRe [ ( 1 + i ω ˜ 1 x ) exp ( i ω ˜ 1 x + i ω ˜ t ) ] ,
(37)

we deduce

C D (t)=2α D Re [ ( 1 + i ω ˜ 1 x ) exp ( i ω ˜ 1 x + i ω ˜ t ) ] Re [ f ˜ ( x , y ) exp ( i ω ˜ t ) ] dxdy.
(38)

We introduce and we calculate the average drag coefficient:

C ˜ D = 1 T 0 T C D ( t ) d t = α D { Re [ f ˜ ( x , y ) exp ( i ω ˜ 1 x ) ] + ω ˜ 1 x Im [ f ˜ ( x , y ) exp ( i ω ˜ 1 x ) ] } d x d y .
(39)

In Figure 3 we present the average drag coefficient (divided by α 2 ) against the reduced frequency for the delta hydrofoil in ground effects. We consider ω ˜ =π/6, ω ˜ 1 =π/2. 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 χ=1 (dash dot line) than for χ=10 (continuous line) i.e. it is bigger when the oscillatory wing is closer to the floor.

Figure 3
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 [5] and [6], the kernel

N(x,y;ξ,η)= x 0 exp(i ω ˜ λ) ( 1 ( λ 2 + y 0 2 ) 3 / 2 + 8 χ 2 λ 2 y 0 2 ( λ 2 + y 0 2 + 4 χ 2 ) 5 / 2 ) dλ
(40)

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

x 0 exp ( i ω ˜ λ ) ( λ 2 + y 0 2 ) 3 / 2 dλ= x 0 exp ( i ω ˜ λ ) 1 ( λ 2 + ϖ 2 y 0 2 ) 3 / 2 dλ+ 1 y 0 2 ( 1 + x 0 x 0 2 + y 0 2 ) ,
(41)
x 0 exp ( i ω ˜ λ ) 1 ( λ 2 + y 0 2 ) 3 / 2 dλ= 0 x 0 exp ( i ω ˜ λ ) 1 ( λ 2 + y 0 2 ) 3 / 2 dλ+ 0 exp ( i ω ˜ λ ) 1 ( λ 2 + y 0 2 ) 3 / 2 dλ,
(42)
0 exp ( i ω ˜ λ ) 1 ( λ 2 + y 0 2 ) 3 / 2 dλ= 1 y 0 2 + 0 cos ω ˜ λ ( λ 2 + y 0 2 ) 3 / 2 dλi 0 sin ω ˜ λ ( λ 2 + y 0 2 ) 3 / 2 dλ.
(43)

The integrals from the right hand part of (43) represent the sine and cosine Fourier transforms of ( λ 2 + y 0 2 ) 3 / 2 and in [20] one shows that

0 cos ω ˜ λ ( λ 2 + y 0 2 ) 3 / 2 dλ= ω ˜ | y 0 | K 1 ( ω ˜ | y 0 | ) ,
(44)
0 sin ω ˜ λ ( λ 2 + y 0 2 ) 3 / 2 dλ= π 2 ω ˜ | y 0 | ( L 1 ( ω ˜ | y 0 | ) I 1 ( ω ˜ | y 0 | ) ) .
(45)

L 1 , L 2 are Strouve functions and I 1 , I 2 , K 1 , K 2 are Bessel functions. We also have

0 x 0 exp ( i ω ˜ λ ) 1 ( λ 2 + y 0 2 ) 3 / 2 d λ = 0 x 0 exp ( i ω ˜ λ ) 1 i ω ˜ λ + ω ˜ 2 λ 2 / 2 ( λ 2 + y 0 2 ) 3 / 2 d λ + i ω ˜ | y 0 | i ω ˜ ( x 0 2 + y 0 2 ) 1 / 2 + ω ˜ 2 x 0 2 ( x 0 2 + y 0 2 ) 1 / 2 ω ˜ 2 2 ln ( x 0 + x 0 2 + y 0 2 ) + ω ˜ 2 2 ln | y 0 | ,
(46)
x 0 exp ( i ω ˜ λ ) 8 χ 2 λ 2 y 0 2 ( λ 2 + y 0 2 + 4 χ 2 ) 5 / 2 d λ = 0 x 0 exp ( i ω ˜ λ ) 8 χ 2 λ 2 y 0 2 ( λ 2 + y 0 2 + 4 χ 2 ) 5 / 2 d λ + 12 χ 2 0 cos ( ω ˜ λ ) ( λ 2 + y 0 2 + 4 χ 2 ) 5 / 2 d λ 12 i χ 2 0 sin ( ω ˜ λ ) ( λ 2 + y 0 2 + 4 χ 2 ) 5 / 2 d λ 0 cos ( ω ˜ λ ) ( λ 2 + y 0 2 + 4 χ 2 ) 3 / 2 d λ + i 0 sin ( ω ˜ λ ) ( λ 2 + y 0 2 + 4 χ 2 ) 3 / 2 d λ ,
(47)
0 cos ( ω ˜ λ ) ( λ 2 + y 0 2 + 4 χ 2 ) 3 / 2 dλ= ω ˜ y 0 2 + 4 χ 2 K 1 ( ω ˜ y 0 2 + 4 χ 2 ) ,
(48)
0 cos ( ω ˜ λ ) ( λ 2 + y 0 2 + 4 χ 2 ) 5 / 2 dλ= 2 ω ˜ 3 y 0 2 + 4 χ 2 K 2 ( ω ˜ y 0 2 + 4 χ 2 ) ,
(49)
0 sin ( ω ˜ λ ) ( λ 2 + y 0 2 + 4 χ 2 ) 3 / 2 d λ = ω ˜ π 2 y 0 2 + 4 χ 2 [ L 1 ( ω ˜ y 0 2 + 4 χ 2 ) I 1 ( ω ˜ y 0 2 + 4 χ 2 ) ] ,
(50)
0 sin ( ω ˜ λ ) ( λ 2 + y 0 2 + 4 χ 2 ) 5 / 2 d λ = ω ˜ 2 π 6 ( y 0 2 + 4 χ 2 ) [ I 2 ( ω ˜ y 0 2 + 4 χ 2 ) L 2 ( ω ˜ y 0 2 + 4 χ 2 ) ] .
(51)

From the previous calculations we deduce that we can split the kernel N(x,y;ξ,η) as follows:

N(x,y;ξ,η)= N 1 (x,y;ξ,η)++ N 14 (x,y;ξ,η),
(52)

where

N 1 = 1 y 0 2 ( x 0 x 0 2 + y 0 2 x 0 | x 0 | ) , N 2 = 1 y 0 2 ( 1 + x 0 | x 0 | ) , N 3 = i ω ˜ x 0 2 + y 0 2 ,
(53)
N 4 = ω ˜ 2 2 ln ( | y 0 | ) ( 1 + x 0 | x 0 | ) , N 5 = ω ˜ 2 2 x 0 | x 0 | ln ( | x 0 | + x 0 2 + y 0 2 ) ,
(54)
N 6 = ω ˜ 2 x 0 2 x 0 2 + y 0 2 , N 7 = ω ˜ | y 0 | K 1 ( ω ˜ | y 0 | ) 1 y 0 2 ω ˜ 2 2 ln ω ˜ | y 0 | 2 + ω ˜ 2 2 ln 1 2 ,
(55)
N 8 (x,y;ξ,η)= 0 x 0 exp ( i ω ˜ u ) 1 i ω ˜ u + ω ˜ 2 u 2 / 2 ( u 2 + y 0 2 ) 3 / 2 du,
(56)
N 9 (x,y;ξ,η)= 0 x 0 exp(i ω ˜ u) 8 χ 2 u 2 y 0 2 ( u 2 + y 0 2 + 4 χ 2 ) 5 / 2 du,
(57)
N 10 = ω ˜ y 0 2 + 4 χ 2 K 1 ( ω ˜ y 0 2 + 4 χ 2 ) ,
(58)
N 11 = 8 ω ˜ χ 2 y 0 2 + 4 χ 2 K 2 ( ω ˜ y 0 2 + 4 χ 2 ) ,
(59)
N 12 (x,y;ξ,η)= ω ˜ π i 2 y 0 2 + 4 χ 2 [ L 1 ( ω ˜ y 0 2 + 4 χ 2 ) I 1 ( ω ˜ y 0 2 + 4 χ 2 ) ] ,
(60)
N 13 (x,y;ξ,η)= 2 ω ˜ 2 π i χ 2 ( y 0 2 + 4 χ 2 ) [ L 2 ( ω ˜ y 0 2 + 4 χ 2 ) I 2 ( ω ˜ y 0 2 + 4 χ 2 ) ] ,
(61)
N 14 (x,y;ξ,η)= i π ω ˜ 2 2 | y 0 | ( I 1 ( ω ˜ | y 0 | ) L 1 ( ω ˜ | y 0 | ) + 2 π ) .
(62)

The integral equation becomes

1 4 π i = 1 14 D f ˜ ( ξ , η ) exp ( i ω ˜ ξ ) N i ( x , y ; ξ , η ) d ξ d η = ( h ( x , y ) x + i ω ˜ h ( x , y ) ) exp ( i ω ˜ x ) .
(63)

The kernels N 1 and N 2 have strong singularities of order 1 y 0 2 . The kernel N 3 has a polar singularity. The kernels N 4 and N 5 have integrable logarithmic singularities. Taking into account the series expansions of L 1 , L 2 , I 1 , I 2 , K 1 , K 2 we may easily prove that the kernels N 7 , N 10 , N 11 , N 12 , N 13 , N 14 have no singularity and they are continuous functions. N 6 , N 8 , N 9 are also continuous functions. We notice that N 2 = N 4 =0 for x 0 <0.

Appendix B: The discretization of the integral equation

We consider the undulatory delta hydrofoil. The equations of the leading edge are

y ± (x)=± x b ;x[0,b].
(64)

In order to ensure the uniqueness of the solution of the integral equation, some analytical results from [2] suggest to presume that there exists a continuous finite function g such that f ˜ (ξ,η)= g ( ξ , η ) ξ 2 b 2 η 2 . We have therefore

D f ˜ ( ξ , η ) exp ( i ω ˜ ξ ) N 1 ( x , y ; ξ , η ) d ξ d η = F P 1 1 1 y 0 2 ( b | η | b g ( ξ , η ) ξ 2 b 2 η 2 exp ( i ω ˜ ξ ) ( x 0 x 0 2 + y 0 2 x 0 | x 0 | ) d ξ ) d η ,
(65)

where FP stands for the finite part of the hypersingular integral as it is introduced by Ch. Fox in [21]. Since the inner integral vanishes for η=1, we assume that

b | η | b g ( ξ , η ) ξ 2 b 2 η 2 exp ( i ω ˜ ξ ) ( x 0 x 0 2 + y 0 2 x 0 | x 0 | ) d ξ = 1 η 2 G ( x , y ; η ) ,
(66)

where G(x,y;η) is finite for η=±1. We consider on D a net consisting of the nodes (grid points, control points) ( x h , y ¯ j )=( h n b, 2 j + 1 2 n ), h=1,,n, j=h,h+1,,h1. For the hypersingular integral occurring in (65) we may use the quadrature formula for equidistant control points given by Dumitrescu [22],

FP 1 1 1 η 2 G ( x k , y ¯ l ; η ) ( y ¯ l η ) 2 dη= j = n n 1 G( x k , y ¯ l ; y ¯ j ) A l j ,
(67)
A l j = arccos ( y j ) A l j = + arccos ( y j + 1 ) + 1 y j 2 y j y ¯ l 1 y j + 1 2 y j + 1 y ¯ l y ¯ l 1 y ¯ l 2 ln | C l ( j + 1 ) C l j | ,
(68)
C l j = 1 y j 1 + y ¯ l 1 + y j 1 y ¯ l 1 y j 1 + y ¯ l + 1 + y j 1 y ¯ l .
(69)

We shall give a quadrature formula for calculating G( x k , y ¯ l ; y ¯ j ). Denoting

x ¯ h j = { x h b 2 n , h < j < h 1 , x h b 4 n , j { h , h 1 } ,
(70)
x ˜ h j = { x h b n , h < j < h 1 , x h b 2 n , j { h , h 1 } ,
(71)
g h j =g( x ¯ h j , y ¯ j ),H(j)= { 1 , j 0 , 0 , j < 0 ,
(72)

we have

1 y ¯ j 2 G ( x k , y ¯ l ; y ¯ j ) b | y ¯ j | b g ( ξ , y ¯ j ) ξ 2 b 2 y ¯ j 2 exp ( i ω ˜ ξ ) ( x 0 x 0 2 + y 0 2 x 0 | x 0 | ) d ξ h = | j | + H ( j ) n [ g h j exp ( i ω ˜ x ¯ h j ) ( x k x ¯ h j ( x k x ¯ h j ) 2 + ( y ¯ l y ¯ j ) 2 x k x ¯ h j | x k x ¯ h j | ) × x ˜ h j x i d ξ ξ 2 b 2 y ¯ j 2 ] ,
(73)

whence it follows that

G( x k , y ¯ l ; y ¯ j ) h = | j | + H ( j ) n g h j B h j k l , B h j k l = E h j D h j k l 1 y ¯ j 2 ,
(74)
E h j =exp(i ω ˜ x ¯ h j ) [ ln ( x h + x h 2 b 2 y ¯ j 2 ) ln ( x ¯ h j + x ˜ h j 2 b 2 y ¯ j 2 ) ] ,
(75)
D h j k l = ( x k x ¯ h j ( x k x ¯ h j ) 2 + ( y ¯ l y ¯ j ) 2 x k x ¯ h j | x k x ¯ h j | ) ,h<j<h1.
(76)

Finally we deduce

D f ˜ ( ξ , η ) exp ( i ω ˜ ξ ) N 1 ( x k , y ¯ l ; ξ , η ) d ξ d η h = 1 n j = h h 1 g h j N h j k l ( 1 )
(77)

with

N h j k l ( 1 ) = A l j B h j k l .
(78)

For N 2 (x,y;ξ,η)= 1 y 0 2 (1+ x 0 | x 0 | ) we get

D f ˜ ( ξ , η ) exp ( i ω ˜ ξ ) N 2 ( x , y ; ξ , η ) d ξ d η = 2 F P x x 1 y 0 2 ( b | η | x g ( ξ , η ) exp ( i ω ˜ ξ ) ξ 2 b 2 η 2 d ξ ) d η .
(79)

Assuming that b | η | x k g ( ξ , η ) exp ( i ω ˜ ξ ) ξ 2 b 2 η 2 dξ= x k 2 / b 2 η 2 G ( k ) ( x k ;η), we have

D f ˜ ( ξ , η ) exp ( i ω ˜ ξ ) N 2 ( x k , y ¯ l ; ξ , η ) d ξ d η = 2 F P x k / b x k / b x k 2 / b 2 η 2 G ( k ) ( x k ; η ) ( y ¯ l η ) 2 d η 2 j = k k 1 G ( k ) ( x k ; y ¯ j ) A l j ( k ) ,
(80)
A l j ( k ) = arccos ( b y j x k ) + arccos ( b y j + 1 x k ) + x k 2 b 2 y j 2 b y j b y ¯ l x k 2 b 2 y j + 1 2 b y j + 1 b y ¯ l b y ¯ l x k 2 b 2 y ¯ l 2 ln | C l ( j + 1 ) ( k ) C l j ( k ) | ,
(81)
C l j ( k ) = x k b y j x k + b y ¯ l x k + b y j x k b y ¯ l x k b y j x k + b y ¯ l + x k + b y j x k b y ¯ l .
(82)

For calculating G ( k ) ( x k ; y ¯ j ) we employ the quadrature formula

G ( k ) ( x k ; y ¯ j ) = 1 x k 2 / b 2 y ¯ j 2 b | y ¯ j | x k g ( ξ , y ¯ j ) exp ( i ω ˜ ξ ) ξ 2 b 2 y ¯ j 2 d ξ = i = | j | + H ( j ) k g i j E i j x k 2 / b 2 y ¯ l 2 .
(83)

At last we find

D f ˜ (ξ,η)exp(i ω ˜ ξ) N 2 ( x k , y ¯ l ;ξ,η)dξdη= h = 1 n j = h h 1 g h j N h j k l ( 2 ) ,
(84)
N h j k l ( 2 ) = { 2 A l j ( k ) E h j x k 2 / b 2 y ¯ l 2 ; h k , 0 ; h > k .
(85)

The singularities of the kernels N 3 , N 4 , N 5 are weaker than 1 y 0 2 . We replace these kernels with N 3 y 0 2 y 0 2 , N 4 y 0 2 y 0 2 , N 5 y 0 2 y 0 2 for obtaining approximation formulas similar to the formulas for N 1 . We get for p=3,4,5:

D f ˜ ( ξ , η ) exp ( i ω ˜ ξ ) N p ( x k , y ¯ l ; ξ , η ) d ξ d η = i = 1 n j = i i 1 g i j N i j k l ( p ) ,
(86)
N i j k l ( p ) = { A l j ( y ¯ l y ¯ j ) 2 E i j 1 y ¯ j 2 N p ( x k , y ¯ l ; x ¯ i j , y ¯ j ) ; i j , 0 ; i = j .
(87)

The kernels N 6 , N 7 , N 8 , N 9 , N 10 , N 11 , N 12 , N 13 , and N 14 are continuous and we utilize the approximation formulas

D f ˜ ( ξ , η ) exp ( i ω ˜ ξ ) N p ( x k , y ¯ l ; ξ , η ) d ξ d η = h = 1 n j = h h 1 g h j N h j k l ( p ) , p = 6 , , 14 ,
(88)
N h j k l ( p ) = E h j N p ( x k , y ¯ l ; x ¯ h j , y ¯ j )/n,p=6,,14.
(89)

For calculating N 7 ( x k , y ¯ l ; x ¯ h j , y ¯ l ) we use the series expansions of the Bessel and Strouve functions and we take into account that

N 7 ( x k , y ¯ l ; x ¯ h j , y ¯ l )= ω ˜ 2 ( ψ ( 1 ) + ψ ( 2 ) ) 4 + π i ω ˜ 4 ,
(90)

where ψ(1)=0.5772, ψ(2)=0.4228. N 8 ( x k , y ¯ l ; x ¯ h j , y ¯ j ) and N 9 ( x k , y ¯ l ; x ¯ h j , y ¯ j ) are integrals which are evaluated numerically with the trapezoidal rule. For calculating the Bessel (MacDonald) functions K 1 and K 2 we may utilize the series expansions. We may also utilize the libraries offered by MATLAB. For calculating the kernels N 12 ( x k , y ¯ l ; x ¯ h j , y ¯ j ) and N 14 ( x k , y ¯ l ; x ¯ h j , y ¯ j ) we use the integral representations

L 1 (x) I 1 (x)= 2 exp ( x ) π + 2 x π 0 π / 2 exp(xcost) ( sin t sin 2 t ) dt,
(91)
L 2 ( x ) I 2 ( x ) = 2 π x + 2 x exp ( x ) 3 π + 2 x 2 3 π 0 π / 2 exp ( x cos t ) ( sin t sin 4 t ) d t .
(92)

The integrals are evaluated numerically with the trapezoidal rule. The approximation formulas for the kernels N 1 , N 2 , N 6 , N 7 , N 8 , N 9 , N 10 , N 11 , N 12 , N 13 , and N 14 were also given in [5] and [6]. The new approximation formulas for the kernels N 3 , N 4 , N 5 are given for the first time herein. Denoting N h j k l = N h j k l ( 1 ) + N h j k l ( 2 ) ++ N h j k l ( 14 ) , we obtain, discretizing the two-dimensional integral equation (32):

1 4 π h = 1 n j = h h 1 g h j N h j k l = ( h ( x k , y ¯ l ) x + i ω ˜ h ( x k , y ¯ l ) ) exp(i ω ˜ x k ).
(93)

References

  1. 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. 10.1007/s11340-009-9289-7

    Google Scholar 

  2. Carabineanu A: Incompressible flow past oscillatory wings of low aspect ratio by the integral equations method. Int. J. Numer. Methods Eng. 1999, 45: 1187-1201. 10.1002/(SICI)1097-0207(19990730)45:9<1187::AID-NME625>3.0.CO;2-Z

    Article  MathSciNet  Google Scholar 

  3. Carabineanu A: Self-propulsion of oscillating wings in incompressible flow. Int. J. Numer. Methods Fluids 2008, 56: 1-21. 10.1002/fld.1513

    Article  MathSciNet  Google Scholar 

  4. Carabineanu A: Self-propulsion of aquatic animals by undulatory or oscillatory motion. Series on Mathematical Modelling of Environmental and Life Sciences Problems. In Proc. 3rd Workshop on Math. Modeling of Environmental and Life Sciences. Edited by: Stelian I, Marinoschi G. Ed. Acad. Romane, Bucharest; 2004:131-145. Constanta, May 2004

    Google Scholar 

  5. Carabineanu A: Self-propulsion of an oscillatory wing including ground effects. Series on Mathematical Modelling of Environmental and Life Sciences Problems. In Proc. 4th Workshop on Math. Modeling of Environmental and Life Sciences. Edited by: Stelian I, Marinoschi G, Popa C. Ed. Acad. Romane, Bucharest; 2004:39-54. Constanta, September 2005

    Google Scholar 

  6. Carabineanu A: Self-propulsion of an oscillatory wing including tunnel effects. An. Univ. Bucur., Mat. 2006, LV: 5-26.

    MathSciNet  Google Scholar 

  7. Dowell EH, Hall KC: Modeling of fluid-structure interaction. Annu. Rev. Fluid Mech. 2001, 33: 445-490. 10.1146/annurev.fluid.33.1.445

    Article  Google Scholar 

  8. Dragoş L: The theory of oscillating thick wings in subsonic flow. Lifting line theory. Acta Mech. 1985, 54: 221-238. 10.1007/BF01184848

    Article  Google Scholar 

  9. Homentcovschi D: Theory of the lifting surface in unsteady motion in an inviscid fluid. Acta Mech. 1977, 27: 205-216. 10.1007/BF01180086

    Article  Google Scholar 

  10. Homentcovschi D: Determination of the vortex associated with the motion of a body in an ideal fluid. Mech. Res. Commun. 1976, 3(3):191-195. 10.1016/0093-6413(76)90009-4

    Article  MathSciNet  Google Scholar 

  11. Lighthill J: Mathematical Biofluiddynamics. SIAM, Philadelphia; 1987.

    Google Scholar 

  12. 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)

    Google Scholar 

  13. Taylor GI: Analysis of the swimming of long and narrow animals. Proc. R. Soc. Lond. Ser. A 1952, 214(117):158-183.

    Article  Google Scholar 

  14. Wu TY: Mathematical biofluiddynamics and mechanophysiology of fish locomotion. Math. Methods Appl. Sci. 2001, 24: 1541-1564. 10.1002/mma.218

    Article  MathSciNet  Google Scholar 

  15. Wu TY: On theoretical modeling of aquatic and aerial animal locomotion. 38. In Advances in Applied Mechanics. Edited by: Giessen E, Wu TY. Elsevier, Amsterdam; 2002:291-353.

    Google Scholar 

  16. Watkins, CE, Runyan, HL, Woolston, DS: On the Kernel function of the integral equation relating the lift and downwash distributions of oscillating finite wings in subsonic flow. NACA T.R. 1234 (1955)

    Google Scholar 

  17. Dragoş L, Carabineanu A: A numerical solution for the equation of the lifting surface in ground effects. Commun. Numer. Methods Eng. 2002, 18: 177-187. 10.1002/cnm.481

    Article  Google Scholar 

  18. Dragoş L, Carabineanu A: The supersonic flow past a thin profile including ground and tunnel effects. Z. Angew. Math. Mech. 2002, 82: 649-652. 10.1002/1521-4001(200209)82:9<649::AID-ZAMM649>3.0.CO;2-2

    Article  Google Scholar 

  19. 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. 2008, 9(3):1-7.

    Google Scholar 

  20. Ditkine V, Proudnikov A: Transformations intégrales et calcul opérationnel. Mir, Moscow; 1978.

    Google Scholar 

  21. Fox C: A generalisation of the Cauchy principal value. Can. J. Math. 1957, 9: 110-115. 10.4153/CJM-1957-015-1

    Article  Google Scholar 

  22. Dumitrescu DF: Three methods for solving Prandtl’s equation. Z. Angew. Math. Mech. 1996, 76(6):1-4.

    Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Adrian Carabineanu.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All authors read and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Răpeanu, E., Carabineanu, A. A Green function approach for the investigation of the incompressible flow past an oscillatory thin hydrofoil including floor effects. Bound Value Probl 2014, 104 (2014). https://doi.org/10.1186/1687-2770-2014-104

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2014-104

Keywords