On considère le problème d'advection-diffusion stationnaire v(∇u, ∇v)+( β*∇u, v) = (f, v) et non stationnaire d/dt (u(t), v) + v(∇u, ∇v)+( β*∇u, v) = (g(t), v), ainsi que le problème d'advection (β*∇u, v) = (f, v) sur un domaine polygonal borné du plan. Le terme de diffusion est approché par des éléments de Crouzeix Raviart et le terme de convection par une méthode upwind sur des volumes barycentriques finis avec un maillage triangulaire. Pour le problème stationnaire d'advection-diffusion, la L²-stabilité (c'est-à-dire indépendante du coefficient de diffusion v) est démontrée pour la solution du problème approché obtenue par cette méthode d'éléments finis et de volumes finis. Pour cela une condition sur la géométrie doit être satisfaite. Des exemples de maillages sont donnés. Toujours avec cette condition géométrique sur le maillage, une inégalité de stabilité (où la discrétisation en temps n'est pas couplée à une condition sur la finesse du maillage) est obtenue pour le cas non-stationnaire. La discrétisation en temps y est faite par un schéma d'Euler implicite. Une majoration de l'erreur, proportionnelle au pas en temps et à la finesse du maillage, est ensuite proposée et exprimée explicitement en fonction des données du problème. Pour le problème d'advection, une approche utilisant la théorie des graphes est utilisée pour obtenir l'existence et l'unicité de la solution, ainsi que le résultat de stabilité. Comme pour la stabilité du problème d'advection-diffusion, une condition géométrique - qui est équivalente pour les points intérieurs du maillage à celle du problème d'advection-diffusion - est nécessaire.

Ce travail est une contribution à la simulation numérique d'écoulements turbulents autour de corps non profilés. Après avoir précisé les ingrédients numériques et les modèles de turbulence utilisés dans nos simulations, nous présentons une étude sur l'évaluation des effets de la procédure dynamique des modèles de sous-maille dans un modèle VMS-LES et une approche hybride RANS/VMS-LES. Des problèmes d'écoulements autour d'un cylindre seul et en tandem sont considérés. Nous étudions ensuite le comportement de modèles de turbulence hybrides pour la simulation d'écoulements en régime sous-critique autour d'un cylindre circulaire. Le calcul de l'écoulement autour d'un cylindre de section rectangulaire par l'approche VMS-LES est aussi présenté. Enfin, dans une dernière partie, après avoir fait une revue des travaux importants sur les schémas d'avancement en temps multirate, nous proposons une nouvelle approche explicite multirate par agglomération de volumes finis que nous appliquons à des calculs d'écoulements turbulents complexes en utilisant un modèle de turbulence hybride. / This work is a contribution to the numerical simulation of turbulent flows around bluff bodies. After specifying the numerical ingredients and the turbulence models used in our simulations, we present a study on the impact of the dynamic sub-grid scale modeling in VMS-LES model and a RANS/VMS-LES hybrid turbulence approach. Simulations of flows around a cylinder and a tandem are performed. Next, we assess the behaviour of some hybrid turbulence models for the simulation of flows around a circular cylinder in the subcritical regime. The computation of the flow around a rectangular cylinder with the VMS-LES approach is also presented. At last, after a review of some important works on multirate time advancing schemes, we propose a new volume-agglomeration explicit multirate approach that is applied to the computation of complex turbulent flows by a hybrid turbulence model.

Dans cette thèse, nous nous intéressons à l'étude d'un modèle de transport de sédiments en nous plaçant sous deux angles d'approche différents. L'un concerne la modélisation numérique du problème et propose une méthode de résolution numérique basée sur un solveur de Riemann approché pour le système de Saint-Venant-Exner qui reste un des systèmes les plus répandus pour traiter le transport sédimentaire par charriage. Ce dernier repose sur un couplage du modèle hydraulique de Saint-Venant et du modèle morphodynamique d'Exner. Le point essentiel de la méthode proposée se situe au niveau du traitement du couplage de ce système. En effet, il existe deux stratégies; la première consiste à découpler la résolution de la partie fluide de la partie solide et les faire interagir à des instants donnés alors que la seconde considère une résolution couplée du système en mettant à jour conjointement les grandeurs hydrauliques et solides aux mêmes instants. Se posera alors la question du choix de la stratégie de résolution pour laquelle nous apporterons des éléments de réponses en comparant les deux approches. L'autre se concentre sur la mise en place d'une méthodologie pour l'étude des incertitudes liées au même modèle. Pour ce faire, nous proposons une formulation stochastique du système de Saint-Venant-Exner et nous cherchons à caractériser la variabilité des sorties par rapport aux paramètres d'entrée naturellement aléatoires. Cette première étude révélera la nécessité de revenir à un système de Saint-Venant avec un fond bruité pour étudier la sensibilité des grandeurs hydrauliques par rapport aux perturbations topographiques. / In this thesis, we are interested on the study of a sediment transport model through two different approaches. One of them concerns the numerical modelling of the problem and proposes a numerical problem-solving method based on an approximate Riemann solver for the Saint-Venant-Exner system which is one of the most common model to deal with sedimentary bed-load transport. This last one is based on a coupling between the hydraulic model of Saint-Venant and the morphodynamic model of Exner. The key point of the proposed method is the treatment of the coupling issue. Indeed, there exists two strategies; the first one consists on decoupling the resolution of the fluid part from the solid part and making them interact at fixed times whereas the second one considers a coupled approach to solve the system by jointly updating the hydraulic and solid quantities at same times. We then raise the issue of the choice of the strategy for which we suggest answers by comparing both approaches. The other one focuses on the development of a methodology to study the uncertainties related to the model previously mentioned. To this end, we propose a stochastic formulation of the Saint-Venant-Exner system and we look for characterizing the variabilities of the outputs in relation to the naturally random input parameters. This first study reveals the need for a return to the Saint-Venant system with a perturbed bed to understand the sensitivity of the hydraulic quantities on the topographical perturbations.

Depuis de nombreuses années, le givrage a été identifié comme un danger dans le domaine de l’aéronautique.L’accrétion de givre se produit lorsque des gouttelettes d’eau surfondue se déposent sur une surface, enparticulier le bord d’attaque d’une aile ou la lèvre d’entrée d’air moteur, et gèlent après l’impact. Ceci peutensuite engendrer une dégradation des performances aérodynamiques, un dysfonctionnement des sondes ouencore un endommagement du moteur. C’est pourquoi cette problématique est étudiée avec attention. Lesessais en vol et en soufflerie étant longs et coûteux, la simulation numérique de l’accrétion de givre est devenueun outil nécessaire dans le processus de conception et de certification des avions. Cette thèse s’inscrit dans le contexte de la modélisation 3D de l’accrétion de givre, et plus particulièrement des couches limites dynamique et thermique qui se développent autour du corps givré. Les outils numériques devant être rapides et robustes, l’approche proposée dans cette thèse pour le calcul aérodynamique est une méthode couplée Euler/couche limite intégrale. Ainsi, un modèle intégral est développé pour représenterle développement de la couche limite dynamique. La partie thermique est modélisée soit par une méthodesimplifiée basée sur des approches algébriques, soit par une méthode intégrale. Cette modélisation des coucheslimites dynamique et thermique est valable sur paroi lisse ou rugueuse et permet de fournir notamment lecoefficient de frottement et le coefficient d’échange thermique qui sont nécessaires pour un calcul d’accrétion degivre. Les équations intégrales de couche limite, associées à leurs relations de fermeture, sont ensuite résoluespar une méthode Volumes-Finis sur maillage surfacique non structuré, qui est bien adaptée pour les géométriescomplexes. De plus, des traitements numériques spécifiques sont mis en œuvre pour améliorer la précision dela méthode au voisinage du point d’arrêt et pour rendre le code robuste au passage du décollement.Après la validation de la méthode de couche limite, le code est utilisé dans les chaînes de givrage 2D et 3Dde l’ONERA pour des applications d’accrétion de givre. Ceci permet de montrer l’intérêt de la méthode entermes de robustesse et de précision par rapport aux codes de couche limite habituellement utilisés dans lescodes de givrage actuels. / Icing has since long been identified as a serious issue in the aeronautical world. Ice accretion occurs whensupercooled water droplets impinge on a surface, particularly the leading edge of a wing or an engine inlet, andfreeze after the impingement. This can lead to degradation of aerodynamic performances, sensor malfunctionor engine damage. This is why this issue is being carefully studied. The lengthy and costly flight and windtunnel tests have made numerical simulation of ice accretion a necessary tool in the aircraft design andcertification process. The present work deals with the 3D numerical modeling of ice accretion, and more particularly the modeling of the dynamic and thermal boundary layers that develop around an iced body. Since numerical tools must befast and robust, the approach proposed in this thesis for aerodynamic computation is a coupled Euler/integralboundary layer method. Thus, an integral model is developed to represent the development of the dynamicboundary layer. The thermal part is modeled either by a simplified method based on algebraic approaches,or by an integral method. This modeling of the dynamic and thermal boundary layers is valid on smoothor rough walls and provides the friction coefficient and heat exchange coefficient that are necessary for thecalculation of ice accretion. The integral boundary layer equations, associated with their closure relations,are then solved by a Finite-Volume method on unstructured surface mesh, that is well suited for complexgeometries. In addition, specific numerical treatments are implemented to improve the accuracy of the methodin the vicinity of the stagnation point and to make the code robust to separated boundary layers.After validation of the boundary layer method, the code is used in ONERA’s 2D and 3D icing tools foricing applications. This demonstrates the value of the method in terms of robustness and accuracy comparedto the boundary layer codes more commonly used in current icing tools.

Le Disjoncteur Basse-Tension (DBT) est un appareil classique de la distribution électrique depuis plus de cinquante ans. Mais aujourd'hui, avec l'arrivée de produits bas-coût fabriqués dans les pays émergents, les industriels sont soumis à une forte pression pour développer de nouveaux systèmes moins encombrants, utilisant d'autres matériaux, ou incorporant davantage de fonctionnalités. Cette recherche est très compliquée dans la mesure où le DBT est un système hautement multi-physique (mécanique, thermique, physique des matériaux, physique des plasmas, ...). De fait, le développement de nouveaux produits passe par un processus empirique long et coûteux. Cet effort pourrait être réduit par l'utilisation de modèles prédictifs permettant d'arriver plus vite à un système fonctionnel. De nos jours, avec l'augmentation des moyens de résolution numérique, de plus en plus de travaux portent sur la description multi-physique en 3D du DBT et notamment sur la chambre de coupure ou l'arc électrique est amorcé, se déplace et doit être éteint, l'objet de nôtre étude. Le travail de cette thèse se divise en deux axes complémentaires : le développement d'un modèle fluide 3D en méthode des volumes finis simulant l'arc électrique et son déplacement dans la chambre de coupure; ainsi que la mise en place d'un dispositif expérimental permettant d'analyser le phénomène physique en œuvre. Pour ces deux points la problématique est abordée dans une configuration simplifiée de DBT où l'arc se déplace entre deux rails parallèles dans une chambre parallélépipédique. Basé sur le savoir-faire du groupe AEPPT, un modèle numérique est établi pour simuler le plasma thermique. Les particularités de ce modèle, du fait de l'application, sont la nécessité d'une résolution précise du champ magnétique en utilisant le calcul de Biot & Savart pour les conditions limites ainsi que l'utilisation de méthodes permettant le déplacement et la commutation de l'arc. La validation de ce modèle se fera à géométrie similaire par confrontation avec l'expérience. En s'inspirant de précédents travaux nous avons réalisé une maquette expérimentale composée d'un réacteur faisant office de chambre de coupure et d'un mécanisme permettant l'amorçage de l'arc dans le réacteur par ouverture rotative du contact à vitesse contrôlée. D'autres paramètres modifiables sont la taille du réacteur ainsi que les matériaux qui le constituent. Les diagnostiques disponibles en plus de la mesure de courant et de tension sont l'imagerie rapide et la mesure de pression en différents points de la chambre de coupure. Notre expérience est utile pour la réalisation d'études paramétriques en découplant facilement les paramètres. En outre, par la mise en évidence des phénomènes prépondérants, notre maquette aide à la mise en place du modèle en plus de permettre sa validation expérimentale. Cette thèse est donc une étape cruciale vers la mise en place d'un modèle prédictif. / Low-Voltage Circuit Breakers (LVCB) are classical apparatuses of electrical distribution since more than fifty years. But nowadays, with the outbreak of low-cost products from the developing countries, industry is under a strong stress in order to improve their devices by making them more compact, using different materials or to implement new functionalities. This research is harsh since LVCB are highly multiphysics systems (mechanics, thermal properties, materials, plasma physics, ...). Therefore, developing new products goes through a long and expensive empirical process. Those efforts could be reduced by using predictive models allowing to get faster to a functional device. With the improvements of the numerical solution capacity, there are more and more works toward the 3D multiphysical description of the LVCB, especially on the extinction chamber where the electrical arc is ignited, moved and must be quenched. This is the subject of our work. The study described here is divided in two complementary parts: development of a 3D fluid model with finite volume method simulating the electrical arc and its movement inside the arc extinction chamber; and the set-up of experimental means to analyse this physical phenomenon. For those two points, we use a simplified LVCB configuration with an arc moving between two parallel rails inside a rectangular box chamber. Based on AEPPT's know-how, a numerical model is established to simulate thermal plasma. Particularities on this model, due to the application, are the resolution of Biot&Savart law to calculate precisely the magnetic field for the boundary condition and the development of methods to model the arc roots movement and commutation of the arc from the moving contact to the rail. Validation of this model will be done with the same geometry by confrontation with the experiment. Inspired by precedent works we designed a test apparatus with a reactor representing the extinction chamber of the LVCB and an opening mechanism allowing arc ignition by contact opening at a specified speed. Other parameters such as size of the chamber and materials can be modified. Measurements will include high speed imaging and pressure acquisition in several points of the reactor in addition to the classical current and voltage measurements. This experiment is useful for parametric studies with its easy uncoupling of the parameters. Moreover, by highlighting the dominating phenomena for arc movement, this set-up helps in the build-up of the model over and above the experimental validation.

Les capacités installées d’énergie éolienne continuent à croître rapidement et prennent une place de plus en plus significative dans le monde. Au fur et à mesure, les études menées sur la conception, la sureté de fonctionnement et la supervision de la chaîne éolienne ont pris progressivement de l’importance. Deux axes de recherche ont été privilégiés dans cette thèse. Le premier concerne la continuité de service d'une éolienne connectée au réseau en présence de défaut de court-circuit entre spires dans une phase du stator de la génératrice asynchrone à cage d'écureuil. L'analyse du défaut ainsi que son impact sur le système éolien et notamment sur la qualité de la puissance produite souligne l'intérêt de développement d'un algorithme de détection et d'isolation rapide, dédié par la suite à la reconfiguration de la commande. Ainsi, une commande tolérante au défaut (CTD) a été conçue de manière à éviter l'arrêt de la production, compenser l'impact de défaut et garder des performances acceptables de la qualité d'énergie produite. Le travail effectué s'est articulé sur les observateurs à mode glissant (OMG), communément connus comme outil puissant pour la supervision et la commande à la fois. Le deuxième axe porte sur la sécurité structurale et la stabilité du système éolien sous contraintes vibratoires. Les travaux se répartissent en deux parties complémentaires : L'établissement d'un modèle numérique tridimensionnel (3-D) sous un logiciel d’analyse par éléments finis (ANSYS) et la réalisation des essais vibratoires sous différentes excitations au sein d'une plateforme vibratoire (TREVISE). Dans ce cadre, un modèle numérique (3-D) d'une éolienne à axe horizontal couplée à un mât et une fondation adéquats a été développé en utilisant la méthode de volumes finis (FVM) afin d'appréhender son comportement vibratoire. Les essais vibratoires expérimentaux valident le modèle numérique et permettent l’identification de la réponse dynamique de la structure d'une manière fine. De plus, nous avons élaboré un modèle expérimental de la tenue de l’éolienne aux contraintes vibratoires de formes aléatoire, sinusoïdale et impulsionnelle. / The wind energy capacity carries on growing quickly and taking an increasingly significant place in the world. Progressively, research studies dealing with designing and supervising wind turbines have become more important. Two areas of research were developed in this thesis. The first one concerns the continuity of service of a wind turbine connected to the grid while an inter-turn short-circuit fault is present in the stator phase of the induction squirrel cage generator. The analysis of the fault as well as its impact on the wind turbine system and mainly on the quality of the produced power highlights the interest of development of a fast detection and isolation algorithm, dedicated to the reconfiguration of the control law. Hence, a fault tolerant control scheme has been established in order to avoid stopping production, compensate the fault impact and maintain acceptable performances of the quality of the produced energy. The carried out work was based on sliding mode observers, commonly known as robust tools for monitoring and controlling at the same time. The second axis concerns the structural modeling and stability checking of the wind system under vibratory stresses. The work is divided into two complementary parts: The establishment of a three-dimensional (3-D) numerical model using a finite element analysis software (ANSYS) and the realization of vibratory tests under different excitations within the platform (TREVISE). In this framework, a numerical (3-D) model of a horizontal axis wind turbine coupled to a suitable tower and foundation was developed basing on the finite volume method (FVM) in order to analyze its vibratory behavior. The experimental vibratory tests validate the numerical model and allow the identification of the dynamic response of the structure in a precise way. In addition, we have developed an experimental model of the behavior of the wind turbine under vibratory stresses of random, sinusoidal and impulse shapes.

Cette thèse est principalement dédiée à la simulation des écoulements massivement décollés dans le domaine spatial. Nous avons restreint notre étude aux écoulements d'arrière-corps, pour lesquels ces décollements sont imposés par des changements brutaux de la géométrie. Dans le domaine spatial, le caractère fortement compressible des écoulements rencontrés impose l'utilisation de schémas numériques robustes. D'un autre coté, la simulation fine de la turbulence impose des schémas d'ordre élevé et peu dissipatifs. Ces deux spécifications, apparemment contradictoires, doivent pourtant coexister au sein d'une même simulation. Les modèles de turbulence ainsi que les schémas de discrétisation sont indissociables et leur couplage doit impérativement être considéré. Les schémas numériques doivent garder leur précision formelle dans des géométries complexes et des maillages très irréguliers imposés par le contexte industriel. Cette étude analyse le schéma de discrétisation utilisé dans le code de calcul FLUSEPA développé par Airbus Defence & Space. Ce schéma est robuste et précis pour des écoulements avec chocs et il présente une faible sensibilité au maillage (l'ordre 3 étant conservé même sur des maillages fortement perturbés). Malheureusement, le schéma possède une trop faible résolvabilité liée à un niveau de dissipation trop élevé pour envisager des simulations hybrides RANS/LES. Pour pallier à cet inconvénient, nous nous sommes penchés vers une solution basée sur un recentrage conditionnel et local : dans les zones dominées par des structures tourbillonnaires, une fonction analytique assure un recentrage local lorsque la stabilité numérique le permet. Cette condition de stabilité assure le couplage entre le schéma et le modèle. De cette manière, les viscosités laminaire et tourbillonnaire sont les seules à jouer un rôle dans les régions dominées par la vorticité et servent aussi à stabiliser le schéma numérique. Cette étude présente de plus une comparaison qualitative et quantitative de plusieurs modèles hybrides RANS/LES, à égalité de maillage et de schéma utilisés Pour cela, un certain nombre d'améliorations (notamment de leur capacité à résoudre les instabilités de Kelvin-Helmohlotz sans retard), proposées dans la littérature ou bien introduites dans cette thèse, sont prises en compte. Les applications numériques étudiées concernent des géométries allant de la marche descendante au lanceur spatial complet à échelle réduite. / This thesis is mainly dedicated to the simulation of massively separated flows in the space domain. We restricted our study to afterbody flows, where the separation is imposed by abrupt geometry changes. In the space domain, highly compressible flows require the use of robust numerical schemes. On the other hand, the simulation of turbulence imposes high-order low dissipative numerical schemes. These two specifications, apparently contradictory, must coexist within the same simulation. The coupling between turbulence models and discretization schemes is of the utmost importance and must be considered. Numerical schemes should keep their formal accuracy on complex geometries and on very irregular meshes imposed by the industrial context. In this research, we analyze the discretization scheme implemented in the FLUSEPA solver, developed by Airbus Defence & Space. Such a scheme is robust and accurate for flows with shocks and exhibits a low sensitivity to the grid (the third order of accuracy being ensured, even on highly irregular grids). Unfortunately, the scheme possesses a too low resolvability related to a too high numerical dissipation for RANS/LES simulations. To circumvent this problem, we considered a conditional and local re-centering strategy: in regions dominated by vortical structures, an analytic function provides local re-centering when a numerical stability condition is satisfied. This stability condition ensures the coupling between the numerical scheme and the model. In this way, only the turbulent and the laminar viscosities play a role in regions dominated by vorticity, and also allow to stabilize the numerical scheme. This study provides also a qualitative and quantitative assessment of several hybrid RANS/LES models, using the same grids and discretization scheme. For this purpose some recent improvements (improving their ability to trigger the Kelvin-Helmohlotz instabilities without delay), proposed in the litterature or suggested in this work, are taken into account. Numerical applications include geometrical configurations ranging from a backward facing step to realistic launcher configurations.

Le calcul d'écoulements dans les générateurs de vapeur des réacteurs à eau pressurisée est un problème complexe, faisant intervenir différents régimes d'écoulement et plusieurs échelles de temps et d'espace. Un scénario accidentel peut être caractérisé par des variations très rapides pour un nombre de Mach de l'ordre de l'unité. A l'inverse en régime nominal l'écoulement peut être stationnaire, à bas nombre de Mach. De plus quelque soit le régime considéré, la complexité de la géométrie d'un générateur de vapeur conduit à modéliser le faisceau de tubes par un milieu poreux, d'où le problème de couplage à l'interface avec le milieu fluide.Un schéma de correction de pression tout-Mach en volumes finis colocalisés a été introduit pour les équations d'Euler et de Navier-Stokes. L'existence d'une solution discrète, la consistance du schéma au sens de Lax et la positivité de l'énergie interne ont été démontrées. Le schéma a été ensuite étendu aux modèles diphasiques homogènes du code GENEPI développé au CEA. Enfin un algorithme Multigrille-AMR a été adaptée pour permettre de mettre en oeuvre notre schéma sur des maillages adaptatifs.Concernant la seconde problématique, une extension de la loi de Beavers-Joseph a été proposée pour le régime convectif. En introduisant un saut d'énergie cinétique à l'interface, on retrouve une loi de type Beavers-Joseph mais avec un coefficient de glissement non-linéaire, qui dépend de la vitesse fluide à l'interface et de la vitesse Darcy. La validité de cette nouvelle condition d'interface a été évaluée en réalisant des calculs de simulation numérique directe à différents nombres de Reynolds. / The numerical simulation of steam generators of pressurized water reactors is a complex problem, involving different flow regimes and a wide range of length and time scales. An accidental scenario may be associated with very fast variations of the flow with an important Mach number. In contrast in the nominal regime the flow may be stationary, at low Mach number. Moreover whatever the regime under consideration, the array of U-tubes is modelled by a porous medium in order to avoid taking into account the complex geometry of the steam generator, which entails the issue of the coupling conditions at the interface with the free-fluid.We propose a new pressure-correction scheme for cell-centered finite volumes for solving the compressible Navier-Stokes and Euler equations at all Mach number. The existence of a discrete solution, the consistency of the scheme in the Lax sense and the positivity of the internal energy were proved. Then the scheme was extended to the homogeneous two-phase flow models of the GENEPI code developed at CEA. Lastly a multigrid-AMR algorithm was adapted for using our pressure-correction scheme on adaptive grids.Regarding the second issue addressed in this work, an extension to the Beavers-Joseph law was proposed for the convective regime. By introducing a jump in the kinetic energy at the interface, we recover an interface condition close to the Beavers-Joseph law but with a non-linear slip coefficient, which depends on the free-fluid velocity at the interface and on the Darcy velocity. The validity of this new transmission condition was assessed with direct numerical simulations at different Reynolds numbers.

La première contribution de cette thèse est l'étude de la modélisation ainsi que de la discrétisation des écoulements multiphasique en géométrie microfluidique. Nous proposons un nouveau schéma d'ordre deux de discrétisation basé sur les méthodes mixtes volumes finis/éléments finis pour le système de Stokes bi-fluides avec tension de surface.Nous présentons alors des comparaison de la précision de ce nouveau schéma avec la discrétisation MAC, en 2D et en 3D-axi. La seconde contribution est relative à l'étude de schémas numériques en temps pour les fluides viscoélastiques. Nous présentons les limites actuelles des modélisations dans ce domaine en étudiant le cas des écoulements de micelles géantes, polymères ayant la capacité de se réorganiser spatialement en fonction du taux de cisaillement. Nous montrons qu'une condition de stabilité liée au ratio de viscosité du polymère et du solvant en temps - très restrictive - existe. Un nouveau schéma est alors proposé pour contourner cette limitation et des études stabilité sont menés pour démontrer nos résultats. / The first contribution of this thesis is the analysis and discretization of multiphase flow in a microfluidic framework.We propose a new scheme wich is second order accurate, based on mixed finite volume / finite element method for the two phases Stokes system with surface tension.We then present the comparison of the accuracy of this new scheme with the MAC discretization in 2D and 3D axi.The second contribution is related to the study of numerical schemes in time for viscoelastic fluids.We present the current limitations in this area by studying the case of flows of wormlike micelles, polymers having the ability to reorganize themselves according to the shear rate.We show that a condition of stability related to the ratio of viscosity of the polymer and solvent in time - very restrictive - exists.A new scheme is then proposed to overcome this limitation and stability studies are conducted to demonstrate our results.

Dans le but d’améliorer le comportement thermo-mécanique des aubes de turbine présentes dans les moteurs d’avion développés par SAFRAN, il est nécessaire de mettre au point un nouveau matériau composite. Le procédé de fabrication de ce matériau est complexe et requiert une densification par voie liquide divisée en deux étapes. Cette thèse s’intéresse à la modélisation numérique de la première étape appelée Slurry Cast/APS.Celle-ci consiste en l’injection et le confinement, dans la préforme fibreuse, de particules préalablement mises en suspension. Pour cela, nous avons développé à l’échelle des fibres,un modèle qui utilise les équations de Navier-Stokes incompressibles et monophasiques ,l’équation de Phillips [Phillips et al., 1992] et une loi rhéologique [Krieger, 1972]. Après validation des résultats numériques par comparaison avec des résultats expérimentaux [Hampton et al., 1997] et théoriques [Belfort et al., 1994], le modèle est utilisé pour simuler l’écoulement autour de géométries de tissage proches du matériau étudié. / In order to improve thermo-mechanical behavior of tubine blades in SAFRAN engines plane, a new composite material is necessary. The manufacturing process to obtain this composite is intricate and requires a two steps fluid densification process. This thesis focuses on numerical simulation of the first one called Slurry Cats/APS. In this step, suspended particles are introducted and captured in the reinforcement. For that purpose,we carry out a model at fiber scale, using Navier-Stokes equations in incompressible and monophasic formulation, Phillips equations [Phillips et al., 1992] and a rheological law [Krieger, 1972]. After validation step consisting in a comparison of computational results with experiments [Hampton et al., 1997] and theorical law [Belfort et al., 1994], this model has been used to simulate flow around geometries similar to those encountered in our composite material.