垂直上升气液环状流流动参数的数值计算

提出了一个新的气核-液膜耦合模型来求解垂直上升气液环状流在充分发展段的流动参数.本模型考虑了液膜、气核以及它们之间的相互影响和作用.模型中基本的气核区域和液膜区域的质量和动量方程由Fluent6.3.26进行求解,而液滴方程以及相界面上的夹带和沉积作用通过用户自定义接口函数UDF(User Defined Functi...     

A model for pressure loss, film thickness, and entrained fraction for gas-liquid annular flow

A two-component (air-water) annular flow model is presented requiring only flow rates, absolute pressure, temperature, and tube diameter. Film thicknesses (base film and wave height) are calculated from a critical film thickness model. Modeled pressure gradient is weighted by wave intermittency to compute average pressure gradient. Film flow rate and wave velocity are estimated using the universal velocity profile in the waves and a piecewise linear profile in the base film. For vertical flow, mean absolute errors for film thickness, wave velocity, and pressure gradient are 9%, 9%, and 19%, respectively. In horizontal flow, mean absolute errors for pressure gradient, base film thickness, and disturbance wave velocity are 17%, 10%, and 14%, respectively, on par with those from single-behavior models that require additional film thickness or other data as inputs.     

Experimental study on axial development of liquid film in vertical upward annular two-phase flow

Accurate measurements of the interfacial wave structure of upward annular two-phase flow in a vertical pipe were performed using a laser focus displacement meter (LFD). The purpose of this study was to clarify the effectiveness of the LFD for obtaining detailed information on the interfacial displacement of a liquid film in annular two-phase flow and to investigate the effect of axial distance from the air–water inlet on the phenomena. Adiabatic upward annular air–water flow experiments were conducted using a 3 m long, 11 mm ID pipe. Measurements of interfacial waves were conducted at 21 axial locations, spaced 110 mm apart in the pipe. The axial distances from the inlet (z) normalized by the pipe diameter (D) varied over z/D = 50–250. Data were collected for predetermined gas and liquid flow conditions and for Reynolds numbers ranging from Re_G = 31,800 to 98,300 for the gas phase and Re_L = 1050 to 9430 for the liquid phase. Using the LFD, we obtained such local properties as the minimum thickness, maximum thickness, and passing frequency of the waves. The maximum film thickness and passing frequency of disturbance waves decreased gradually, with some oscillations, as flow developed. The flow development, i.e., decreasing film thickness and passing frequency, persisted until the end of the pipe, which means that the flow might never reach the fully developed state. The minimum film thickness decreased with flow development and with increasing gas flow rate. These results are discussed, taking into account the buffer layer calculated from Karman’s three-layer model. A correlation is proposed between the minimum film thickness obtained in relation to the interfacial shear stress and the Reynolds number of the liquid.     

Prediction of amount of entrained droplets in vertical annular two-phase flow

Prediction of amount of entrained droplets or entrainment fraction in annular two-phase flow is essential for the estimation of dryout condition and analysis of post dryout heat transfer in light water nuclear reactors and steam boilers. In this study, air–water and organic fluid (Freon-113) annular flow entrainment experiments have been carried out in 9.4 and 10.2 mm diameter test sections, respectively. Both the experiments covered three distinct pressure conditions and wide range of liquid and gas flow conditions. The organic fluid experiments simulated high pressure steam–water annular flow conditions. In each experiment, measurements of entrainment fraction, droplet entrainment rate and droplet deposition rate have been performed by using the liquid film extraction method. A simple, explicit and non-dimensional correlation developed by Sawant [Sawant, P.H., Ishii, M., Mori, M., 2008. Droplet entrainment correlation in vertical upward co-current annular two-phase flow. Nucl. Eng. Des. 238 (6), 1342–1352] for the prediction of entrainment fraction is further improved in this study in order to account for the existence of critical gas and liquid flow rates below which no entrainment is possible.Additionally, a new correlation is proposed for the estimation of minimum liquid film flow rate at the maximum entrainment fraction condition. The improved correlation successfully predicted the newly collected air–water and Freon-113 entrainment fraction data. Furthermore, the correlations satisfactorily compared with the air–water, helium–water and air–genklene experimental data measured by Willetts [Willetts, I.P., 1987. Non-aqueous annular two-phase flow. D.Phil. Thesis, University of Oxford]. However, comparison of the correlations with the steam–water data available in literature showed significant discrepancies. It is proposed that these discrepancies might have been caused due to the inadequacy of the liquid film extraction method used to measure the entrainment fraction or due to the change in mechanism of entrainment under high liquid flow conditions.     

Turbulence modification in vertical upward annular flow passing through a throat section

Experimental studies on the turbulence modification in annular two-phase flow passing through a throat section were carried out. The turbulence modification in multi-phase flow due to the interactions between two-phases is one of the most interesting scientific issues and has attracted research attention. In this study, the gas-phase turbulence modification in annular flow due to the gas–liquid phase interaction is experimentally investigated. The annular flow passing through a throat section is under the transient state due to the changing cross sectional area of the channel and resultantly the superficial velocities of both phases are changed compared with a fully developed flow in a straight pipe. The measurements for the gas-phase turbulence were precisely performed by using a constant temperature hot-wire anemometer, and made clear the turbulence structure such as velocity profiles, fluctuation velocity profiles. The behavior of the interfacial waves in the liquid film flow such as the ripple or disturbance waves was also observed. The measurements for the liquid film thickness by the electrode needle method were also performed to measure the base film thickness, mean film thickness, maximum film thickness and wave height of the ripple or the disturbance waves.     

Three-phase upward flow in a vertical pipe

In this paper we develop an approach to design a three-phase, gas–solid–liquid flow system that transports pneumatically scarified solid particles, including sticky ones, through a vertical pipe. The proposed system permits the introduction and maintenance of a liquid film that coats the pipe’s inner wall and acts as a lubricant that ensures sticky particles continue to move upward without permanently adhering to the pipe wall. The system’s operating conditions fall within the boundaries of the annular dispersed region on a typical flow pattern map of vertical flow of a gas–liquid mixture. High gas superficial velocities combined with low liquid superficial velocities characterize such a region. A combination of a modified one-dimensional, two-fluid annular dispersed flow model and a one-dimensional pneumatic conveying model is shown to describe this transport process satisfactorily. Solution of the combined models produces all the necessary design parameters including power requirements and superficial velocities of the two-fluid media needed to transport a given amount of solid particles. Results of model calculations are compared with rare three-phase flow data obtained prior to the development of the present model, by an independent experimental team that used the physical conditions of the present approach. Reasonable agreement justifies the use of the combined model for engineering design purposes.     

A model for droplet entrainment in heated annular flow

The ability to accurately predict droplet entrainment in annular two-phase flow is required to effectively calculate the interfacial mass, momentum, and energy transfer, which characterizes nuclear reactor safety, system design, analysis, and performance. Most annular flow entrainment models in the open literature are formulated in terms of dimensionless groups, which do not directly account for interfacial instabilities. However, many researchers agree that there is a clear presence of interfacial instability phenomena having a direct impact on droplet entrainment. The present study proposes a model for droplet entrainment, based on the underlying physics of droplet entrainment from upward co-current annular film flow that is characteristic to light water reactor safety analysis. The model is developed based on a force balance and stability analysis that can be implemented into a transient three-field (continuous liquid, droplet, and vapor) two-phase heat transfer and fluid flow systems analysis computer code.     

A separated flow model for predicting two-phase pressure drop and evaporative heat transfer for vertical annular flow

A separated flow model has been developed that is applicable to vertical annular two-phase flow in the purely convective heat transfer regime. Conservation of mass, momentum, and energy are used to solve for the liquid film thickness, pressure drop, and heat transfer coefficient. Closure relationships are specified for the interfacial friction factor, liquid film eddy-viscosity, turbulent Prandtl number, and entrainment rate. Although separated flow models have been reported previously, their use has been limited, because they were tested over a limited range of flow and thermal conditions. The unique feature of this model is that it has been tested and calibrated against a vast array of two-phase pressure drop and heat transfer data, which include upflow, downflow, and microgravity flow conditions. The agreements between the measured and predicted pressure drops and heat transfer coefficients are, on average, better or comparable to the most reliable empirical correlations. This separated flow model is demonstrated to be a reliable and practical predictive tool for computing two-phase pressure drop and heat transfer rates. All of the datasets have been obtained from the open literature.     

An improved semi-empirical friction model for gas-liquid two-phase flow in horizontal and near horizontal pipes

Pressure drop and liquid hold-up are two very important fluid flow parameters in design and control of multiphase flow pipelines. Friction factors play an important role in the accurate calculation of pressure drop. Various empirical and semi-empirical closure relations exist in the literature to calculate the liquid-wall, gas-wall and interfacial friction in two-phase pipe flow.However most of them are empirical correlations found under special experimental conditions. In this paper by modification of a friction model available in the literature, an improved semiempirical model is proposed. The proposed model is incorporated in the two-fluid correlations under equilibrium conditions and solved. Pressure gradient and velocity profiles are validated against experimental data. Using the improved model, the pressure gradient deviation from experiments diminishes by about 3%; the no-slip condition at the interface is satisfied and the velocity profile is predicted in better agreement with the experimental data.     

A multiphase,micro-scale PIV measurement technique for liquid film velocity measurements in annular two-phase flow

Prediction methods for two-phase annular flow require accurate knowledge of the velocity profile within the liquid film flowing at its perimeter as the gradients within this film influence to a large extent the overall transport processes within the entire channel. This film, however, is quite thin and variable and traditional velocimetry methods have met with only very limited success in providing velocity data. The present work describes the application of Particle Image Velocimetry (PIV) to the measurement of velocity fields in the annular liquid flow. Because the liquid is constrained to distances on the order of a millimeter or less, the technique employed here borrows strategies from micro-PIV, but micro-PIV studies do not typically encounter the challenges presented by annular flow, including very large velocity gradients, a free surface that varies in position from moment to moment, the presence of droplet impacts and the passage of waves that can be 10 times the average thickness of the base film. This technique combines the seeding and imaging typical to micro-PIV with a unique lighting and image processing approach to deal with the challenges of a continuously varying liquid film thickness and interface. Mean velocity data are presented for air–water in two-phase co-current upward flow in a rectangular duct, which are the first detailed velocity profiles obtained within the liquid film of upward vertical annular flow to the authors’ knowledge. The velocity data presented here do not distinguish between data from waves and data from the base film. The resulting velocity profiles are compared with the classical Law of the Wall turbulent boundary layer model and found to require a decreased turbulent diffusivity for the model to predict well. These results agree with hypotheses previously presented in the literature.     

Studying disturbance waves in vertical annular flow with high-speed video

In many annular two-phase gas–liquid flows, large disturbance waves propagate liquid mass. These waves are important for modeling of gas-to-liquid momentum transfer and liquid film behavior. High-speed videos of vertical upflow have been analyzed to extract individual and average wave data. Two types of structures, coherent waves and piece waves, have been identified in these flows. Velocities, lengths, and temporal spacings of individual waves and average velocities, lengths, frequencies, and intermittencies have been studied as functions of both gas and liquid flow rates. Velocity and frequency increase with liquid and gas flow rates, length decreases with increasing gas flow and increases with increasing liquid flow, and intermittency is predominantly an increasing function of liquid flow.     

Planar laser-induced fluorescence (PLIF) measurements of liquid film thickness in annular flow. Part I: Methods and data

Most approaches to the modeling of annular flow require information regarding the thin liquid film surrounding the central gas core. This film is hypothesized to present a rough surface to the gas core, enhancing interfacial shear and pressure loss, with the roughness closely linked to the height of the film. This height is typically obtained from conductance probe measurements. The present work used planar laser-induced fluorescence to provide direct visualization of the liquid film in upward vertical air–water annular flow. Images were processed to produce the distribution of film heights. The standard deviation and average film thickness are found to be an increasing function of liquid flow and a decreasing function of gas flow, with the standard deviation approaching 0.4 times the average at sufficient liquid flow.     

Planar laser-induced fluorescence (PLIF) measurements of liquid film thickness in annular flow. Part II: Analysis and comparison to models

Film thickness distributions in upward vertical air–water annular flow have been determined using planar laser-induced fluorescence (PLIF). Film thickness data are frequently used to estimate interfacial shear and pressure loss. This film roughness concept has been used in a number of models for annular flow of varying complexity. The PLIF data are presently applied to the single-zone interfacial shear correlation of Wallis; the more detailed model of Owen and Hewitt; and the two-zone (base film and waves) model of Hurlburt, Fore, and Bauer. For the present data, these models all under-predict the importance of increasing liquid flow on pressure loss and interfacial shear. Since high liquid flow rates in annular flow induce disturbance wave and entrainment activity, further modeling in these areas is advised.     

Flow regime development analysis in adiabatic upward two-phase flow in a vertical annulus

In this work radial and axial flow regime development in adiabatic upward air-water two-phase flow in a vertical annulus has been investigated. Local flow regimes have been identified using conductivity probes and neural networks techniques. The inner and outer diameters of the annulus are 19.1 mm and 38.1 mm, respectively. The equivalent hydraulic diameter of the flow channel, D_H, is 19.0 mm and the total length is 4.37 m. The flow regime map includes 1080 local flow regimes identifications in 72 flow conditions within a range of 0.01 m/s < 〈j_g〉 < 30 m/s and 0.2 m/s < 〈j_f〉 < 3.5 m/s where 〈j_g〉 and 〈j_f〉 are, respectively, superficial gas and liquid velocities. The local flow regime has been classified into four categories: bubbly, cap-slug, churn-turbulent and annular flows. In order to study the radial and axial development of flow regime the measurements have been performed at five radial locations. The three axial positions correspond to z/D_H = 52, 149 and 230, where z represents the axial position. The flow regime indicator has been chosen as some statistical parameters of local bubble chord length distributions and self-organized neural networks have been used as mapping system. This information has been also used to compare the results given by the existing flow regime transition models. The local flow regime is characterized basically by the void fraction and bubble chord length. The radial development of flow regime shows partial and complete local flow regime combinations. The radial development is controlled by axial location and superficial liquid velocity. The radial flow regime transition is always initiated in the center of the flow channel and it is propagated towards the channel boundaries. The axial development of flow regime is observed in all the flow maps and it is governed by superficial liquid velocity and radial location. The prediction results of the models are compared for each flow regime transition.     

Algebraic turbulence modeling in adiabatic and evaporating annular two-phase flow

The study considers algebraic turbulence modeling in adiabatic and evaporating annular two-phase flow, focusing in particular on momentum and heat transfer (so-called ‘convective boiling’) through the annular liquid film. In contrast with single-phase wall-bounded flow theory, shear-driven annular liquid films are assumed here to behave as fluid-bounded flows, mostly interacting with the shearing gas-entrained droplets core flow. Besides providing velocity and temperature profiles through the liquid film, the turbulence model proposed here predicts key parameters such as the average liquid film thickness, the void fraction and the convective boiling heat transfer coefficient with accuracies comparable or better than those of leading design correlations. This turbulence model is part of a unified annular flow modeling suite that includes methods to predict the entrained liquid fraction and the axial frictional pressure gradient. The underlying heat transfer database covers nine fluids (water, two hydrocarbons and six refrigerants) for vertical and horizontal tubes of 1.03-14.4 mm diameter and pressures of 0.1-7.2 MPa. Importantly, this study shows that there appears to be no macro-to-microscale transition when it comes to annular flow. Simply better physical modeling is required to span this range.     

Time and spatially resolved measurements of interfacial waves in vertical annular flow

Film thickness measurements have been performed in a vertical air/water annular flow in a pipe of 0.05 m diameter. A sensor has been built which allows to measure the film thickness evolution in time at 320 positions, such that the interface of the vertical annular flow can be reconstructed. The large-scale structures moving on the interface are described statistically, with a special attention to the disturbance waves. Probability density functions and mean statistics are given for the height, length, velocity, frequency and spatial distribution of the disturbance waves. In particular, it is shown that the disturbance waves are three-dimensional structures with large height fluctuations in the circumferential and axial direction, giving a meandering path between the maximum height around the circumference. It is also shown that the disturbance waves can flow with a slight inclination with respect to the axial direction. Finally, the disturbance waves are shown to be located randomly in space, within a Gamma distribution whose order only depends on the liquid superficial velocity. Due to the nature of the Gamma distribution, it could indicate that the spatial distribution of the disturbance waves results from a cascade of coalescence processes between the original disturbance waves on the film.     

A film flow model for analysing gravity-driven,thin wavy fluid films

To analyse the physics underlying gravity-driven runoff of thin wavy films, a film flow model is developed, and is solved with computational fluid dynamics. This model is based on the lubrication theory, and takes into account the gravitational, wall shear and surface tension forces. A key characteristic of the model is that it assumes only one computational cell over the film height, which enables studying film flow on larger computational domains. A main aim of this study is to perform a detailed validation of the numerical model. The film flow model is validated against several experiments of gravity-driven, thin fluid films on smooth surfaces. The time-averaged film thickness and the fluid speed profiles predicted by the model show very good agreement with experimental results. Similarly, the film flow model is able to predict the wave speeds with sufficient accuracy. The energy spectra of the waves, where higher frequency waves are present in film flows at higher Reynolds numbers, show an exponentially decaying trend at these high frequencies. The model performs better than the Nusselt equation for film flows, which under-predicts the time-averaged film thickness and over-predicts the time-averaged fluid speeds, even for flows at low Reynolds numbers. The film flow model is compared qualitatively for fingering behaviour. This model also allows to investigate film flows on large surfaces, which can be rough, curved and of complex geometrical shape.     

Numerical simulation of flow boiling for organic fluid with high saturation temperature in vertical porous coated tube

A semi-analytical model is developed for the prediction of flow boiling heat transfer inside vertical porous coated tubes. The model assumes that the forced convection and nucleate boiling coexist together in the annular flow regime. Conservations of mass, momentum, and energy are used to solve for the liquid film thickness and temperature. The heat flux due to nucleate boiling consists of those inside and outside micro-tunnels. To close the equations, a detailed analysis of various forces acting on the bubble is presented to predict its mean departure diameter. The active nucleation site density of porous layer is determined from the pool boiling correlation by introducing suppression factor. The flow boiling heat transfer coefficients of organic fluid (cumene) with high saturation temperature in a vertical flame-spraying porous coated tube are studied numerically. It is shown that the present model can predict most of the experimental values within ±20%. The numerical results also indicate that the nucleate boiling contribution to the overall heat transfer coefficient decreases from 50% to 15% with vapor quality increasing from 0.1 to 0.5.     

Counter current gas-liquid vertical flow, model for flow pattern and pressure drop

Counter-current flow pattern transition and pressure drop are modeled. Emphasis is placed on the understanding of the transition mechanisms from a mechanistic point of view.

Unlike the case of co-current flow, in counter-current flow, the situations of “no solution” as well as “multiple solutions” for the flow pattern and pressure drop exist. These possibilities are discussed and criteria for the actual flow pattern that will take place are suggested.

Some of the results are supported by data (from the literature), others are somewhat tentative suggesting future experimental verification is needed.