Numerical simulation of a pulsatile flow through a flexible channel

Laboratoire de Mathématiques, Informatique et Applications, Université de Haute-Alsace, 4, rue des Frères Lumière, 68093 Mulhouse Cedex, France; cornel.murea@uha.fr

Received: 15 January 2005
Revised: 1 June 2006
Revised: 15 November 2006

Abstract

An algorithm for approximation of an unsteady fluid-structure interaction problem is proposed. The fluid is governed by the Navier-Stokes equations with boundary conditions on pressure, while for the structure a particular plate model is used. The algorithm is based on the modal decomposition and the Newmark Method for the structure and on the Arbitrary Lagrangian Eulerian coordinates and the Finite Element Method for the fluid. In this paper, the continuity of the stresses at the interface was treated by the Least Squares Method. At each time step we have to solve an optimization problem which permits us to use moderate time step. This is the main advantage of this approach. In order to solve the optimization problem, we have employed the Broyden, Fletcher, Goldforb, Shano Method where the gradient of the cost function was approached by the Finite Difference Method. Numerical results are presented.