This thesis provides a comprehensive investigation of flow-sound-structure coupling in spring-loaded valves subjected to air flow. While they are commonly used in a multitude of applications, these types of valves have been found to experience severe vibrations when interaction is present among the structure, the hydrodynamic field, and the acoustic field for a range of operational valve structural characteristics, flow parameters, and connected piping length. The first part of this investigation was aimed at characterizing experimentally the valve’s dynamic behaviour and the parameters affecting the onset of self-excited instability. The occurrence of instability was mainly driven by the presence of acoustic feedback: the connected length of piping had to be sufficiently long, with a longer pipe correlating to more severe vibrations. In addition, it was found that the valve’s oscillation frequency depends on the modal characteristics of the combined valve piping system, rather than the structural natural frequency alone. Furthermore, an increase in the valve’s spring stiffness caused the vibrations to become more severe. Meanwhile, other parameters such as initial spring preload force and valve plate area only had moderate effects on the stability behaviour of the valve. The second part of the investigation sought to develop a theoretical model that could simulate the valve’s response when subjected to air flow while considering the effects of acoustic feedback and impact on the seat and limiter. Thus, a structural model of the valve was developed based on a single-degree-of-freedom model of the system with impact computed based on a pseudo-force method. The hydrodynamic field relied on a one dimensional unsteady Bernoulli description of the flow. Finally, the acoustic interaction was accounted for using the one-dimensional wave equation resolved using a finite difference scheme. The model has demonstrated remarkable agreement with the experimental results. It has shown an ability to predict the modal characteristics of the system as well as correctly predict the effect of increased stiffness or increased piping length on vibration amplitude. The final part of the investigation consisted in designing countermeasures to mitigate the effects of this self-excited instability mechanism. A concentric Helmholtz-type cavity resonator, an orifice plate, and an anechoic termination are placed at the downstream side of a model valve which were seen to be unstable in the experimental and modelling phases of the investigation. All tested devices were able to eliminate the self excited instability mechanism. The applicability and robustness of each of these methods were discussed. / Thesis / Doctor of Philosophy (PhD)

A Computational Fluid Dynamics (CFD) tool is developed in this thesis to solve complex fluid-structure interaction (FSI) problems. The fluid domain is based on Smoothed Particle Hydro-dynamics (SPH) and the structural domain employs large-deformation Finite Element Method (FEM). Validation tests of SPH and FEM are first performed individually. A loosely-coupled SPH-FEM model is then proposed for solving FSI problems. Validation results of two benchmark FSI problems are illustrated (Antoci et al., 2007; Souto-Iglesias et al., 2008). The first test case is flow in a sloshing tank interacting with an elastic body and the second one is dam-break flow through an elastic gate. The results obtained with the SPH-FEM model show good agreement with published results and suggest that the SPH-FEM model is a viable and effective numerical tool for FSI problems. This research is then applied to simulate a two-dimensional free-stream flow interacting with a deformable, pressurized surface, such as an ACV/SES bow seal. The dynamics of deformable surfaces such as the skirt/seal systems of the ACV/SES utilize the large-deformation FEM model. The fluid part including the air inside the chamber and water are simulated by SPH. A validation case is performed to investigate the application of SPH-FEM model in ACV/SES via comparison with experimental data (Zalek and Doctors, 2010). The thesis provides the theory of the SPH and FEM models incorporated and the derivation of the loosely-coupled SPH-FEM model. The validation results have suggested that this SPH-FEM model can be readily applied to skirt/seal dynamics of ACV/SES interacting with free-surface flow. / Ph. D.

Underwater explosions (UNDEX) and a ship's vulnerability to them are problems of interest in early-stage ship design. A series of events occur sequentially in an UNDEX scenario in both the fluid and structural domains and these events happen over a wide range of time and spatial scales. Because of the complexity of the physics involved, it is a common practice to separate the description of UNDEX into early-time and late-time, and far-field and near-field. The research described in this dissertation is focused on the simulation of near-field and early-time UNDEX. It assembles a hybrid framework of algorithms to provide results while maintaining computational efficiency. These algorithms include Runge-Kutta, Discontinuous Galerkin, Level Set, Direct Ghost Fluid and Embedded Boundary methods. Computational fluid dynamics (CFD) solvers are developed using this framework of algorithms to demonstrate the computational methods and their ability to effectively and efficiently solve UNDEX problems. Contributions, made in the process of satisfying the objective of this research include: the derivation of eigenvectors of flux Jacobians and their application to the implementation of the slope limiter in the fluid discretization; the three-dimensional extension of Direct Ghost Fluid Method and its application to the multi-fluid treatment in UNDEX flows; the enforcement of an improved non-reflecting boundary condition and its application to UNDEX simulations; and an improvement to the projection-based embedded boundary method and its application to fluid-structure interaction simulations of UNDEX problems. / Doctor of Philosophy / Underwater explosions (UNDEX) and a ship's vulnerability to them are problems of interest in early-stage ship design. A series of events occur sequentially in an UNDEX scenario in both the fluid and structural domains and these events happen over a wide range of time and spatial scales. Because of the complexity of the physics involved, it is a common practice to separate the description of UNDEX into early-time and late-time, and far-field and near-field. The research described in this dissertation is focused on the simulation of near-field and early-time UNDEX. It assembles a hybrid framework of algorithms to provide results while maintaining computational efficiency. These algorithms include Runge-Kutta, Discontinuous Galerkin, Level Set, Direct Ghost Fluid and Embedded Boundary methods. Computational fluid dynamics (CFD) solvers are developed using this framework of algorithms to demonstrate these computational methods and their ability to effectively and efficiently solve UNDEX problems.

<p dir="ltr">The CFD-CAA method combines computational fluid dynamics (CFD) and computational aeroacoustics (CAA) techniques to analyze the interaction between fluid flow and the generation and propagation of sound. CFD is primarily concerned with simulating fluid flow patterns, while CAA focuses on predicting noise generation and its propagation in fluids. The CFD-CAA method provides a powerful tool for understanding and predicting the acoustic behavior of turbulent flows. By combining the strengths of CFD and CAA, this approach provides more precise and comprehensive analyses across various fields, thereby contributing to enhanced designs and noise control strategies.</p><p dir="ltr">Within industrial applications, a primary concern is noise source identification. This process enables engineers to locate and quantify the strength of noise sources within a system, facilitating the implementation of more effective strategies during the design process. A novel methodology, computational statistically optimized near-field acoustic holography (C-SONAH), is proposed to virtually identify aeroacoustic sources. Initially, sound pressure is obtained using the CFD-CAA method, followed by the application of the SONAH algorithm to locate acoustic sources and predict the sound field. C-SONAH offers computational advantages over direct CAA methods for simulating sound produced by systems with rotating elements, as CAA analyzes sources on the moving elements, making sound field calculation computationally expensive. The SONAH procedure converts these rotating sources into a series of equivalent stationary planar or cylindrical waves, reducing the number of sources and the time required to compute the sound field from each source. This methodology was demonstrated by characterizing the aerodynamic noise produced by a bladeless fan. The sound pressure level obtained by C-SONAH method was validated by the data predicted by the direct CFD-CAA method. Acoustic maps were reconstructed at different locations and frequencies, revealing that the C-SONAH method can predict noise sources generated by airflow and rotating components within the fan. Thus, it serves as an effective tool for understanding the aeroacoustic noise generation mechanism and guiding the design optimization of similar products.</p><p dir="ltr">The CFD-CAA method is also a powerful tool for design optimization. Computational simulations are typically less expensive and time-consuming than building and maintaining experimental setups, particularly for large or complex projects. Additionally, simulations reduce the need for multiple physical prototypes, which can shorten the development cycle. CFD-CAA simulations provide detailed flow and acoustic field data, including variables that may be difficult or impossible to measure experimentally, such as pressure distributions, velocity fields, and turbulent structures. In this dissertation, aeroacoustic characteristics and flow field information of vortex whistles were investigated using the CFD-CAA method. The simulation results clearly illustrate the swirling motion created in the vortex whistle cylinder and also demonstrate the linear frequency versus flow rate relationship characteristic of the whistle. The design of the vortex whistle was optimized based on the acoustic response and flow resistance by both simulations and experiments. The results suggest that the whistle with a thin inlet exhibits the best performance at high flow rates, while the whistle with a scale of 0.5 is the most sensitive to low flow rates, making it suitable for pediatric applications.</p><p dir="ltr">In CFD-CAA simulations, the time step typically cannot be too small due to limited computational resources. This constraint results in an aliasing error in spectral analysis. Consequently, an anti-aliasing operation prior to sampling is necessary to eliminate such errors from the acoustic source terms. In the present study, an anti-aliasing filter based on the compact finite difference formulation was designed within a time-domain, compact filter scheme. This filter was directly applied to the Navier-Stokes solver prior to sampling for CAA analysis. A cavity flow case was simulated to validate this mitigation strategy. The results indicate that the artificial spectral peak induced by aliasing error is removed without affecting other signature peaks. The anti-aliasing filter was also applied to more complex cases for predicting the acoustic field of a vortex whistle. The acoustic field around the vortex whistle, with both constant and variable inlet flow rates, was simulated, and the aliasing peak was successfully removed. Although the peak magnitudes decreased slightly due to the filter, the signature frequencies remained unchanged. Thus, the simulation with anti-aliasing operation can predict acoustic features without introducing aliasing errors, even if the time step is not sufficiently small, thereby significantly reducing simulation time.</p><p dir="ltr">In engineering applications, once noise sources are identified, the subsequent concern is noise reduction. An effective strategy for noise reduction involves acoustical absorbing materials to minimize noise emissions from components. Traditionally, experiments in engineering applications have focused on surface treatments to explore noise control techniques. However, the CFD-CAA method commonly assumes smooth and purely reflective wall surfaces. Consequently, there is growing interest in incorporating impedance boundary conditions into the CFD-CAA method. Since impedance boundary conditions are defined in the frequency domain, while CFD-CAA simulations operate in the time domain, direct implementation is not feasible. To address this issue, several methods have been proposed to define time-domain impedance boundary conditions in simulations. In the present study, a wall softness model was implemented in the CFD-CAA method and to examine a vortex whistle featuring an acoustically permeable surface. In simulations, an impedance boundary condition representing the properties of melamine foam was defined over the surface of a cylindrical cavity. The simulation results were validated against experimental data obtained from a vortex whistle with melamine foam. The findings revealed that the impedance of the melamine foam contributed to noise reduction at high frequencies. Additionally, at low airflow rates, the impedance boundary condition enhanced the signal-to-noise ratio for the low-frequency peak, which is advantageous in clinical applications.</p>

[pt] O escoamento de cápsulas suspensas em uma fase líquida através de canais e capilares micrométricos representa um problema complexo que ocorre em diferentes aplicações, de glóbulos vermelhos em hemodinâmica até escoamento em meios porosos. Em aplicações de meios porosos, a compreensão da dinâmica na microescala é fundamental para avaliar o comportamento macroscópico do escoamento. Canais e capilares com constrição podem ser usados para modelar uma garganta conectando dois poros adjacentes. O escoamento de uma cápsula suspensa através de tais modelos foi analisado para avaliar as características do escoamento considerando os efeitos inerciais (isto é, número de Reynolds finito), incluindo a máxima diferença de pressão necessária para empurrar uma cápsula através da constrição em função do raio da cápsula, a tensão inicial e o material da membrana, geometria do canal e do capilar, assim como as condições de escoamento. De fato, neste estudo, a resposta da pressão é fundamental para avaliar o efeito de bloqueio da cápsula. As fases líquidas internas e externas foram descritas pelas equações de Navier-Stokes, enquanto que a dinâmica da membrana da cápsula foi modelada por uma estrutura flexível 1-D tipo mola. O problema de interação fluido-estrutura foi resolvido usando o método de elementos finitos acoplado ao método de fronteira imersa. Os resultados mostraram a redução da mobilidade da fase contínua devido à presença da cápsula através da constrição. Tais resultados podem ser usados para projetar microcápsulas para bloquear caminhos preferenciais de fluxo da água no processo de deslocamento de óleo em meios porosos. / [en] The flow of capsules suspended in a liquid phase through small channels and capillaries poses a complex problem presented in different applications, from red blood cells on hemodynamics to flow in porous media. In applications of porous media, the understanding of microscale dynamics is fundamental to assess the macroscopic flow behavior. Constricted channels and capillaries can be used to model a pore throat connecting two adjacent pore bodies. The flow of a suspended capsule through such models was analyzed to evaluate the flow characteristics considering inertial effects (i.e. finite Reynolds numbers), including the maximum pressure difference required to push a capsule through the constriction as a function of capsule radius, initial membrane tension, membrane material, channel and capillary geometries, as well as flow conditions. In fact, in this study, the pressure response is fundamental in order to assess the capsule blocking mechanism. Inner and outer liquid phases were described by the Navier-Stokes equations and capsule membrane dynamics was modeled by a 1-D spring-like flexible structure. The fluid-structure interaction problem was solved using the finite element method coupled with the immersed boundary method. Results showed the mobility reduction of the continuous phase due to the presence of a capsule as it flows through the constriction. Such results can be used to design microcapsules to block preferential water flow paths in oil displacement process in porous media.

L'Interaction Fluide-Structure (IFS) décrit une classe très générale de problème physique, ce qui explique la nécessité de développer une méthode numérique capable de simuler le problème FSI. Pour cette raison, un solveur IFS est développé qui peut traiter un écoulement de fluide incompressible en interaction avec des structures différente: élastique ou rigide. Dans cet aspect, le solveur peut couvrir une large gamme d'applications.La méthode proposée est développée dans le cadre d'une formulation monolithique dans un contexte Eulérien. Cette méthode consiste à considérer un seul maillage et résoudre un seul système d'équations avec des propriétés matérielles différentes. La fonction distance permet de définir la position et l'interface de tous les objets à l'intérieur du domaine et de fournir les propriétés physiques pour chaque sous-domaine. L'adaptation de maillage anisotrope basé sur la variation de la fonction distance est ensuite appliquée pour assurer une capture précise des discontinuités à l'interface fluide-solide.La formulation monolithique est assurée par l'ajout d'un tenseur supplémentaire dans les équations de Navier-Stokes. Ce tenseur provient de la présence de la structure dans le fluide. Le système est résolu en utilisant une méthode élément fini et stabilisé suivant la formulation variationnelle multiéchelle. Cette formulation consiste à décomposer les champs de vitesse et pression en grande et petite échelles. La particularité de l'approche proposée réside dans l'enrichissement du tenseur de l'extra contraint.La première application est la simulation IFS avec un corps rigide. Le corps rigide est décrit en imposant une valeur nul du tenseur des déformations, et le mouvement est obtenu par la résolution du mouvement de corps rigide. Nous évaluons le comportement et la précision de la formulation proposée dans la simulation des exemples 2D et 3D. Les résultats sont comparés avec la littérature et montrent que la méthode développée est stable et précise.La seconde application est la simulation IFS avec un corps élastique. Dans ce cas, une équation supplémentaire est ajoutée au système précédent qui permet de résoudre le champ de déplacement. Et la contrainte de rigidité est remplacée par la loi de comportement du corps élastique. La déformation et le mouvement du corps élastique sont réalisés en résolvant l'équation de convection de la Level-Set. Nous illustrons la flexibilité de la formulation proposée par des exemples 2D. / Numerical simulations of fluid-structure interaction (FSI) are of first interest in numerous industrial problems: aeronautics, heat treatments, aerodynamic, bioengineering... Because of the high complexity of such problems, analytical study is in general not sufficient to understand and solve them. FSI simulations are then nowadays the focus of numerous investigations, and various approaches are proposed to treat them. We propose in this thesis a novel monolithic approach to deal with the interaction between an incompressible fluid flow and rigid/ elastic material. This method consists in considering a single grid and solving one set of equations with different material properties. A distance function enables to define the position and the interface of any objects with complex shapes inside the volume and to provide heterogeneous physical properties for each subdomain. Different anisotropic mesh adaptation algorithms based on the variations of the distance function or on using error estimators are used to ensure an accurate capture of the discontinuities at the fluid-solid interface. The monolithic formulation is insured by adding an extra-stress tensor in the Navier-Stokes equations coming from the presence of the structure in the fluid. The system is then solved using a finite element Variational MultiScale (VMS) method, which consists of decomposition, for both the velocity and the pressure fields, into coarse/resolved scales and fine/unresolved scales. The distinctive feature of the proposed approach resides in the efficient enrichment of the extra constraint. In the first part of the thesis, we use the proposed approach to assess its accuracy and ability to deal with fluid-rigid interaction. The rigid body is prescribed under the constraint of imposing the nullity of the strain tensor, and its movement is achieved by solving the rigid body motion. Several test case, in 2D and 3D with simple and complex geometries are presented. Results are compared with existing ones in the literature showing good stability and accuracy on unstructured and adapted meshes. In the second, we present different routes and an extension of the approach to deal with elastic body. In this case, an additional equation is added to the previous system to solve the displacement field. And the rigidity constraint is replaced with a corresponding behaviour law of the material. The elastic deformation and motion are captured using a convected level-set method. We present several 2D numerical tests, which is considered as classical benchmarks in the literature, and discuss their results.

Une fistule artérioveineuse (FAV) est un accès vasculaire permanent créé par voie chirurgicale en connectant une veine et une artère chez le patient en hémodialyse. Cet accès vasculaire permet de mettre en place une circulation extracorporelle partielle afin de remplacer les fonctions exocrines des reins. En France, environ 36000 patients sont atteint d’insuffisance rénale chronique en phase terminale, stade de la maladie le plus grave qui nécessite la mise en place d’un traitement de suppléance des reins : l’hémodialyse. La création et présence de la FAV modifient significativement l’hémodynamique dans les vaisseaux sanguins, au niveau local et systémique ainsi qu’à court et à plus long terme. Ces modifications de l’hémodynamiques peuvent induire différents pathologies vasculaires, comme la formation d’anévrysmes et de sténoses. L’objectif de cette étude est de mieux comprendre le comportement mécanique et l’hémodynamique dans les vaisseaux de la FAV. Nous avons étudié numériquement les interactions fluide-structure (IFS) au sein d’une FAV patient-spécifique, dont la géométrie a été reconstruite à partir d’images médicales acquises lors d’un précédent doctorat. Cette FAV a été créée chez le patient en connectant la veine céphalique du patient à l’artère radiale et présente une sténose artérielle réduisant de 80% la lumière du vaisseau. Nous avons imposé le profil de vitesse mesuré sur le patient comme conditions aux limites en entrée et un modèle de Windkessel au niveau des sorties artérielle et veineuse. Nous avons considéré des propriétés mécaniques différentes pour l’artère et la veine et pris en compte le comportement non-Newtonien du sang. Les simulations IFS permettent de calculer l’évolution temporelle des contraintes hémodynamiques et des contraintes internes à la paroi des vaisseaux. Nous nous sommes demandées aussi si des simulations non couplées des équations fluides et solides permettaient d’obtenir des résultats suffisamment précis tout en réduisant significativement le temps de calcul, afin d’envisager son utilisation par les chirurgiens. Dans la deuxième partie de l’étude, nous nous sommes intéressés à l’effet de la présence d’une sténose artérielle sur l’hémodynamique et en particulier à ses traitements endovasculaires. Nous avons dans un premier temps simulé numériquement le traitement de la sténose par angioplastie. En clinique, les sténoses résiduelles après angioplastie sont considérées comme acceptables si elles obstruent moins de 30% de la lumière du vaisseau. Nous avons donc gonflé le ballonnet pour angioplastie avec différentes pressions de manière à obtenir des degrés de sténoses résiduelles compris entre 0 et 30%. Une autre possibilité pour traiter la sténose est de placer un stent après l’angioplastie. Nous avons donc dans un deuxième temps simulé ce traitement numériquement et résolu le problème d’IFS dans la fistule après la pose du stent. Dans ces simulations, la présence du stent a été prise en compte en imposant les propriétés mécaniques équivalentes du vaisseau après la pose du stent à une portion de l’artère. Dans la dernière partie de l’étude nous avons mis en place un dispositif de mesure par PIV (Particle Image Velocimetry). Un moule rigide et transparent de la géométrie a été obtenu par prototypage rapide. Les résultats expérimentaux ont été validés par comparaison avec les résultats des simulations numériques. / An arteriovenous fistula (AVF) is a permanent vascular access created surgically connecting a vein onto an artery. It enables to circulate blood extra-corporeally in order to clean it from metabolic waste products and excess of water for patients with end-stage renal disease undergoing hemodialysis. The hemodynamics results to be significantly altered within the arteriovenous fistula compared to the physiological situation. Several studies have been carried out in order to better understand the consequences of AVF creation, maturation and frequent use, but many clinical questions still lie unanswered. The aim of the present study is to better understand the hemodynamics within the AVF, when the compliance of the vascularwall is taken into account. We also propose to quantify the effect of a stenosis at the afferent artery, the incidence of which has been underestimated for many years. The fluid-structure interactions (FSI) within a patient-specific radio-cephalic arteriovenous fistula are investigated numerically. The considered AVF presents an 80% stenosis at the afferent artery. The patient-specific velocity profile is imposed at the boundary inlet, and a Windkessel model is set at the arterial and venous outlets. The mechanical properties of the vein and the artery are differentiated. The non-Newtonian blood behavior has been taken into account. The FSI simulation advantageously provides the time-evolution of both the hemodynamic and structural stresses, and guarantees the equilibrium of the solution at the interface between the fluid and solid domains. The FSI results show the presence of large zones of blood flow recirculation within the cephalic vein, which might promote neointima formation. Large internal stresses are also observed at the venous wall, which may lead to wall remodeling. The fully-coupled FSI simulation results to be costly in computational time, which can so far limit its clinical use. We have investigated whether uncoupled fluid and structure simulations can provide accurate results and significantly reduce the computational time. The uncoupled simulations have the advantage to run 5 times faster than the fully-coupled FSI. We show that an uncoupled fluid simulation provides informative qualitative maps of the hemodynamic conditions in the AVF. Quantitatively, the maximum error on the hemodynamic parameters is 20%. The uncoupled structural simulation with non-uniform wall properties along the vasculature provides the accurate distribution of internal wall stresses, but only at one instant of time within the cardiac cycle. Although partially inaccurate or incomplete, the results of the uncoupled simulations could still be informative enough to guide clinicians in their decision-making. In the second part of the study we have investigated the effects of the arterial stenosis on the hemodynamics, and simulated its treatment by balloon-angioplasty. Clinically, balloon-angioplasty rarely corrects the stenosis fully and a degree of stenosis remains after treatment. Residual degrees of stenosis below 30% are considered as successful. We have inflated the balloon with different pressures to simulate residual stenoses ranging from 0 to 30%. The arterial stenosis has little impact on the blood flow distribution: the venous flow rate remains unchanged before and after the treatment and thus permits hemodialysis. But an increase in the pressure difference across the stenosis is observed, which could cause the heart work load to increase. To guarantee a pressure drop below 5 mmHg, which is considered as the threshold stenosis pressure difference clinically, we find that the residual stenosis degree must be 20% maximum.

Ce travail de thèse s’appuie sur une étude expérimentale en tube à choc, plus précisément, c’est une approche expérimentale de l'étude de l'interaction fluide-structure. Considérons un solide indéformable auquel on laisse un degré de liberté en rotation autour d'un axe. Cette structure ferme un espace clos. Si le contenu de l'espace clos subit le passage d'une onde de choc, ce solide va être mis mouvement et tourner autour de son axe. Concrètement, l'onde de choc va augmenter les caractéristiques physiques, en particulier sa pression, du fluide en contact avec la face impactée de cette porte. La face opposée de la porte ne subissant pas ou que très peu l'influence de l'onde de choc, une seule de ses faces est soumise à la surpression. Au moment de l'impact, le déséquilibre ainsi créé impose une action mécanique sur la porte qui va la faire accélérer et tourner autour de son axe de rotation. Jusqu'à ce stade tout est relativement simple. La difficulté intervient à l'instant où la porte commence à s'ouvrir, car les frontières du volume dans lequel le fluide évolue sont modifiées. Des fuites apparaissent et le gaz qui était maintenu dans un volume clos peut maintenant s'écouler vers un milieu libre. Une communication entre les gaz agissant de chaque coté de la porte est créée modifiant leurs propriétés et par conséquent la pression agissant sur chaque côté de la porte. Les actions mécaniques qui s'appliquent sur la porte ne sont plus les mêmes, et par conséquent l'accélération que la porte subit aussi. Au fur et à mesure que la porte change de position, le problème fluide continue d'être modifié et change en retour son action sur la porte. Cette interaction perdure soit jusqu'à ce que les limites du problème cessent d'être modifiées, la porte ne peut plus bouger, ou bien lorsque les actions mécaniques agissant sur la porte s'équilibrent, les fluides de chaque côté de la porte étant dans le même état physique. Le travail présenté ici est une étude des paramètres du fluide ou du solide en mouvement qui sont les acteurs de la loi comportementale gérant ce système complexe. Pour ce faire, nous avons réalisé une maquette expérimentale mettant en action la physique que nous venons de décrire et nous l'avons adaptée à un tube à choc. En éprouvant de nombreuses configurations expérimentales, nous avons pu déterminer comment l'écoulement interne d'un tube à choc évolue lorsqu'il est plus ou moins ouvert à son extrémité. Comment une porte fermée réagit-elle à l'impact d'une onde de choc et quelles en sont les conséquences sur l'évolution des fluides mis en jeu? Quelles sont les conséquences d'une position différente de la porte au moment de l'impact avec l'onde de choc? Ou encore, quel rôle joue l'intensité de l'onde de choc incidente ou l'inertie de la porte sur toute cette dynamique? / This thesis is based on an experimental study carried out in shock tube; in particular, this is an experimental approach to the study of fluid-structure interaction. Consider a rigid body which is allowed to rotate only around an axis and which closes a confined space. If a shock wave crosses the content of the confined space, the body will accelerate and rotate around its axis. Specifically, the shock wave will increase the physical characteristics, especially its pressure, of the fluid acting on the impacted face of the door. The opposite side of the door is not influenced by the incident shock wave, only one of its faces is subjected to overpressure. Following the first impact, the resulting imbalance imposes a mechanical action on the door that will increase its speed and make it turn around its rotation axis. The difficulty comes when the door begins to open: the volume boundaries in which the fluid is contained are modified. Leaks occur and the gas kept in this closed volume can now flow to the atmosphere. Communication between the gas acting on each side of the door is created modifying their properties and consequently the pressure acting on each side of the door.The mechanical actions that apply to the door are no more the same with time, and therefore the acceleration of the door is changing. As the door moves, the fluid problem continues to be changed and in turn it changes its action on the door. This interaction process continues until either the limits of the problem ceases to be changed, the door cannot move, or when the mechanical actions acting on the door are in equilibrium, fluids on each side of the door are in the same physical state. The presented work is a study of the parameters of the fluid or the solid motion which are main actors in the behavioral law managing this complex system. In this aim, we designed an experimental device involving the physics that we have described and we have adapted it to a shock tube. Testing many experimental configurations, we could determine how the internal flow of a shock tube evolves when the end of this shock tube is more or less open.How a closed door reacts to the impact of a shock wave and what are the implications for the evolution of the involved fluids? What are the consequences of a different position of the door at the instant of the impact with the incident shock wave? What role plays the intensity of the incident shock wave or the inertia of the door on this dynamic?

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%

Le sujet de la thèse porte sur l'étude (existence, unicité, régularité) et le contrôle de problèmes fluide-structure possédant des conditions limites sur la pression. Le système étudié couple une partie fluide, décrite par les équations de Navier-Stokes incompressibles dans un domaine 2D et une partie structure, décrite par une équation 1D de poutre amortie située sur une partie du bord du domaine fluide. Dans le Chapitre 2, on étudie l'existence de solutions fortes pour ce modèle. Nous démontrons des résultats de régularité optimale pour le système de Stokes avec conditions de bord mixtes sur un domaine non régulier. Ces résultats sont ensuite utilisés pour prouver l'existence et l'unicité de solutions fortes, locales en temps, pour le système fluide-structure sans hypothèse de petitesse sur les données initiales. Le Chapitre 3 réutilise l'analyse précédente dans le cadre de solutions périodiques en temps. Nous développons un critère d'existence de solutions périodiques pour un problème parabolique abstrait. Ce critère est ensuite appliqué au système fluide-structure et nous obtenons l'existence de solutions strictes, périodiques et régulières en temps, pour des termes sources périodiques suffisamment petits. Le quatrième volet de la thèse porte sur la stabilisation du système fluide-structure au voisinage d'une solution périodique. Le système linéarisé sous-jacent est décrit à l'aide d'un opérateur A(t) dont le domaine dépend du temps. Nous démontrons l'existence d'un opérateur parabolique d'évolution pour ce système linéaire. Cet opérateur est ensuite utilisé, dans le cadre de la théorie de Floquet, pour étudier le comportement asymptotique du système. Nous adaptons la théorie existante pour des opérateurs à domaine constant au cas de domaine non constant. Nous obtenons la stabilisation exponentielle du système linéaire à l'aide d'un contrôle sur la frontière du domaine fluide. / In this thesis we study the well-posedness (existence, uniqueness, regularity) and the control of fluid-structure system with boundary conditions involving the pressure. The fluid part of the system is described by the incompressible Navier- Stokes equations in a 2D rectangular type domain coupled with a 1D damped beam equation localised on a boundary part of the fluid domain. In Chapter 2 we investigate the existence of strong solutions for this model. We prove optimal regularity results for the Stokes system with mixed boundary conditions in non-regular domains. These results are then used to obtain the local-in-time existence and uniqueness of strong solutions for the fluid-structure system without smallness assumption on the initial data. Chapter 3 uses the previous analysis in the framework of periodic (in time) solutions. We develop a criteria for the existence of periodic solutions for an abstract parabolic system. This criteria is then used on the fluid- structure system to prove the existence of a periodic and regular in time strict solution, provided that the periodic source terms are small enough. In Chapter 4 we study the stabilisation of the fluid-structure system in a neighbourhood of a periodic solution. The underlying linear system involves an operator A(t) with a domain which depends on time. We prove the existence of a parabolic evolution operator for this linear system. This operator is then used to apply the Floquet theory and to describe the asymptotic behaviour of the system. We adapt the known results for an operator with constant domain to the case of operators with non constant domain. We obtain the exponential stabilisation of the linear system with control acting on a part of the boundary of the fluid domain.