EDP Sciences Journals List
Advanced Search
Subscriber Authentication Point
Home All issues Volume 14 Article
Free access article

Issue ESAIM: Proc.
Volume 14, 2005
CEMRACS 2004 - Mathematics and applications to biology and medicine
Page(s) 115 - 123
DOI http://dx.doi.org/10.1051/proc:2005010

ESAIM: Proc., September 2005, Vol. 14, pp. 115-123
DOI: 10.1051/proc:2005010

A penalty method for the simulation of fluid - rigid body interaction

João Janela1, Aline Lefebvre2 and Bertrand Maury2

1  CEMAT, Instituto Superior Técnico and Departamento de Matemática, ISEG, Lisbon, Portugal;
2  Laboratoire de Mathématiques, Université Paris Sud, Orsay, France;

jjanela@iseg.utl.pt
Aline.Lefebvre@math.u-psud.fr
Bertrand.Maury@math.u-psud.fr

Abstract
We present here a method to simulate the motion of a rigid body in a fluid. The method is based on a variational formulation on the whole fluid/solid domain, with some constraints on the unknown and the test functions. These constraints are relaxed by introducing a penalty term, which leads to a minimization problem over unconstrained functional spaces. This makes the method straightforward to implement from any finite element Stokes/Navier-Stokes solver. It is shown that, as the penalty parameter goes to infinity, we recover the coupled fluid-solid equations. We apply this approach to a simplified 2D model of the aortic valve.


Résumé

Nous présentons ici une méthode de simulation du mouvement d'un corps rigide dans un fluide. Cette méthode est basée sur une formulation variationnelle sur tout le domaine fluide/solide, avec des contraintes sur l'inconnue et sur les fonctions test. On relaxe ces contraintes en introduisant un terme de pénalisation, qui conduit à un problème de minimisation sur des espaces fonctionnels non contraints. Ainsi, tout solveur éléments finis pour Stokes/Navier Stokes permet de programmer aisément cette méthode. On montre que, quand le paramètre de pénalisation tend vers l'infini, on retrouve le système d'équations couplées fluide/solide. Cette approche est appliquée à un modèle 2D simplifié de la valve aortique.


Mathematics Subject Classification. 35Q30, 65M25, 65M60

Key words: Navier-Stokes equation, rigid motion, finite element, penalty method, valve


© EDP Sciences, ESAIM 2005


What is OpenURL?

The OpenURL standard is a protocol for transmission of metadata describing the resource that you wish to access. An OpenURL link contains article metadata and directs it to the OpenURL server of your choice. The OpenURL server can provide access to the resource and also offer complementary services (specific search engine, export of references...). The OpenURL link can be generated by different means.
  • If your librarian has set up your subscription with an OpenURL resolver, OpenURL links appear automatically on the abstract pages.
  • You can define your own OpenURL resolver with your EDPS Account. In this case your choice will be given priority over that of your library.
  • You can use an add-on for your browser (Firefox or I.E.) to display OpenURL links on a page (see http://www.openly.com/openurlref/). You should disable this module if you wish to use the OpenURL server that you or your library have defined.