# Fluid-Structure Interaction Effects on the Propulsion of an Flexible Composite Monofin

Authors: Adil El Baroudi; Fulgence Razafimahery
Journal: Journal of Engineering
 (2014)

---

## Abstract

Finite element method has been used to analyze the propulsive efficiency of a swimming fin. Fluid-structure interaction model can be used to study the effects of added mass on the natural frequencies of a multilayer anisotropic fin oscillating in a compressible fluid. Water by neglecting viscidity effects has been considered as a surrounding fluid and the frequency response of the fin has been compared with that of vacuum conditions. It has been shown that because of the added mass effects in water environment, the natural frequencies of the fin decrease.


---

## Body

## 1. Introduction

Multilayer anisotropic structure has wide applications in areas such as modern construction engineering, biomechanical engineering, aerospace industries, aircraft construction, and the components of nuclear power plants. It is therefore very important that the modal and dynamic analysis of multilayer anisotropic structure when subjected to different loading conditions be clearly understood so that they may be safely used in these industrial applications.It is well known that the natural frequencies of structures in contact with fluid are different from those in vacuum. Therefore, the prediction of natural frequency changes due to the presence of the fluid is important for designing structures which are in contact with or immersed in fluid. In general, the effect of the fluid force on the structure is represented as added mass, which lowers the natural frequency of the structure from that which would be measured in a vacuum. This decrease in the natural frequency of the fluid-structure system is caused by increasing the kinetic energy of the coupled system without a corresponding increase in strain energy.In this paper the propulsive efficiency of a swimming fin has been studied. Dynamic analysis of aquatic locomotion is a fundamental parameter in the performance search. In the case of swimming with fins, the propulsive efficiency depends on several factors. Most models suggested aim to evaluate the dynamic performances, including drag and lift which are the two relevant parameters relevant to quantifing the propulsive efficiency of a fin. Some proposed models are essentially of discrete type [1, 2], while others, by being inspired by organs of propulsion of marine cetaceans, use continuous models [3, 4]. Most of these authors do not account for the highly coupled nature of the problem. In fact, for rate of stresses observed in actual swimming, the coupling between the fluid and the fin becomes stronger.


## 2. Governing Equations

The numerical formulations used include the displacement formulation [5], the potential formulation [6], the pressure formulation [7], and the combination of some of them [8]. Finite element method is used to extract the natural frequencies and modal shapes. To compute the natural vibration modes of a fluid alone, the fluid is typically described either by pressure or by displacement potential variables. When the fluid is coupled with a solid, standard methods to solve (1) and (2) consist in eliminating either the pressure or the displacement potential [9]. However, in both cases nonsymmetric eigenvalue problems are obtained (see, e.g., [10]). To avoid this drawback, Morand and Ohayon introduce in [6] an alternative procedure which consists in using pressure and displacement potential simultaneously. In this section we summarize their approach; further details and discussions can be found in their book [11].In this work, we assume an amateur swimmer, where the scale of velocityU0 is supposed to be very small compared to the compression wave velocities cL in the fin. Indeed, some amateur swimmers note that, when making foot movements at low frequency, the resonance phenomenon and buckling phenomena appear. And we cannot explain why these phenomena tend to occur, because the natural frequencies of the fin which would be measured in the vacuum are higher than of the beat frequency of an ankle, for example. In this study, we assume that the swimmer does not disturb the free surface of the fluid domain. This leads to neglect the gravity effects.Dimensional analysis of coupled equations (Navier-Stokes equations and governing equations of nonlinear elasticity) of a fluid-structure interaction model [12] reveals several dimensionless parameters. One of its dimensionless parameters, α=U0/cL, is called displacement parameter. The displacement parameter α allows characterizing the nature of the coupling problem considered in this work. In the amateur swimmer hypothesis where U0≪cL, we can set the parameter α to a very low value and we can show that the convective terms and viscosity terms can be neglected for the fluid model [12]. We can also assume the assumption of small deformations for the fin. The resulting model is called inertial coupling [13]. The real shape of the fin is given in Figure 1, but for the sake of the simplicity the problem is bidimensional (Figure 2) and the fin is immersed in a large pool. The fin is modeled by a multilayer linear elastic transverse anisotropic material. The different layers constituting the fin are denoted by Ωi and have the density ρi. We denote by ui the displacement field in the fin and p the pressure field in the fluid. c0 and ρ0 denote the sound celerity and density of the fluid, respectively. The longitudinal axis of the fin is denoted by x. The force F given in (1) is used to describe the motion of the fin. θi denotes the orientation of fibers relative to the longitudinal axis x on the fin and takes 0∘ or 90∘. Here, each layer is made of either fiberglass or carbon fiber.Figure 1
Geometry of 3D fin.Figure 2
2D computational domain.The use of the ALE method is not essential in this study because the material is assumed linear. In the frame attached to the fin, the problem is to find(ui,p) solutions:(i)
solid domain(Ωi):
(1)ρi∂2ui∂t2=∇·σui+ρiF,σui=Kθiɛui;(ii)
fluid domain(Ωf):
(2)1ρ0c02∂2p∂t2=∇·1ρ0∇p-ρ0F;(iii)
fluid-solid interaction(Γ):
(3)σun=-pn,∇p-ρ0F·n=-ρ0∂2u∂t2·n;(iv)
other boundary conditions:(4)u=0Γ0,∇p-ρ0F·n=0Γf,p=0Γe∪Γg∪Γs.


## 3. Modal Analysis of Composite Monofin

Modal analysis of elastic submerged structures is needed in every modern construction and has wide engineering application especially in ocean engineering. In this study, modal analysis is important to predict the dynamic behavior of the submerged fin. It is well known that the natural frequencies of the submerged elastic structures are different from those in vacuum. The effect of fluid forces on the submerged fin is represented as added mass, which decreases the natural frequencies of the submerged fin from those which would be measured in the vacuum. This decrease in the natural frequencies of the submerged structures is caused by the increase of the kinetic energy of the fluid-fin system without a corresponding increase in strain energy. This step seems important to calculate the variations of natural frequencies of the fin for different situations. For this, we looked at the modes of the fin in the vacuum and water. Indeed, to test the quality of a fin, it is usual to search its quasi-static deformed shape and dynamic response in air. The aim is to check if results of tests carried out of the water are strongly influenced by the presence of the surrounding fluid. In addition, frequencies can have accurate information in the dynamic behavior of the system. By introducing the spaces of test functionsV={v∈H1(Ωs),v=0(Γ0)} and ϕ∈Q=H1(Ωf), the weak formulation of (1) and (2) holds:
(5)∫Ωs‍σ(u):ɛ(v)dx-ω2∫Ωs‍ρu·vdx+∫Γ‍pv·ndΓ=0,∫Ωf‍1ρ0∇p·∇ϕdx-ω2∫Ωf‍pϕρ0c02dx-ω2∫Γ‍u·nϕdΓ=0,
where
(6)∫Ωs‍σu:ɛvdx=∑i=1NL‍∫Ωi‍σui:ɛvidx,∫Ωs‍ρu·vdx=∑i=1NL‍∫Ωi‍ρiui·vidx.NL is a number of layers. Using Lagrange finite elements, where uh∈P2×P2 and ph∈P1, discretization of the weak formulation (5) induces a nonsymmetrical system:
(7)KsBOKpUP=ω2MsOMaMpUP,
where U and P are the vectors of nodal values for u and p, respectively. The submatrices of (7) are defined by
(8)VTKsU=∫Ωs‍σu:ɛvdx,ΦTKpP=∫Ωf‍pϕc02dx,VTMsU=∫Ωs‍ρu·vdx,VTBP=∫Γ‍pv·ndΓ,ΦTMpP=∫Ωf‍∇p·∇ϕdx,ΦTMaU=∫Γ‍ρ0u·nϕdΓ,
where V and Φ are the vectors of nodal values for v and ϕ, respectively. Ma is the added mass matrix (symmetric and positive definite [11]). The nonsymmetric system (7) was solved using the commercial software Comsol Multiphysics.Two types of calculations were carried out. The first is when the palm is plunged into the vacuum and the second into water. We give below the results for a model of up to five layers (NL=5) and the natural frequencies in vacuum and water. The fibers of each layer are arranged alternately along the two directions orthogonal axes x and y of the mean plane of the fin. Tables 2 and 4 show that arrangement of layers has a strong influence on the natural frequencies, and the added mass decreases the natural frequencies. Figures 3 and 4 show that arrangement of layers has no influence on the coupled modal shapes (Tables 1 and 3).Table 1
Natural frequenciesf (Hz) of the fin in the case of 0°/90°/0°/90°/0°.
Oriented fibers
In vacuum
f
1
f
2
f
3
0°/90°
12.6
76.0
205.3
0°/90°/0°
29.3
122.9
265.0
0°/90°/0°/90°
22.6
90.1
192.0
0°/90°/0°/90°/0°
28.8
79.7
185.1Table 2
Natural frequenciesf (Hz) of the fin in the case of 0°/90°/0°/90°/0°.
Oriented fibers
In water
f
1
f
2
f
3
0°/90°
1.7
46.0
103.1
0°/90°/0°
3.5
21.7
55.1
0°/90°/0°/90°
2.4
14.5
37.3
0°/90°/0°/90°/0°
2.7
11.8
33.4Table 3
Natural frequenciesf  (Hz) of the fin in the case of 90°/0°/90°/0°/90°.
Oriented fibers
In vacuum
f
1
f
2
f
3
90°/0°
12.6
73.5
158.4
90°/0°/90°
9.5
51.3
105.5
90°/0°/90°/0°
22.1
56.2
115.5
90°/0°/90°/0°/90°
17.4
44.8
91.6Table 4
Natural frequenciesf (Hz) of the fin in the case of 90°/0°/90°/0°/90°.
Oriented fibers
In a water
f
1
f
2
f
3
90°/0°
1.7
13.3
39.4
90°/0°/90°
1.1
8.7
23.8
90°/0°/90°/0°
2.2
9.5
23.8
90°/0°/90°/0°/90°
1.6
6.8
17.2The colors pertain to the pressure field and the arrows to the velocity field in the case of0∘/90∘/0∘/90∘/0∘.
(a)
Coupled mode 1 (f1=2.7 (Hz))
(b)
Coupled mode 2 (f2=11.8 (Hz))The colors pertain to the pressure field and the arrows to the velocity field in the case of90∘/0∘/90∘/0∘/90∘.
(a)
Coupled mode 1 (f1=1.6 (Hz))
(b)
Coupled mode 2 (f2=6.8 (Hz))


## 4. Dynamic Analysis of Composite Monofin

The dynamic problem was conducted using the data proposed in [14]. For this, the fin is subjected to combined translational and rotation motions. In this case, the quantity F introduced in the model problem (1) has the expression
(9)F=xω˙2(t)+yω¨(t)-h¨(t)sin[ω(t)]yω˙2(t)-xω¨(t)-h¨(t)cos⁡[ω(t)],
where
(10)ωt=θ0sin2πft,ht=h0sin2πft-ψ,θ0=40∘, ψ=π/2, h0=1c,  f=0.225[Hz], and c=0.7 is the chord of the profile, that is to say, the length of the fin. The phase ψ is introduced to model the muscle dissymmetry.To avoid a resonant frequency, the excitation frequency is taken far enough from the first natural frequency of the coupled system. The hydrodynamic parameters that seem most relevant are the total forceR=∫Γ‍σ(u)ndΓ exerted on the fin during the movement phase. The two components of R are, respectively, the drag (D) and lift (L) of the fin. The quantity T=-D is called thrust. Different types of layers exist in the manufacture of fins. Throughout the model the thickness of the fin is fixed in advance. We use the same physical characteristics as in the case of modal analysis. Using the same notation as before, the weak formulation of boundary value problem (1) and (2) is then written:
(11)d2dt2∫Ωs‍ρu·vdx+∫Ωs‍σ(u):ɛ(v)dx+∫Γ‍pv·ndΓ=-∫Ωs‍ρF·vdx,d2dt2∫Ωf‍pϕc02dx+∫Γ‍ρ0u·nϕdΓ+∫Ωf‍∇p·∇ϕdx=∫Ωf‍ρ0F·∇ϕdx.
In this section, we use a particular kinematics proposed in [14, 15], even if our models are not exactly similar. Indeed, the kinematics will allow us in the future to develop a new experimental protocol for measuring various hydrodynamic parameters of a fin. As the model problem (1)-(2) is linear, it is interesting to see the different contributions of each elementary movement in the dynamic response of the fin.


### 4.1. Dynamic Response in the Case of Translational Motion

The rotationω(t) is canceled and the movement is then sinusoidal along the direction y. Figure 5 shows that the two-layer model seems to give a greater thrust than the other models. This is consistent with the results of modal analysis, where this is the first model that has the lowest frequency. This type of movement is not interesting for the propulsive efficiency. Indeed, it leads to a zero mean propulsive efficiency. On the other hand, we see a greater amplitude for thrust, compared to the lift.Thrust of fin in the case of translation motion.
(a)
0
∘
/
9
0
∘
/
0
∘
/
9
0
∘
/
0
∘
(b)
9
0
∘
/
0
∘
/
9
0
∘
/
0
∘
/
9
0
∘


### 4.2. Dynamic Response in the Case of Rotation Motion

The functionh(t) is canceled and the movement will be a sinusoidal rotation around foot. According to Figure 7, the two-layer model always gives a greater thrust than the other models. But by eliminating this model in the response curves, we can see that the five-layer model gives the best performance. The three-layer model gives a better lift compared to other models. Thus, this type of movement provides a propulsive efficiency rather interesting. This phenomenon is also well observed in the movement of marine mammals. On the other hand, such movement can be interesting if you want to stay stationary at one position.


### 4.3. Dynamic Response in the Case of Combined Rotational-Translational Motion

In order to have a reasonable performance of the system, we must combine both translational and rotation motion and take the full expression of the excitation forceF. According to Figure 9, the two-layer model always gives a greater thrust than the other models. In general, the three-layer model seems to give a better compromise. Indeed, its thrust remains positive all the time, while its lift is of negative value and has less importance than other models. It is possible that by varying some physical parameters, we can significantly reduce hydrodynamic quantities, such as the moment and lift.


## 4.1. Dynamic Response in the Case of Translational Motion

The rotationω(t) is canceled and the movement is then sinusoidal along the direction y. Figure 5 shows that the two-layer model seems to give a greater thrust than the other models. This is consistent with the results of modal analysis, where this is the first model that has the lowest frequency. This type of movement is not interesting for the propulsive efficiency. Indeed, it leads to a zero mean propulsive efficiency. On the other hand, we see a greater amplitude for thrust, compared to the lift.Thrust of fin in the case of translation motion.
(a)
0
∘
/
9
0
∘
/
0
∘
/
9
0
∘
/
0
∘
(b)
9
0
∘
/
0
∘
/
9
0
∘
/
0
∘
/
9
0
∘


## 4.2. Dynamic Response in the Case of Rotation Motion

The functionh(t) is canceled and the movement will be a sinusoidal rotation around foot. According to Figure 7, the two-layer model always gives a greater thrust than the other models. But by eliminating this model in the response curves, we can see that the five-layer model gives the best performance. The three-layer model gives a better lift compared to other models. Thus, this type of movement provides a propulsive efficiency rather interesting. This phenomenon is also well observed in the movement of marine mammals. On the other hand, such movement can be interesting if you want to stay stationary at one position.


## 4.3. Dynamic Response in the Case of Combined Rotational-Translational Motion

In order to have a reasonable performance of the system, we must combine both translational and rotation motion and take the full expression of the excitation forceF. According to Figure 9, the two-layer model always gives a greater thrust than the other models. In general, the three-layer model seems to give a better compromise. Indeed, its thrust remains positive all the time, while its lift is of negative value and has less importance than other models. It is possible that by varying some physical parameters, we can significantly reduce hydrodynamic quantities, such as the moment and lift.


## 5. Conclusions

Finally, the above results allow us to draw some conclusions.(i)
The presence of layers provides some flexibility as indicated by the results of modal analysis. The first mode is flexural type, which justifies the use of models proposed in [1].(ii)
Fins with anisotropic material structures allow implementing a technique of layers parameterization to improve performance. It is quite possible now to bring special attention to the structure of the layers and types of constituent materials thereof.(iii)
The sensitivity of the dynamic behavior of the model with respect to the materials used and the boundary conditions for the fluid domain should be noted. Indeed, the presence or absence of rigid walls alters significantly the natural modes of the coupled system. Thus, the dynamic behavior of a swimmer depends on the localization in the pool where it is at the given moment. To obtain a better thrust, the fin has to be elastic and has to be sought at least in rotation. The amplitude of the vertical translation must be controlled to avoid a too high lift, in order to remain at a constant depth. The use of multilayer fins allows controlling an excessive variation of lift (Figures6, 8, and 10).(iv)
Most experimental results we know [16, 17] are mainly interested in kinematic aspect of the mechanical problem. Nevertheless, the results obtained in the case of a rigid fin [14] allowed us to have a basis for comparison. We found that dynamic responses curves are similar but with different amplitudes. These differences in results obtained can be explained by the type of models used (rigid fin and flexible composite fin).Lift of fin in the case of translation motion.
(a)
0
∘
/
9
0
∘
/
0
∘
/
9
0
∘
/
0
∘
(b)
9
0
∘
/
0
∘
/
9
0
∘
/
0
∘
/
9
0
∘Thrust of fin in the case of rotation motion.
(a)
0
∘
/
9
0
∘
/
0
∘
/
9
0
∘
/
0
∘
(b)
9
0
∘
/
0
∘
/
9
0
∘
/
0
∘
/
9
0
∘Lift of fin in the case of rotation motion.
(a)
0
∘
/
9
0
∘
/
0
∘
/
9
0
∘
/
0
∘
(b)
9
0
∘
/
0
∘
/
9
0
∘
/
0
∘
/
9
0
∘Thrust of fin in the case of combined rotational and translational motion.
(a)
0
∘
/
9
0
∘
/
0
∘
/
9
0
∘
/
0
∘
(b)
9
0
∘
/
0
∘
/
9
0
∘
/
0
∘
/
9
0
∘Lift of fin in the case of combined rotational and translational motion.
(a)
0
∘
/
9
0
∘
/
0
∘
/
9
0
∘
/
0
∘
(b)
9
0
∘
/
0
∘
/
9
0
∘
/
0
∘
/
9
0
∘In this paper the modal and dynamic analysis is proposed to understand the behavior of a flexible composite fin with a good accuracy. The publications in the literatures deal with the behavior of fins; few authors have not studied the case of coupled boundary conditions [18]. It is for this reason this study was conducted on the request of a company specializing in the design of fins in order to improve the propulsion of a flexible composite fin.



---

*Source: 541953-2014-12-14.xml*