Abonnement à la biblothèque: Guest
Portail numérique Bibliothèque numérique eBooks Revues Références et comptes rendus Collections
Computational Thermal Sciences: An International Journal
ESCI SJR: 0.249 SNIP: 0.434 CiteScore™: 1.4

ISSN Imprimer: 1940-2503
ISSN En ligne: 1940-2554

Computational Thermal Sciences: An International Journal

DOI: 10.1615/ComputThermalScien.2015014329
Forthcoming Article

Simulations of two-phase flows with a multifield approach

Solene Fleau
EDF R&D, Université Paris-Est
Stéphane Mimouni
Nicolas Mérigoux
Stéphane Vincent
Laboratoire de Modélisations et Simulations Multi-Echelle, Université Paris-Est, France


Safety issues in nuclear power plant involve complex bubbly flows. To predict the behavior of these flows, the two-fluid approach is often used. Nevertheless, this model induces a numerical diffusion of interfaces, which results in a poor accuracy in the calculation of the local parameters. Therefore, to simulate large interfaces such as slugs or free surfaces, located methods have been developed using the single-fluid model. In this paper, the two approaches have been coupled in the CMFD code NEPTUNE_CFD to simulate adiabatic separated flows. The averaged momentum balance equations are solved for each field and are followed by an artificial compression step, which fixes the interface thickness and ensures mass conservation. Moreover, since the two-fluid model allows the existence of relative velocities at the interface, a drag force is used to cancel them. This article proposes also a new formulation for this force, to take into account the physical properties of the flow. To validate this approach, an analytical test case with a static bubble has been simulated with a mesh refinement test. Then, the simulation of the Kelvin-Helmholtz instability has been performed to highlight the effect of the modification of the drag force. The sheared interface is particularly sensible to this force, which has an important influence on the flow parameters such as the interface velocity. Finally, these developments have been compared to three other codes by simulating the Rayleigh-Taylor instability.