Flow Rate Measurement in Multiphase Flow Using Turbine Flow-meter

对涡轮流量传感器进行了理论分析,给出了涡轮流量计仪表常数的计算方法,讨论了获得较大固有仪表常数K_0时涡轮传感器结构参数（如叶片数、涡轮半径、口径等）的优化组合问题,通过多相流动实验,总结出K_0与流动密度之间的实验关系,由此给出用涡轮流量计测量多相流的半理论半经验公式,并在油井多相流量测量中得到了实际应用,符合较好。

Multiphase Flow Measurement in the Offshore Oil Industry

In recent years, considerable research has been conducted into the development of a three-phase flowmeter suitable for use in an offshore environment, and oil/gas/water three-phase metering becomes an important aspect in multiphase flow measurement. This paper discusses the importance of three-phase flow measurement in offshore oil industry, describes the current development in this area, and points out the principal strategies which may be used to meter three-phase flow.

A coupled model for multiphase fluid flow and sedimentation deformation in oil reservoir and its numerical simulation

A mathematical model for coupled multiphase fluid flow and sedimentation deformation is developed based on fluid-solid interaction mechanism. A finite difference-finite element numerical approach is presented. The results of an example show that the fluid-solid coupled effect has great influence on multiphase fluid flow and reservoir recovery performances, and the coupled model has practical significance for oilfield development.

Experimental study and simulation principles of an oil-gas multiphase transportation system

Presented is an experimental study on the performance of an oil-gas multiphase transportation system, especially on the multiphase flow patterns, multiphase pumping and multiphase metering of the system. A dynamic simulation analysis is conducted to deduce simulation parameters of the system and similarity criteria under simplified conditions are obtained. The reliability and feasibility of two-phase flow experiment with oil and natural gas simulated by water and air are discussed by using the similarity criteria.

The kinetic-theory for dilute solid liquid 2-phase flow

On the basis of a brief review of the continuum theory for macroscopic descriptions and the kinetic theory for microscopic descriptions in solid/liquid two-phase flows, some suggestions are presented, i.e. the solid phase may be described by the Boltzmann equation and the liquid phase still be described by conservation laws in the continuum theory. Among them the action force on the particles by the liquid fluid is a coupling factor which connects the phases. For dilute steady solid/liquid two-phase flows, the particle velocity distribution function can be derived by analogy with the procedures in the kinetic theory of gas molecules for the equilibrium state instead of being assumed, as previous investigators did. This done, more detailed information, such as the velocity probability density distribution, mean velocity distribution and fluctuating intensity etc. can be obtained directly from the particle velocity distribution function or from its integration. Experiments have been performed for dilute solid/liquid two-phase flow in a 4 x 6 cm2 sized circulating square pipe system by means of laser Doppler anemometry so that the theories can be examined. The comparisons show that the theories agree very well with all the measured data.

Numerical study on heavy rigid particle motion in a plane wake flow by spectral element method

Experimental particle dispersion patterns in a plane wake flow at a high Reynolds number have been predicted numerically by discrete vortex method (Phys. Fluids A 1992; 4:2244-2251; Int. J. Multiphase Flow 2000; 26:1583-1607). To address the particle motion at a moderate Reynolds number, spectral element method is employed to provide an instantaneous wake flow field for particle dynamics equations, which are solved to make a detail classification of the patterns in relation to the Stokes and Froude numbers. It is found that particle motion features only depend on the Stokes number at a high Froude number and depend on both numbers at a low Froude number. A ratio of the Stokes number to squared Froude number is introduced and threshold values of this parameter are evaluated that delineate the different regions of particle behavior. The parameter describes approximately the gravitational settling velocity divided by the characteristic velocity of wake flow. In order to present effects of particle density but preserve rigid sphere, hollow sphere particle dynamics in the plane wake flow is investigated. The evolution of hollow particle motion patterns for the increase of equivalent particle density corresponds to that of solid particle motion patterns for the decrease of particle size. Although the thresholds change a little, the parameter can still make a good qualitative classification of particle motion patterns as the inner diameter changes.

Two-phase flow and pool boiling heat transfer in microgravity

Researches on two-phase flow and pool boiling heat transfer in microgravity, which included groundbased tests, flight experiments, and theoretical analyses, were conducted in the National Microgravity Laboratory/CAS. A semi-theoretical Weber number model was proposed to predict the slug-to-annular flow transition of two-phase gas–liquid flows in microgravity, while the influence of the initial bubble size on the bubble-to-slug flow transition was investigated numerically using the Monte Carlo method. Two-phase flow pattern maps in microgravity were obtained in the experiments both aboard the Russian space station Mir and aboard IL-76 reduced gravity airplane. Mini-scale modeling was also used to simulate the behavior of microgravity two-phase flow on the ground. Pressure drops of two-phase flow in microgravity were also measured experimentally and correlated successfully based on its characteristics. Two space experiments on pool boiling phenomena in microgravity were performed aboard the Chinese recoverable satellites. Steady pool boiling of R113 on a thin wire with a temperature-controlled heating method was studied aboard RS-22, while quasi-steady pool boiling of FC-72 on a plate was studied aboard SJ-8. Ground-based experiments were also performed both in normal gravity and in short-term microgravity in the drop tower Beijing. Only slight enhancement of heat transfer was observed in the wire case, while enhancement in low heat flux and deterioration in high heat flux were observed in the plate case. Lateral motions of vapor bubbles were observed before their departure in microgravity. The relationship between bubble behavior and heat transfer on plate was analyzed. A semi-theoretical model was also proposed for predicting the bubble departure diameter during pool boiling on wires. The results obtained here are intended to become a powerful aid for further investigation in the present discipline and development of two-phase systems for space applications.

Seed bubbles stabilize flow and heat transfer in parallel microchannels

Seed bubbles are generated on microheaters located at the microchannel upstream and driven by a pulse voltage signal, to improve flow and heat transfer performance in microchannels. The present study investigates how seed bubbles stabilize flow and heat transfer in micro-boiling systems. For the forced convection flow, when heat flux at the wall surface is continuously increased, flow instability is self-sustained in microchannels with large oscillation amplitudes and long periods. Introduction of seed bubbles in time sequence improves flow and heat transfer performance significantly. Low frequency (similar to 10 Hz) seed bubbles not only decrease oscillation amplitudes of pressure drops, fluid inlet and outlet temperatures and heating surface temperatures, but also shorten oscillation cycle periods. High frequency (similar to 100 Hz or high) seed bubbles completely suppress the flow instability and the heat transfer system displays stable parameters of pressure drops, fluid inlet and outlet temperatures and heating surface temperatures. Flow visualizations show that a quasi-stable boundary interface from spheric bubble to elongated bubble is maintained in a very narrow distance range at any time. The seed bubble technique almost does not increase the pressure drop across microsystems, which is thoroughly different from those reported in the literature. The higher the seed bubble frequency, the more decreased heating surface temperatures are. A saturation seed bubble frequency of 1000-2000 Hz can be reached, at which heat transfer enhancement attains the maximum degree, inferring a complete thermal equilibrium of vapor and liquid phases in microchannels. Benefits of the seed bubble technique are the stabilization of flow and heat transfer, decreasing heating surface temperatures and improving temperature uniformity of the heating surface.

Multi-channel effect of condensation flow in a micro triple-channel condenser

Multi-channel effect is important to understand transport phenomenon in phase change systems with parallel channels. In this paper, visualization studies were performed to study the multi-channel effect in a silicon triple-channel condenser with an aspect ratio of 0.04. Saturated water vapor was pumped into the microcondenser, which was horizontally positioned. The condenser was cooled by the air natural convention heat transfer in the air environment. Flow patterns are either the annular flow at high inlet vapor pressures, or a quasi-stable elongated bubble at the microchannel upstream followed by a detaching or detached miniature bubble at smaller inlet vapor pressures. The downstream miniature bubble was detached from the elongated bubble tip induced by the maximum Weber number there. It is observed that either a single vapor thread or dual vapor threads are at the front of the elongated bubble. A miniature bubble is fully formed by breaking up the vapor thread or threads. The transient vapor thread formation and breakup process is exactly symmetry against the centerline of the center channel. In side channels, the Marangoni effect induced by the small temperature variation over the channel width direction causes the vapor thread formation and breakup process deviating from the side channel centerline and approaching the center channel. The Marangoni effect further forces the detached bubble to rotate and approach the center channel, because the center channel always has higher temperatures, indicating the multi-channel effect. 

Investigations of phase inversion and frictional pressure gradients in upward and downward oil–water flow in vertical pipes

The present study has attempted to investigate phase inversion and frictional pressure gradients during simultaneous vertical flow of oil and water two-phase through upward and downward pipes. The liquids selected were white oil (44 mPa s viscosity and 860 kg/m3 density) and water. The measurements were made for phase velocities varying from 0 to 1.24 m/s for water and from 0 to 1.87 m/s for oil, respectively. Experiments were carried either by keeping the mixture velocity constant and increasing the dispersed phase fraction or by keeping the continuous phase superficial velocity constant and increasing the dispersed phase superficial velocity. From the experimental results, it is shown that the frictional pressure gradient reaches to its lower value at the phase inversion point in this work. The points of phase inversion are always close to an input oil fraction of 0.8 for upward flow and of 0.75 for downward flow, respectively. A few models published in the literature are used to predict the phase inversion point and to compare the results with available experimental data. Suitable methods are suggested to predict the critical oil holdup at phase inversion based on the different viscosity ratio ranges. Furthermore, the frictional pressure gradient is analyzed with several suitable theoretical models according to the existing flow patterns. The analysis reveals that both the theoretical curves and the experimental data exhibit the same trend and the overall agreement of predicted values with experimental data is good, especially for a high oil fraction.

Shock wave reflection from a wedge in a dusty gas

The present paper contains a detailed study of shock wave reflection from a wedge placed in various suspensions. In past works, the incident shock propagated initially in pure gas and the suspension started only at the leading edge of the deflecting wedge. However, in the present case the entire flow field is filled with a gas-dust suspension and the initial shock wave has steady-state structure relative to the shock front. In former studies the transmitted shock wave starts its propagation into the suspension and is reflected from the wedge at the same time. It is therefore obvious that the two unrelated processes of (2D) reflection and (1D) "transitional" relaxation occur simultaneously. In the present case the suspension behind the incident shock wave has reached steady state (i.e., it is a traveling wave) before the shock reaches the wedge leading edge. The reflection process from the deflecting wedge is studied for different dust mass loadings and different dust-particle diameter. It is shown that when the dust loading is low and the dust particle diameter is small the wave reflection pattern is similar to that observed in a similar pure gas case. In addition, an equilibrium state is reached, behind the evolved waves, very quickly. On the other hand, when the dust loading is relatively high and/or the dust particle diameter is relatively large, the observed reflection wave pattern is very different from that seen in a similar pure gas case. In such cases it takes much longer time to reach an equilibrium state behind the reflecting waves. It is also shown that the dust presence significantly affects the (gas) pressure on the wedge surface. The higher the dust loading is, the higher the pressure on the wedge surface. Suspensions composed of solid particle of different size, but having the same dust mass loading, will approach the same equilibrium pressure. However, it will take longer time to reach an equilibrium state for suspensions having large diameter particles. (C) 2004 Elsevier Ltd. All rights reserved.

Studies on two-phase co-current air/non-Newtonian shear-thinning fluid flows in inclined smooth pipes

In this work. co-current flow characteristics of air/non-Newtonian liquid systems in inclined smooth pipes are studied experimentally and theoretically using transparent tubes of 20, 40 and 60 turn in diameter. Each tube includes two 10 m lone pipe branches connected by a U-bend that is capable of being inclined to any angle, from a completely horizontal to a fully vertical position. The flow rate of each phase is varied over a wide range. The studied flow phenomena are bubbly, plug flow, slug flow, churn flow and annular flow. These are observed and recorded by a high flow. stratified flow. -speed camera over a wide range of operating conditions. The effects of the liquid phase properties, the inclination angle and the pipe diameter on two-phase flow characteristics are systematically studied. The Heywood-Charles model for horizontal flow was modified to accommodate stratified flow in inclined pipes, taking into account the average void fraction and pressure drop of the mixture flow of a gas/non-Newtonian liquid. The pressure drop gradient model of Taitel and Barnea for a gas/Newtonian liquid slug flow was extended to include liquids possessing shear-thinning flow behaviour in inclined pipes. The comparison of the predicted values with the experimental data shows that the models presented here provide a reasonable estimate of the average void fraction and the corresponding pressure drop for the mixture flow of a gas/ non-Newtonian liquid. (C) 2007 Elsevier Ltd. All rights reserved.

Experimental Study on Subcooled Pool Boiling in Microgravity Utilizing

Turbulent collision of inertial particles: point-particle based, hybrid simulations and beyond

Point-particle based direct numerical simulation (PPDNS) has been a productive research tool for studying both single-particle and particle-pair statistics of inertial particles suspended in a turbulent carrier flow. Here we focus on its use in addressing particle-pair statistics relevant to the quantification of turbulent collision rate of inertial particles. PPDNS is particularly useful as the interaction of particles with small-scale (dissipative) turbulent motion of the carrier flow is mostly relevant. Furthermore, since the particle size may be much smaller than the Kolmogorov length of the background fluid turbulence, a large number of particles are needed to accumulate meaningful pair statistics. Starting from the relative simple Lagrangian tracking of so-called ghost particles, PPDNS has significantly advanced our theoretical understanding of the kinematic formulation of the turbulent geometric collision kernel by providing essential data on dynamic collision kernel, radial relative velocity, and radial distribution function. A recent extension of PPDNS is a hybrid direct numerical simulation (HDNS) approach in which the effect of local hydrodynamic interactions of particles is considered, allowing quantitative assessment of the enhancement of collision efficiency by fluid turbulence. Limitations and open issues in PPDNS and HDNS are discussed. Finally, on-going studies of turbulent collision of inertial particles using large-eddy simulations and particle- resolved simulations are briefly discussed.

Large-eddy simulation of turbulent collision of heavy particles in isotropic turbulence

The small-scale motions relevant to the collision of heavy particles represent a general challenge to the conventional large-eddy simulation (LES) of turbulent particle-laden flows. As a first step toward addressing this challenge, we examine the capability of the LES method with an eddy viscosity subgrid scale (SGS) model to predict the collision-related statistics such as the particle radial distribution function at contact, the radial relative velocity at contact, and the collision rate for a wide range of particle Stokes numbers. Data from direct numerical simulation (DNS) are used as a benchmark to evaluate the LES using both a priori and a posteriori tests. It is shown that, without the SGS motions, LES cannot accurately predict the particle-pair statistics for heavy particles with small and intermediate Stokes numbers, and a large relative error in collision rate up to 60% may arise when the particle Stokes number is near St_K=0.5. The errors from the filtering operation and the SGS model are evaluated separately using the filtered-DNS (FDNS) and LES flow fields. The errors increase with the filter width and have nonmonotonic variations with the particle Stokes numbers. It is concluded that the error due to filtering dominates the overall error in LES for most particle Stokes numbers. It is found that the overall collision rate can be reasonably predicted by both FDNS and LES for St_K>3. Our analysis suggests that, for St_K<3, a particle SGS model must include the effects of SGS motions on the turbulent collision of heavy particles. The spectral analysis of the concentration fields of the particles with different Stokes numbers further demonstrates the important effects of the small-scale motions on the preferential concentration of the particles with small Stokes numbers.