Un nouveau procédé propre et performant de préparation de surface reposant sur l'utilisation d'un jet d'azote sous très haute pression et basse température est développé. Comparativement aux techniques classiques de préparation de surface, ce nouveau procédé apporte des avantages environnementaux (pas de génération de déchets additionnels), techniques (modes d'action sur le revêtement à traiter) et énergétiques (consommation d'énergie). Cependant, le fonctionnement du procédé de traitement de surface en question, qui existe sous forme de prototype à l'heure actuelle, repose sur des bases empiriques et les phénomènes physiques qu'il met en jeu sont encore bien mal appréhendés. Pour comprendre ces phénomènes et les modéliser, une analyse numérique est présentée. Cette analyse s'inclut dans une stratégie d'étude qui vise à augmenter progressivement le degré de complexité de la modélisation. Un premier modèle vise à décrire l'évolution du fluide en temps et en espace en supposant un écoulement compressible visqueux et axisymétrique. Dans ce modèle, le fluide est considéré comme parfait. Une extension directe de ce modèle est ensuite présentée, elle propose de tenir compte des effets de fluide réel dans l'écoulement compressible monophasique. Ces modèles sont implémentés dans un code de calculs par volumes finis. Des cas tests sont étudiés afin de valider les modèles numériques. Une étude sur des configurations de type industrielles, représentatives des conditions d'utilisation du procédé de traitement de surface par jet d'azote, est ensuite menée / A new and efficient process of surface treatment is developed and exists in prototype form at present. The process aims at injecting an inert gas, such as nitrogen, at supercritical conditions through a nozzle. The jet resulting from the expansion of the fluid at high pressure and supercritical temperature impinges normally a flat surface. Alternatively to water jet technologies, which need expansive purification of water after use, and to other classical surface treatment process, this process provides environmental (no generation of additional waste), technical (action on the coating) and energetical (efficiency) benefits. However, the physical phenomena involved in the jet are still poorly understood at present. To understand and model these phenomena, numerical analysis is presented. This analysis is part of a study strategy that aims at gradually increasing the complexity of modeling. A first model aims at describing the evolution of the fluid in time and space, assuming a compressible axisymmetric viscous flow. In this model, the fluid is assumed perfect. A direct extension of this model is then presented where real fluid effects are taken into account in the compressible flow. These models are implemented in finite volume CFD code. Test cases are studied to validate the numerical models. A study of industrial-type configurations, representative of the conditions of use of the process of surface treatment by nitrogen is then conducted

Dans cette thèse on s'est intéressé à l'étude de problèmes concernant la théorie des écoulements compressibles. Dans une première partie on a traité le problème de transport instationnaire avec un champ de vitesse peu régulier, on a établi un résultat d'existence en passant à la limite sur des schémas numériques volumes finis avec un choix décentré amont qui garantie la positivité de la masse volumique. Pour le problème de Stokes, le résultat est démontré par deux approches : une approche par schéma numérique et une approche par régularité visqueuse.Dans la première méthode on propose une discrétisation qui combine la méthode des éléments finis et la méthode des volumes finis qui repose sur les espaces Crouzeix-Raviart. Une première difficulté de ce travail est de démontrer les estimations sur la solution discrète, en particulier à cause de la présence de la gravité dans le terme source de l'équation de quantité de mouvement. Le fait de considérer une loi d'état très générale conduit des difficultés supplémentaires en particulier dans le passage à la limite sur cette équation.Dans la deuxième méthode, le résultat d'existence est démontré en utilisant une approximation par viscosité. Ceci consiste essentiellement en deux parties : l'étude du problème de convection diffusion (qui apparait dans le problème régularisé) où on démontre l'existence et l'unicité de solution et en deuxième partie le passage à la limite sur le problème régularisé. / This thesis is concerned with the study of problems relating in the theory of compressible flows . We prove the existence of the considered problems in a first part by passing to the limit on the numerical schemes proposed for the discretisation of these problems. In the second part, the existence result is obtained by passing to the limit on the approximate solutions given by a corresponding regularized problem.The main result is to prove the existence of a solution of the stationnary compressible Stokes problem with a general equation of state.We first prove this result by passing to the limit on the numerical scheme as the mesh size tends to zero. The fact to consider a general E.O.S induces some additional difficulties in particular to get estimates on the discrete solution (which comes also from the presence of the gravity in the momentum equation) and in the passage to the limit on the E.O.S.We also prove the existence result by passing to the limit on a regularized problem. We first treat the convection-diffusion problem (which appears in the regularized problem), we give an existence and uniqueness result, and we then prove estimates on the approwimate solutions and pass to the limit on the regularized problem.

La topographie d'une surface se compose généralement de plusieurs échelles, depuis l'échelle macroscopique (sa géométrie physique), jusqu'aux échelles microscopiques ou atomiques appelées rugosité. L'évolution spatiale et géométrique de la rugosité fournit une description plus complète de la surface, et une interprétation physique de certains problèmes importants tels que le frottement et les mécanismes d'usure pendant le contact mécanique entre deux surfaces. La topographie d'une surface rugueuse est de nature aléatoire, ce qui traduit par des altitudes spatialement corrélées, appelées pics et vallées. La relation entre leurs densités de probabilité et leurs propriétés géométriques sont les aspects fondamentaux qui ont été développés dans cette thèse, en utilisant la théorie des champs aléatoires et la géométrie intégrale. Un modèle aléatoire approprié pour représenter une surface rugueuse a été mis en place et étudié au moyen des paramètres les plus significatifs, dont les changements influencent la géométrie des ensembles de niveaux (excursion sets) de cette surface. Les ensembles de niveaux ont été quantifiés par des fonctionnelles connues sous le nom de fonctionnelles de Minkowski, ou d'une manière équivalente sous le nom de volumes intrinsèques. Dans un premier temps, les volumes intrinsèques des ensembles de niveaux sont calculés analytiquement sur une classe de modèles mixtes, qui sont définis par la combinaison linéaire d'un champ aléatoire Gaussien et d'un champ de t-student (t-field), et ceux d'une classe de champs aléatoires asymétriques appelés skew-t. Ces volumes sont comparés et testés sur des surfaces produites par des simulations numériques. Dans un second temps, les modèles aléatoires proposés ont été appliqués sur des surfaces réelles acquises à partir d'une cupule d'UHMWPE (provenant d’une prothèse totale de hanche) avant et après les processus d'usure. Les résultats ont montré que le champ aléatoire skew-t est un modèle mieux approprié pour décrire la rugosité de surfaces usées, contrairement aux modèles adoptés dans la littérature. Une analyse statistique, basée sur le champ aléatoire skew-t, est ensuite proposée pour détecter les niveaux des pics/vallées de la surface usée et pour décrire le comportement et la fonctionnalité de la surface usée. / Surface topography is, generally, composed of many length scales starting from its physical geometry, to its microscopic or atomic scales known by roughness. The spatial and geometrical evolution of the roughness topography of engineering surfaces avail comprehensive understanding, and interpretation of many physical and engineering problems such as friction, and wear mechanisms during the mechanical contact between adjoined surfaces. Obviously, the topography of rough surfaces is of random nature. It is composed of irregular hills/valleys being spatially correlated. The relation between their densities and their geometric properties are the fundamental topics that have been developed, in this research study, using the theory of random fields and the integral geometry.An appropriate random field model of a rough surface has been defined by the most significant parameters, whose changes influence the geometry of its excursion. The excursion sets were quantified by functions known as intrinsic volumes. These functions have many physical interpretations, in practice. It is possible by deriving their analytical formula to estimate the parameters of the random field model being applied on the surface, and for statistical analysis investigation of its excursion sets. These subjects have been essentially considered in this thesis. Firstly, the intrinsic volumes of the excursion sets of a class of mixture models defined by the linear combination of Gaussian and t random fields, then for the skew-t random fields are derived analytically. They have been compared and tested on surfaces generated by simulations. In the second stage, these random fields have been applied to real surfaces measured from the UHMWPE component, involved in application of total hip implant, before and after wear simulation process. The primary results showed that the skew-t random field is more adequate, and flexible for modelling the topographic roughness. Following these arguments, a statistical analysis approach, based on the skew-t random field, is then proposed. It aims at estimating, hierarchically, the significant levels including the real hills/valleys among the uncertain measurements. The evolution of the mean area of the hills/valleys and their levels enabled describing the functional behaviour of the UHMWPE surface over wear time, and indicating the predominant wear mechanisms.

Cette thèse s'intéresse à la conception de méthodes de discrétisation non-conforme pour un modèle de poromécanique. Le but de ce travail est de simplifier les couplages liant la géomécanique d'un milieu poreux à l'écoulement polyphasique compositionnel ayant cours en son sein tels qu'ils sont réalisés actuellement dans l'industrie pétrolière, en discrétisant sur un même maillage, typiquement non-conforme car à l'image de la lithologie, la mécanique et l'écoulement. La nouveauté consiste donc à traiter la mécanique par une méthode d'approximation non-conforme sur maillages généraux. Dans cette thèse, nous nous concentrons sur un modèle d'élasticité linéaire. Les difficultés inhérentes à son approximation non-conforme sont son manque de coercivité (se traduisant par la nécessité de satisfaire une inégalité de Korn sur un espace discret discontinu), ainsi que le phénomène de verrouillage numérique lorsque le matériau tend à devenir incompressible. Dans une première partie, nous construisons un espace d'approximation sur maillages généraux, s'apparentant à une extension de l'espace de Crouzeix-Raviart. Nous explicitons ses propriétés d'approximation et de conformité, et montrons que ce dernier est adapté à une discrétisation primale coercive et robuste au locking du modèle d'élasticité sur maillages généraux. La méthode proposée est moins coûteuse que son équivalent éléments finis (en termes de propriétés) P2. Nous nous intéressons dans une deuxième partie à l'approximation non-conforme d'un modèle couplé de poroélasticité. Nous étudions la convergence d'une famille de schémas numériques dont la discrétisation en espace utilise le formalisme des schémas Gradient, auquel appartient la méthode développée pour la mécanique. Nous prouvons la convergence de telles approximations vers la solution de régularité minimale du problème continu, indépendamment des paramètres physiques du système / This manuscript focuses on the conception of nonconforming discretization methods for a poromechanical model. The aim of this work is to ease the coupling between the geomechanics and the multiphase compositional Darcy flow in porous media by discretizing mechanics and flow on the same mesh, typically nonconforming as it represents the lithology. Hence, the novelty hinges on a nonconforming treatment of mechanics on general meshes. In this work, we focus on a linear elasticity model. The nonconforming approximation of such a model is not straightforward owing to its lack of coercivity (meaning that a discrete Korn's inequality must hold on a discontinuous discrete space) and to the numerical locking phenomenon occurring as the material becomes incompressible. In a first part, we design an approximation space on general meshes, which can be viewed as an extension of the so-called Crouzeix-Raviart space. We study its approximation and conformity properties, and prove that this latter is well-adapted to the design of a primal, coercive, and locking-free discretization of the elasticity model on general meshes. The proposed method is less costly than its finite element equivalent (in terms of properties) P2. In a second part, we tackle the nonconforming approximation of a coupled poroelasticity model. We study the convergence of a family of numerical schemes whose space discretization relies on the Gradient schemes framework, to which belongs the method developed for mechanics. We prove the convergence of such approximations toward the minimal regularity solution of the continuous problem, and independently of the choice of physical parameters

Dans cette thèse, nous avons étudié la simulation numérique des phénomènes d'interaction fluide-structure entre un fluide compressible et une structure déformable. En particulier, nous nous sommes intéressés au couplage par une approche partitionnée entre une méthode de Volumes Finis pour résoudre les équations de la mécanique des fluides compressibles et une méthode d'Eléments discrets pour le solide, capable de prendre en compte la fissuration. La revue des méthodes existantes de domaines fictifs ainsi que des algorithmes partitionnés couramment utilisés pour le couplage conduit à choisir une méthode de frontières immergées conservative et un schéma de couplage explicite. Il est établi que la méthode d'Eléments Discrets utilisée permet de retrouver le comportement macroscopique du matériau et que le schéma symplectique employé assure la préservation de l'énergie du solide. Puis nous avons développé un algorithme de couplage explicite entre un fluide compressible non-visqueux et un solide indéformable. Nous avons montré des propriétés de conservation exacte de masse, de quantité de mouvement et d'énergie du système ainsi que de consistance du schéma de couplage. Cet algorithme a été étendu au couplage avec un solide déformable, sous la forme d'un schéma semi-implicite. Cette méthode a été appliquée à l'étude de problèmes d'écoulements non-visqueux autour de structures mobiles : les comparaisons avec des résultats numériques et expérimentaux existants démontrent la très bonne précision de notre méthode / This work aims at the numerical simulation of compressible fluid/deformable structure interactions. In particular, we have developed a partitioned coupling algorithm between a Finite Volume method for the compressible fluid and a Discrete Element method capable of taking into account fractures in the solid. A survey of existing fictitious domain methods and partitioned algorithms has led to choose an Embedded Boundary method and an explicit coupling scheme. We first showed that the Discrete Element method used for the solid yielded the correct macroscopic behaviour and that the symplectic time-integration scheme ensured the preservation of energy. We then developed an explicit coupling algorithm between a compressible inviscid fluid and an undeformable solid. Mass, momentum and energy conservation and consistency properties were proved for the coupling scheme. The algorithm was then extended to the coupling with a deformable solid, in the form of a semi-implicit scheme. Finally, we applied this method to unsteady inviscid flows around moving structures: comparisons with existing numerical and experimental results demonstrate the excellent accuracy of our method

L’élimination des inclusions indésirables par flottation en cuve agitée est un procédé largement utilisé dans l’industrie de l’aluminium. Ce procédé consiste en l’injection de gaz au sein de la cuve par le biais d’un rotor ; les bulles ainsi formées captent, lors de leur ascension dans le métal liquide, les inclusions et les drainent en surface où elles sont éliminées mécaniquement.Dans le but de prédire l’évolution numérale de la population d’inclusions dans la cuve de flottation au cours du temps, un modèle 2D a été développé. Afin de rendre compte du transport convectif des inclusions au sein de l’écoulement diphasique métal liquide – bulles, des phénomènes de décantation, d’agglomération et de flottation, ce modèle couple la Mécanique des Fluides Numérique et la technique des Bilans de Population. Une attention particulière a été portée aux modèles d’agglomération et de flottation pour qu’ils soient applicables sur une grande plage de conditions de turbulence.Les simulations numériques réalisées, nous ont notamment permis de mettre en évidence l’importance de veiller à une bonne répartition des bulles dans la cuve, en particulier dans les zones de forte turbulence. En parallèle de ce travail de modélisation, une analyse de la littérature et des travaux expérimentaux sur métal liquide ont été réalisés, afin de mieux comprendre le rôle du chlore dans le procédé de flottation. Un effet de seuil a ainsi pu être mis en évidence : une certaine quantité de chlore injectée dans le métal est nécessaire pour que la flottation devienne efficace. La compréhension du rôle du chlore reste cependant partielle, mais des pistes de recherche prometteuses ont put être identifiées / Inclusions removal by flotation in mechanically agitated vessels is widely used in liquid aluminium treatments. This process consists of gas injection in a molten aluminium reactor using an impeller. Inclusions are attached to the bubbles while rising into the liquid bulk and released in the froth at the vessel surface.With the aim of both a better understanding of the physical processes acting during flotation and the optimization of the refining process, a 2D mathematical modelling of the flotation tank was set up. Transport phenomena, sedimentation, agglomeration and flotation of inclusions were considered here. The model couples Population Balance with convective transport of the inclusions, in the multiphase flow field. Simulations showed that a good dispersion of the bubbles into the tank, especially in the high turbulence zones, is essential for an efficient process. This situation was achieved with a strong agitation. Gas flow affected the process efficiency weakly.Theoretical and experimental studies on chlorine effects were done. It was noticed that a critical amount of chlorine into the melt is required for the flotation to be efficient. Complete understanding of chlorine role is not achieved yet, but promising ideas for future research have been exhibited

Ce travail est consacré au développement de schémas numériques pour approcher les solutions de modèles d'écoulement type “depth averaged”. Dans un premier temps, nous détaillons la construction d'approches Volumes Finis pour le système Shallow Water avec termes sources sur maillages non structurés. En se basant sur une reformulation appropriée des équations, nous mettons en place un schéma équilibré et préservant la positivité de la hauteur d'eau, et suggérons des extensions MUSCL adaptées. La méthode est capable de gérer des topographies irrégulières et exhibe de fortes propriétés de stabilité. L'inclusion des termes de friction fait l'objet d'une analyse poussée, aboutissant à l'établissement d'une propriété type “Asymptotic Preserving” à travers l'amélioration d'un autre récent schéma Volumes Finis. La seconde composante de cette étude concerne les méthodes Elements Finis type Galerkin discontinu. Certaines des idées avancées dans le contexte Volumes Finis sont employées pour aborder le système Shallow Water surmaillages triangulaires. Des résultats numériques sont exposés et la méthode se révèle bien adaptée à la description d'une large variété d'écoulements. Partant de ces observations nous proposons finalement d'exploiter ces caractéristiques pour étendre l'approche à une nouvelle famille d'équations type Green-Nadghi. Des validations numériques sont également proposées pour valider le modèle numérique. / This work is devoted to the development of numerical schemes to approximatesolutions of depth averaged flow models. We first detail the construction of Finite Volume approaches for the Shallow Water system with source terms on unstructured meshes. Based on a suitable reformulation of the equations, we implement a well-balanced and positive preserving approach, and suggest adapted MUSCL extensions. The method is shown to handle irregular topography variations and demonstrates strong stabilities properties. The inclusion of friction terms is subject to a thorough analysis, leading to the establishment of some Asymptotic Preserving property through the enhancement of another recent Finite Volume scheme.The second aspect of this study concerns discontinuous Galerkin Finite Elementmethods. Some of the ideas advanced in the Finite Volume context areemployed to broach the Shallow Water system on triangular meshes. Numericalresults are exposed and the method turns out to be well suited to describe a large variety of flows. On these observations we finally propose to exploit its features to extend the approach to a new family of Green-Nadghi equations. Numerical experiments are also proposed to validate this numerical model.

[pt] Neste trabalho apresentamos um algoritmo baseado em agrupamento de dados para o mapeamento automático de horizontes e de falhas sísmicas a partir de dados sísmicos 3D. Apresentamos uma técnica para quantizar o volume sísmico de entrada a partir dos neurônios do grafo resultante do processo de treinamento de uma instância do algoritmo Growing Neural Gas (GNG). No conjunto de amostras de entrada utilizadas pelo GNG, cada amostra representa um voxel do volume de entrada, e retém informações da vizinhança vertical desse voxel. Depois da etapa de treinamento, a partir do grafo gerado pelo GNG um novo volume quantizado é gerado, e nesse volume possíveis ambigüidades e imperfeições existentes no volume de entrada tendem a ser minimizadas. A partir do volume quantizado descrevemos uma nova técnica de extração de horizontes, desenvolvida com o objetivo de que seja possível mapear horizontes na presença de estruturas geológicas complexas, como por exemplo horizontes que possuam porções completamente desconectadas por uma ou mesmo diversas falhas sísmicas. Também iniciamos o desenvolvimento de uma abordagem de mapeamento de falhas sísmicas utilizando informações presentes no volume quantizado. Os resultados obtidos pelo processo de mapeamento de horizontes, testado em volumes diferentes, foram bastante promissores. Além disso, os resultados iniciais obtidos pelo processo de extração de falhas sugerem que a técnica pode vir a ser uma boa alternativa para a tarefa. / [en] In this work we present a clusterization-based method to map seismic horizons and faults from 3D seismic data. We describe a method used to quantize an initial seismic volume using a trained instance of the Growing Neural Gas (GNG) algorithm. To accomplish this task we create a training set where each sample corresponds to an entry volume voxel, retaining its vertical neighboring information. After the training procedure, the resulting graph is used to create a quantized version of the original volume. In this quantized volume both horizons and faults are more evidenced in the data, and we present a method that uses the created volume to map seismic horizons, even when they are completely disconnected by seismic faults. We also present another method that uses the quantized version of the volume to map the seismic faults. The horizon mapping procedure, tested in different volume date, yields good results. The preliminary results presented for the fault mapping procedure also yield good results, but needs further testing.

Cette thèse est centrée autour du développement et de l'analyse des schémas volumes finis robustes afin d'approcher les solutions du modèle diphasique compressible en milieux poreux hétérogènes et anisotropes. Le modèle à deux phases compressibles comprend deux équations paraboliques dégénérées et couplées dont les variables principales sont la saturation du gaz et la pression globale. Ce système est discrétisé à l'aide de deux méthodes différentes (CVFE et DDFV) qui font partie de la famille des volumes finis. La première classe à laquelle on s'intéresse consiste à combiner la méthode des volumes finis et celle des éléments finis. Dans un premier temps, on considère un schéma volume finis upwind pour la partie convective et un schéma de type éléments finis conformes pour la diffusion capillaire. Sous l'hypothèse que les coefficients de transmissibilités sont positifs, on montre que la saturation vérifie le principe du maximum et on établit des estimations d'énergies permettant de démontrer la convergence du schéma. Dans un second temps, on a mis en place un schéma positif qui corrige le précédent. Ce schéma est basé sur une approximation des flux diffusifs par le schéma de Godunov. L'avantage est d'établir la bornitude des solutions approchées ainsi que les estimations uniformes sur les gradients discrets sans aucune contrainte ni sur le maillage ni sur la perméabilité. En utilisant des arguments classiques de compacité, on prouve rigoureusement la converge du schéma. Chaque schéma est validé par des simulations numériques qui montrent bien le comportement attendu d'une telle solution. Concernant la deuxième classe, on s'intéressera tout d'abord à la construction et à l'étude d'un nouveau schéma de type DDFV (Discrete Duality Finite Volume) pour une équation de diffusion non linéaire dégénérée. Cette méthode permet d' avantage de prendre en compte des maillages très généraux et des perméabilités quelconques. L'idée clé de cette discrétisation est d'approcher les flux dans la direction normale par un schéma centré et d'utiliser un schéma décentré dans la direction tangentielle. Par conséquent, on démontre que la solution approchée respecte les bornes physiques et on établit aussi des estimations d'énergie. La convergence du schéma est également établie. Des résultats numériques confirment bien ceux de la théorie. Ils exhibent en outre que la méthode est presque d'ordre deux. / The objective of this thesis is the development and the analysis of robust and consistent numerical schemes for the approximation of compressible two-phase flow models in anisotropic and heterogeneous porous media. A particular emphasis is set on the anisotropy together with the geometric complexity of the medium. The mathematical problem is given in a system of two degenerate and coupled parabolic equations whose main variables are the nonwetting saturation and the global pressure. In view of the difficulties manifested in the considered system, its cornerstone equations are approximated with two different classes of the finite volume family. The first class consists of combining finite elements and finite volumes. Based on standard assumptions on the space discretization and on the permeability tensor, a rigorous convergence analysis of the scheme is carried out thanks to classical arguments. To dispense with the underlined assumptions on the anisotropy ratio and on the mesh, the model has to be first formulated in the factional flux formulation. Moreover, the diffusive term is discretized by a Godunov-like scheme while the convective fluxes are approximated using an upwind technique. The resulting scheme preserves the physical ranges of the computed solution and satisfies the coercivity property. Hence, the convergence investigation holds. Numerical results show a satisfactory qualitative behavior of the scheme even if the medium of interest is anisotropic. The second class allows to consider more general meshes and tensors. It is about a new positive nonlinear discrete duality finite volume method. The main point is to approximate a part of the fluxes using a non standard technique. The application of this ideato a nonlinear diffusion equation yields surprising results. Indeed,not only is the discrete maximum property fulfilled but also the convergence of the scheme is established. Practically, the proposed method shows great promises since it provides a positivity-preserving and convergent scheme with optimal convergence rates.

O fenômeno de injeção de vapor em tanques de alívio é de relevância no projeto de reatores de água leve, sejam eles do tipo reator de água pressurizada (PWR) ou reator de água fervente (BWR). Este fenômeno permite a rápida absorção do vapor injetado em massa de água, por meio de sua condensação, uma vez que este vapor pode conter contaminantes químicos ou radiológicos que não permitem o seu descarte diretamente no ambiente. Desta forma, facilita-se a coleta do vapor produzido por descarga de vapor da água do resfriamento do reator, radiologicamente contaminada, e evita-se o que projeto de dispositivos e equipamentos necessite considerar a elevada pressão do vapor. A rapidez com que se dá a condensação é fruto de processos físicos que ocorrem na interface de vapor e água e que ainda não possuem modelo analítico e numérico definido. Em 1972 um modelo semi-empírico foi proposto, o qual, desde então, vem evoluindo. Não obstante, até o presente momento, não há modelo definitivo que se proponha a abranger toda extensão das condições experimentais. Estes modelos são fortemente dependentes do fluxo de massa que atravessa a interface de vapor e água, entretanto, até a presente data, não há expressão que determine este fluxo de massa, de tal forma que o valor de 275 Kg/m2/s vem sendo assumido como \"representativo da ordem de grandeza do fenômeno\" até o presente momento. Neste trabalho, é proposto um método de cálculo analítico do fluxo de massa, considerando-se como premissa a isentropia da injeção, e o desenvolvimento da 1ª e 2ª leis da Termodinâmica. Ainda, o fenômeno é analisado experimentalmente, por meio da análise dos dados produzidos no experimento do Circuito Termo Hidráulico de 150 bar (Loop 150), realizado nas dependências do CENTRO TECNOLÓGICO DA MARINHA EM SÃO PAULO. Por fim, um modelo numérico em software comercial foi desenvolvido para complementar a análise. Os resultados obtidos comprovam que a formulação isentrópica do fluxo de massa corrige de maneira satisfatória o fluxo de massa constante utilizado até então nos modelos semi-empíricos. Tal comprovação se deu através de análise numérica e da confrontação com dados experimentais obtidos na literatura. / The phenomenon of vapor injection in relief tanks presents relevance in the design of light water reactors, be they of the type pressurized water reactor (PWR) or boiling water reactor (BWR). This phenomenon allows the rapid absorption of the vapor injected in a mass of water, by condensation. Since this vapor may contain chemical or radiological contaminants that do not allow its discharge directly in the environment, it must be collected. The condensation avoids the design of devices and equipment, which need to consider the high vapor pressure, and allows the vapor to be collected. The rapidity with which the condensation occurs is the result of physical processes that occur at the interface of steam and water. These processes do not yet have a defined numerical and analytical model. In 1972 a semi-empirical model was proposed, which has, since then, evolved. Nevertheless, up to the present moment, there is no definitive model that intends to cover every extension of the experimental conditions. These models are strongly dependent on the mass flow through the steam and water interface, however, up to date, there is no expression that determines this mass flow. For the sake of this, the value of 275 kg / m2 / s has been assumed as \"representative of the order of magnitude of the phenomenon\" up to the present moment. In this work, a method of analytical calculation of mass flow is proposed, considering as premise the isotropy of the injection, and the development of the 1st and 2nd laws of thermodynamics. Still, the phenomenon is analyzed experimentally, by means of the analysis of the data produced in the experiment of Loop 150, realized in dependencies of the CENTRO TECNOLÓGICO DA MARINHA EM SÃO PAULO. Finally, a numerical model in commercial software was developed to complement the analysis. The result is proven with an isentropic mass flow formulation, which satisfactorily corrected the mass used in the former semi-empirical models. Such verification was performed through a series of data and confrontation with experimental data in the literature.