Experimental Study on Interfacial Area Transport of Two-Phase Flow under Vibration Conditions

An experimental study on air-water two-phase flow under vibration condition has been conducted using double-sensor conductivity probe. The test section is an annular geometry with hydraulic diameter of 19.1mm. The vibration frequency ranges from 0.47Hz to 2.47Hz. Local measurements of void fraction, interfacial area concentration (IAC), and Sauter mean diameter have been performed along one radius in the vibration direction. The result shows that local parameters fluctuate continuously around the base values in the vibration cycle. Additional bubble force due to inertia is used to explain lateral bubble motions. The fluctuation amplitudes of local void fraction and IAC increase significantly with vibration frequency.The radial distribution of local parameters at the maximum vibration displacement is specifically analyzed. In the void fraction and IAC profiles, the peak near the inner wall is weakened or even disappearing and a strong peak skewed to outer wall is gradually observed with the increase of vibration frequency. The nondimensional peak void fraction can reach a maximum of 49% and the mean relative variation of local void fraction can increase to more than 29% as the vibration frequency increases to 2.47Hz. But the increase of vibration frequency does not bring significant change to bubble diameter.


Introduction
Two-phase flows are frequently encountered in industrial applications.In the present state of the arts, the two-fluid model is considered to be the most detailed and accurate formulation in analyzing the thermofluid dynamics of twophase systems.In two-fluid model, liquid and gas phases are separately described by using two sets of conservation equations for mass, momentum, and energy.The existence of the interfacial transfer terms which represent the mass, momentum, and energy transfer at the liquid-vapor interface is one of the most important characteristics of the two-fluid model formulation.To dynamically predict the interfacial transport phenomena, Kocamustafaogullari and Ishii [1] proposed the interfacial area transport equation (IATE).The IATE is based on statistical mechanics and its concept has been fully established.But the sink and source terms of interfacial area due to bubble interactions are still being developed.
Recent advances in two-phase measurement techniques have stimulated further research on interfacial area transport models.Particularly, the improvements in probe methods [2][3][4][5][6] allow accurate measurements of local flow parameters, such as void fraction, interfacial area concentration, and interfacial velocity.With the help of these techniques, a rigorous database has been collected for two-phase flow in various flow geometries [7][8][9][10].
However, all the above-mentioned experiments were performed under static condition.Earthquakes struck to nuclear reactors frequently in recent years and some even resulted in scram in nuclear power plants [11,12].Many studies have been conducted in consideration of the nuclear reactor safety when an earthquake occurs, including the structural safety design of nuclear reactors and siting of nuclear power plants.But there is little knowledge of thermal fluid characteristics in the reactor system during the earthquake.Two-phase flow is an important phenomenon in reactor core and its behavior is closely related to the operation of nuclear reactors.Due to the different properties of the two phases and the complex interactions at phase boundary, the seismic vibration may significantly change the two-phase flow structures in reactor core.This change may result in deterioration in local heat transfer, even critical heat flux (CHF).In addition, the changes in void fraction and interfacial area concentration may further bring reactivity feedback to reactor core and change neutron flux as well as power distribution in the reactor.Therefore, in order to understand the effect of earthquake vibration on nuclear reactor precisely, the study on two-phase flow behaviors under vibration condition is required.Nariai and Tanaka [13] performed experimental study of void fraction in subcooled boiling two-phase flow under oscillating conditions.The void fraction decreased drastically when the vibration frequency reached 10 Hz or higher.Hibiki and Ishii [14] investigated the effects of flow-induced vibration on local parameters of upward air-water two-phase flow.The flow-induced vibration promoted bubble coalescence and led to increase of Sauter mean diameter and decrease of interfacial area concentration for relatively low superficial liquid velocity.The void fraction profiles also changed from "wall peak" to "core peak" or "transition."Lee et al. [15] performed CHF experiments with vibrating test section.The results showed that the CHF was enhanced by increasing vibration amplitude and frequency.Chen et al. [16] studied the adiabatic two-phase flow under vibration conditions with impedance meters.It was found that the low-frequency vibration had a limited impact on the time-averaged area-averaged void fraction.However, the instantaneous void fraction signals were affected and the results were shown by fast Fourier transformation (FFT) analysis.Arai et al. [17] performed visualization experiments for vertical bubbly flow in rectangular water tank on which sine wave vibration was applied.The bubble behaviors were affected by oscillation and the cycle of bubble inclination angel was almost the same as the cycle of oscillation table velocity.Xiao et al. [18] conducted preliminary experiment of vibration effect on bubbly two-phase flow.Important phenomena including the variations of local void fraction and interfacial area concentration were proposed.But the detailed influence of vibration parameters on two-phase flow structure was not addressed.
As described above, some studies have been conducted to investigate the vibration effect on two-phase flow.The results indicate that the flow behaviors will be affected by structure vibration.However, experimental data concerning about the detailed variations of local interfacial parameters under vibration condition is still lacking.In the current research, an experiment has been performed to measure the void fraction, interfacial area concentration, and Sauter mean diameter locally using double-sensor conductivity probe under different vibration conditions.

Experimental Apparatus
2.1.Experimental Facility.The schematic of the experimental facility is shown in Figure 1.The test section is an annular geometry that is formed by a Pyrex glass tube on the outside and a stainless steel rod on the inside.The inner rod has a 19.1 mm outer diameter and the glass tube has a 38.1 mm inner diameter.The overall length of the flow channel is 2.32 m.A two-phase mixture injection system is located at the bottom of the test section, which consists of a sparger and water flow channel.The sparger is made up of porous material with an average pore size of 10 m and can generate almost uniform bubbles of approximately 1-2 mm in diameter.During the experiment, filtered and chemically treated water with conductivity of about 300 S/cm was held in the main tank.Water was driven by a centrifugal pump and controlled by a ball valve and a frequency converter.Air was supplied by a compressor and a pressure regulator was employed to maintain a constant pressure.After flowing through the test section, the air was released into the atmosphere through a separator while water was circulated by the centrifugal pump.The superficial liquid velocity ⟨  ⟩ and superficial gas velocity ⟨  ⟩ ranged from 0.52 m/s to 1.03 m/s and from 0.195 m/s to 0.403 m/s, respectively.Here, the symbol ⟨ ⟩ represents the area-averaged value.The water flow rate was measured using an electromagnetic flow meter with an accuracy of ±1%.Two rotameters with different measuring ranges (0.08-0.57m 3 /h and 0.28-2.83m 3 /h) were used to measure the gas flow rate with an accuracy of ±3% of full scale.The water temperature was kept at 25 ± 0.5 ∘ C.
The annulus test section is connected to a vibration module through a rectangular vibration beam.The vibration module mainly consists of a fixed I-beam, a vibration beam, a ball bearing, a sliding track, three damping springs, a 3-HP motor, a frequency controller, and an eccentric cam.The vibration beam is driven back and forth by the eccentric cam and damping springs.By controlling the motor speed (ms, in rpm) and using different eccentric wheels, various vibration conditions can be acquired.In view of the earthquake intensity scale and safety operation of the vibration module [16], an eccentric cam with eccentricity of 9.5 mm is used.The motor speed ranges from 30 rpm to 150 rpm.An accelerometer is installed on the vibration beam to measure the vibration acceleration () of the test section.The accuracy of the accelerometer is ±1% of full scale (2, where  is the gravitational acceleration and 1 = 9.8 m/s 2 ).

2.2.
Instrumentation.Double-sensor conductivity probe was used to measure local two-phase flow parameters including void fraction, interfacial area concentration, and Sauter mean diameter under both static and vibration conditions.The local time-averaged void fraction was simply the accumulated time during which the sensor was exposed to the gas phase divided by the total sampling time.The bubble interface velocity was calculated by taking into account the distance between the tips of the upstream and downstream sensors to the time difference between the upstream and downstream signal.The interfacial area concentration can be calculated from interface velocity by using some mathematical formulations and statistical considerations.The detailed measurement methodology can be referred to the previous works [3,19].The measurement port was located at / ℎ = 106, where  is the distance from the measurement port to the inlet port and  ℎ is the hydraulic diameter of the flow channel.The port was designed to support conductivity probe and the probe can be radially traversed in the annulus.Probe measurements were taken along one of the radius in the vibration direction as shown in Figure 2. Local flow measurements using doublesensor conductivity probe were performed at eleven radial locations of /(  −   ) = 0.13, 0.18, 0.23, 0.33, 0.43, 0.52, 0.62, 0.72, 0.82, 0.86, and 0.91, where  is the radial distance from the inner rod wall and (  −   ) is the radial clearance of the annulus.Probe signals were collected at a high sampling frequency of 40 kHz.The data acquisition time was adjusted depending on specific flow condition to make sure that at least 1500 bubbles were detected by the leading sensor of the probe in every vibration interval to be defined in Section 3.3.Differential pressure measurement was also taken at the measurement port.The accuracy of the pressure transducer was within ±0.375% of full measurement range.It is important to note that the probe signals must be synchronized with the vibration signal in order to reveal the vibration effect on local flow parameters in time.To make sure of this, the same data acquisition board was used to acquire output signals from the probe and the accelerometer simultaneously.
In order to verify the accuracy of double-sensor probe measurement, the area-averaged quantities obtained by integrating the local parameters over the flow channel were compared with those measured by other cross-calibration methods such as differential pressure gauge for void fraction and rotameter for superficial gas velocity.The comparison results are shown in Table 1.The area-averaged superficial gas velocity from probe in Table 1 was calculated based on the local void fraction and local gas velocity measured by the double-sensor probe.Good agreements were obtained   with averaged relative deviations of 10% and 11.5% for area-averaged void fraction and superficial gas velocity, respectively.The benchmark experiments for double-sensor conductivity probe have also been performed in earlier studies [20,21].The relative difference of interfacial area concentration between the image processing method and doublesensor probe was within 10%.Sauter mean diameter depends on void fraction and interfacial area concentration.The propagated uncertainty for Sauter mean diameter is 19.7%.

Acceleration Analysis.
The vibration frequency of the test section was controlled by the motor speed.Five motor speeds which ranged from 30 rpm to 150 rpm were adopted during the vibration experiment.Vibration acceleration was measured by the accelerometer.Signal filtering and FFT methods were used to process the acceleration signals.The vibration acceleration changes linearly along the length of the test section.In the following, the vibration parameters stand for the vibration condition at the measurement location of / ℎ = 106.Figure 3 shows the vibration acceleration and displacement profiles under different motor speeds.Here, the test section is defined to vibrate around -axis and counterclockwise around -axis is defined as the positive vibration direction.It can be seen from Figure 3 that the vibration parameters vary periodically with the vibration motion.The Sauter mean diameter, formation of big cap bubbles decreases the interfacial area concentration.As a result, the profile of the IAC is different from the void fraction profile.The Sauter mean diameter has a core peak distribution, resulting in a concave profile of the interfacial area concentration near the channel center for this flow condition.
The vibration experiments were then carried out in sequence from low to high frequency.Structure vibration may stimulate phase interactions and interfacial transport in two-phase flow.In order to verify the vibration effect on local time-averaged flow parameters, the parameters measured under vibration conditions are compared with those measured under static condition.The mean relative variation (MRV), which represents the average of local time-averaged flow parameters' variation in the 11 radial measurement locations caused by vibration, is defined.Here, the MRV of void fraction, interfacial area concentration, and Sauter mean diameter have the following expressions: where  is the total number of radial measurement locations and  denotes the th radial measurement location; the subscripts  and  = 0 denote the local time-averaged values measured under vibration condition with vibration frequency of  and under static condition, respectively.In general, the vibration effect on local time-averaged flow parameters is very limited with MRV less than 8%.Particularly, according to the analysis in every radial measurement location, the variation of local time-averaged parameters is within ±10%.This can be explained from two aspects.On one hand, the pressure fluctuation caused by the vibration of the test section is still low compared to the hydrostatic head and friction loss.On the other hand, the instantaneous change in local flow parameters resulting from structure vibration is balanced out by the time averaging process.Therefore, an adequate analytical method is needed to reveal the transient variations of local flow parameters.

Variation of Local Flow Parameters under Vibration
Condition.The test section vibrates periodically with preset frequency during vibration experiment.In order to obtain the transient variations of local two-phase flow parameters under vibration condition, the vibration interval is defined.Specifically, the vibration interval represents a continuously finite interval where the flow behavior is assumed to remain the same in the vibration cycle.Local flow parameters measured by conductivity probe are averaged values over a sufficient amount of time for statistical convergence.Therefore, the same vibration intervals in every vibration cycle are grouped together and the time-averaged flow parameters in such a combined interval are used to stand for the parameters at the instantaneous vibration position.There are two considerations involved in selection of the vibration interval.One is that the interval cannot be too large to keep the above assumption reasonable.The second is that the interval cannot be too small.The total sampling time of each combined interval should be sufficiently large to minimize statistical errors [22].Based on these two considerations, the vibration cycle is divided into twelve vibration intervals with each vibration interval covering /6 vibration phase angle.In the present experiment, the vibration parameter is measured by the accelerometer.The acceleration signal as shown in Figure 3 is first analyzed to provide information for separating vibration cycle into twelve vibration intervals.Since the probe signals are collected simultaneously with the acceleration signal, the probe signals are then separated and regrouped into twelve parts according to the above information.As a result, the local flow parameters at each vibration interval as well as the variation of local parameters in the vibration cycle can be determined.Figure 6 shows the variation of local flow parameters at five representative radial locations in the vibration cycle.As a comparison, the respective local time-averaged flow parameters measured under static condition are presented as dotted lines in this figure.Since the probe measurements were taken only along one radius as indicated in Figure 2, the variation of local parameters along this radius is mainly discussed in this study.Plus, the inner wall and outer wall in the following denote the radial locations where /(  −   ) = 0 and 1, respectively.It can be seen from Figure 6 that the local parameters fluctuate continuously around the dotted lines in the vibration cycle.In the flow region near the inner wall (/(  −   ) = 0.13 and 0.33), the void fraction and interfacial area concentration decrease as the test section moves from  2 to  1 and increase as the test section moves from  1 to  2 .Opposite trend can be observed in the flow region near the outer wall (/(  −   ) = 0.72 and 0.91).The fluctuation amplitude of void fraction and interfacial area concentration around the channel center (/(  −   ) = 0.52) is relatively small.This finding can be attributed to the significant radial migration of bubbles under vibration, which mainly results from the coaction of forces acting on bubbles.Some researchers [23][24][25] have reviewed the forces acting on bubbles in lateral direction under static condition.The lateral motion of bubbles is mainly determined by three forces: the wall force, the lift force, and the turbulent dispersion force.Bubble coalescence and breakup process can also influence bubble motions.In the vibration experiment, there are two additional forces acting on bubbles.One is the lateral component of buoyancy force.It is directed toward outer wall and inner wall, respectively, when the test section is in the positive and negative vibration directions.In fact, this force is negligible since the inclination angle, which is defined as the ratio of vibration displacement to the total length of the test section, is rather small.The other one is the lateral force due to inertia.This force plays an important role in radial bubble motion and its magnitude and direction vary continuously with vibration acceleration.The sum of bubble forces projected to the radial direction directly determines the variation of local flow parameters.
As the test section moves from  = 0 to  =  1 , the vibration acceleration of the test section has a negative direction.The liquid phase tends to accumulate near the inner wall under inertia effect and the bubbles are driven toward outer wall.As a result, the local void fraction decreases in the radial locations of /(  −  ) = 0.13 and 0.33 and increases in the radial locations of /(  −   ) = 0.72 and 0.91 during this stage.With the bubbles moving to the outer wall, the radial gradient of void fraction becomes higher and the turbulent dispersion force is enhanced.As the test section moves from  =  1 to  = 0, the vibration acceleration keeps its negative direction but its value decreases from | 2 | to almost 0.During this stage, the turbulent dispersion force gradually exceeds the additional bubble force caused by inertia and bubbles begin to reverse their direction to move toward inner wall.Therefore, the local void fraction starts to increase near the inner wall region and decrease near the outer wall at a certain moment in this stage.The opposite trend can be expected during the other half period of the vibration cycle.The variation of bubble Sauter mean diameter in the vibration cycle is similar to that of the void fraction and interfacial area concentration, but the variation amplitude is much smaller.Thus, the vibration condition does not lead to much bubble coalescence or breakup in the present experiment.It can also be observed from Figure 7 that the fluctuation amplitude of the local flow parameters increases with vibration frequency.The absolute change of void fraction near the wall can reach more than 17% under vibration frequency of  = 2.473 Hz.This indicates that the flow structure may change from finely dispersed bubbly flow to closely packed bubbly flow and the interfacial transport terms in the two-phase flow may be significantly changed.
It is important to note that, according to the design mechanism of the vibration module, the vibration of the test section is axisymmetric and the symmetry axis is defined as -axis.The vibration acceleration of the test section is approximately symmetric in the positive and negative vibration directions, as shown in Figure 3.As a result, the lateral forces acting on bubbles including the additional force due to inertia and other interphase forces are almost symmetric with respect to -axis.Therefore, the flows can be expected as axisymmetric under vibration condition and the measurement of two-phase flow along one radius is sufficient to reveal the two-phase flow behaviors in the vibration direction during structure vibration.

Vibration Frequency Effect on Variation of Local Flow
Parameters.Two-phase flow structure at the vibration displacement of  =  1 is used to investigate the vibration frequency effect on variation of local flow parameters.Vibration frequency  ranges from 0.472 Hz to 2.473 Hz.The radial distributions of local parameters under different vibration frequency are shown in Figure 7.The local time-averaged flow parameters measured under static condition are known as base values.As a comparison, these parameters are also shown in Figure 7 as the profiles indicted by  = 0.0 Hz.As can be seen from the figure, the void fraction and interfacial area concentration are increased near the outer wall while decreased near the inner wall compared to the base values.This is consistent with the foregoing analysis and the phenomenon is enhanced with the increase of vibration frequency.When  = 0.5 Hz, the vibration acceleration of the test section is within 0.014 and the vibration velocity is less than 0.03 m/s.The additional bubble force due to inertia is rather small under this vibration condition.Therefore, the void fraction profile is almost the same as that in the stationary experiment.As the vibration frequency increases, the liquid inertia effect becomes significant and the additional bubble force increases greatly.As a result, radial bubble motions with bubbles moving toward outer wall become prominent as the test section moves from  = 0 to  =  1 .When the vibration frequency increases to 2.473 Hz, the peak near the inner wall is weakened (⟨  ⟩ = 1.03 m/s and ⟨  ⟩ = 0.195 m/s) or even disappearing (⟨  ⟩ = 0.52 m/s and ⟨  ⟩ = 0.195 m/s) and a strong peak skewed to outer wall is formed in the void fraction profile.From Figure 7 it can also be seen that the vibration effect on bubble Sauter mean diameter is not marked as compared with void fraction and interfacial area concentration.Therefore, bubble coalescence or breakup caused by structure vibration rarely occurs in bubbly two-phase flow.
where  -NonD ,  ,-NonD , and  sm,-NonD are the nondimensional peak values with respect to void fraction, interfacial area concentration, and Sauter mean diameter.The subscripts  and 0 denote the local peak values along the radius under vibration and static conditions, respectively.It can be seen from Figure 8 that the nondimensional peak void fraction and interfacial area concentration increase significantly with vibration frequency.Compared to the peak values measured in the stationary experiment, the peak void fraction can increase by 49% as the vibration frequency  increases to 2.473 Hz.But for the flow condition of ⟨  ⟩ = 1.03 m/s and ⟨  ⟩ = 0.403 m/s, the nondimensional peak void fraction is much smaller.This may be attributed to strong mixing of two-phase flow due to bubble-induced turbulence.The nondimensional peak Sauter mean diameter is within ±5%, which means that structure vibration does not change bubble size greatly in the present experiment.The average variation of local flow parameters throughout the 11 radial measurement locations is evaluated in this part through comparing the local parameters at vibration displacement of  =  1 with base values in order to further reveal the vibration frequency effect.Similar to MRV, the mean relative variation at  =  1 (MRV  ) is defined for evaluation.The MRV  has the same equations as MRV (i.e., (1)) except that the subscript  is replaced by subscript - in the expressions of MRV  .Here, the subscript - represents the local time-averaged values measured at the vibration displacement of  =  1 with vibration frequency of .The change of MRV  with respect to vibration frequency is shown in Figure 9.As the vibration frequency  increases from 0.472 Hz to 2.473 Hz, the average change of local void fraction can reach more than 29% for all the flow conditions under consideration.The change in local void fraction is most obvious for low liquid velocity of ⟨  ⟩ = 0.52 m/s with MRV  of void fraction reaching up to 41% when  = 2.473 Hz.The MRV  of Sauter mean diameter also increases slightly with , but the variation range is within 7%.

Conclusions
An experimental database has been collected in an annulus using double-sensor conductivity probe under vibration condition.The vibration amplitude of the test section is 9.5 mm and the vibration frequency ranges from 0.472 Hz to 2.473 Hz.Local measurements have been performed along one radius in the vibration direction.The local time-averaged flow parameters including void fraction, interfacial area concentration, and Sauter mean diameter measured in the stationary experiment are used as base values in the analysis of vibration effect.Based on the experimental data, the vibration effect especially the vibration frequency effect on local parameters has been discussed and the following results have been obtained.
The vibration effect on local time-averaged flow parameters increases slightly with vibration frequency but the overall effect is very limited.Compared to the base values, the mean variation of local time-averaged flow parameters is within 8% and the vibration effect is further weakened as the void fraction increases.Local flow parameters fluctuate continuously around the respective base values in the vibration cycle and the fluctuation amplitude increases with vibration frequency.The absolute change of void fraction near the channel wall can reach 17% during the vibration cycle under vibration frequency of  = 2.473 Hz.Additional bubble force due to inertia is found to play an important role in lateral bubble motions and the coaction of bubble forces determines the radial distribution of local parameters under vibration condition.With the increase of vibration frequency, the peak near the inner wall is weakened or even disappearing, and a strong peak skewed to outer wall is gradually observed in the void fraction and interfacial area concentration profiles at the vibration displacement of  =  1 .The nondimensional peak void fraction and interfacial area concentration increase significantly with vibration frequency and  -NonD can reach up to 49% at  = 2.473 Hz.The average variations of local void fraction and interfacial area concentration along the radius also increase with vibration frequency.MRV  of void fraction and interfacial area concentration can both increase to more than 25% for the flow conditions under consideration.But the increase of vibration frequency does not change bubble diameter significantly in the present work condition.
The current study shows that the vibration parameters have great influence on local parameters distribution under structure vibration and the vibration impact also varies with specific flow condition.Therefore, further study on vibration amplitude effect and the vibration effect on flow structure beyond bubbly flow is suggested.

Figure 3 :
Figure 3: Vibration acceleration and displacement profiles under different motor speeds.

Figure 4 :
Figure 4: Profiles of local time-averaged flow parameters under static condition.

Figure 5
shows the vibration effect on time-averaged flow parameters.For the flow condition with relatively low void fraction (i.e., ⟨  ⟩ = 1.03 m/s and ⟨  ⟩ = 0.195 m/s), the change in local time-averaged void fraction due to vibration is bigger than other flow conditions and the MRV of void fraction increases with vibration frequency .This phenomenon becomes insignificant with the increase of void fraction.The MRV of interfacial area concentration shows similar profiles.

Figure 5 :
Figure 5: Vibration frequency effect on MRV of time-averaged flow parameters.

Figure 6 :Figure 7 :
Figure 6: Variation of local flow parameters in the vibration cycle.

Figure 8 :
Figure 8: Vibration frequency effect on peak values of local flow parameters at  =  1 .

Figure 9 :
Figure 9: Vibration frequency effect on MRV  of local flow parameters.

Table 2 :
Vibration parameters under five motor speeds.The increase in ⟨  ⟩ also augments the void fraction at the peak and makes the void fraction peak sharp.For the flow condition of ⟨  ⟩ = 1.03 m/s, bubble coalescence occurs as ⟨  ⟩ increases to 0.403 m/s.The radial migration of bubbles plays an important role in the evolution of local flow structure.
fluctuation amplitude of the vibration acceleration increases with the motor speed and the vibration acceleration can reach up to 0.264 as the motor speed increases to 150 rpm.Table2summarizes the main vibration parameters under five motor speeds.The subscripts 1 and 2 represent the maximum values in the positive and negative vibration directions, respectively.The vibration of the test section was driven by eccentric cam and damping springs.The force exerted by damping springs may change its direction during vibration cycle, resulting in the difference between | 1 | and | 2 |.3.2.Local Time-Averaged Flow Parameters.Local two-phase flow measurements using double-sensor conductivity probe were first performed under static condition.Figure4shows the radial distributions of local time-averaged void fraction (), interfacial area concentration (  ), and Sauter mean diameter ( sm ).For the flow condition of ⟨  ⟩ = 0.52 m/s and ⟨  ⟩ = 0.195 m/s, intermediate peak characterized as broad peak in void fraction near the channel wall and plateau with medium void fraction around the channel center is observed.When ⟨  ⟩ is increased to 1.03 m/s, the radial position of the void fraction peak is moved toward the channel wall.tration and void fraction are very similar if the Sauter mean diameter is radially uniform.For ⟨  ⟩ = 0.195 m/s, the Sauter mean diameter profiles are almost uniform along the radius with some decrease near the channel wall.This yields a greater peak-to-average ratio in the interfacial area concentration profiles.But for ⟨  ⟩ = 1.03 m/s and ⟨  ⟩ = 0.403 m/s, the