Le delta du fleuve Sénégal est le théâtre de crues importantes, le plus souvent catastrophiques et la dernière en date a nécessité l'ouverture d'une brèche dans la Langue de Barbarie qui est une fine bande de sable séparant le delta du fleuve de la mer. Depuis, cette brèche ne cesse de s'agrandir sous l'action conjuguée des eaux du fleuve et de la mer. Modéliser ce phénomène d'élargissement nécessite de connaître à tout moment les caractéristiques de l'écoulement des eaux dans le delta, y compris en période de crue. Nous avons opté pour la réalisation de deux algorithmes de simulation. C'est ce travail que je présente ici. D'une part, nous considérons l'équation de Navier-Stokes tridimensionnelle et nous faisons l'hypothèse d'une faible épaisseur d'eau, ce qui est acceptable comparé aux dimensions du delta. Cette hypothèse nous permet de faire, en chaque point, une intégration selon la verticale pour obtenir une équation de Saint-Venant bidimensionnelle simplifiée dans un plan xoy sans les hypothèses approximatives (fond plat ou conditions aux limites trop généreuses dans l'un ou l'autre, qui ne reproduisent pas fidèlement le phénomène étudié). Cette équation ne nous permet d'obtenir les champs de vitesses et de pressions que dans un plan horizontal. Pour prendre en compte les apports extérieurs responsables d'une crue, nous introduisons une équation 1D de conservation de la masse d'eau. Le couplage entre l'équation de Saint-Venant 2D et l'équation de conservation 1D conduit à une modélisation (2D1/2 et non 3D) du phénomène étudié. D'autre part, la phase d'élargissement de la brèche résulte de la problématique fluide-structure, une intéraction entre un écoulement et une structure fortement déformable. Il s'agit d'un problème multi-physique reposant sur plusieurs équations d'état (Darcy, poroélastique ...) et impliquant divers couplages. Sous l'hypothèse que les déformations de la structure dues à l'hydrodynamique restent à l'échelle microscopique, nous obtenons des résultats analytiques et numériques de l'évolution de l'interface à long terme.

Cette thèse est motivée par la modélisation et la simulation numérique des phénomènes d'interaction fluide-structure autour de valves cardiaques. L'interaction avec la paroi des vaisseaux est traitée avec une formulation Arbitraire Lagrange Euler (ALE), tandis que l'interaction avec les valves est traitée à l'aide de multiplicateurs de Lagrange, dans une formulation de type Domaines Fictifs (FD). Après une présentation de synthèse des diverses méthodes utilisées en interaction fluide-structure dans les écoulements sanguins, nous décrivons une méthode permettant de simuler la dynamique d'une valve immergée dans un écoulement visqueux incompressible. L'algorithme de couplage est partionné, ce qui permet de conserver des solveurs fluides et structures indépendants. Le maillage du fluide est mobile pour suivre la paroi des vaisseaux, mais indépendant du maillage des valves. Ceci autorise des très grands déplacements sans nécessiter de remaillage. Nous proposons une stratégie pour gérer le contact entre plusieurs valves. L'algorithme est totalement indépendant des solveurs de structures et est bien adapté au couplage fluide-structure partionné. Enfin, nous proposons un schéma de couplage semi-implicite permettant de mêler efficacement les formulations ALE et FD. Toutes les méthodes considérées sont accompagnées de nombreux tests numériques en 2D et 3D.

Le propos de cette étude est la modélisation numérique de la guitare acoustique dans le domaine temporel. La méthode consiste en l'élaboration d'un modèle qui s'attache à décrire les phénomènes vibratoires et acoustiques mis en jeu depuis le pincer de corde jusqu'au rayonnement 3D du son. La corde est modélisée par une équation des ondes amortie 1D. Elle est couplée à la table d'harmonie via le chevalet. Le mouvement de la table est régi par le modèle de plaque mince amortie de Kirchhoff--Love pour un matériau orthotrope et hétérogène, percée d'un trou et encastrée sur son bord externe. Le reste du corps de la guitare (fond, bords, manche...) est supposé rigide. La table rayonne à l'intérieur et à l'extérieur de la cavité. La modélisation complète du champ acoustique rayonné est une approche originale comparativement aux études antérieures portant sur la guitare. On obtient un système d'équations aux dérivées partielles que l'on résout numériquement dans le domaine temporel. On utilise une méthode spectrale spécifique pour la résolution de l'équation de plaque dynamique de Kirchhoff-Love. Pour l'équation de corde et l'équation des ondes acoustiques, on utilise une méthode mixte standard pour l'approximation spatiale et des différences finies centrées en temps. Le problème d'interaction fluide-structure est résolu par une méthode de domaines fictifs qui permet d'approcher finement la géométrie de la guitare tout en utilisant un maillage cubique régulier pour le calcul du champ sonore 3D. L'originalité du schéma de résolution du modèle est un couplage stable entre une méthode de résolution exacte en temps et une méthode discrète. Un nombre important de simulations numériques est réalisées, montrant la validité de la méthode et les très riches potentialités d'une telle approche.

Les haubans de pont peuvent vibrer sous l'effet conjoint du vent et de la pluie (Rain-Wind-Induced Vibrations, RWIV). On observe toujours deux filets d'eau ruisselant le long du hauban sujet aux RWIV. De plus, une expérience récente (Alam et Zhou 2005), réalisée sur un cylindre fixe, montre que le sillage d'air est considérablement perturbé par la présence des filets. Nous présentons un modèle de film mince entourant un cylindre soumis au vent qui permet de prévoir l'apparition de filets, leur position et leur forme. Nous montrons en outre que l'existence des filets n'est possible que dans un intervalle de nombre de Weber, qui compare effet du vent et tension de surface. Un modèle d'interaction entre les filets et l'écoulement d'air environnant est ensuite exposé. Pour un intervalle de nombre de Weber, le détachement tourbillonnaire dévie de la loi de Strouhal et la portance fluctuante s'exerçant sur le cylindre est amplifiée, ce qui est en accord avec les expériences d'Alam. Pour les conditions de vent et de pluie qui président à l'apparition des RWIV, l'intervalle d'existence des filets et l'intervalle d'amplification de portance se recouvrent. Est-ce à dire que les RWIV peuvent être expliquées par un accrochage filets-sillage ?

Dans cette thèse on a développé des algorithmes semi-implicites procédures partagées et monolithiques pour l'interaction entre un fluide gouverné par le modèle de Navier Stokes et une structure. Dans le premier chapitre, on présente un algorithme semi-implicite procédures partagées pour l'interaction entre un fluide et une structure gouvernée soit par les équations d'élasticité linéaire ou soit par le modèle de Saint-Venant Kirchhoff non linéaire. Dans le second chapitre, on propose un algorithme semi-implicite procédures partagées pour l'interaction entre un fluide et une structure de modèle linéaire et on montre un résultat de stabilité inconditionnelle en temps de l'algorithme. Un problème d'optimisation est résolu dans les deux algorithmes précédents, afin de satisfaire les conditions de continuité des vitesses et d'égalité des contraintes à l'interface. Durant les itérations de BFGS pour résoudre le problème d'optimisation, le maillage fluide reste fixe et la matrice fluide n'est factorisée qu'une seule fois, ce qui réduit l'effort de calcul. Dans le troisième chapitre, un algorithme semi-implicite monolithique pour l'interaction entre un fluide et une structure de modèle linéaire est proposé. L'algorithme utilise un maillage global pour le domaine fluide structure. La condition de continuité des vitesses à l'interface est automatiquement satisfaite et celle de l'égalité des contraintes n'apparaît pas explicitement dans la formulation faible. A chaque pas de temps on résout un système monolithique d'inconnues vitesse et pression définies sur le domaine global. Le temps CPU est réduit quand l'approche monolithique est utilisée à la place des procédures partagées.

Dans cette thèse nous nous intéressons a l'approximation numérique des équations incompressibles de Navier-Stokes évoluant dans un domaine en mouvement par la méthode des éléments spectraux et des intégrateurs en temps d'ordre élève. Dans une première phase, nous présentons la méthode des éléments spectraux et les outils de base pour effectuer des discrétisations spectrales du type Galerkin ou Galerkin avec intégration numérique (G-NI). Nous couvrons un large éventail de possibilités concernant les éléments de référence, fonctions de base, points d'interpolation et points de quadrature. Dans cette approche, l'intégration et la différentiation des fonctions polynomiales est faite numériquement grâce a l'aide d'ensembles de points convenables. En ce qui concerne la différenciation, nous présentons une étude numérique des points qui doivent être utilisés pour atteindre une meilleure stabilité numérique (parmi les choix que nous avons actuellement). Deuxièmement, nous introduisons les équations incompressibles stationnaires et non-stationnaires de Stokes et de Navier-Stokes et son approximation spectrale. Dans le cas non-stationnaire, nous introduisons une combinaison de la méthode Backward Différentiation Formula (BDF) et une formule d'extrapolation du même ordre pour l'intégration par rapport au temps. Une fois les équations discrétisées, un système linéaire doit être résolu pour obtenir la solution approchée. Dans ce contexte, nous resolvons ce système avec un préconditionneur par blocs. Nous montrons que le préconditionneur est optimal par rapport au nombre d'iterations utilisées par la méthode GMRES dans le cas stationnaire, mais pas dans le cas non-stationnaire. Une autre alternative est d'utiliser les méthodes de factorization algébrique de type Yosida et séparer le calcul de la vitesse et de la pression. Un cas test est présente pour déterminer les proprietes de convergence de ce type de méthodes dans notre contexte. Troisièmement, nous 'tendons les algorithmes développés dans le cas ou le domaine est fixe au cadre de la formulation Arbitraire Lagrange-Euler (ALE). La question de la définition d'une carte ALE d'ordre élevé est aborée. Cela permet de construire un domaine de calcul qui est d'ecrit avec des éléments courbes. Un cas test utilisant une méthode directe et les méthodes Yosida-q pour resoudre le système linéaire est présente pour montrer les ordres de convergence de la méthode proposée. Finalement, nous appliquons la méthode développée pour résoudre une un problème d'interaction fluide-structure pour un exemple simple bidimensionnel d'hémodynamique. Nous considérons deux approches: une implicite entièrement couplée et une semi-implicite.

Ce travail est consacré à l'étude de problèmes de couplage fluide-structure par des méthodes de frontière immergée et de pénalisation.<br />Dans une première partie nous abordons les différents aspects de la simulation numérique de modèles couplant un fluide et une structure (formulation de l'écoulement, localisation de la structure et prise en compte des interactions), puis nous présentons les modèles sur lesquels nous nous concentrons dans la suite: fluide incompressible/membrane élastique et fluide incompressible/solide rigide.<br />Dans la deuxième partie nous étudions la stabilité numérique d'une formulation level-set de la méthode de frontière immergée pour le couplage fluide/membrane élastique. Une nouvelle condition de stabilité pour ce modèle est dérivée par analyse d'un modèle linéarisé 1D, puis validée numériquement en 2D.<br />Enfin la troisième partie est consacrée à l'analyse numérique d'une méthode de pénalisation pour le couplage fluide/solide rigide, dans laquelle la vitesse rigide est calculée par projection. Nous montrons la convergence du problème pénalisé vers une formulation faible du couplage, et illustrons ce résultat sur des tests numériques en sédimentation.

Cette thèse est motivée par la modélisation et la simulation numérique des phénomènes d'interaction fluide-structure autour de valves cardiaques. L'interaction avec la paroi des vaisseaux est traitée avec une formulation Arbitraire Lagrange Euler (ALE), tandis que l'interaction avec les valves est traitée à l'aide de multiplicateurs de Lagrange, dans une formulation de type Domaines Fictifs (FD). Après une présentation de synthèse des diverses méthodes utilisées en interaction fluide-structure dans les écoulements sanguins, nous décrivons une méthode permettant de simuler la dynamique d'une valve immergée dans un écoulement visqueux incompressible. L'algorithme de couplage est partitionné, ce qui permet de conserver des solveurs fluides et structures indépendants. Le maillage du fluide est mobile pour suivre la paroi des vaisseaux, mais indépendant du maillage des valves. Ceci autorise des très grands déplacements sans nécessiter de remaillage. Nous proposons une stratégie pour gérer le contact entre plusieurs valves. L'algorithme est totalement indépendant des solveurs de structures et est bien adapté au couplage fluide-structure partitionné. Enfin, nous proposons un schéma de couplage semi-implicite permettant de mêler efficacement les formulations ALE et FD. Toutes les méthodes considérées sont accompagnées de nombreux tests numériques en 2D et 3D.

Les travaux résumés dans ce mémoire s'articulent, essentiellement, autour des deux thématiques suivantes: la modélisation et la simulation numériques de systèmes couplés (Chapitres 1-3) et les méthodes d'éléments finis stabilisées pour des problèmes transitoires (Chapitre 4). Ces travaux sont essentiellement motivés par l'étude de la stabilité aéroélastique de structures du génie civil et la simulation numérique de l'écoulement du sang et de l'électrophysiologie cardiaque. Dans le cadre de l'interaction fluide-structure, nous couplons les équations de Navier-Stokes en domaine mobile avec l'équation de l'élastodynamique non-linéaire. Nous étudions la stabilité des états d'équilibre du système à partir de l'analyse des solutions harmoniques d'un problème linéaire spécifique. Dans le contexte de la simulation temporelle, nous proposons une méthode de Newton exacte pour la résolution des schémas de couplage implicite. Puis nous nous intéressons à la question suivante: comment éviter le couplage fort sans compromettre la stabilité? Cette question est abordée de deux points de vue différents: via le couplage semi-implicite avec projection et par un traitement faible approprié des conditions d'interface au niveau discret. Nous abordons aussi la simulation numérique des ECG en utilisant un modèle mathématique 3D complet, entièrement basée sur des EDP/EDO. Les principaux ingrédients de ce modèle sont: dynamique phénoménologique au niveau cellulaire, équation bidomaine (dans le cœur) et équation de Laplace généralisée (dans le torse). D'autres aspects essentiels à la modélisation sont élucidés, ce qui nous permet de simuler des ECGs complets réalistes. Quelques schémas de discrétisation en temps pour l'équation bidomaine et le système couplé cœur-torse sont analysés. Enfin, nous généralisons la méthode de pénalisation intérieure conforme au problème d'Oseen et aux équations de Navier-Stokes transitoires. Des estimations d'erreur a priori (uniformes par rapport à la viscosité) sont fournies pour des approximations vitesse/pression du même ordre. Une analyse d'erreur abstraite pour des méthodes de stabilisation symétriques est présentée pour l'équation de Stokes et l'équation de réaction-advection-diffusion transitoires. Dans le cas de Stokes, nous montrons que l'instabilité des petits pas de temps peut être éliminée par un choix judicieux de l'approximation de la vitesse initiale. Pour l'équation de réaction-advection-diffusion, nous contournons le problème de la réduction de la structure creuse de la matrice (due à l'opérateur de stabilisation) par un traitement explicite de la stabilisation.

Les travaux de cette thèse contribuent à la vérification et à la validation d'un modèle diphasique instationnaire, le mo- dèle bifluide de Baer-Nunziato, pour modéliser les phénomènes de transitoires hydrauliques tels que les coups de bélier et les marteaux d'eau, qui peuvent apparaître dans les réseaux de tuyauteries industrielles. Il s'agit d'abord de modéliser les écoulements de transitoires hydrauliques avec le modèle bifluide en représentation eulérienne, puis d'étendre ce modèle en formalisme ALE (Arbitrary Lagrangian Eulerian) pour prendre en compte l'interaction fluide-structure (IFS). Pour mo- déliser les écoulements, des lois de fermetures du modèle bifluide concernant les termes interfaciaux, les termes sources et les lois thermodynamiques ont d'abord été étudiées. Ensuite, le système complet a été simulé avec une méthode à pas fractionnaires qui admet deux étapes, l'une pour la résolution de la partie convective, l'autre pour la prise en compte des termes sources. Pour la partie convective, le schéma de Rusanov a d'abord été vérifié. Des problèmes de stabilité ayant été observés, un nouveau schéma plus stable à pas fractionnaires a été proposé et vérifié. En ce qui concerne les termes sources, quatre schémas de relaxation non-instantanés qui représentent respectivement le retour à l'équilibre de pression, les transferts de quantité de mouvement, de chaleur et de masse, sont appliqués successivement. Ces schémas ont été étendus aux lois thermodynamiques 'Stiffened Gas généralisées' afin de représenter le changement de phase eau-vapeur. Après avoir retrouvé certains phénomènes typiques associés aux écoulements de transitoires hydrauliques le modèle bi- fluide a été confronté aux résultats de l'expérience de Simpson, qui est un cas classique de coup de bélier et à ceux de l'expérience Canon, dédiée à la décompression rapide d'un fluide à haute pression dans une tuyauterie. Par ailleurs, le modèle bifluide a été comparé avec deux modèles homogènes sur ces deux expériences. Enfin, une version ALE du mo- dèle bifluide a été mise en œuvre et vérifiée sur un cas de propagation d'ondes de pression dans une conduite flexible, en écoulement 'quasi-monophasique', ou diphasique. La variation de la célérité des ondes de pression dans le fluide liée au couplage entre le fluide et la structure a été bien retrouvée. La validation a été effectuée sur une expérience qui étudie la réponse d'une tuyauterie remplie d'eau soumise à un pic de pression violent (140 bar). Les simulations sont en bon accord avec les données expérimentales.