Le récent développement des menaces terroristes renforce l'effort de recherche du CEA et d'EDF pour la protection des citoyens et des installations. De nombreux scénarios doivent être envisagés comme, par exemple, la chute d'un avion de ligne sur une structure de génie civil. La dispersion du carburant dans la structure, son embrasement sous forme de boule de feu et les effets thermiques associés sont des éléments essentiels du problème. L'utilisation de modèles numériques est indispensable car des expériences seraient difficiles à mettre en œuvre, coûteuses et dangereuses. Le problème type que l'on cherche à modéliser est donc l'impact d'un réservoir rempli de fluide, sa déchirure et la dispersion de son contenu. C'est un problème complexe qui fait intervenir une structure mince avec un comportement fortement non-linéaire allant jusqu'à rupture, un fluide dont la surface libre peut varier drastiquement et des interactions fluide-structure non permanentes. L'utilisation des méthodes numériques traditionnelles pour résoudre ce problème semble difficile, essentiellement parce qu'elles reposent sur un maillage. Cela complique la gestion des grandes déformations, la modélisation des interfaces variables et l'introduction de discontinuités telles que les fissures. Afin de s'affranchir de ces problèmes, la méthode sans maillage SPH (\og Smoothed Particle Hydrodynamics \fg) a été utilisée pour modéliser le fluide et la structure. Ce travail, inscrit dans la continuité de recherches précédentes, a permis d'étendre un modèle de coque SPH à la modélisation des ruptures. Un algorithme de gestion des interactions fluide-structure a également été adapté à la topologie particulière des coques. Afin de réduire les coûts de calcul importants liés à ce modèle, un couplage avec la méthode des éléments finis a également été élaboré. Il permet de n'utiliser les SPH que dans les zones d'intérêt où la rupture est attendue. Finalement, des essais réalisés par l'ONERA sont étudiés pour valider la méthode. Ces travaux ont permis de doter le logiciel de dynamique rapide Europlexus d'un outil original et efficace pour la simulation des impacts de structures minces en interaction avec un fluide. Un calcul démonstratif montre enfin la pertinence de l'approche et sa mise en œuvre dans un cadre industriel.

Cette thèse a pour objet l'étude du comportement asymptotique des solutions des équations des fluides de grades 2 et 3. Dans le premier chapitre, on étudie les profils asymptotiques au premier ordre des solutions des équations des fluides de grade 2 en dimension 3. On démontre que les solutions des équations des fluides de grade 2 convergent vers des solutions particulières et explicites des équations de la chaleur, lorsque le temps tend vers l'infini. Ce résultat montre en particulier que les fluides de grade 2 se comportent asymptotiquement comme les fluides newtoniens régis par les équations de Navier-Stokes. Pour cette étude, on utilise les variables d'échelles (ou variables autosimilaires), et on effectue des estimations d'énergies dans divers espaces fonctionnels, en particulier dans des espaces de Sobolev à poids polynomiaux. La description des profils asymptotiques est obtenue sous des conditions de petitesse sur les données initiales de l'équation.Le second chapitre de cette thèse traite des profils asymptotiques à l'ordre 1 des solutions des équations des fluides de grade 3 en dimension 2. À l'instar des résultats du premier chapitre, on obtient ici aussi la convergence des solutions de ces équations vers des solutions explicites des équations de la chaleur. Les outils utilisés pour cette étude sont semblables à ceux utilisés pour les fluides de grade 2 en dimension 3, à savoir les variables autosimilaires et des estimations d'énergies. Dans ce cas aussi, on conclut que les fluides de grade 3 se comportent asymptotiquement comme les fluides newtoniens.Dans le dernier chapitre, on étudie l'existence d'un attracteur pour les équations des fluides de grade 3 en dimension 2 avec des conditions périodiques. On considère donc les solutions faibles de ces équations à données initiales dans l'espace de Sobolev H¹. Ces solutions faibles définissent un semi-groupe généralisé. Ensuite, on montre que les solutions à données initiales dans H² possèdent un attracteur global pour la topologie H¹. Pour ce travail, on utilise un schéma de Galerkin, des estimations a priori et une méthode de monotonie. Les principales difficultés que l'on rencontre sont liées au peu de régularité des données initiales et au fait que l'on ne sait par si les solutions des équations des fluides de grade 3 à données H¹ sont uniques.

Afin de contribuer à l'amélioration des simulations numériques de foyers industriels avec recycle de gaz brûlés et combustible liquide, la modélisation de la combustion turbulente non prémélangée est abordée sous deux de ses aspects : le traitement des problèmes diphasiques et la prise en compte des configurations complexes de mélange (dilution, recirculation interne).Une flamme spray éthanol/oxygène diluée par du dioxyde de carbone est d'abord calculée en LES dans un formalisme Euler-Lagrange et une chimie détaillée pour une résolution atteignant 250 µm. Les conditions d'injection du brouillard de gouttes sont déterminées à partir de mesures expérimentales de granulométrie et d'anémométrie phase Doppler.Dans une seconde partie, un nouveau formalisme de description de la combustion turbulente, basé sur l'introduction de temps caractéristiques de l'histoire du mélange est développé, puis validé en LES sur le cas de référence d'une flamme jet dans un écoulement co-courant vicié.

Cette thèse est consacrée à l'étude numérique des interactions non-linéaires entre des vagues et un corps rigide perçant la surface libre. La méthode développée repose sur un modèle d'éléments de frontière qui réduit la dimensionnalité du problème d'une dimension. Dans un premier temps, un modèle2D est appliqué à des géométries simples et permet de démontrer la pertinence de l'approche envisagée pour la prédiction des mouvements d'une structure flottante soumise à des vagues monochromatiques régulières. Dans un second temps, en nous inspirant d'un modèle potentiel non-linéaire 3D développé par Grilli textit{et al.}~cite{grilli2001fully}, nous proposons une généralisation de la méthode pour des maillages triangulaires non-structurés de surfaces 3D. Le modèle développé permet de traiter des configurations arbitraires de plusieurs cylindres verticaux en interaction avec les vagues. Nous présentons des cas de validation de nature académique qui permettent d'apprécier le comportement du modèle numérique. Puis nous nous tournons vers l'application visée par EDF R&D, qui concerne le dimensionnement d'éoliennes off-shore flottantes. Un flotteur de type semi-submersible est évalué à l'aide du modèle non-linéaire / This PhD work is devoted to the study of nonlinear interactions between waves and floating rigid structures. The developed model relies on a boundary element method which reduces the dimensionality of the problem by one. First, a 2D model is applied to basic geometries and allows us to demonstrate the validity of the method for predicting the motion of a floating structrure subject to incoming monochromatic regular waves. Secondly, getting inspired by the 3D fully nonlinear potential flow model of Grilli textit{et al.}~cite{grilli2001fully}, we propose a novel model which generalizes the method for unstructured triangular meshes of 3D surfaces. The proposed model is able to deal with arbitrary configurations of multiple vertical cylinders interacting with the waves. We present academic validation test cases which show how the model works and behaves. Finally, we study situations of interest for EDF R&D related to floating off-shore wind turbines. A semi-submersible platform is evaluated with the nonlinear model

Mon travail de thèse a porté sur la compréhension «systèmes dynamiques de la dynamique à grande échelle dans l’écoulement pleinement développé de cisaillement turbulent. Dans le plan écoulement de Couette, simulation des grandes échelles (LES) est utilisée pour modéliser petits mouvements d’échelle et de ne résoudre mouvements à grande échelle afin de calculer non linéaire ondes progressives (SNT) et orbites périodiques relatives (RPO). Artificiel sur-amortissement a été utilisé pour étancher une gamme croissante de petite échelle motions et prouvent que les motions grande échelle sont auto-entretenue. Les solutions d’onde inférieure branche itinérantes qui se trouvent sur le bassin laminaire turbulent limite sont obtenues pour ces simulation sur-amortie et continue encore dans l’espace de paramètre à des solutions de branche supérieure. Cette approche ne aurait pas été possible si, comme supposé dans certains enquêtes précédentes, les mouvements à grande échelle dans le mur bornées flux de cisaillement sont forcée par un mécanisme fondé sur l’existence de structures actives à plus petite échelle. En flux Poseuille, orbites périodiques relatives à décalage réflexion symétrie sur la limite du bassin laminaire turbulent sont calculés en utilisant DNS. Nous montrons que le RPO trouvé sont connectés à la paire de voyager vague (TW) solution via bifurcation mondiale (noeud-col-période infinie bifurcation). La branche inférieure de cette solution TW évoluer dans un état de l’envergure localisée lorsque le domaine de l’envergure est augmentée. La solution de branche supérieure développe plusieurs stries avec un espacement de l’envergure compatible avec des mouvements à grande échelle en régime turbulent. / My thesis work focused on ‘dynamical systems’ understanding of the large-scale dynamics in fully developed turbulent shear flow. In plane Couette flow, large-eddy simulation (L.E.S) is used to model small scale motions and to only resolve large-scale motions in order to compute nonlinear traveling waves (NTW) and relative periodic orbits (RPO). Artificial over-damping has been used to quench an increasing range of small-scale motions and prove that the motions in large-scale are self-sustained. The lower-branch traveling wave solutions that lie on laminar-turbulent basin boundary are obtained for these over-damped simulation and further continued in parameter space to upper branch solutions. This approach would not have been possible if, as conjectured in some previous investigations, large-scale motions in wall bounded shear flows are forced by mechanism based on the existence of active structures at smaller scales. In Poseuille flow, relative periodic orbits with shift-reflection symmetry on the laminar-turbulent basin boundary are computed using DNS. We show that the found RPO are connected to the pair of traveling wave (TW) solution via global bifurcation (saddle-node-infinite period bifurcation). The lower branch of this TW solution evolve into a spanwise localized state when the spanwise domain is increased. The upper branch solution develops multiple streaks with spanwise spacing consistent with large-scale motions in turbulent regime.

Le bruit émis par les nouvelles architectures de moteurs aéronautiques a été considérablement réduit dans les dernières années. Les différentes sources de bruit ont été identifiées et pour la plupart réduites. Cependant, la contribution relative du bruit de combustion au bruit global a augmenté progressivement avec la décroissance des autres sources. Deux mécanismes de génération de bruit de combustion ont été identifiés : le bruit direct qui est produit par des fluctuations du dégagement de chaleur dû à la combustion, et le bruit indirect qui est généré par l’accélération des spots d’entropie. Dans ce travail, les mécanismes de génération et propagation du bruit entropique sont étudiés par des simulations numériques aux grandes échelles (en anglais LES) et par des modèles analytiques. Dans un premier temps, une configuration simplifiée du phénomène est étudiée : des spots d’entropie sont créés par des résistances chauffantes et ensuite accélérés par une tuyère pour générer du bruit indirect. Cette configuration a été simulée et ses résultats validés par des campagnes expérimentales. Ensuite, la simulation numérique est utilisée pour mieux comprendre les mécanismes de génération du bruit indirect et ses interactions avec des effets visqueux et non visqueux. Dans une seconde partie, une configuration de turbine haute pression à un seul étage est utilisée pour étudier le bruit indirect d’une façon plus réaliste. Dans les deux parties de cette thèse, les résultats numériques sont comparés à des théories analytiques pour mieux comprendre les avantages et inconvénients d’une méthode par rapport à l’autre. / Combustion noise is increasing its relative contribution to aircraft noise, while other sources are being reduced and new low-NOx emission combustion chambers being built. Two mechanisms are responsible for this noise source: direct noise in which acoustic waves are generated by the flame and propagate to the outlet of the aero-engine, and indirect noise, where entropy waves generate noise as they are accelerated and decelerated in the turbine stages. In this work, the analytical models used for the propagation of waves through non-homogeneous flows, including the generation of indirect noise, are revised and extended. In the first part, the quasi-1D case is studied, extending the analytical method to non-zero frequencies and validating the results with numerical methods and experimental data. In the second part, the 2D method for the case of compact turbine blades is studied and validated using numerical simulations of a rotating blade and of a complete turbine stage. Finally, in the third part of this thesis, these models are combined with reactive and compressible Large Eddy Simulations (LES) of combustion chambers to build a hybrid approach, named CHORUS, able to predict combustion noise.

Les préoccupations relatives à la consommation de l’énergie et notamment des pertes provoquent une demande d’optimisation toujours plus forte des procédés. La recherche de l’efficacité maximale dans les échangeurs de chaleur est d’autant plus forte lorsque ceux-ci ne fonctionnent pas au régime nominal. L’étude présentée dans ce mémoire se concentre sur une problématique de l’efficacité des échangeurs thermiques : la distribution diphasique. Le but de ce travail a été de montrer la faisabilité d’un dispositif améliorant la répartition des deux phases dans un échangeur à mini canaux. La méthode choisie a été l’introduction d’ultrasons par le biais de générateurs placés dans le distributeur produisant une fontaine. Cette technologie présente l’avantage d’être facilement modulable en ne changeant que la tension d’entrée des appareils. Ceci permet donc de pouvoir adapter le dispositif aux différentes conditions opératoires. Des essais expérimentaux ont été réalisés dans différentes conditions, en changeant notamment le débit de chaque phase. La densité de flux massique a été variée de 60 kg.m-2.s-1 à 450 kg.m-2.s-1 et le titre massique de moins de 1% à plus de 23% en entrée de section d’essais. Les comparaisons entre les différents essais se sont faite par la mesure des débits de chaque phase sortant des canaux. L’introduction des ultrasons a également été étudiée en jouant sur le placement et le nombre des générateurs d’ultrasons. Cette technologie s’est alors montré viable car la distribution diphasique est améliorée dans une grande majorité des cas testés. Des phénomènes intéressant ont même été observés, permettant lors d’une prochaine étude, des approches différentes mais complémentaires. / Concerns about energy consumption including thermal losses cause a need in greater optimization of the processes. The goal of the maximum efficiency in heat exchangers is even stronger when they do not work at nominal conditions. The study presented in this paper focuses on the two-phase flow distribution problematic. The aim of this work was to demonstrate the feasibility of a device that can improve the distribution of the two phases in a mini-channel heat exchanger. The method chosen was the introduction of ultrasound through generators placed in the distributor where they can produce an ultrasonic fountain. This technology has the advantage of being easily adjustable by changing the input voltage of the devices. This therefore allows the device to adapt to different operating conditions. Experimental tests have been done under different conditions, in particular by changing the flow rate of each phase. The mass flow density was varied from 60 kg.m-2.s-1 to 450 kg.m-2.s-1 and the mass quality of less than 1% to over 23% at the tests section inlet. Comparisons between different tests were made by measuring flow rates of each phase at the outlet of the channels. The introduction of ultrasound was also studied by varying the placement and number of the generators of ultrasound. This technology has been shown as viable two-phase distribution enhancement in a majority of cases tested. Interesting phenomena were even observed, which allows, for a future study, different but complementary approaches.

Résumé : Les nanotubes de carbone mono-parois (SWCNTs) sont très peu dispersibles dans les solvants et ils ont besoin d'être chimiquement modifiés avant leur utilisation dans beaucoup d'applications. Ce travail se concentre sur la synthèse du matériau des SWCNTs chimiquement modifié par une approche in situ. Les objectifs principaux de cette recherche sont : I) explorer le procédé chimique in situ pendant la synthèse des SWCNTs et 2) examiner de manière approfondie l'effet de l'environnement réactif sur les SWCNTs. Les effets du type de catalyseur et son contenu sur le produit fini des SWCNTs, synthétisé par plasma thermique inductif (PTI), ont été étudiés pour remplacer le cobalt (Co) toxique dans la matière première. À cet égard, trois mélanges de catalyseurs différents (c.-à-d. Ni-Y203, Ni-Co-Y203, et Ni-Mo-Y203) ont été utilisés. Les résultats expérimentaux ont montré que le type de catalyseur affecte la qualité des SWCNTs. Une qualité similaire peut être produite lorsque la même quantité de Co est remplacée par le Ni. En outre, des résultats observés dans les travaux expérimentaux ont été explicités par les résultats des calculs thermodynamiques. La therrnogravimétrie (TG) a été utilisée tout au long du travail pour caractériser les échantillons de SWCNTs. La TG a tout d'abord été normalisée par l'étude des effets des trois principaux paramètres instrumentaux (rampe de température, RT, la masse initiale de l'échantillon, MI, et le débit de gaz, D) sur le T, et largeur à mi-hauteur (LMH) obtenu à partir de graphiques TG et TG dérivés de noir de carbone, respectivement. Par conséquent, un plan factoriel à deux niveaux a été prévu. L'analyse statistique a montré que l'effet de RT, MI, et à un degré moindre D est significatif sur la LMH et négligeable sur Tonss. Une méthodologie a ensuite été développée sur la base de la synthèse des SWCNT en utilisant le système PTI, à travers une approche chimique in situ. L'ammoniac (NH3) a été choisi et injecté à contre-courant dans le réacteur PTI à trois débits différents et en utilisant quatre types de buses différentes. La simulation numérique a indiqué un meilleur mélange du NH3 dans le réacteur PTI lorsqu'une buse particulière a été utilisée. Les résultats expérimentaux montrent l'augmentation d'intensité de D-bande dans les spectres Raman d'échantillons SWCNTs lors de l'injection du NH3. Le NH3 pourrait augmenter la teneur en azote du produit fini de SWCNTs jusqu'à 10 fois. L'échantillon des SWCNTs traitée avec 15% vol de NH3 a montré une dispersion accrue dans le diméthylformamide et l'isopropanol. Les nanostructures de carbone en forme d'oignon et plane, ont aussi été observées. Une caractérisation complémentaire sur l'échantillon des SWCNTs traités par NH3 à 15% vol., a indiqué une modification de la surface des nanotubes, où des tubes métalliques ont montré une plus grande réactivité avec NH3 que les semi-conducteurs. Le modèle, y compris le champ d'écoulement thermique du réacteur et la cinétique de décomposition thermique de NH3 a suggéré une modification de surface des SWCNTs en deux étapes dans laquelle les nanotubes réagissent premièrement avec les espèces intermédiaires de H et de NH2. Le NH3 s'adsorbe ensuite chimiquement sur les nanotubes. Le modèle a également suggéré que les espèces intermédiaires comme le NNH et le N2H2 jouent un rôle principalement en conduisant la décomposition du NH3 plutôt que la modification chimique des SWCNTs. // Abstract : Pristine single-walled carbon nanotubes (SWCNTs) are poorly dispersible and insoluble in many solvents and need to be chemically modified prior to their use in many applications. This work is focused on the investigation of the synthesis of chemically modified SWCNTs material through an in situ approach. The main objectives of the presented research are: I) to explore the in situ chemical process during the synthesis of SWCNT and 2) to closely examine the effect of a reactive environment on SWCNTs. Effects of the catalyst type and content on the SWCNTs final product, synthesized by induction thermal plasma (1TP), were studied to replace toxic cobalt (Co) in the feedstock. In this regard, three different catalyst mixtures (i.e. Ni-Y203, Ni-Co-Y203, and Ni-Mo-Y203) were used. Experimental results showed that the catalyst type affects the quality of the SWCNT final product. Similar quality SWCNTs can be produced when the same amount of Co was replaced by Ni. Moreover, the results observed in this experimental work were further explained by thermodynamic calculation results. Thermogravimetry (TG) was used throughout the work to characterize the SWCNTs product. TG was firstly standardized by studying the effects of three main instrumental parameters (temperature ramp, TR, initial mass of the sample, 1M, and gas flow rate, FR) on the Lise, and full-width half maximum (FWHM) obtained from TG and derivative TG graphs of carbon black, respectively. Therefore, a two-level factorial statistical design was performed. The statistical analysis showed that the effect of TR, IM, and to a lower extent, FR, is significant on FWHM and insignificant on T01, 1. A methodology was then developed based upon the SWCNTs synthesis using the 1TP system, through an in situ chemistry approach. Ammonia (NH3) was selected and counter-currently injected into the ITP reactor at three different flow rates and by four different nozzle designs. Numerical simulation indicated a better mixing of NH3 in the ITP reactor when a certain nozzle was used. The experimental results showed the increase of D-band intensity in the Raman spectra of SWCNT samples upon the NH3 injection. NH3 could increase the nitrogen content of the SWCNTs final product up to 10 times. The SWCNTs sample treated with 15 vol% NH3 showed an enhanced dispersibility in Dimethylformamide and Isopropanol. Onion-like and planar carbon nanostructures were also observed. Complementary characterization on the SWCNT samples treated by 15 vol% NH3 indicated the surface modification of nanotubes. Metallic tubes showed a higher reactivity with NH3 than semiconducting ones. The model including the reactor thermo-flow field and NH3 thermal decomposition kinetics suggested a two-step SWCNT surface modification in which nanotubes firstly react with H and NH2 intermediates and later, NH3 chemisorbs on the nanotubes. The model also suggested that the intermediate species, like NNH and N2H2, play a rote primarily in driving the NH3 decomposition rather than the chemical modification of SWCNTs. [symboles non conformes]

La coalescence intervient dans de nombreuses situations physiques, naturelles ou industrielles, de la microphysique des nuages à la stabilité des émulsions ou l’assèchement des pétroles. Dans toutes ces situations, il est crucial de comprendre les mécanismes physiques en jeu, de manière pouvoir influencer la coalescence, la favoriser ou au contraire l’inhiber, selon les besoins. Dans cette thèse, nous étudions la coalescence dans l’air entre deux gouttes attachées et décomposons le processus global en quatre étapes : l’approche avec drainage du film d’air entre les gouttes, le perçage des interfaces, l’ouverture du pont résultant de ce perçage, les oscillations amorties conduisant à l’équilibre de la goutte résultante. Les théories décrivant les étapes 1, 2 et 4 font intervenir des modèles hydrodynamiques continus, se plaçant à une échelle macroscopique. Cependant, à l’articulation entre les deux premières étapes, intervient le perçage des interfaces, processus gouverné par des forces dont la portée correspond à une échelle de quelques dizaines de nanomètres. Une des difficultés les plus importantes dans l’étude de la coalescence est celle de l’intégration des processus ayant lieu à un niveau moléculaire, dans une théorie du continuum dont l’échelle caractéristique est bien supérieure. L’objectif est de faire le lien entre les différentes échelles : y a-t-il des interactions entre les processus se produisant à ces différentes échelles ? Pour répondre à cette question, nous développons trois axes de travail, engageant chacun une échelle caractéristique. L’un est l’étude, au niveau macroscopique du micromètre, de l’ouverture du pont liquide. Grâce à une caméra rapide, plusieurs régimes d’écoulement sont mis en évidence. Les modèles théoriques existants concernent essentiellement le régime visqueux, et aucun modèle complet ne décrit le régime purement inertiel. Nous explorons expérimentalement ce régime et décrivons la forme et la longueur du pont, à l’aide d’ondes capillaires. Nous mettons en évidence l’existence de deux lignes de très forte courbure, que nous appelons singularités, qui naissent sur le lieu de perçage des interfaces et se propagent presque sans déformation de part et d’autre. Ces singularités, conditionnées par la tension superficielle, moteur de la coalescence, façonnent la forme du pont liquide et donc l’écoulement dans ce dernier. Nous proposons un modèle simple d’écoulement inertiel, basé sur la forme du pont liée à ces singularités. Ce modèle permet de mieux comprendre les rôles des forces hydrodynamiques et de la courbure dans l’évolution temporelle de la largeur du pont. Un autre axe est une étude expérimentale par Microscope à Force Atomique, qui permet de décrire les forces responsables de la coalescence à l’échelle nanométrique, les déformations des gouttes intervenant à cette échelle et leur rôle dans la rupture des interfaces. Les mesures de forces entre goutte et flaque, puis entre deux gouttes sont effectuées avec un AFM principalement en mode dynamique de Modulation de Fréquence. Elles permettent de mettre en évidence une distance seuil de déclenchement de l’instabilité hydrodynamique responsable de la coalescence et de mesurer cette distance en fonction des propriétés physiques du liquide et du rayon des gouttes. Un diagramme de coalescence est proposé, qui permet de prévoir la valeur de la distance de déclenchement de la coalescence et le rôle des déformations d’interfaces à l’échelle nanométrique. Enfin, les oscillations du pont liquide, générées par la coalescence, sont étudiées, les modes et fréquences propres sont calculés numériquement par la méthode des éléments finis, puis comparés aux valeurs expérimentales mesurées à partir des films acquis par caméra rapide.

Les méthodes de domaines fictifs permettent de simuler numériquement des écoulements autour de structures complexes et/ou mobiles à l’aide de maillages simples. L’objet solide est alors « immergé » dans un domaine de calcul englobant le fluide et le solide. Dans un premier temps, on étudie une méthode de pénalisation, qui consiste à ajouter un terme dans l’équation de conservation de la quantité de mouvement du fluide afin d’imposer la vitesse du solide. Grâce à des développements asymptotiques, on obtient une estimation de l’erreur induite par cette approche lorsque le solide est en mouvement. Ce procédé est ensuite couplé avec un schéma de projection vectorielle permettant d’imposer la contrainte d’incompressibilité. La convergence du schéma ainsi obtenu, vers les équations de Navier-Stokes, est établie. Dans un second temps,une approche originale capable de traiter des écoulements multiphasiques est développée : la méthode de porosité variable. L’idée principale est de considérer le solide comme un milieu sans masse. La discrétisation des bilans massiques de chaque phase est alors modifiée, de sorte que le volume total occupé par l’ensemble des phases fluides soit égal au volume laissé libre parle solide. Cette méthode est validée numériquement sur un ensemble varié de cas test comprenantd es écoulements monophasiques incompressibles et compressibles ainsi que des écoulements diphasiques. / Fictitious domain methods allow to simulate flows around complex and/or moving bodies with simple meshes. The object is "immersed" in a domain that contains fluid and solid volumes. The penalization method, which consists in adding a term in the momentum balance equation, in order to impose the solid velocity, is studied in a first part. Thanks to asymptotic expansions, the order of the error induced by this method is computed for moving bodies. This approach is then coupled with a Vector Penalty Projection scheme that permits to impose the incompressibility constraint. The convergence of the penalized scheme towards the Navier-Stokes equations is established. In a second part, an original approach, able to treat multiphase flowsis presented: the Time and Space Dependent Porosity method. The key idea is to consider the solid as a medium without mass. The discretization of the mass balance equation is modified,so that the total volume occupied by all fluid phases and the solid is equal to the total volume.This method is numerically validated on a set of various test cases including incompressible or compressible single phase flows and two-phase flows.