The onset of slugging in horizontal stratified air-water flow

A simple criterion for the transition from stratified flow to the slug or plug flow regime in horizontal rectangular ducts is presented. A theoretical model is developed and shown to be consistent with data obtained by blowing air over stationary water in channels with a wide range of geometrical parameters. Some data for cocurrent and countercurrent flow are also correlated successfully.     

Modulation of an internal gravity wave packet in a stratified shear flow

Modulations of an internal gravity wave packet in a stratified shear flow are discussed in the weakly nonlinear and weakly dispersive context. It is shown that the modulations are described by a variable coefficient nonlinear Schrödinger equation when the modulations are confined to the direction of wave propagation. Transverse modulations couple the nonlinear Schrödinger equation to the mean flow equations. For long waves, it is shown that the modulation equations may be somewhat simplified. An Appendix describes the equations governing long wave resonance.     

Characteristics of air-water flow in an emptying tank under different conditions

This paper experimentally studied the features of air-water flow during the emptying of a water-filled prismatic tank with a bottom orifice under different conditions. The experiments were conducted with both circular and elliptical orifices, with and without ventilation. The evolution of bubbles, water pressure variation, and water level change with time were recorded in the experiments and analyzed. Based on the results, the evolution of bubbles could be mainly divided into three stages of formation, deformation, and decomposition. Ventilation was found important to the emptying process, with which the drainage efficiency was much higher than that under the unventilated condition. Additionally, under the unventilated condition, the drainage efficiency with the circular orifice was slightly higher than that with the elliptical orifice.     

Effect of pipe diameter on flow patterns for air-water flow in horizontal pipes

New results are presented on interfacial patterns observed for air and water flowing in horizontal 2.54 and 9.53 cm pipelines close to atmospheric conditions. This work differs from previous studies in that measurements of pressure fluctuations at two locations separated in the streamwise direction are used to detect slugs. The liquid flow needed to initiate slugs at low gas velocities is strongly affected by pipe diameter and appears to depend on a linear instability. At high gas velocities the transition is approximately independent of pipe diameter and is explained by a nonlinear mechanism associated with the coalescence of roll waves. The initiation of slugs in the annular flow regime is determined to occur at much lower liquid flows than had been reported by previous investigators. The transition from stratified to annular flow is different in smaller-diameter pipes than in larger pipes because wave wetting plays a more important role.     

Effect of flow-obstruction geometry on pressure drops in horizontal air-water flow

Experiments were carried out to investigate the pressure drop due to flow obstructions in horizontal air-water flow. The axial pressure distribution along a 25.4 mm inside diameter tube, with and without flow obstructions was measured using multi-tube manometers. Various obstruction shapes and sizes were investigated. Pressure-loss coefficients and two-phase multipliers were derived for twelve different flow obstructions.Estimates of the kinetic energy and momentum of the flow were also obtained from radial void distribution measurements. The pressure drop through the obstructions in two-phase flow was found to depend strongly on the kinetic energy and momentum of the liquid intercepted by the flow obstruction. Buoyancy-induced flow stratification caused a strongly non-symmetrical effect, and resulted in large pressure drops for flow obstructions located in the bottom part of the channel.     

Experimental investigation of horizontal forced-vibration effect on air-water two-phase flow

In order to investigate the potential seismic vibrations effect on two-phase flow in an annular channel, experimental tests with air-water two-phase flow under horizontal vibrations were carried out. A low-speed eccentric-cam vibration module capable of operating at motor speed of 45–1200 rpm (f = 0.75–20 Hz) was attached to an annular channel, which was scaled down from a prototypic BWR fuel sub-channel with inner and outer diameters of 19.1 mm and 38.1 mm, respectively. The two-phase flow was operated in the ranges of 〈j_f〉 = 0.25–1.00 m/s and 〈j_g〉 = 0.03–1.46 m/s with 27 flow conditions, and the vibration amplitudes controlled by cam eccentricity (E) were designed for the range of 0.8–22.2 mm. Ring-type impedance void meters were utilized to detect the area-averaged time-averaged void fraction under stationary and vibration conditions. A systematic experimental database was built and analyzed with effective maps in terms of flow conditions (〈j_g〉-〈j_f〉) and vibration conditions (E-f and f-a), and the potential effects were expressed by regions on the maps. In the 〈j_g〉-〈j_f〉 maps, the void fraction was found to potentially decrease under vibrations in bubbly flow regime and relatively lower liquid flow conditions, which may be explained by the increase of distribution parameter. Whereas and the void fraction may increase at the region closed to bubbly-to-slug transition boundary under vibrations, which may be explained by the changes of drift velocity due to flow regime change from bubbly to slug flows. No significant change in void fraction was found in slug flow regime under the present test conditions.     

A charnock-based estimate of interfacial resistance and roughness for internal,fully-developed,stratified, two-phase horizontal flow

Charnock's (1955) relation between interfacial roughness and drag is utilized in order to determine these two quantities for fully-developed flows in tubes of arbitrary cross-section.     

Heat transfer and interfacial drag in countercurrent steam-water stratified flow

Local condensation heat transfer coefficients and interfacial shear stresses have been measured for countercurrent stratified flow of steam and subcooled water in rectangular channels over a wide range of inclination angles (4–87°) at two aspect ratios. Dimensionless correlations for the interfacial friction factor have been developed that show that it is a function of the liquid Reynolds number only. Empirical correlations of the heat transfer coefficient, based upon the bulk flow properties, have also been set up for the whole body of data encompassing the different inclination angles and aspect ratios. These indicate that the Froude number as a dimensionless gas velocity is a better correlating parameter than the gas Reynolds number. As an alternative approach, a simple dimensionless relationship for the beat transfer coefficient was obtained by analogy between heat and momentum transfer through the interface. Finally, a turbulence-centered model has been modified by using measured interfacial parameters for the turbulent velocity and length scales, resulting in good agreement with the data.     

Electro-magnetic-field-induced flow and interfacial instabilities in confined stratified liquid layers

The electro-magneto-hydrodynamic (EMHD) flow and instabilities engendered by the Lorenz force arising from interaction between externally applied perpendicular electric and magnetic fields are investigated in layers of two immiscible liquids in a channel. A new finite wave-number EMHD instability mode is uncovered by the Orr–Sommerfeld analysis, in addition to the interfacial and shear modes which also arise in the pressure-driven flows. Thus, EMHD can be controlled for micro-channel transport, heat and mass transfer, mixing, micro-emulsion generation, etc.     

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.     

Structural and interfacial stability of multiple solutions for stratified flow

Prediction of the liquid level in stratified two-phase upwards flow shows that one may have multiple solutions. In this case it is necessary to determine which solutions will actually occur and whether hysteresis is possible, namely whether it is possible to have two or more solutions for the same operating conditions. In this work the stability of the solutions for stratified flow is considered using two types of stability analyses: (1) structural stability analysis; and (2) interfacial stability analysis (Kelvin—Helmholtz, K—H). For the K—H stability analysis we used two methods: an approximate simplified method suggested by Taitel & Dukler; and a more rigorous method suggested by Barnea, which is based on a combination of the viscous K—H and inviscid K—H analyses. The results show that whenever three solutions exist only the first, i.e. the solution with the thinnest liquid level, is stable. The middle solution is always structurally unstable (linearly), whereas the third solution is structurally unstable to large disturbances (non-linear stability). The third solution is usually also unstable to the K—H type of instability. As a result it is concluded that hysteresis is not possible and that only the thinnest solution will be observed practically.     

Prediction of liquid level distribution in horizontal gas-liquid stratified flows with interfacial level gradient

A one-dimensional momentum equation has been derived based on a two-fluid model and used to predict the axial distribution of liquid level or void fraction in steady, cocurrent, gas-liquid stratified flows in horizontal circular pipes and rectangular channels. The equation is carefully formulated to incorporate the effect of interfacial level gradient. Two different critical liquid levels are found from the momentum equation and are adopted as a boundary condition to calculate the liquid level or void fraction distribution upstream of the channel exit. The predicted void fraction distributions are compared with the experimental data obtained in a rectangular channel in this work and other data reported for large-diameter pipes. Good agreement is shown for air-kerosene, air-water and stream-water stratified flows with a smooth gas-liquid interface.     

The stratified flow of gas and non-newtonian liquid in horizontal pipes

Predictions of pressure drop and holdup are presented for the stratified flow of gas and non-Newtonian liquid obeying the Ostwald-de Waele power law model. The model of Taitel & Dukler (1976) for gas/Newtonian liquid flow is extended to liquids possessing either shear-thinning or shear-thickening laminar flow behaviour and computed results are given for flow behaviour indices in the range $\text{0.1 ≤ n ≤ 2}$. In particular, conditions are defined for drag reduction of the liquid flow by the presence of the gas. It is concluded that drag reduction occurs over the largest ranges of liquid and gas flow rates at the lowest $\text{n}$ values, provided that liquid flow remains laminar, but that maximum drag reduction may be expected for shear-thickening liquids with $\text{n}$ values of 2 or greater. Ratios of the liquid flow rate in the presence of gas to that for liquid flow alone under a constant pressure gradient are also presented. These ratios frequently exceed unity and are greatest for highly shear-thinning liquids.Although the Taitel & Dukler approach is consistent with experiments on gas/Newtonian liquid flow, and, in addition, appears to be valid for immiscible Newtonian liquid-liquid systems, provided that the viscosity ratio of the two phases is at least five, experiments are required to confirm its applicability for gas/non-Newtonian systems.     

An investigation into the interfacial shear stress contribution in two-phase stratified flow

This paper provides a combined theoretical and experimental investigation into the contribution of interfacial shear stress in certain co and counter-current flows in circular pipes. Based on momentum balance two equations were developed for such flows then two fluid systems of significantly different density ratio were experimentally tested to quantify these equations.     

Non-newtonian liquid-air stratified flow through horizontal tubes—ii

Non-Newtonian liquid gas stratified flow data were obtained using 0.052 and 0.025 m dia horizontal circular ducts. Unless the liquid velocity was very low, the flow pattern generally observed was non-uniform stratified flow having an interfacial level gradient between the two phases. The Heywood-Charles model is valid for predicting the pressure drop and liquid holdup in pseudoplastic (shear thinning) non-Newtonian liquid-gas uniform stratified flow. Two-phase drag reduction, which is predicted by the Heywood-Charles model did not occur because there was a transition to semi-slug flow before the model criteria were reached. Interfacial liquid and gas shear stresses were compared.     

Numerical analysis for slugging of steam-water stratified two-phase flow in horizontal duct

Thermalhydraulic transient phenomena of a steam-water two-phase flow was calculated numerically in order to investigate the onset of slugging from a stratified flow in a horizontal duct. Conservation equations were solved by the finite difference method using a two-phase flow analyzer ‘MINCS’. The analysis was performed to investigate the initiation of slugging with and without phase change, or condensation. The present instability criteria for the onset of slugging with no condensation agreed well with that of the Mishima–Ishii relation while it was much lower than that defined by the Kelvin–Helmholtz instability criteria. However, as the temperature difference between phases increased, steam velocity became higher for the onset of slugging condition. The characteristics of flow reversal and water hammering which were the consequences of slugging with condensation, were investigated and described. It is expected that this modeling could be well applied to complicated thermalhydraulic phenomena accompanied by flow reversal and water hammering in power plants.