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.

Ce travail consiste en une étude des phénomènes de coalescence dans un nuage de gouttes, par la simulation numérique directe d'un écoulement turbulent gazeux, couplée avec une approche de suivi Lagrangien pour la phase dispersée. La première étape consiste à développer et valider une méthode de détection des collisions pour une phase polydispersée. Elle est ensuite implémentée dans un code couplé de simulation directe et de suivi Lagrangien existant. Des simulations sont menées pour une turbulence homogène isotrope de la phase continue et pour des phases dispersées en équilibre avec le fluide. L'influence de l'inertie des gouttes et de la turbulence sur le taux de coalescence des gouttes est discutée dans un régime de coalescence permanente. Un aperçu est donné de la prise en compte d'autres régimes de collision et de coalescence entre gouttes. Ces simulations sont la base de développement et de validation des approches utilisées dans les calculs à l'échelle industrielle. En particulier, les résultats des simulations sont comparés avec les prédictions d'une approche Lagrangienne de type Monte-Carlo et de l'approche Eulerienne 'Direct Quadrature Method of Moments' (DQMOM). Différents types de fermeture des termes de coalescence sont validés. Les uns sont basés sur l'hypothèse de chaos-moléculaire, les autres sont capables de prendre en compte des corrélations de vitesses des gouttes avant la collision. Il est montré que cette derniere approche prédit beaucoup mieux le taux de coalescence par comparaison avec les résultats des simulations déterministes. / Coalescence in a droplet cloud is studied in this work by means of direct numerical simulation of the turbulent gas flow, which is coupled with a Lagrangian tracking of the disperse phase. In a first step, a collision detection algorithm is developed and validated, which can account for a polydisperse phase. This algorithm is then implemented into an existing code for direct numerical simulations coupled with a Lagrangian tracking scheme. Second, simulations are performed for the configuration of homogeneous isotropic turbulence of the fluid phase and a disperse phase in local equilibrium with the fluid. The influence of both droplet inertia and turbulence intensity on the coalescence rate of droplets is discussed in a pure permanent coalescence regime. First results are given, if other droplet collision outcomes than permanent coalescence (i.e. stretching and reflexive separation) are considered. These results show a strong dependence on the droplet inertia via the relative velocity of the colliding droplets at the moment of collision. The performed simulations serve also as reference data base for the development and validation of statistical modeling approaches, which can be used for simulations of industrial problems. In particular, the simulation results are compared to predictions from a Lagrangian Monte-Carlo type approach and the Eulerian 'Direct Quadrature Method of Moments' (DQMOM) approach. Different closures are validated for the coalescence terms in these approaches, which are based either on the assumption of molecular-chaos, or based on a formulation, which allows to account for the correlation of droplet velocities before collision by the fluid turbulence. It is shown that the latter predicts much better the coalescence rates in comparison with results obtained by the performed deterministic simulations.

Les aquifères fissurés de type karst contiennent d'importantes ressources en eau. Ces aquifères sont complexes et hétérogènes sur une gamme d'échelles importantes. Leur gestion nécessite l'utilisation d'outils et de méthodologies adaptés. Dans le cadre de cette étude, différents outils et méthodologies numériques d'étude ont été développés pour la modélisation des aquifères karstiques, et plus généralement, des milieux poreux très fissurés 2D et 3D - en mettant l'accent sur la morphologie et sur le comportement hydrodynamique du milieu à travers la notion de changement d'échelle ("second changement d'échelle", reposant sur un modèle d'écoulement local de type Darcy et/ou Poiseuille avec quelques généralisations). Plusieurs axes sont explorés concernant la morphologie du milieu poreux fissuré (milieux aléatoires, milieux booléens avec réseaux statistiques de fissures, mais aussi, modèles morphogénétiques). L'étude du changement d'échelle hydrodynamique tourne autour du concept de macro perméabilité. Dans un premier temps, l'étude porte sur un modèle de perte de charge linéaire darcien. Les perméabilités effectives sont calculées numériquement en termes des fractions volumiques de fissures et du contraste de perméabilité matrice/fissures. Elles sont analysées et comparées à des modèles théoriques (analytiques). Une étude particulière des effets de quasi-percolation pour les grands contrastes aboutit à la définition de trois fractions critiques liées à des seuils de percolation. Pour tenir compte des effets inertiels dans les fissures, l'étude est étendue au cas d'une loi locale comprenant un terme quadratique en vitesse (Darcy/Ward-Forchheimer). Une perméabilité macroscopique équivalente non linéaire est définie et analysée à l'aide d'un modèle inertiel généralisé (linéaire/puissance). Enfin, l'anisotropie hydraulique à grande échelle du milieu fissuré est étudiée, en termes de perméabilités directionnelles, à l'aide d'une méthode numérique d'immersion. / Karstic aquifers contain large subsurface water resources. These aquifers are complex and heterogeneous on a large range of scales. Their management requires appropriate numerical tools and approaches. Various tools and numerical methodologies have been developed to characterize andmodel the geometry and hydraulic properties of karstic aquifers, more generally, of highly fissured 2D and 3D porous media. In this study, we emphasize morphological characterization, and we analyze hydrodynamic behavior through the concept of upscaling ("second upscaling"). Concerning the morphology of fissured porous media, several axes are explored : random media, composite random Boolean media with statistical properties, and morphogenetic models. Hydrodynamic upscaling is developed using the macro-permeability concept. This upscaling method is based on either Darcy's linear law, or on a linear/quadratic combination of Darcy's and Ward-Forchheimer's quadratic law (inertial effects). First, the study focuses on Darcy's linear head loss law, and Darcian effective permeabilities are calculated numerically in terms of volume fractions of fissures and "fissure/matrix" permeability contrasts. The results are analysed and compared with analytical results and bounds. A special study of percolation and quasi-percolation effects, for high contrasts, leads to defined three critical fractions. These critical fractions are "connected" to percolation thresholds. Secondly, in order to consider inertial effect in fissures, the study is extended to a local law with a quadratic velocity term (Darcy/Ward-Forchheimer). Then, an equivalent nonlinear macroscopic permeability is defined and analysed using a generalized inertial model (linear/power). Finally, the large scale hydraulic anisotropy of fissured medium is studied, in terms of directional permeabilities, using an "immersion" numerical method.

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.

Dans ce travail, on s'intéresse au développement d'une méthode hybride qui couple spatialement les schémas d'intégration temporelle explicite et implicite. Afin de répondre aux contraintes du solveur industriel FLUSEPA, les schémas explicite Heun et implicite Crank-Nicolson ont été hybridés via un paramètre de transition : l'approche mise en place est appelée schéma AION. Cette dernière est étudiée en détails avec une attention particulière sur son comportement spectral et sa capacité à maintenir l'ordre de précision. On montre que le traitement hybride a d'intéressants comportements dissipatif et dispersif tout en empêchant la réflexion d'ondes parasites et en maintenant la précision attendue. De plus, l'approche hybride est validée sur plusieurs cas académiques à la fois pour les flux convectifs et pour les flux diffusifs. Et comme espéré, la méthode est plus intéressante en terme de temps de calcul que les méthodes standards d'intégration temporelle. Pour l'extension de cette approche à la méthode temporelle adaptative présente dans FLUSEPA, il a été nécessaire d'améliorer le traitement qui permet à la méthode d’être conservative tout en obtenant des propriétés spectrales acceptables. Finalement l'approche hybride a été aussi étendue pour la modélisation RANS/LES de la turbulence avec des temps de calcul intéressants tout en capturant la physique de l'écoulement / This work deals with the design of a hybrid time integrator that couples spatially explicit and implicit time integrators. In order to cope with the industrial solver of Ariane Group called FLUSEPA, the explicit scheme of Heun and the implicit scheme of Crank-Nicolson are hybridized using the transition parameter : the whole technique is called AION time integration. The latter is studied into details with special focus on spectral behaviour and on its ability to keep the accuracy. It is shown that the hybrid technique has interesting dissipation and dispersion properties while maintaining precision and avoiding spurious waves. Moreover, this hybrid approach is validated on several academic test cases for both convective and diffusive fluxes. And as expected the method is more interesting in term of computational time than standard time integrators. For the extension of this hybrid approach to the temporal adaptive method implemented in FLUSEPA, it was necessary to improve some treatments in order to maintain conservation and acceptable spectral properties. Finally the hybrid time integration was also applied to a RANS/LES turbulent test case with interesting computational time while capturing the flow physics.

Ce travail de thèse concerne l’étude des écoulements de fluides gazeux compressibles laminaires subsonique-supersonique dans une tuyère de type convergent-divergent. Les écoulements étudiés sont à nombres de Reynolds modérés et s’affranchissent de l’hypothèse de condition adiabatique de paroi couramment utilisée afin de mieux prendre en compte les phénomènes de transfert de chaleur par convection et rayonnement avec le milieu extérieur. Cette étude des phénomènes de transferts conjugués a permis de déterminer le comportement dynamique simultané du fluide et de la paroi de la tuyère. Enfin, compte tenu des niveaux élevés de températures mis en jeu dans ces systèmes, une analyse concernant le comportement thermomécanique de l’ensemble de la structure de paroi avec des matériaux monocouches et multicouches a été réalisé. De nombreuses configurations géométriques, propriétés physiques et conditions aux limites sur le fluide et la paroi ont été analysées. Les résultats présentés montrent, la structure des écoulements à travers les iso-contours de vitesses, des nombres de Mach, des pressions dans le fluide, des températures dans le fluide et dans la paroi ainsi que les déformations et les contraintes de la paroi qui résultent des couplages thermomécaniques. Une analyse des performances de la tuyère, en termes de force de poussée et de coefficient de débit spécifique, est largement discutée dans ce travail. / This work concerns the study of flows of compressible gaseous laminar subsonic-supersonic nozzle in a convergent-divergent type. The flows are studied to moderate Reynolds numbers and free themselves from the assumption of adiabatic wall conditions commonly used to better take into account the phenomena of heat transfer by convection and radiation with the external environment. This study combined transfer phenomena was determined simultaneously the dynamic behavior of the fluid and the wall of the nozzle. Finally, given high levels of temperatures at stake in these systems, an analysis of the thermomechanical behavior of the entire wall structure with monolayer and multilayer materials was performed. Many geometric configurations, physical properties and boundary conditions on the fluid and the wall were analyzed. The results presented show the structure of the flow through the iso-contours of speed, Mach numbers, pressures in the fluid, temperatures in the fluid and in the wall as well as the deformations and stresses resulting from the wall thermomechanical couplings. A performance analysis of the nozzle, in terms of thrust coefficient and specific yield, is widely discussed in this work.

Le refroidissement seculaire du manteau terrestre et de la lithosphere continentale se traduit par l'augmentation de la surface de terres emergees. L'augmentation corollaire des processus d'alteration et d'erosion des silicates a des consequences majeures pour les enveloppes externes. Nous avons developpe un modele physique qui permet d'evaluer la surface de terres emergees en fonction de la temperature du manteau, de la surface totale de continents, et de la distribution des altitudes continentales. Nos resultats numeriques montrent qu'a la fin de l'Archeen, moinsde 15% de la surface terrestre etaient emergee, en accord avec nombre d'observations geologiques et geochimiques. Pour estimer le refroidissement seculaire de la lithosphere continentale, nous avons combine des modeles thermo-mecaniques avec des observations de terrain. Nos resultats, contraints par des donnees geologiques, suggerent que la temperature au Moho a diminue de ~ 200ºC en 2,7 Ga dans le craton des Pilbaras. Pour evaluer l'effet de la croissance continentale sur l'evolution de la surface de terres emergees, nous avons developpe un modele base sur un modele d'evolution thermique publie. Nos resultats suggerent que la surface emergee, de moins de 5% de la surface terrestre a l'Archeen, depend peu de la croissance continentale. Ceci permet de reconcilier l'evolution du 87Sr/86Sr oceanique avec une croissance continentale precoce. Les continents sont enrichis en phosphate, element essentiel a la biosphere. Leur emergence aurait donc engendre une augmentation de la production d'oxygene par des micro-organismes photosynthetiques, contribuant ainsi a l'oxydation de l'atmosphere il y a 2,4 Ga.

L'objet de cette étude est d'utiliser la décomposition en modes propres (POD) pour établir un modèle de la composante déterministe et aléatoire d'un écoulement turbulent présentant une instationnarité à grande échelle pseudo périodique.<br />Une étude expérimentale de l'écoulement bidimensionnel en moyenne en aval d'un cylindre semi-circulaire, par vélocimétrie par image de particules en deux temps (PIV2T) caractérise l'écoulement<br />Une analyse POD du champ de vitesse permet d'extraire les modes spatiaux et de définir un paramètre de phase décrivant l'instationnarité à grande échelle qui régit la partie déterministe. La modélisation de l'évolution temporelle des coefficients associés aux modes s'effectue par des fonctions soit harmoniques pour la partie déterministe, soit stochastiques pour la partie aléatoire.<br />La modélisation est en accord avec les mesures expérimentales des premiers moments statistiques en un point et des fonctions de corrélation spatio-temporelle du champ de vitesse.

L'écoulement de von Karman est l'écoulement produit par la rotation de deux disques coaxiaux placés en vis-à-vis. Dans la première partie de cette thèse, nous présentons un certain nombre de résultats expérimentaux visant à la caractérisation des fluctuations à grande échelle d'un écoulement d'allongement sensiblement unitaire. Dans certaines circonstances, nous mettons notamment à jour une "bifurcation globale" ainsi qu'un comportement "intermittent" de la structure de l'écoulement. Nous présentons par ailleurs une relation qui relie les statistiques des fluctuations de l'écoulement aux valeurs des couples fournis par les moteurs. Enfin, nous ébauchons l'étude de l'influence de la force de Coriolis sur l'écoulement. La deuxième partie de cette thèse présente des résultats numériques et expérimentaux portant sur l'étude des caractéristiques magnéto-hydrodynamiques de l'écoulement, et notamment sur sa capacité à produire un champ magnétique par "effet dynamo". Les résultats expérimentaux ont été obtenus à haut nombre de Reynolds magnétique dans le montage "VKS", qui utilise du sodium liquide.