On two-phase flow solvers in irregular domains with contact line

We present numerical methods that enable the direct numerical simulation of two-phase flows in irregular domains. A method is presented to account for surface tension effects in a mesh cell containing a triple line between the liquid, gas and solid phases. Our numerical method is based on the level-set method to capture the liquid–gas interface and on the single-phase Navier–Stokes solver in irregular domain proposed in [35]to impose the solid boundary in an Eulerian framework. We also present a strategy for the implicit treatment of the viscous term and how to impose both a Neumann boundary condition and a jump condition when solving for the pressure field. Special care is given on how to take into account the contact angle, the no-slip boundary condition for the velocity field and the volume forces. Finally, we present numerical results in two and three spatial dimensions evaluating our simulations with several benchmarks

An efficient mass-preserving interface-correction level set/ghost fluid method for droplet suspensions under depletion forces

Aiming for the simulation of colloidal droplets in microfluidic devices, we present here a numerical method for two-fluid systems subject to surface tension and depletion forces among the suspended droplets. The algorithm is based on an efficient solver for the incompressible two-phase Navier–Stokes equations, and uses a mass-conserving level set method to capture the fluid interface. The four novel ingredients proposed here are, firstly, an interface-correction level set (ICLS) method; global mass conservation is achieved by performing an additional advection near the interface, with a correction velocity obtained by locally solving an algebraic equation, which is easy to implement in both 2D and 3D. Secondly, we report a second-order accurate geometric estimation of the curvature at the interface and, thirdly, the combination of the ghost fluid method with the fast pressurecorrection approach enabling an accurate and fast computation even for large density contrasts. Finally, we derive a hydrodynamic model for the interaction forces induced by depletion of surfactant micelles and combine it with a multiple level set approach to study short-range interactions among droplets in the presence of attracting forces

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

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

Numerical Study of Interfacial Flow using Algebraic Coupled Level Set-Volume of Fluid (A-CLSVOF) Method

Solving interfacial flows numerically has been a challenge due to the lack of sharpness and the presence of spurious currents at the interface. Two methods, Algebraic Coupled Level Set-Volume of Fluid (A-CLSVOF) method and Ghost Fluid Method (GFM) have been developed in the finite volume framework and employed in several interfacial flows such as Rayleigh-Taylor instability, rising bubble, impinging droplet and cross-flow oil plume. In the static droplet simulation, A-CLSVOF substantially reduces the spurious currents. The capillary wave relaxation shows that this method delivers results comparable to those of more rigorous methods such as Front Tracking methods for fine grids. The results for the other interfacial flows also compared well with the experimental results. Next, interfacial forces are implemented by enlisting the finite volume discretization of Ghost Fluid Method. To assess the A-CLSVOF/GFM performance, four cases are studied. In the case of the static droplet in suspension, the combined A-CLSVOF/GFM produces a sharp and accurate pressure jump compared to the traditional CSF (continuum surface force) implementation. For the linear two-layer shear flow, GFM sharp treatment of the viscosity captured the velocity gradient across the interface. For a gaseous bubble rising in a viscous fluid, GFM outperforms CSF by almost 10%. Also, a Decoupled Pressure A-CLSVOF/GFM method (DPM) has been developed which separates pressure into two pressure components, one accounting for interfacial forces such as surface tension and another representing the rest of flow pressure. It is proven that the DPM implementation results in more efficiency in PISO (Pressure Implicit with Splitting of Operators) loop. A two-phase solver is used to study buoyant oil discharge in quiescent and cross-flow ambient. Different modes of breakup including dripping, jetting (axisymmetric and asymmetric) and atomization for cross-flow oil jet are captured

Fluids-membrane interaction with a full Eulerian approach based on the level set method

A fully Eulerian approach to predict fluids-membrane behaviours is presented in this paper. Based on the numerical model proposed by Ii et al. (2012), we present a sharp methodology to account for the jump conditions due to hyperelastic membranes. The membrane is considered infinitely thin and is represented by the level set method. Its deformations are obtained from the transport of the components of the left Cauchy-Green tensor throughout time. Considering the linear or a hyperelastic material law, the surface stress tensor is computed and gives the force exerted by the membrane on the surrounding fluids. The membrane force is taken into account in the Navier-Stokes equations as jump conditions on the pressure and on the velocity derivatives by imposing suitable singular source terms in cells crossed by the interface. To prevent stability issues, an extension algorithm has been developed to remove the normal derivatives of the scalar fields specific to the membrane. In particular, a subcell resolution at the interface of the extrapolated variable is proposed for increasing the accuracy of the extension algorithm. These improvements are validated by comparing our numerical results with benchmarks from the literature. Moreover, a new benchmark is proposed for fluids with both different viscosities and different densities to target applications where a gas and a liquid phase are separated by a membrane

Convection in the Melt

A physical problem involving the melting/freezing of a phase-change material (PCM) is the applied setting of this research. The development of models that couple the partial differential equations for energy transport and fluid motion with phases of differing densities is a primary goal of the research. In Chapter 2, a general framework is developed for the formulation of conservation laws that admit interfaces. A notion of weak solution is developed and its relation with classical and other weak formulations is discussed. Conditions that hold across various kinds of interfaces are also developed. The formulation is examined for the conservation of mass, momentum and energy in Chapter 3. In Chapter 4, a numerical method for the solution of conservation law equations is given. The method uses a Crank-Nicolson time discretization and solves the implicit equations with a Newton/Approximate Factorization technique. The method captures interfaces and is consistent with the control volume weak formulations of Chapter 2. The numerical solution converges to the distributional solution of the conservation law. In Chapter 5, three applications of the theory are developed and numerical computations are presented. First, a one dimensional problem is studied involving conservation of mass. momentum and energy in a phase-change material with a liquid density larger than that of the solid. The second application is a suction problem in two dimensions. The bulk movement of a liquid and void are simulated with and without the effects of surface tension. The third application is to a three-dimensional simulation of the heating of a cylindrical canister of PCM in 1-g and 0-g. For this simulation the Marangoni stress is the important driving force on the flow