Numerical modeling of slug flow initiation in a horizontal channels using a two-fluid model

This paper presents a methodology for modeling slug initiation and growth in horizontal ducts. Transient two-fluid equations are solved numerically using a class of high-resolution shock capturing methods. The advantage of this method is that slug formation and growth in a stratified regime can be calculated directly from the solutions to the flow field differential equations. In addition, by using high-resolution shock capturing methods that do not contain numerical diffusion, the discontinuity generated by slugging in the flow field can be modeled with good accuracy. The two-fluid model is shown to be well-posed mathematically only under certain conditions. Under these circumstances, the two-fluid model is capable of correctly predicting and modeling the flow physics. When ill-posed, an unbounded instability occurs in the flow field solution, and the instability amplitude increases exponentially with decreasing mesh sizes. This work shows that there are three zones associated with slug formation. In addition, long wavelength slugs are shown to initiate from short wavelength waves. These short waves are generated at the interface of the two phases by the Kelvin-Helmholtz hydrodynamic instability. The results obtained through numerical modeling show good agreement with experimental results.     

Numerical simulation of compressible two-phase flow using a diffuse interface method

In this article, a high-resolution diffuse interface method is investigated for simulation of compressible two-phase gas–gas and gas–liquid flows, both in the presence of shock wave and in flows with strong rarefaction waves similar to cavitations. A Godunov method and HLLC Riemann solver is used for discretization of the Kapila five-equation model and a modified Schmidt equation of state (EOS) is used to simulate the cavitation regions. This method is applied successfully to some one- and two-dimensional compressible two-phase flows with interface conditions that contain shock wave and cavitations. The numerical results obtained in this attempt exhibit very good agreement with experimental results, as well as previous numerical results presented by other researchers based on other numerical methods. In particular, the algorithm can capture the complex flow features of transient shocks, such as the material discontinuities and interfacial instabilities, without any oscillation and additional diffusion. Numerical examples show that the results of the method presented here compare well with other sophisticated modeling methods like adaptive mesh refinement (AMR) and local mesh refinement (LMR) for one- and two-dimensional problems.     

Numerical simulation of internal gravity waves using a lattice gas model

Internal waves are modelled in two different circumstances: in a continuously stratified fluid and at the interface between two immiscible fluids. This is done using the lattice gas approach. The standard single phase model and an immiscible two-phase model are both modified to incorporate gravitational interactions. Standing internal waves are set up in both models and are seen to oscillate under the action of the gravitational interaction. The results obtained suggest that the lattice gas approach can be a useful tool in the modelling of such phenomena. © 1998 John Wiley & Sons, Ltd.     

Scaling laws for gas-solid riser flow through two-fluid model simulation

Scale up of gas-solid circulating fluidized bed (CFB) risers poses many challenges to researchers.In this paper,CFD investigation of hydrodynamic scaling laws for gas-solid riser flow was attempted on the basis of two-fluid model simulations,in particular,the recently developed empirical scaling law of Qi,Zhu,and Huang (2008).A 3D computational model with periodic boundaries was used to perform numerical experiments and to study the effect of various system and operating parameters in hydrodynamic scaling o...     

One-dimensional two-fluid equations for horizontal stratified two-phase flow

Advanced computer codes for water reactor loss-of-coolant analysis are based on the use of the two-fluid model of two-phase flow, in which conservation equations are solved for the gas and liquid phases separately. The standard two-fluid equations, however, sometimes predict the growth of instabilities in the flow, and occasionally become improperly posed. These difficulties have in the past led to the proposal of several different forms for the conservations equations.To help resolve these uncertainties a widely accepted form of the one-dimensional two-fluid equations is used to calculate wave propagation speeds, and stability limits, for the illustrative case of a frictionless horizontal stratified gas-liquid flow. Calculated propagation velocities are shown to agree with the appropriate limit of an exact solution, and the predicted stability limits are found consistent with available observations on the stability of the stratified flow regime.These comparisons help improve confidence in the ability of the two-fluid equations to analyse more complex problems in transient two-phase flow.     

Numerical simulation and homogenization of two-phase flow in heterogeneous porous media

A mathematically rigorous method of homogenization is presented and used to analyze the equivalent behavior of transient flow of two incompressible fluids through heterogeneous media. Asymptotic expansions and H-convergence lead to the definition of a global or effective model of an equivalent homogeneous reservoir. Numerical computations to obtain the homogenized coefficients of the entire reservoir have been carried out via a finite element method. Numerical experiments involving the simulation of incompressible two-phase flow have been performed for each heterogeneous medium and for the homogenized medium as well as for other averaging methods. The results of the simulations are compared in terms of the transient saturation contours, production curves, and pressure distributions. Results obtained from the simulations with the homogenization method presented show good agreement with the heterogeneous simulations.     

Numerical simulations of counter-current two-phase flow experiments in a PWR hot leg model using an interfacial area density model

In order to improve the understanding of counter-current two-phase flows and to validate new physical models, CFD simulations of 1/3rd scale model of the hot leg of a German Konvoi PWR with rectangular cross section was performed. Selected counter-current flow limitation (CCFL) experiments at the Helmholtz–Zentrum Dresden–Rossendorf (HZDR) were calculated with ANSYS CFX 12.1 using the multi-fluid Euler–Euler modeling approach. The transient calculations were carried out using a gas/liquid inhomogeneous multiphase flow model coupled with a k-ω turbulence model for each phase. In the simulation, the surface drag was approached by a new correlation inside the Algebraic Interfacial Area Density (AIAD) model. The AIAD model allows the detection of the morphological form of the two phase flow and the corresponding switching via a blending function of each correlation from one object pair to another. As a result this model can distinguish between bubbles, droplets and the free surface using the local liquid phase volume fraction value. A comparison with the high-speed video observations shows a good qualitative agreement. The results indicated that quantitative agreement of the CCFL characteristics between calculation and experimental data was obtained. The goal is to provide an easy usable AIAD framework for all Code users, with the possibility of the implementation of their own correlations.     

Numerical calculations for two-phase jet flow of a mixture of pulverized-coal and gas

In this paper, Euler-Lagrange type equations are used to describe the jet flow of a mixture of pulverized-coal and gas, which is an unsteady axisymmetric two-phase flow. By means of the finite-difference method, the coal particle's distribution, velocity and trajectory in the flow field are obtained. The coal particles are represented by a finite number of computational particles. Each particle's diameter is randomly assigned according to a given distribution. The states of the computational particles are different from each other. Turbulence is accounted for in a stochastic model. Explicit time-splitting scheme is used to calculate the strongly coupling interphase term. The numerical results are reasonable. The comparison between the numerical results and the experiment data for the case of the oil droplet injection shows good agreement. This numerical technique can be extended to the calculation of other two-phase flows of dilute particles or a droplet system. Mr. Mei Renwei also participated in the work of this paper.     

考虑活塞回复的侧向后喷武器两相流数值模拟

为研究活塞回复运动对火药燃气流动的影响,基于两相流理论对活塞控制侧向后喷武器的发射过程进行了数值模拟研究。考虑控制侧向后喷通道开闭的活塞-弹簧系统的往复运动,建立了结合膛内气固两相流、活塞腔内流固耦合和侧向排气管内气体瞬态流动的武器发射过程数学模型,并将数值模拟结果与相关文献进行了比较验证。得到了该武器发射过程中膛内流场分布与稀疏波传播特性,并与普通武器的膛内流场进行了对比分析。进一步研究了活塞回复运动对火药燃气流动和减后坐效率的影响。结果表明:相对于不考虑活塞的回复运动,在弹丸初速都降低1.52%的情况下,因为活塞回复关闭后喷通道,其减后坐效率由38.86%下降到32.88%,说明在此类武器研究中,不可忽视活塞回复运动。     

Numerical simulation of drag-reducing channel flow by using bead-spring chain model

Numerical simulations of the drag-reducing turbulent channel flow caused by polymer addition are performed. A bead-spring chain model is employed as a model of polymer aggregation. The model consists of beads and springs to represent the polymer dynamics. Three drag-reduction cases are studied with different spring constants that correspond to the relaxation time of the polymer. The energy budget is mainly focused upon to discuss the drag-reduction mechanism. Our results show that a decreasing pressure-strain correlation mainly contributes to strengthening the anisotropy of the turbulence. Furthermore, energy transport by the polymer models attenuates the turbulence. These viscoelastic effects on the drag-reducing flow are intensified with decreasing spring constant. By visualizing the flow field, it is found that this polymer energy transport is related to the orientation of the polymer.     

MPS -DEM 方法数值模拟超长管道固液混输流动

近些年来,深海矿产开发成为热点,其中涉及的固液混输问题逐渐引起学者的关注。如矿粒在海床上的输运,涉及水平管道的输运问题;矿粒从海床输运到海上浮式建筑物,涉及垂直管道的输运问题;矿粒在不同倾角管道内的输运,涉及倾斜管道的输运问题。对于该类问题,由于流体和固体之间的物理差异较大,对于固液两相往往采用不同的数值处理方法。本文使用MPS-DEM耦合方法模拟管道内的固液混输运动。用MPS方法模拟管道内流体的运动,用DEM方法追踪管道内颗粒的运动。在管道的进出口使用周期性边界条件,将多段管道之间进行拼接,从而实现了对超长管道内高雷诺数的颗粒输运现象的模拟。最后,分析管道颗粒和流场的速度及粒子分布等信息,分析流动特征。得到的结果会为超长管道内固液混输现象的研究提供一些借鉴意义及参考价值。     

Closure relations effects on the prediction of the stratified two-phase flow stability via the two-fluid model

The possibility of predicting the exact long wave linear stability boundary via the two-fluid (TF) model for horizontal and inclined stratified two-phase flow is examined. The application of the TF model requires the introduction of empirical closure relations for the velocity profile shape factors and for the wave induced wall and interfacial shear stresses. The latter are recognized as the problematic closure laws. In order to explore the closure relations effects and to suggest the necessary modifications that can improve the stability predictions of the TF model, the results are compared with the exact long wave solution of the Orr–Sommerfeld equations for the two-plate geometry. It is demonstrated that with the shape factors corrections and the inclusion of wave induced stresses effects, the TF model is able to fully reproduce the exact long wave neutral stability curves. The wave induced shear stresses in phase with the wave slope, which give rise to the so called “sheltering force”, were found to have a remarkable destabilizing effect in many cases of horizontal and inclined flows. In such cases, the sheltering effects must be included in the TF model, otherwise the region of smooth stratified flow would be significantly over predicted. Based on the results of the exact analysis, a simple closure relation for the sheltering term in the TF model is provided.     

Numerical simulation of a low-emission gas turbine combustor using KIVA-II

A numerical study was performed to investigate chemically reactive flows with sprays inside a staged turbine combustor (STC) using a modified version of the KIVA-II code. This STC consists of a fuel nozzle (FN), a rich-burn (RB) zone, a converging connecting pipe, a quick-quench (QQ) zone, a diverging connecting pipe and a lean-combustion (LC) zone. From the computational viewpoint, it is more efficient to split the STC into two subsystems, called FN/RB zone and QQ/LC zones, and the numerical solutions were obtained separately for each subsystem. This paper addresses the numerical results of the STC which is equipped with an advanced airblast fuel nozzle. The airblast nozzle has two fuel injection passages and four air flow passages. The input conditions used in this study were chosen similar to those encountered in advanced combustion systems. Preliminary results generated illustrate some of the major features of the flow and temperature fields inside the STC. Velocity, temperature and some critical species information inside the FN/RB zone are given. Formation of the co- and counter-rotating bulk flow and the sandwiched-ring-shaped temperature field, typical of the confined inclined jet-in-cross-flow, can be seen clearly in the QQ/LC zones. The calculations of the mass-weighted standard deviation and the pattern factor of temperature revealed that the mixing performance of the STC is very promising. The temperature of the fluid leaving the LC zone is very uniform. Prediction of the NO_x emission shows that there is no excessive thermal NO_x produced in the QQ/LC zones for the case studied. From the results obtained so far, it appears that the modified KIVA-II code can be used to guide the low-emission combustion experiments.     

长药室装药中多点点传火一维两相流数值模拟

针对药室长度超过3m的长药室装药设计问题,设计了以电点火具为点火源的多点点火系统,在点传火模拟实验装置上进行了多点点火的实验研究.同时,建立了点传火模拟实验装置中多点点传火过程的数学物理模型,采用MacCormack差分格式进行了数值求解,得到了点传火模拟装置中传火管的压力分布、固相和气相速度以及空隙率等数值解,分析了...     

Numerical simulation of annular two-phase flow considering the four involved mass transfers

In the present study, a numerical model is developed for simulation of annular two-phase flow considering bubbly flow regime in the liquid film along with the four involved mechanisms of mass transfer those are evaporation, entrainment, deposition and condensation. In the numerical approach, liquid film accompanied by fine nucleated bubbles are simulated with innovative model named suction model, the whole domain containing liquid film and the vapor core is simulated by volume of fluid model. While the vapor and the entrained droplets are treated as homogeneous flow. The interface between the liquid and the vapor is traced by level set formulation. The model is then validated by experimental models of Lee & Lee and Stevanovic et al. and shows a good precision such that it predicts the experimental results of Stevanivic et al. Better than their own numerical model. This issue is due to the least possible simplifying assumptions along with considering the effect of boiling in liquid film and all mechanisms of mass transfer in the fluid flow.     

A two-fluid model for critical vapour-liquid flow

One-dimensional two-fluid equations are used to calculate the mass flux of initially saturated or subcooled water discharging from a pipe in critical flow. The model allows in a general way for thermal non-equilibrium between the liquid and vapour bubbles, and for interphase relative motion. The theory is shown to be in good agreement with the measured critical flow-rates of pressurised water over a wide pressure range for a single choice of parameters characterising (i) the density of nucleation sites in the liquid, (ii) the liquid superheat required to cause bubble nucleation. Predictions are made of the critical mass flux of initially saturated water in pipes of the range of sizes of interest in water-reactor blowdown safety analysis. Results indicate that for pipes up to ten diameters in length flows will be significantly higher than values obtained from conventional homogeneous thermal equilibrium flow theory.     

An investigation of pressure transients in pipelines with two-phase bubbly flow

This paper presents a two-dimensional model for the analysis of the pressure transient of a two-phase homogeneous bubbly mixture flowing in a pipeline and the numerical integration using the centre implicit method (CIM). Experiments were conducted to confirm the proposed sonic speed equation of an air–water mixture for an air concentration of less than 1%. The 2D CIM model is compared with the method of characteristics (MoC) for a two-phase bubbly flow in a pipeline. The comparisons show that the proposed 2D CIM model generally gives good agreement with the method of characteristics.     

Numerical investigation of two-phase rotating flow

Finite difference solutions of the two-fluid equations of motion for a particle (droplet)-fluid mixture in a rotating finite axisymmetric cylinder are presented. The numerical method, which can be regarded as an extension of the Harlow & Amsden approach, employs forward time and centred space discretization and treats implicitly the pressure, Coriolis and volume flux terms. The computed flow fields are examined via a detailed comparison to previous analytic approximations, which illuminates both the physical and numerical aspects and the validity of these approximations.