Lattice Boltzmann Simulation of Surface Tension Dominated Vortices Behaviour in Two-Phase Mixing Layer

Author(s):  
Y. Y. Yan ◽  
Y. Q. Zu

Surface tension dominating mixings and interfacial interactions are major phenomena of multiphase flow in microchannels and a variety of micro mixers. Such phenomena are concerned with interfacial interactions not only at fluid-solid interface but also at different fluids/phases interfaces. In this paper, vortices behaviours in a mixing layer of two immiscible fluids are studied numerically. The lattice Boltzmann method (LBM) is employed to simulat surface tension dominated mixing process. As a mesoscopic numerical method, the LBM has many advantages, which include the ability of incorporating microscopic interactions, the simplicity of programming and the nature of parallel algorithm and is therefore ideal for simulating multiphase flow. In this article, the index function methodology of the LBM is employed to simulate surface tension dominated vertices behaviour in a two-dimensional immiscible two-phase mixing layer. The initial interface between two-fluids is evenly distributed around the midpoint in vertical direction. Different velocity perturbations which consist of a basic wave and a series of sub-harmonic waves are forced at the entrance of a rectangular mixing layer of the flow field. By changing the strength of surface tension and the combinations of perturbation waves, the effects of the surface tension and the velocity perturbation on vortices merging are investigated. The vortices contours and frequency spectrums are used to analyse the mechanism of vortices merging. Some interesting phenomena, which do not take place in a single-phase mixing layer, are observed and the corresponding mechanism is discussed in details.

2017 ◽  
Vol 28 (09) ◽  
pp. 1750120 ◽  
Author(s):  
Yong Peng ◽  
Yun Fei Mao ◽  
Bo Wang ◽  
Bo Xie

Equations of State (EOS) is crucial in simulating multiphase flows by the pseudo-potential lattice Boltzmann method (LBM). In the present study, the Peng and Robinson (P–R) and Carnahan and Starling (C–S) EOS in the pseudo-potential LBM with Exact Difference Method (EDM) scheme for two-phase flows have been compared. Both of P–R and C–S EOS have been used to study the two-phase separation, surface tension, the maximum two-phase density ratio and spurious currents. The study shows that both of P–R and C–S EOS agree with the analytical solutions although P–R EOS may perform better. The prediction of liquid phase by P–R EOS is more accurate than that of air phase and the contrary is true for C–S EOS. Predictions by both of EOS conform with the Laplace’s law. Besides, adjustment of surface tension is achieved by adjusting [Formula: see text]. The P–R EOS can achieve larger maximum density ratio than C–S EOS under the same [Formula: see text]. Besides, no matter the C–S EOS or the P–R EOS, if [Formula: see text] tends to 0.5, the computation is prone to numerical instability. The maximum spurious current for P–R is larger than that of C–S. The multiple-relaxation-time LBM still can improve obviously the numerical stability and can achieve larger maximum density ratio.


2013 ◽  
Vol 300-301 ◽  
pp. 1062-1066
Author(s):  
Yang Yu ◽  
Li Chen ◽  
Jian Hua Lu ◽  
Guo Xiang Hou

Free-surface model with surface tension and wall adhesion(wetting) is a very efficient technique to simulate two-phase flows with high density and viscosity ratios, such as etching and casting processes. In this paper, a conservative surface tension and wall adhesion model based on lattice Boltzmann single-phase free-surface method is proposed. The effectiveness of the model is demonstrated by simulating the flows induced by wall adhesion and surface tension, and filling processes in a 2D cavity.


2011 ◽  
Vol 9 (5) ◽  
pp. 1347-1361 ◽  
Author(s):  
Yoshito Tanaka ◽  
Masato Yoshino ◽  
Tetsuo Hirata

AbstractA thermal lattice Boltzmann method (LBM) for two-phase fluid flows in nucleate pool boiling process is proposed. In the present method, a new function for heat transfer is introduced to the isothermal LBM for two-phase immiscible fluids with large density differences. The calculated temperature is substituted into the pressure tensor, which is used for the calculation of an order parameter representing two phases so that bubbles can be formed by nucleate boiling. By using this method, two-dimensional simulations of nucleate pool boiling by a heat source on a solid wall are carried out with the boundary condition for a constant heat flux. The flow characteristics and temperature distribution in the nucleate pool boiling process are obtained. It is seen that a bubble nucleation is formed at first and then the bubble grows and leaves the wall, finally going up with deformation by the buoyant effect. In addition, the effects of the gravity and the surface wettability on the bubble diameter at departure are numerically investigated. The calculated results are in qualitative agreement with other theoretical predictions with available experimental data.


2005 ◽  
Vol 14 (1) ◽  
pp. 41-47 ◽  
Author(s):  
Wenchun Li ◽  
Guilin Hu ◽  
Zhe Zhou ◽  
Jianren Fan ◽  
Kefa Cen

1998 ◽  
Vol 09 (08) ◽  
pp. 1393-1404 ◽  
Author(s):  
D. J. Holdych ◽  
D. Rovas ◽  
J. G. Georgiadis ◽  
R. O. Buckius

Lattice-Boltzmann (LB) models provide a systematic formulation of effective-field computational approaches to the calculation of multiphase flow by replacing the mathematical surface of separation between the vapor and liquid with a thin transition region, across which all magnitudes change continuously. Many existing multiphase models of this sort do not satisfy the rigorous hydrodynamic constitutive laws. Here, we extend the two-dimensional, seven-speed Swift et al. LB model1 to rectangular grids (nine speeds) by using symbolic manipulation (MathematicaTM) and compare the LB model predictions with benchmark problems, in order to evaluate its merits. Particular emphasis is placed on the stress tensor formulation. Comparison with the two-phase analogue of the Couette flow and with a flow involving shear and advection of a droplet surrounded by its vapor reveals that additional terms have to be introduced in the definition of the stress tensor in order to satisfy the Navier–Stokes equation in regions of high density gradients. The use of Mathematica obviates many of the difficulties with the calculations "by-hand," allowing at the same time more flexibility to the computational analyst to experiment with geometrical and physical parameters of the formulation.


2012 ◽  
Vol 21 (12) ◽  
pp. 124703 ◽  
Author(s):  
Hai-Qiong Xie ◽  
Zhong Zeng ◽  
Liang-Qi Zhang ◽  
Gong-You Liang ◽  
Hiroshi Mizuseki ◽  
...  

2013 ◽  
Vol 10 (02) ◽  
pp. 335-370 ◽  
Author(s):  
PANTERS RODRÍGUEZ-BERMÚDEZ ◽  
DAN MARCHESIN

We study Riemann solutions for a system of two nonlinear conservation laws that models buoyancy-driven flow of three immiscible fluids in a porous medium, which do not exchange mass. We also assume that the fluids are incompressible and the flow occurs in the vertical spatial dimension. We consider the simplified case in which two of the three fluids have equal densities, obtaining the Riemann solutions by the wave curve method. As expected, the solutions contain waves traveling both upwards and downwards. The sequences of waves contain rarefactions, shocks (sometimes traveling with characteristic speed), and constant states. The shocks found in this work are proper or generalized Lax shock waves. The solutions we found are L1-continuous with respect to the initial data. Waves involving only two fluids often take part in three-phase flow Riemann solutions; this is the basis of a useful tool (the wedge construction) to obtain shocks separating states in distinct two-phase regimes having a common fluid. This tool is similar to fractional flow theory, or Oleinik's convex construction. In this investigation, the wave curve method from the theory of conservation laws is combined with numerical calculations.


2019 ◽  
Vol 863 ◽  
Author(s):  
Alessio Roccon ◽  
Francesco Zonta ◽  
Alfredo Soldati

We propose a physically sound explanation for the drag reduction mechanism in a lubricated channel, a flow configuration in which an interface separates a thin layer of less-viscous fluid (viscosity $\unicode[STIX]{x1D702}_{1}$) from a main layer of a more-viscous fluid (viscosity $\unicode[STIX]{x1D702}_{2}$). To single out the effect of surface tension, we focus initially on two fluids having the same density and the same viscosity ($\unicode[STIX]{x1D706}=\unicode[STIX]{x1D702}_{1}/\unicode[STIX]{x1D702}_{2}=1$), and we lower the viscosity of the lubricating layer down to $\unicode[STIX]{x1D706}=\unicode[STIX]{x1D702}_{1}/\unicode[STIX]{x1D702}_{2}=0.25$, which corresponds to a physically realizable experimental set-up consisting of light oil and water. A database comprising original direct numerical simulations of two-phase flow channel turbulence is used to study the physical mechanisms driving drag reduction, which we report between 20 and 30 percent. The maximum drag reduction occurs when the two fluids have the same viscosity ($\unicode[STIX]{x1D706}=1$), and corresponds to the relaminarization of the lubricating layer. Decreasing the viscosity of the lubricating layer ($\unicode[STIX]{x1D706}<1$) induces a marginally decreased drag reduction, but also helps sustaining strong turbulence in the lubricating layer. This led us to infer two different mechanisms for the two drag-reduced systems, each of which is ultimately controlled by the outcome of the competition between viscous, inertial and surface tension forces.


1985 ◽  
Vol 150 ◽  
pp. 381-394 ◽  
Author(s):  
Yuriko Renardy ◽  
Daniel D. Joseph

We consider the flow of two immiscible fluids lying between concentric cylinders when the outer cylinder is fixed and the inner one rotates. The interface is assumed to be concentric with the cylinders, and gravitational effects are neglected. We present a numerical study of the effect of different viscosities, different densities and surface tension on the linear stability of the Couette flow. Our results indicate that, with surface tension, a thin layer of the less-viscous fluid next to either cylinder is linearly stable and that it is possible to have stability with the less dense fluid lying outside. The stable configuration with the less-viscous fluid next to the inner cylinder is more stable than the one with the less-viscous fluid next to the outer cylinder. The onset of Taylor instability for one-fluid flow may be delayed by the addition of a thin layer of less-viscous fluid on the inner wall and promoted by a thin layer of more-viscous fluid on the inner wall.


Sign in / Sign up

Export Citation Format

Share Document