L'imagerie par Résonance Magnétique permet de mesurer l'écoulement sanguin. Au niveau cardiovasculaire, elle permet d'acquérir non seulement des images anatomiques du cœur et des gros vaisseaux mais aussi des images fonctionnelles de vitesse par contraste de phase. Cette technique offre des perspectives dans l'étude de la dynamique des fluides et dans la caractérisation des artères, en particulier pour les grosses artères systémiques comme l'aorte dont le rôle est primordial dans la circulation sanguine. Par ailleurs, l'un des paramètres qui entrent en jeu dans la détermination de la fonction cardiaque et du comportement vasculaire est la pression artérielle. La méthode de référence de la mesure de pression dans l'aorte étant le cathétérisme, plusieurs méthodes combinant la modélisation à l'imagerie ont été proposées afin d'estimer un gradient de pression de façon non invasive. Ce travail de thèse propose de mesurer la pression dans un segment d'aorte grâce à un modèle 1D simplifié et en utilisant les données mesurées par IRM et un modèle 0D représentant le réseau vasculaire périphérique comme conditions aux limites. Aussi, afin d'adapter le modèle à l'aorte du patient, une loi de pression exprimant une relation entre la section aortique à la pression et basée sur la compliance a été utilisée. Cette dernière, liée à la vitesse d'onde de pouls (VOP), a été mesurée en IRM sur les ondes de vitesse.Par ailleurs, les séquences de codage de vitesse et d'accélération sont longues et ponctuées d'artéfacts dus au mouvement du patient. Une apnée est requise afin de limiter le mouvement respiratoire. Cependant, la durée de l'apnée atteint 25 à 30 secondes pour de telles séquences, ce qui est souvent impossible à tenir pour les malades. Une technique d'optimisation de séquences dynamiques par réduction du champ de vue est proposée et étudiée. La technique décrit un dépliement des régions repliées par différence complexe de deux images, l'une codée et l'autre non codée en vitesse. Cette méthode réalise une réduction de plus de 25% de la durée d'apnée. / Magnetic Resonance Imaging (MRI) is used to measure blood flow. It allows assessing not only dynamic images of the heart and the large arteries, but also functional velocity images by means of Phase Contrast. This promising technique is important for studying fluid dynamics and characterizing the arteries, especially the large systemic arteries that play a prominent role in the blood circulation. One of the parameters used for determining the cardiac function and the vascular behavior is the arterial pressure. The reference technique for measuring the aortic pressure is catheterism, but several methods combining imaging and mathematical modeling have been proposed in order to non-invasively estimate a pressure gradient. This work proposes to measure pressure in an aortic segment through a simplified 1D model using MRI measured flow and 0D model representing the peripheral vascular system as boundary conditions. To adapt the model to the aorta of a patient, a pressure law was used forming a relation between the aortic section area and pressure, based on compliance, which is linked to pulse wave velocity (PWV) estimated on MRI measured flow waves.Scan duration was optimized, as it is often a limitation during image acquisition. Velocity and acceleration sequences require a long time and may cause artifacts. Hence, they are acquired during apnea to avoid respiratory motion. However, for such acquisitions, a subject would have to hold their breath for more than 25 seconds which can pose difficulties for some patients. A technique that allows dynamic acquisition time optimization through field of view reduction was proposed and studied. The technique unfolds fold-over regions by complex difference of two images, one of which is motion encoded and the other acquired without an encoding gradient. By implementing this method, we decrease the acquisition time by more than 25%

Dans le cadre de la lutte en mer, la détection acoustique des structures immergées s’effectue généralement sur des signaux stationnaires. Une nouvelle génération de sonars permet de détecter sur des signaux transitoires. Ceci implique de compléter le processus de conception des projets industriels qui ne tient compte d’exigences qu'en matière de bruits rayonnés en régime stationnaire. Il est donc nécessaire de comprendre les mécanismes de transfert des sources de bruit transitoires sur les structures immergées. Cette thèse s’inscrit dans ce cadre et consiste à étudier les mécanismes de transfert vibratoire et de rayonnement acoustique qui peuvent intervenir sur ces structures lorsque l’excitation est transitoire. L’analyse porte sur différents éléments de la chaine de transfert : le rayonnement dans l’eau du bordé, la diffraction des ondes par les raidisseurs, et le comportement résonnant des structures internes supportant les matériels. Le premier chapitre présente une analyse bibliographique autour de l’étude des phénomènes vibroacoustiques transitoires des structures immergées, de l’influence d’un fluide lourd sur le comportement vibroacoustique des plaques, et des méthodes de calcul vibroacoustiques en régime transitoire. Dans le second chapitre nous étudions la réponse transitoire d’une plaque infinie immergée soumise à une force impulsionnelle ponctuelle. La méthode de calcul s’appuie sur les calculs spectraux fréquences-nombre d’onde. Les réponses temporelles sont obtenues par transformées de Fourier inverses. L’analyse des spectres et des réponses temporelles de l’accélération vibratoire de la plaque et de la pression rayonnée, met en évidence l’influence de la présence du fluide. La prise en compte de l’inertie rotationnelle et du cisaillement à travers le modèle de plaque de Mindlin-Timoshenko est également étudiée. Ces résultats sont confrontés à une expérimentation présentée dans le troisième chapitre. La structure étudiée est une plaque rectangulaire posée horizontalement à la surface d’une cuve remplie d’eau. Deux types de sources transitoires sont utilisés : marteau de choc, lâché d’une bille. La comparaison des résultats numériques et expérimentaux montre que l’on retrouve certains phénomènes évoqués précédemment. L’effet des raidisseurs sur le rayonnement acoustique fait l’objet du quatrième chapitre. Une plaque raidie périodiquement dans une direction est considéré. L’influence des ondes de Bloch-Floquet sur la réponse temporelle est étudiée. Les résultats sont comparés à des mesures effectuées sur une barge d’essais. Dans le cinquième chapitre, l’effet des structures internes est étudié à partir d’un modèle de plaque couplé à un système résonnant constitué d’un assemblage poutre-plaque. La méthode des inertances est utilisée pour obtenir les forces de couplage entre les différents éléments. Les signaux temporels sont étudiés en fonction de l’importance de la rupture d’inertance entre la plaque et l’assemblage. / In the sea, the acoustic detection of other battle engines is done by detecting mostly stationary signals. However, new types of detection systems are being developed, and are able to detect and to analyze transient signals. Therefore, the industrial conception process needs to be improved, so that the underwater vehicles transient noises can be taken in account. In order to do so, the mechanism of vibroacoustic transfer of transient sources of submerged structures has to be understood. The object of this thesis is then the study of the mechanism of vibration transfer and acoustic radiation of those structures when they are excited by transient sources. The shell radiation in the water, the wave diffraction by circumferential stiffeners and the resonant behavior of internal substructures are analyzed. The first chapter presents the bibliographical study of three themes: the study of transient phenomenon of submerged structures, the influence of heavy fluid coupling on vibroacoustic behavior of plates, and the different calculation methods in transient vibroacoustics. In the second chapter, we study the transient response of a submerged infinite plate excited by an impulsively point force. First, the calculations are done in the wavenumber-frequency domain. Then the spatio-temporal responses are obtained using inverse Fourier transforms. The discretization of wavenumber and frequency domains and the damping model are studied. The analysis of frequency and time responses of the plate vibration and the radiated pressure enable us to observe the influence of heavy fluid coupling. Besides the Mindlin-Timoshenko plate model is also used and the effect of rotation inertia and shear stress are studied. In the third chapter, these numerical results are confronted to experimental data, obtained experimentally. The studied structure is a rectangular plate lying on the surface of a water tank. Two different excitations are used: an impact hammer and the free fall of a steel ball. The study of the correlation between numerical and experimental results showed that some phenomena are observed in both cases. The influence of stiffeners on the acoustic radiation is the theme of the fourth chapter. An infinite plate which is periodically stiffened through one direction is considered. The effect of Bloch-Floquet waves on time response is studied. Numerical results are compared to measurements data obtained on an industrial submerged structure. In the fifth chapter, the effect of internal structures is analyzed by modelling an infinite plate coupled to a resonant system made of a beam and a rectangular finite plate. The inertance coupling method is used to obtain the coupling forces between the different substructures. Influence of inertance difference between the substructures is illustrated by the time signals.

Le contrôle d'écoulements turbulents est un enjeu majeur en aérodynamique. Cependant, la présence d'un grand nombre de degrés de libertés et d'une dynamique complexe rend délicat la modélisation dynamique de ces écoulements qui est pourtant nécessaire à la conception d'un contrôle efficace. Au cours de cette thèse, différentes directions ont été suivies afin de développer des modèles réduits dans des configurations réalistes d'écoulements et d'utiliser ces modèles pour le contrôle.Premièrement, la décomposition en modes dynamiques (DMD), et certaines de ses variantes, ont été exploitées en tant que base réduite afin d'extraire au mieux le comportement dynamique de l'écoulement. Par la suite, nous nous sommes intéressés à l'assimilation de données 4D-Var qui permet de combiner des informations inhomogènes provenant d'un modèle dynamique, d'observations et de connaissances a priori du système. Nous avons ainsi élaboré des modèles réduits POD et DMD d'un écoulement turbulent autour d'un cylindre à partir de données expérimentales PIV. Finalement, nous avons considéré le contrôle d'écoulement dans un contexte d'interaction fluide/structure. Après avoir montré que les mouvements de solides immergés dans le fluide pouvaient être représentés comme une contrainte supplémentaire dans le modèle réduit, nous avons stabilisé un écoulement de sillage de cylindre par oscillation verticale. / Control of turbulent flows is still today a challenge in aerodynamics. Indeed, the presence of a high number of active degrees of freedom and of a complex dynamics leads to the need of strong modelling efforts for an efficient control design. During this PhD, various directions have been followed in order to develop reduced-order models of flows in realistic situations and to use it for control. First, dynamic mode decomposition (DMD), and some of its variants, have been exploited as reduced basis for extracting at best the dynamical behaviour of the flow. Thereafter, we were interested in 4D-variational data assimilation which combines inhomogeneous informations coming from a dynamical model, observations and an a priori knowledge of the system. POD and DMD reduced-order models of a turbulent cylinder wake flow have been successfully derived using data assimilation of PIV measurements. Finally, we considered flow control in a fluid-structure interaction context. After showing that the immersed body motion can be represented as an additional constraint in the reduced-order model, we stabilized a cylinder wake flow by vertical oscillations.

Dans cette thèse, nous étudions des modèles mathématiques concernant certains problèmes d'écoulement de fluide à densité variable. Le premier chapitre résume l'ensemble de la thèse et se concentre sur les résultats obtenus, la nouveauté et la comparaison avec la littérature existante. Dans le deuxième chapitre, nous étudions la stabilisation locale des équations non homogènes de Navier-Stokes dans un canal 2d autour du flot de Poiseuille. Nous concevons un contrôle feedback de la vitesse qui agit sur l'entrée du domaine de sorte que la vitesse et la densité du fluide soient stabilisées autour du flot de Poiseuille, à condition que la densité initiale soit donnée par une constante additionnée d'une perturbation dont le support se situe loin du bord latéral du canal. Dans le troisième chapitre, nous étudions un système couplant les équations de Navier-Stokes compressibles à une structure élastique située à la frontière du domaine fluide. Nous prouvons l'existence locale de solutions solides pour ce système couplé. Dans le quatrième chapitre, notre objectif est d'étudier la nulle- contrôlabilité d'un problemè d'interaction fluide-structure linéarisé dans un canal bi dimensional. L'écoulement du fluide est ici modélisé par les équations de Navier-Stokes compressibles. En ce qui concerne la structure, nous considérons une poutre de type Euler-Bernoulli amortie située sur une partie du bord. Dans ce chapitre, nous établissons une inégalité d'observabilité pour le problème considéré d'interaction fluid-structure linéarisé qui constitue le premier pas vers la preuve de la nulle contrôlabilité du système. / In this thesis we study mathematical models concerning some fluid flow problems with variable density. The first chapter is a summary of the entire thesis and focuses on the results obtained, novelty and comparison with the existing literature. In the second chapter we study the local stabilization of the non-homogeneous Navier-Stokes equations in a 2d channel around Poiseuille flow. We design a feedback control of the velocity which acts on the inflow boundary of the domain such that both the fluid velocity and density are stabilized around Poiseuille flow provided the initial density is given by a constant added with a perturbation, such that the perturbation is supported away from the lateral boundary of the channel. In the third chapter we prove the local in time existence of strong solutions for a system coupling the compressible Navier-Stokes equations with an elastic structure located at the boundary of the fluid domain. In the fourth chapter our objective is to study the null controllability of a linearized compressible fluid structure interaction problem in a 2d channel where the structure is elastic and located at the fluid boundary. In this chapter we establish an observability inequality for the linearized fluid structure interaction problem under consideration which is the first step towards the direction of proving the null controllability of the system.

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. / Our aim was to develop some partitioned procedures and monolithic semi-implicit algorithms for solving the interaction between a fluid governed by Navier Stokes equations and a structure. In the first chapter, we propose a partitioned procedures semi-implicit algorithm for solving fluid-structure interaction problems, with a structure governed either by linear elasticity equations or by the non-linear Saint-Venant Kirchhoff model. In the second chapter, we present a partitioned procedures semi-implicit algorithm for solving fluid-structure interaction problem with a linear model for the structure and we prove an unconditional stability result of the algorithm. In the above algorithms, an optimization problem must be solved in order to get the continuity of the velocity as well as the continuity of the stress at the interface. During the iterations of BFGS for solving the optimization problem, the fluid mesh does not move and the fluid matrix is only factorized once, which reduces the computational effort. In the fast chapter, we present a monolithic semi-implicit algorithm for solving fluid-structure interaction problem with linear model for the structure. The algorithm uses one global mesh for the fluid-structure domain. The continuity of velocity at the interface is automatically satisfied and the continuity of stress does not appear explicitly in the global weak form due to the action and reaction principle. At each time step, we have to solve a monolithic system of unknowns velocity and pressure defined on the global fluid-structure domain. When the monolithic approach is used the CPU time is reduced compared to a particular partitioned procedures strategy.

L'exploration pétrolière en eaux profondes crée des nouveaux défis technologiques. Certains de ces problèmes sont liés à des très longues structures servant à relier la plate-forme à la tête de puits. L'augmentation de la profondeur implique des structures plus longues et, par conséquent, plus souple. Deux problèmes liés à ce type de structure sont étudiés dans ce rapport : le positionnement de la structure pour sa connexion à la tête de puits (opération de rentrée), et les vibrations induites par vortex. La première partie de ce rapport est consacrée au contrôle de l'opération de rentrée. Deux stratégies différentes de planification du mouvement sont proposées. La première est basée sur un modèle obtenu par approche modale. L'autre stratégie approxime le comportement de la structure par l'équation d'un câble de Bernoulli avec amortissement. Les solutions approximées de cette EDP sont directement utilisées dans la conception du contrôleur. Deux types de contrôle en boucle fermée sont proposées: un système de suivi de trajectoire en utilisant une fonction de Lyapunov, et un autre qui utilise l'inverse du système. La deuxième partie de ce rapport est consacrée à l'étude des vibrations induites par vortex (VIV) et à leur contrôle. Elle présente une première stratégie de contrôle afin de réduire les VIV. Cette stratégie est basée sur une analyse modale des équations du système. La loi de commande génère un déplacement en opposition de phase par rapport au VIV, l'atténuation résultant d'un effet d'antirésonance. Les résultats de simulation sont confirmés par des expériences réalisées sur une maquette en modèle réduit.

Ce travail de recherche traite du couplage entre un liquide incompressible, irrotationnel et son contenant : une structure élastique. Cette interaction fluide-structure est traitée dans le cadre des petites déformations autour d'un état d'équilibre.Dans un premier temps, on présente une méthode d'introduction des sources dissipatives visqueuses dans le liquide à partir des équations du système couplé conservatif en s'appuyant sur une approche de type fluide potentiel généralement utilisée pour traiter les problèmes de couplage fluide-structure linéarisés non amortis. Un modèle d'amortissement diagonal est alors choisi pour le liquide et les effets dissipatifs de celui-ci sont pris en compte en calculant les coefficients d'amortissement modaux. Seuls les effets dissipatifs liées à la viscosité du liquide sont alors pris en compte. Le système couplé dissipatif obtenu possède une matrice d'amortissement non symétrique. Une résolution de ce système à amortissement non classique est alors présentée et les expressions des réponses fréquentielle et temporelle linéarisées sont données pour différents types d'excitations.Dans un deuxième temps, le liquide est supposé non visqueux et les forces de tension surfacique sont prises en compte. Cette configuration concerne principalement les satellites où le système couplé est en situation de microgravité. Une formulation du problème conservatif permettant de prendre en compte l'incompressibilité du fluide, la condition de continuité à l'interface fluide structure, les effets de capillarité du fluide ainsi que les effets éventuels de précontraintes statiques est alors établie. On se propose pour cela d'utiliser une méthode énergétique via le Principe de Moindre Action. La démarche est alors décomposée en deux étapes : une étude statique afin de déterminer la position de référence, puis une étude dynamique linéarisée autour de cette position d'équilibre. Cette formulation forme notamment une base pour l'introduction des sources dissipatives liées aux effets de capillarité via la méthode précédemment introduite.

Mémoire numérisé par la Division de la gestion de documents et des archives de l'Université de Montréal

Les travaux présentés dans ce mémoire s'articulent essentiellement autour de deux thèmes : l'étude mathématique et numérique des phénomènes d'interaction fluide-structure et la modélisation de l'appareil respiratoire humain. C'est au cours de ma thèse que j'ai commencé à travailler sur les phénomènes d'interaction fluide-structure. Ces travaux se sont poursuivis par la suite au CEREMADE (Université Paris Dauphine) puis au sein du projet REO à l'INRIA. Je me suis essentiellement intéressée au cas des fluides newtoniens, visqueux, incompressibles satisfaisant aux équations de Navier-Stokes. La structure est, quant à elle, mobile : rigide ou déformable. On se place dans le cas ou les déplacements de la structure ont une amplitude telle que l'on ne peut pas les supposer infinitésimaux, ce qui va faire apparaitre des non linéarités géométriques. J'ai étudie ces problèmes aussi bien d'un point de vue théorique (existence de solutions faibles ou fortes pour des systèmes stationnaires ou instationnaires) que d'un point de vue numérique (étude de la stabilité et de la convergence des schémas, conception d'algorithmes performants). Les applications ou de tels phénomènes apparaissent sont nombreuses : tant en aérodynamique qu'en biomécanique. Nous pourrions citer beaucoup d'exemples : écoulement autour d'un bateau, écoulement sanguin dans les artères, pour l'hydrodynamique (fluide en phase liquide), écoulement autour d'ailes d'avion, étude de l'influence des vents sur le tablier d'un pont, pour l'aéroélasticité (fluide en phase gazeuse). En particulier, ils apparaissent dans l'étude des écoulements biologiques : écoulement sanguin dans les artères, écoulement de l'air dans l'appareil respiratoire. Ces dernières années je me suis particulièrement intéressée à ce type d'applications et plus spécifiquement à la modélisation mathématique et numérique de l'appareil respiratoire, tout d'abord dans le cadre d'une ACI nouvelle interface des mathématiques le poumon vous dis-je, puis au sein du projet REO. Je m'intéresse à la modélisation de l'écoulement de l'air dans les voies aériennes, la modélisation du transport et dépôt de particules ainsi qu'à la modélisation des tissus élastiques pulmonaire. Ces applications nécessitent le développement de modèles mathématiques nouveaux et de méthodes numériques spécifiques et adaptées.

Dans le cadre de ce travail, une technique non-intrusive est proposée pour coupler la méthode Smoothed Particle Hydrodynamics (SPH) à la méthode des Eléments Finis afin de résoudre numériquement des problèmes dynamiques et non-linéaires d’interaction fluide-structure en permettant l’utilisation des pas de temps différents dans les deux domaines de calcul (fluide et solide). Ces développements sont motivés par le besoin de simuler numériquement des phénomènes rapides et très non-linéaires qui prennent en compte des impacts en se servant des intégrateurs temporels explicites dans chaque sous-domaine de calcul (Newmark explicite pour le solide et Runge-Kutta 2 pour le fluide). De ce fait, le pas de temps de stabilité est limité par des caractéristiques intrinsèques au modèle numérique du phénomène étudié et en conséquence, il devient important de pouvoir intégrer chaque sous-domaine numérique avec un pas de temps proche de son pas de temps de stabilité. Pour permettre d’utiliser un pas de temps proche du pas de temps de stabilité pour chaque sous-domaine, des méthodes de décomposition de domaines dual-Schur sont implémentées et validées pour des cas en 1-D, 2-D, et 3-D. Des simulations numériques d’impacts de cailloux sur des aubes des turbines hydrauliques sont aussi effectue´es afin de prédire le dommage que cet évènement peut engendrer. / A method to couple smoothed particle hydrodynamics and finite elements methods for nonlinear transient fluid–structure interaction simulations by adopting different time-steps depending on the fluid or solid sub-domains is proposed. These developments were motivated by the need to simulate highly non-linear and sudden phenomena that take into acount solid impacts and hence require the use of explicit time integrators on both sub-domains (explicit Newmark for the solid and Runge–Kutta 2 for the fluid). However, due to critical time-step required for the stability of the explicit time integrators in, it becomes important to be able to integrate each sub-domain with a different time-step while respecting the features that a previously developed mono time-step coupling algorithm offered. For this matter, a dual-Schur decomposition method originally proposed for structural dynamics was considered, allowing to couple time integrators of the Newmark family with different time-steps with the use of Lagrange multipliers.