Depuis plusieurs années à nos jours, les paliers à gaz sont utilisés avec succès sur une large gamme de turbo-machines. Certains de ces systèmes sont utilisés dans des environnements de contrôle de l’environnement par gaz réfrigérant. Dans ce travail, nous présentons un modèle théorique et numérique qui tient compte de la transition du lubrifiant vapeur / liquide, la transition de l’écoulement laminaire / turbulent et les variations 3D de la viscosité et de la température dans le fluide et les solides pour les deux situations statiques et dynamiques. Ce modèle comporte : la résolution de l'équation de Reynolds généralisée pour les fluides compressibles à viscosité variable en 3D, la description des effets de la turbulence en utilisant l'approche phénoménologique de Elrod, en utilisant un champ de viscosité turbulente 3D, la résolution de l’équation d’état non linéaire du lubrifiant, capable de décrire la transition vapeur / liquide et une approche thermique local afin d'obtenir une estimation 3D de la température du fluide, grâce à l'équation d'énergie pour film mince. La prise en compte également des effets thermiques dans les solides. Dans cette étude, nous avons montré l'importance d'une description précise des paramètres du film fluide, dont les variations influencent largement le comportement du palier. Parmi les principales théories, il y a: lubrifiant compressible, avec un comportement non-linéaire près de la transition vapeur / liquide, la transition vapeur / liquide et le calcul des paramètres équivalents du mélange, un écoulement turbulent du fluide pour le palier GFB à grande vitesse en utilisant un modèle 3D de la viscosité turbulente, un comportement 3D pour la viscosité, en particulier les variations inter-films (dépendant de la température), et un comportement 3D pour la température, en particulier dans le sens transversal du film afin d'être compatible avec la viscosité, mais également dans la direction axiale afin de tenir en compte du gradient de température potentiel qui modifie considérablement le profil 3D de la température du palier. Ces deux comportements statiques et dynamiques du palier compliant GFB sont analysés. / For years now, gas bearings are successfully used over a large panel of turbo-machineries. Some of these systems are bound to be run in controlled environments such as refrigerating gas. In this work we present a theoretical and numerical model which takes into account the vapor/liquid lubricant transition, the laminar/turbulent flow transition and both temperature and viscosity 3D variations in the fluid and the solids for both static and dynamic situations. This model involves: the resolution of the generalized Reynolds equation for compressible fluids with 3D variable viscosity, the description of the turbulence effects by the phenomenological approach of Elrod, using a 3D eddy viscosity field, the resolution of a non-linear equation of state for the lubricant, able to describe the vapor/liquid transition and a local thermal approach to obtain a 3D estimation of the fluid temperature, thanks to the thin-film energy equation. The thermal effects in solids are also taken into account. In this study, we showed the importance of an accurate description of the film parameters, which variations largely influence the bearing behaviour. Among the principal theories, there are: compressible lubricant, with an appropriate non-linear behaviour when close to the vapor/liquid transition, vapor/liquid transition and calculation of the mixture equivalent parameters, turbulent flow for high-speed GFBs with a 3D eddy viscosity mode, a 3D behaviour for viscosity, particularly the cross-film variations, (temperature dependent)and a 3D behaviour for temperature, particularly in cross-film direction in order to be consistent with viscosity, but also in the axial direction in order to account for potential temperature gradient which considerably modifies the bearing 3D temperature profile. Both static and dynamic behaviours of GFBs are analysed.

Dans un système de machines rotatives comme un moteur à turbine, les paliers d'éléments roulants à grande vitesse jouent un rôle important dans le support de l'arbre ou du rotor rotatif, et ont besoin d'une lubrification pour assurer leur fonction. Sauf qu'une petite quantité d'huile est nécessaire pour former le film lubrifiant élastohydrodynamique dans la zone de contact, la plus grande partie du lubrifiant reste en suspension dans l'air, formant un mélange huile/air. Ce phénomène entraîne des pertes hydrauliques parasitaires excessives lorsque les éléments roulants se translatent et tournent dans l'environnement fluide, ce qui peut constituer une partie relativement importante de la perte de puissance totale du roulement, appelée traînée d'enroulement et pertes de barattage. Pour une vitesse de rotation jusqu'à 3× 106 Ndm, la contribution de la traînée/dérive au total peut atteindre 50%. Cependant, jusqu'à présent, il existe peu d'approches utilisées directement pour l'estimation des pertes par traînage, qui ne pouvait fournir qu'une approximation plutôt grossière. Dans cette thèse, la méthode CFD est utilisée pour analyser d'abord l'écoulement autour d'un cylindre circulaire de longueur finie avec deux extrémités libres dans un espace ouvert. Ensuite, le modèle est remplacé par plusieurs cylindres circulaires en ligne pris en sandwich par deux parois plates, ce qui représente une approche simplifiée. Le fluide est ici considéré comme incompressible, représentant un fluide monophasé équivalent pour l'écoulement diphasique huile/air à l'intérieur de la cavité de palier avec des propriétés de fluide spécifiées. Les résultats indiquent que l'écoulement autour de l'élément de rouleau de longueur finie est perturbé par ses deux extrémités libres, les anneaux environnants, la cage et d'autres éléments roulants. Il est proposé une relation entre le coefficient de traînée et le nombre de Reynolds approprié pour un cylindre circulaire dans les roulements à rouleaux (1<L/D<6), ainsi qu'une formulation pour la prévision des pertes de barattage. L'écoulement diphasique huile/air à l'intérieur de la cavité de palier avec lubrification sous la course est également étudié dans ce travail. Le volume couplé de niveau de fluide (CLS-VOF) est utilisé pour démontrer la distribution du lubrifiant le long de la circonférence du palier. L'effet de divers facteurs est étudié, par ex. la vitesse d'injection d'huile, le diamètre de la buse, les propriétés de l'huile et l'angle d'injection de l'huile. La vitesse de rotation de tous les composants du palier est étudiée en particulier pour quantifier leur influence sur la fraction du volume d'huile à l'intérieur de la cavité du palier. Les résultats démontrent que non seulement la vitesse de rotation relative de l'anneau interne, mais la vitesse de la cage elle-même pourrait changer la distribution d'huile. / In a rotating machinery system like turbine engine, high speed rolling element bearings play an important role in supporting the rotating shaft or rotor, and need lubrication to insure their function. Except a small quantity of oil is needed to form the elastohydrodynamic lubricant film in the contact zone, most of lubricant remains in suspension in air, forming an oil/air mixture. This phenomenon leads to excessive parasitic hydraulic losses when rolling elements translate and rotate into the fluid environment, which may constitute a relatively large portion of the bearing's total power loss, named windage drag and churning losses. For high speed applications, i.e. for rotational speed up to 3× 10^6 Ndm, the contribution of drag/windage loss to the total one may reach up to 50%. However, so far there are few approaches used directly for drag and churning losses estimation, which could only provide a rather gross approximation. In this thesis, the Computational Fluid Dynamics (CFD) method is employed to analyze first the flow around one finite-length circular cylinder with two free ends in an open space. Then the model is changed to several in-line circular cylinders sandwiched by two flat walls, which represents a simplified approach. The fluid here is regarded as incompressible, representing an equivalent one-phase fluid for the oil/air two-phase flow inside the bearing cavity with specified fluid properties. The results indicate that the flow around the finite length roller element is perturbed by its two free ends, the surrounding rings, the cage and other rolling elements. A relationship between the drag coefficient and the Reynolds number suitable for circular cylinder in roller bearings (1<L/D<6) is proposed, as well as a formulation for churning losses prediction. The oil/air two phase flow inside the bearing cavity with under-race lubrication is also studied in this work. The coupled level-set volume of fluid (CLS-VOF) method is employed to demonstrate the lubricant distribution along the bearing circumference. The effect of various factors is studied, e.g. the oil injection velocity, the nozzle diameter, the oil properties, and the oil injecting angle. Rotational speed of all the bearing components are studied particularly to quantify their influence to the oil volume fraction inside the bearing cavity. The results demonstrate that not only the inner-ring relative rotational speed, but the cage speed itself could change the oil distribution. The results can be used for the precise lubrication design to optimate the oil distribution inside the bearing.

Cette thèse propose une modélisation numérique 2D des écoulements en zone de jet de rive avec un code Volume-Of-Fluid. Dans un premier temps, le détail de la structure interne de l’écoulement de jet de rive généré par l’effondrement d’un front d’onde turbulent sur une plage imperméable lisse est étudié. Le modèle numérique est ensuite étendu à la simulation des écoulements en milieu poreux internes à la plage. L’utilisation d’une unique équation de quantité de mouvement (VARANS) et de la méthode 1-fluide, permet de résoudre simultanément les écoulements de l’eau et de l’air à la surface et dans une plage perméable. Ce modèle a été confronté à une série de cas tests analytiques et à de récentes mesures expérimentales. Les résultats numériques montrent l’aptitude du modèle VOF-VARANS à reproduire les écoulements en zone de jet de rive sur une plage imperméable fixe. / A 2D numerical modeling of flows in the swash zone is proposed using a Volume-Of- Fluid code. The detailed flow structure of a bore-driven swash event over an impermeable beach is first studied. The numerical model is then developed to account for porous media flow within the beach. The unique VARANS momentum equation and 1-fluid method used allow to solve simultaneously both surface and subsurface flows of air and water phases in the swash zone. This model is validated against a series of analytical tests cases and confronted to recent experimental measurements. The numerical results highlight the ability of the VOF-VARANS model to reproduce swash flows over and within a permeable beach.

Les boucles diphasiques à pompage thermo-capillaire de type LHP (pour Loop Heat Pipe, en anglais), dont le fonctionnement s’apparente à celui d’un caloduc, permettent un transfert de chaleur particulièrement efficace et entièrement passif entre une source chaude et une source froide. Ce transfert s’effectue au moyen d’un fluide diphasique, mû grâce à la force motrice capillaire générée par un matériau poreux contenu dans l’évaporateur/réservoir de la LHP. Outre son rôle de barrière hydraulique entre les phases liquide et vapeur, ce matériau doit assurer une fonction de barrière thermique afin de favoriser l’évaporation du liquide. L’aptitude du matériau à remplir ses fonctions dépend étroitement de sa microstructure, elle-même liée à la méthode de fabrication. Dès lors, l’objectif de la thèse est d’associer la science des matériaux à celle de la thermique, pour améliorer les procédures de fabrication de structures capillaires existantes ou tester de nouvelles méthodes, et aboutir à des structures dont les caractéristiques sont en adéquation avec celles qui sont recherchées. / The capillary pumped loops (CPL) or loop heat pipes (LHP), whom the operating principle is similar to classic heat pipes, enable an efficient heat transfer between a hot source and a cold source without additional energy sources. Indeed, a porous structure provides a capillary force that enables a two-phase fluid to circulate around the loop, transferring the heat from the evaporator to the condenser. The porous structure acts as a hydraulic barrier between the two phases and as a thermal barrier enabling the liquid evaporation. The ability of the capillary structure to fulfil its mission depends on its microstructure, and thus on the manufacturing process. Therefore, the objective of the present thesis is to join the thermal sciences with the material sciences in order to improve the existing manufacturing procedure or even to test new ones. It aims at obtaining capillary structures corresponding to heat transfer applications.

Cette thèse s’inscrit dans le cadre du projet FUI THERMOFLUIDE-RT impliquant des Grands Groupes (Zodiac DS, Safran Hispano, MBDA), des PME (Atmostat, ADR, ControlSys) et cinq laboratoires (CRIStAL, LML Arts et Métiers Paris Tech, LEGI Grenoble, LMT ENS Cachan, CEA-Liten Grenoble). Le but est d’étudier un nouveau système de refroidissement de l’électronique. La technologie retenue est celle d’une boucle fluide diphasique à pompage mécanique (BFDPM). La thèse traite la modélisation dynamique et la validation expérimentale des composants de la boucle. Ceci permet de prévoir l’efficacité du système à partir de ses paramètres d’entrée, d’analyser les problèmes de régimes transitoires, et de proposer un outil de dimensionnement. La méthodologie bond graph est retenue à cause du caractère multi-physique des composants. D’abord, la problématique de base et le contexte sont présentés. Ceci permet d’introduire la solution retenue, celle des BFDPM. Les objectifs de la thèse sont décrits. Ensuite, une description du banc expérimental développé au cours de cette thèse est proposée. Les trois chapitres suivant sont consacrés à l’étude théorique et expérimentale des équipements de la boucle. Chacun de ces chapitres commence par l’état de l’art sur les travaux de modélisation et les corrélations des coefficients d’échange et des pertes de charge. Une seconde partie décrit les phénomènes et les équations. Une troisième partie est réservée à la validation des modèles. Un dernier chapitre récapitule les travaux de couplage des modèles dynamiques validés séparément. En conclusion, un récapitulatif des contributions est effectué. Des perspectives à court et moyen terme sont proposées / This thesis is part of the collaborative project FUI THERMOFLUIDE-RT involving major groups (Zodiac DS, Safran Hispano, and MBDA), SMEs (Atmostat Alcen, ADR, AER, ControlSys) and five laboratories (CRIStAL Ecole Centrale de Lille, LML Arts et Métiers Paris Tech, LEGI Grenoble, LMT ENS Cachan, CEA-Liten Grenoble). The main purpose is to study a new electronic cooling system. The technology chosen consists of a two-phase fluid loop mechanically pumped (TPLMP). The thesis deals with the dynamic modeling and experimental validation of the cooling components. The developed dynamic model allows to predict the efficiency of the cooling loop, to conduct the study of transitional regimes, and provides an original tool dedicated to design the loop components. The bond graph methodology is adopted because of the multi physics character of the studied components. First, the basic issues and the industrial context are presented. This allows to introduce the chosen solution (TPLMP). The objectives of the thesis are described. Then, a description of the rig test is proposed. The following three chapters are devoted to a theoretical and experimental study of the loop equipment. Each chapter begins with a state of the art on modeling and correlations of the heat exchange coefficients and losses. A second part of the chapter describes phenomena and equations. A third part is dedicated to the experimental validation. A final chapter presents the coupling works of dynamic models validated separately. Finally, a summary of all contributions is made. Prospects for future developments in short and medium term are proposed.

France Télécom possède des armoires de télécommunication dont la puissance est limitée à cause de la dissipation thermique des équipements actifs qui entraîne une augmentation de leur température interne. La puissance des équipements limite le nombre de clients qu'il est possible de connecter aux services des réseaux à hauts débits. En plus de cette contrainte, les armoires sont soumises à des effets liés au climat (ensoleillement) qui peuvent être sévères et difficiles à maîtriser. Ceci nécessite l’intégration de systèmes de refroidissement permettant de maintenir la température des composants en dessous de la limite imposée (55°C). C’est dans cet objectif que ce travail de thèse a été mené au sein du laboratoire FEMTO-ST en collaboration avec le service R&D de France Télécom à Lannion. Le refroidissement par changement de phase est favorisé pour maintenir la température de fonctionnement du système stable et pour être utilisé dans les systèmes à haute densité de puissance. Les boucles diphasiques sont des systèmes de refroidissement pour le contrôle thermique et fonctionnent passivement sans pompage mécanique du fluide caloporteur. Après une étude bibliographique sur les boucles de refroidissement diphasiques et leurs applications, on a constaté que les boucles thermosiphons sont particulièrement adaptées aux applications où le faible coût, l'efficacité énergétique et la fiabilité d’entretien sont souhaités. Cette étude a été conduite en suivant un cahier de charge proposé par France Télécom qui consiste à : (i) développer un modèle numérique permettant de modéliser les transferts échangés entre l’armoire de télécommunication et le milieu ambiant, (ii) mener une étude expérimentale en vue de concevoir une boucle thermosiphon pour le refroidissement d’armoires de télécommunication.Le mémoire de cette thèse montre la limitation des systèmes de refroidissement classiques utilisant des écoulements d’air en convection forcée ou autre fluides sans changement de phase. Un modèle numérique est développé afin de permettre la prédiction des températures à l’entrée des boitiers chauffants pour différentes conditions climatiques. Le choix est porté sur l’utilisation d’une modélisation par réseau nodal. La modélisation est effectuée en tridimensionnel et en régime transitoire. Nous avons également modélisé le rayonnement solaire auquel est soumise l’armoire de télécommunication. Le modèle développé a été validé en effectuant une comparaison entre les résultats issus de la modélisation et ceux obtenus à partir des expériences menées au laboratoire et à la plateforme CLIMA chez France Télécom. Les essais sont effectués en régime transitoire en imposant une puissance électrique et en faisant varier la température ambiante ou la densité de flux thermique solaire. L’ensemble des résultats obtenus ont permis de constituer une base de données. Le deuxième objectif fixé dans le cadre de ce travail de thèse est la conception d’un système de refroidissement sous forme d’une boucle thermosiphon. La contrainte principale qui a guidée cette conception était le fait que la boucle doit refroidir l’armoire et assurer une température d’air à l’entrée des équipements inférieure à la limite imposée par la norme ETSI. Ceci nous a mené à concevoir un prototype de boucle thermosiphon dont la puissance thermique qu’il doit dissiper est imposée. On a montré que ce prototype permet de dissiper des puissances thermiques allant jusqu’à 470 W en utilisant une petite charge de npentane. Nous avons effectué des essais sur le refroidissement du prototype d’armoire de télécommunication en utilisant la boucle thermosiphon légèrement modifiée. On montre que les performances thermiques obtenues en utilisant un mode de refroidissement en boucle thermosiphon sont meilleures. Les boucles thermosiphons semblent intéressantes pour un refroidissement passif de matériels déployés dans un réseau de télécommunication... / France Telecom owns telecommunication cabinets whose power is limited because of the heat dissipation of active devices which leads to increased internal temperature. Power equipment limits the number of clients that can connect to networks services with high data rates. In addition to this constraint, the cabinets are subject to climate-related impacts (sunlight) that can be severe and difficult to master. This requires the integration of cooling systems to maintain the temperature of components below the limit (55 ° C). It is with this aim that this work was conducted in the laboratory Femto-ST in collaboration with the R & D department of France Telecom in Lannion.Cooling the phase change is promoted to maintain the operating temperature of the stable and system for use in systems with high power density. The loops are two-phase cooling systems for thermal control and operate passively without mechanical pumping of the coolant.After a literature review on two-phase cooling loops and their applications, it was found that the thermosyphon loops are particularly suitable for applications where low cost, energy efficiency and reliability maintenance are desired. This study was conducted by following a set of specifications proposed by France Telecom which involves: (i) develop a numerical model to model transfers exchanged between the cabinet and the telecommunications environment, (ii) conduct an experimental study to design a thermosyphon loop for cooling telecommunication cabinets.The memory of this thesis shows the limitation of conventional cooling systems using air flow forced convection or other fluids without phase change. A numerical model is developed to enable the prediction of temperatures at the inlet of heated enclosures for different climatic conditions. The choice is focused on the use of a nodal network modeling. The modeling is done by three-dimensional and transient. We also modeled the solar radiation, which applies to the telecommunications closet. The developed model was validated by comparison between the results of modeling and those obtained from experiments in the laboratory and platform CLIMA at France Telecom. The tests are performed by imposing transient electrical power and varying the temperature or heat flux density solar. All the results obtained allowed to establish a database.The second goal as part of this thesis is the design of a cooling system as a thermosyphon loop. The main constraint has guided this design was that the loop needs to cool the cabinet and provide air temperature at the inlet of the equipment below the limit imposed by the ETSI. This led us to design a prototype of thermosyphon loop with a heat output that must be dissipated is imposed. We showed that this prototype is used to dissipate the heat ratings up to 470 W using a small load of npentane.We conducted tests on the prototype cooling telecommunication cabinet using slightly modified thermosyphon loop. We show that the thermal performance obtained by using a cooling mode loop thermosyphon are better. Thermosyphon loops seem interesting for passive cooling of equipment deployed in a telecommunications network. Indeed, being able to use an air conditioning system independent and requires no energy should be promoted in a reduction of overall energy consumption.

Une étude tridimensionnelle du déferlement plongeant des vagues est proposée sur la base de résultats obtenus par Simulation des Grandes Echelles des équations de Navier-Stokes en formulation diphasique. Après une étude bibliographique situant le problème, la première partie du travail est consacrée au développement et à la validation des méthodes numériques employées pour la résolution des équations. Dans une seconde partie de travail, l'outil numérique est exploité afin d'étudier différents problèmes bidimensionnels de propagation d'ondes solitaires. La dernière partie du travail consiste à appliquer les méthodes numériques au déferlement tridimensionnel d'ondes sinusoïdales périodiques de larges amplitudes. Une description complète du phénomène est proposée, incluant la génération du splash-up ainsi que la dissipation de l'énergie par les différentes structures turbulentes présentes dans l'écoulement.

Ce travail de thèse est consacré à l'étude numérique et théorique de la modulation de la turbulence par des particules. Cette étude s'appuie sur des résultats issus de simulations de type Euler/Lagrange qui résolvent directement les équations instantanées de la phase gazeuse et effectuent un suivi de trajectoires des particules. La configuration étudiée représente une nappe de particules injectées à haute vitesse dans une turbulence homogène isotrope décroissante. Le mouvement des particules est supposé uniquement gouverné par la force de traînée visqueuse. Le chargement en particules est suffisamment important pour que les particules influent sur la phase gazeuse (couplage inverse) mais suffisamment faible pour pouvoir négliger les collisions interparticulaires. Une analyse des équations de transport des principales grandeurs moyennes de l'écoulement est menée pour déterminer les effets directs et indirects des particules sur la turbulence fluide. L'étude des transferts d'énergie entre phases montre que la présence des particules tend à détruire la turbulence gazeuse au centre de la nappe et à l'augmenter à la périphérie. Ce dernier effet est causé par la forte corrélation entre la distribution de particules et la vitesse instantanée du gaz. Le modèle k-epsilon est ensuite étudié et la validité de ses hypothèses de fermeture en écoulement diphasique est éprouvée à l'aide de tests a priori. une nouvelle formulation de type viscosité turbulente, fonction des paramètres diphasiques, est utilisée pour modéliser le tenseur de Reynolds du gaz. Une équation de Langevin diphasique est également testée pour modéliser les équations de vitesse de dérive et de covariance des fluctuations de vitesse fluide-particules.

La thèse est une contribution expérimentale et théorique à la modélisation d'écoulements de bruts pétroliers produits par dépressurisation du réservoir. Il s'agit de produire l'huile en dessous du point de bulle et de profiter de l'expansion de la phase gaz pour déplacer l'huile (procédé appelé "Solution Gas Drive"). Jusqu'à présent les simulateurs de réservoir permettaient la simulation d'un tel mécanisme pour les huiles légères. Un problème se pose avec les huiles lourdes dont la forte viscosité et la physico-chimie intensifient l'état hors d'équilibre dû au mécanisme de changement de phase. Les simulateurs numériques classiques ne prennent pas en compte ce type de mécanisme hors équilibre, ni l'écoulement du gaz sous forme dispersée. Bien qu'un grand nombre de travaux aient déjà été réalisés essentiellement à l'échelle des pores, l'étude bibliographique montre qu'il est nécessaire d'avoir un modèle décrit par des équations continues (approche de Darcy). L'étude expérimentale des propriétés physico-chimiques (tension dynamique, élasticité de surface) des principaux constituants des bruts montre que les asphaltènes ont un effet tensioactif favorisant la formation de mousse au dessus d'une concentration seuil. Enfin, l'étude théorique des différents processus physiques permet d'établir un modèle d'écoulement continu, ne faisant intervenir que des grandeurs macroscopiques mesurables. La nucléation et le transfert de masse sont modélisés par une fonction volumique de transfert, la nucléation étant basée sur le modèle de bulles préexistantes. Deux phases gaz peuvent coexister : l'une continue à l'échelle du milieu, l'autre dispersée dans l'huile sous forme de bulles. Les simulations de plusieurs expériences, réalisées avec des roches et des fluides différents, montrent que le modèle permet d'interpréter et de prédire les productions de gaz et d'huile.

L'objectif de cette thèse est de développer des méthodes numériques dans un code de simulation numérique direct contenant une méthode de suivi Lagrangien de phase disperse. Ces méthodes numériques ont permis d'optimiser le suivi Lagrangien parallèle des gouttes afin de pouvoir simuler de très grands domaines. Un algorithme de gestion des collisions entre les gouttes a aussi été développé et mis en place ainsi qu'une procédure permettant d'intégrer des corps immergés de formes quelconques. Ces méthodes numériques ont permis d'effectuer des simulations de référence et d'obtenir des statistiques détaillées de diverses configurations telles que le croisement de jets de gouttes qui se collisionnent ou bien la propagation d'un front de flamme dans un spray évaporant. La comparaison de ces caractéristiques avec les mêmes simulations où la prise en charge du spray était effectuée par une méthode Eulérienne multi-fluide ont été faites. Les très bonnes concordances qualitatives et quantitatives entre les prises en charge Lagrangiennes et Eulériennes ont permis de montrer les capacités du modèle Eulérien à rendre compte de phénomènes physiques complexes. Ce modèle Eulérien permettra de bien meilleurs scaling lors de son portage sur des architectures massivement parallèles.