skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A Ghost Fluid/Level Set Method for boiling flows and liquid evaporation: Application to the Leidenfrost effect

Abstract

The development of numerical methods for the direct numerical simulation of two-phase flows with phase change, in the framework of interface capturing or interface tracking methods, is the main topic of this study. We propose a novel numerical method, which allows dealing with both evaporation and boiling at the interface between a liquid and a gas. Indeed, in some specific situations involving very heterogeneous thermodynamic conditions at the interface, the distinction between boiling and evaporation is not always possible. For instance, it can occur for a Leidenfrost droplet; a water drop levitating above a hot plate whose temperature is much higher than the boiling temperature. In this case, boiling occurs in the film of saturated vapor which is entrapped between the bottom of the drop and the plate, whereas the top of the water droplet evaporates in contact of ambient air. The situation can also be ambiguous for a superheated droplet or at the contact line between a liquid and a hot wall whose temperature is higher than the saturation temperature of the liquid. In these situations, the interface temperature can locally reach the saturation temperature (boiling point), for instance near a contact line, and be cooler in other places.more » Thus, boiling and evaporation can occur simultaneously on different regions of the same liquid interface or occur successively at different times of the history of an evaporating droplet. Standard numerical methods are not able to perform computations in these transient regimes, therefore, we propose in this paper a novel numerical method to achieve this challenging task. Finally, we present several accuracy validations against theoretical solutions and experimental results to strengthen the relevance of this new method.« less

Authors:
; ; ;
Publication Date:
OSTI Identifier:
22572332
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Computational Physics; Journal Volume: 316; Other Information: Copyright (c) 2016 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; ACCURACY; BOILING; BOILING POINTS; CALCULATION METHODS; CAPTURE; COMPUTERIZED SIMULATION; DROPLETS; EVAPORATION; FILMS; HEAT EXCHANGERS; INTERFACES; LIQUIDS; SATURATION; TRANSIENTS; TWO-PHASE FLOW; VALIDATION; VAPORS; WALLS; WATER

Citation Formats

Rueda Villegas, Lucia, Alis, Romain, Lepilliez, Mathieu, and Tanguy, Sébastien, E-mail: tanguy@imft.fr. A Ghost Fluid/Level Set Method for boiling flows and liquid evaporation: Application to the Leidenfrost effect. United States: N. p., 2016. Web. doi:10.1016/J.JCP.2016.04.031.
Rueda Villegas, Lucia, Alis, Romain, Lepilliez, Mathieu, & Tanguy, Sébastien, E-mail: tanguy@imft.fr. A Ghost Fluid/Level Set Method for boiling flows and liquid evaporation: Application to the Leidenfrost effect. United States. doi:10.1016/J.JCP.2016.04.031.
Rueda Villegas, Lucia, Alis, Romain, Lepilliez, Mathieu, and Tanguy, Sébastien, E-mail: tanguy@imft.fr. 2016. "A Ghost Fluid/Level Set Method for boiling flows and liquid evaporation: Application to the Leidenfrost effect". United States. doi:10.1016/J.JCP.2016.04.031.
@article{osti_22572332,
title = {A Ghost Fluid/Level Set Method for boiling flows and liquid evaporation: Application to the Leidenfrost effect},
author = {Rueda Villegas, Lucia and Alis, Romain and Lepilliez, Mathieu and Tanguy, Sébastien, E-mail: tanguy@imft.fr},
abstractNote = {The development of numerical methods for the direct numerical simulation of two-phase flows with phase change, in the framework of interface capturing or interface tracking methods, is the main topic of this study. We propose a novel numerical method, which allows dealing with both evaporation and boiling at the interface between a liquid and a gas. Indeed, in some specific situations involving very heterogeneous thermodynamic conditions at the interface, the distinction between boiling and evaporation is not always possible. For instance, it can occur for a Leidenfrost droplet; a water drop levitating above a hot plate whose temperature is much higher than the boiling temperature. In this case, boiling occurs in the film of saturated vapor which is entrapped between the bottom of the drop and the plate, whereas the top of the water droplet evaporates in contact of ambient air. The situation can also be ambiguous for a superheated droplet or at the contact line between a liquid and a hot wall whose temperature is higher than the saturation temperature of the liquid. In these situations, the interface temperature can locally reach the saturation temperature (boiling point), for instance near a contact line, and be cooler in other places. Thus, boiling and evaporation can occur simultaneously on different regions of the same liquid interface or occur successively at different times of the history of an evaporating droplet. Standard numerical methods are not able to perform computations in these transient regimes, therefore, we propose in this paper a novel numerical method to achieve this challenging task. Finally, we present several accuracy validations against theoretical solutions and experimental results to strengthen the relevance of this new method.},
doi = {10.1016/J.JCP.2016.04.031},
journal = {Journal of Computational Physics},
number = ,
volume = 316,
place = {United States},
year = 2016,
month = 7
}
  • While Eulerian schemes work well for most gas flows, they have been shown to admit nonphysical oscillations near some material interfaces. In contrast, Lagrangian schemes work well at multimaterial interfaces, but suffer from their own difficulties in problems with large deformations and vorticity characteristic of most gas flows. The authors believe that the most robust schemes will combine the best properties of Eulerian and Lagrangian schemes. In this paper, they propose a new numerical method for treating interfaces in Eulerian schemes that maintains a Heaviside profile of the density with no numerical smearing along the lines of earlier work andmore » most Lagrangian schemes. They use a level set function to track the motion of a multimaterial interface in an Eulerian framework. In addition, the use of ghost cells (actually ghost nodes in the finite difference framework) and a new isobaric fix technique allows them to keep the density profile from smearing out, while still keeping the scheme robust and easy to program with simple extensions to multidimensions and multilevel time integration, e.g., Runge-Kutta methods. In contrast, previous methods used ill-advised dimensional splitting for multidimensional problems and suffered from great complexity when used in conjunction with multilevel time integrators.« less
  • This paper presents a novel methodology for simulating incompressible two-phase flows by combining an improved version of the conservative level set technique introduced in [E. Olsson, G. Kreiss, A conservative level set method for two phase flow, J. Comput. Phys. 210 (2005) 225-246] with a ghost fluid approach. By employing a hyperbolic tangent level set function that is transported and re-initialized using fully conservative numerical schemes, mass conservation issues that are known to affect level set methods are greatly reduced. In order to improve the accuracy of the conservative level set method, high order numerical schemes are used. The overallmore » robustness of the numerical approach is increased by computing the interface normals from a signed distance function reconstructed from the hyperbolic tangent level set by a fast marching method. The convergence of the curvature calculation is ensured by using a least squares reconstruction. The ghost fluid technique provides a way of handling the interfacial forces and large density jumps associated with two-phase flows with good accuracy, while avoiding artificial spreading of the interface. Since the proposed approach relies on partial differential equations, its implementation is straightforward in all coordinate systems, and it benefits from high parallel efficiency. The robustness and efficiency of the approach is further improved by using implicit schemes for the interface transport and re-initialization equations, as well as for the momentum solver. The performance of the method is assessed through both classical level set transport tests and simple two-phase flow examples including topology changes. It is then applied to simulate turbulent atomization of a liquid Diesel jet at Re=3000. The conservation errors associated with the accurate conservative level set technique are shown to remain small even for this complex case.« less
  • Abstract not provided.