Ce travail porte sur la modélisation de la formation et la dispersion d'un nuage de gouttes, par déconfinement d'un liquide: agression extérieure ou situation accidentelle. Le but est la construction d'un modèle apte à reproduire simultanément les conditions génératrices de la formation du nuage et l'évolution de ce nuage dans le temps (dispersion). La principale difficulté réside en la différence des modèles adaptés à la description d'écoulements caractérisant chaque étape du phénomène global : modèle d'écoulement multiphasique à phases compressibles (milieux continus) initialement, puis fragmentation et formation du nuage de gouttes dispersées dans une phase porteuse (modèle d'écoulements dilués). En l'absence de modèle analytique unique apte à décrire l'ensemble de ces processus, on propose une approche originale pour réaliser un couplage effectif entre ces deux modèles. La problématique de formation et de dispersion de liquide implique la prise en compte de plusieurs phénomènes physiques: fragmentation, transferts de chaleur et de masse ainsi que la traînée entre les phases. Ces différents phénomènes sont introduits dans le modèle global via des termes d'interactions présents dans les systèmes d'équations. La construction de ce modèle complet à permis la réalisation de calculs décrivant la formation et la dispersion d'un nuage de gouttes pouvant intervenir lors de situations accidentelles sur des sites industriels par exemple. / This work focuses on modeling the formation and the dispersion of a cloud of droplets, induced by ejection of a liquid, resulting from an external aggression or an accidental situation. The goal is to build a model able to reproduce simultaneously the conditions which generate the cloud formation and the cloud evolution in time (dispersion). The main difficulty lies in the differences between the already existing models adapted to the description of flows which are able to characterize each stage of the global phenomenon: initially a multiphase flow model with compressible phases (Continuum), then the atomization and the formation of a cloud of droplets dispersed in a carrier phase (dilute flow model). We propose a new approach to achieve an effective coupling between these two models. The problem of the formation and the dispersion of the liquid requires to take into account several physical phenomena: atomization, heat and mass transfers and drag between phases. These phenomena are included in the global model through interaction terms involved in the systems of equations. The construction of this model has permited the realization of calculations describing the formation and dispersion of a cloud of droplets which may occur during, for axample, in accidental situations at industrial sites.

Le bon fonctionnement et les performances des turbines à vapeur sont liés à l’état de la vapeur et notamment au taux d’humidité qu’elle contient. EDF souhaite pouvoir maîtriser les phénomènes spécifiques à ces problématiques afin d’améliorer l’utilisation et l’évolution de ses turbines. Le sujet de recherche concerne la modélisation de la formation de l’humidité dans un corps de turbine et l’étude des couplages entre la phase liquide et les instationnarités. Dans ce contexte, la démarche adoptée est la suivante : la présence d’humidité est prise en compte à l’aide d’un modèle homogène, couplé à des modèles de condensation permettant de prendre en compte les phénomènes hors-équilibre thermodynamique : le grossissement et la nucléation des gouttes d’eau dans la vapeur. Pour mener à bien les calculs, des méthodes numériques adaptées aux gaz réels ont été utilisées et testées à l’aide d’un code monodimensionnel avant d’être intégrées dans le code 3D elsA. Deux types de modèles de condensation ont été mis en œuvre, considérant ou non la polydispersion des gouttes dans la vapeur. Les couplages instationnaires entre la condensation et l’écoulement principal ont été étudiés à différents niveaux d’observations (1D, 1D − 3D, 3D). Il a été montré que la méthode des moments apporte une richesse supplémentaire par rapport à un modèle mono-dispersé, et permet de mieux capter les couplages instationnaires entre l’humidité et le champ principal. / In addition to conventional turbomachinery problems, both the behavior and performances of steam turbines are highly dependent on the vapour thermodynamic state and the presence of a liquid phase. EDF, the main French electricity producer, is interested in further developing its’ modelling capabilities and expertise in this area to allow for operational studies and long-term planning. This PhD thesis explores the modelling of wetness formation and growth in a steam turbine and an analysis of the coupling between the liquid phase and the main flow unsteadiness. To this end, the work in this thesis took the following approach. Wetness was accounted for using a homogeneous model coupled with transport equations to take into account the effects of non-equilibrium phenomena, such as the growth of the liquid phase and nucleation. The real gas attributes of the problem demanded adapted numerical methods. Before their implementation in the 3D elsA solver, the accuracy of the chosen models was tested using a developed one-dimensional nozzle code. In this manner, various condensation models were considered, including both polydispersed and monodispersed behaviours of the steam. Finally, unsteady coupling effects were observed from several perspectives (1D, 1D − 3D, 3D), demonstrating the ability of the method of moments to sustain unsteady phenomena which were not apparent in a simple monodispersed model.

Une méthode de frontières immergées est développée pour la simulation d’écoulements compressibles et validée au travers de cas-tests spécifiques (réflexion d’ondes acoustiques et quantification de la conservation de la masse dans des canaux inclinés). La simulation aux grandes échelles (LES) d’une cavité transsonique est ensuite présentée. Le bouclage aéro-acoustique, très sensible aux conditions aux limites, est reproduit avec précision par la LES dans le cas où les parois sont immergées dans un maillage structurée. La comparaison des stratégies de modélisation de sous-maille pour cet écoulement transsonique et l’adaptation des filtres en présence de frontières immergées sont également discutées. Le rôle, souvent sous-estimé, du schéma de viscosité artificiel, est quantifié.Dans la dernière partie du manuscrit, des études sont réalisées pour aider au dimensionnement d’un nouveau concept de chambre de combustion où la flamme est stabilisée par la recirculation de gaz brûlés dans une cavité (chambre TVC pour Trapped Vortex Combustor). La modélisation de la combustion turbulente est basée sur une chimie tabulée, couplée à une fonction densité de probabilité présumée (PCM-FPI). L’étude de la dynamique de la flamme est réalisée pour diverses conditions de fonctionnement (débit de l’écoulement principal et présence ou non d’un swirl). Les spécificités de mise en œuvre de la simulation d’un écoulement de ce type sont discutées et un soin particulier est apporté au traitement de la condition de sortie, qui constitue un point sensible de la chaîne de modélisation. Les phénomènes d’instabilités et de retour de la flamme sont mis en évidence ainsi que les modifications à apporter au dispositif afin de minimiser ces effets. L’existence d’un cycle limite acoustique est souligné et une formule permettant d’anticiper le niveau des fluctuations de pression est proposée et validée. Une correction au modèle PCM-FPI est présentée afin de préserver la vitesse de flamme et d’assurer une reproduction plus précise de la dynamique de flamme. / An immersed boundary method has been developed for the simulation of compressible flow and validated with reference test cases (pressure wave reflection and quantification of mass conservation for various inclined channels). Large Eddy Simulation (LES) of a transonic cavity is then presented. The aeroacoustic feedback loop, which is highly sensitive to the boundary conditions, was accurately reproduced where the walls are immersed inside a structured grid. The comparison between the modeling approaches for this transonic flow and the correction of the filtering operation near immersed boundaries are also discussed. The often underestimated role of the numerical artificial dissipation is also quantified.In the last part of this manuscript, many studies are realized to help in the design of a new combustion chamber for Trapped Vortex Combustor (TVC). The turbulent combustion model is based on tabulated chemistry and a presumed probability density function (PCM-FPI) method.The flame dynamics is studied for various operating conditions (flowrate of the main flow and presence of swirl motion). Details concerning the realization of such a flow are discussed and special care is taken for the treatment of the most sensitive outlet boundary condition. The phenomena of combustion instabilities and of flame backflow are highlighted along with the modifications to be made for the device to minimize these effects. The existence of a acoustic limit cycle is emphasized and a formula is proposed and validated to anticipate the level of pressure fluctuations. Finally a correction to the PCM-FPI model is suggested to preserve the flame front speed and to ensure a more accurate description of the flame dynamics.

Les travaux de thèse sont axés sur les méthodes numériques pour les écoulements diphasiques, compressibles, à faible vitesse, avec apparition soudaine de forts gradients de pression. La vitesse matérielle de chacune des phases étant très petite devant la célérité des ondes acoustiques, le régime d'écoulement est dit à faible nombre de Mach. Dans ce travail, la loi d'état de la phase considérée contient toujours une information mesurant sa plus ou moins grande compressibilité. Ainsi, la faible compressibilité de l'eau peut produire un régime d'écoulement où des sauts de pression importants apparaissent même si le nombre de Mach est très faible. La première partie de la thèse s'est focalisée sur un modèle diphasique dit homogène-équilibré. Les deux phases de l'écoulement ont alors la même vitesse, pression, température et même potentiel chimique. Un premier travail a été la construction de solveurs de Riemann approchés dits tout-nombre-de-Mach. En l'absence de transitoire rapide, ces solveurs basent leur contrainte de pas de temps sur la vitesse des ondes matérielles lentes et sont donc précis pour suivre ces dernières. En revanche, lorsqu'une onde de choc rapide traverse l'écoulement, ces solveurs s'adaptent automatiquement afin de la capturer. La seconde partie de la thèse s'est focalisée sur la prise en compte du couplage convection-source dans le cadre des modèles en approche bi-fluides avec effets de relaxation pression-vitesse. Dans ces modèles, les deux phases de l'écoulement possèdent leur propre jeu de variables. Dans ce travail, un schéma implicite à mailles décalées, basé sur l'influence des termes sources dans des problèmes de Riemann linéaires, a été proposé / The present work focuses on numerical methods for low-material velocity compressible two-phase flows with high pressure jumps. In this context, the material velocity of both phases is small compared with the celerity of the acoustic waves. The flow is said to be a low-Mach number flow. In this work, the equation of state of the considered phase always contains information relative to its compressibility. For example, the low-compressibility of liquid water may lead to fast transients in which high pressure jumps are produced even if the flow Mach number is low. The first part of this work has leaned on two-phase homogeneous-equilibrium models. Thus, both phases have the same velocity, pressure, temperature and the same chemical potential. The construction of what is called an all-Mach-number approximate Riemann solver has been conducted. When no fast transients come through the flow, the above solvers enable computations with CFL conditions based on low-material velocities. As a result, they remain accurate to follow slow material interfaces, or subsonic contact discontinuities. However, when fast shock waves propagate, these solvers automatically adapt in order to capture them. The second part of the thesis has been dedicated to the design of numerical methods enhancing the coupling between convection and relaxation for two-fluid models containing pressure-velocity relaxation effects. In such models, both phases have their own set of variables. A time-implicit staggered scheme, based on the influence of relaxation source terms on linear Riemann problems has been proposed.

Cette thèse contribue à la simulation numérique des écoulements d'un fluide compressible modélisé par les équations de Euler et Navier-Stokes: étude d'un schéma d'ordre élévé basé sur une matrice de masse, modélisation des écoulement turbulents compressibles à très haut Reynolds, développement des estimateurs et correcteurs d'erreurs a posteriori et a priori, et adaptation de maillage anisotrope pour les fonctionnelles d'observation. Concernant la prédiction des écoulements turbulents, on s'est intéressé aux modèles hybrides de type RANS/LES comportant les nouveautés suivantes: traitement des tourbillons de grande échelle utilisant la formulation VMS (Variational Multi-Scale) et du RANS employé sur la paroi sur une distance imposée via une zone de protection conçue afin d'éviter le phénomène assez commun apelé "grid induced model depletion". Le niveau de viscosité du modèle VMS-LES est de plus controllé par un procédé de double filtre dynamique. La seconde partie concerne l'adaptation de maillage anisotrope pour mieux observer une fontionnelle d'observation. Les estimations a priori sont réalises pour le modèle des équations d'Euler et Navier-Stokes en instationnaire en 2D et 3D. A partir de ces estimations on sait définir les maillages optimaux au cours du calcul instationnaire, en fonction de l'état et de l'état adjoint. Le système d'optimalité est discrétisé et résolu à l'aide d'une méthode de point fixe instationnaire global, comportant une stratégie de stockage/recalcul pour le couplage état/ état adjoint. Des applications à la propagation d'ondes de choc et d'ondes acoustiques sont présentées.

Lors que nous examinons numériquement des phénomènes multiphasiques suite à un accidentgrave dans le réacteur nucléaire, la dimension caractéristique des zones multi-fluides(non-réactifs et réactifs) s’avère beaucoup plus petite que celle du bâtiment réacteur, cequi fait la Simulation Numérique Directe de la configuration à peine réalisable. Autrement,nous proposons de considérer la zone de mélange multiphasique comme une interface infinimentfine. Puis, le solveur de Riemann réactif est inséré dans la Méthode des ÉquationsDiscrètes Réactives (RDEM) pour calculer le front de combustion à grande vitesse représentépar une interface discontinue. Une approche anti-diffusive est ensuite couplée avec laRDEM afin de précisément simuler des interfaces réactives. La robustesse et l’efficacité decette approche en calculant tant des interfaces multiphasiques que des écoulements réactifssont à la fois améliorées grâce à la méthode ici proposée : upwind downwind-controlled splitting(UDCS). UDCS est capable de résoudre précisément des interfaces avec les maillagesnon-structurés multidimensionnels, y compris des fronts réactifs de détonation et de déflagration. / When numerically investigating multiphase phenomena during severe accidents in a reactorsystem, characteristic lengths of the multi-fluid zone (non-reactive and reactive) are foundto be much smaller than the volume of the reactor containment, which makes the directmodeling of the configuration hardly achievable. Alternatively, we propose to consider thephysical multiphase mixture zone as an infinitely thin interface. Then, the reactive Riemannsolver is inserted into the Reactive Discrete Equations Method (RDEM) to compute highspeed combustion waves represented by discontinuous interfaces. An anti-diffusive approachis also coupled with RDEM to accurately simulate reactive interfaces. Increased robustnessand efficiency when computing both multiphase interfaces and reacting flows are achievedthanks to an original upwind downwind-controlled splitting method (UDCS). UDCS is capableof accurately solving interfaces on multi-dimensional unstructured meshes, includingreacting fronts for both deflagration and detonation configurations.

Cette thèse présente la construction de schémas distribuant le résidu (RD) d'ordre très élevés, pour la discrétisation d'équations d'advection-diffusion multidimensionnelles et stationnaires sur maillages non structurés. Des schémas linéaires ainsi que des schémas non linéaires sont considérés. Une approximation de la solution polynomiale par morceaux et continue sur chaque élément est adoptée, de plus une procédure de reconstruction du gradient que celle de la solution numérique est utilisée afin d'avoir une représentation continue de la solution numérique et de son gradient. Il est montré que le gradient doit être reconstruit avec la même précision de la solution, sans quoi la précision formel du schéma numérique est perdue dans les cas où les effets de diffusion prévalent sur les effets d'advection, et aussi quand l'advection et la diffusion sont également importants. Ensuite, la méthode est étendue à des systèmes d'équations, en particulier aux équations de Navier-Stokes et aux équations RANS. La précision, l'efficacité et la robustesse du solveur RD implicite sont démontrées sur plusieurs cas tests. / The construction of compact high-order Residual Distribution schemes for the discretizationof steady multidimensional advection-diffusion problems on unstructuredgrids is presented. Linear and non-linear scheme are considered. A piecewise continuouspolynomial approximation of the solution is adopted and a gradient reconstructionprocedure is used in order to have a continuous representation of both thenumerical solution and its gradient. It is shown that the gradient must be reconstructedwith the same accuracy of the solution, otherwise the formal accuracy ofthe numerical scheme is lost in applications in which diffusive effects prevail overthe advective ones, and when advection and diffusion are equally important. Thenthe method is extended to systems of equations, with particular emphasis on theNavier-Stokes and RANS equations. The accuracy, efficiency, and robustness of theimplicit RD solver is demonstrated using a variety of challenging aerodynamic testproblems.

Cette thèse présente la construction d’un schéma aux Résidus Distribués p-adaptatif pour la discrétisation des équations d’Euler ainsi qu’un schéma aux Résidus Distribués hp-adaptatif pour les équations de Navier- Stokes pénalisées. On rappelle tout d’abord les équations d’Euler et de Navier-Stokes ainsi que leurs versions non dimensionnelles. Les définitions et propriétés de base des schémas aux Résidus Distribués sont ensuite présentées. On décrit alors la construction d’un schéma aux Résidus Distribués p-adaptatif pour les équations d’Euler. La construction du schéma p-adaptatif est basée sur la possibilité d’exprimer le résidu total d’un élément K de degré k (au sens où l’élément fini (K; P; Sigma ) est un élément fini de degré k) comme une somme pondérée des résidus totaux de ses sous-éléments de degré 1. La solution discrète ainsi obtenue est en général discontinue à l’interface entre un élément subdivisé et un élément non subdivisé. Ceci contredit l’hypothèse de continuité de la solution qui est utilisée pour démontrer le théorème de Lax-Wendroff discret pour les schémas aux Résidus Distribués. Cependant, on montre que cette hypothèse peut être assouplie. La conséquence pratique est que si l’on emploie des quadratures particulières dans l’implémentation numérique, on peut quand même démontrer le théorème de Lax-Wendroff discret, ce qui garantit la convergence du schéma numérique vers une solution faible des équations d’origine. Les formules qui permettent d’exprimer le résidu total comme une somme pondérée des résidus totaux des sous-éléments sont à la base de la méthode de p-adaptation présentée ici. Dans le cas quadratique, la formule est obtenue avec les classiques fonctions de base de Lagrange en dimension deux et avec des fonctions de base de Bézier en dimension trois. Ces deux formules sont ensuite généralisées à des degrés polynomiaux quelconques en dimension deux et trois avec des fonctions de base de Bézier. Dans la deuxième partie de la thèse, on présente l’application du schéma p-adaptatif aux équations pénalisées de Navier-Stokes avec adaptation de maillage anisotrope. . En pratique, on combine le schéma p-adaptatif avec la méthode IBM-LS-AUM (Immersed Boundary Method with Level Sets and Adapted Unstructured Meshes). La méthode IBM-LS-AUM permet d’imposer les conditions aux bords grâce à la méthode de pénalisation et l’adaptation anisotrope du maillage à la solution numérique et à la level-set augmente la précision de la solution et de la représentation de la surface. Une fois la méthode IBM-LS-AUM combinée avec le schéma p-adaptatif, il est alors possible d’utiliser des éléments d’ordre élevés en-dehors de la zone où la pénalisation est appliquée. La méthode est robuste comme le montrent les diverses expérimentations numériques à des vitesses faibles à élevées et à différents nombres de Reynolds. / This thesis presents the construction of a p-adaptive Residual Distribution scheme for the steady Euler equations and a hp-adaptive Residual Distribution scheme for the steady penalized Navier-Stokes equations in dimension two and three. The Euler and Navier-Stokes equations are recalled along with their non dimensional versions. The basis definitions and properties of the steady Residual Distribution schemes are presented. Then, the construction of a p-adaptive Residual Distribution scheme for the Euler equations is considered. The construction of the p-adaptive scheme is based upon the expression of the total residual of an element of a given degree k (in the Finite Element sense) into the total residuals of its linear sub-elements. The discrete solution obtained with the p-adaptive scheme is then a one degree polynomial in the divided elements and a k-th degree polynomial in the undivided ones. Therefore, the discrete solution is in general discontinuous at the interface between a divided element and an undivided one. This is in apparent contradiction with the continuity assumption used in general to demonstrate the discrete Lax-Wendroff theorem for Residual Distribution schemes. However, as we show in this work, this constrain can be relaxed. The consequence is that if special quadrature formulas are employed in the numerical implementation, the discrete Lax-Wendroff theorem can still be proved, which guaranties the convergence of the p-adaptive scheme to a weak solution of the governing equations. The formulas that express the total residual into the combination of the total residuals of the sub-elements are central to the method. In dimension two, the formula is obtained with the classical Lagrange basis in the quadratic case and with the Bézier basis in dimension three. These two formulas are then generalized to arbitrary polynomial degrees in dimension two and three with a Bézier basis. In the second part of the thesis the application of the p-adaptive scheme to the penalized Navier-Stokes equations with anisotropic mesh adaptation is presented. In practice, the p-adaptive scheme is used with the IBM-LS-AUM (Immersed Boundary Method with Level Sets and Adapted Unstructured Meshes) method. The IBM-LS-AUM allows to impose the boundary conditions with the penalization method and the mesh adaptation to the solution and to the level-set increases the accuracy of the representation of the surface and the solution around walls. When the IBM-LSAUM is combined with the p-adaptive scheme, it is possible to use high-order elements outside the zone where the penalization is applied. The method is robust as shown by the numerical applications at low to large Mach numbers and at different Reynolds in dimension two and three.

L’objectif de ce travail de thèse est de proposer la Méthodes des éléments finis de Galerkin discontinus (DGFEM) à la discrétisation des équations compressibles de Navier-Stokes. Plusieurs challenges font l’objet de ce travail. Le premier aspect a consisté à montrer l’ordre de convergence optimal de la méthode DGFEM en utilisant les polynômes d’interpolation d’ordre élevé. Le deuxième aspect concerne l’implémentation de méthodes de ‘‘shock-catpuring’’ comme les limiteurs de pentes et les méthodes de viscosité artificielle pour supprimer les oscillations numériques engendrées par l’ordre élevé (lorsque des polynômes d’interpolation de degré p>0 sont utilisés) dans les écoulements transsoniques et supersoniques. Ensuite nous avons implémenté des estimateurs d’erreur a posteriori et des procédures d ’adaptation de maillages qui permettent d’augmenter la précision de la solution et la vitesse de convergence afin d’obtenir un gain de temps considérable. Finalement, nous avons montré la capacité de la méthode DG à donner des résultats corrects à faibles nombres de Mach. Lorsque le nombre de Mach est petit pour les écoulements compressibles à la limite de l’incompressible, la solution souffre généralement de convergence et de précision. Pour pallier ce problème généralement on procède au préconditionnement qui modifie les équations d’Euler. Dans notre cas, les équations ne sont pas modifiées. Dans ce travail, nous montrons la précision et la robustesse de méthode DG proposée avec un schéma en temps implicite de second ordre et des conditions de bords adéquats. / The aim of this thesis is to deal with compressible Navier-Stokes flows discretized by Discontinuous Galerkin Finite Elements Methods. Several aspects has been considered. One is to show the optimal convergence of the DGFEM method when using high order polynomial. Second is to design shock-capturing methods such as slope limiters and artificial viscosity to suppress numerical oscillation occurring when p>0 schemes are used. Third aspect is to design an a posteriori error estimator for adaptive mesh refinement in order to optimize the mesh in the computational domain. And finally, we want to show the accuracy and the robustness of the DG method implemented when we reach very low mach numbers. Usually when simulating compressible flows at very low mach numbers at the limit of incompressible flows, there occurs many kind of problems such as accuracy and convergence of the solution. To be able to run low Mach number problems, there exists solution like preconditioning. This method usually modifies the Euler. Here the Euler equations are not modified and with a robust time scheme and good boundary conditions imposed one can have efficient and accurate results.

Cette thèse s'intéresse aux simulations dépendantes du temps impliquant des géometries fixes ou mobiles. Ce type de simulations est l'objet d'attentes grandissantes de la part des industriels, qui souhaiteraient voir réaliser ce type de calculs de façon systématique au sein de leurs centres de recherche, ce qui n'est clairement pas le cas à l'heure actuelle. Ce travail tente de satisfaire en partie cette demande et vise notamment à améliorer la précision ainsi que l'efficacité en termes de temps de calcul des algorithmes actuellement utilisés dans ce contexte. Les méthodes d'adaptation de maillage anisotrope par prescription d'un champ de métriques, qui ont aujourd'hui atteint une certaine maturité, notamment dans leur application aux simulations stationnaires, constituent une piste très prometteuse pour l'amélioration des calculs évoluant en temps, mais leur extension dans ce contexte est loin d'être triviale. Quant à leur utilisation sur les simulations en géométries mobiles, seules quelques tentatives peuvent être répertoriées, et très peu portent sur des problèmes réalistes en trois dimensions. Cette étude présente plusieurs nouveautés sur ces questions, notamment l'extension de l'adaptation de maillage multi-échelles par champ de métriques aux problèmes instationnaires en géométries fixes et mobiles. Par ailleurs, essentiellement dans une optique de réduction des temps de calculs, une stratégie originale à été adoptée pour réaliser des calculs impliquant des maillages mobiles. Notamment, il est démontré par la pratique dans cette thèse qu'il possible de déplacer des objets en trois dimensions sur de grandes distances en maintenant le nombre de sommets du maillage constant, c'est-à-dire en limitant les types d'opérations de modification de maillage autorisés. Il en résulte un gain conséquent en terme de temps de calcul aussi bien au niveau du déplacement de maillage qu'au niveau de la résolution numérique. Par ailleurs, un nouveau schéma est proposé qui permet de gérer les changements de connectivité du maillage de manière cohérente avec la description Arbitrary-Lagrangian-Eulerian des équations physiques. La plupart de ces nouvelles méthodes ont été appliquées à la simulation d'écoulements fluides compressibles autour de géometries complexes en deux et trois dimensions d'espace.