Method of Successive Approximations for a Fluid Structure Interaction Problem

Abstract

In this paper, we present a method for solving coupled problem. This method is mainly based on the successive approximations method. The external force acting on the structure is replaced by λ = p (x1, H + u (x1, λ)). Then we have a nonlinear equation of unknown λ to solve by successive approximations method. By this method, we obtain easily the analytic expression of the displacement. In addition, good results are obtained with only a few iterations.

Share and Cite:

Sow, A. and Mbaye, I. (2014) Method of Successive Approximations for a Fluid Structure Interaction Problem. Applied Mathematics, 5, 2299-2304. doi: 10.4236/am.2014.515223.

1. Introduction

Problem involved in fluid structure interactions occurs in a wide variety of engineering problem and therefore has attracted the interest of many investigations from different engineering disciplines. As a result, much effort has gone into the development of general computational method for fluid structure system by Osses, Fernandez, quarteroni, Blouza, Mbaye [1] -[7] .

In this paper, successive approximations method is applied to solve a fluid-structure interaction problem. We replace the external force acting on the interface between fluid and structure by.

Then we introduce a nonlinear equation to solve by successive approximations such that the coupled problem is achievable.

By this method, we obtain good approximate solutions. In addition, the analytic solution of beam equation can be computed easily.

The fluid is modeled by two dimensional Stokes equations for steady flow and the structure is represented by the one dimensional beam equation.

2. Presentation of the Problem

We denote by the two dimensional domain occupied by the fluid, the elastic interface between fluid and structure and be the remaining external boundaries of the fluid as depicted in Figure 1.

defines the reference domain.

3. Governing Equations

3.1. Structure Equation

We start from the simple equation that governs the structure. The simplified beam equation is:

, (1)

. (2)

where is the displacement of the structure, is the external force of the structure, , is the young modulus, is the thickness of the structure This equation is good representation of the structure for small deformation.

3.2. Fluid Equation

We suppose that the fluid is governing by the Stokes equations for steady flow in:

, (3)

, (4)

, (5)

, (6)

, (7)

, (8)

, (9)

where denotes the fluid velocity, denotes the pressure,  denotes the volume force of the fluid,  the viscosity of the fluid and denotes the inflow velocity profile of the fluid, is the unit outward normal vector, is the identity matrix, on the symmetric axis we have the non-penetration condition.

4. Formulation of Coupled Problem

The problem is to find, and such that:

(10)

, (11)

, (12)

, (13)

, (14)

(15)

, (16)

, (17)

, (18)

We have a fluid structure interaction problem. The domain of the fluid depends on the displacement and the displacement depends on the velocity and the pressure of the fluid.

5. Successive Approximations Method

We assume that. Corresponding to each, we consider the coupled problem:

(19)

, (20)

, (21)

, (22)

, (23)

, (24)

, (25)

, (26)

, (27)

To solve this coupled problem, we need to solve a nonlinear equation of unknown define as: by the successive approximations method. Then we will find , and .

Description of the Method

The weak formulation of the fluid and the structure is given by Grandmont, Murea [8] [9] .

We summarize step by step our computational method to find such that:

(28)

, (29)

  Step 1: We give, the initial displacement and the fluid domain are compute.

  Step 2: We solve the Stokes equation by finite elements methods in the reference domain. We find   Step 3: we fine  Step 4: Do:

- - compute and- solve the stokes equation in, we find

- - we compute- While ()  Step 5: Give, ,.

6. Numerical Results

For each, the analytic solution of the beam equations is

.

We assume that the velocity on is such that and for all in Murea [10] .

The parameter values of the fluid and the structure are:

Parameter related to fluid: The fluid velocity is , the fluid density is , the channel length is , the channel width .

Parameter related to structure: The structure thickness , Young’s modulus is ,.

We use the Lagrange finite element to approach the velocities and Lagrange finite element is used to approach the pressure. FreeFem++ Hecht [11] is using for the numerical tests.

7. Conclusion

In this work, we applied successive approximations method to solve fluid structure interaction problem. This method gives good results when the displacement is small. After 11 iterations, we found a good approximate solution of the nonlinear equation and also we obtained the solution of coupled problem.

Table 1 and Figure 2 show that decreases to zero when iterations increase. Figure 3 and Figure 4 display the behavior of the fluid flow and the pressure wave respectively after 11.

Figure 1. Sets appearing to the fluid structure problem.

Figure 2. Error with respect to number iterations.

Figure 3. Fluid velocity and structure displacement.

Figure 4. Fluid pressure.

Table 1. Different values of after 11 iterations.

In a forthcoming work, we will be showed the theoretical convergence of and also successive approximations method will be used to solve an unsteady fluid structure interaction problem.

Conflicts of Interest

The authors declare no conflicts of interest.

 [1] Osses, A. and Puel, J.P. (2000) Approximate Controllability for a Linear Model of Fluid Structure Interactiont. ESAIM Control, Optimization and Calculus of Variation, 4, Article ID: 4976513. [2] Fernandez, M.A. and Moubachir, M. (2003) And Exact Block Newton Algorithm for Solving Fluid Structure Interaction Problem. Compte rendu de l’académie des sciences, paris, Serie I, 336, 681-686. [3] Fernandez, M.A. and Moubachir, M. (2005) A Newton Method Using Exact Jacobians for Solving Fluid-Structure Coupling. Computers and Structures, 83, 127-142. [4] Quateroni, A., Tuveri, M. and Veneziani, A. (2000) Computational Vascular Fluid Dynamics, Problem Models and Methods. Computing and Visualization in Science, 2, 163-197. http://dx.doi.org/10.1007/s007910050039 [5] Blouza, A., Dumas, L. and Mbaye, I. (2008) Multiobjective Optimization of a Stent in a Fluid-Structure Context. Proceedings of Genetic and Evolutionary Computation Conference (GECCO), 2056-2060.http://dx.doi.org/10.1145/1388969.1389021 [6] Mbaye, I. (2012) Controllability Approach for a Fluid-Structure Interaction Problem. Applied Mathematics, 3, 213-216.http://dx.doi.org/10.4236/am.2012.33034 [7] Mbaye, I. and Murea, C. (2008) Numerical Procedure with Analytic Derivative for Unsteady Fluid-Structure Interaction. Communications in Numerical Methods in Engineering, 24, 1257-1275. [8] Grandmont, C. (1998) Existence et Unicité de Solution d’un Problème de Couplage Fluide-Structure Bidimensionnel stationnaire. Compte rendu de l’académie des sciences, Paris. [9] Murea, C. and Maday, C. (1996) Existence of an Optimal Control for a Nonlinear Fluid Cable Interaction Problem. Rapport de Recherche CEMRACS, C.I.R.M., Luminy. [10] Murea, C. (2005) The BFGS Algorithm for Nonlinear Least Squares Problem Arising from Blood Flow in Arteries. Computers Mathematics with Applications, 49, 171-186. http://dx.doi.org/10.1016/j.camwa.2004.11.002 [11] Hecht, F. and Pironneau, O. (2003) A Finite Element Software for PDE FreeFem++.www.freefem.org/ff++/