Characteristics of two-phase one-component flow with slip

All loss-of-coolant-accident analysis codes solve the one-dimensional hydraulic conservation equations for a two-phase compressible fluid. Some codes allow slip between the phases, this slip being specified by a slip model. In this paper the characteristics of these conservation equations with well-known slip models (Jones, Moody, Fauske, Smith and CISE) are evaluated as a function of pressure, quality and mass velocity. Regions in which two of the characteristics are complex are shown on contour maps. The significance of these regions is discussed.     

Phase idstribution in horizontal gas-liquid two-phase bubbly flow

An investigation on phase distribution in air-water two-phase flow in horizontal circular channel was conducted by using the double-sensor resistivity probe.The variations of phase distribution with variations of gas and liquid volumetric fluxes were analyzed and the present data were com0pared with some of other researcher‘s data and existing models,It was fund there exists more complicated phase distribution pattern in horizontal flow system than invertical flow.The radial local void fraction profiles are similar at the same measurement angle with various gas and liquid flow rates.However,an asymmetric profile can be observed at a given slice of the pipe cross-section.     

Phase distribution in horizontal gas-liquid two-phase bubbly flow

1IntroductionTwo-phasebubblyflowisencounteredinawidevaxietyofindustrialapplications;Suchasintheproductionofelectricalpowerandpetrochendcajsindustry.Oneofthemosttwortantandyetleastunderstandingaspectsoftwo-phaseflowislateralphasedistribution.Theinformationonphasedistributionisarequiredparameterforhydrodynalincandthermalcalculationsinmanypracticalapplications.Considerableexperimentalandseal-theoreticalresearchhasbeenconductedbymanyresearchersinverticalupwaxdsordownwaxdsbubblynow.[1-7]However,L…     

Internal structure and interfacial velocity development for bubbly two-phase flow

This paper describes an experimental study of the internal structure of air-water flowing horizontally. The double-sensor resistivity probe technique was applied for measurements of local interfacial parameters, including void fraction, interfacial area concentration, bubble size distributions, bubble passing frequency and bubble interface velocity. Bubbly flow patterns at several flow conditions were examined at three axial locations, L/D = 25, 148 and 253, in which the first measurement represents the entrance region where the flow develops, and the second and third may represent near fully developed bubbly flow patterns. The experimental results are presented in three-dimensional perspective plots of the interfacial parameters over the cross-section. These multi-dimensional presentations showed that the local values of the void fraction, interfacial area concentration and bubble passing frequency were nearly constant over the cross-section at L/D = 25, with slight local peaking close to the channel wall. Although similar local peakings were observed at the second and third locations, the internal flow structure segregation due to buoyancy appeared to be very strong in the axial direction. A simple comparison of profiles of the interfacial parameters at the three locations indicated that the flow pattern development was a continuous process. Finally, it was shown that the so-called “fully developed” bubbly two-phase flow pattern cannot be established in a horizontal pipe and that there was no strong correspondence between void fraction and interface velocity profiles.     

Pressure wave phenomena in bubbly liquid

Propagation of shock waves in dilute bubbly liquids is investigated numerically taking into account internal phenomena inside the bubbles. Governing equations for the bubbly liquid are formulated with emphasis on the radial and transverse motions of the bubbles. A numerical method, where individual bubbles are traced to estimate the effects of transverse motions and volumetric changes on the wave phenomena, is developed. Numerical results under several conditions reveal that the radial motion of the bubbles, which is affected by the internal phenomena, such as thermal conduction through the bubble wall, have significant influences on the change of propagation velocity of the shock wave and relaxation phenomena behind the wave. Slippage between bubbles and liquid do not have so much influence on the wave phenomena as the thermal conditions inside the bubble. The nondimensional thermal diffusivity,

is one of the most effective parameters to be correlated with the wave propagation processes, such as the propagation velocity and the waveform change with the wave propagation.     

Interfacial area density, mean radius and number density measurements in bubbly two-phase flow

This paper presents a new experimental method for measuring the local interfacial area density, average bubble radius and number density for bubbly two-phase flow. Experiments were performed using a KfK dual-sensor resistivity probe. The results obtained for interfacial area density using this method agreed with other methods previously proposed, except for high void fractions, close to the transition to slug flow, where prior methods are known to no longer be valid.     

基于图像法的气液两相稀疏泡状流气泡参数分析

采用图像法对垂直上升管气液两相流中稀疏上升气泡进行实验测量研究。该方法使用高速摄像机拍摄气泡运动图像,经图像处理后,提取气泡的特征参数,分别绘制稀疏气泡上升过程中速度变化曲线和单个气泡上升过程中面积变化及形心位置变化曲线,分析气泡参数,总结运动规律。实验结果表明,采用图像法可以很好地完成对气泡参数的分析。     

Study of incident water hammer in an engineering loop under two-phase flow experiment

Water hammer can occur in case of an inflow of sub-cooled water into pipes or other parts of the equipment filled with steam or steam-water mixture. The shock loading due to dynamic pressure during water hammer induces high stresses in the walls of piping especially in the bends. A study was carried out to analyze real life incident of water hammer wherein an attempt was made to recreate water hammer involving steam-water interaction in an instrumented engineering loop to capture pressure, temperature, flow rate and fluid levels in the loop. The arrival time of multiple shocks and the time lag between them is examined in the paper with the support of energy balance equation for variable flow processes.     

Experiment and numerical simulation of bubbly two-phase flow across horizontal and inclined rod bundles

Experimental and numerical analyses were carried out on vertically upward air-water bubbly two-phase flow behavior in both horizontal and inclined rod bundles with either in-line or staggered array. The inclination angle of the rod bundle varied from 0 to 60° with respect to the horizontal. The measured phase distributions indicated non-uniform characteristics, particularly in the direction of the rod axis when the rods were inclined. The mechanisms for this non-uniform phase distribution is supposed to be due to: (1) Bubble segregation phenomenon which depends on the bubble size and shape; (2) bubble entrainment by the large scale secondary flow induced by the pressure gradient in the horizontal direction which crosses the rod bundle; (3) effects of bubble entrapment by vortices generated in the wake behind the rods which travel upward along the rod axis; and (4) effect of bubble entrainment by local flows sliding up along the front surface of the rods. The liquid velocity and turbulence distributions were also measured and discussed. In these speculations, the mechanisms for bubble bouncing at the curved rod surface and turbulence production induced by a bubble were discussed, based on visual observations. Finally, the bubble behaviors in vertically upward bubbly two-phase flow across horizontal rod bundle were analyzed based on a particle tracking method (one-way coupling). The predicted bubble trajectories clearly indicated the bubble entrapment by vortices in the wake region.     

On the development of a model for predicting phase separation phenomena in dividing two-phase flow

A model to predict the dividing flow characteristics for annular flow in a T-junction is proposed consisting of mixture and vapour phase continuity equations, two pressure change correlations and a closure relationship. The pressure change from the inlet through the run of the T is modelled by way of a balance of axial momentum at the junction based on a separated flow assumption. The branch pressure change is modelled using a balance of mechanical energy for the branching flow consisting of reversible and irreversible components. The closure relationship links the phase separation characteristics with the junction pressure changes. It involves a balance between pressure and inertia forces within the junction volume defining a dividing surface for each phase between the run and branch flows. The branch quality is then determined using a well-defined inlet flow distribution. The model is capable of predicting the experimentally observed phase separation characteristics from three independent studies of annular/steam—water and air-water flow in dividing T-junctions.     

The measurement of void waves in bubbly two-phase flows

The propagation of void fraction disturbances (i.e. void waves) in bubbly two-phase air-oil flow was investigated experimentally.The void wave data often showed two simultaneous void wave propagations; in particular, the classical kinematic void wave and a faster void wave propagation associated with bubble clusters. Significantly, the bubbly-to-slug flow regime transition was found to be associated with amplification of the void waves associated with the propagating bubble clusters.These data should be valuable for the assessment of the closure laws used in two-fluid models and for the development of mechanistic models for the bubbly-slug flow regime transition.     

Experimental investigation and CFD simulation of horizontal stratified two-phase flow phenomena

For the investigation of stratified two-phase flow, two horizontal channels with rectangular cross-section were built at Forschungszentrum Dresden-Rossendorf (FZD). The channels allow the investigation of air/water co-current flows, especially the slug behaviour, at atmospheric pressure and room temperature. The test-sections are made of acrylic glass, so that optical techniques, like high-speed video observation or particle image velocimetry (PIV), can be applied for measurements. The rectangular cross-section was chosen to provide better observation possibilities. Moreover, dynamic pressure measurements were performed and synchronised with the high-speed camera system.CFD post-test simulations of stratified flows were performed using the code ANSYS CFX. The Euler–Euler two fluid model with the free surface option was applied on grids of minimum 4 × 10⁵ control volumes. The turbulence was modelled separately for each phase using the k–ω-based shear stress transport (SST) turbulence model. The results compare very well in terms of slug formation, velocity, and breaking. The qualitative agreement between calculation and experiment is encouraging and shows that CFD can be a useful tool in studying horizontal two-phase flow.     

Interfacial area concentration in bubbly flow

Interfacial area concentration is one of the most important parameters in two-phase flow. It is related to mass, momentum and energy transfer at the interface. In bubbly flow, it has close relation to bubble diameter. Interfacial area concentration was measured with particular attention to bubble diameter at inlet. The relationship between bubble diameter and interfacial area was theoretically considered. The possibilities of measurement methods of interfacial area concentration with multiprobes were discussed and a preliminary experiment for one of these methods was carried out.     

CFD studies on the phenomena around counter-current flow limitations of gas/liquid two-phase flow in a model of a PWR hot leg

In order to improve the understanding of counter-current two-phase flow and to validate new physical models, CFD simulations of a 1/3rd scale model of the hot leg of a German Konvoi pressurized water reactor (PWR) with rectangular cross section were performed. Selected counter-current flow limitation (CCFL) experiments conducted at Helmholtz-Zentrum Dresden-Rossendorf (HZDR) were calculated with ANSYS CFX using the multi-fluid Euler–Euler modelling approach. The transient calculations were carried out using a gas/liquid inhomogeneous multiphase flow model coupled with a shear stress transport (SST) turbulence model.In the simulation, the drag law was approached by a newly developed correlation of the drag coefficient (Höhne and Vallée, 2010) in the Algebraic Interfacial Area Density (AIAD) model. The model can distinguish the 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 indicate also a quantitative agreement between calculations and experimental data for the CCFL characteristics and the water level inside the hot leg channel.     

Liquid mean velocity and turbulence in a horizontal air-water bubbly flow

The liquid phase turbulent structure of an air-water bubbly horizontal flow in a circular pipe has been investigated experimentally. Three-dimensional measurements were implemented with two "X" type probes oriented in different planes, and local liquid-phase velocities and turbulent stresses were simultaneously obtained. Systematic measurements were conducted covering a range of local void fraction from 0 to 11.7%. The important experiment results and parametric trends are summarized and discussed.     

Numerical simulation of lateral phase distribution in turbulent upward bubbly two-phase flows

New constitutive models for the interfacial forces acting on bubbles were developed for accurately predicting the lateral phase distribution in turbulent bubbly two-phase flow in vertical channels. Several experimental measurements have revealed that the lateral void profile in bubbly two-phase flow varies from the void peaking near the wall to the almost flat distributions as the liquid velocity increases. However, within the authors' knowledge, the effect of liquid velocity on the void profile has not been successfully predicted by the existing models; this would indicate the strong limitation of the existing multidimensional two-phase flow models. In view of these, the validity of the present constitutive models was tested in varied conditions of the liquid velocity as well as the bubble size. Since several assumptions were required in the models mainly due to the insufficient knowledge of the bubble motion, further improvements should still be needed. Nevertheless, the predicted lateral phase distributions were found to be in reasonably good agreement with available experimental data. It is hence expected that the present constitutive models can effectively be used in the practical applications and also be the base of the more sophisticated ones.     

Numerical simulation of secondary flow in bubbly turbulent flow in sub-channel

Secondary flow in bubbly turbulent flow in sub-channel was simulated by using an algebraic turbulence stress model. The mass, momentum, turbulence energy and bubble diffusion equations were used as fundamental equation. The basis for these equations was the two-fluid model: the equation of liquid phase was picked up from the equation system theoretically derived for the gas-liquid two-fluid turbulent flow. The fundamental equation was transformed onto a generalized coordinate system fitted to the computational domain in sub-channel. It was discretized for the SIMPLE algorithm using the finite-volume method. The shape of sub-channel causes a distortion of the computational mesh, and orthogonal nature of the mesh is sometimes broken. An iterative method to satisfy a requirement for the contra-variant velocity was introduced to represent accurate symmetric boundary condition. Two-phase flow at a steady state was simulated for different magnitudes of secondary flow and void fraction. The secondary flow enhanced the momentum transport in sub-channel and accelerated the liquid phase in the rod gap. This effect was slightly mitigated when the void fraction increased. The acceleration can contribute to effective cooling in the rod gap. The numerical result implied a phenomenon of industrial interest. This suggested that experimental approach is necessary to validate the numerical model and to identify the phenomenon.