Numerical Simulation Of Two-phase Flow Research Articles

Numerical analysis of tsunami-triggered oil spill from industrial parks in Osaka Bay

Oil storage tanks in industrial parks were heavily damaged by a large-scale tsunami, and large-scale oil spill occurred, which led to a fire. Moreover, it was found that mud samples in these regions had a high level of oil contamination after the disaster, showing that the tsunami triggered high turbid water mixed with spill oil and deposited on the sea bottom. For industrial zones around the coastal line with the potential for tsunamis, the building of a scenario of a tsunami-triggered oil spill from these industrial parks is urgently needed for planning the ship evacuating routes and the evacuation of residences nearby. This was carried out by numerical simulations of two-phase flow around the industrial parks. The results could build a scenario of a tsunami-triggered oil spill from industrial zones and could help to add the oil spill effect in reviewing the risk management of industrial zones around the coastal line.

A weighted multiple-relaxation-time lattice Boltzmann method for multiphase flows and its application to partial coalescence cascades

We present a lattice Boltzmann method (LBM) with a weighted multiple-relaxation-time (WMRT) collision model and an adaptive mesh refinement (AMR) algorithm for direct numerical simulation of two-phase flows in three dimensions. The proposed WMRT model enhances the numerical stability of the LBM for immiscible fluids at high density ratios, particularly on the D3Q27 lattice. The effectiveness and efficiency of the proposed WMRT–LBM–AMR is validated through simulations of (a) buoyancy-driven motion and deformation of a gas bubble rising in a viscous liquid; (b) the bag-breakup mechanism of a falling drop; (c) crown splashing of a droplet on a wet surface; and (d) the partial coalescence mechanism of a liquid drop at a liquid–liquid interface. The numerical simulations agree well with available experimental data and theoretical approximations where applicable.

Simulation of two-phase flow and mass transfer with deformable interface

Two-phase flow and interphase mass transfer are very popular in chemical processes. Rapid development of numerical simulation plays a crucial role in deep understanding of the mechanism of momentum, mass and heat transfer in two-phase systems. However, establishing and solving the mathematical models of two-phase flow and transport processes, especially the accurate description of deformable moving interface, are critical and difficult. This review addresses the level set method and the boundary element method, which are involved with the numerical simulation of two-phase flow of deformable surfaces, and presents a detailed discussion of the limits and advantages of these typical methods. We show the research status and prospect of the mathematical methods of deformation interface of fluid particle coalescence and breakage, micro-fluidics and droplets in two-phase flow field. In the end, we briefly introduce the numerical simulation methods of interphase mass transfer across deformable interface.

A new approach to solve mixture Multi-phase Flow model using Time Splitting Projection Method

This paper was aimed at introducing a new approach to numerical simulation of two-phase flows by solving the mixture model equations. In this approach, an effective numerical algorithm was developed based on the time splitting projection method to solve the mixture model equations for an incompressible two-phase flow. In this study, a new technique was developed to solve the mixture continuity equation based on the time splitting projection method. One of the advantages of the mixture model is that, in addition to determining the velocities of each phase, it solves only one set of momentum equations and calculates the velocity differences between the phases by solving the relative velocity equation. Accordingly, an extended finite volume vertical two-dimensional numerical model was used to determine the pressure distribution, velocity field and volume fraction of each phase at each time step. The model has been used in various tests to simulate unsteady flow problems. Comparison between numerical results, analytical solutions and experimental data demonstrated a satisfactory performance.

Physical and numerical simulation of two-phase flow in the area of mini-channel

The results of physical and numerical modeling to determine the rate of evaporation of ethanol, at a flow of air flow velocities in the range of 0.0139 - 0.139 m/s and temperature in the range of 10 - 45 °C are placed. The time of the establishment of the diffusion equilibrium in the work area and differences in the results of numerical and physical modeling are have determined.

Applying Different Strategies within OpenFOAM to Investigate the Effects of Breakup and Collision Model on the Spray and in-Cylinder Gas Mixture Attribute

In the current study, a 3-D numerical simulation of two-phase flow has been conducted in a direct injection CI engine using the Eularian-Lagrangian approach and a new breakup model. The newly modified breakup scheme has been implemented for simulating the ultra-high pressure diesel injection. The effects of droplet breakup and collision model on the spray and in-cylinder gas characteristics have been examined using the open source code OpenFOAM. Spray penetration and cone angle are investigated as spray properties and surrounding gas motion are studied by in-cylinder gas velocity and pressure distribution for non-evaporating conditions. In addition, vapor penetration of the evaporating spray is presented to study the effects of current scheme on the evaporating condition. The continuous field is described by RANS equations and dynamics of the dispersed droplet is modeled by Lagrangian tracking scheme. Results of the proposed modified KHRT model are compared against other default methods in OpenFOAM and favorable agreement is achieved. Robustness and accuracy of different breakup schemes and collision models are also verified using the published experimental data. It is demonstrated that the proposed breakup scheme and Nordin collision model display very accurate results in the case of ultra-high pressure injection.

Numerical simulation of two-phase flow in fractured porous media using streamline simulation and IMPES methods and comparing results with a commercial software

Streamline simulation is developed to simulate waterflooding in fractured reservoirs. Conventional reservoir simulation methods for fluid flow simulation in large and complex reservoirs are very costly and time consuming. In streamline method, transport equations are solved on one-dimensional streamlines to reduce the computation time with less memory for simulation. First, pressure equation is solved on an Eulerian grid and streamlines are traced. Defining the “time of flight”, saturation equations are mapped and solved on streamlines. Finally, the results are mapped back on Eulerian grid and the process is repeated until the simulation end time. The waterflooding process is considered in a fractured reservoir using the dual porosity model. Afterwards, a computational code is developed to solve the same problem by the IMPES method and the results of streamline simulation are compared to those of the IMPES and a commercial software. Finally, the accuracy and efficiency of streamline simulator for simulation of two-phase flow in fractured reservoirs has been proved.

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.

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.

Numerical study on the wettability dependent interaction of a rising bubble with a periodic open cellular structure

A phase-field method for interface resolving numerical simulations of two-phase flows with OpenFOAM® is validated for the buoyancy-driven rise of a single air bubble through a viscous stagnant liquid using experimental data from literature. The validation encompasses the terminal bubble rise velocity and the instantaneous cutting of the bubble by a solid horizontal cylinder. In the latter process, the numerical method takes into account the equilibrium contact angle of the three-phase system. The numerical method is then used to study the behavior of a single air bubble rising through a representative subdomain of a periodic open cellular structure (POCS) with cubic cell geometry filled with stagnant water. The results indicate that the bubble shape and path do significantly depend on the structure wettability. In the industrial application of POCS for enhancing mass transfer and as catalytic supports, the utilization of structures with high wettability (low contact angles) is expected to be beneficial.

A mass-conserving lattice Boltzmann method with dynamic grid refinement for immiscible two-phase flows

A mass-conserving lattice Boltzmann method (LBM) for multiphase flows is presented in this paper. The proposed LBM improves a previous model (Lee and Liu, 2010 [21] ) in terms of mass conservation, speed-up, and efficiency, and also extends its capabilities for implementation on non-uniform grids. The presented model consists of a phase-field lattice Boltzmann equation (LBE) for tracking the interface between different fluids and a pressure-evolution LBM for recovering the hydrodynamic properties. In addition to the mass conservation property and the simplicity of the algorithm, the advantages of the current phase-field LBE are that it is an order of magnitude faster than the previous interface tracking LBE proposed by Lee and Liu (2010) [21] and it requires less memory resources for data storage. Meanwhile, the pressure-evolution LBM is equipped with a multi-relaxation-time (MRT) collision operator to facilitate attainability of small relaxation rates thereby allowing simulation of multiphase flows at higher Reynolds numbers. Additionally, we reformulate the presented MRT-LBM on nonuniform grids within an adaptive mesh refinement (AMR) framework. Various benchmark studies such as a rising bubble and a falling drop under buoyancy, droplet splashing on a wet surface, and droplet coalescence onto a fluid interface are conducted to examine the accuracy and versatility of the proposed AMR-LBM. The proposed model is further validated by comparing the results with other LB models on uniform grids. A factor of about 20 in savings of computational resources is achieved by using the proposed AMR-LBM. As a more demanding application, the Kelvin–Helmholtz instability (KHI) of a shear-layer flow is investigated for both density-matched and density-stratified binary fluids. The KHI results of the density-matched fluids are shown to be in good agreement with the benchmark AMR results based on the sharp-interface approach. When a density contrast between the two fluids exists, a typical chaotic structure in the flow field is observed at a Reynolds number of 10000, which indicates that the proposed model is a promising tool for direct numerical simulation of two-phase flows.

One-dimensional comparison of numerical approaches on two-phase flow in the membrane electrode assembly of PEMFC

Numerical simulation of two-phase flow is important for the performance evaluation of proton exchange membrane fuel cell (PEMFC). Exchange of water between vapor and liquid phases in the gas diffusion layer (GDL) was considered by Nguyen and co-workers (2000) using the evaporation and condensation rate constants. Hsuen and Yin (2011), in which a pseudo-phase-equilibrium function approximates the vapor/liquid phase equilibrium, proposed another two-phase flow approach. Both methods are capable of predicting the cell performance and the species concentration profiles in the membrane electrode assembly (MEA) without the need to track liquid water front in the GDL explicitly. The present study compares the numerical results of these two approaches in detail. In addition, one-dimensional simplified semi-analytical solution based on the instantaneous phase equilibrium assumption is used to validate the numerical calculations. The evaporation/condensation approach allows the existence of super-saturated vapor near the interface of cathode GDL and membrane. Thus, membrane can be better hydrated near the interface once liquid water emerges in the GDL as compared to the pseudo-phase equilibrium method. Consequently, smaller ohmic resistance and hence better cell performance are observed with Nguyen's prediction once liquid water appears in the GDL.

Choice of measure source terms in interface coupling for a model problem in gas dynamics

This paper is devoted to the mathematical and numerical analysis of a coupling procedure for one-dimensional Euler systems. The two systems have different closure laws and are coupled through a thin fixed interface. Following the work of [5], we propose to couple these systems by a bounded vector-valued Dirac measure, concentrated at the coupling interface, which in the applications may have a physical meaning. We show that the proposed framework allows to control the coupling conditions and we propose an approximate Riemann solver based on a relaxation approach preserving equilibrium solutions of the coupled problem. Numerical experiments in constrained optimization problems are then presented to assess the performances of the present method. 1. Introduction The study of large-scale and complex problems exhibiting a wide range of physical space and time scales (see for instance [62, 35, 14]), usually requires separate solvers adapted to the resolution of specific scales. This is the case of many industrial flows. Let us quote, for example, the numerical simulation of two-phase flows applied to the burning liquid oxygen-hydrogen gas in rocket engines [58]. This kind of flow contains both separated and dispersed two-phase flows, due to atomization and evaporation phenomena. This requires appropriate models and solvers for separated and dispersed phases that have to be appropriately coupled. Another example concerns turbomachine flows which can be modeled by the Euler equations of gas dynamics with different closure laws between the stages of the turbine, where the conditions of temperature and pressure are strongly heterogeneous. The coupling of these different systems is thus necessary to give a complete description of the flow inside the whole turbine. The method of interface coupling allows to represent the evolution of such flows, where different models are separated by fixed interfaces. First, coupling conditions are specified at the interface to exchange information between the systems. The definition of transmission conditions generally results from physical consideration, e.g. the conservation or the continuity of given variables. Then, the transmission conditions are represented at the discrete level. The study of interface coupling for nonlinear hyperbolic systems has received attention for several years. In [43], the authors study the scalar case from both mathematical and numerical points of view.

New Leakage Current Particulate Matter Sensor for On-Board Diagnostics

Structure and principle of the new leakage current particulate matter (PM) sensor are introduced and further study is performed on the PM sensor with the combination of numerical simulation and bench test. High voltage electrode, conductive shell, and heaters are all built-in. Based on the principle of Venturi tube and maze structure design, this sensor can detect transient PM concentrations. Internal flow field of the sensor and distribution condition of PM inside the sensor are analyzed through gas-solid two-phase flow numerical simulation. The experiment was also carried out on the whole sensor system (including mechanical and electronic circuit part) and the output signals were analyzed. The results of simulation and experiment reveal the possibility of PM concentration (mass) detection by the sensor.

Numerical simulation of two-phase flows in complex geometries by using the volume-of-fluid/immersed-boundary method

The numerical modeling of two-phase flows in complex geometries is of special relevance in many scientific and industrial applications. In this paper, we describe a numerical method to perform three-dimensional simulations of such flow problems. This method adopts a volume-of-fluid (VOF) approach to capture and advance the fluid interface, and it integrates the fluid solver with the immersed boundary (IB) modeling of arbitrary-shape walls and moving bodies. The shape and movement of solid geometries are efficiently represented by an auxiliary signed distance function (SDF) with local coordinate transformation. Validation tests have been conducted using the present method, and the computational results are in good agreements with reference solutions and experimental data. We further present its application to the simulation of the air–water flow in a twin screw kneader. Hence, the adequacy and suitability of the present VOF–IB method are shown to successfully simulate complicated two-phase flows interacting with general geometries.

A stable and linear time discretization for a thermodynamically consistent model for two-phase incompressible flow

A new time discretization scheme for the numerical simulation of two-phase flow governed by a thermodynamically consistent diffuse interface model is presented. The scheme is consistent in the sense that it allows for a discrete in time energy inequality. An adaptive spatial discretization is proposed that conserves the energy inequality in the fully discrete setting by applying a suitable post processing step to the adaptive cycle. For the fully discrete scheme a quasi-reliable error estimator is derived which estimates the error both of the flow velocity, and of the phase field. The validity of the energy inequality in the fully discrete setting is numerically investigated.

Pressure and volume changes of an air bubble in a liquid water flow through a heated micro-channel

Purpose – The purpose of this paper is to study a gas bubble flowing in a micro-channel filled with liquid and to quantify the compressibility effects induced in the bubble by a heat supply at the walls of the channel. Design/methodology/approach – The paper presents a model and its numerical implementation. A hybrid method combining front-tracking techniques and a Heaviside step function is introduced to ensure an accurate satisfaction of the mass and energy conservation laws. Findings – Compressibility effects in the bubble are quantified. Test cases for numerical simulations of two-phase flows involving heat transfer are proposed. Originality/value – The authors present original test cases in which expansion or compression of a gas bubble flowing in a liquid are induced by heat transfer at the wall.

Numerical analysis of evaporation in microchannel under capillary pumping

Abstract The paper presents numerical simulation of two-phase flow in a heated capillary with evaporation on the meniscus. To solve the problem, a model of evaporation from meniscus was developed in which the dynamics of liquid-vapour interface is described by the Cahn-Hilliard equation. The numerical simulations were performed using commercial software for 2D axially symmetric case. The flow evolution was analysed for different values of heat transfer coefficient at the capillary wall and inlet liquid mass flow rate.

The response of two-phase hydrothermal systems to changing magmatic heat input at mid-ocean ridges

We develop numerical simulations of two-phase flow in a NaCl–H2O fluid subject to time varying basal heat flux boundary conditions in order to understand the response of hydrothermal vent temperature and salinity to changing magmatic heat input. The results show that periodic changes in basal heat input on a time scale of several years will not be detected in a continuous time series record of temperature measurements. Fluctuations in vent salinity may be recorded, however. For models with monotonic decay of the magmatic heat flux, a decline in vent temperature may not be observed for several years; however, once single phase conditions are established at the base of the system, a pulse of brine-derived fluid is expected to appear at the surface, followed by a gradual decline of salinity to the seawater value. The pulse of brine salinity is expected to occur before an observed decline in vent temperature. Observed rapid changes in vent temperature and salinity associated with either eruptive or non-eruptive magmatic events are not likely a result of changes in basal heat flux.

Numerical Simulations of Two-Phase Flow in a Self-Aerated Flotation Machine and Kinetics Modeling

A new boundary condition treatment has been devised for two-phase flow numerical simulations in a self-aerated minerals flotation machine and applied to a Wemco 0.8 m3 pilot cell. Airflow rate is not specified a priori but is predicted by the simulations as well as power consumption. Time-dependent simulations of two-phase flow in flotation machines are essential to understanding flow behavior and physics in self-aerated machines such as the Wemco machines. In this paper, simulations have been conducted for three different uniform bubble sizes (db = 0.5, 0.7 and 1.0 mm) to study the effects of bubble size on air holdup and hydrodynamics in Wemco pilot cells. Moreover, a computational fluid dynamics (CFD)-based flotation model has been developed to predict the pulp recovery rate of minerals from a flotation cell for different bubble sizes, different particle sizes and particle size distribution. The model uses a first-order rate equation, where models for probabilities of collision, adhesion and stabilization and collisions frequency estimated by Zaitchik-2010 model are used for the calculation of rate constant. Spatial distributions of dissipation rate and air volume fraction (also called void fraction) determined by the two-phase simulations are the input for the flotation kinetics model. The average pulp recovery rate has been calculated locally for different uniform bubble and particle diameters. The CFD-based flotation kinetics model is also used to predict pulp recovery rate in the presence of particle size distribution. Particle number density pdf and the data generated for single particle size are used to compute the recovery rate for a specific mean particle diameter. Our computational model gives a figure of merit for the recovery rate of a flotation machine, and as such can be used to assess incremental design improvements as well as design of new machines.