Contrairement aux liquides ordinaires, les fluides complexes comme le sang exhibent des comportements étranges qui dépendent essentiellement des structures sous-jacentes qui les composent. La simulation des écoulements sanguins continue de poser un formidable défi pour les modélisations théoriques et numériques dont l'intérêt est de développer des méthodes et des outils de simulation pour la communauté médicale. Nous proposons dans cette thèse une contribution à ce projet qui sera majoritairement centré sur les aspects numériques et informatiques. Nous nous sommes particulièrement intéressés à l'interaction entre le sang et la paroi vasculaire, qui joue un rôle important dans les grandes artères comme l'aorte. Nous nous sommes aussi investis dans la simulation du transport des cellules sanguines dans le sang. Pour la résolution des équations aux dérivées partielles décrivant nos modèles d'hémodynamique, nous avons choisi d'utiliser des méthodes numériques dont la précision pourra être accrue de manière arbitraire. Dans ce but, les principaux ingrédients qui ont été mis en oeuvre sont (i) la méthode des éléments finis basée sur des approximations de Galerkin d'ordre arbitraire en espace et géométrie, (ii) la méthode ALE pour la prise en compte de la mobilité des domaines pour des déplacements d'ordre arbitraire, (iii) les couplages implicites et semi-implicites pour l'interaction fluide-structure. Nous proposons également une nouvelle formulation de la méthode de la frontière élargie visant à modéliser le transport de particules déformables immergées dans un fluide. Nos simulations numériques se sont appuyées sur la librairie de calcul Feel++, spécialisée dans la résolution d'EDP. Outre l'implémentation des modèles physiques, nous y avons développé diverses fonctionnalités nécessaires à la mise en oeuvre de nos méthodes : interpolation, méthode de Galerkin non standard, méthode ALE, environnement pour l'interaction fluide-structure. De plus, de par la taille des géométries et la complexité des modèles mis en jeu, le passage au calcul parallèle a été indispensable pour pouvoir réaliser nos simulations. Ainsi, nous avons décrit le développement qui a été effectué dans cette librairie pour permettre le déploiement de nos programmes sur des architectures parallèles. / Unlike ordinary liquids, complex fluids such as blood exhibit strange behavior mainly dependent underlying structures that compose them. Simulation of blood flow remains a formidable challenge for theoretical and numerical modeling whose interest is to develop methods and simulation tools for the medical community. We propose in this work a contribution to this project will be mainly focused on numerical and computational aspects. We are particularly interested in the interaction between the blood and the vascular wall, which plays an important role in the large arteries as the aorta. We also invested in transport simulation of blood cells in the blood. For solving partial differential equations describing our hemodynamic models we chose to use numerical methods whose accuracy can be increased arbitrarily. For this purpose, the main ingredients that have been used are textit {(i)} the finite element method based on Galerkin approximations of arbitrary order in space and geometry, (i) the ALE method for taking into account the mobility of areas for movements of arbitrary order, (ii) the implicit and semi-implicit coupling for fluid-structure interaction. We also propose a new formulation of the method of boundary extended to model the transport of deformable particles immersed in a fluid. Our numerical simulations were based on the library Feel++ specializing on the EDP resolution. In addition to implementing physical models, we have developed various features are necessary for the implementation our methods: interpolation, Galerkin non-standard method, ALE method, framework for fluid-structure interaction. In addition, the size and geometry of the complexity of the models involved, the transition to parallel computing was essential in order to achieve our simulations. Thus, we have described the development was carried out in the library to allow deployment of our programs on parallel architectures.

Dans la conception des turbomachines modernes, la prédiction des phénomènes aéroélastiques est devenue un point clé. La tendance à réduire la masse et à augmenter la charge des composants aérodynamiques accroit le risque de rupture. Dans un tel contexte, la compréhension et la bonne prédiction des diverses instabilités constituent un enjeu industriel et scientifique majeur. Le présent travail de recherche a pour objectif d’améliorer la prédiction des phénomènes instationnaires intervenant dans les problèmes d’aéroélasticité en turbomachines. Cette thèse est plus particulièrement axée sur la simulation de l’interaction onde de choc/couche limite. Le support d’étude est une tuyère transsonique présentant un écoulement avec des zones décollées. L’oscillation forcée de l’onde de choc est simulée grâce à une méthode de petites perturbations instationnaires couplée avec une hypothèse de turbulence variable. Cette approche est validée par comparaison avec des mesures. Elle permet une prédiction tout à fait satisfaisante du premier harmonique de pression sur la paroi de la tuyère. Ce travail a montré la nécessité de linéariser le modèle de turbulence. Le besoin de dériver le modèle de turbulence nous a amené à investiguer la modélisation faite pour prédire l’interaction onde de choc/couche limite. Un modèle de turbulence à deux équations complété par une équation de « retard » est implémenté afin de capter un déséquilibre de la turbulence. Les résultats obtenus en tuyère sont cohérents avec la théorie mais une surproduction d’énergie turbulente en présence de bord d’attaque rend le modèle inefficace pour des configurations de turbomachines. Au final, l’introduction d’un limiteur de viscosité turbulente dans un modèle de turbulence à deux équations s’avère donner de bons résultats. La méthode de dérivation du modèle est alors présentée sur le modèle de Wilcox proposé en 2008. Enfin, la technique de linéarisation est étendue à la problématique aéroélastique. Une approche de couplage fluide-structure faible est adoptée. L’oscillation structurelle des aubages suivant les modes propres est considérée mais en laissant la fréquence évoluer au cours du couplage. La nouvelle méthode utilisée s’appuie sur la construction d’un méta-modèle du comportement dynamique du fluide afin de résoudre directement le système fluide-structure couplé. Cette technique est validée sur une configuration de grille annulaire de turbine en haut subsonique et présente l’avantage d’un temps de calcul réduit. / In modern turbomachinery design, predicting aerolastic phenomena has become a key point. The development of highly loaded components, while reducing their weight, increases the risk of failure. In this context, good understanding and prediction of various instabilities are a major industrial and scientific challenge. This research work aims to improve the prediction of unsteady phenomena involved in turbomachinery aeroelasticity. This study focuses especially on the simulation of shock wave/boundary layer interaction. To begin with, a transonic nozzle separated flow is investigated. Forced oscillation of the shock wave system is simulated through a small unsteady perturbation method combined with the assumption of variable turbulence. This approach is validated against exprimental measurements. The first harmonic of pressure on the wall of the nozzle is predicted quite satisfactorily. The need to linearize the turbulence model was shown of high importance. Deriving the turbulence model, leads us to investigate the turbulence modeling performed to predict the shockwave/boundary layer interaction. A two equations turbulence model supplemented by a "time-lagged" equation is implemented to capture non-equilibrium effects of turbulence. All achieved results for a nozzle are consistent with theory, but overproduction of turbulent kinetic energy at leading edge makes the model useless for turbomachinery configurations. However, the introduction of an eddy viscosity stress limiter inside a two-equation turbulence model proves to give good results. The derivation method is thus presented on such a model, precisely on Wilcox model proposed in 2008. Finally, the linearization technique is extended to aeroelastic problems. A loose fluid-structure coupling strategy is adopted. The structural oscillation of the blades is considered for eigen-modes but frequency is free to change during coupling resolution. The new approach is based on the building of a meta-model to describe the fluid dynamic behavior in order to solve directly the coupled fluid-structure system. This technique is validated on a standard high subsonic turbine configuration and takes advantage of a reduced computation time.

Ce travail est consacré à l’étude numérique de l’interaction entre un fluide compressible et une structure indéformable, en adaptant une famille récente de schémas d’ordre très élevé à la prise en compte de conditions aux bords particulières entre le fluide et la structure. Plus précisément,on évalue l’apport de schémas d’ordre strictement supérieur à 3 par rapport à des stratégies plus classiques dans la littérature restreintes aux ordres 1 et 2. Un résultat important est qu’il est possible de réaliser le couplage à tout ordre et qu’il existe des configurations pour lesquelles on observe un gain important pour les ordres élevés. Une revue bibliographique est faite rappelant les résultats théoriques concernant les systèmes hyperboliques et décrivant les méthodes utilisées dans la littérature pour la simulation de la dynamique des gaz et la prise en compte des conditions aux bords. Un schéma sur grilles cartésiennes décalées et d’ordre très élevé est proposé pour la résolution des équations d’Euler en 1D/2D. Ce schéma est basé sur le formalisme Lagrange-projection et bien que formulé en énergie interne assure conservation et consistance faible grâce à un correctif en énergie interne. Parallèlement, l’étude pour les systèmes hyperboliques linéaires de discrétisation à l’ordre très élevé des conditions aux bords est faite. Elle met en évidence la nécessité pour l’ordre élevé de s’intéresser à la stabilité des schémas ainsi obtenus. À partir de ces travaux, la prise en compte de conditions aux bords en vitesse normale imposée est réalisée pour les équations d’Euler en 1D et 2D. Enfin, une procédure de couplage entre fluide compressible et structure indéformable est proposée. / This work is devoted to the construction of stable and high-order numerical methods in order to simulate fluid - rigid body interactions. In this manuscript, a bibliographic overview is done, which highlights theoretical results about hyperbolic system of conservation laws, as well as the methods available in the literature for the hydrodynamics simulation and the numericalboundary treatment. A high-order accurate scheme is proposed on staggered Cartesian grids to approximate the solution of Euler equations in 1D and 2D. The scheme relies on Lagrange-remap formalism, and although formulated in internal energy, ensures both conservation and weak consistency thanks to an internal energy corrector. In the same time, the study of high-order numerical boundary treatment for linear hyperbolic system is done. It highlights the necessity to focus especially on the linear stability of the effective scheme. Starting from the linear results, the numerical boundary treatment with imposed normal velocity is done for Euler equations in 1D and 2D. Last, the coupling between a compressible fluid and a rigid body is realized, using the designed procedure for numerical boudary treatment.

La sûreté des réacteurs nucléaires constitue une condition primordiale et nécessaire à leur industrialisation. Pour les réacteurs à neutrons rapides refroidis au sodium, cette sûreté passe par la possibilité de détection d’une fuite d’eau dans le sodium au niveau du générateur de vapeur à des stades très précoces de leur déclenchement. Le présent travail de thèse s’inscrit dans le cadre du développement d’une technique non intrusive pour la détection d’une réaction sodium-eau dans le générateur de vapeur d’un réacteur nucléaire refroidi au sodium. On désire identifier le bruit vibratoire de cette réaction à partir de mesures d’accélération sur la virole externe du générateur. Cependant, les vibrations dues à la fuite peuvent être noyées dans le bruit généré par l’écoulement du sodium ou par d’autres sources d’excitation. Afin d’augmenter, le rapport signal à bruit, on s’intéresse à l’apport d’un filtrage spatial par formation de voies pour ce type de système. Ce dernier se caractérise notamment par la présence d’un couplage structure-fluide lourd (c.-à-d. virole-sodium) et des modes vibroacoustiques qui apparaissent aux mêmes fréquences d’intérêt pour la détection de la fuite. Pour réaliser cette étude, nous considérons une maquette expérimentale composée d’une conduite cylindrique reliée à un circuit hydraulique par deux brides très rigides. La source à identifier est simulée par un hydrophone de taille réduite en mode émetteur positionné à l’intérieur de la conduite par un dispositif mécanique dédié à cet effet. Le bruit perturbateur est induit par l’écoulement de l’eau, supposé turbulent à partir d’un certain débit. L’antenne est composée de vingt-cinq accéléromètres placés sur la conduite d’essai. L’objectif consiste à traiter simultanément les signaux accéléromètriques pour faire ressortir le signal de la fuite tout en rejetant le bruit dû à l’écoulement. Deux types de traitement par formation de voies sont considérés : la formation de voies conventionnelle (dites de Bartlett) et une formation de voies optimisée qui vise à maximiser le rapport signal à bruit de la fuite. / The safety of nuclear reactors represents a necessary and vital condition for the exploitation of nuclear plants with liquid-cooled cores. This safety passes by the ability to detect and anticipate the earliest stages of a water leak into sodium within the heat generator. The study detailed in this manuscript has been initiated in a framework aiming to develop nonintrusive detection techniques for sodium water chemical reactions. Its main goal is to uncover the signal of a water leak into sodium, based on the measurements of the vibratory field recorded by the means of accelerometers externally mounted on the heat generator shell. However, such a spectrum is masked by the plant general background noise, especially that generated by heavy fluid flow (i.e. sodium flow) during actual operating times of the power plant. In order to increase the signal-to-noise ratio of the leak, beamforming technique for the acceleration measurements of the mechanical system is considered. The aforementioned system is characterized by a strong non-linear coupling between the the heat generator’s cylindrical shell and the heavy fluid in motion. In particular, fluid motion and the acoustic emission of the leak seem to excite some eigen modes of the cylinder, at these same frequencies where the acoustic signature of the leak is at its highest amplitude. For the purpose of our study, a cylindrical mock-up connected by some very rigid links to a hydraulic circuit is considered. A hydrophone emission excites the mock-up from within and is being accounted for the acoustic leak. Water flows inside the cylinder at turbulent Reynolds number. An array of twenty-five accelerometers is mounted on the mock-up shell using ceramic insulators. The main goal is to combine the twenty-five signals in such a way that allows the increase of the SNR for the acoustic source while rejecting water flow noise. Two beamforming techniques are applied and compared: classical Bartlett beamforming as well as optimized beamforming for SNR maximization (Max SN

Le but du travail présenté est la compréhension puis la simulation numérique des mécanismes d'érosion des augets de turbine Pelton par impacts répétés de gouttes d'eau dans le but de prédire la durée de vie des composants. Tout d'abord, les phénomènes de propagation d'ondes dans les milieux fluide et solide sont étudiés. Cela permet de mettre en lumière l'évolution temporelle et la distribution spatiale des pressions de contact, et l'apparition de microjets par éjection supersonique du fluide au contact. Les études expérimentales de l'érosion par gouttes d'eau traduisent un dommage basé sur la fissuration par fatigue. Des simulations numériques en dynamique rapide couplées fluide-structure sont alors effectuées. Le domaine solide est discrétisé par la Méthode des éléments Finis (MEF), et le domaine fluide par la méthode Smoothed Particle Hydrodynamics (SPH), qui est une méthode particulaire (sans maillage) particulièrement adaptée aux grandes distorsions et au suivi des surfaces libres. L'analyse des états de contraintes vient corroborer la nature cyclique de l'endommagement. La simulation d'érosion est alors réalisée à l'aide de critères de fatigue multiaxiaux. Le choix se porte vers un premier critère général de l'American Society of Mechanical Engineers (ASME), utilisant les valeurs principales des différences de contraintes au cours du temps. Le second choix concerne un critère à plan critique : le critère de Dang Van 2. Il traite séparément la contrainte hydrostatique et le cisaillement alterné maximal local. Ces critères permettent de définir les régions érodées du solide au bout d'un nombre d'impact donné, ce qui fait de cette démarche une méthode prédictive. Une étude paramétrique pour différentes tailles de gouttes et vitesses d'impact est ensuite réalisée, puis on évalue l'influence de la présence d'une couche de coating. / The goal of this work is the comprehension and the numerical simulation of erosion caused by repeated droplet impact on Pelton turbine buckets, to predict the lifetime of these components. First, waves propagation phenomenon inside fluid and solid domains are presented, which allows determining the time evolution and spatial distribution of contact pressure, and the birth of lateral microjets by supersonic ejection of the fluid on the contact. Experimental studies of erosion by droplet impact highlight a fatigue cracking-based erosion mechanism. Then, coupled FSI computation are performed. The solid subdomain is discretized by the Finite Element Method (FEM), and the fluid subdomain by the Smoothed Particle Hydrodynamics (SPH), which is a particle method (meshless) effectively recommended for large distortions and free surface tracking. Stress analysis confirms the cyclic nature of the damage mechanism, and erosion simulation is performed using multiaxial fatigue criteria. The first selected criterion is a general one from the American Society of Mechanical Engineers (ASME) using principal values of stress differences over time. The second one is the Dang van 2 criterion, belonging to the family of critical plane criteria. This criterion considers separately the effects due to hydrostatic stress on one hand, and the ones induced by maximum local shear on the other. These two criteria are used to building the equivalent eroded zones of the solid subdomain for a given number of impacts, which allows to qualify this procedure as a predictive predictive. Finally, a parametric study for different droplet sizes and velocites is computed, and the effects of a coating layer are investigated.

La méthode des mobilités énergétiques a été développée dans le domaine des moyennes fréquences, afin de prévoir le comportement vibratoire des structures par bandes de fréquence. Dans ce travail, on étend l'approche du cas du couplage vibroacoustique entre un fluide et une structure. Les mobilités acoustiques et structurales sont étudiées pour le milieu fluide et la structure respectivement. La surface du couplage est maillée en pavés, afin de réduire le problème en mobilités énergétiques ponctuelles. Une simulation numérique a été réalisée dans le cas d'une plaque couplée à une cavité parallélépipédique, où les influences du nombre de pavés et de l'amortissement des sous-systèmes sont étudiées. L'énergie cinétique de la plaque, la puissance transmise au milieu acoustique et le facteur de rayonnement de la structure sont calculés et comparés avec un calcul de référence (méthode modale) et avec des résultats expérimentaux. Les critères de maillage (taille de pavé par rapport aux longueurs d'ondes de la structure et du milieu fluide) sont déterminés avant et après la fréquence critique, pour assurer une modélisation correcte du couplage dans une bande de fréquence donnée. Afin de tenir compte d'une structure compliquée à modéliser, un calcul mixte (utilisation des mobilités mécaniques mesurées et des mobilités acoustiques calculées) est également présenté.

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.

Cette thèse consiste à mettre en place un couplage externe aéro-thermo-mécanique, sur la base d'un schéma partitionné, entre les codes de recherche CEDRE (mécanique des fluides, volumes finis) et Z-set (modules indépendants pour la mécanique des structures et la thermique du solide, éléments finis). Les résultats numériques sont confrontés à ceux de l'expérience (une campagne de mesures a été menée dans le cadre de cette étude), dans le cas d'un problème complexe lié au domaine de l'aérospatial : l'interaction flamme-paroi. Ce phénomène est piloté par la thermique, à travers le flux de chaleur pariétal généré par la flamme. A cause de la disparité des temps caractéristiques thermiques entre les milieux fluide et solide, la partie aéro-thermique du couplage est traitée de façon simplifiée, en considérant le fluide comme une suite d'états stationnaires. L'échauffement de la plaque métallique provoque sa déformation (la loi de comportement mécanique du matériau est de type élasto-visco-plastique). Le déplacement de l'interface fluide-structure est propagé sur le maillage fluide. En se basant sur les similitudes entre jets non réactifs et réactifs (de type flamme) dans le cas de l'impact, des calculs couplés sont menés dans des configurations 2D et 3D de l'impact d'un jet chaud non réactif. / This thesis consists in setting up an external fluid-thermal-structural coupling, based on a partitionned scheme, between the research codes CEDRE (fluid mechanics, finite volumes) and Z-set (independent solvers for structural mechanics and heat transfer through the solid). The numerical results are compared with experimental data, to study a complex problem related to the aerospace certification process: the flame-wall interaction. This phenomenon is is driven by the heat flux generated by the flame close to the wall. Because of the disparity of thermal characteristic times between the fluid and the solid, the aero-thermal part of the coupling is simplified by considering the fluid as a sequence of steady states. The heating of the metallic plate causes its deformation (the material has a viscoplastic behavior). The displacement of the fluid-structure interface is propagated through the fluid mesh. Based on similitudes between impinging reacting jets (flames) and non-reacting jets, coupled computations are performed in 2D and 3D configurations with an equivalent non-reacting hot jet.

Le projet CLIC (City Lightweight Innovative Cab) dans lequel s’inscrit cette thèse de doctorat vise à développer une cabine de camion allégée sans dégrader les performances vibratoires et acoustiques. Pour cela il est nécessaire d’établir dans un premier temps un modèle de prédiction vibroacoustique du système couplé structure/espace intérieur incluant l’influence des matériaux dissipatifs (amortissement ou absorption) dans le domaine des moyennes fréquences. Les méthodes basées sur les éléments finis et les approches statistiques les plus couramment utilisées étant peu adaptées pour ce domaine de fréquence (coût de calcul important, méthodes peu flexibles), nous utiliserons le formalisme de la méthode SmEdA (Statistical modal Energy distribution Analysis). L’objectif principal de cette thèse de doctorat est dès lors, d’étendre cette méthode à la prise en compte de l’effet d’amortissement induit par des matériaux dissipatifs. La méthodologie se divise en trois étapes: 1. Les modèles équivalents des matériaux dissipatifs sont établis: (1) un modèle de plaque équivalent pour décrire la plaque amortie par un ou plusieurs patch(s) viscoélastique(s) et (2) un modèle de fluide équivalent pour décrire un matériau poreux agissant dans la cavité. 2. Chaque sous-système amorti est modélisé par éléments finis. Les méthodes MSE (Modal Strain Energy) et MSKE (Modal Strain Kinetic Energy) sont ensuite utilisées pour estimer les facteurs de perte modaux de chaque sous-système. 3. Le calcul SmEdA est effectué sur le système couplé en prenant en compte les facteurs de pertes modaux de chaque sous-système estimés dans la deuxième étape. Le point d’excitation est appliquée à la plaque, en supposant la force stationnaire et large bande. Afin de valider la méthodologie proposée un cas semi-complexe composé d’une plaque rectangulaire couplée à une cavité parallélépipédique est considéré. Ce système peut être utilisé pour étudier l’interaction vibroacoustique entre la structure de la cabine et l’intérieur de l’habitacle. Deux cas d’amortissement sont étudiés pour le système semi-complexe plaque-cavité: (1) un cas où la plaque est amortie avec un (ou plusieurs) patch(s) viscoélastique(s) et (2) un cas où un matériau poreux est placé dans la cavité. Le problème vibroacoustique est pour chaque cas modélisé suivant les trois étapes proposées et analyses dans le formalisme de la méthode SmEdA. Les résultats sont ensuite comparés au cas de référence (sans matériau dissipatif). La dernière partie de la thèse porte sur la validation expérimentale pour chaque cas test de la méthodologie numérique proposée. a mobilité d’éntrée, la puissance injectée et les énergies des sous-systèmes sont comparées aux prédictions numériques. Enfin les facteurs de pertes modaux des sous-systèmes estimés par les méthodes MSE et MSKE sont comparés aux résultats expérimentaux obtenus par la méthode d’analyse modale à haute résolution (méthode ESPRIT). / The project CLIC (City Lightweight Innovative Cab) aims to develop a lighter-weighted truck that maintains NVH performances of the initial design. This PhD research is then to establish a vibroacoustic prediction model of a complex structure-bounded fluid system (cabin structure coupled to cabin space) including dissipative treatments (damping or absorbing materials) for the mid-frequency domain. Since most commonly used element based and statistical methods are not suitable for this frequency domain, a proper prediction tool, which should be flexible in modeling capabilities and feasible in computational cost, must be implemented. The SmEdA (Statistical modal Energy distribution Analysis) method is considered in this thesis to comply with these requirements. The main objective of this research is to extend this method for taking account of the damping effect induced by dissipative materials. Development and validation of the methodology are carried out. 1. Dissipative materials are represented by simplified equivalent models: (1) the equivalent single layer model for describing the plate covered with a viscoelastic layer and (2) the equivalent fluid model for describing a porous material into the cavity. 2. Each subsystem including the equivalent models of the dissipative materials is modeled with FEM(Finite Element Model). The FE matrices including the energy dissipation are then computed. The MSE (Modal Strain Energy) and MSKE (Modal Strain Kinetic Energy) methods are used to estimate the modal damping loss factor of each subsystem mode. 3. The SmEdA calculation is performed on a whole system considering the modal damping loss factors estimated in the second step for each subsystem. The power is injected into the plate at a localized point by the stationary white noise force and subsequently, the SmEdA parameters are computed. To validate the proposed methodology, laboratory test cases of the structure-fluid problem composed of a rectangular plate coupled to a parallelepipedic cavity are considered. Such system can be used to study the vibroacoustic interaction between structure and fluid. Two damped test cases of the plate-cavity system are studied: (1) a system with a viscoelastic damping pad on the plate and (2) a system with a composite fibre in the cavity. The damped test cases are modeled following the three steps and are analyzed in the framework of SmEdA. The results are then compared to the original case with no damping treatment. The last part of the thesis presents an experimental validation of the numerical computation results on each test case. Measured quantities such as input mobility, injected power and subsystem energies are compared to the numerical predictions. The modal damping loss factors of the damped subsystems estimated with MSE and MSKE methods are compared to the experimental results estimated by a high-resolution modal analysis method (ESPRIT method).

La simulation de cas réalistes d'écoulements ou de transferts thermiques implique souvent l'utilisation d'obstacles ou d'interfaces de forme complexe. De part leur manque de flexibilité, les maillages structurés ne sont pas initialement adaptés au traitement d'interfaces irrégulières, ces dernières coïncidant rarement avec les lignes du maillage. Afin de permettre à l'approche structurée de traiter des interfaces complexes avec précision, des méthodes dites de domaines fictifs sont nécessaires. La première contribution de cette thèse est une nouvelle méthode de travail sur maillage curviligne structuré qui permet de réutiliser de nombreuses méthodes fonctionnant initialement sur des maillages cartésiens sur maillages curvilignes. Nous avons ensuite mis au point deux nouvelles méthodes de domaines fictifs : la méthode de pénalisation de sous-maille (PSM) pour la gestion des frontières immergées pour les équations elliptiques et de Navier-Stokes et la méthode d'interface immergée algébrique (IIA) pour les problèmes d'interfaces immergées pour les équations elliptiques. L'un des intérêts de ces deux méthodes à l'ordre deux en espace est leur simplicité. Ces différents développements ont finalement été appliqués à des cas de couplage fluide-structure académiques et réalistes (sédimentation d'un cylindre, hydroplanage d'un pneu, écoulements dans une tête de forage et convection naturelle dans la grotte de Lascaux).