EXPERIMENTAL STUDY OF AIR-WATER TWO-PHASE FLOW IN PARALLEL HELICALLY COILED PIPES
Panella, Bruno
2012-01-01
The air-water two-phase flow in a 12 mm inner diameter parallel helically coiled pipes is investigated with three different coils diameters. Void fraction, flow rate distribution and two-phase pressure drops along the pipes in the parallel channels are measured. The test two-phase pressure drops are compared with theoretical ones, in terms of multipliers and friction factors. The instabilities arisen during the experimental tests are investigated and are related to the void fraction and flow ...
Two-phase air-water stratified flow measurement using ultrasonic techniques
Fan, Shiwei; Yan, Tinghu; Yeung, Hoi
2014-04-01
In this paper, a time resolved ultrasound system was developed for investigating two-phase air-water stratified flow. The hardware of the system includes a pulsed wave transducer, a pulser/receiver, and a digital oscilloscope. The time domain cross correlation method is used to calculate the velocity profile along ultrasonic beam. The system is able to provide velocities with spatial resolution of around 1mm and the temporal resolution of 200μs. Experiments were carried out on single phase water flow and two-phase air-water stratified flow. For single phase water flow, the flow rates from ultrasound system were compared with those from electromagnetic flow (EM) meter, which showed good agreement. Then, the experiments were conducted on two-phase air-water stratified flow and the results were given. Compared with liquid height measurement from conductance probe, it indicated that the measured velocities were explainable.
Numerical simulation of air-water two-phase flow over stepped spillways
CHENG; Xiangju; CHEN; Yongcan
2006-01-01
Stepped spillways for significant energy dissipation along the chute have gained interest and popularity among researchers and dam engineers. Due to the complexity of air-water two-phase flow over stepped spillways, the finite volume computational fluid dynamics module of the FLUENT software was used to simulate the main characteristics of the flow. Adopting the RNG k-ε turbulence model, the mixture flow model for air-water two-phase flow was used to simulate the flow field over stepped spillway with the PISO arithmetic technique. The numerical result successfully reproduced the complex flow over a stepped spillway of an experiment case, including the interaction between entrained air bubbles and cavity recirculation in the skimming flow regime, velocity distribution and the pressure profiles on the step surface as well. The result is helpful for understanding the detailed information about energy dissipation over stepped spillways.
Developments in the research of air-water two-phase flows in turbomachinery
Recently, engineering problems associated with two-phase flows in turbomachinery have become increasingly important in relation to the safety analysis of nuclear reactors or the usage of low quality energy resources; the research on this subject has been promoted. It is a really knotty problem caused by the multiform flow patterns as well as the variety of its applications. However, the mechanics in two-phase machines may involve similar phenomena. In this paper, developments of the research of air-water mixtures in turbomachinery will be briefly reviewed, and the mechanics of two-phase flows in rotating flow fields and the prediction methods of the performance of turbomachinery based on some analytical models are discussed. (author)
Air-water two-phase flow pressure drop across various components of AHWR fuel bundle
Single-phase (water) and two-phase (air-water) experiments were carried out for the measurement of pressure drops across various components of a prototype full scale 54-rod fuel bundle of proposed AHWR (Advanced Heavy Water Reactor). From the measured values of pressure drops, the friction factor for fuel bundle and the loss coefficients for the tie plates and spacers were estimated. The single-phase experimental data were compared with different existing correlations. Correlations have been proposed based on the data generated with the air-water mixture which can be used for prediction of pressure drop across fuel channel (with 54 rod fuel bundle) of AHWR under normal operating conditions with appropriate correction factor for steam-water flow. Also a heuristic approach to predict the Lockhart-Martinelli parameter has been presented. Further, a new correlation for two-phase friction multiplier applicable to 54-rod cluster geometry has been developed based on two-phase experimental pressure drop data. The effect of mixture mass flux on the two-phase friction multiplier has been probed and the assessment of existing friction multiplier correlations has also been carried out with the test data. (author)
Interfacial structures of confined air-water two-phase bubbly flow
Kim, S.; Ishii, M.; Wu, Q.; McCreary, D.; Beus, S.G.
2000-08-01
The interfacial structure of the two-phase flows is of great importance in view of theoretical modeling and practical applications. In the present study, the focus is made on obtaining detailed local two-phase parameters in the air-water bubbly flow in a rectangular vertical duct using the double-sensor conductivity probe. The characteristic wall-peak is observed in the profiles of the interracial area concentration and the void fraction. The development of the interfacial area concentration along the axial direction of the flow is studied in view of the interfacial area transport and bubble interactions. The experimental data is compared with the drift flux model with C{sub 0} = 1.35.
Structure of air-water two-phase flow in helically coiled tubes
Air-water two-phase flow in helically coiled tubes is investigated experimentally to elucidate the effects of centrifugal acceleration on the flow regime map and the spatial and the temporal flow structure distribution. Three kinds of test tubes with 20 mm inner diameters including a straight tube are used to compare the turbulent flow structure. Superficial velocities up to 6 m/s are tested so that the centrifugal Froude number covers a range from 0 to 3. The interfacial structure is photographed from two directions by a high-speed video system with synchronized measurement of local pressure fluctuations. The results reveal that the flow transition line alters due to centrifugal force acting on the liquid phase in the tube. In particular, the bubbly flow regime is narrowed significantly. The pressure fluctuation amplitude gets large relatively to the average pressure loss as void fraction increases. The frequency spectra of the pressure fluctuation have plural peaks in the case of strong curvature, implying that the periodicity of slugging two-phase flow is collapsed by an internal secondary flow activated inside the liquid phase. Moreover, under large Froude number conditions, the substantial velocity of the gas phase that biases to the inner side of the helical coil is slower than the total superficial velocity because the liquid flow is allowed to pass through the outer side and so resembles a radial stratified flow
Time-resolved Fast Neutron Radiography of Air-water Two-phase Flows
Zboray, Robert; Dangendorf, Volker; Mor, Ilan; Tittelmeier, Kai; Bromberger, Benjamin; Prasser, Horst-Michael
Neutron imaging, in general, is a useful technique for visualizing low-Z materials (such as water or plastics) obscured by high-Z materials. However, when significant amounts of both materials are present and full-bodied samples have to be examined, cold and thermal neutrons rapidly reach their applicability limit as the samples become opaque. In such cases one can benefit from the high penetrating power of fast neutrons. In this work we demonstrate the feasibility of time-resolved, fast neutron radiography of generic air-water two-phase flows in a 1.5 cm thick flow channel with Aluminum walls and rectangular cross section. The experiments have been carried out at the high-intensity, white-beam facility of the Physikalisch-Technische Bundesanstalt, Germany. Exposure times down to 3.33 ms have been achieved at reasonable image quality and acceptable motion artifacts. Different two-phase flow regimes such as bubbly slug and churn flows have been examined. Two-phase flow parameters like the volumetric gas fraction, bubble size and bubble velocities have been measured.
Experimental study of flow monitoring instruments in air-water, two-phase downflow
The performance of a turbine meter, target flow meter (drag disk), and a gamma densitometer was studied in air-water, two-phase vertical downflow. Air and water were metered into an 0.0889-m-ID (3.5-in.) piping system; air flows ranged from 0.007 to 0.3 m3/sec (16 to 500 scfm) and water flows ranged from 0.0006 to 0.03 m3/sec (10 to 500 gpm). The study included effects of flow rate, quality, flow regime, and flow dispersion on the mean and fluctuating components of the instrument signals. Wire screen flow dispersers located at the inlet to the test section had a significant effect on the readings of the drag disk and gamma densitometer, but had little effect on the turbine. Further, when flow dispersers were used, mass flow rates determined from the three instrument readings and a two-velocity, slip flow model showed good agreement with actual mass flow rate over a three-fold range in quality; mass flows determined with the drag disk and densitometer readings assuming homogeneous flow were nearly as accurate. However, when mass flows were calculated using the turbine and densitometer or turbine and drag disk readings assuming homogeneous flow, results were scattered and relatively inaccurate compared to the actual mass flows. Turbine meter data were used with a two-velocity turbine model and continuity relationships for each phase to determine the void fraction and mean phase velocities in the test section. The void fraction was compared with single beam gamma densitometer results and fluid momentum calculated from a two-velocity model was compared with drag disk readings
An ultrasonic methodology is proposed for the measurement for two phase (air-water) flow parameters. Ultrasonic backscattered signals were used to analyze the following parameters: average number of bubbles, interfacial area and void fraction. The results show a strong correlation between the parameters and the ultrasonic power signal obtained. (author)
Mehta, Hemant B.; Jyotirmay Banerjee
2014-01-01
The experimentally developed flow pattern maps for micro-scale channels reported by various researchers differ significantly. Also, no theoretical models effectively predict the flow regime transition boundaries in micro-scale channel. The present work proposes an empirical model for air-water two-phase flow pattern transition boundaries for minichannel diameters between 2 to 5mm. Moreover, experiments are conducted with 2.5 mm diameter horizontal circular minichannel to develop a flow regime...
Kim, T. H.; Yun, B. J.; Jeong, J. H. [Pusan National University, Geunjeong-gu, Busan (Korea, Republic of)
2015-05-15
Studies were mostly about flow in upward flow in medium size circular tube. Although there are great differences between upward and downward flow, studies on vertical upward flow are much more active than those on vertical downward flow in a channel. In addition, due to the increase of surface forces and friction pressure drop, the pattern of gas-liquid two-phase flow bounded to the gap of inside the rectangular channel is different from that in a tube. The downward flow in a rectangular channel is universally applicable to cool the plate type nuclear fuel in research reactor. The sub-channel of the plate type nuclear fuel is designed with a few millimeters. Downward air-water two-phase flow in vertical rectangular channel was experimentally observed. The depth, width, and length of the rectangular channel is 2.35 mm, 66.7 mm, and 780 mm, respectively. The test section consists of transparent acrylic plates confined within a stainless steel frame. The flow patterns of the downward flow in high liquid velocity appeared to be similar to those observed in previous studies with upward flow. In downward flow, the transition lines for bubbly-slug and slug-churn flow shift to left in the flow regime map constructed with abscissa of the superficial gas velocity and ordinate of the superficial liquid velocity. The flow patterns observed with downward flow at low liquid velocity are different from those with upward flow.
Two-phase air-water flows:Scale effects in physical modeling
PFISTER Michael; CHANSON Hubert
2014-01-01
Physical modeling represents probably the oldest design tool in hydraulic engineering together with analytical approaches. In free surface flows, the similitude based upon a Froude similarity allows for a correct representation of the dominant forces, namely gravity and inertia. As a result fluid flow properties such as the capillary forces and the viscous forces might be incorrectly reproduced, affecting the air entrainment and transport capacity of a high-speed model flow. Small physical models operating under a Froude similitude systematically underestimate the air entrainment rate and air-water interfacial properties. To limit scale effects, minimal values of Reynolds or Weber number have to be respected. The present article summarizes the physical background of such limitations and their combination in terms of the Morton number. Based upon a literature review, the existing limits are presented and discussed, resulting in a series of more conservative recommendations in terms of air concentration scaling. For other air-water flow parameters, the selection of the criteria to assess scale effects is critical because some parameters (e.g., bubble sizes, turbulent scales) can be affected by scale effects, even in relatively large laboratory models.
Air-water two-phase vertical upward flow regime identification with cross-sectional visualization
The paper presents flow regime identification in two-phase vertical pipe with 30 mm diameter using a high-speed camera and a wire-mesh sensor. In experiments, the wire-mesh sensor which consists of 16x16 wires is installed 2800 mm away from the entrance for capturing the local instantaneous phase distribution, and the high-speed camera which can monitor two-phase flow at the maximum 500,000 fps is set up in order to obtain the axial images at that time. The experimental observations for flow regime identification are compared with existing semi-theoretical model. The comparison results show a good agreement between them. (author)
Hemant B. Mehta
2014-01-01
Full Text Available The experimentally developed flow pattern maps for micro-scale channels reported by various researchers differ significantly. Also, no theoretical models effectively predict the flow regime transition boundaries in micro-scale channel. The present work proposes an empirical model for air-water two-phase flow pattern transition boundaries for minichannel diameters between 2 to 5mm. Moreover, experiments are conducted with 2.5 mm diameter horizontal circular minichannel to develop a flow regime map. The proposed empirical model is found to provide good agreement with the experimental data. Comparisons are also shown with the work of Mandhane et al. (1974, Taitel and Dukler (1976, Barnea et al. (1983, Damianides and Westwater (1988, Coleman and Garimella (1999, Yang and Shieh (2001, Venkatesan et al. (2010.
Two-phase air-water flows: Scale effects in physical modelling
Pfister, Michael; Chanson, Hubert
2014-01-01
Physical modeling represents probably the oldest design tool in hydraulic engineering together with analytical approaches. In free surface flows, the similitude based upon a Froude similarity allows for a correct representation of the dominant forces, namely gravity and inertia. As a result fluid flow properties such as the capillary forces and the viscous forces might be incorrectly reproduced, affecting the air entrainment and transport capacity of a high-speed model flow. Small physical mo...
Droplet velocity in horizontal two-phase air/water free jet flow in stagnant ambient
Surma, R.; Friedel, L. [Wissenschaftlicher Arbeitsbereich Stroemungsmechanik, TU Hamburg-Harburg, Eissendorferstr. 38, D-21073 Hamburg (Germany)
2004-08-01
For the case of an isothermal free jet flow the velocity and the size of the drops were measured with a Phase-Doppler-Anemometer. The axial velocity decay downstream from the nozzle outlet and the radial drop velocity profile are reproducible within reasonable limits, if a so-called homogeneous flow mixture density in the nozzle outlet is introduced. In case of low mass flow qualities an adjustment of the entrainment coefficient and a virtual jet origin as well as a reduced exponent for the Gaussian function type radial velocity profile is required. (Abstract Copyright [2004], Wiley Periodicals, Inc.)
The slug flow transitions and related phenomena for horizontal two-phase flows were studied for a better prediction of two-phase flows that typically appear during the reactor loss-of-coolant accidents (LOCAs). For better representation of the flow conditions experimentally, two large-scaled facility: TPTF for high-pressure steam/water two-phase flows and large duct test facility for air/water two-phase flows, were used. The visual observation of the flow using a video-probe was performed in the TPTF experiments for good understanding of the phenomena. The currently-used models and correlations based mostly on the small-scale low-pressure experiments were reviewed and improved based on these experimental results. The modified Taitel-Dukler model for prediction of transition into slug flow from wavy flow and the modified Steen-Wallis correlation for prediction of onset of liquid entrainment from the interfacial waves were obtained. An empirical correlation for the gas-liquid interfacial friction factor was obtained further for prediction of liquid levels at wavy flow. The region of slug flow regime that is generally under influences of the channel height and system pressure was predicted well when these models and correlations were applied together. (author). 90 refs
Numerical simulation of slug flow regime for an air water two-phase flow in horizontal pipes
Slug flow is a quite common multiphase flow regime in horizontal pipelines and channels, which can be potentially hazardous to the structure of the pipe system or to apparatus and processes following the slug flow pipe section due to the strong oscillating pressure levels formed behind liquid slugs. Areas of application are in the chemical and process industry as well as in safety research and thermo-hydraulic engineering for nuclear power plants. The intended paper deals with the feasibility and accuracy of CFD simulations for an air-water slug flow in a horizontal circular pipe of diameter D = 0.054 m and a pipe length of up to 8 m. In the past most investigations of the slug flow regime in horizontal pipelines and channels have been carried out on experimental test rigs. Due to the transient and three-dimensional character of slug flow regime and the resulting numerical effort only a few attempts of numerical simulation have been made. In principal three different computational approaches can be applied for the simulation of horizontal slug flows: - 'frozen slug' in a domain with moving wall boundaries, where the absolute value of the prescribed wall velocity is equal to the slug propagation velocity in the pipe. The slug propagation velocity and the slug length/period has to be known in advance. - Transient 3-D simulation in a short computational domain with periodic boundary conditions. A driving pressure force has to be prescribed to compensate the kinetic energy losses due to wall friction. Furthermore it has to be ensured, that the geometrical dimensions of the computational domain do not affect the computed slug flow length and time scales. - Transient 3. simulation of slug flow in a long pipe segment with inlet/outlet boundary conditions. The later of the three computational approaches provides the highest predictive capability, also it is the most computational intensive approach. The presented paper will discuss the general aspects of feasibility
To predict the behavior of air-water two-phase flows in centrifugal pumps, we have proposed a three-dimensional numerical method on the basis of an inviscid bubbly flow model with slippage between two phases. The void fractions calculated distribute unevenly and their maximum exceeds an applicability of the model. To extend its applicability, a newly modified model is proposed in this paper by assuming that the bubbles in such high void fraction regions coalesce with each other and adhere to the neighboring impeller walls so as to form a fixed cavity. Using this model, the flows in a radial-flow pump are solved. The cavity obtained increases progressively from the shroud to the hub in the section just after the impeller inlet when the inlet void fraction exceeds a critical value and finally fills the section, showing close relation with the experiments when the pump loses its function due to an air-filled blockade. (author)
To predict the behavior of air-water two-phase flows in a centrifugal pump impeller, a three-dimensional numerical method is proposed based on a bubbly flow model. If it is assumed that the mixtures are homogeneous bubbly flow containing fine bubbles compared with the characteristic length of the impeller channel, then the equations of motion of the mixtures are represented by those of liquid phase and its velocity is expressed as a potential for the quasi-harmonic equation. The equations are solved by use of the finite element method to obtain the velocities and pressures, and the equation of motion of an air bubble is integrated numerically on this flow field to obtain the void fraction. These calculations are repeated until the solutions converge. The results obtained show good agreement with experiments within the range of bubbly flow regime. (author)
Image processing analysis on the air-water slug two-phase flow in a horizontal pipe
Dinaryanto, Okto; Widyatama, Arif; Majid, Akmal Irfan; Deendarlianto, Indarto
2016-06-01
Slug flow is a part of intermittent flow which is avoided in industrial application because of its irregularity and high pressure fluctuation. Those characteristics cause some problems such as internal corrosion and the damage of the pipeline construction. In order to understand the slug characteristics, some of the measurement techniques can be applied such as wire-mesh sensors, CECM, and high speed camera. The present study was aimed to determine slug characteristics by using image processing techniques. Experiment has been carried out in 26 mm i.d. acrylic horizontal pipe with 9 m long. Air-water flow was recorded 5 m from the air-water mixer using high speed video camera. Each of image sequence was processed using MATLAB. There are some steps including image complement, background subtraction, and image filtering that used in this algorithm to produce binary images. Special treatments also were applied to reduce the disturbance effect of dispersed bubble around the bubble. Furthermore, binary images were used to describe bubble contour and calculate slug parameter such as gas slug length, gas slug velocity, and slug frequency. As a result the effect of superficial gas velocity and superficial liquid velocity on the fundamental parameters can be understood. After comparing the results to the previous experimental results, the image processing techniques is a useful and potential technique to explain the slug characteristics.
Local gas- and liquid-phase measurements for air-water two-phase flows in a rectangular channel
Local gas- and liquid-phase measurements of various gas-liquid two-phase flows, including bubbly, cap-bubbly, slug, and churn-turbulent flows, were performed in an acrylic vertical channel with a rectangular cross section of 30 mm x 10 mm and height of 3.0 m. All the measurements were carried out at three measurement elevations along the flow channel, with z/Dh = 9, 72, and 136, respectively, to study the flow development. The gas-phase velocity, void fraction, and bubble number frequency were measured using a double-sensor conductivity probe. A high-speed imaging system was utilized to perform the flow regime visualization and to provide additional quantitative information of the two-phase flow structure. An image processing scheme was developed to obtain the gas-phase velocity, void fraction, Sauter mean diameter, bubble number density, and interfacial area concentration. The liquid-phase velocity and turbulence measurements were conducted using a particle image velocimetry-planar laser-induced fluorescence (PIV-PLIF) system, which enables whole-field and high-resolution data acquisition. An optical phase separation method, which uses fluorescent particles and optical filtration technique, is adopted to extract the velocity information of the liquid phase. An image pre-processing scheme is imposed on the raw PIV images acquired to remove noises due to the presence of bubble residuals and optically distorted particles in the images captured by the PIV-PLIF system. Due to the better light access and less bubble distortion in the narrow rectangular channel, the PIV-PLIF system were able to perform reasonably well in flows of even higher void fractions as compared to the situations with circular pipe test sections. The flow conditions being studied covered various flow regime transitions, void fractions, and liquid-phase flow Reynolds numbers. The obtained experimental data can also be used to validate two-phase CFD results. (author)
Zboray, Robert; Dangendorf, Volker; Mor, Ilan; Bromberger, Benjamin; Tittelmeier, Kai
2015-01-01
In a previous work we have demonstrated the feasibility of high-frame-rate, fast-neutron radiography of generic air-water two-phase flows in a 1.5 cm thick, rectangular flow channel. The experiments have been carried out at the high-intensity, white-beam facility of the Physikalisch-Technische Bundesanstalt, Germany, using an multi-frame, time-resolved detector developed for fast neutron resonance radiography. The results were however not fully optimal and therefore we have decided to modify ...
Analysis of Air-Water Two Phase Flow for K-HERMES-HALF Experiment using RELAP5/MOD3
Park, Rae Joon; Ha, Kwang Soon; Kim, Sang Baik; Hong, Seong Wan [Korea Atomic Energy Research Institute, Daejeon (Korea, Republic of); Heo, Sun [KHNP Nuclear Engineering and Technology Institute, Daejeon (Korea, Republic of)
2011-05-15
The IVR (In-Vessel corium Retention) through the ERVC (External Reactor Vessel Cooling) is known to be an effective means for maintaining the integrity of the reactor pressure vessel during a severe accident in a nuclear power plant. This measure has been adopted in some low-power reactors such as the AP600, AP1000, and the Loviisa nuclear power plants as a design feature, and in the high-power reactor of the APR (Advanced Power Reactor) 1400 as an accident management strategy for severe accident mitigation. As part of a study on two-phase flow in the reactor cavity under external reactor vessel cooling in the APR1400, K-HERMES-HALF experiment (Hydraulic Evaluation of Reactor cooling Mechanism by External Self-induced flow-HALF scale) had performed at KAERI. This large-scale experiment using a half-height and half-sector model of the APR1400 uses the non-heating method of the air injection. In this research, K-HERMES-HALF test results had been evaluated by using RELAP5/MOD3 computer code to observe and evaluate the two-phase natural circulation phenomena through the annulus gap between the outer reactor vessel and the vessel insulation material
Minemura, Kiyoshi; Kinoshita, Katsuhiko [Nagoya Univ. (Japan). School of Informatics and Sciences; Ihara, Masaru; Furukawa, Hironori; Egashira, Kazuyuki [Japan National Oil Corp., Chiba (Japan). Petroleum Engineering Lab.
1995-12-31
To establish the optimum design parameters of offshore oil well centrifugal pumps, which should deliver crude oil containing a large amount of gas, various shapes of pump impeller with different outlet blade angles, locations of leading-edge and numbers of impeller blades as the design parameters were tested with various rotating speeds and suction pressures under air-water two-phase flow conditions. The greater the outlet blade angle, the less the degradation of the pump performance becomes, showing the optimum blade angle approximately equals to 90{degree}.
Zboray, Robert [Paul Scherrer Institute, PSI Villigen 5232 (Switzerland); Dangendorf, Volker; Bromberger, Benjamin; Tittelmeier, Kai [Physikalisch-Technische Bundesanstalt (PTB), Braunschweig 38116 (Germany); Mor, Ilan [Soreq NRC, Yavne 81800 (Israel)
2015-07-15
In a previous work, we have demonstrated the feasibility of high-frame-rate, fast-neutron radiography of generic air-water two-phase flows in a 1.5 cm thick, rectangular flow channel. The experiments have been carried out at the high-intensity, white-beam facility of the Physikalisch-Technische Bundesanstalt, Germany, using an multi-frame, time-resolved detector developed for fast neutron resonance radiography. The results were however not fully optimal and therefore we have decided to modify the detector and optimize it for the given application, which is described in the present work. Furthermore, we managed to improve the image post-processing methodology and the noise suppression. Using the tailored detector and the improved post-processing, significant increase in the image quality and an order of magnitude lower exposure times, down to 3.33 ms, have been achieved with minimized motion artifacts. Similar to the previous study, different two-phase flow regimes such as bubbly slug and churn flows have been examined. The enhanced imaging quality enables an improved prediction of two-phase flow parameters like the instantaneous volumetric gas fraction, bubble size, and bubble velocities. Instantaneous velocity fields around the gas enclosures can also be more robustly predicted using optical flow methods as previously.
Zboray, Robert; Dangendorf, Volker; Mor, Ilan; Bromberger, Benjamin; Tittelmeier, Kai
2015-07-01
In a previous work, we have demonstrated the feasibility of high-frame-rate, fast-neutron radiography of generic air-water two-phase flows in a 1.5 cm thick, rectangular flow channel. The experiments have been carried out at the high-intensity, white-beam facility of the Physikalisch-Technische Bundesanstalt, Germany, using an multi-frame, time-resolved detector developed for fast neutron resonance radiography. The results were however not fully optimal and therefore we have decided to modify the detector and optimize it for the given application, which is described in the present work. Furthermore, we managed to improve the image post-processing methodology and the noise suppression. Using the tailored detector and the improved post-processing, significant increase in the image quality and an order of magnitude lower exposure times, down to 3.33 ms, have been achieved with minimized motion artifacts. Similar to the previous study, different two-phase flow regimes such as bubbly slug and churn flows have been examined. The enhanced imaging quality enables an improved prediction of two-phase flow parameters like the instantaneous volumetric gas fraction, bubble size, and bubble velocities. Instantaneous velocity fields around the gas enclosures can also be more robustly predicted using optical flow methods as previously.
In a previous work, we have demonstrated the feasibility of high-frame-rate, fast-neutron radiography of generic air-water two-phase flows in a 1.5 cm thick, rectangular flow channel. The experiments have been carried out at the high-intensity, white-beam facility of the Physikalisch-Technische Bundesanstalt, Germany, using an multi-frame, time-resolved detector developed for fast neutron resonance radiography. The results were however not fully optimal and therefore we have decided to modify the detector and optimize it for the given application, which is described in the present work. Furthermore, we managed to improve the image post-processing methodology and the noise suppression. Using the tailored detector and the improved post-processing, significant increase in the image quality and an order of magnitude lower exposure times, down to 3.33 ms, have been achieved with minimized motion artifacts. Similar to the previous study, different two-phase flow regimes such as bubbly slug and churn flows have been examined. The enhanced imaging quality enables an improved prediction of two-phase flow parameters like the instantaneous volumetric gas fraction, bubble size, and bubble velocities. Instantaneous velocity fields around the gas enclosures can also be more robustly predicted using optical flow methods as previously
An electromagnetic velocity meter has been applied to measure the liquid velocity distribution along a large vertical pipe (Inner diameter: 0.48 m, Length of flow path: 2 m) under air-water two-phase flow. This measurement is performed to examine the flow structure along the large vertical pipe where the flow structure has not been fully understood yet. The experiment was performed under atmospheric pressure and superficial air and water velocities in the test section were 0.02--0.87 m/s and 0.01--0.2 m/s, respectively. The accuracy of the electromagnetic velocity meter was firstly checked and was confirmed to be within the error of ±10% for the local liquid velocity up to about 2 m/s under bubbly two-phase flow. The velocity meter was used to measure the radial distribution of local liquid velocity including flow direction in the large vertical pipe. With increasing air flow rate, the axial liquid velocity at the center of the pipe becomes higher, the direction of axial liquid flow near the wall becomes downward and the degree of anisotropy of liquid velocity fluctuation becomes larger. A developing region exists below about 1 m from the bottom of the test section and the flow structure above the elevation is considered to be almost developed based on the measurement of the radial distribution of axial liquid velocity
Zboray, Robert; Mor, Ilan; Bromberger, Benjamin; Tittelmeier, Kai
2015-01-01
In a previous work we have demonstrated the feasibility of high-frame-rate, fast-neutron radiography of generic air-water two-phase flows in a 1.5 cm thick, rectangular flow channel. The experiments have been carried out at the high-intensity, white-beam facility of the Physikalisch-Technische Bundesanstalt, Germany, using an multi-frame, time-resolved detector developed for fast neutron resonance radiography. The results were however not fully optimal and therefore we have decided to modify the detector and optimize it for the given application, which is described in the present work. Furthermore, we managed to improve the image post-processing methodology and the noise suppression. Using the tailored detector and the improved post-processing significant increase in the image quality and an order of magnitude lower exposure times, down to 3.33 ms, have been achieved with minimized motion artifacts. Similar to the previous study, different two-phase flow regimes such as bubbly slug and churn flows have been e...
The Effect of Sudden Change in Pipe Diameter on Flow Patterns of Air-Water Two-Phase
无
2001-01-01
Flow patterns upstream and downstream of a sudden-contraction cross-section in a vertical straight pipe were presented. By comparing with flow patterns in uniform croes-section vertical tubes, the effect of the sudden change in pipe diameter on flow patterns was analyzed. Flow pattern transition mechanisms were discussed and transition criteria for flow pattern transitions were deduced accordingly using the dimensional analysis.
Two-phase compressible flow fields of air-water are investigated numerically in the fixed Eulerian grid framework. The phase interface is captured via volume fractions of ech phase. A way to model two phase compressible flows as a single phase one is found based on an equivalent equation of states of Tait's type for a multiphase cell. The equivalent single phase field is discretized using the Roe's approximate Riemann solver. Two approaches are tried to suppress the pressure oscillation phenomena at the phase interface, a passive advection of volume fraction and a direct pressure relaxation with the compressible form of volume fraction equation. The direct pressure equalizing method suppresses pressure oscillation successfully and generates sharp discontinuities, transmitting and reflecting acoustic waves naturally at the phase interface. In discretizing the compressible form of volume fraction equation, phase interfaces are geometrically reconstructed to minimize the numerical diffusion of volume fraction and relevant variables. The motion of a projectile in a water-filled tube which is fired by the release of highly pressurized air is simulated presuming the flow field as a two dimensional one, and several design factors affecting the projectile movement are investigated
Full text of publication follows: Within the framework of the nuclear power plant lifetime issue, the assessment of the French 900 MWe (3-loops) series reactor pressure vessel (RPV) integrity has been performed. A simplified analysis has shown that the most severe loading conditions are given by the small break loss of coolant accidents due to the pressurized injection of cold water (9 deg. C) into the cold leg and down comer of the RPV. During these transient scenarios, single or two-phase (uncovered cold leg) flows have been shown in the cold leg, depending on the crack size and RPV model (900 MWe or 1300 MWe). An experimental study has been carried out, on the one hand, to consolidate the numerical results obtained with CFD home code (Code-Saturne) which mainly showed the stratified flow in the cold leg and the fluid flow separation and its oscillations in the down comer during a single phase scenario. These physical phenomena are important for the thermal RPV loading assessment. On the other hand, the absence of experimental two-phase data necessitated to carry out an experimental study around the mixing area behavior (free surface, stratified flow) during an ECC injection with an uncovered cold leg. The new EDF R and D mock up, called HYBISCUS, is a facility which is made out of Plexiglas (atmosphere pressure) and represents a half scale CP0 geometry with one cold leg and part of the down comer. The mock up modularity allows us to insert representative ECC nozzles and a thermal shield. In reference to the reactor scenarios, the experimental operating conditions are derived from the conservation of the density effects (Froude number). For that, a heated salted water flow is used to represent the ECC injection whereas water represents the cold leg fluid. This mock up has been defined in order to represent single phase flow (cold leg and down comer full of water) or two-phase flow (uncovered cold leg) ECC scenarios. This paper reports experimental results
Experimental investigation about the onset flooding was made on an air-water two-phase flow in inclined pipe. The pipe inclination angles examined were 30deg, 45deg and 60deg from horizontal. The tube diameter was 26 mm I.D., and superficial velocities of air and water ranged from 2.2 - 25 m/sec and from 0.018 - 0.141 m/sec, respectively. As a result, it was found that (1) there were two characteristic regimes before and after the onset of flooding, pre-flooding and post-flooding, (2) unstable movement in the liquid film with roll wave was occurred at the onset of flooding, (3) pressure gradient increased suddenly at the onset of flooding, while it decreased gradually in post-flooding after violent fluctuation, (4) the fluctuation in post-flooding was the biggest among the three regimes, (5) the flooding gas velocity increased with the pipe inclination due to the increment of gravity component which has to be sustained by the gas flow. (author)
To predict the behavior of gas-liquid two-phase flows in a centrifugal pump impeller, a three-dimensional numerical method is proposed on the basis of a bubbly flow model. Under the assumption of homogeneous bubbly flow entraining fine bubbles, the equation of motion of the mixture is represented by that of liquid-phase and the liquid velocity is expressed as a potential for a quasi-harmonic equation. This equation is solved with a finite element method to obtain the velocities, and the equation of motion of an air bubble is integrated numerically in the flow field to obtain the void fraction. These calculations are iterated to obtain a converged solution. The method has been applied to a radial-flow pump, and the results obtained have been confirmed by experiments within the range of bubbly flow regime
Two-phase flow characterisation by nuclear magnetic resonance
The results presented in this paper demonstrate the performance of the PFGSE-NMR to obtain a complete characterisation of two-phase flows. Different methods are proposed to characterise air-water flows in different regimes: stationary two-phase flows and flows in transient condition. Finally a modified PFGSE is proposed to analyse the turbulence of air-water bubbly flow. (author)
无
2001-01-01
Flow patterns in upstream and downstream straight tubes of sudden-changedareas in ahorizontal straight pipe were experimentally examined. Both sudden-expansioncross-section (SECS)and sudden-contraction cross-section (SCCS) were investigated. The flow pattern mapsupstream anddownstream were delineated and compared with those in straight tubes with uniformcross-sections.The effects of the SECS and SCCS on flow patterns were discussed and analyzed.Furthermore, flowpattern transition mechanisms resulting in occurrences of different flow patternswere simplydiscussed and some transition criteria for the flow pattern transitions were deduced byusing the non-dimensionlized analysis method.
刘雪敏; 李舟航; 吴玉新; 吕俊复
2012-01-01
The air-water two phase flow pattern in vertical tube was experimentally investigated with different tube inner diameter of 20 mm and 8 mm under atmosphere condition. The bubbly flow, slug flow, annular flow and mist flow are observed for two size tubes. Most of the experimental points agree well with Hewitt and Roberts flow pattern map when putting them on the map. With the experimental results, the range line between the flow patterns is suggested for the tubes of 20 mm inner diameter as well as 8 mm. For the vertical flow up, the air superficial velocity when the annular flow appears is independent of the water superficial velocity. The two phase flow patterns range line is similar in tendency and scope for different tubes. The range line between mist flow and annular flow for different tubes is almost overlapping while that between bubbly flow and slug flow is inconsistent for different tubes. The confusion error between bubbly flow and slug flow is of no great importance for heat transfer because the heat transfer for bubbly flow or for slug flow is same. It is believed that the effect of the tube diameter on the two-phase flow in vertical tubes could be neglected, and the Hewitt and Roberts flow pattern map could be adopted without any consideration of tube diameter.%以空气和水作为介质,在管径分别为20mm与8mm的垂直上升管内进行了常压下气液两相流流型的实验研究,得到了这两种管径下泡状流、弹状流、乳状流和环状流等流型,得到的绝大部分实验点与Hewitt和Roberts流型图相符合,并根据实验结果修正了流型图的转化边界.对于气液两相流垂直上升流动,环状流发生所需的气相折算速度几乎不随液相折算速度的变化而变化.不同管径条件下,各种气液两相流流型发生的范围和转换趋势基本一致,乳状流向环状流的转换界限基本重叠,而泡状流与弹状流的界限变化大一些.由于弹状流的换热与泡状流的换热
Sato, S. (Tsuyama National College of Technology, Okayama (Japan)); Furukawa, A. (Kyushu University, Fukuoka (Japan). Faculty of Engineering); Takamatsu, Y. (Ariake National College of Technology, Fukuoka (Japan))
1993-11-25
An air/water two-phase flow experiment was carried out on the impellers of a centrifugal pump to study the lifting performance and the flow aspect of gas phase. In the experiment, pressurized air is fed to the blowing pipe of a vertical-shaft type pump through a compressor. Transparent acryl resin was used to form the side wall, etc. of the casing through which a video picture of the flowing aspect of the gas phase was taken. The results showed the flowing aspect suddenly changes due to the increase of air flow rate in a low flow rate region where the angle of incidence of flow is large, and the negative pressure sides of impellers were covered with gas. At this time, the lift lowers sharply and discontinuously and then gradually with extension of the gas residence region. This effect appeared more clearly as the outlet angle of the impeller increased. The experimental result agrees roughly with that of the air-bubble calculation in the region where the gas-liquid ratio is so low that a fine air-bubble flow is maintained, but not in the region where the lift lowers sharply, approaching the result of separate flow calculation. The lift after the gas residence region occurs decreases gradually with the increase of air flow rate, showing the same tendency as the result of separate-flow calculation. 13 refs., 14 figs., 1 tab.
Francisco García
2009-01-01
Full Text Available Se realizó un estudio experimental de patrones de flujo bifásico aire-agua en tuberías horizontales y ligeramente inclinadas. Se realizaron 493 experimentos de flujo bifásico aire-agua de los cuales 191 corresponden a tubería horizontal y 302 corresponden a flujo ascendente. Las distribuciones espaciales de los experimentos incluyen los patrones de flujo estratificado liso y ondulado, tapón, anular y burbuja dispersa. Se desarrollan mapas de patrones de flujo experimentales para cada ángulo de inclinación y se evalúa la capacidad de predicción de cuatro modelos mecanicistas y dos modelos de correlación utilizados comúnmente en la literatura para determinar patrones de flujo. Comparando con los valores experimentales se encuentra que los modelos seleccionados tienen un porcentaje de acierto mayor al 75 %.An experimental study in horizontal and slightly inclined pipelines was performed. A total of 493 air-water two-phase flow experiments were carried out, of which 191 correspond to horizontal pipelines and 302 correspond to upward flow. The space distributions of the experiments include the smooth and wavy stratified, slug, annular and dispersed bubble flow patterns. Experimental maps of the flow patterns for each inclination angle are developed and the prediction capability of four mechanistic models and two correlating models commonly used in the literature for determining flow patterns are evaluated. Comparison between calculated and experimental values indicates that the selected models have a success percentage greater than 75%.
无
2001-01-01
In this companion paper, flow patterns in the upstream and downstream tubes of a sudden-expansion cross-section (SECS) in a vertical straight pipe were presented. The effect of SECS on flow patterns upstream and downstream was analyzed by comparing with flow patterns in uniform cross-section vertical tubes. It is found the effect is great. There exist great instabilities of two-phase flow in the neighboring areas of the SECS both downstream and upstream.
Two-Phase Flow in Heterogeneous Media
Ghaffari, Hamed O
2009-01-01
In this study, we investigate the appeared complexity of two-phase flow (air-water) in a heterogeneous soil where the supposed porous media is non-deformable media which is under the time-dependent gas pressure. After obtaining of governing equations and considering the capillary pressure-saturation and permeability functions, the evolution of the models unknown parameters were obtained. In this way, using COMSOL (FEMLAB) and fluid flow-script Module, the role of heterogeneity in intrinsic permeability was analysed. Also, the evolution of relative permeability of wetting and non-wetting fluid, capillary pressure and other parameters were elicited.
New concept of analytical method for two-phase flow
The authors are developing a new analytical method for vertical upward two-phase flow based on a concept that two-phase flow with minimum pressure energy consumption rate is the most stable and easily flowable two-phase flow for the given boundary conditions and, thus, such two-phase flow should be realized actually. Although this concept is applied basically one-dimensionally in the analytical method, gravity convection effect due to density difference between liquid film on the channel wall and two-phase flow core in the central region of the channel is taken into account through a two-dimensional turbulent flow analysis. An air-water two-phase flow experiment was performed to verify the proposed analytical method. In the present paper, results of the experimental analysis with the proposed method are reported. (author)
Interfacial Area and Interfacial Transfer in Two-Phase Flow Systems (Volume IV. Chapters 15-19)
Guo, T.; Park, J.; Kojasoy, G.
2003-03-15
Experiments were performed on horizontal air-water bubbly two-phase flow, axial flow, stratified wavy flow, and annular flow. Theoretical studies were also undertaken on interfacial parameters for a horizontal two-phase flow.
Interfacial Area and Interfacial Transfer in Two-Phase Flow Systems (Volume II. Chapters 6-10)
Guo, T.; Park, J.; Kojasoy, G.
2003-03-15
Experiments were performed on horizontal air-water bubbly two-phase flow, axial flow, stratified wavy flow, and annular flow. Theoretical studies were also undertaken on interfacial parameters for a horizontal two-phase flow.
Interfacial Area and Interfacial Transfer in Two-Phase Flow Systems (Volume III. Chapters 11-14)
Guo, T.; Park, J.; Kojasoy, G.
2003-03-15
Experiments were performed on horizontal air-water bubbly two-phase flow, axial flow, stratified wavy flow, and annular flow. Theoretical studies were also undertaken on interfacial parameters for a horizontal two-phase flow.
Interfacial Area and Interfacial Transfer in Two-Phase Flow Systems (Volume I. Chapters 1-5)
Guo, T.; Park, J.; Kojasoy, G.
2003-03-15
Experiments were performed on horizontal air-water bubbly two-phase flow, axial flow, stratified wavy flow, and annular flow. Theoretical studies were also undertaken on interfacial parameters for a horizontal two-phase flow.
Wallis, Graham B.
1989-01-01
Some features of two recent approaches of two-phase potential flow are presented. The first approach is based on a set of progressive examples that can be analyzed using common techniques, such as conservation laws, and taken together appear to lead in the direction of a general theory. The second approach is based on variational methods, a classical approach to conservative mechanical systems that has a respectable history of application to single phase flows. This latter approach, exemplified by several recent papers by Geurst, appears generally to be consistent with the former approach, at least in those cases for which it is possible to obtain comparable results. Each approach has a justifiable theoretical base and is self-consistent. Moreover, both approaches appear to give the right prediction for several well-defined situations.
Development of methods for mass flow measurement of non-stationary two-phase flows
The topics of the present work are: 1) functioning and testing of a True Mass Flow Meter (TMFM), 2) development of an infra-red absorption measuring process to determine the single components and the mass flows in a non-stationary air/water two-phase flow, 3) presentation of a radionuclide measuring method to measure two-phase mass flow, 4) description of a test stand to calibrate various two-phase mass flow measuring methods. (RW/LH)
Flow pattern maps in two phase flow: present panorama
In this work is presented a general panorama on the condition that watch over the related understanding to the pattern maps of flux regimes in the two-phase flow. The revision that has been done no exhaustive treat of flux patterns observed in vertical and horizontal ducts. As resulting of this investigation, it has been to make evident the necessity of lighting up with precision the use of flux pattern maps that they are not framed respect to really two-phase flow, but that they correspond really to the simultaneous flux of a gas and a liquid un miscible flowing in adiabatic conditions. The case more common of late these is the relative to the air-water mixture. The observed necessity has generated in the Thermo fluids Department of National Institute of Nuclear Research the restlessness of realizing experimental studies in this area. This in spite of being motive of research over 40 years and also of counting with a vast reported bibliography, on one the hand it has not conveyed to obtain representations of general character. And on the other hand it has origined a great confusion about the applicability of available information. In the same way it is described the advances developed in the experimental studies in the field of forced convection, as to only phase as one in two phases. (Author)
Flow instabilities in two-phase flow system with and without phase change
The gas-liquid two-phase flow of various types, such as single component or multiple components, and boiling two-phase flow or insulated two-phase flow, exist in piping systems, and the undesirable phenomena for the operation of systems such as the large scale pulsation of flow rate and the uneven distribution of flow may occur according to the condition. Generally these phenomena are called unstable flow. The author has carried out the research on unstable flow with air-water two-phase flow system, but a question arose to what extent the results in air-water system are applicable to boiling system. The unstable flow is explained with some examples. In this study, the similarity of pulsation in boiling system and insulated system was clarified, using the examples of pressure drop oscillation and flow rate distribution, and the theory to treat them in unified way was presented. The range of discussion is limited to the phenomena that do not depend on the microstructure of flow. The experimental setups were Freon boiling system, air-water capillary system and air-water vertical tube system. The characteristics of pressure drop oscillation and the fundamental mechanism, the theoretical analysis of pressure drop oscillation, the uneven distribution of flow rate in parallel tubes, the stability of flow rate distribution, and the numerical simulation are reported. (Kako, I.)
Håland, Gaute
2010-01-01
The use of multiphase flow meters (MPFMs) increase each year. The main reasons for this are the cost benefits connected to the use of such a meter by eliminating exclusive lines to the test separator and the test separator itself and the demand for a meter that can measure the produced fluids in real time without separating the phases. The MPFMs measures the flow rate by use of different measurement principles and techniques, and in combination of these, e.g. by means of Venturi and void frac...
Two-phase flow in refrigeration systems
Gu, Junjie; Gan, Zhongxue
2013-01-01
Two-Phase Flow in Refrigeration Systems presents recent developments from the authors' extensive research programs on two-phase flow in refrigeration systems. This book covers advanced mass and heat transfer and vapor compression refrigeration systems and shows how the performance of an automotive air-conditioning system is affected through results obtained experimentally and theoretically, specifically with consideration of two-phase flow and oil concentration. The book is ideal for university postgraduate students as a textbook, researchers and professors as an academic reference book, and b
Coherent Calculation for Air-Water Flow and Boiling Flow by Using CUPID Code
The Korea Atomic Energy Research Institute has been developing a three-dimensional thermal-hydraulic code, called CUPID, which was motivated from practical needs for the realistic simulation of two-phase flows in nuclear reactor components. This paper presents coherent simulation of an air-water flow test and a sub-cooled boiling flow test, and the model implementation of related to them. The closure relations for the air-water flow and sub-cooled boiling flow are turbulence model, interfacial non-drag force, interfacial condensation, wall evaporation model, interfacial area transport equation, and so on
Entrainment 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 of the experiments, measurements of entrainment fraction, droplet entrainment rate and droplet deposition rate have been performed by using a liquid film extraction method. A simple, explicit and non-dimensional correlation developed by Sawant et al. (2008a) 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 (1987). (author)
A research on the mechanisms of transition from annular flow in two-phase pipeline flow
Various kinds mechanisms of transitions from two-phase annular flow in tubes were studied and modelled, and the affection factors on the transitions were also discussed. Some mathematical equations and transition criteria for every mechanisms presented were derived, and an unified general criterion for the annular flow transitions in whole range of pipe inclinations was recommended. The boundaries predicted show good agreement with the air-water two-phase experimental data
Laser diagnostics in two phase flows
Krueger, S.
2001-06-01
The existence of a huge lack of experimental data from both technical and fundamental two phase flows was mentioned. The development of laser based non-intrusive measurement techniques to overcome this problem were the task of this work. An optical flow algorithm was adapted for the determination of the velocity fields of continuous and dispersed phase in flow systems. It was used as data reduction method for the newly developed gaseous imaging velocimetry (GIV) technique. The measurement technique including the data reduction has been validated by comparing it to the well-established particle image velocimetry (PIV). Its applicability on scalar data from 2D two-phase flows and reacting gaseous flows was demonstrated. Laser based measurement techniques concerning 3D two-phase flows have also been developed. Solutions for the measurement of the velocity field of the gaseous phase in between the droplets as well as of the liquid phase in an automotive DI spray have been given. (orig.)
Two-phase flow in fractured rock
This report gives the results of a three-day workshop on two-phase flow in fractured rock. The workshop focused on two-phase flow processes that are important in geologic disposal of nuclear waste as experienced in a variety of repository settings. The goals and objectives of the workshop were threefold: exchange information; describe the current state of understanding; and identify research needs. The participants were divided into four subgroups. Each group was asked to address a series of two-phase flow processes. The following groups were defined to address these processes: basic flow processes; fracture/matrix interactions; complex flow processes; and coupled processes. For each process, the groups were asked to address these four issues: (1) describe the two-phase flow processes that are important with respect to repository performance; (2) describe how this process relates to the specific driving programmatic issues given above for nuclear waste storage; (3) evaluate the state of understanding for these processes; and (4) suggest additional research to address poorly understood processes relevant to repository performance. The reports from each of the four working groups are given here
Two-phase flow measurement by pulsed neutron activation techniques
The Pulsed Neutron Activation (PNA) technique for measuring the mass flow velocity and the average density of two-phase mixtures is described. PNA equipment can be easily installed at different loops, and PNA techniques are non-intrusive and independent of flow regimes. These features of the PNA technique make it suitable for in-situ measurement of two-phase flows, and for calibration of more conventional two-phase flow measurement devices. Analytic relations governing the various PNA methods are derived. The equipment and procedures used in the first air-water flow measurement by PNA techniques are discussed, and recommendations are made for improvement of future tests. In the present test, the mass flow velocity was determined with an accuracy of 2%, and average densities were measured down to 0.08 g/cm3 with an accuracy of 0.04 g/cm3. Both the accuracy of the mass flow velocity measurement and the lower limit of the density measurement are functions of the injected activity and of the total number of counts. By using a stronger neutron source and a larger number of detectors, the measurable density can be decreased by a factor of 12 to .007 g/cm3 for 12.5 cm pipes, and to even lower ranges for larger pipes
Study on flooding in two-phase flow
In a countercurrent two-phase flow, where gas phase flows in the upward direction against a gravity-driven liquid downflow, the liquid downflow rate begins to be limited when the gas flow rate exceeds a certain threshold value. This phenomenon, termed 'flooding', may occur during a loss-of-coolant accident (LOCA) at such locations in reactor coolant system as steam generator (SG) U-tubes in a pressurized water reactor (PWR). Flooding generally tends to reduce the amount of water available for core cooling in emergency situations. Flooding has been studied for various flow conditions and geometries, in particular for vertical channels. Most of these studies were concerned with those situations where the lower entry of the channel is exposed to the gas phase or a gas-continuous two-phase flow, and scarcely dealt with such situations where the liquid is the continuous phase at the channel lower entry. However, in a PWR small-break LOCA, where the reactor coolant inventory is depleted only slowly, the latter situations would be encountered more frequently than the former. The present study is concerned with flooding in a vertical channel whose lower entry is facing to a liquid-continuous two-phase flow. Experiments were conducted using Freon R-113 as a simulant of high-pressure steam-water two-phase flow. Experimental results indicate that flooding for this situation initiates when the two-phase mixture swell level in the channel, which indicates large fluctuations with time, reaches the channel top entry at the peaks of level fluctuations. It was also found that the flooding correlation developed formerly by the authors for air-water flows can be applied to the present R-113 case if the difference in fluid properties are considered appropriately. (author)
Characteristics of low-mass-velocity vertical gas-liquid two-phase flow
Low-mass-velocity two-phase flow in a vertical pipe shows lower void fraction than high-mass-velocity two-phase flow even though their qualities are the same. In order to clarify the flow characteristics of the low-mass-velocity two-phase flow, air-water two-phase flow experiments were conducted under the froth or annular flow conditions. Experimental results show that wall shear stress is positive even though both gas and liquid superficial velocities are positive. Measured water film average velocity is negative under this condition. These results indicate that local flow reversal should exist along the channel wall. This local flow reversal gives to the low void fraction in low-mass-velocity two-phase flow. It is also clarified that the drift flux model can be applied to the low-mass-velocity two-phase flow with local reversal. (author)
Development of gas-liquid two-phase flow interfacial structure in a confined bubbly flow
In gas-liquid two-phase flow systems, the interfacial structure specifies the geometric capability of the interfacial transfer of mass, momentum, and energy between the two phases. In view of this, extensive experiments have been carried out in an air-water upward two-phase flow through a test section of 20-cm in width and 1-cm in gap. In it, the local two-phase flow parameters were acquired by the double-sensor conductivity probe at three different elevations in a wide range of the bubbly flow conditions. The acquired local parameters include void fraction ( α), interfacial area concentration (ai), bubble velocity (Ub), bubble Sauter mean diameter (Dsm) and bubble frequency. By taking advantage of the transparent two-dimensional flow path, the flow regime map was constructed through flow visualization. Examination of the measured parameters reveals the development of the interfacial structure due to bubble interactions. (author)
Two-phase flow and heat transfer under low gravity
Frost, W.
1981-11-01
Spacelab experiment to investigate two-phase flow patterns under gravity uses a water-air mixture experiment. Air and water are circulated through the system. The quality or the mixture or air-water is controlled. Photographs of the test section are made and at the same time pressure drop across the test section is measured. The data establishes a flow regime map under reduced gravity conditions with corresponding pressure drop correlations. The test section is also equipped with an electrical resistance heater in order to allow a flow boiling experiment to be carried out using Freon II. High-speed photographs of the test section are used to determine flow patterns. The temperature gradient and pressure drop along the duct can be measured. Thus, quality change can be measured, and heat transfer calculated.
Research on one-dimensional two-phase flow
In the Part I, the author describes about the fundamental form of the hydraulic basic equations for a one-dimensional two-phase flow (two fluid model). Most of the discussions are concentrated on the treatment of phase change inertial force terms in the equations of motion and the author's equations of motion have a strong uniqueness on the following three points in comparison with conventional equations of motion. (1) To express force balance of unit mass two-phase fluid instead of that of unit volume two-phase fluid. (2) To pick up the unit existing mass and the unit flowing mass as the unit mass of two-phase fluid. (3) To apply the kinetic energy principle instead of the momentum low in the evaluation of stational inertia force term. In these three, the item (1) is for excluding a part of momentum change or kinetic energy change due to mass change of the elementary part of fluid, which is independent of force. The item (2) is not to introduce a phenomenological physical model into the evaluation of phase change inertial force term. And the item (3) is for correctly applying the momentum law taking into account the difference of representative velocities between the main flow fluid (vapor phase or liquid phase) and the phase change part of fluid. In the Part II, characteristics of various kinds of high speed two-phase flow are clarified theoretically by using the basic equations derived in the Part I. It is demonstrated that the steam-water two-phase critical flow with violent flashing and the air-water two-phase critical flow without phase change can be described with fundamentally the same basic equations. Furthermore, by comparing the experimental data from the two-phase critical discharge test and the author's theoretical prediction, the two-phase discharge coefficient, CD, for large sharp-edged orifice is determined as the value which is not affected by the experimental facility characteristics, etc.. (J.P.N.)
An introduction to two-phase flows
This course aims at proposing the necessary background for a rational approach to two-phase flows which are notably present in numerous industrial devices and equipment designed to perform energy transfer or mass transfer. The first part proposes a phenomenological approach to main two-phase flow structures and presents their governing variables. The second part presents some proven measurement techniques. The third part focuses on modelling. It recalls the equation elaboration techniques which are based on basic principles of mechanics and thermodynamics and on the application of different averaging operators to these principles. Some useful models are then presented such as models of pressure loss in a duct. The last chapter addresses some fundamental elements of heat transfers in ebullition and condensation
Modeling of two-phase slug flow
When gas and liquid flow in a pipe, over a range of flow rates, a flow pattern results in which sequences of long bubbles, almost filling the pipe cross section, are successively followed by liquid slugs that may contain small bubbles. This flow pattern, usually called slug flow, is encountered in numerous practical situations, such as in the production of hydrocarbons in wells and their transportation in pipelines; the production of steam and water in geothermal power plants; the boiling and condensation in liquid-vapor systems of thermal power plants; emergency core cooling of nuclear reactors; heat and mass transfer between gas and liquid in chemical reactors. This paper provides a review of two phase slug flow modeling
Stability of oscillatory two phase Couette flow
Coward, Adrian V.; Papageorgiou, Demetrios T.
1993-01-01
We investigate the stability of two phase Couette flow of different liquids bounded between plane parallel plates. One of the plates has a time dependent velocity in its own plane, which is composed of a constant steady part and a time harmonic component. In the absence of time harmonic modulations, the flow can be unstable to an interfacial instability if the viscosities are different and the more viscous fluid occupies the thinner of the two layers. Using Floquet theory, we show analytically in the limit of long waves, that time periodic modulations in the basic flow can have a significant influence on flow stability. In particular, flows which are otherwise unstable for extensive ranges of viscosity ratios, can be stabilized completely by the inclusion of background modulations, a finding that can have useful consequences in many practical applications.
Coupling two-phase fluid flow with two-phase darcy flow in anisotropic porous media
Chen, J.
2014-06-03
This paper reports a numerical study of coupling two-phase fluid flow in a free fluid region with two-phase Darcy flow in a homogeneous and anisotropic porous medium region. The model consists of coupled Cahn-Hilliard and Navier-Stokes equations in the free fluid region and the two-phase Darcy law in the anisotropic porous medium region. A Robin-Robin domain decomposition method is used for the coupled Navier-Stokes and Darcy system with the generalized Beavers-Joseph-Saffman condition on the interface between the free flow and the porous media regions. Obtained results have shown the anisotropic properties effect on the velocity and pressure of the two-phase flow. 2014 Jie Chen et al.
Two-Phase Flow Complexity in Heterogeneous Media
Ghaffari, Hamed O
2009-01-01
In this study, we investigate the appeared complexity of two-phase flow (air/water) in a heterogeneous soil where the supposed porous media is non-deformable media which is under the timedependent gas pressure. After obtaining of governing equations and considering the capillary pressuresaturation and permeability functions, the evolution of the model unknown parameters were obtained. In this way, using COMSOL (FEMLAB) and fluid flow/script Module, the role of heterogeneity in intrinsic permeability was analysed. Also, the evolution of relative permeability of wetting and non-wetting fluid, capillary pressure and other parameters were elicited. In the last part, a complex network approach to analysis of emerged patterns will be employed.
Refrigeration. Two-Phase Flow. Flow Regimes and Pressure Drop
Knudsen, Hans-Jørgen Høgaard
2002-01-01
The note gives the basic definitions used in two-phase flow. Flow regimes and flow regimes map are introduced. The different contributions to the pressure drop are stated together with an imperical correlation from the litterature....
Two Phase Flow Simulation Using Cellular Automata
The classical mathematical treatment of two-phase flows is based on the average of the conservation equations for each phase.In this work, a complementary approach to the modeling of these systems based on statistical population balances of aut omata sets is presented.Automata are entities defined by mathematical states that change following iterative rules representing interactions with the neighborhood.A model of automata for two-phase flow simulation is presented.This model consists of fie lds of virtual spheres that change their volumes and move around a certain environment.The model is more general than the classical cellular automata in two respects: the grid of cellular automata is dismissed in favor of a trajectory generator, and the rules of interaction involve parameters representing the actual physical interactions between phases.Automata simulation was used to study unsolved two-phase flow problems involving high heat flux rates. One system described in this work consists of a vertical channel with saturated water at normal pressure heated from the lower surface.The heater causes water to boil and starts the bubble production.We used cellular automata to describe two-phase flows and the interaction with the heater.General rule s for such cellular automata representing bubbles moving in stagnant liquid were used, with special attention to correct modeling of different mechanisms of heat transfer.The results of the model were compared to previous experiments and correlations finding good agreement.One of the most important findings is the confirmation of Kutateladze's idea about a close relation between the start of critical heat flux and a change in the flow's topology.This was analyzed using a control volume located in the upper surface of the heater.A strong decrease in the interfacial surface just before the CHF start was encountered.The automata describe quite well some characteristic parameters such as the shape of the local void fraction in the
MHD Generators Operating with Two-Phase Liquid Metal Flows
with a sequence of tests in which the generator performance will be studied at mixture, qualities up to 50%. The trends of the data accumulated to date have been verified by analysis. A second generator in which the two-phase mixture is passed directly through the generator with no. deliberate attempt to form a film is also being studied. Experimental studies of air-water mixtures have shown that the conductivity of two-phase mixtures does not deteriorate excessively until mixture qualities of 30% are reached. A continuous functional relationship between the quality and conductivity exists even though the flow pattern traverses the gamut of two-phase flow regimes, from bubble to dispersed annular. From these tests, it appears that the electrical conduction is also through an annular film which exists in the various flow regimes. Such a generator has been operated successfully at very low qualities and flow rates, < 5%; currently its performance characteristics in the higher quality regime are being investigated. (author)
Two-phase flow simulation of aeration on stepped spillway
CHENG Xiangju; LUO Lin; ZHAO Wenqian; LI Ran
2004-01-01
Stepped spillways have existed as escape works for a very long time. It is found that water can trap a lot of air when passing through steps and then increasing oxygen content in water body, so stepped spillways can be used as a measure of re-aeration and to improve water quality of water body. However, there is no reliable theoretical method on quantitative calculation of re-aeration ability for the stepped spillways. By introducing an air-water two-phase flow model, this paper used k-ε turbulence model to calculate the characteristic variables of free-surface aeration on stepped spillway. The calculated results fit with the experimental results well. It supports that the numerical modeling method is reasonable and offers firm foundation on calculating re-aeration ability of stepped spillways. The simulation approach can provide a possible optimization tool for designing stepped spillways of more efficient aeration capability.
A software technique for flow-rate measurement in horizontal two-phase flow
Darwich, T.D.; Toral, H.; Archer, J.S. (Imperial Company (EG))
1991-08-01
This paper presents a software technique for measuring individual phase flow rates in two-phase flow. The technique is based on the extraction, classification, and identification of stochastic features from turbulent pressure and void-fraction waveforms. Experiments in a horizontal air/water loop showed that a set of stochastic features is uniquely related to the individual phase flow rates. The software flowmeter is calibrated in situ by compilation of feature sets related to individual phase flow rates in a data base. On-line flow-rate measurement is made by a pattern recognition technique that identifies the best match to the measured feature vector from the calibration data base.
Air-water upward flow in prismatic channel of rectangular base
Experiments had carried out to investigate the two-phase upward air-water flow structure, in a rectangular test section, by using independent measuring techniques, which comprise direct viewing and photography, electrical probes and gamma-ray attenuation. Flow pattern maps and correlations for flow pattern transitions, void fraction profiles, liquid film thickness and superficial average void fraction are proposed and compared to available data. (Author)
Tracer Partitioning in Two-Phase Flow
Sathaye, K.; Hesse, M. A.
2012-12-01
The concentration distributions of geochemical tracers in a subsurface reservoir can be used as an indication of the reservoir flow paths and constituent fluid origin. In this case, we are motivated by the origin of marked geochemical gradients in the Bravo Dome natural CO2 reservoir in northeastern New Mexico. This reservoir contains 99% CO2 with various trace noble gas components and overlies the formation brine in a sloping aquifer. It is thought that magmatic CO2 entered the reservoir, and displaced the brine. This displacement created gradients in the concentrations of the noble gases. Two models to explain noble gas partitioning in two-phase flow are presented here. The first model assumes that the noble gases act as tracers and uses a first order non-linear partial differential equation to compute the volume fraction of each phase along the displament path. A one-way coupled partial differential equation determines the tracer concentration, which has no effect on the overall flow or phase saturations. The second model treats each noble gas as a regular component resulting in a three-component, two-phase system. As the noble gas injection concentration goes to zero, we see the three-component system behave like the one-way coupled system of the first model. Both the analytical and numerical solutions are presented for these models. For the process of a gas displacing a liquid, we see that a noble gas tracer with greater preference for the gas phase, such as Helium, will move more quickly along the flowpath than a heavier tracer that will more easily enter the liquid phase, such as Argon. When we include partial miscibility of both the major and trace components, these differences in speed are shown in a bank of the tracer at the saturation front. In the three component model, the noble gas bank has finite width and concentration. In the limit where the noble gas is treated as a tracer, the width of the bank is zero and the concentration increases linearly
Two-Phase Flow Pressure Drop in Superhydrophobic Channels
Stevens, Kimberly; Crockett, Julie; Maynes, Daniel R.; Iverson, Brian C.
2015-01-01
Superhydrophobic surfaces promote dropwise condensation, which increases the rate of thermal transport, making them desirable for use in condensers. Adiabatic two-phase flow loops have been constructed to gain insight into the hydrodynamics of two-phase systems, laying the groundwork for further study of condensing flow on superhydrophobic surfaces. A two-phase flow loop to measure pressure drop and visualize the flow patterns of two-phase flow in superhydrophobic channels relative to classic...
Numerical calculation of two-phase flows
The theoretical study of time-varying two-phase flow problems in several space dimensions introduces such a complicated set of coupled nonlinear partial differential equations that numerical solution procedures for high-speed computers are required in almost all but the simplest examples. Efficient attainment of realistic solutions for practical problems requires a finite- difference formulation that is simultaneously implicit in the treatment of mass convection, equations of state, and the momentum coupling between phases. Such a method is described, the equations on which it is based are discussed, and its properties are illustrated by means of examples. In particular, the capability for calculating physical instabilities and other time-varying dynamics, at the same time avoiding numerical instability is emphasized. The computer code is applicable to problems in reactor safety analysis, the dynamics of fluidized dust beds, raindrops or aerosol transport, and a variety of similar circumstances, including the effects of phase transitions and the release of latent heat or chemical energy. (U.S.)
Advances in two-phase flow instrumentation
Multiphase flow measurements have become increasingly in a number of process and power systems. However, the need to predict system behavior under transient and accident conditions in nuclear reactors has given impetus to research in this area. Since moving internal interfaces make theoretical predictions difficult, much information for design and supporting analyses is based on experimental observation. The simplest models involving parameters representing mixture density and mixture mass flux, assume thermal equillibrium of the two phases, and are applicable only to a limited number of situations. Most of the parameters, such as interface area and local mixture density, needed for more sophisticated models, are particularly difficult to measure. At present, there are no truly direct methods for measuring local void fraction or mass flux. Local measurements can be taken for a cross-section using, for example, a system of simultaneously quick-closing valves. These valves obtained for the cross section can be integrated, and the result compared with direct measurements for an entire pipeline. Consistent results tend to support the response-model used
Hasan, Abbas; Lucas, Gary
2007-01-01
In two phase flow, differential pressures technique can be used to measure the volume fraction of the gas phase. In the case where no restriction is available in the pipeline, the differential pressure technique can be used only in vertical or inclined pipelines. Two phase air-water pressure drop across a Venturi meter may change its sign from positive to negative due to change in the compressibility of the gas phase. In other words, the inlet of the venturi (upstream section) is not...
Two-phase flow measurements with advanced instrumented spool pieces and local conductivity probes
A series of two-phase, air-water and steam-water tests performed with instrumented spool pieces and with conductivity probes obtained from Atomic Energy of Canada, Ltd. is described. The behavior of the three-beam densitometer, turbine meter, and drag flowmeter is discussed in terms of two-phase models. Application of some two-phase mass flow models to the recorded spool piece data is made and preliminary results are shown. Velocity and void fraction information derived from the conductivity probes is presented and compared to velocities and void fractions obtained using the spool piece instrumentation
Measurements of local two-phase flow parameters in a boiling flow channel
Local two-phase flow parameters were measured lo investigate the internal flow structures of steam-water boiling flow in an annulus channel. Two kinds of measuring methods for local two-phase flow parameters were investigated. These are a two-conductivity probe for local vapor parameters and a Pitot cube for local liquid parameters. Using these probes, the local distribution of phasic velocities, interfacial area concentration (IAC) and void fraction is measured. In this study, the maximum local void fraction in subcooled boiling condition is observed around the heating rod and the local void fraction is smoothly decreased from the surface of a heating rod to the channel center without any wall void peaking, which was observed in air-water experiments. The distributions of local IAC and bubble frequency coincide with those of local void fraction for a given area-averaged void fraction. (author)
Experimental study on flow pattern transitions for inclined two-phase flow
In this paper, experimental data on flow pattern transition of inclination angles from 0-90 are presented. A test section is constructed 2 mm long and I.D 1inch using transparent material. The test section is supported by aluminum frame that can be placed with any arbitrary inclined angles. The air-water two-phase flow is observed at room temperature and atmospheric condition using both high speed camera and void impedance meter. The signal is sampled with sampling rate 1kHz and is analyzed under fully-developed condition. Based on experimental data, flow pattern maps are made for various inclination angles. As increasing the inclination angels from 0 to 90, the flow pattern transitions on the plane jg-jf are changed, such as stratified flow to plug flow or slug flow or plug flow to bubbly flow. The transition lines between pattern regimes are moved or sometimes disappeared due to its inclined angle
Air-water countercurrent annular flow
Bharathan, D.
1979-09-01
Countercurrent annular flow of air and water in circular tubes of diameters ranging from 6.4 to 152 mm is investigated. Experimental measurements include liquid fraction, pressure gradients and countercurrent gas and liquid fluxes. Influences of tube end geometries on the countercurrent fluxes are isolated. Analogies between countercurrent flow, open channel flow, and compressible flow are established. Interfacial momentum transfer between the phases are characterized by empirical friction factors. The dependence of interfacial friction factors on tube diameter is shown to yield a basis for extending the present results to larger tubes.
Air-water countercurrent annular flow
Countercurrent annular flow of air and water in circular tubes of diameters ranging from 6.4 to 152 mm is investigated. Experimental measurements include liquid fraction, pressure gradients and countercurrent gas and liquid fluxes. Influences of tube end geometries on the countercurrent fluxes are isolated. Analogies between countercurrent flow, open channel flow, and compressible flow are established. Interfacial momentum transfer between the phases are characterized by empirical friction factors. The dependence of interfacial friction factors on tube diameter is shown to yield a basis for extending the present results to larger tubes
Flow visualization on flooding phenomena in two phase flow
Flow pattern on flooding phenomena in counter-current two phase flow in vertical tube is experimentally studied by means of dye tracer technique for the following condition: the length -to-diameter ratio L/D=30. Just before flooding, oscillatory motion of the liquid upwards and downwards in the tube is observed. Under flooding conditions, churn flow in the tube is observed. (author)
Chi, Changqing
2016-01-01
Ferrofluids currently are the only type of magnetic liquid materials with wide practical use. The theory on ferrofluids is an example of success to apply statistics to science. Ferrofluids are two-phase liquids consisting of dispersed nanoscale ferromagnetic particles suspended in a carrier fluid. Due to their tiny size, individual ferromagnetic particles clearly exhibit Brownian motions. Only when a large number of randomly-moving particles are subject to an external magnetic field, can they...
Next steps in two-phase flow: executive summary
DiPippo, R.
1980-09-01
The executive summary includes the following topics of discussion: the state of affairs; the fundamental governing equations; the one-dimensional mixture model; the drift-flux model; the Denver Research Institute two-phase geothermal flow program; two-phase flow pattern transition criteria; a two-fluid model under development; the mixture model as applied to geothermal well flow; DRI downwell instrumentation; two-phase flow instrumentation; the Sperry Research Corporation downhole pump and gravity-head heat exchanger systems; and the Brown University two-phase flow experimental program. (MHR)
Aspects of two-phase gas--liquid flow
A wide range of topics related to current research on liquid-gas flow is reviewed, and the relevance of these topics to the design of heat exchangers is discussed. Information is included on flow patterns; system variables; mathematical models for parallel flow and non-parallel flow; critical two-phase flow; unsteady flow; and types of two-phase flow equipment used in industry. (U.S.)
Two-phase flow dynamics in ECC
The present report summarizes the achievements within the project ''Two-phase Systems and ECC''. The results during 1978 - 1980 are accounted for in brief as they have been documented in earlier reports. The results during the first half of 1981 are accounted for in greater detail. They contain a new model for the Basset force and test runs with this model using the test code RISQUE. Furthermore, test runs have been performed with TRAC-PD2 MOD 1. This code was implemented on Edwards Pipe Blowdown experiment (a standard test case) and UC-Berkeley Reflooding experiment (a non-standard test case.) (Auth.)
Numerical method for two-phase flow discontinuity propagation calculation
In this paper, we present a class of numerical shock-capturing schemes for hyperbolic systems of conservation laws modelling two-phase flow. First, we solve the Riemann problem for a two-phase flow with unequal velocities. Then, we construct two approximate Riemann solvers: an one intermediate-state Riemann solver and a generalized Roe's approximate Riemann solver. We give some numerical results for one-dimensional shock-tube problems and for a standard two-phase flow heat addition problem involving two-phase flow instabilities
Program determines two-phase flow
When a mixture of a gas and a liquid flows along a horizontal pipe, it is possible to have up to seven different flow patterns. These flow patterns are: 1. Dispersed. When nearly all the liquid is entrained as spray by the gas; 2. Annular. The liquid forms a film around the inside wall of the pipe, and the gas flows at a high velocity as a central core; 3. Bubble. Bubbles of gas move along at about the same velocity as the liquid; 4. Stratified. The liquid flows along the bottom of the pipe and the gas flows above over a smooth gas-liquid interface; 5. Wave. Is similar to stratified except the interface is disturbed by waves moving in the direction of flow; 6. Slug. Waves are picked up periodically in the form of frothy slugs that move at a much greater velocity than the average liquid velocity; 7. Plug. Alternate plugs of liquid and gas move along the pipe
Structural developments of turbulent two-phase flow in large pipes
In connection with the thermohydraulic problems of two-phase flow that may be encountered under certain operating conditions in piping systems containing heat sources and sinks such as a CANDU reactor heat transport system, this study investigates some of the turbulent characteristics of both cocurrent air-water two-phase flow and single phase flow in large pipes with horizontal orientation. Pitot tubes together with hot film anemometry have been shown to be an adequate measurement system in turbulent dispersed two-phase flow. A practical semi-empirical formula has been developed to predict local mixture velocity as a function of differential head read by Pitot tube, local void fraction, flow pattern constant, gas-liquid properties, momentum transfer factor and two-phase flow quality. The structural developments of the dispersed mixture velocity was studied along a straight horizontal PVC run and expressed in terms of the radial distance and the pipeline length. A correlation is introduced to determine the local mixture velocity in terms of radial and streamwise distance, two-phase flow quality, gas and liquid densities. A similar correlation is presented to predict the local developments of the void fraction. In addition to those normalized correlations, hypothetical interpretations of the experienced phenonema are presented. It was found that the mixture velocity is significantly influenced by the volumetric mixing ratio of both phases. Conclusions are drawn in the special cases of turbulent single and two-phase flow
Visualization and measurement of two-phase flow in tight rod bundle by neutron tomography
Neutron tomography thermal-hydraulic measurement technique is originally developed based on the neutron radiography, computed tomography and two-phase flow measurement techniques. The purpose of the developing is to measure the void fraction distribution in the Reduced-Moderation Water Reactor which is a water-cooled breeder reactor designed by the JAERI as a future reactor. We have visualized and measured the void distribution of air/water two-phase flow and boiling flow in tight-lattice rod bundles. We used the research reactor JRR-3 as a neutron source. Three-dimensional data can be obtained in order to evaluate the numerical analysis codes. In this manuscript, the neutron tomography system, comparison between the reconstruction methods of computed tomography and examples of the measured two-phase flow data which was taken in the 7 rod bundle with a gap between rods of 1.0 mm. (author)
Acoustics of two-phase pipe flows
van Dijk
2005-01-01
Acoustic signals that are recorded in oil pipelines contain information about the flow. In order to extract this information from the pressure recordings, detailed knowledge about the transmission properties of sound waves in the pipes is required.
INFLUENCE OF SURFACTANT ON TWO-PHASE FLOW REGIME AND PRESSURE DROP IN UPWARD INCLINED PIPES
XIA Guo-dong; CHAI Lei
2012-01-01
The influence of a surfactant on the two-phase flow regime and the pressure drop in upward inclined pipes is investigated for various gas/liquid flow rates.The air/water and air/100 ppm sodium dodecyl sulphate aqueous solution are used as the working fluids.The influence of the surfactant on the two-phase flow regime in upward inclined pipes is investigated using the electrical tomographic technique.For 0°,2.5° and 5° pipe inclinations,the surfactant has obvious effect on the transition from the stratified wavy flow to the annular flow,and the range of the stratified smooth flow regime is also extended to higher gas velocities.For 10°pipe inclination,no stratified flow regime is observed in the air/water flow.In the air/surfactant solution system,however,the stratified flow regime can be found in the range of USG =10m/s-28m/s and USL =0.07 m/s-0.2 m/s.For all inclination angles,the changes of the pressure gradient characteristics are accompanied with the flow pattern transitions.Adding surfactant in a two-phase flow would reduce the pressure gradient significantly in the slug flow and annular flow regimes.In the annular flow regime,the pressure gradient gradually becomes free of the influence of the upward inclined angle,and is only dependent on the property of the two-phase flow.
Stochastic modelling of two-phase flows including phase change
Stochastic modelling has already been developed and applied for single-phase flows and incompressible two-phase flows. In this article, we propose an extension of this modelling approach to two-phase flows including phase change (e.g. for steam-water flows). Two aspects are emphasised: a stochastic model accounting for phase transition and a modelling constraint which arises from volume conservation. To illustrate the whole approach, some remarks are eventually proposed for two-fluid models. (authors)
Thermo-Fluid Dynamics of Two-Phase Flow
Ishii, Mamrou
2011-01-01
"Thermo-fluid Dynamics of Two-Phase Flow, Second Edition" is focused on the fundamental physics of two-phase flow. The authors present the detailed theoretical foundation of multi-phase flow thermo-fluid dynamics as they apply to: Nuclear reactor transient and accident analysis; Energy systems; Power generation systems; Chemical reactors and process systems; Space propulsion; Transport processes. This edition features updates on two-phase flow formulation and constitutive equations and CFD simulation codes such as FLUENT and CFX, new coverage of the lift force model, which is of part
Two-phase-flow models and their limitations
An accurate prediction of transient two-phase flow is essential to safety analyses of nuclear reactors under accident conditions. The fluid flow and heat transfer encountered are often extremely complex due to the reactor geometry and occurrence of transient two-phase flow. Recently considerable progresses in understanding and predicting these phenomena have been made by a combination of rigorous model development, advanced computational techniques, and a number of small and large scale supporting experiments. In view of their essential importance, the foundation of various two-phase-flow models and their limitations are discussed in this paper
Sensitivity study of poisson corruption in tomographic measurements for air-water flows
Munshi, P. (Fraunhofer Institute for Nondestructive Testing, Saarbrucken (Germany)); Vaidya, M.S. (Indian Institute of Technology, Kanpur (India))
1993-01-01
An application of computerized tomography (CT) for measuring void fraction profiles in two-phase air-water flows was reported earlier. Those attempts involved some special radial methods for tomographic reconstruction and the popular convolution backprojection (CBP) method. The CBP method is capable of reconstructing void profiles for nonsymmetric flows also. In this paper, we investigate the effect of corrupted CT data for gamma-ray sources and aCBP algorithm. The corruption in such a case is due to the statistical (Poisson) nature of the source.
Sensitivity study of poisson corruption in tomographic measurements for air-water flows
An application of computerized tomography (CT) for measuring void fraction profiles in two-phase air-water flows was reported earlier. Those attempts involved some special radial methods for tomographic reconstruction and the popular convolution backprojection (CBP) method. The CBP method is capable of reconstructing void profiles for nonsymmetric flows also. In this paper, we investigate the effect of corrupted CT data for gamma-ray sources and aCBP algorithm. The corruption in such a case is due to the statistical (Poisson) nature of the source
Review on two-phase flow instabilities in narrow spaces
Instabilities in two-phase flow have been studied since the 1950s. These phenomena may appear in power generation and heat transfer systems where two-phase flow is involved. Because of thermal management in small size systems, micro-fluidics plays an important role. Typical processes must be considered when the channel hydraulic diameter becomes very small. In this paper, a brief review of two-phase flow instabilities encountered in channels having hydraulic diameters greater than 10 mm are presented. The main instability types are discussed according to the existing experimental results and models. The second part of the paper examines two-phase flow instabilities in narrow spaces. Pool and flow boiling cases are considered. Experiments as well as theoretical models existing in the literature are examined. It was found that several experimental works evidenced these instabilities meanwhile only limited theoretical developments exist in the literature. In the last part of the paper an interpretation of the two-phase flow instabilities linked to narrow spaces are presented. This approach is based on characteristic time scales of the two-phase flow and bubble growth in the capillaries
Gravity Independence of Microchannel Two-Phase Flow Project
National Aeronautics and Space Administration — Most of the amassed two-phase flow and heat transfer knowledge comes from experiments conducted in Earth’s gravity. Space missions span varying gravity...
Research on one-dimensional two-phase flow
In Part I the fundamental form of the hydrodynamic basic equations for a one-dimensional two-phase flow (two-fluid model) is described. Discussions are concentrated on the treatment of phase change inertial force terms in the equations of motion and the author's equations of motion which have a remarkable uniqueness on the following three points. (1) To express force balance of unit mass two-phase fluid instead of that of unit volume two-phase fluid. (2) To pick up the unit existing mass and the unit flowing mass as the unit mass of two-phase fluid. (3) To apply the kinetic energy principle instead of the momentum low in the evaluation of steady inertial force term. In these three, the item (1) is for excluding a part of momentum change or kinetic energy change due to mass change of the examined part of fluid, which is independent of force. The item (2) is not to introduce a phenomenological physical model into the evaluation of phase change inertial force term. And the item (3) is for correctly applying the momentum law taking into account the difference of representative velocities between the main flow fluid (vapor phase or liquid phase) and the phase change part of fluid. In Part II, characteristics of various kinds of high speed two-phase flow are clarified theoretically by the basic equations derived. It is demonstrated that the steam-water two-phase critical flow with violent flashing and the airwater two-phase critical flow without phase change can be described with fundamentally the same basic equations. Furthermore, by comparing the experimental data from the two-phase critical discharge test and the theoretical prediction, the two-phase discharge coefficient, CD, for large sharp-edged orifice is determined as the value which is not affected by the experimental facility characteristics, etc. (author)
Adiabatic boiling of two-phase coolant in upward flow
A mathematical model of the process of adiabatic boiling (self-condensation) of a two-phase coolant in upward (downward) flow is developed. The model takes account of changes in phase properties with static pressure decrease. The process is investigated numerically. Approximate analytical formulas for design calculations are obtained. It is shown that effects of adiabatic boiling (self-condensation) should be taken into account when calculating two-phase coolant flow in stretched vertical channels
Mechanistic multidimensional analysis of horizontal two-phase flows
The purpose of this paper is to discuss the results of analysis of two-phase flow in horizontal tubes. Two flow situations have been considered: gas/liquid flow in a long straight pipe, and similar flow conditions in a pipe with 90 deg. elbow. The theoretical approach utilizes a multifield modeling concept. A complete three-dimensional two-phase flow model has been implemented in a state-of-the-art computational multiphase fluid dynamics (CMFD) computer code, NPHASE. The overall model has been tested parametrically. Also, the results of NPHASE simulations have been compared against experimental data for a pipe with 90 deg. elbow.
Two-Phase flow instrumentation for nuclear accidents simulation
Monni, G.; De Salve, M.; Panella, B.
2014-11-01
The paper presents the research work performed at the Energy Department of the Politecnico di Torino, concerning the development of two-phase flow instrumentation and of models, based on the analysis of experimental data, that are able to interpret the measurement signals. The study has been performed with particular reference to the design of power plants, such as nuclear water reactors, where the two-phase flow thermal fluid dynamics must be accurately modeled and predicted. In two-phase flow typically a set of different measurement instruments (Spool Piece - SP) must be installed in order to evaluate the mass flow rate of the phases in a large range of flow conditions (flow patterns, pressures and temperatures); moreover, an interpretative model of the SP need to be developed and experimentally verified. The investigated meters are: Turbine, Venturi, Impedance Probes, Concave sensors, Wire mesh sensor, Electrical Capacitance Probe. Different instrument combinations have been tested, and the performance of each one has been analyzed.
Numerical simulation of two phase flows in heat exchangers
The report presents globally the works done by the author in the thermohydraulic applied to nuclear reactors flows. It presents the studies done to the numerical simulation of the two phase flows in the steam generators and a finite element method to compute these flows. (author)
Two phase discharge flow prediction in safety valves
Safety relief valves (SRV) are necessary elements in the protection of any pressurised system and the prediction of the expected discharge flows is an important consideration for the valve sizing to ensure that rupture pressures do not occur. The high speed flows that occur inside the SRV are complex particularly when a two-phase flow is involved and lead to a less capable protection device which result in larger valves compared to single phase flows. In this paper the ability of a CFD based two phase mixture model to predict the critical flows of air and water through a safety valve is examined. An industrial refrigeration safety relief valve of ¼″ inlet bore size has been tested experimentally over a pressure range of 6–15 barg and air mass qualities from 0.1 to 1 when discharging to near atmospheric conditions for a fully open condition. A two-dimensional mixture model consisting of mixture mass, momentum, and energy equations, combined with a liquid mass equation and the standard k–ε turbulence model for mixture turbulent transport has been used to predict the two phase flows through the valve. The mixture model results have been compared with the Homogenous Equilibrium Model (HEM) commonly used for in valve sizing in non flashing two phase flow conditions. The accuracy of the models over the two phase flow range are quantified and discussed
Definition of two-phase flow behaviors for spacecraft design
Reinarts, Thomas R.; Best, Frederick R.; Miller, Katherine M.; Hill, Wayne S.
1991-01-01
Data for complete models of two-phase flow in microgravity are taken from in-flight experiments and applied to an adiabatic flow-regime analysis to study the feasibility of two-phase systems for spacecraft. The data are taken from five in-flight experiments by Hill et al. (1990) in which a two-phase pump circulates a freon mixture and vapor and liquid flow streams are measured. Adiabatic flow regimes are analyzed based on the experimental superficial velocities of liquid and vapor, and comparisons are made with the results of two-phase flow regimes at 1 g. A motion analyzer records the flow characteristics at a rate of 1000 frames/sec, and stratified flow regimes are reported at 1 g. The flow regimes observed under microgravitational conditions are primarily annular and include slug and bubbly-slug regimes. The present data are of interest to the design and analysis of two-phase thermal-management systems for use in space missions.
Multidimensional two-phase flow regime distribution in a PWR downcomer during an LBLOCA refill phase
The multidimensional countercurrent two-phase flow regimes that occur in a pressurized-water reactor (PWR) vessel downcomer during the refill phase of a large-break loss-of-coolant accident are studied using a transparent 1/10 scale model of a PWR vessel. The various flow regimes and their distribution in the downcomer have been identified and mapped for a range of air-water flooding experiments. The two-phase flow patterns that are identified in the downcomer included various types of film flows, droplet flows, countercurrent churn flows and cocurrent flows depending on the flooding condition. Through observation of the two-phase flow dynamics it was deduced that the physical mechanisms associated with the flooding processes could be separated into a liquid entrainment process and a film flow reversal process. In addition to the above exercise, the effect of non-uniform injection of water into the downcomer via different combinations of cold leg was studied similarly by determining flooding curves and flow pattern maps. It was found that differences in the flooding characteristic were noticeable for various water inlet configurations when compared with the uniform injection case. The differences could be explained qualitatively in terms of the flooding mechanisms identified previously by examining the flow patterns in the downcomer for the non-uniform injection tests. ((orig.))
Two-phase flow characteristics analysis code: MINCS
Two-phase flow characteristics analysis code: MINCS (Modularized and INtegrated Code System) has been developed to provide a computational tool for analyzing two-phase flow phenomena in one-dimensional ducts. In MINCS, nine types of two-phase flow models-from a basic two-fluid nonequilibrium (2V2T) model to a simple homogeneous equilibrium (1V1T) model-can be used under the same numerical solution method. The numerical technique is based on the implicit finite difference method to enhance the numerical stability. The code structure is highly modularized, so that new constitutive relations and correlations can be easily implemented into the code and hence evaluated. A flow pattern can be fixed regardless of flow conditions, and state equations or steam tables can be selected. It is, therefore, easy to calculate physical or numerical benchmark problems. (author)
Study of the Gas-Liquid Two-Phase Flow Measuring Method Based on the V-Cone Flow Meter
Dong, Feng; Hu, Jun
2007-06-01
The paper presents a proper correlation for measuring air-water two-phase flow based on the V-Cone flow meter (β=0.65). The calibrating methods of discharge coefficient and the corrective coefficient θ of V-Cone were introduced. On the basis of theoretical and experimental studies, the Lin Zonghu correlation was studied. The relationship between θ and static pressure (P) or gas-liquid density ratio was obtained. The root-square-errows of gas mass flow and that of liquid mass flow were 2.25% and 2.23% respectively according to the modified Lin Zonghu correlation.
Analysis of water hammer in two-component two-phase flows
The water hammer phenomena caused by a sudden valve closure in air-water two-phase flows must be clarified for the safety analysis of LOCA in reactors and further for the safety of boilers, chemical plants, pipe transport of fluids such as petroleum and natural gas. In the present paper water hammer phenomena caused by a sudden valve closure in two-component two-phase flow are investigated theoretically and experimentally. The phenomena are more complicated that in single phase-flow due to the fact of the presence of compressible component. Basic partial differential equations based on a one-dimensional homogeneous flow model are solved by the method of characteristic. The analysis is extended to include friction in a two-phase mixture depending on the local flow pattern. The profiles of the pressure transients, the propagation velocity of pressure waves and the effect of valve closure on the transient pressure are found. Different two-phase flow pattern and frictional pressure drop correlations were used including Baker, Chesholm and Beggs ampersand Bril correlations
Designing piping systems for two-phase flow
A wide range of industrial systems, such as thermosiphon reboilers and chemical reactors, involve two-phase gas-liquid flow in conduits. Design of these systems requires information about the flow regime, pressure drop, slug velocity and length, and heat transfer coefficient. An understanding of two-phase flow is critical for the reliable and cost-effective design of such systems. The successful design of a pipeline in two-phase flow, for example, is a two-step process. The first step is the determination of the flow regime. If an undesirable flow regime, such as slug flow, is not anticipated and adequately designed for, the resulting flow pattern can upset a tower control system or cause mechanical failures of piping components. The second step is the calculation of flow parameters such as pressure drop and density to size lines and equipment. Since the mechanism of fluid flow (and heat transfer) depends on the flow pattern, separate flow models are required for different flow patterns
Structure of two-phase slug flow in vertical channels
Based on a suggested model of two-phase slug flow in a wide range of varying regime parameters the available literary data and experimental data obtained by the authors on the intensity of void fraction fluctuations, lengths of slugs and liquid plugs are generalized. It is shown that the magnitude of the void fraction is the determining parameter in the formation of the flow regime and structure. The technique for calculating the characteristic frequency of fluctuations of a two-phase flow is suggested
In view of the need to determine void fraction and flow regime of vapor-liquid two-phase flow in the steam generator test model, domestic made optical probe was applied on a small-scale freon two-phase flow test rig. Optical probe signals were collected at a sampling rate up to 500 Hz and converted into digital form. Both the time signal, and the amplitude probability density function and FFT spectrum function calculated thereof were analysed in the time and frequency domains respectively. The threshold characterizing vapor or liquid contact with the probe tip was determined from the air-water two-phase flow pressure drop test results. Then, the boiling freon two-phase flow void fraction was determined by single threshold method, and compared with numerical heat transfer computation. Typical patterns which were revealed by the above-mentioned time signal and the functions were found corresponding to distinct flow regimes, as corroborated by visual observation. The experiment shows that the optical probe was a promising technique for two-phase flow void fraction measurement and flow regime identification (3 refs., 15 figs., 1 tab.)
Two-Phase Slug Flow Experiments with Viscous Liquids
Diaz, Mariana J.C.
2016-01-01
The challenges behind the multiphase transport of oil and gas mixtures are increasing as the oil and gas industry is moving towards production from non-conventional reservoirs and in remote locations. Transport of high viscosity fluids in long multiphase pipelines is a particular challenge. Previous experiments have shown that gas-liquid slug flow is a frequent two-phase flow pattern at high liquid viscosities. The slug flow regime is an unstable flow, which may lead to operati...
Experimental investigation two phase flow in direct methanol fuel cells
Direct methanol fuel cells (DMFC) have received many attentions specifically for portable electronic applications since it utilize methanol which is in liquid form in atmospheric condition and high energy density of the methanol. Thus it eliminates the storage problem of hydrogen. It also eliminates humidification requirement of polymeric membrane which is a problem in PEM fuel cells. Some electronic companies introduced DMFC prototypes for portable electronic applications. Presence of carbon dioxide gases due to electrochemical reactions in anode makes the problem a two phase problem. A two phase flow may occur at cathode specifically at high current densities due to the excess water. Presence of gas phase in anode region and liquid phase in cathode region prevents diffusion of fuel and oxygen to the reaction sites thus reduces the performance of the system. Uncontrolled pressure buildup in anode region increases methanol crossover through membrane and adversely effect the performance. Two phase flow in both anode and cathode region is very effective in the performance of DMYC system and a detailed understanding of two phase flow for high performance DMFC systems. Although there are many theoretical and experimental studies available on the DMFC systems in the literature, only few studies consider problem as a two-phase flow problem. In this study, an experimental set up is developed and species distributions on system are measured with a gas chromatograph. System performance characteristics (V-I curves) is measured depending on the process parameters (temperature, fuel ad oxidant flow rates, methanol concentration etc)
Simulation of Two-Phase Flow in Sloshing Tanks
Luppes, Roel; Veldman, Arthur; Wemmenhove, Rik; Kuzmin, A
2011-01-01
The CFD simulation tool ComFLOW is applied to study the effect of tank motions on two-phase flow phenomena inside a sloshing tank. An improved VOF method is used to assure an accurate description of the fluid displacement. With a novel “gravity-consistent” density averaging method, spurious velociti
Pigging analysis for gas-liquid two phase flow in pipelines
A new method to analyze transient phenomena caused by pigging in gas-liquid two-phase flow is developed. During pigging, a pipeline is divided into three sections by two moving boundaries, namely the pig and the leading edge of the liquid slug in front of the pig. The basic equations are mass, momentum and energy conservation equations. The boundary conditions at the moving boundaries are determined from the mass conservation across the boundaries, etc. A finite difference method is used to solve the equations numerically. The method described above is also capable of analyzing transient two-phase flow caused by pressure and flow rate changes. Thus the over-all analysis of transient two-phase flow in pipelines becomes possible. A series of air-water two-phase flow pigging experiments was conducted using 105.3 mm diameter and 1436.5 m long test pipeline. The agreement between the measured and the calculated results is very good
Full text of publication follows: A new type instrumentation, average bidirectional flow tube, was suggested to apply to the single and two phase flow condition. Its working principle is similar to that of the Pitot tube. The pressure measured at the front of the flow tube is equal to the total pressure, while that measured at the rear tube is slightly less than static pressure of flow field due to the suction effect at the downstream. It gives an amplification effect of measured pressure difference at the flow tube. The proposed instrumentation has the characteristics that it could be applicable to low flow condition and measure bidirectional flow. It was tested in the air-water vertical and horizontal test sections which have 0.08 m inner diameter. The pressure difference across the average bidirectional flow tube, system pressure, average void fraction and injection phasic mass flow rates were measured on the measuring plane. Test was performed primarily in the single phase water and air flow condition to get the amplification factor k of the flow tube. The test was also performed in the air-water two phase flow condition and the covered flow regimes were bubbly, slug, churn turbulent flow in the vertical pipe and stratified flow in the horizontal pipe. In order to calculate the phasic and total mass flow rates from the measured differential pressure, Chexal drift-flux correlation and momentum exchange factor between the two phases were introduced. The test result shows that the suggested instrumentation with the measured void fraction, Chexal drift-flux correlation and Bosio and Malnes' momentum exchange model can predict the phasic mass flow rates within 15% error compared to the true values. A new momentum exchange model was also suggested and it gives up to 5% improvement of the measured mass flow rate compared to combination of Bosio and Malnes' momentum exchange model. (authors)
Multiparticle imaging velocimetry measurements in two-phase flow
The experimental flow visualization tool, Particle Image Velocimetry (PIV), is being extended to determine the velocity fields in two and three-dimensional, two-phase fluid flows. In the past few years, the technique has attracted quite a lot of interest. PIV enables fluid velocities across a region of a flow to be measured at a single instant in time in global domain. This instantaneous velocity profile of a given flow field is determined by digitally recording particle (microspheres or bubbles) images within the flow over multiple successive video frames and then conducting flow pattern identification and analysis of the data. This paper presents instantaneous velocity measurements in various two and three- dimensional, two-phase flow situations. (author)
Advanced Conceptual Models for Unsaturated and Two-Phase Flow in Fractured Rock
Harihar Rajaram; Robert J. Glass; Michael J. Nicholl; Thomas R. Wood
2007-06-24
The Department of Energy Environmental Management Program is faced with two major issues involving two-phase flow in fractured rock; specifically, transport of dissolved contaminants in the Vadose Zone, and the fate of Dense Nonaqueous Phase Liquids (DNAPLs) below the water table. Conceptual models currently used to address these problems do not correctly include the influence of the fractures, thus leading to erroneous predictions. Recent work has shown that it is crucial to understand the topology, or 'structure' of the fluid phases (air/water or water/DNAPL) within the subsurface.
Effect of entrained liquid on turbulent mixing rate between subchannels in annular two-phase flows
Turbulent mixing rates of gas and liquid phases between the subchannels have been measured for various air-water two-phase annular flows in a multiple channel consisting of the two identical circular subchannels. In order to study effect of entrained liquid in the gas core on the turbulent mixing rates, experiments were conducted for two types of liquid injection method, i.e., a small bore nozzle placed in the subchannel center and a porous wall, at a fixed gas injection method. The result showed that the effect of entrained liquid on the turbulent mixing rates of both phases is negligibly small. (author)
Slug flooding in air-water countercurrent vertical flow
This paper is to study slug flooding in the vertical air-water countercurrent flow loop with a porous liquid injector in the upper plenum. More water penetration into the bottom plenum in slug flooding is observed than the annular flooding because the flow regime changes from the slug flow regime or periodic slug/annular flow regime to annular flow regime due to the hysteresis between the onset of flooding and the bridging film. Experiments were made tubes of 0.995 cm, 2.07 cm, and 5.08 cm in diameter. A mechanistic model for the slug flooding with the solitary wave whose height is four time of the mean film thickness is developed to produce relations of the critical liquid flow rate and the mean film thickness. After fitting the critical liquid flow rate with the experimental data as a function of the Bond number, the gas flow rate for the slug flooding is obtained by substituting the critical liquid flow rate to the annular flooding criteria. The present experimental data evaluate the slug flooding condition developed here by substituting the correlations for mean film thickness models in the literature. The best prediction was made by the correlation for the mean film thickness of the present study which is same as Feind's correlation multiplied by 1.35. (author)
Mathematical modeling of disperse two-phase flows
Morel, Christophe
2015-01-01
This book develops the theoretical foundations of disperse two-phase flows, which are characterized by the existence of bubbles, droplets or solid particles finely dispersed in a carrier fluid, which can be a liquid or a gas. Chapters clarify many difficult subjects, including modeling of the interfacial area concentration. Basic knowledge of the subjects treated in this book is essential to practitioners of Computational Fluid Dynamics for two-phase flows in a variety of industrial and environmental settings. The author provides a complete derivation of the basic equations, followed by more advanced subjects like turbulence equations for the two phases (continuous and disperse) and multi-size particulate flow modeling. As well as theoretical material, readers will discover chapters concerned with closure relations and numerical issues. Many physical models are presented, covering key subjects including heat and mass transfers between phases, interfacial forces and fluid particles coalescence and breakup, a...
Measuring the volumetric flow rate of each of the flowing components is required to be monitored in production logging applications. Hence it is necessary to measure the flow rates of gas, oil and water in vertical and inclined oil wells. An increasing level of interest has been shown by the researchers in developing system for the flow rate measurement in multiphase flows. This paper describes the experimental methodology using a miniature, local four-sensor probe for the measurement of dispersed flow parameters in bubbly two-phase flow for spherical bubbles. To establish interdependent among different parameters corresponding to dispersed flow, the available model has been used to experimentally obtain different parameters such as volume fraction, velocity and bubble shape of the dispersed phase in the bubbly air-water flow.
Flow behavior and pressure drop of two-phase flow through C-shaped bend in vertical plane, (1)
Experimental results are presented on the flow behavior, average void fraction and pressure drop in air-water two-phase flow mixture flowing upward through a C-shaped bend curved in vertical plane. The curved test section of transparent acrylic resin tubing was varied in four versions of (a) 90 mm, (b) 132.5 mm and (c) 180 mm radii of curvature with 16 mm inside diameter tube, and (d) 135 radius with 24 mm diameter tube. The combined action of gravity and centrifugal force acting on the two-phase flow is expressed in terms of a modified Froude number representing the balance of radial forces between those acting on the liquid and the gaseous phases of the flow passing through the reclined vertical U-bend. The average void fraction in the curved test section was determined, and empirically correlated to the pressure drop, by means of a series of nondimensional numbers. (author)
Stability of interfacial waves in two-phase flows
The influence of the interfacial pressure and the flow distribution in the one-dimensional two-fluid model on the stability problems of interfacial waves is discussed. With a proper formulation of the interfacial pressure, the following two-phase phenomena can be predicted from the stability and stationary criteria of the interfacial waves: onset of slug flow, stationary hydraulic jump in a stratified flow, flooding in a vertical pipe, and the critical void fraction of a bubbly flow. It can be concluded that the interfacial pressure plays an important role in the interfacial wave propagation of the two-fluid model. The flow distribution parameter may enhance the flow stability range, but only plays a minor role in the two-phase characteristics. (author). 20 refs., 3 tabs., 4 figs
Safety relief valve performance for two-phase flow
The performance of main steam safety relief valve has been evaluated with respect only to the steam. In the present study, two-phase flow and subcooled water blow-out tests with model valves were performed in order to evaluate the valve's characteristics and performance. From the test results, it was made clear that not only for the steam but also for the two-phase flow the measurement data were hardly affected by scaling and also that the reaction force of the fluid to the valve stem was hardly dependent upon the void fraction. Analytical study was performed using the two-phase flow model in the valve. The results of the analysis showed good agreement with the test data. It was shown from the test and analysis results that the reaction force of the two-phase flow and subcooled water to the valve stem was almost as much as that of the steam flow, and the integrity of the safety relief valve could be maintained. (author)
Wire-mesh sensors for two-phase flow investigations
In the annual report 1996 a new wire-mesh sensor for gas-liquid flows was presented. It was used to visualise the cavitation bubble behind a fast acting shut-off valve in a pipeline with a time resolution of over 1000 frames per second for the first time. In the last two years the sensor was applied to an air-water flow in a vertical pipeline (inner diameter D=51.2 mm) to study the flow structure in a wide range of superficial velocities. Besides the void fraction distributions, the high resolution of the sensor allows to calculate bubble size distributions from the primary measuring data. It was possible to study the evolution of the bubble size distribution along the flow path with growing distance from the gas injection (inlet length, L). (orig.)
Velocity and energy relaxation in two-phase flows
Meyapin, Yannick; Gisclon, Marguerite
2009-01-01
In the present study we investigate analytically the process of velocity and energy relaxation in two-phase flows. We begin our exposition by considering the so-called six equations two-phase model [Ishii1975, Rovarch2006]. This model assumes each phase to possess its own velocity and energy variables. Despite recent advances, the six equations model remains computationally expensive for many practical applications. Moreover, its advection operator may be non-hyperbolic which poses additional theoretical difficulties to construct robust numerical schemes |Ghidaglia et al, 2001]. In order to simplify this system, we complete momentum and energy conservation equations by relaxation terms. When relaxation characteristic time tends to zero, velocities and energies are constrained to tend to common values for both phases. As a result, we obtain a simple two-phase model which was recently proposed for simulation of violent aerated flows [Dias et al, 2010]. The preservation of invariant regions and incompressible li...
Takenaka, N; Fujii, T; Mizubata, M; Yoshii, K
1999-01-01
Three-dimensional void fraction distribution of air-water two-phase flow in a 4x4 rod-bundle near a spacer was visualized by fast neutron radiography using a CT method. One-dimensional cross sectional averaged void fraction distribution was also calculated. The behaviors of low void fraction (thick water) two-phase flow in the rod bundle around the spacer were clearly visualized. It was shown that the void fraction distributions were visualized with a quality similar to those by thermal neutron radiography for low void fraction two-phase flow which is difficult to visualize by thermal neutron radiography. It is concluded that the fast neutron radiography is efficiently applicable to two-phase flow studies.
Electrical Capacitance Probe Characterization in Vertical Annular Two-Phase Flow
Grazia Monni
2013-01-01
Full Text Available The paper presents the experimental analysis and the characterization of an electrical capacitance probe (ECP that has been developed at the SIET Italian Company, for the measurement of two-phase flow parameters during the experimental simulation of nuclear accidents, as LOCA. The ECP is used to investigate a vertical air/water flow, characterized by void fraction higher than 95%, with mass flow rates ranging from 0.094 to 0.15 kg/s for air and from 0.002 to 0.021 kg/s for water, corresponding to an annular flow pattern. From the ECP signals, the electrode shape functions (i.e., the signals as a function of electrode distances in single- and two-phase flows are obtained. The dependence of the signal on the void fraction is derived and the liquid film thickness and the phase’s velocity are evaluated by means of rather simple models. The experimental analysis allows one to characterize the ECP, showing the advantages and the drawbacks of this technique for the two-phase flow characterization at high void fraction.
Controlling two-phase flow in microfluidic systems using electrowetting
Gu, Hao
2011-01-01
Electrowetting (EW)-based digital microfluidic systems (DMF) and droplet-based two-phase flow microfluidic systems (TPF) with closed channels are the most widely used microfluidic platforms. In general, these two approaches have been considered independently. However, integrating the two technologie
Dynamic Modeling of Phase Crossings in Two-Phase Flow
Madsen, Søren; Veje, Christian; Willatzen, Morten
2012-01-01
here a numerical implementation and novel study of a fully distributed dynamic one-dimensional model of two-phase flow in a tube, including pressure drop, heat transfer, and variations in tube cross-section. The model is based on a homogeneous formulation of the governing equations, discretized by a...
Determination of bubble parameters in two-phase flow
A development of a probe-detector system for measurement of bubble parameters like size, rise velocity and void fraction in two-phase flow is presented. The method uses an electro resistivity probe and a compact electronic circuit has been developed for obtain this purpose. (author)
Fluid dynamics of cryogenic two-phase flows
The objective of this study was to examine the flow behavior of a methane hydrate/methane-liquid hydrogen dispersed two-phase fluid through a given design of a moderator chamber for the ESS target system. The calculations under simplified conditions, e.g., taking no account of heat input from outside, have shown that the computer code used, CFX, was able to simulate the behavior of the two-phase flow through the moderator chamber, producing reasonable results up to a certain level of the solid phase fraction, that allowed a continuous flow process through the chamber. Inlet flows with larger solid phase fractions than 40 vol% were found to be a ''problem'' for the computer code. From the computer runs based on fractions between 20 and 40 vol%, it was observed that with increasing solid phase fraction at the inlet, the resulting flow pattern revealed a strong tendency for blockage within the chamber, supported by the ''heavy weight'' of the pellets compared to the carrying liquid. Locations which are prone to the development of such uneven flow behavior are the areas around the turning points in the semispheres and near the exit of the moderator. The considered moderator chamber with horizontal inlet and outlet flow for a solid-liquid two-phase fluid does not seem to be an appropriate design. (orig.)
Flow regime transition criteria for two-phase flow in a vertical annulus
Julia, J. Enrique, E-mail: bolivar@emc.uji.es [Departamento de Ingenieria Mecanica y Construccion, Universitat Jaume I., Campus de Riu Sec, 12071 Castellon (Spain); Hibiki, Takashi [School of Nuclear Engineering, Purdue University, 400 Central Dr., West Lafayette, IN 47907-2017 (United States)
2011-10-15
Highlights: > Flow regime transition model is presented for two-phase flows in a vertical annulus. > The transition criteria is easy to be implemented in computational codes. > Final equations do not need experimental input. > New developed model shows better predicting capabilities than existing correlations. > New developed model shows good predicting capabilities in boiling flow. - Abstract: In this work, a new flow regime transition model is proposed for two-phase flows in a vertical annulus. Following previous works, the flow regimes considered are bubbly (B), slug (S) or cap-slug (CS), churn (C) and annular (A). The B to CS transition is modeled using the maximum bubble package criteria of small bubbles. The S to C transition takes place for small annulus perimeter flow channels and it is assumed to occur when the mean void fraction over the entire region exceeds that over the slug-bubble section. If the annulus perimeter is larger that the distorted bubble limit the cap-slug flow regime will be considered since in these conditions it is not possible to distinguish between cap and partial-slug bubbles. The CS to C transition is modeled using the maximum bubble package criteria. However, this transition considers the coalescence of cap and spherical bubbles in order to take into account the flow channel geometry. Finally, the C to A transition is modeled assuming two different mechanisms, (a) flow reversal in the liquid film section along large bubbles; (b) destruction on liquid slugs or large waves by entrainment or deformation. In the S to C and C to A flow regime transitions the annulus flow channel is considered as a rectangular flow channel with no side walls. In all the modeled transitions the drift-flux model is used to obtain the final correlations. The final equations for every flow regime transition are easy to be implemented in computational codes and not experimental input is needed. The prediction accuracy of the newly developed model has been
High speed motion neutron radiography of two-phase flow
Current research in the area of two-phase flow utilizes a wide variety of sensing devices, but some limitations exist on the information which can be obtained. Neutron radiography is a feasible alternative to ''see'' the two-phase flow. A system to perform neutron radiographic analysis of dynamic events which occur on the order of several milliseconds has been developed at Oregon State University. Two different methods have been used to radiograph the simulated two-phase flow. These are pulsed, or ''flash'' radiography, and high speed movie neutron radiography. The pulsed method serves as a ''snap-shot'' with an exposure time ranging from 10 to 20 milliseconds. In high speed movie radiography, a scintillator is used to convert neutrons into light which is enhanced by an optical intensifier and then photographed by a high speed camera. Both types of radiography utilize the pulsing capability of the OSU TRIGA reactor. The principle difficulty with this type of neutron radiography is the fogging of the image due to the large amount of scattering in the water. This difficulty can be overcome by using thin regions for the two-phase flow or using heavy water instead of light water. The results obtained in this paper demonstrate the feasibility of using neutron radiography to obtain data in two-phase flow situations. Both movies and flash radiographs have been obtained of air bubbles in water and boiling from a heater element. The neutron radiographs of the boiling element show both nucleate boiling and film boiling. (Auth.)
Two phase flow instabilities in horizontal straight tube evaporator
Liang, Nan; Shuangquan, Shao; Tian, Changqing; Yan, Y. Y.
2010-01-01
Abstract It is essential to ensure the stability of a refrigeration system if the oscillation in evaporation process is the primary cause for the whole system instability. This paper is concerned with an experimental investigation of two phase flow instabilities in a horizontal straight tube evaporator of a refrigeration system. The relationship between pressure drop and mass flow with constant heat flux and evaporation pressure is measured and determined. It is found that there is...
Simulation of two-phase flow with varying surface tension.
Lervåg, Karl Yngve
2008-01-01
This thesis is a study on the effects of varying surface tension along an interface separating two fluids. Varying surface tension leads to tangential forces along the interface. This is often called the Marangoni effect. These forces are discussed in detail, and two test cases are considered to analyse the Marangoni effect, and to verify the present implementation. The first test studies steady-state two-phase flow where the fluids are separated with plane interfaces and the flow is driv...
Two-phase Flow Distribution in Heat Exchanger Manifolds
Vist, Sivert
2004-01-01
The current study has investigated two-phase refrigerant flow distribution in heat exchange manifolds. Experimental data have been acquired in a heat exchanger test rig specially made for measurement of mass flow rate and gas and liquid distribution in the manifolds of compact heat exchangers. Twelve different manifold designs were used in the experiments, and CO2 and HFC-134a were used as refrigerants.
Two-Phase flow instrumentation for nuclear accidents simulation
The paper presents the research work performed at the Energy Department of the Politecnico di Torino, concerning the development of two-phase flow instrumentation and of models, based on the analysis of experimental data, that are able to interpret the measurement signals. The study has been performed with particular reference to the design of power plants, such as nuclear water reactors, where the two-phase flow thermal fluid dynamics must be accurately modeled and predicted. In two-phase flow typically a set of different measurement instruments (Spool Piece – SP) must be installed in order to evaluate the mass flow rate of the phases in a large range of flow conditions (flow patterns, pressures and temperatures); moreover, an interpretative model of the SP need to be developed and experimentally verified. The investigated meters are: Turbine, Venturi, Impedance Probes, Concave sensors, Wire mesh sensor, Electrical Capacitance Probe. Different instrument combinations have been tested, and the performance of each one has been analyzed
Design and construction of two phases flow meter
This paper deals with design of the gamma ray correlometer and flow loop system for measuring the velocity between two parallel cross-sections of a pipeline. In the laboratory, the radioisotope source and detector were collimated by brass with small beam slit respectively. The flow loop system consists of transparent pipeline, adjustable frequency pump and water container. As a result, when the construction of the flow loop and correlometer is completed, the velocity of two phases flow can be measured by the cross-correlation techniques. (Author)
Hydrodynamics of single- and two-phase flow in inclined rod arrays
Required inputs for thermal-hydraulic codes are constitutive relations for fluid-solid flow resistance, in single-phase flow, and interfacial momentum exchange (relative phase motion), in two-phase flow. An inclined rod array air-water experiment was constructed to study the hydrodynamics of multidimensional porous medium flow in rod arrays. Velocities, pressures, and bubble distributions were measured in square rod arrays of P/d = 1.5, at 0, 30, 45, and 90 degree inclinations to the vertical flow direction. Constitutive models for single-phase flow resistance are reviewed, new comprehensive models developed, and an assessment with previously published and new data made. The principle of superimposing one-dimensional correlations proves successful for turbulent single-phase inclined flow. For bubbly two-phase incline flow a new flow separation phenomena was observed and modeled. A two-region liquid velocity model is developed to explain the experimentally observed phenomena. Fundamental data for bubbles rising in rod arrays were also taken
Two-phase flow instabilities in a vertical annular channel
Babelli, I.; Nair, S.; Ishii, M. [Purdue Univ., West Lafayette, IN (United States)
1995-09-01
An experimental test facility was built to study two-phase flow instabilities in vertical annular channel with emphasis on downward flow under low pressure and low flow conditions. The specific geometry of the test section is similar to the fuel-target sub-channel of the Savannah River Site (SRS) Mark 22 fuel assembly. Critical Heat Flux (CHF) was observed following flow excursion and flow reversal in the test section. Density wave instability was not recorded in this series of experimental runs. The results of this experimental study show that flow excursion is the dominant instability mode under low flow, low pressure, and down flow conditions. The onset of instability data are plotted on the subcooling-Zuber (phase change) numbers stability plane.
Hydrodynamics of single- and two-phase flow in inclined rod arrays
Required inputs for thermal-hydraulic codes are constitutive relations for fluid-solid flow resistance, in single-phase flow, and interfacial momentum exchange (relative phase motion), in two-phase flow. An inclined rod array air-water experiment was constructed to study the hydrodynamics of multidimensional porous medium flow in rod arrays. Velocities, pressures, bubble distributions, and void fractions were measured in inline and rotational square rod arrays of P/d = 1.5, at 0, 30, 45, and 90 degree inclinations to the vertical flow direction. Constitutive models for single-phase flow resistance are reviewed, new comprehensive models developed, and an assessment with previously published and new data made. The principle of superimposing one-dimensional correlations proves successful for turbulent single-phase inclined flow. For bubbly two-phase yawed flow through incline rod arrays a new flow separation phenomena was observed and modeled. Bubbles of diameters significantly smaller than the rod diameter travel along the rod axis, while larger diameter bubbles move through the rod array gaps. The outcome is a flow separation not predictable with current interfacial momentum exchange models. This phenomenon was not observed in rotated square rod arrays. Current interfacial momentum exchange models were confirmed for this rod arrangement. Models for the two phase flow resistance multiplier for cross flow were reviewed and compared with data from cross and yawed flow rod arrays. Both drag and lift components of the multiplier were well predicted by the homogenous model. Other models reviewed overpredicted the data by a factor of two
The additional force arising from rolling motion causes the change of two phase flow resistance. The characteristics of air-water mixture flow resistance in narrow rectangular channel (40 mm × 1.6 mm) in rolling motions were experimentally studied with the rolling periods of 8 s, 12 s, 16 s and the angles of 10°, 15°, 30°, respectively. The experimental results show that the transient frictional factor changes periodically. The larger the mass quality and the maximum rolling angle and the smaller the rolling period, the larger the fluctuation amplitude of frictional pressure drop. A new correction for predicting the transient frictional factor of two-phase flow in rolling motion is proposed with a good accuracy. (authors)
Analysis of water hammer in two-component two-phase flows
The water hammer phenomena caused by a sudden valve closure in air-water two-phase flows must be clarified for the safety analysis of LOCA in reactors and further for the safety of boilers, chemical plants, pipe transport of fluids such as petroleum and natural gas. In the present work water hammer phenomena caused by sudden valve closure in two-component two-phase flows are investigated theoretically and experimentally. The phenomena are more complicated than in single phase-flows due to the fact of the presence of compressible component. Basic partial differential equations based on a one-dimensional homogeneous flow model are solved by the method of characteristic. The analysis is extended to include friction in a two-phase mixture depending on the local flow pattern. The profiles of the pressure transients, the propagation velocity of pressure waves and the effect of valve closure on the transient pressure are found. Different two-phase flow pattern and frictional pressure drop correlations were used including Baker, Chesholm and Beggs and Bril correlations. The effect of the flow pattern on the characteristic of wave propagation is discussed primarily to indicate the effect of void fraction on the velocity of wave propagation and on the attenuation of pressure waves. Transient pressure in the mixture were recorded at different air void fractions, rates of uniform valve closure and liquid flow velocities with the aid of pressure transducers, transient wave form recorders interfaced with an on-line pc computer. The results are compared with computation, and good agreement was obtained within experimental accuracy
Flow regime development analysis in adiabatic upward two-phase flow in a vertical annulus
Julia, J. Enrique [Departamento de Ingenieria Mecanica y Construccion, Universitat Jaume I, Campus de Riu Sec, Castellon 12071 (Spain); Ozar, Basar [School of Nuclear Engineering, Purdue University, 400 Central Dr., West Lafayette, IN 47907-2017 (United States); Jeong, Jae-Jun [Korea Atomic Energy Research Institute, 150 Dukjin, Yuseong, Daejeon 305-353 (Korea, Republic of); Hibiki, Takashi [School of Nuclear Engineering, Purdue University, 400 Central Dr., West Lafayette, IN 47907-2017 (United States); Ishii, Mamoru, E-mail: ishii@purdue.ed [School of Nuclear Engineering, Purdue University, 400 Central Dr., West Lafayette, IN 47907-2017 (United States)
2011-02-15
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{sub 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 <
Dynamic modelling for two-phase flow systems
Several models for two-phase flow have been studied, developing a thermal-hydraulic analysis code with one of these models. The program calculates, for one-dimensional cases with variable flow area, the transient behaviour of system process variables, when the boundary conditions (heat flux, flow rate, enthalpy and pressure) are functions of time. The modular structure of the code, eases the program growth. In fact, the present work is the basis for a general purpose accident and transient analysis code in nuclear reactors. Code verification has been made against RETRAN-02 results. Satisfactory results have been achieved with the present version of the code. (Author)
Two-fluid model for two-phase flow
Ishii, M.
1987-06-01
The two-fluid model formulation is discussed in detail. The emphasis of the paper is on the three-dimensional formulation and the closure issues. The origin of the interfacial and turbulent transfer terms in the averaged formulation is explained and their original mathematical forms are examined. The interfacial transfer of mass, momentum, and energy is proportional to the interfacial area and driving force. This is not a postulate but a result of the careful examination of the mathematical form of the exact interfacial terms. These two effects are considered separately. Since all the interfacial transfer terms involve the interfacial area concentration, the accurate modeling of the local interfacial area concentration is the first step to be taken for a development of a reliable two-fluid model closure relations. The interfacial momentum interaction has been studied in terms of the standard-drag, lift, virtual mass, and Basset forces. Available analytical and semi-empirical correlations and closure relations are reviewed and existing shortcomings are pointed out. The other major area of importance is the modeling of turbulent transfer in two-phase flow. The two-phase flow turbulence problem is coupled with the phase separation problem even in a steady-state fully developed flow. Thus the two-phase turbulence cannot be understood without understanding the interfacial drag and lift forces accurately. There are some indications that the mixing length type model may not be sufficient to describe the three-dimensional turbulent and flow structures. Although it is a very difficult challenge, the two-phase flow turbulence should be investigated both experimentally and analytically with long time-scale research.
Measurement of liquid layer in two-phase flow by neutron radiography
The main objective of applying NRG to heat-transferring flow phenomena is to visualize without external disturbance the phenomena in metal pipes and so on which are difficult to visualize. Besides, from the information on the luminance of images, quantitative physical information can be obtained. In Japan Atomic Energy Research Institute, in fiscal year 1993, the experiment on the measurement of liquid films in gas-liquid two-phase flow was carried out. The steam-water experiment was carried out for the purpose of obtaining the systematic data on the thinning of liquid films by the spacers in the simulated fuel assembly channels. The air-water experiment was carried out for the purpose of obtaining the data on the film thickness of falling water against rising air flow in a rectangular channel. The experimental setup and the results are reported. The achievements so far and the plan for hereafter are described. (K.I.)
Experimental Study on Two-Phase Flow in Horizontal Rectangular Minichannel with Y-Junction
Agus Santoso
2016-03-01
Full Text Available An experimental study was conducted to investigate two-phase air-water flow characteristics, in horizontal rectangular minichannel with Y-junction. The width (W, the height (H and the hydraulic diameter (DH of the rectangular cross section for the upstream side of the junction are 4.60 mm, 2.50 mm and 3.24 mm, while those for the downstream side are 2.36 mm, 2.50 mm and 2.43 mm. The entire test section was machined from transparent acrylic block, so that the flow structure could be visualized. Liquid single-phase and air-liquid twophase flow experiments were conducted at room temperature. The flow pattern, the bubble velocity, the bubble length, and the void fraction were measured with a high-speed video camera. Pressure profile upstream and downstream from the junction was also measured for the respective flows, and the pressure loss due to the contraction at the junction was determined from the pressure profiles. Two flow patterns, i.e., slug and annular flows, were observed in the fully-developed region apart from the junction. In the analysis, the frictional pressure drop data, the two-phase frictional multiplier data, bubble velocity data, bubble length data and void fraction data were compared with calculations by some correlations in literatures. In addition, new pressure loss coefficient correlations for the pressure drop at the junction has been proposed. Results of such experiment and analysis are described in the present paper.
Assessment of intermittent two-phase flow using a high-speed visualization technique
The intermittent two-phase horizontal flow is recognized by many authors to be one of the most highly complex inherently unsteady flow pattern. Great efforts have been made in recent years to develop theoretical models to predict the intermittent flow interfacial parameters. However, as far as the flow structure itself is concerned, these models are not capable to predicting them yet. The liquid slug and gas bubble velocity, length and evolution along the flow are examples of structural parameters that require experimental specialised work to their evaluation. In this paper, we describe the application of a visualization technique to assess the intermittent two-phase air-water horizontal flow. The technique consists in employing a high-speed digital camera to take a series of pictures of the bubbles moving in the upper side of a horizontal pipe, where a gas-liquid mixture flows. From the images obtained the bubbles can be tracked along. Analyzing frame by frame it is possible to extract the gas bubbles lengths and velocities. This investigation is restricted to liquid superficial velocities ranging from 0.3 to 2.0 m/s and to gas superficial velocities ranging from 0.1 to 1.6 m/s. (author)
Reactor vessel and core two-phase flow ultrasonic densitometer
A local ultrasonic density (LUD) detector has been developed by EG and G Idaho, Inc., at the Idaho National Engineering Laboratory (INEL) for the Loss-of-Fluid Test (LOFT) reactor vessel and core two-phase flow density measurements. The principle of operating the sensor is the change in propagation time of a torsional ultrasonic wave in a metal transmission line as a function of the density of the surrounding media. A theoretical physics model is presented which represents the total propagation time as a function of the sensor modulus of elasticity and polar moment of inertia. Separate effects tests and two-phase flow tests have been conducted to characterize the detector. Tests show the detector can perform in a 3430C pressurized water reactor environment and measure the average density of the media surrounding the sensor
Two-phase flow experiments through intergranular stress corrosion cracks
Experimental studies of critical two-phase water flow, through simulated and actual intergranular stress corrosion cracks, were performed to obtain data to evaluate a leak flow rate model and investigate acoustic transducer effectiveness in detecting and sizing leaks. The experimental program included a parametric study of the effects of crack geometry, fluid stagnation pressure and temperature, and crack surface roughness on leak flow rate. In addition, leak detection, location, and leak size estimation capabilities of several different acoustic transducers were evaluated as functions of leak rate and transducer position. This paper presents flow rate data for several different cracks and fluid conditions. It also presents the minimum flows rate detected with the acoustic sensors and a relationship between acoustic signal strength and leak flow rate
Controlling two-phase flow in microfluidic systems using electrowetting
Gu, Hao
2011-01-01
Electrowetting (EW)-based digital microfluidic systems (DMF) and droplet-based two-phase flow microfluidic systems (TPF) with closed channels are the most widely used microfluidic platforms. In general, these two approaches have been considered independently. However, integrating the two technologies into one allows to combine the advantages of both worlds: (i) high throughput (from TPF) and (ii) precise control over each individual drop (from EW). Thus the aim of this thesis was to investiga...
Phase appearance or disappearance in two-phase flows
Cordier, Floraine; Degond, Pierre; Kumbaro, Anela
2011-01-01
This paper is devoted to the treatment of specific numerical problems which appear when phase appearance or disappearance occurs in models of two-phase flows. Such models have crucial importance in many industrial areas such as nuclear power plant safety studies. In this paper, two outstanding problems are identified: first, the loss of hyperbolicity of the system when a phase appears or disappears and second, the lack of positivity of standard shock capturing schemes such as the Roe scheme. ...
Controllability and observability in two-phase porous media flow
Van Doren, J.F.M.; Van den Hof, P.M.J.; Bosgra, O.H.; Jansen, J. D.
2013-01-01
Reservoir simulation models are frequently used to make decisions on well locations, recovery optimization strategies etc. The success of these applications is, among other aspects, determined by the controllability and observability properties of the reservoir model. In this paper it is shown how the controllability and observability of two-phase flow reservoir models can be analyzed and quantified with aid of generalized empirical Gramians. The empirical controllability Gramian can be inter...
A diffuse interface model for two-phase ferrofluid flows
Nochetto, Ricardo H.; Salgado, Abner J.; Tomas, Ignacio
2016-01-01
We develop a model describing the behavior of two-phase ferrofluid flows using phase field-techniques and present an energy-stable numerical scheme for it. For a simplified, yet physically realistic, version of this model and the corresponding numerical scheme we prove, in addition to stability, convergence and as by-product existence of solutions. With a series of numerical experiments we illustrate the potential of these simple models and their ability to capture basic phenomenological feat...
Fluctuation model of a nonequilibrium two-phase channel flow
An ill-posed Cauchy problem for a model of a nonequilibrium two-phase flow in the barotropic approximation is transformed into a well-posed problem by changing the type of the initial hyperbolic equations. Approximation of fluctuations of the phase velocities by a random delta-correlated process and averaging of the equations over its realizations generate a system of parabolic equations. Results of numerical integration of this system are compared with experiment and calculations by well-known models
Fluctuation model of a nonequilibrium two-phase channel flow
Krivoshei, F.A. [Inst. of Engineering Thermophysics, Kiev (Ukraine)
1994-12-01
An ill-posed Cauchy problem for a model of a nonequilibrium two-phase flow in the barotropic approximation is transformed into a well-posed problem by changing the type of the initial hyperbolic equations. Approximation of fluctuations of the phase velocities by a random delta-correlated process and averaging of the equations over its realizations generate a system of parabolic equations. Results of numerical integration of this system are compared with experiment and calculations by well-known models.
Recent advances in two-phase flow numerics
The authors review three topics in the broad field of numerical methods that may be of interest to individuals modeling two-phase flow in nuclear power plants. The first topic is iterative solution of linear equations created during the solution of finite volume equations. The second is numerical tracking of macroscopic liquid interfaces. The final area surveyed is the use of higher spatial difference techniques
Recent advances in two-phase flow numerics
Mahaffy, J.H.; Macian, R. [Pennsylvania State Univ., University Park, PA (United States)
1997-07-01
The authors review three topics in the broad field of numerical methods that may be of interest to individuals modeling two-phase flow in nuclear power plants. The first topic is iterative solution of linear equations created during the solution of finite volume equations. The second is numerical tracking of macroscopic liquid interfaces. The final area surveyed is the use of higher spatial difference techniques.
Two-phase flow induced vibrations in CANDU steam generators
The U-Bend region of nuclear steam generators tube bundles have suffered from two-phase cross flow induced vibrations. Tubes in this region have experienced high amplitude vibrations leading to catastrophic failures. Turbulent buffeting and fluid-elastic instability has been identified as the main causes. Previous investigations have focused on flow regime and two-phase flow damping ratio. However, tube bundles in steam generators have vapour generated on the surface of the tubes, which might affect the flow regime, void fraction distribution, turbulent intensity levels and tube-flow interaction, all of which have the potential to change the tube vibration response. A cantilevered tube bundle made of electric cartridges heaters was built and tested in a Freon-11 flow loop at McMaster University. Tubes were arranged in a parallel triangular configuration. The bundle was exposed to two-phase cross flows consisting of different combinations of void from two sources, void generated upstream of the bundle and void generated at the surface of the tubes. Tube tip vibration response was measured optically and void fraction was measured by gamma densitometry technique. It was found that tube vibration amplitude in the transverse direction was reduced by a factor of eight for void fraction generated at the tube surfaces only, when compared to the upstream only void generation case. The main explanation for this effect is a reduction in the correlation length of the turbulent buffeting forcing function. Theoretical calculations of the tube vibration response due to turbulent buffeting under the same experimental conditions predicted a similar reduction in tube amplitude. The void fraction for the fluid-elastic instability threshold in the presence of tube bundle void fraction generation was higher than that for the upstream void fraction generation case. The first explanation of this difference is the level of turbulent buffeting forces the tube bundle was exposed to
Two-phase boundary layer prediction in upward boiling flow
In the present work, the numerical modelling of the two-phase turbulent boundary layer in upward boiling flow was investigated. First, non-dimensional liquid velocity and temperature profiles in the two-phase boundary layer were validated on the one-dimensional section of a pipe with prescribed radial void fraction profiles. Simulations were performed on a fine grid with a commercial code CFX-5 using the k-ω turbulence model. A significant deviation of results from the analytical single-phase and two-phase wall functions from the literature was observed. Second, a wall boiling model in a vertical heated pipe was simulated (CFX-5) on the coarse grid. Here the prediction of the two-phase thermal boudary layer was compared to the experimental data, k-ω calculation on the fine grid and against the singlephase analytical wall function. Again a major deviation against single-phase temperature wall function was obtained. Presented analyses suggest that the existing analytical velocity and temperature wall functions cannot be valid for the boiling boundary layer with the high void fraction on the wall. (author)
An improved scaling model of buffeting lift forces in air-water flows
This paper presents the results of a series of experiments to study the influence of diameter on the loading of a single rigid cylinder subjected to air-water cross-flow. Five rigid cylinders of same length and different diameters (12.15 * 10-3 m to 31.9 * 10-3 m) were tested over void fractions ranging from 10% to 80%. The fluctuating lift forces on the cylinder are measured and represented in the form of power spectral density. A scaling model of these forces previously developed from one series of experiments with one tube diameter (12.15 * 10-3 m) is tested on these new results by investigating the effect of tube diameter D. Unlike single phase results where the force spectra vary as D(3), it is shown that for two-phase flows, the force spectra vary as D2. The experimental data collapse remarkably well. Both local void fraction and flow regime appear to be sensitive parameters. It confirms the importance of a precise knowledge of the local characteristics of two-phase flows in the study of buffeting forces mechanisms. (authors)
Characteristics of spatiotemporal intermittency in two phase flows
The characteristics of temporal and spatiotemporal intermittent, or slug flow, in two phase flashing flows are analyzed in the context of intermittency transition to chaos. Particularly, the possibility of occurence of generic intermittency routes (Types I, II and III) is investigated in the presence of density wave instability in the system and conclusive evidence for the occurrence of Types I and III is presented. Particularly important is the prediction of the mean slug length and the slug length distribution since the slug flow represents an unfavorable flow regime for gas-liquid transportation in pipes. Identification of two generic intermittency routes automatically gives a quantitative prediction concerning the length distribution of laminar and turbulent (slug) phases. Spatiotemporal analysis, based on the bi-othogonal decoposition and concepts from information theory, provides quantitative characterization and prediction of slug flows and possible mechanism of transition from spatiotemporal intermittency to spatiotemporal chaos (turbulence) is outlined. The analysis of the intermittency in two-phase flows was performed on the experimental data obtained in the study of natural circulation instabilities during small break loss-of-coolant accident. 7 refs., 14 figs
As part of a study on a two-phase natural circulation flow between the outer reactor vessel and the insulation material in the reactor cavity under an external reactor vessel cooling of APR (Advanced Power Reactor) 1400, a K-HERMES-HALF (Hydraulic Evaluation of Reactor cooling Mechanism by External Self-induced flow-HALF scale) experiment was performed at KAERI (Korea Atomic Energy Research Institute) using an air injection method. This experiment was analyzed to verify and evaluate the experimental results using the RELAP5/MOD3 computer code. In addition, the geometry scaling on full height & full sector, and a material scaling between air-water and steam-water two phase natural circulation flow, have been performed for an application of the experimental results to an actual APR1400. The RELAP5/MOD3 results on the water circulation mass flow rate are very similar to the experimental results, in general. The water circulation mass flow rate of the full height & full sector case is approximately 7.6-times higher than that of the K-HERMEL-HALF case. The water circulation mass flow rate of the air injection case is 20-50 % higher than that of the steam injection case at 20 % of the injection rate. (author)
High-frame rate, fast neutron imaging of two-phase flow in a thin rectangular channel
Zboray, R; Dangendorf, V; Stark, M; Tittelmeier, K; Cortesi, M; Adams, R
2015-01-01
We have demonstrated the feasibility of performing high-frame-rate, fast neutron radiography of air-water two-phase flows in a thin channel with rectangular cross section. The experiments have been carried out at the accelerator facility of the Physikalisch-Technische Bundesanstalt. A polychromatic, high-intensity fast neutron beam with average energy of 6 MeV was produced by 11.5 MeV deuterons hitting a thick Be target. Image sequences down to 10 millisecond exposure times were obtained using a fast-neutron imaging detector developed in the context of fast-neutron resonance imaging. Different two-phase flow regimes such as bubbly slug and churn flows have been examined. Two phase flow parameters like the volumetric gas fraction, bubble size and bubble velocities have been measured. The first results are promising, improvements for future experiments are also discussed.
Numerical simulation of two phase flows in heat exchangers
The author gives an overview of his research activity since 1981. He first gives a detailed presentation of properties and equations of two-phase flows in heat exchangers, and of their mathematical and numerical investigation: semi-local equations (mass conservation, momentum conservation and energy conservation), homogenized conservation equations (mass, momentum and enthalpy conservation, boundary conditions), equation closures, discretization, resolution algorithm, computational aspects and applications. Then, he reports the works performed in the field of turbulent flows, hyperbolic methods, low Mach methods, the Neptune project, and parallel computing
Interfacial area measurements in two-phase flow
A thorough understanding of two-phase flow requires the accurate measurement of the time-averaged interfacial area per unit volume (also called the time-averaged integral specific area). The so-called 'specific area' can be estimated by several techniques described in the literature. These different methods are reviewed and the flow conditions which lead to a rigourous determination of the time-averaged integral specific area are clearly established. The probe technique, involving local measurements seems very attractive because of its large range of application
Modelling experimental investigation of two-phase gas-liquid flow patterns in horizontal pipes
In two-phase flow systems the effect of flow patterns on the behavior of the dynamic flow and heat transfer is very important. The study of flow patterns belongs in the fundamental research that have strong practicality. Flow patterns and their transitions in horizontal pipes have been up to now experimentally studied by a number of investigators, and many flow pattern maps and their transition correlations have been proposed. Besides the complexity of flow patterns themselves, the following account for the occurrence of this impracticability: most of these studies were confined to air-water (i.e. two components) systems; most of these studies were confined to normal atmospheric condition; most of these studies were confined to adiabatic condition; there have not common principles on the classification of flow patterns and common definitions of flow patterns yet; and the better coordinate parameters of flow pattern map have not been found as yet. In order to reveal the effect of parameters and components on flow patterns and their transitions thoroughly, in this study the modelling experiment and analysis of horizontal heated channel flow patterns for single component fluid medium (F-12) at higher pressure (the pressures of F-12 were 10-15 bar, corresponding to water system pressures 65-95 bar) were made
Modeling and numerical study of two phase flow
This thesis describes the modelization and the simulation of two-phase systems composed of droplets moving in a gas. The two phases interact with each other and the type of model to consider directly depends on the type of simulations targeted. In the first part, the two phases are considered as fluid and are described using a mixture model with a drift relation (to be able to follow the relative velocity between the two phases and take into account two velocities), the two-phase flows are assumed at the equilibrium in temperature and pressure. This part of the manuscript consists of the derivation of the equations, writing a numerical scheme associated with this set of equations, a study of this scheme and simulations. A mathematical study of this model (hyperbolicity in a simplified framework, linear stability analysis of the system around a steady state) was conducted in a frame where the gas is assumed baro-tropic. The second part is devoted to the modelization of the effect of inelastic collisions on the particles when the time of the simulation is shorter and the droplets can no longer be seen as a fluid. We introduce a model of inelastic collisions for droplets in a spray, leading to a specific Boltzmann kernel. Then, we build caricatures of this kernel of BGK type, in which the behavior of the first moments of the solution of the Boltzmann equation (that is mass, momentum, directional temperatures, variance of the internal energy) are mimicked. The quality of these caricatures is tested numerically at the end. (author)
Development of Interfacial Structure in a Confined Air-Water Cap-Turbulent and Churn-Turbulent Flow
The objective of the present work is to study and model the interfacial structure development of air-water two-phase flow in a confined test section. Experiments of a total of 9 flow conditions in a cap-turbulent and churn-turbulent flow regimes are carried out in a vertical air-water upward two-phase flow experimental loop with a test section of 20-cm in width and 1-cm in gap. The miniaturized four-sensor conductivity probes are used to measure local two-phase parameters at three different elevations for each flow condition. The bubbles captured by the probes are categorized into two groups in view of the two-group interfacial area transport equation, i.e., spherical/distorted bubbles as Group 1 and cap/churn-turbulent bubbles as Group 2. The acquired parameters are time-averaged local void fraction, interfacial velocity, bubble number frequency, interfacial area concentration, and bubble Sauter mean diameter for both groups of bubbles. Also, the line-averaged and area-averaged data are presented and discussed. The comparisons of these parameters at different elevations demonstrate the development of interfacial structure along the flow direction due to bubble interactions
On drag reduction in a two-phase flow
Gatapova, E. Ya.; Ajaev, V. S.; Kabov, O. A.
2015-02-01
Bubbles collected on a local hydrophobic surface with nanocoating in a two-phase flow in a minichannel have been detected experimentally. It has been proposed to use the effect of concentration of gas bubbles on hydrophobic segments of the surface of the channel with contrast wettability for ensuring drag reduction. A two-dimensional flow model with the Navier slip condition in the region of the bubble layer gives criteria of drag reduction, depending on the slip length, dimension of bubbles, and dimension of the segment with nanocoating. The presence of the bubble layer on half of the surface of the channel can increase the flow rate of a liquid flowing through the channel by 40% at a fixed pressure gradient.
Solutal Marangoni instability in layered two-phase flows
Picardo, Jason R; Pushpavanam, S
2015-01-01
In this paper, the instability of layered two-phase flows caused by the presence of a soluble surfactant (or a surface active solute) is studied. The fluids have different viscosities, but are density matched to focus on Marangoni effects. The fluids flow between two flat plates, which are maintained at different solute concentrations. This establishes a constant flux of solute from one fluid to the other in the base state. A linear stability analysis is performed, using a combination of asymptotic and numerical methods. In the creeping flow regime, Marangoni stresses destabilize the flow, provided a concentration gradient is maintained across the fluids. One long wave and two short wave Marangoni instability modes arise, in different regions of parameter space. A well-defined condition for the long wave instability is determined in terms of the viscosity and thickness ratios of the fluids, and the direction of mass transfer. Energy budget calculations show that the Marangoni stresses that drive long and shor...
A real two-phase submarine debris flow and tsunami
Pudasaini, Shiva P.; Miller, Stephen A. [Department of Geodynamics and Geophysics, Steinmann Institute, University of Bonn Nussallee 8, D-53115, Bonn (Germany)
2012-09-26
The general two-phase debris flow model proposed by Pudasaini is employed to study subaerial and submarine debris flows, and the tsunami generated by the debris impact at lakes and oceans. The model, which includes three fundamentally new and dominant physical aspects such as enhanced viscous stress, virtual mass, and generalized drag (in addition to buoyancy), constitutes the most generalized two-phase flow model to date. The advantage of this two-phase debris flow model over classical single-phase, or quasi-two-phase models, is that the initial mass can be divided into several parts by appropriately considering the solid volume fraction. These parts include a dry (landslide or rock slide), a fluid (water or muddy water; e.g., dams, rivers), and a general debris mixture material as needed in real flow simulations. This innovative formulation provides an opportunity, within a single framework, to simultaneously simulate the sliding debris (or landslide), the water lake or ocean, the debris impact at the lake or ocean, the tsunami generation and propagation, the mixing and separation between the solid and fluid phases, and the sediment transport and deposition process in the bathymetric surface. Applications of this model include (a) sediment transport on hill slopes, river streams, hydraulic channels (e.g., hydropower dams and plants); lakes, fjords, coastal lines, and aquatic ecology; and (b) submarine debris impact and the rupture of fiber optic, submarine cables and pipelines along the ocean floor, and damage to offshore drilling platforms. Numerical simulations reveal that the dynamics of debris impact induced tsunamis in mountain lakes or oceans are fundamentally different than the tsunami generated by pure rock avalanches and landslides. The analysis includes the generation, amplification and propagation of super tsunami waves and run-ups along coastlines, debris slide and deposition at the bottom floor, and debris shock waves. It is observed that the
A real two-phase submarine debris flow and tsunami
The general two-phase debris flow model proposed by Pudasaini is employed to study subaerial and submarine debris flows, and the tsunami generated by the debris impact at lakes and oceans. The model, which includes three fundamentally new and dominant physical aspects such as enhanced viscous stress, virtual mass, and generalized drag (in addition to buoyancy), constitutes the most generalized two-phase flow model to date. The advantage of this two-phase debris flow model over classical single-phase, or quasi-two-phase models, is that the initial mass can be divided into several parts by appropriately considering the solid volume fraction. These parts include a dry (landslide or rock slide), a fluid (water or muddy water; e.g., dams, rivers), and a general debris mixture material as needed in real flow simulations. This innovative formulation provides an opportunity, within a single framework, to simultaneously simulate the sliding debris (or landslide), the water lake or ocean, the debris impact at the lake or ocean, the tsunami generation and propagation, the mixing and separation between the solid and fluid phases, and the sediment transport and deposition process in the bathymetric surface. Applications of this model include (a) sediment transport on hill slopes, river streams, hydraulic channels (e.g., hydropower dams and plants); lakes, fjords, coastal lines, and aquatic ecology; and (b) submarine debris impact and the rupture of fiber optic, submarine cables and pipelines along the ocean floor, and damage to offshore drilling platforms. Numerical simulations reveal that the dynamics of debris impact induced tsunamis in mountain lakes or oceans are fundamentally different than the tsunami generated by pure rock avalanches and landslides. The analysis includes the generation, amplification and propagation of super tsunami waves and run-ups along coastlines, debris slide and deposition at the bottom floor, and debris shock waves. It is observed that the
A real two-phase submarine debris flow and tsunami
Pudasaini, Shiva P.; Miller, Stephen A.
2012-09-01
The general two-phase debris flow model proposed by Pudasaini [1] is employed to study subaerial and submarine debris flows, and the tsunami generated by the debris impact at lakes and oceans. The model, which includes three fundamentally new and dominant physical aspects such as enhanced viscous stress, virtual mass, and generalized drag (in addition to buoyancy), constitutes the most generalized two-phase flow model to date. The advantage of this two-phase debris flow model over classical single-phase, or quasi-two-phase models, is that the initial mass can be divided into several parts by appropriately considering the solid volume fraction. These parts include a dry (landslide or rock slide), a fluid (water or muddy water; e.g., dams, rivers), and a general debris mixture material as needed in real flow simulations. This innovative formulation provides an opportunity, within a single framework, to simultaneously simulate the sliding debris (or landslide), the water lake or ocean, the debris impact at the lake or ocean, the tsunami generation and propagation, the mixing and separation between the solid and fluid phases, and the sediment transport and deposition process in the bathymetric surface. Applications of this model include (a) sediment transport on hill slopes, river streams, hydraulic channels (e.g., hydropower dams and plants); lakes, fjords, coastal lines, and aquatic ecology; and (b) submarine debris impact and the rupture of fiber optic, submarine cables and pipelines along the ocean floor, and damage to offshore drilling platforms. Numerical simulations reveal that the dynamics of debris impact induced tsunamis in mountain lakes or oceans are fundamentally different than the tsunami generated by pure rock avalanches and landslides. The analysis includes the generation, amplification and propagation of super tsunami waves and run-ups along coastlines, debris slide and deposition at the bottom floor, and debris shock waves. It is observed that the
Characteristics of gas-liquid two-phase flow in a vertical small diameter tube at a medium pressure
Most of correlations for calculating two-phase flow parameters, such as flow pattern transitions, void fraction and pressure drop, have been developed based on the experimental data on tubes greater than 10 mm in diameter at near atmospheric pressures. For that, the applicability of such correlations is doubtful to the flow in small diameter tubes at a medium pressure as seen in compact heat exchangers like residential room conditioners. In this connection, the purpose of this study is to provide experimental data for gas-liquid two-phase vertical flows in a small diameter tube at medium pressures since the published data for such flows is limited to examine existing correlations and/or develop a new one. Experiments have been conducted on air-water two-phase flows in a vertical circular tube of 9.48 mm internal diameter. In the experiment, system pressure in the channel has been systematically changed from 0.2 to 0.7 MPa (absolute) to study the effect of the pressure on two-phase flow parameters, i.e., two-phase flow pattern transitions, bubble size in bubble flow, void fraction, interfacial shear force, frictional pressure drop and static pressure fluctuations. Furthermore, the respective data obtained have been compared with existing correlations. (author)
Interfacial area transport equation of gas-liquid two-phase flow across a horizontal tube
It is important to improve accuracy of the two-phase flow simulation codes, since these codes are used for the performance estimation and the safety designs of the industrial products which two-phase flow occur inside. In the governing equations of the two-phase flow simulation model, interfacial area concentration is one of the most important parameters to improve the model. Some researcher developed a transport equation of the interfacial area concentration that include source and sink terms due to bubble coalescence and breakup. However, these equations use numerous parameters to obtain the source and sink terms. In this study, we have successfully expressed the source and sink term as functions of a single variable to simplify and enhance their usefulness. The coalescence of bubbles was expressed as a function of the local void fraction, and the breakup of bubbles was expressed as a function of the time averaged local liquid velocity. Using this model, we can simulate the local void fraction and interfacial area concentration. In order to validate the simulation model, we also carried out experiments to measure a time averaged local interfacial area concentration and a time averaged local void fraction in a vertical upward air-water two-phase flow across horizontal tube with a double sensor electrical conductivity probe. The test section comprised a 50mm × 50 mm square channel partially obstructed by a single horizontal tube inside. The superficial gas velocity and superficial liquid velocity ranged from 0.019 to 0.215 m/s and from 0.133 to 0.533 m/s, respectively. Measurements of the void fraction and the interfacial area concentration were performed at three axial locations of z = 110, 340, and 640 mm and transverse locations from r = 0 to 25 mm. The experimental data set were compared to the result of the numerical simulations. Good agreement is obtained between these results. (author)
Experimental and numerical investigation on two-phase flow instabilities
Ruspini, Leonardo Carlos
2013-03-01
Two-phase flow instabilities are experimentally and numerically studied within this thesis. In particular, the phenomena called Ledinegg instability, density wave oscillations and pressure drop oscillations are investigated. The most important investigations regarding the occurrence of two-phase flow instabilities are reviewed. An extensive description of the main contributions in the experimental and analytical research is presented. In addition, a critical discussion and recommendations for future investigations are presented. A numerical framework using a hp-adaptive method is developed in order to solve the conservation equations modelling general thermo-hydraulic systems. A natural convection problem is analysed numerically in order to test the numerical solver. Moreover, the description of an adaptive strategy to solve thermo-hydraulic problems is presented. In the second part of this dissertation, a homogeneous model is used to study Ledinegg, density wave and pressure drop oscillations phenomena numerically. The dynamic characteristics of the Ledinegg (flow excursion) phenomenon are analysed through the simulation of several transient examples. In addition, density wave instabilities in boiling and condensing systems are investigated. The effects of several parameters, such as the fluid inertia and compressibility volumes, on the stability limits of Ledinegg and density wave instabilities are studied, showing a strong influence of these parameters. Moreover, the phenomenon called pressure drop oscillations is numerically investigated. A discussion of the physical representation of several models is presented with reference to the obtained numerical results. Finally, the influence of different parameters on these phenomena is analysed. In the last part, an experimental investigation of these phenomena is presented. The designing methodology used for the construction of the experimental facility is described. Several simulations and a non
Flow regimes of adiabatic gas-liquid two-phase under rolling conditions
Yan, Chaoxing; Yan, Changqi; Sun, Licheng; Xing, Dianchuan; Wang, Yang; Tian, Daogui
2013-07-01
Characteristics of adiabatic air/water two-phase flow regimes under vertical and rolling motion conditions were investigated experimentally. Test sections are two rectangular ducts with the gaps of 1.41 and 10 mm, respectively, and a circular tube with 25 mm diameter. Flow regimes were recorded by a high speed CCD-camera and were identified by examining the video images. The experimental results indicate that the characteristics of flow patterns in 10 mm wide rectangular duct under vertical condition are very similar to those in circular tube, but different from the 1.41 mm wide rectangular duct. Channel size has a significant influence on flow pattern transition, boundary of which in rectangular channels tends asymptotically towards that in the circular tube with increasing the width of narrow side. Flow patterns in rolling channels are similar to each other, nevertheless, the effect of rolling motion on flow pattern transition are significantly various. Due to the remarkable influences of the friction shear stress and surface tension in the narrow gap duct, detailed flow pattern maps of which under vertical and rolling conditions are indistinguishable. While for the circular tube with 25 mm diameter, the transition from bubbly to slug flow occurs at a higher superficial liquid velocity and the churn flow covers more area on the flow regime map as the rolling period decreases.
Objective characterization of interfacial structures in two-phase flow
In view of establishing a detailed and reliable measurement technique for characterizing the interfacial structures and identifying flow regimes in two-phase flow, two objective approaches are presented. First, the state-of-the-art four-sensor conductivity probe technique is presented to obtain the detailed local information. The newly designed four-sensor conductivity probe accommodates the double-sensor probe capability. Hence, it can be applied in a wide range of two-phase flow regimes spanning from bubbly to churn-turbulent flows with a measurement error of approximately ± 10%. The signal processing scheme is developed such that it accounts for the missing bubbles and defective signals. Furthermore, it categorizes the acquired parameters into two groups based on bubble cord length information. Local information on the void fraction, interfacial area concentration, Sauter mean diameter, interface velocity for each group of bubbles was obtained successfully. Second, a global measurement technique using the non-intrusive impedance voidmeter and neural networks is presented. In this method, an advanced non-intrusive impedance voidmeter provides global interfacial structure information to neural networks which are used to identify the flow regimes. Both supervised and self-organizing neural network learning paradigms performed flow regime identification successfully. In the application of this global method, two approaches are presented, namely: One based on the Probability Density Function (PDF input method), and another based on the ordered set of void fraction measurements which were acquired in a very short time period (instantaneous direct signal input method). The direct signal input method minimizes the time required for identifying the flow regime
One-dimensional transient unequal velocity two-phase flow by the method of characteristics
An understanding of two-phase flow is important when one is analyzing the accidental loss of coolant or when analyzing industrial processes. If a pipe in the steam generator of a nuclear reactor breaks, the flow will remain critical (or choked) for almost the entire blowdown. For this reason the knowledge of the two-phase maximum (critical) flow rate is important. A six-equation model--consisting of two continuity equations, two energy equations, a mixture momentum equation, and a constitutive relative velocity equation--is solved numerically by the method of characteristics for one-dimensional, transient, two-phase flow systems. The analysis is also extended to the special case of transient critical flow. The six-equation model is used to study the flow of a nonequilibrium sodium-argon system in a horizontal tube in which the nonequilibrium sodium-argon system in a horizontal tube in which the critical flow condition is at the entrance. A four-equation model is used to study the pressure-pulse propagation rate in an isothermal air-water system, and the results that are found are compared with the experimental data. Proper initial and boundary conditions are obtained for the blowdown problem. The energy and mass exchange relations are evaluated by comparing the model predictions with results of void-fraction and heat-transfer experiments. A simplified two-equation model is obtained for the special case of two incompressible phases. This model is used in the preliminary analysis of batch sedimentation. It is also used to predict the shock formation in the gas-solid fluidized bed
A study on the two-phase flow characteristics of nanofluids
While a considerable body of research exists regarding enhancements of the heat transfer using nanofluids, the basic hydraulic phenomenon of a nanofluid has not been investigated as much. Several studies were reported related to the pressure drop of nanofluid flow and a few researches on the hydraulic characteristics of two-phase nanofluid flow were conducted. Two-phase Flow Analysis in a Helical Wire Inserted Tube using CFD Code An analysis on the two-phase flow in a helical wire inserted tube using commercial CFD code, CFX11.0, was performed in bubbly flow and annular flow regions. The analysis method was validated with the experimental results of Keishi Takeshima. Bubbly and annular flows in a 10 mm inner diameter tube with varying pitch lengths and inserted wire diameters were simulated using the same analysis methods after validation. The geometry range of p/D was 1-4 and e/D was 0.08-0.12. The results show that the inserted wire with a larger diameter increased swirl flow generation. An increasing swirl flow was seen as the pitch length increased. Regarding pressure loss, smaller pitch lengths and inserted wires with larger diameters resulted in larger pressure loss. The average liquid film thickness increased as the pitch length and the diameter of the inserted wire increased in the annular flow region. Both in the bubbly flow and annular flow regions, the effect of pitch length on swirl flow generation and pressure loss was more significant than that of the inserted wire diameters. An Experimental Study on the Two-Phase Flow Characteristics of Nanofluids The main objective of this study is to investigate the basic hydraulic phenomenon of the nanofluid in the two-phase flow region. For the accomplishment of this objective, a series of experiments have been performed. The first one is the pressure drop and pressure fluctuation measurements in a vertically upward air-water two-phase flow. The air and the water based nanofluid were used as working fluids under
Effect of flow direction on two-phase flow distribution of refrigerants at a T-junction
The present study experimentally investigated the effect of flow direction and other flow parameters on two-phase flow distribution of refrigerants at a T-junction, and also suggested a prediction model for refrigerant in a T-junction by modifying previous model for air-water flow. R-22, R-134a, and R-410A were used as test refrigerants. As geometric parameters, the direction of the inlet or branch tube and the tube diameter ratio of branch to inlet tube were chosen. The measured data were compared with the values predicted by the models developed for air-water or steam-water mixture in the literature. We propose a modified model for application to the reduced T-junction and vertical tube orientation. Among the geometric parameters, the branch tube direction showed the biggest sensitivity to the mass flow rate ratio for the gas phase, while the inlet quality showed the biggest sensitivity to the mass flow rate ratio among the inlet flow parameters
Two-phase flow instability in a parallel multichannel system
HOU Suxia
2009-01-01
The two-phase flow instabilities observed in through parallel multichannel can be classified into three types, of which only one is intrinsic to parallel multichannel systems. The intrinsic instabilities observed in parallel multichannel system have been studied experimentally. The stable boundary of the flow in such a parallel-channel system are sought, and the nature of inlet flow oscillation in the unstable region has been examined experimentally under various conditions of inlet velocity, heat flux, liquid temperature, cross section of channel and entrance throttling. The results show that parallel multichannel system possess a characteristic oscillation that is quite independent of the magnitude and duration of the initial disturbance, and the stable boundary is influenced by the characteristic frequency of the system as well as by the exit quality when this is low, and upon raising the exit quality and reducing the characteristic frequency, the system increases its instability, and entrance throttling effectively contributes to stabilization of the system.
Mathematical model of two-phase flow in accelerator channel
О.Ф. Нікулін
2010-01-01
Full Text Available The problem of two-phase flow composed of energy-carrier phase (Newtonian liquid and solid fine-dispersed phase (particles in counter jet mill accelerator channel is considered. The mathematical model bases goes on the supposition that the phases interact with each other like independent substances by means of aerodynamics’ forces in conditions of adiabatic flow. The mathematical model in the form of system of differential equations of order 11 is represented. Derivations of equations by base physical principles for cross-section-averaged quantity are produced. The mathematical model can be used for estimation of any kinematic and thermodynamic flow characteristics for purposely parameters optimization problem solving and transfer functions determination, that take place in counter jet mill accelerator channel design.
Particle-fluid two-phase flow modeling
Mortensen, G. A.; Trapp, J. A.
This paper describes a numerical scheme and computer program, DISCON, for the calculation of two-phase flows that does not require the use of flow regime maps. This model is intermediate between-thermal instantaneous and the averaged two-fluid model. It solves the Eulerian continuity, momentum, and energy equations for each liquid control volume, and the Lagrangian mass, momentum, energy, and position equations for each bubble. The bubbles are modeled individually using a large representative number of bubbles, thus, avoiding the numerical diffusion associated with Eulerian models. DISCON has been used to calculate the bubbling of air through a column of water and the subcooled boiling of water in a flow channel. The results of these calculations are presented.
Two phase flow models in DxUNSp code platform
Catalin NAE
2011-09-01
Full Text Available The aim of this work is to find an efficient implementation for a two phase flow model in an existing URANS CFD code platform (DxUNSp, initially based on unsteady URANS equations with a k- turbulence model and various other extensions, ranging from a broad selection of wall laws up to a very efficient LES model. This code has the capability for development for nonreacting/reacting multifluid flows for research applications and is under continuous progress. It is intend to present mainly three aspects of this implementation for unstructured mesh based solvers, for high Reynolds compressible flows: the importance of the 5/7 equation model, performance with respect to a basic test cases and implementation details of the proposed schemes. From a numerical point of view, we propose a new approximation schemes of this system based on the VFRoe-ncv.
Stability of stratified two-phase flows in horizontal channels
Barmak, Ilya; Ullmann, Amos; Brauner, Neima; Vitoshkin, Helen
2016-01-01
Linear stability of stratified two-phase flows in horizontal channels to arbitrary wavenumber disturbances is studied. The problem is reduced to Orr-Sommerfeld equations for the stream function disturbances, defined in each sublayer and coupled via boundary conditions that account also for possible interface deformation and capillary forces. Applying the Chebyshev collocation method, the equations and interface boundary conditions are reduced to the generalized eigenvalue problems solved by standard means of numerical linear algebra for the entire spectrum of eigenvalues and the associated eigenvectors. Some additional conclusions concerning the instability nature are derived from the most unstable perturbation patterns. The results are summarized in the form of stability maps showing the operational conditions at which a stratified-smooth flow pattern is stable. It is found that for gas-liquid and liquid-liquid systems the stratified flow with smooth interface is stable only in confined zone of relatively lo...
The PDF method for Lagrangian two-phase flow simulations
A recent turbulence model put forward by Pope (1991) in the context of PDF modelling has been used. In this approach, the one-point joint velocity-dissipation pdf equation is solved by simulating the instantaneous behaviour of a large number of Lagrangian fluid particles. Closure of the evolution equations of these Lagrangian particles is based on stochastic models and more specifically on diffusion processes. Such models are of direct use for two-phase flow modelling where the so-called fluid seen by discrete inclusions has to be modelled. Full Lagrangian simulations have been performed for shear-flows. It is emphasized that this approach gives far more information than traditional turbulence closures (such as the K-ε model) and therefore can be very useful for situations involving complex physics. It is also believed that the present model represents the first step towards a complete Lagrangian-Lagrangian model for dispersed two-phase flow problems. (authors). 21 refs., 6 figs
Modeling transient two-phase stratified flow in pipelines
An isothermal, two-fluid model, comprised of separate mass and linear momentum conservation equations for the gas and liquid phases was formulated. Interfacial mass transfer effects were modeled via the black oil method. Both equal and unequal phase pressure formulations were evaluated. The model was used to investigate transient two-phase stratified flow in pipelines. An explicit numerical scheme was used to solve the system of equations. Experimental data were collected in an existing 425 m long, 76.2 mm diameter horizontal pipeline. Good agreement was observed between experimental and predicted results
Strongly coupled dispersed two-phase flows; Ecoulements diphasiques disperses fortement couples
Zun, I.; Lance, M.; Ekiel-Jezewska, M.L.; Petrosyan, A.; Lecoq, N.; Anthore, R.; Bostel, F.; Feuillebois, F.; Nott, P.; Zenit, R.; Hunt, M.L.; Brennen, C.E.; Campbell, C.S.; Tong, P.; Lei, X.; Ackerson, B.J.; Asmolov, E.S.; Abade, G.; da Cunha, F.R.; Lhuillier, D.; Cartellier, A.; Ruzicka, M.C.; Drahos, J.; Thomas, N.H.; Talini, L.; Leblond, J.; Leshansky, A.M.; Lavrenteva, O.M.; Nir, A.; Teshukov, V.; Risso, F.; Ellinsen, K.; Crispel, S.; Dahlkild, A.; Vynnycky, M.; Davila, J.; Matas, J.P.; Guazelli, L.; Morris, J.; Ooms, G.; Poelma, C.; van Wijngaarden, L.; de Vries, A.; Elghobashi, S.; Huilier, D.; Peirano, E.; Minier, J.P.; Gavrilyuk, S.; Saurel, R.; Kashinsky, O.; Randin, V.; Colin, C.; Larue de Tournemine, A.; Roig, V.; Suzanne, C.; Bounhoure, C.; Brunet, Y.; Tanaka, A.T.; Noma, K.; Tsuji, Y.; Pascal-Ribot, S.; Le Gall, F.; Aliseda, A.; Hainaux, F.; Lasheras, J.; Didwania, A.; Costa, A.; Vallerin, W.; Mudde, R.F.; Van Den Akker, H.E.A.; Jaumouillie, P.; Larrarte, F.; Burgisser, A.; Bergantz, G.; Necker, F.; Hartel, C.; Kleiser, L.; Meiburg, E.; Michallet, H.; Mory, M.; Hutter, M.; Markov, A.A.; Dumoulin, F.X.; Suard, S.; Borghi, R.; Hong, M.; Hopfinger, E.; Laforgia, A.; Lawrence, C.J.; Hewitt, G.F.; Osiptsov, A.N.; Tsirkunov, Yu. M.; Volkov, A.N.
2003-07-01
flow, current distribution and mass transfer along a vertical gas evolving electrode; a two-way coupled model for dilute multiphase flows. Topic 3: turbulence modulation by particles, droplets or bubbles in dense systems: influence of particles on the transition to turbulence in pipe flow; comparison between a point particle model and a finite-diameter-model for the particle turbulence interaction in a suspension; the effect on turbulence by bubbles rising through it under buoyancy; the physical mechanisms of modifying the structure of turbulent homogeneous shear flows by dispersed particles; influence of hydrodynamic interactions between particles on the turbulent flow in a suspension; review of relationships between Lagrangian and Eulerian scales; a two-point PDF for modelling turbulent dispersed two-phase flows and derivation of a two field model; mathematical and numerical modeling of two-phase compressible flows with micro-inertia. Topic 4: collective effects in dispersed two-phase flows clustering and phase distribution: hydrodynamic structure of downward bubbly flow; influence of gravity on the dynamics of a turbulent bubbly pipe flow; experimental study of two-phase flows; particle clusters formed in dispersed gas-solid flows: simulations and experiments; experimental study of the turbulence in bubbly flows at high void fraction; first step in the study of the correlation between air/water flow fluctuations and random buffering forces; clustering and settling velocity of micro-droplets in a grid turbulence. Topic 5: large scale instabilities and gravity driven dispersed flows: new 'non-isothermal' linear instability modes in fluidized beds and bubbly flows; large scale instability in a confined buoyant shear layer; convective instability in uniform dispersed layers; structures in gravity driven bubbly flows; effects of concentration profiles on velocity profiles in sewer; pyroclastic density currents viewed as mammoth scale two-phase flows; mixing and
Construction of the two-phase critical flow test facility
The two-phase critical test loop facility has been constructed in the KAERI engineering laboratory for the simulation of small break loss of coolant accident entrained with non-condensible gas of SMART. The test facility can operate at 12 MPa of pressure and 0 to 60 C of sub-cooling with 0.5 kg/s of non- condensible gas injection into break flow, and simulate up to 20 mm of pipe break. Main components of the test facility were arranged such that the pressure vessel containing coolant, a test section simulating break and a suppression tank inter-connected with pipings were installed vertically. As quick opening valve opens, high pressure/temperature coolant flows through the test section forming critical two-phase flow into the suppression tank. The pressure vessel was connected to two high pressure N2 gas tanks through a control valve to control pressure in the pressure vessel. Another N2 gas tank was also connected to the test section for the non-condensible gas injection. The test facility operation was performed on computers supported with PLC systems installed in the control room, and test data such as temperature, break flow rate, pressure drop across test section, gas injection flow rate were all together gathered in the data acquisition system for further data analysis. This test facility was classified as a safety related high pressure gas facility in law. Thus the loop design documentation was reviewed, and inspected during construction of the test loop by the regulatory body. And the regulatory body issued permission for the operation of the test facility
Droplets Formation and Merging in Two-Phase Flow Microfluidics
Hao Gu
2011-04-01
Full Text Available Two-phase flow microfluidics is emerging as a popular technology for a wide range of applications involving high throughput such as encapsulation, chemical synthesis and biochemical assays. Within this platform, the formation and merging of droplets inside an immiscible carrier fluid are two key procedures: (i the emulsification step should lead to a very well controlled drop size (distribution; and (ii the use of droplet as micro-reactors requires a reliable merging. A novel trend within this field is the use of additional active means of control besides the commonly used hydrodynamic manipulation. Electric fields are especially suitable for this, due to quantitative control over the amplitude and time dependence of the signals, and the flexibility in designing micro-electrode geometries. With this, the formation and merging of droplets can be achieved on-demand and with high precision. In this review on two-phase flow microfluidics, particular emphasis is given on these aspects. Also recent innovations in microfabrication technologies used for this purpose will be discussed.
Droplets formation and merging in two-phase flow microfluidics.
Gu, Hao; Duits, Michel H G; Mugele, Frieder
2011-01-01
Two-phase flow microfluidics is emerging as a popular technology for a wide range of applications involving high throughput such as encapsulation, chemical synthesis and biochemical assays. Within this platform, the formation and merging of droplets inside an immiscible carrier fluid are two key procedures: (i) the emulsification step should lead to a very well controlled drop size (distribution); and (ii) the use of droplet as micro-reactors requires a reliable merging. A novel trend within this field is the use of additional active means of control besides the commonly used hydrodynamic manipulation. Electric fields are especially suitable for this, due to quantitative control over the amplitude and time dependence of the signals, and the flexibility in designing micro-electrode geometries. With this, the formation and merging of droplets can be achieved on-demand and with high precision. In this review on two-phase flow microfluidics, particular emphasis is given on these aspects. Also recent innovations in microfabrication technologies used for this purpose will be discussed. PMID:21731459
Computer simulation of two-phase flow in nuclear reactors
Two-phase flow models dominate the requirements of economic resources for the development and use of computer codes which serve to analyze thermohydraulic transients in nuclear power plants. An attempt is made to reduce the effort of analyzing reactor transients by combining purpose-oriented modelling with advanced computing techniques. Six principles are presented on mathematical modeling and the selection of numerical methods, along with suggestions on programming and machine selection, all aimed at reducing the cost of analysis. Computer simulation is contrasted with traditional computer calculation. The advantages of run-time interactive access operation in a simulation environment are demonstrated. It is explained that the drift-flux model is better suited than the two-fluid model for the analysis of two-phase flow in nuclear reactors, because of the latter's closure problems. The advantage of analytical over numerical integration is demonstrated. Modeling and programming techniques are presented which minimize the number of needed arithmetical and logical operations and thereby increase the simulation speed, while decreasing the cost. (orig.)
Study on interfacial area transport of two-phase flow in vertical and large circular pipe
Local distribution of the interfacial area concentration (IAC) and void fraction of air-water two-phase flow in a vertical pipe were investigated by using the measurement of optical fiber probe. The inner diameter of the circular pipe was 100 mm, and superficial gas and liquid velocities ranged from 0 to 0.1 m/s and from 0 to 1.0 m/s. accordingly. It is found that the bubble passing through frequency dominates the radial distribution profile of IAC. The axial interfacial area concentration was calculated based on the interfacial area transport model of Ishii-Kim. In the mean time, the analysis of four bubble interaction mechanisms determining the change of IAC shows that it is the local pressure inside the pipe that dominates the value of IAC in the axial direction. A new correlation including the effect of pressure was given for calculating the axial IAC. (authors)
Two-phase counter-current flow limitations in complex piping systems
Experiments have been performed to investigate two-phase air-water counter-current flow limitations in two complex geometry test sections. One of the test sections was chosen to represent the complexity of CANDU feeders, at a 1:4 scale. In the second configuration, the effect of declining the horizontal parts of the test section was investigated. For the tested geometry, the results showed that the onset of flooding occurs at lower gas velocity as compared to those required to initiate flooding in vertical tubes or single 90 deg vertical bends. However, the critical gas velocity at the zero liquid penetration limit was found to be comparable to that of a 90 deg vertical bend. (author) 18 refs., 78 figs
Dynamic and spectroscopic characteristics of atmospheric gliding arc in gas-liquid two-phase flow
In this study, an atmospheric alternating-current gliding arc device in gas-liquid two-phase flow has been developed for the purpose of waste water degradation. The dynamic behavior of the gas-liquid gliding arc is investigated through the oscillations of electrical signals, while the spatial evolution of the arc column is analyzed by high speed photography. Different arc breakdown regimes are reported, and the restrike mode is identified as the typical fluctuation characteristic of the hybrid gliding arc in air-water mixture. Optical emission spectroscopy is employed to investigate the active species generated in the gas-liquid plasma. The axial evolution of the OH (309 nm) intensity is determined, while the rotational and vibrational temperatures of the OH are obtained by a comparison between the experimental and simulated spectra. The significant discrepancy between the rotational and translational temperatures has also been discussed.
Non-equilibrium effects in transient two-phase flow
Depressurisation tests were carried out on Refrigerant 113 liquid flowing in a horizontal pipeline, under conditions where vapour was formed by flash evaporation. The tests covered a range of initial velocities up to 2.1 m/s in a 51 mm diameter glass pipeline at starting pressure around 1.5 bar and with varying rates of depressurisation. Measurements were made of local pressure and temperature, circulation rate, pressure difference and void fraction variation over a test section length of 2m. The local pressure and temperature measurements give a direct indication of non-equilibrium effects. The vapour formed during the flash evaporation process quickly formed a stratified type flow and a theoretical model was developed on this basis. The model includes the transient two phase low conservation equations allied to a heat transfer equation. Satisfactory agreement between theoretical predictions and experimental results was obtained. (author)
Flooding in counter-current two-phase flow
Flooding is a phenomenon which is best described as the transition from counter-current to co-current flow. Early notice was taken of this phenomenon in the chemical engineering industry. Flooding also plays an important role in the field of two-phase heat transfer since it is a limit for many systems involving counter-current flow. Practical applications of flooding limited processes include wickless thermosyphons and the emergency core cooling system (ECCS) of pressurized water nuclear reactors. The phenomenon of flooding also is involved in the behavior of nuclear reactor core materials during severe accident conditions where flooding is one of the mechanisms governing the motion of the molten fuel pin cladding
Geometry and Material Scaling on Two Phase Natural Circulation Flow for K-HERMES-HALF Experiment
Park, Rae Joon; Ha, Kwang Soon; Kim, Sang Baik; Hong, Seong Wan [Korea Atomic Energy Research Institute, Daejeon (Korea, Republic of); Heo, Sun [KHNP Nuclear Engineering and Technology Institute, Daejeon (Korea, Republic of)
2011-05-15
As part of a study on two-phase natural circulation mass flow rate between the outer reactor vessel and vessel insulation in the reactor cavity under the IVR (In-Vessel corium retention) through the ERVC (External Reactor Vessel Cooling) in APR1400, K-HERMES-HALF experiment (Hydraulic Evaluation of Reactor cooling Mechanism by External Self-induced flow) had performed at KAERI. This large-scale experiment using a half-height and half-sector model of the APR1400 uses the non-heating method of the air injection. For this reason, it is necessary to evaluate the geometry scaling on full height and full sector and a material scaling between air-water and steam-water two phase natural circulation flow for an application of the experimental results to an actual APR1400. In the geometry scaling, two cases, such as a half height and half sector and a full height and full sector, had performed by using RELAP5/MOD3 computer code. In the material scaling, two cases, such as an air injection and a steam injection, had performed to compare the air injection experimental results with the steam injection case
CFD analyses for water / air tests to investigate the RPV exterior two-phase flow behavior
The full text follows. Siemens / KWU develops a new boiling water reactor called SWR 1000. It's safety concept will consist of passive safety equipment combined with active systems, and through this diversity, meets the goal of reducing the probability of core damage compared to existing nuclear plants. Siemens / KWU performs - in co-operation with VTT - tests to quantify the safety margins of the exterior cooling concept for the SWR 1000, which are supposed to be very high, by measuring the critical heat fluxes (CHFs). The following stepwise procedure will be applied for the investigation of the CHFs, whereas the first two steps will be needed to design a model in such a way, that it represents the flow around the RPV: -) Water / air experiments with a 1:10 scaled global model; -) Water / air experiments with a 1:10 section model; -) Water / steam experiments with a 1:1-section model. The CHF will be determined by heating the reproduced reactor pressure vessel wall and measuring wall temperatures. FLUENT4 and CFX -codes have been successfully used for validating the two-phase flow in 1:10 air-water experiments and predicting behavior is 1:10 tests. The results are summarized in the paper. The capability of CFD codes for modeling multidimensional two-phase flow is discussed
Numerical prediction of horizontal two phase flow using an interfacial area density model
During a Lost of Coolant Accident (LOCA) cold Emergency Core Cooling (ECC) water can be injected into the cold leg of the primary loop of PWR. A relevant problem occurring in this situation is the development of wavy stratified flows which can lead to slug generation. Because slug flow cannot be predicted with the required accuracy and spatial resolution by the one-dimensional system codes, the stratified flows are increasingly modelled with computational fluid dynamics (CFD) codes. In CFD, closure models are required that must be validated. The recent improvements of the multiphase flow modelling in the ANSYS CFX code make it now possible to simulate these mechanisms in detail. In order to validate existing and further developed multiphase flow models, high-resolution measurement data is needed in time and also in space. For the experimental investigation of co-current air/water flows, the HAWAC (Horizontal Air/Water Channel) was built. The channel allows in particular the study of air/water slug flow under atmospheric pressure. Parallel to the experiments, CFD calculations were carried out. A picture sequence recorded during slug flow was compared with the equivalent CFD simulation made. The two-fluid model was applied with a special turbulence damping procedure at the free surface. An Algebraic Interfacial Area Density (AIAD) model on the basis of the implemented mixture model was introduced, which 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. The behaviour of slug generation and propagation at the experimental setup was qualitatively reproduced by the simulation, while local deviations require a continuation of the work. The creation of small instabilities due to pressure surge or an increase of
Studies of Two-Phase Flow Dynamics and Heat Transfer at Reduced Gravity Conditions
Witte, Larry C.; Bousman, W. Scott; Fore, Larry B.
1996-01-01
The ability to predict gas-liquid flow patterns is crucial to the design and operation of two-phase flow systems in the microgravity environment. Flow pattern maps have been developed in this study which show the occurrence of flow patterns as a function of gas and liquid superficial velocities as well as tube diameter, liquid viscosity and surface tension. The results have demonstrated that the location of the bubble-slug transition is affected by the tube diameter for air-water systems and by surface tension, suggesting that turbulence-induced bubble fluctuations and coalescence mechanisms play a role in this transition. The location of the slug-annular transition on the flow pattern maps is largely unaffected by tube diameter, liquid viscosity or surface tension in the ranges tested. Void fraction-based transition criteria were developed which separate the flow patterns on the flow pattern maps with reasonable accuracy. Weber number transition criteria also show promise but further work is needed to improve these models. For annular gas-liquid flows of air-water and air- 50 percent glycerine under reduced gravity conditions, the pressure gradient agrees fairly well with a version of the Lockhart-Martinelli correlation but the measured film thickness deviates from published correlations at lower Reynolds numbers. Nusselt numbers, based on a film thickness obtained from standard normal-gravity correlations, follow the relation, Nu = A Re(sup n) Pr(exp l/3), but more experimental data in a reduced gravity environment are needed to increase the confidence in the estimated constants, A and n. In the slug flow regime, experimental pressure gradient does not correlate well with either the Lockhart-Martinelli or a homogeneous formulation, but does correlate nicely with a formulation based on a two-phase Reynolds number. Comparison with ground-based correlations implies that the heat transfer coefficients are lower at reduced gravity than at normal gravity under the same
The Condensation effect on the two-phase flow stability
considering riser condensation and of correcting the localized friction due to the presence of the two-phase mixture in the two-phase region.These effects are more important for high heating power and high inlet subcooling. CAREM 25 nuclear power reactor is investigated to get the stability boundary map. The flow instability regions are appeared at low and high core power. In the low heat flux range, the trends of the thermal equilibrium - equal velocity (homogeneous) model and the thermal non equilibrium - non equal velocity model are the same because the steam quality is small.In the high heat flux range, for the subcooled boiling number and the phase change number, the marginal stability boundaries are crossed in a point, determining tow different regions, of high and low inlet subcooling.For the first region, the steam quality calculation of the first model is greater and has the effect of stabilizing the system more than the second one.For the second region, the two-phase region length calculation of the first model is smaller and has the effect of stabilizing the system less than the second one. In general, the model predicts a more stable system with an increase in inlet restriction or riser condensation or system pressure or a decrease in exit restriction
Tsunami Generated by a Two-Phase Submarine Debris Flow
Pudasaini, S. P.
2012-04-01
The general two-phase debris flow model proposed by Pudasaini (2011) is employed to study subaerial and submarine debris flows, and the tsunami generated by the debris impact at lakes and oceans. The model includes several essential physical aspects, including Mohr-Coulomb plasticity for the solid stress, while the fluid stress is modelled as a solid volume fraction gradient enhanced non-Newtonian viscous stress. The generalized interfacial momentum transfer includes the viscous drag, buoyancy, and the virtual mass. The generalized drag covers both the solid-like and fluid-like contributions, and can be applied to linear to quadratic drags. Strong couplings exist between the solid and the fluid momentum transfer. The advantage of the real two-phase debris flow model over classical single-phase or quasi-two-phase models is that by considering the solid (and/or the fluid) volume fraction appropriately, the initial mass can be divided into several (even mutually disjoint) parts; a dry (landslide or rock slide), a fluid (water or muddy water; e.g., dams, rivers), and a general debris mixture material as needed in real flow simulations. This offers a unique and innovative opportunity within a single framework to simultaneously simulate (a) the sliding debris (or landslide), (b) the water lake or ocean, (c) the debris impact at the lake or ocean, (d) tsunami generation and propagation, (e) mixing and separation between the solid and the fluid phases, and (f) sediment transport and deposition process in the bathymetric surface. The new model is applied to two-phase subaerial and submarine debris flows. Benchmark numerical simulations reveal that the dynamics of the debris impact induced tsunamis are fundamentally different than the tsunami generated by pure rock avalanche and landslides. Special attention is paid to study the basic features of the debris impact to the mountain lakes or oceans. This includes the generation, amplification and propagation of the multiple
Propagation of density disturbances in air-water flow
Nassos, G. P.
1969-01-01
Study investigated the behavior of density waves propagating vertically in an atmospheric pressure air-water system using a technique based on the correlation between density change and electric resistivity. This information is of interest to industries working with heat transfer systems and fluid power and control systems.
Solutal Marangoni instability in layered two-phase flow
Picardo, Jason; Radhakrishna, T. G.; Pushpavanam, S.
2015-11-01
In this work, the instability of layered two-phase flow caused by the presence of a surface-active solute is studied. The fluids are density matched to focus on surfactant effects. The fluids flow between two flat plates, which are maintained at different solute concentrations. This establishes a constant flux of soluble surfactant from one fluid to the other, in the base state. A linear stability analysis is carried out, supported by energy budget calculations. The flow is first analyzed in the creeping flow regime. Long wave as well as short wave Marangoni instabilities are identified, each with a distinct energy signature. The short wave instability manifests as two distinct modes, characterized by the importance of interfacial deformations or lack thereof. The primary instability switches between these different modes as parameters are varied. The effect of small but finite inertia on these solutal Marangoni modes is then examined. The effect of soluble surfactant on a finte inertia flow is also studied, with focus on the transition from the viscosity-induced instability to solutal Marangoni instability. This analysis is relevant to microfluidic applications, such as solvent extraction, in which mass transfer is carried out between stratified immiscible fluids.
Wall function approach for boiling two-phase flows
One of the important goals of the NURESIM project is to assess and improve the simulation capability of the three-dimensional two-fluid codes for prediction of local boiling flow processes. The boiling flow is strongly affected by local mechanisms in the turbulent boundary layer near the heated wall. Wall-to-fluid transfer models for boiling flow with the emphasis on near-wall treatment are being addressed in the paper. Since the computational grid of the 3D two-fluid models is too coarse to resolve the variable gradients in the near-wall region, the use of wall functions is a common approach to model the liquid velocity and temperature profile adjacent to the heated wall. The wall function model for momentum, based on the surface roughness analogy has been discussed and implemented in the NEPTUNECFD code. The model has been validated on several upward boiling flow experiments, differing in the geometry, working fluid and operating conditions. The simulations with the new wall function model show an improved prediction of flow parameters over the boiling boundary layer. Furthermore, a wall function model for the energy equation, based on enhanced two-phase wall friction has been derived and validated.
Stability of stratified two-phase flows in horizontal channels
Barmak, I.; Gelfgat, A.; Vitoshkin, H.; Ullmann, A.; Brauner, N.
2016-04-01
Linear stability of stratified two-phase flows in horizontal channels to arbitrary wavenumber disturbances is studied. The problem is reduced to Orr-Sommerfeld equations for the stream function disturbances, defined in each sublayer and coupled via boundary conditions that account also for possible interface deformation and capillary forces. Applying the Chebyshev collocation method, the equations and interface boundary conditions are reduced to the generalized eigenvalue problems solved by standard means of numerical linear algebra for the entire spectrum of eigenvalues and the associated eigenvectors. Some additional conclusions concerning the instability nature are derived from the most unstable perturbation patterns. The results are summarized in the form of stability maps showing the operational conditions at which a stratified-smooth flow pattern is stable. It is found that for gas-liquid and liquid-liquid systems, the stratified flow with a smooth interface is stable only in confined zone of relatively low flow rates, which is in agreement with experiments, but is not predicted by long-wave analysis. Depending on the flow conditions, the critical perturbations can originate mainly at the interface (so-called "interfacial modes of instability") or in the bulk of one of the phases (i.e., "shear modes"). The present analysis revealed that there is no definite correlation between the type of instability and the perturbation wavelength.
Local measurement of interfacial area, interfacial velocity and liquid turbulence in two-phase flow
Double sensor probe and hotfilm anemometry methods were developed for measuring local flow characteristics in bubbly flow. The formulation for the interfacial area concentration measurement was obtained by improving the formulation derived by Kataoka and Ishii. The assumptions used in the derivation of the equation were verified experimentally. The interfacial area concentration measured by the double sensor probe agreed well with one by the photographic method. The filter to validate the hotfilm anemometry for measuring the liquid velocity and turbulent intensity in bubbly flow was developed based on removing the signal due to the passing bubbles. The local void fraction, interfacial area concentration, interfacial velocity, Sauter mean diameter, liquid velocity, and turbulent intensity of vertical upward air-water flow in a round tube with an inner diameter of 50.8 mm were measured by using these methods. A total of 54 data sets were acquired consisting of three superficial gas flow rates, 0.015-0.076 m s-1, and three superficial liquid flow rates, 0.600, 1.00, and 1.30 m s-1. The measurements were performed at the three locations: L/D=2, 32, and 62. This data is expected to be used for the development of reliable constitutive relations which reflect the true transfer mechanisms in two-phase flow. (orig.)
Local measurement of interfacial area, interfacial velocity and liquid turbulence in two-phase flow
Double sensor probe and hotfilm anemometry methods were developed for measuring local flow characteristics in bubbly flow. The formulation for the interfacial area concentration measurement was obtained by improving the formulation derived by Kataoka and Ishii. The assumptions used in the derivation of the equation were verified experimentally. The interfacial area concentration measured by the double sensor probe agreed well with one by the photographic method. The filter to validate the hotfilm anemometry for measuring the liquid velocity and turbulent intensity in bubbly flow was developed based on removing the signal due to the passing bubbles. The local void fraction, interfacial area concentration, interfacial velocity, Sauter mean diameter, liquid velocity, and turbulent intensity of vertical upward air-water flow in a round tube with inner diameter of 50.8 mm were measured by using these methods. A total of 54 data sets were acquired consisting of three superficial gas flow rates, 0.039, 0.067, and 0.147 m/s, and three superficial liquid flow rates, 0.60, 1.00, and 1.30 m/s. The measurements were performed at the three locations: L/D=2, 32, and 62. This data is expected to be used for the development of reliable constitutive relations which reflect the true transfer mechanisms in two-phase flow. (author)
Bubble Generation in a Flowing Liquid Medium and Resulting Two-Phase Flow in Microgravity
Pais, S. C.; Kamotani, Y.; Bhunia, A.; Ostrach, S.
1999-01-01
The present investigation reports a study of bubble generation under reduced gravity conditions, using both a co-flow and a cross-flow configuration. This study may be used in the conceptual design of a space-based thermal management system. Ensuing two-phase flow void fraction can be accurately monitored using a single nozzle gas injection system within a continuous liquid flow conduit, as utilized in the present investigation. Accurate monitoring of void fraction leads to precise control of heat and mass transfer coefficients related to a thermal management system; hence providing an efficient and highly effective means of removing heat aboard spacecraft or space stations. Our experiments are performed in parabolic flight aboard the modified DC-9 Reduced Gravity Research Aircraft at NASA Lewis Research Center, using an air-water system. For the purpose of bubble dispersion in a flowing liquid, we use both a co-flow and a cross-flow configuration. In the co-flow geometry, air is introduced through a nozzle in the same direction with the liquid flow. On the other hand, in the cross-flow configuration, air is injected perpendicular to the direction of water flow, via a nozzle protruding inside the two-phase flow conduit. Three different flow conduit (pipe) diameters are used, namely, 1.27 cm, 1.9 cm and 2.54 cm. Two different ratios of nozzle to pipe diameter (D(sub N))sup * are considered, namely (D(sub N))sup * = 0.1 and 0.2, while superficial liquid velocities are varied from 8 to 70 cm/s depending on flow conduit diameter. It is experimentally observed that by holding all other flow conditions and geometry constant, generated bubbles decrease in size with increase in superficial liquid velocity. Detached bubble diameter is shown to increase with air injection nozzle diameter. Likewise, generated bubbles grow in size with increasing pipe diameter. Along the same lines, it is shown that bubble frequency of formation increases and hence the time to detachment of a
The purpose of this project is to investigate how changes in interfacial chemical properties affect two-phase transport relationships. Specifically, the objective is to develop a quantitative means that will enable the prediction of changes in the capillary pressure-saturation relationship, a fundamental constitutive relationship in multiphase flow, from changes in interfacial properties, such as adsorption and electrophoretic mobility, through a knowledge of their effect on wettability. The information presented here summarizes the progress we have made in the past eight months of the second project period. Working with a system composed of air-water-silica-cetyltrimethylammonium bromide (CTAB), we have obtained a relationship between degree of adsorption and the surface charge of silica (as measured by electrophoretic mobility), and the drainage and imbibition capillary pressure relationships of system. The bulk of this report describes the completed set of measurements for the air-water-silica-CTAB system at pH 6. We are currently working on a comparable set of measurements for the xylene-water-silica-CTAB system at pH 6. Described here are the interfacial tension, contact angle and preliminary drainage capillary pressure measurements. Our work to date shows a dependence of surface properties on pH. Consequently, in the coming year, we will also complete a set of measurements at another pH value to show the effect of pH on capillary pressure relationships
Buffeting lift forces and local air-water flow aspects around a rigid cylinder
A new experimental programme is conducted in order to relate the characteristics of two-phase flow around a rigid cylinder with the resulting lift forces. The local characteristics of air-water flow measured in the vicinity of the cylinder provide a useful source of information about the effects of flow on the excitation mechanisms. In particular, a selection of relevant parameters has been identified which, with the help of a standard dimensional analysis, may explain the energetic contents of buffeting forces. Among the parameters effective in reducing the data are the flow regime, bubble frequency and gravity forces. In addition, in the range of bubbly regimes, the magnitude of the random forces is found to be related to the local fluctuations of void fraction. Finally, a new formulation is proposed to collapse the dimensionless spectra of the buffeting lift forces in a single characteristic curve. This analysis shows a marked improvement over the collapse of data in comparison with previous normalized models. (authors)
Final Report - Advanced Conceptual Models for Unsaturated and Two-Phase Flow in Fractured Rock
Nicholl, Michael J.
2006-07-10
The Department of Energy Environmental Management Program is faced with two major issues involving two-phase flow in fractured rock; specifically, transport of dissolved contaminants in the Vadose Zone, and the fate of Dense Nonaqueous Phase Liquids (DNAPLs) below the water table. Conceptual models currently used to address these problems do not correctly include the influence of the fractures, thus leading to erroneous predictions. Recent work has shown that it is crucial to understand the topology, or ''structure'' of the fluid phases (air/water or water/DNAPL) within the subsurface. It has also been shown that even under steady boundary conditions, the influence of fractures can lead to complex and dynamic phase structure that controls system behavior, with or without the presence of a porous rock matrix. Complicated phase structures within the fracture network can facilitate rapid transport, and lead to a sparsely populated and widespread distribution of concentrated contaminants; these qualities are highly difficult to describe with current conceptual models. The focus of our work is to improve predictive modeling through the development of advanced conceptual models for two-phase flow in fractured rock.
Two-phase flow heat transfer in nuclear reactor systems
Complete text of publication follows: Heat transfer and phase change phenomena in two-phase flows are often encountered in nuclear reactor systems and are therefore of paramount importance for their optimal design and safe operation.The complex phenomena observed especially during transient operation of nuclear reactor systems necessitate extensive theoretical and experimental investigations. This special issue brings seven research articles of high quality. Though small in number, they cover a wide range of topics, presenting high complexity and diversity of heat transfer phenomena in two-phase flow. In the last decades a vast amount of research has been devoted to theoretical work and computational simulations, yet the experimental work remains indispensable for understanding of two-phase flow phenomena and for model validation purposes. This is reflected also in this issue, where only one article is purely experimental, while three of them deal with theoretical modelling and the remaining three with numerical simulations. The experimental investigation of the critical heat flux (CHF) phenomena by means of photographic study is presented in the paper of J. Park et al. They have used a high-speed camera system to observe the transient boiling characteristics on a thin horizontal cylinder submerged in a pool of water or highly wetting liquid. Experiments show that the initial boiling process is strongly affected by the properties and wettability of the liquid. The authors have stressed the importance of the local scale observation leading to better understanding of the transient CHF phenomena. In the article of G. Espinosa-Paredes et al. a theoretical work concerning the derivation of transport equations for two-phase flow is presented. The author proposes a novel approach based on derivation of nonlocal volume averaged equations which contain new terms related to nonlocal transport effects. These non-local terms act as coupling elements between the phenomena
Phase appearance or disappearance in two-phase flows
Cordier, Floraine; Kumbaro, Anela
2011-01-01
This paper is devoted to the treatment of specific numerical problems which appear when phase appearance or disappearance occurs in models of two-phase flows. Such models have crucial importance in many industrial areas such as nuclear power plant safety studies. In this paper, two outstanding problems are identified: first, the loss of hyperbolicity of the system when a phase appears or disappears and second, the lack of positivity of standard shock capturing schemes such as the Roe scheme. After an asymptotic study of the model, this paper proposes accurate and robust numerical methods adapted to the simulation of phase appearance or disappearance. Polynomial solvers are developed to avoid the use of eigenvectors which are needed in usual shock capturing schemes, and a method based on an adaptive numerical diffusion is designed to treat the positivity problems. An alternate method, based on the use of the hyperbolic tangent function instead of a polynomial, is also considered. Numerical results are presente...
Equations of two-phase flow in spray chamber
李新禹; 张志红; 金星; 徐杰
2009-01-01
The downstream water-air heat and moisture transfer system in a moving coordinate was studied. The relationship between the diameter of the misted droplets and the spray pressure was determined. Based on the theory of the relative velocity,the two-phase flow mode of the spray chamber and the efficiency equation for heat and moisture exchange were established. Corrections were carried out for the efficiency equation with spray pressure of 157 kPa. The results show that the pressure plays an important part in determining the efficiency of heat and moisture exchange. When the spray pressure is less than 157 kPa,better coincidence is noticed between the theoretical analysis and the test results with the error less than 6%. Greater error will be resulted in the case when the spray pressure is beyond 157 kPa. After the correction treatment,the coincidence between the theoretical and the experimental results is greatly improved.
Unsteady interfacial coupling of two-phase flow models
The primary coolant circuit in a nuclear power plant contains several distinct components (vessel, core, pipes,...). For all components, specific codes based on the discretization of partial differential equations have already been developed. In order to obtain simulations for the whole circuit, the interfacial coupling of these codes is required. The approach examined within this work consists in coupling codes by providing unsteady information through the coupling interface. The numerical technique relies on the use of an interface model, which is combined with the basic strategy that was introduced by Greenberg and Leroux in order to compute approximations of steady solutions of non-homogeneous hyperbolic systems. Three different coupling cases have been examined: (i) the coupling of a one-dimensional Euler system with a two-dimensional Euler system; (ii) the coupling of two distinct homogeneous two-phase flow models; (iii) the coupling of a four-equation homogeneous model with the standard two-fluid model. (author)
Thermo-fluid dynamic theory of two-phase flow
A detailed discussion on the formulation of various mathematical models of two-phase flows based on the conservation laws of mass, momentum and energy is presented. Special emphasises have been put on the local instant formulation and on the time-averaged macroscopic models. Two important models have been presented: 1) the two-fluid model which is formulated by considering each phase separately, and 2) the diffusion model which is formulated by considering the mixture as a whole, thus it is expressed in terms of three mixture conservation equations of mass, momentum and energy with one additional diffusion equation. The present formulation can be used to obtain the wall laws by analyses, i.e., the friction factors and heat transfer coefficients. The high-light of the report is the proposed constitutive laws with special emphasis on the interfacial exchange terms as well as on the turbulent transfer terms
Measurement of two-phase flow momentum with force transducers
This paper reports on two strain-gage-based drag transducers developed to measure two-phase flow in simulated pressurized water reactor (PWR) test facilities. One transducer, a drag body (DB), was designed to measure the bidirectional average momentum flux passing through an end box. The second drag sensor, a break through detector (BTD), was designed to sense liquid downflow from the upper plenum to the core region. After prototype sensors passed numerous acceptance tests, transducers were fabricated and installed in two experimental test facilities, one in Japan and one in West Germany. High-quality data were extracted from both the DBs and BTDs for a variety of loss-of-coolant (LOCA) scenarios. The information collected from these sensors has added to the understanding of the thermohydraulic phenomena that occur during the refill/reflood stage of a LOCA in PWR
Numerical modeling of two-phase transonic flow
Halama, Jan; Benkhaldoun, F.; Fořt, Jaroslav
2010-01-01
Roč. 80, č. 88 (2010), s. 1624-1635. ISSN 0378-4754 Grant ostatní: GA ČR(CZ) GA201/08/0012 Institutional research plan: CEZ:AV0Z20760514 Keywords : two-phase flow * condensation * fractional step method Subject RIV: BK - Fluid Dynamics Impact factor: 0.812, year: 2010 http://www.sciencedirect.com/science?_ob=MImg&_imagekey=B6V0T-4VNK68X-2-R&_cdi=5655&_user=640952&_pii=S0378475409000421&_origin=search&_coverDate=04%2F30%2F2010&_sk=999199991&view=c&wchp=dGLzVlb-zSkWb&md5=5ba607428fac339a3e5f67035d3996d0&ie=/sdarticle.pdf
Test of the EG and G two-phase mass flow rate instrumentation at Kernforschungszentrum Karlsruhe
This report presents the data analyses of experiments designed to understand the behavior of a free field drag disc turbine transducer (DTT) and a three beam γ densitometer in steady-state horizontal steam-water and air-water flow. The pressure was varied between 2 and 75 bars, the experiments were made at a mass flow rate and void fraction range where various quite separated flow regimes occurred. Two different test sections with 103 mm ID (5' pipe) and 66 mm ID (3' pipe) were used. Information on flow regime and phase distribution in the cross section was obtained with local impedance probes, measurements of the axial distribution of phase velocities in the test section piping were made with the radiotracer technique. The best overall accuracy of mass flow rate determined by combining two of the three available instruments is obtained by the combination of γ-densitometer and drag disc. From the experiments, single calibration factors are determined which depend only on the γ-densitometer reading. A time averaged separated two-phase model for the DTT is postulated which shows that the DTT measures the local parameters. To obtain the pipe averaged mass flux, a density correction is proposed. For some experiments the radiotracer technique combined with the γ-densitometer for measuring the mass flow rate was tested. (orig./HP) 891 HP/orig.- 892 HIS
Transition from slug to annular flow in horizontal air-water flow
The transition from slug to annular flow in horizontal air-water and steam-water flow was investigated. Test sections of 50; 66.6 and 80 mm ID were used. The system pressure was 0.2 and 0.5 MPa in the air-water experiments and 2.5; 5; 7.5 and 10 MPa in the steam-water experiments. For flow pattern detection local impedance probes were used. This method was compared in a part of the experiments with differential pressure and gamma-beam measurements. The flow regime boundary is shifting strongly to smaller values of the superficial gas velocity with increasing pressure. Correlations from literature fit unsatisfactorily the experimental results. A new correlation is presented. (orig.)
Unsteady flow analysis of a two-phase hydraulic coupling
Hur, N.; Kwak, M.; Lee, W. J.; Moshfeghi, M.; Chang, C.-S.; Kang, N.-W.
2016-06-01
Hydraulic couplings are being widely used for torque transmitting between separate shafts. A mechanism for controlling the transmitted torque of a hydraulic system is to change the amount of working fluid inside the system. This paper numerically investigates three-dimensional turbulent flow in a real hydraulic coupling with different ratios of charged working fluid. Working fluid is assumed to be water and the Realizable k-ɛ turbulence model together with the VOF method are used to investigate two-phase flow inside the wheels. Unsteady simulations are conducted using the sliding mesh technique. The primary wheel is rotating at a fixed speed of 1780 rpm and the secondary wheel rotates at different speeds for simulating different speed ratios. Results are investigated for different blade angles, speed ratios and also different water volume fractions, and are presented in the form of flow patterns, fluid average velocity and also torques values. According to the results, blade angle severely affects the velocity vector and the transmitted torque. Also in the partially-filled cases, air is accumulated in the center of the wheel forming a toroidal shape wrapped by water and the transmitted torque sensitively depends on the water volume fraction. In addition, in the fully-filled case the transmitted torque decreases as the speed ration increases and the average velocity associated with lower speed ratios are higher.
Particle migration in two-phase, viscoelastic flows
Jaensson, Nick; Hulsen, Martien; Anderson, Patrick
2014-11-01
Particles suspended in creeping, viscoelastic flows can migrate across stream lines due to gradients in normal stresses. This phenomenon has been investigated both numerically and experimentally. However, particle migration in the presence of fluid-fluid interfaces is hardly studied. We present results of simulations in 2D and 3D of rigid spherical particles in two-phase flows, where either one or both of the fluids are viscoelastic. The fluid-fluid interface is assumed to be diffuse and is described using Cahn-Hilliard theory. The particle boundary is assumed to be sharp and is described by a boundary-fitted, moving mesh. The governing equations are solved using the finite element method. We show that differences in normal stresses between the two fluids can induce a migration of the particle towards the interface in a shear flow. Depending on the magnitude of the surface tension and the properties of the fluids, particle migration can be halted due to the induced Laplace pressure, the particle can be adsorbed at the interface, or the particle can cross the interface into the other fluid. Dutch Polymer Institute (DPI), P.O. Box 902, 5600 AX Eindhoven, The Netherlands.
Critical two-phase flow through rough slits
The knowledge of the two-phase mass flow rate through a crack in the wall of nuclear or chemical reactor components is very important under the leak-before-break criterion point of view. For providing a qualified analytical tool for calculating critical mass flow rates through such a crack a detailed test program was carried out using subcooled water up to pressures of 14 MPa. A real crack and several simulated cracks (rough slits) were examined experimentally. The important parameters such as inlet pressure, subcooling temperature of water, slit width, and inner surface roughness were varied in a wide range and the measured data compared with calculated values from different models. The data comparison indicates that the model published by Pana leads to predictions which agree best with the observed data. First calculations were carried out using the friction coefficient ζ, which results from single phase flow measurements. A correlation has been developed to calculate ζ from the geometrical dimensions of the crack and was integrated into an advanced version of the Pana model. The modified Pana model was qualified against some hundreds of test values. The measured data were predicted with a relative standard deviation of less than 20%. (orig.)
Uncertainty analysis of two-phase flow pressure drop calculations
Siqueira, Cezar A.M.; Costa, Bruno M.P.; Fonseca Junior, Roberto da; Gonalves, Marcelo de A.L. [PETROBRAS, Rio de Janeiro, RJ (Brazil)
2004-07-01
The simulation of multiphase flow in pipes is usually performed by petroleum engineers with two main purposes: design of new pipelines and production systems; diagnosis of flow assurance problems in existing systems. The tools used for this calculation are computer codes that use published pressure drop correlations developed for steady-state two-phase flow, such as Hagedorn-Brown, Beggs and Brill and others. Each one of these correlations is best suited for a given situation and the engineer must find out the best option for each particular case, based on his experience. In order to select the best correlation to use and to analyze the results of the calculation, the engineer must determine the reliability of computed values. The uncertainty of the computation is obtained by considering uncertainties of the correlation adopted, of the calculation algorithm and the input data. This paper proposes a method to evaluate the uncertainties of this type of calculation and presents an analysis of these uncertainties. The uncertainty analysis also allows the identification of the parameters that are more significant for the final uncertainty of the simulation. Therefore it makes possible to determine which are the input parameters that must be determined with higher accuracy and the ones that may have lower accuracy, without reducing the reliability of the results. (author)
Test of the EG and G two-phase mass flow rate instrumentation at Kernforschungszentrum Karlsruhe
For many experiments which investigate the Loss-of-Coolant Accident (LOCA) in nuclear reactors, proper measurement of the two-phase mass flow is of great importance. This report presents the experimental description and the data of experiments designed to understand the behaviour of a free field drag disc turbine transducer (DTT) and a three beam gamma densitometer in steady-state horizontal steam-water and air-water flow. The pressure was varied between 2 and 75 bars, the experiments were made at a mass flow rate and void fraction range where various quite separated flow regimes occurred. Two different test sections with 103 mm ID (5 pipe) and 66 mm ID (3 pipe) were used. Information on flow regime and phase distribution in the cross section was obtained with local impedance probes, measurements of the axial distribution of phase velocities in the test section piping were made with the radiotracer technique. These techniques are of great help for the physical interpretation of the single instrument readings. The results of detailed data analyses are given in another report. (orig.)
Simulation of two-phase flows by domain decomposition
This thesis deals with numerical simulations of compressible fluid flows by implicit finite volume methods. Firstly, we studied and implemented an implicit version of the Roe scheme for compressible single-phase and two-phase flows. Thanks to Newton method for solving nonlinear systems, our schemes are conservative. Unfortunately, the resolution of nonlinear systems is very expensive. It is therefore essential to use an efficient algorithm to solve these systems. For large size matrices, we often use iterative methods whose convergence depends on the spectrum. We have studied the spectrum of the linear system and proposed a strategy, called Scaling, to improve the condition number of the matrix. Combined with the classical ILU pre-conditioner, our strategy has reduced significantly the GMRES iterations for local systems and the computation time. We also show some satisfactory results for low Mach-number flows using the implicit centered scheme. We then studied and implemented a domain decomposition method for compressible fluid flows. We have proposed a new interface variable which makes the Schur complement method easy to build and allows us to treat diffusion terms. Using GMRES iterative solver rather than Richardson for the interface system also provides a better performance compared to other methods. We can also decompose the computational domain into any number of sub-domains. Moreover, the Scaling strategy for the interface system has improved the condition number of the matrix and reduced the number of GMRES iterations. In comparison with the classical distributed computing, we have shown that our method is more robust and efficient. (author)
Development of two-phase Flow Model, 'SOBOIL', for Sodium
The objective of this research is to develop a sodium two-phase flow analysis model, 'SOBOIL', for the assessment of the initial stage of the KALIMER HCDA (Hypotherical Core Disruptive Accident). The 'SOBOIL' is basically similar to the multi-bubble slug ejection model used in SAS2A[1]. When a bubble is formed within the liquid slug, the bubble fills the whole cross section of the coolant channel except for a film left on the cladding or on the structure. Up to nine bubbles, separated by the liquid slugs, are allowed in the channel at any time. Each liquid slug flow rate in the model is performed in 2 steps. In the first step, the preliminary flow rate in the liquid slug is calculated neglecting the effect of changes in the vapor bubble pressures over the time step. The temperature and pressure distributions, and interface velocity at the interface between the liquid slug and vapor bubble are also calculated during this process. The new vapor temperature and pressure are then determined from the balance between the net energy transferred into the vapor and the change of the vapor energy. The liquid flow is finally calculated considering the change of the vapor pressure over a time step and the calculation is repeated until specified elapsed time is met. Continuous effort, therefore, must be made on the examination and improvement for the model to become reliable. To this end, much interest must be concentrated in the relevant international collaborations for access to a reference model or test data for the verification
Characterization and modeling of annular two-phase flows
Three aspects of annular two-phase flow are studied: (a) wave motion on falling films, (b) flow transition from downflow to upflow, and (c) the upflow. For the case of wave motion on falling films, it is shown that the assumption of the Nusselt velocity profile for finite-amplitude waves is solution of the wave profile, wave velocity, and velocity components within the wave is developed. An algorithm based on collocation methods is also detailed and can be applied to extend the model to solve for higher order terms in the velocity profile. Comparisons with experimental studies show good agreement. Flow transition and the upflow experiments are conducted in a 5.08 x 10-2m inner diameter, 6.5m long Plexiglas column. The liquid rates are varied from 0 to 0.126 kg/s and the gas rates from 0 to 0.0524 kg/s. At four measuring stations along the length of the column, an electrical conductance technique which employs two electrodes mounted flush with the wall is utilized to measure film thickness and pressure transducers are used to make the pressure measurements. Flow visualization studies indicate that flooding takes place as a result of entrainment from the crests of large waves. The effect of column length and pore size of the feed device on flooding velocities is studied. No previous correlation or theory is found to be fully adequate. A speculative interaction among system parameters is proposed to form a basis for a physical model for flooding phenomena
Numerical methods for two-phase flow with contact lines
Walker, Clauido
2012-07-01
This thesis focuses on numerical methods for two-phase flows, and especially flows with a moving contact line. Moving contact lines occur where the interface between two fluids is in contact with a solid wall. At the location where both fluids and the wall meet, the common continuum descriptions for fluids are not longer valid, since the dynamics around such a contact line are governed by interactions at the molecular level. Therefore the standard numerical continuum models have to be adjusted to handle moving contact lines. In the main part of the thesis a method to manipulate the position and the velocity of a contact line in a two-phase solver, is described. The Navier-Stokes equations are discretized using an explicit finite difference method on a staggered grid. The position of the interface is tracked with the level set method and the discontinuities at the interface are treated in a sharp manner with the ghost fluid method. The contact line is tracked explicitly and its dynamics can be described by an arbitrary function. The key part of the procedure is to enforce a coupling between the contact line and the Navier-Stokes equations as well as the level set method. Results for different contact line models are presented and it is demonstrated that they are in agreement with analytical solutions or results reported in the literature.The presented Navier-Stokes solver is applied as a part in a multiscale method to simulate capillary driven flows. A relation between the contact angle and the contact line velocity is computed by a phase field model resolving the micro scale dynamics in the region around the contact line. The relation of the microscale model is then used to prescribe the dynamics of the contact line in the macro scale solver. This approach allows to exploit the scale separation between the contact line dynamics and the bulk flow. Therefore coarser meshes can be applied for the macro scale flow solver compared to global phase field simulations
CFD code validation against stratified air-water flow experimental data
Pressurized Thermal Shock (PTS) modelling has been identified as one of the most important industrial needs related to nuclear reactor safety. A severe PTS scenario limiting the Reactor Pressure Vessel (RPV) lifetime is the cold water Emergency Core Cooling (ECC) injection into the cold leg during a Loss of Coolant Accident (LOCA). Since it represents a big challenge for numerical simulations, this scenario was selected within the NURESIM (European Platform for Nuclear Reactor Simulations) Integrated Project as a reference two-phase problem for CFD code validation. This paper presents a CFD analysis of a stratified air-water flow experimental investigation performed at the Institut de Mecanique des Fluides de Toulouse in 1985 [1], which shares some common physical features with the ECC injection in PWR cold leg. Numerical simulations have been carried out with two commercial codes (Fluent and Ansys CFX), and a research code NEPTUNECFD (developed by EDF and CEA). The aim of this work, carried out at the University of Pisa within the NURESIM IP, is to validate the free surface flow model implemented in the codes against the available experimental data, and to perform code to code benchmarking. Obtained results suggest the relevance of three-dimensional effects and stress the importance of a suitable interface drag coefficient modelling. A relevant improvement of results has been achieved with 3D simulations, even if the air velocity profile was still significantly underestimated. (author)
Study on two-phase flow dynamics in steam injectors
Analytical and experimental studies have been conducted on large-scale steam injectors for a next-generation reactor. The steam injectors are simple, compact, passive steam jet pumps for a steam-injector-driven passive core injection system (SI-PCIS) or steam-injector-driven primary loop recirculation system (SI-PLR). In order to check the feasibility of such large-scale steam injectors, we developed the separate-two-phase flow models installed in the PHOENICS Code, and scale-model tests were conducted for both SI-PCIS and SI-PLR. A 1/2 scale SI-PCIS model achieved a discharge pressure of almost 8 MPa with 7 MPa steam and 0.4 MPa water, and a 1/5 scale SI-PLR model attained a discharge pressure of 12.5 MPa with 3 MPa steam and 7 MPa water. Both results are in good agreement with the analysis, confirming the feasibility of both systems. The systems will help to simplify the next generation of BWRs. (author)
DNS and LES of two-phase flows with cavitation
Hickel, Stefan
2014-01-01
We report on recent progress in the physical and numerical modeling of compressible two-phase flows that involve phase transition between the liquid and gaseous state of the fluid. The high-speed dynamics of cavitation bubbles is studied in well-resolved simulations (DNS) with a sharp-interface numerical model on a micro scale. The underlying assumption of the employed evaporation/condensation model is that phase change occurs in thermal non-equilibrium and that the associated timescale is larger than that of the wave dynamics. Results for the collapse of a spherical vapor bubble close to a solid wall are discussed for three different bubble-wall configurations. The major challenge for such numerical investigations is to accurately reproduce the dynamics of the interface between liquid and vapor during the entire collapse process, including the high-speed dynamics of the late stages, where compressibility of both phases plays a decisive role. Direct interface resolving simulations are intractable for real wor...
Hybrid dynamic modeling for two phase flow condensers
In this paper, a hybrid modeling approach is proposed to describe the dynamic behavior of the two phase flow condensers used in air-conditioning and refrigeration systems. The model is formulated based on fundamental energy and mass balance governing equations, and thermodynamic principles, while some constants and less important variables that change very little during normal operation, such as cross-sectional areas, mean void fraction, the derivative of the saturation enthalpy with respect to pressure, etc., are lumped into several unknown parameters. These parameters are then obtained by experimental data using least squares identification method. The proposed modeling method takes advantages of both physical and empirical modeling approaches, can accurately predict the transient behaviors in real-time and significantly reduce the computational burden. Other merits of the proposed approach are that the order of the model is very low and all the state variables can be easily measured. These advantages make it easy to be applied to model based control system design. The model validation studies on an experimental system show that the model predicts the system dynamic well. -- Highlights: • A hybrid modeling approach is proposed to describe the dynamic behavior of condensers. • This modeling approach balances the trade-offs between complexity and accuracy. • The model order is very low and all the state variables are available for measurement. • The model validation studies show that the model predicts the system dynamic well. • The model is suitable for dynamic analysis and model-based controller design
Experimental and Theoretical Studies on Two-Phase Flows.
Koh, Christopher James
This thesis, comprised of two parts, deals with the flow of suspensions. Part I concerns specifically with the stability of a single drop translating through a quiescent, unbounded suspending fluid at low Reynolds number. The evolution of the shape of an initially nonspherical drop as it translates is studied numerically and experimentally. For finite capillary numbers, it is shown that the drop reverts to a sphere provided that the initial deformation is small enough. However, beyond certain critical initial deformation, the drop deforms continuously. For initially prolate shapes, the drop elongates with the formation of a tail; for initially oblate shapes, the drop flattens with the formation of a cavity at its rear. Experiments extend the numerical results. It is found that initially unstable prolate drops break up into multiple droplets, while initially unstable oblate drops deform in double-emulsion drops. Part II of this thesis considers the flow of high concentration solid suspensions through a rectangular channel. By adapting the well-known Laser Doppler Anemometry, an experimental technique is developed to measure the velocity as well as particle volume fraction of the suspension. A crucial element in these experiments is the reduction of the optical turbidity of the suspension. To accomplish this goal, a systematic method based on refractive-index-matching of the two phases is employed. Experimental results show that the velocity profile is blunted while the concentration profile has a maximum near the center. The qualitative features of the experimental data compare reasonably well with theoretical predictions based on the shear-induced particle migration theory.
Dynamic modeling strategy for flow regime transition in gas-liquid two-phase flows
In modeling gas-liquid two-phase flows, the concept of flow regime has been widely used to characterize the global interfacial structure of the flows. Nearly all constitutive relations that provide closures to the interfacial transfers in two-phase flow models, such as the two-fluid model, are often flow regime dependent. Currently, the determination of the flow regimes is primarily based on flow regime maps or transition criteria, which were developed for steady-state, fully-developed flows and have been widely applied in nuclear reactor system safety analysis codes. As two-phase flows are dynamic in nature (fully-developed two-phase flows generally do not exist in real applications), it is of importance to model the flow regime transition dynamically to be able to predict two-phase flows more accurately. The present work aims to develop a dynamic modeling strategy to determine flow regimes in gas-liquid two-phase flows through introduction of interfacial area transport equations (IATEs) within the framework of a two-fluid model. The IATE is a transport equation that models the interfacial area concentration by considering the creation and destruction of the interfacial area, such as the fluid particle (bubble or liquid droplet) disintegration, boiling and evaporation; and fluid particle coalescence and condensation, respectively. For the flow regimes beyond bubbly flows, a two-group IATE has been proposed, in which bubbles are divided into two groups based on their size and shapes, namely group-1 and group-2 bubbles. A preliminary approach to dynamically identify the flow regimes is discussed, in which discriminators are based on the predicted information, such as the void fraction and interfacial area concentration. The flow regime predicted with this method shows good agreement with the experimental observations. (author)
Turbulent mixing rates of gas and liquid phases between adjacent subchannels have been measured for various air-water two-phase annular flows in a multiple channel consisting of two identical circular subchannels. In order to study effect of liquid entrainment in the gas core on the turbulent mixing rates, experiments were conducted for two types of liquid injection method, i.e., a small bore nozzle placed in the subchannel center and a porous wall. The result showed that the effect of liquid entrainment on the turbulent mixing rates of both phases is negligibly small. (author)