Effect of Entrainment on the Liquid Film Behavior in Pipe Elbows

: Multiphase flow entrainment in natural gas engineering significantly influences the safety and efficiency of oil companies since it affects both the flow and the heat transfer process, but its mechanisms are not fully understood. Additionally, current computational fluid dynamics (CFD) methodologies seldom consider entrainment behavioral changes in pipe elbows. In this article, a verified CFD method is used to study the entrainment behavior, mechanism, and changes in an elbow. The results show that droplet diameter in a developed annular flow follows a negative skewness distribution; as the radial distance (from the wall) increases, the fluctuation in the droplets becomes stronger, and the velocity difference between the gas and the droplets increases linearly. Turbulence bursts and vortices sucking near the wall jointly contribute to droplet entrainment. As the annular flow enters the elbow, the secondary flow promotes the film expansion to the upper and lower parts of the pipe. Droplets re-occur near the elbow exit intrados, and their size is much smaller than those in the upstream pipe. Vortices sucking under low gas velocity play an important role in this process. These findings provide guidelines for safety and flow assurance issues in natural gas production and transportation and bridge the gap between multiphase flow theory and natural gas engineering.


Introduction
Currently, gas-dominant multiphase flow (i.e., churn flow and annular flow) has many natural gas engineering applications.Typically, in annular flow, high-speed gas moves through the pipe center while a liquid film flows around the pipe wall [1][2][3].In this process, the strong interaction between the gas core and liquid film shears part of the liquid rollwave crests and leads to droplet entrainment [4].Entrainment is critical in gas-dominant multiphase flow since it affects the mass and heat exchanges between phases as well as phase distribution and pressure drop.These changes destabilize the flow and make it hard to predict, which negatively affects the safety and flow assurance in natural gas production and transportation.Droplet entrainment characteristics are also important for developing an accurate multiphase flow model [5,6].However, the entrainment mechanism is not fully understood, and the empirical models available are based on simplified assumptions in specific conditions that are not accurate in many cases.
Some researchers have argued that the interfacial wave instability between two phases affects droplet entrainment [7][8][9].Depending on the liquid flow rate, two types of waves exist on the liquid film: ripple waves and disturbance waves.When the liquid flow rate is low, ripple waves (small-scale) occur, which exist for a short time and move slowly.Therefore, droplets cannot be extracted from the liquid film.For high gas and liquid flow rates, the disturbance waves (long-length) become dominant.Since disturbance waves move faster than ripple waves, their amplitude is much higher, and their lifespan is longer [10,11].Usually, for pipelines with a small diameter, the disturbance waves occur coherently on the circumference [12].According to the literature, high-speed gas shears off the disturbance wave crest to form droplets; then, the droplets enter the gas core [13,14].In this process, the gas drag force, as well as gravity, deforms and stretches the crest against the retaining liquid surface tension.As a result, the liquid ligament breaks up, and the droplet enters the gas core.However, the liquid volume detached from the wave crest is not taken into consideration in this mechanism, which makes it impossible to calculate the entrainment rate.
There are many studies on liquid sampling, flow patterns, and factors that influence liquid distribution as well as erosion using both experimental and simulation methods.Some of these are summarized in Table 1.Regarding the empirical modeling, for vertical flow, Oliemans et al. [25] developed a correlation based on their experimental database that covers a wide range of flow conditions.Later, Ishii and Mishima [26] derived a verified correlation based on the gas Weber number and the liquid Reynolds number.Researchers also studied droplet sizes entrained in the gas core and proposed correlations based on their experimental database.Most of the correlations [27-29] considered the Weber number (the ratio between inertia and surface forces) as a key parameter since it indicates the interaction between the external stress force and the surface forces acting on droplets, an interaction which is important for the droplet breakup mechanism and droplet size.Kolev [27] derived a correlation to calculate the droplet diameter from Kelvin-Helmholtz instability.Fore [30] put forward a modified expression to predict the volume median diameter in vertical and horizontal flows in a pipe.Ishii [31] proposed a semi-empirical correlation to estimate the average maximum droplet size for roll-wave entrainment.
Summarizing the entrainment models above, Wicks's model and Paleev's model are based on experimental data in horizontal annular flow, so they are unlikely to provide a precise prediction for vertical or inclined annular flow.In contrast, Oliemans' and Ishii's models can provide accurate entrainment in vertical annular flow (superficial gas velocities of 15 to 40 m/s, superficial liquid velocities of 0.06 to 0.2 m/s, and pipeline diameter of 6 to 32 mm) since they use a vertical annular flow database.Magrini's experiments [32] show that Paleev's model provides the best prediction for horizontal annular flow.For the entrained droplet size calculation, most of the correlations choose the Weber number as the key factor.Additionally, these semi-empirical correlations are based on experimental data in specific conditions, which limits their applications.
In brief, the entrainment mechanism in annular flow is not fully understood; related models and correlations are empirical and inaccurate, which limits their application in practice.Many details on droplet entrainment are not included in the current CFD methodologies available in the literature, and no studies on how entrainment affects liquid film behavior in the elbow were found.
In this article, a verified CFD method is used to study entrainment behavior in pipelines.Firstly, entrainment characteristics like droplet diameter and velocity distribution are analyzed in detail.Then, the entrainment mechanism is interpreted based on a vortex theory.Lastly, entrainment changes in elbows are studied.These detailed analyses of liquid film and droplet entrainment help to evaluate elbow safety, with elbows being sensitive to erosion as well as corrosion, and solve flow assurance problems in natural gas production and transportation.

Entrainment Modeling
The correlation established by Wicks and Dukler [15] is shown below: where R is the entrainment group (unit ft 3 /lbf-hr); q L and q G are the mass flow rates of liquid and gas, respectively; N Wecr is the critical Weber number ranging from 13 to 22; W LE is the droplet mass flow rate; and dp/dL is the pressure gradient.
The correlation developed by Paleev [16] is listed below: where ρ C is the mixture density of the core defined as where F E is the entrainment fraction; ρ L and ρ G are the density of the liquid and gas, respectively; µ L is the dynamic viscosity of the liquid; v SL and v SG are the superficial liquid velocity and superficial gas velocity, respectively; and σ is the surface tension.Oliemans' model [25] is given below: Ishii's model [26] is represented below: We SG N 0.25 where

Droplet Size
The critical Weber number is the most widely used criterion to estimate the droplet size; for the condition of high density difference (droplet in gas core), the empirical value of the Weber number is 12.The expression is as follows [27]: where We crit is the critical Weber number, u g is the gas velocity, and ∅ d is the droplet diameter.
For annular flow, the volume median diameter of the droplet derived by Fore [30] can be calculated using the following equation: where ∅ vm is the droplet volume median diameter, and We g and Re g are the gas Weber number and the gas Reynolds number, respectively.
This correlation is simple and accurate enough to predict the droplet volume median diameter in horizontal and vertical pipelines.
The semi-empirical correlation to estimate the average maximum droplet size proposed by Ishii [26] is as follows: where ∅ max is the average maximum droplet size.

Experimental Background
The experimental data used for liquid entrainment verification in this study were taken from Magrini [32].In the experiment, air from the compressor and water from the tank were measured separately and then mixed to develop a horizontal annular flow in a 76.2 mm ID pipeline.Then, the fluid entered the test section after it was fully developed.The test section included two parts: the isokinetic sampling section and the film removal section.The isokinetic sampling section measured the liquid entrainment in the gas, and the film removal section determined the liquid film flow rate.Finally, the fluid flowed out of the system from the outlet.
The isokinetic sampling section included an L-shape isokinetic probe, two valves, a container, and a gas flowmeter.The probe inlet was placed against the gas flow direction so that the droplets in the gas core could enter the probe easily.The probe could be located at different positions of the pipe cross-section.Valves were applied to control the sampling rate and keep the gas velocity in both the probe and the pipeline the same.The sampled droplets were collected in the container, and the sampled gas was vented out through the flowmeter.Then, the container was weighed to calculate the entrainment.
The film removal section includes an inserted porous sleeve that allows the liquid film to pass through, and the cavity between the sleeve and outer shell is used to store the sampled liquid.When passing the film removal section, liquid film passes through the porous section and enters the cavity.The entrained droplets in the pipe center move out directly without being removed.The collected liquid is used to calculate the film flow rate.
In these experiments, the superficial gas velocity ranged from 40 m/s to 80 m/s, and the superficial liquid velocity ranged from 0.0035 m/s to 0.04 m/s.Each test lasted for 5 min.The entrainment flux at each location can be calculated as The entrainment fraction can be calculated as where E X is the entrainment flux, V E is the collected liquid volume, t S is the sample time, and A is the cross-sectional area.

CFD Simulation
It is difficult to acquire details of the entrainment process since it is random and changes fast.Fortunately, CFD can provide abundant information about the whole domain at any specific time.For liquid entrainment simulation, a volume of fluid (VOF) model based on the Euler method is adopted, in which the interface between immiscible fluids can be obtained.The main steps in the liquid entrainment simulation are as follows: (1) solve the single set of momentum equations shared by all the phases; (2) track the volume fraction in each computational cell through the domain; (3) solve the volume fraction equation through implicit time discretization.A flow chart representing the methodology is shown in Figure 1.

Geometry
The geometry used in this investigation is shown in Figure 2. The upstream length, L, is 10D, the curvature radius of the elbow is 1.5D, and the downstream length is 3D.

Meshing
A high-quality mesh is critical for a successful simulation.Since the entrainment process occurs at the interface of the two phases near the wall, elements in this region were refined to obtain the details. Figure 3 provides the overview and cross-section details of the mesh.

Meshing
A high-quality mesh is critical for a successful simulation.Since the entrainment process occurs at the interface of the two phases near the wall, elements in this region were refined to obtain the details. Figure 3 provides the overview and cross-section details of the mesh.It can be seen that the domain is meshed by polygonal cells.Polygonal cells can reduce the skewness, improve the parallel solver load balancing, and minimize communication by decreasing the number of partition interface cells [33]. Figure 3c shows an enlarged view of the near-wall region, where the first-layer thickness of the mesh is 20 µm and the maximum cell size is 1 mm, the growth ratio is 1.2, and the row number of the boundary is 5.The length of the cells in both circumferential and axial directions is 1 mm, and the total number of elements is 11,338,564.Since the maximum cell size is in the same order as the droplet size, this mesh can reveal detailed information during the entrainment process.

Simulation Settings and Boundary Conditions
In this study, the annular flow simulation is performed first to analyze the entrainment process in the pipeline.The settings of the simulation are shown in Table 2.The boundary conditions of the flow simulation are shown in Table 2. Water (liquid phase) enters the pipe circumferentially along the pipe wall, and the water entrance is ring-shaped.Air (gas phase) enters the pipe in the center of the cross-section, and the shape of the air entrance is a circle.The superficial gas velocity is 24.6 m/s, and the superficial liquid velocity is 0.1 m/s; the initial values of the corresponding air and water mass flow rates are listed in Table 3, and the time step is 0.000005 s.Once the mass flow rate difference between the inlet and outlet decreases below 0.01 kg/s and the flow state in the pipe becomes steady, the simulation can be stopped.

Entrainment Fraction Validation
In this section, experimental data from the literature [32] is adopted to verify the entrainment fraction in a CFD simulation.The working conditions in both experiments and simulations are the same.The superficial gas and liquid velocities are 40 m/s and 0.01 m/s, respectively.The entrainment fractions measured by isokinetic sampling and film removal methods in a vertical orientation are 0.42 and 0.51, respectively.
Based on Equations ( 18)-( 20), the entrainment fractions calculated by Oliemans' and Ishii's models are 0.53 and 0.55, respectively.As for the simulation, at the position 1D before the elbow entrance, the area-weighted average void fraction of the gas core and pipe cross-section are extracted and converted to the average entrainment fraction using the following equations: ) where A d and A f are the area of droplets and liquid film in the cross-section, respectively; A c is the area of the gas core; h gc and h P are the liquid holdup of the gas core and pipe cross-section, respectively; and F ECFD is the entrainment fraction calculated using CFD.
The entrainment fractions calculated above use Equations ( 4)-( 7) and are summarized in Table 4.As shown in Table 4, the entrainment fraction measured using isokinetic sampling is smaller than that of the film removal method.Droplets adhering to the small diameter tube surface may cause this difference.As for the film removal method, the liquid collected may contain the redeposited droplets, which makes the measured value larger.Thus, the real value should be between 0.42 and 0.51.Both Oliemans' model and Ishii's model give larger values than the experiments, indicating that the two models are not accurate enough.The entrainment fraction predicted using CFD is 0.47, just within the range (0.42 to 0.51), so the entrainment fraction predicted with CFD is reasonable.

Entrainment Behavior in the Pipeline
After the flow simulation stops, the phase contours as well as the velocity vectors of the pipeline longitudinal section and representative cross-sections are extracted to analyze the entrainment behavior in the pipeline.

Entrainment Distribution
Figure 4 shows the phase distribution and velocity vector in the pipeline.It can be seen that as the upstream length increases, the gas in the pipe center moves much faster than that near the wall, the gas velocity fluctuation in the pipe becomes obvious, and the uniform liquid film starts waving and then breaks up into droplets (see Figure 4b).Then, the droplets are gradually entrained into the pipe center by the high-speed gas.When the upstream length reaches 6D, the droplet distribution seldom changes, and the entrainment process reaches dynamic equilibrium.After entering the elbow, most of the droplets hit the elbow extrados and coalesce into liquid film.Then, the film flows along the extrados side to the downstream pipe (see Figure 4c).Representative cross-sections of the upstream pipe are extracted to reveal more information (see Figure 4d).At an upstream length of 2D, there is a film wave at the right side of the cross-section (circled).The corresponding velocity vectors indicate that a turbulence burst here contributes to the wave.When the upstream length reaches 2.4D, the film becomes discontinuous, and "liquid bridges" (circled) occur, which are also related to the turbulence bursts near the wall.As the flow develops, turbulence bursts become stronger and break the liquid bridges into droplets, and then, the droplets are entrained into the pipe center by the gas.At an upstream length beyond 6D, the droplets are distributed uniformly in the pipe cross-section.Droplet coalescence (circled) is also found in the flow at an upstream length of 7D.
Other researchers like Zahedi [34] and Farokhipour [21] also studied liquid distribution in annular flow pipelines.Zahedi [34] only reported the liquid distribution in the middle of the elbow (elbow angle = 45 • ), and no liquid was found at the elbow extrados; however, this seems inaccurate, because droplets carried by the gas core will hit the extrados and form a liquid film.Farokhipour's simulation [21] showed the liquid film distribution in the elbow well, but no droplets can be observed in the figures.Compared with their studies, the simulation in this study provides more details of droplet distribution, like droplet coalescence, which can be seen as a remedy for entrainment simulation.Another difference is the grid.Other researchers use hexahedral cells in their simulations, but polygonal cells are adopted in this study to reduce skewness to obtain more accurate results.Additionally, droplet distribution and entrainment changes in the elbow are also investigated, which have not featured in other simulations.

Entrainment Characteristics
Further, the void fraction (area-weighted average) of the cross-sections (every 0.25D) in the upstream section is extracted and summarized in Figure 5.It can be seen that for the upstream length of 2D, the void fraction is steady but fluctuates slightly around 0.95.As the upstream length increases to 6D, the void fraction rises to 0.978.Meanwhile, its fluctuation becomes stronger, and droplets start to take shape and enter the gas core.Then, after 6D, the void fraction becomes steady again, with its value fluctuating more vigorously between 0.978 and 0.981.It seems that for the annular liquid inlet, the liquid droplet entrainment and re-deposition process reach dynamic equilibrium when the upstream length is beyond 6D.
According to Ishii and Kataoka [5,35], for the case where the liquid is injected smoothly as a film at the inlet, the entrained fraction expression takes the form of exponential relaxation.For the case in this study, the trend of the curve in Figure 5 is also exponential relaxation, which means the CFD results here are reasonable.
Then, the droplet information (droplet number, droplet diameter, droplet velocity, and droplet position (radial distance from the wall)) in the pipeline's longitudinal section is collected to further analyze the entrainment.Droplet diameters and positions are measured by the scale of the phase distribution contour; droplet velocity is determined by the velocity vector in the center of the droplet.
Figure 6 shows a histogram of the relationship between droplet diameter and droplet number distribution in the longitudinal section in the developed annular flow.It can be seen that the droplet number roughly follows a negative skew distribution.Most droplet diameters are distributed in a 1 mm to 2.5 mm range, and the droplet number in the range of 2 mm to 2.5 mm is the largest, while in the range of 0 to 0.5 mm, it is the smallest.In contrast to one study [36], for the conditions of low-pressure (below 3 atm) and low-viscosity fluid (water), the droplet diameter distribution in this study is similar to the upper limit log-normal (ULLN) distribution given by Mugele and Evans [37].The maximum size of the droplet can be calculated based on the Kelvin-Helmholtz instability and critical Weber number.
where ∅ dmax is the maximum diameter of the droplet, and We crit is the critical Weber number, which equals 12.For this study, the maximum diameter of the droplet calculated using this equation is 2.7 mm, which is represented by the last bar in Figure 6.For this study, the calculated volume median diameter of the droplet is 1.1 mm, while the volume median droplet diameter predicted using CFD simulation is 1.7 mm (54.5% higher than the calculated value).Further improvement may be achieved by refining the mesh, which is limited by the current computational resources.
Figure 7 shows the relationship between the radial distance (from the wall) and the droplet diameter in a developed annular flow.As shown in this figure, the average value of all the droplet diameters is 1.81 mm.Droplet diameter fluctuates around the average value, and its amplitude increases with the radial distance.Table 5 presents the variance distribution of droplet diameters.It can be seen that as the radial distance (from the wall) increases, the expectation of the droplet diameter varies very slightly around 1.81 mm, by less than ±0.05 mm.However, for the variance, it increases from 0.19 to 0.63.This means that the change in droplet diameter becomes increasingly unstable the closer the droplet is to the gas core.This may be because droplets closer to the pipe center are more likely to be broken up or coalesced by the gas core turbulence, which strengthens the fluctuation.According to Kolev [27], the droplet breakup mechanism can be expressed as a balance between the external stress force and the surface force.During the breakup, external stress force tries to disrupt the droplets, while surface tension force tries to avoid droplet deformation.From the critical Weber number, Equation ( 11), it can be seen that the droplet diameter varies with the inverse square of the gas velocity.As the radial distance from the wall increases, the gas velocity becomes higher, so the minimum value of droplet diameter decreases.On the other hand, according to the Kelvin-Helmholtz instability, as the gas velocity increases, the gas-liquid interface becomes more unstable.Additionally, droplet diffusion becomes stronger, so droplets are more likely to coalesce, which increases the maximum droplet diameter.
Figure 8 shows velocity distribution (both gas and droplet) according to radial distance (from the wall).It can be seen that both gas and droplet velocities rise and fluctuate as radial distance increases, while the velocity difference between gas and droplets also increases, which means that the slippage becomes more obvious.Moreover, the fluctuations in both gas and droplets become stronger.In order to clearly express the droplet slippage, the droplet slippage ratio can be calculated with the following equation.
where S d is the droplet slippage ratio, and u d is the droplet velocity.It can be found that the slippage ratio near the wall (27.3%) is smaller than that in the pipeline center (33.8%).

Entrainment Mechanism
In this subsection, the mechanism of the entrainment process in annular flow is interpreted gradually, as shown in Figure 9.In the initial state (see Figure 9a), the liquid film and gas are uniform and there is no fluctuation.Since there is a velocity difference between gas and liquid film, vortices occur at the two-phase interface under the effect of shear and friction, which produce waves in the film and spread downward (see Figure 9b).As the flow develops, the crest stretched by the vortices becomes higher, and its bottom is necking under the effect of surface tension (see Figure 9c,d).Then, the necking breaks, the liquid entering the gas becomes droplets, and the rest of the liquid becomes film (see Figure 9e).The high-speed gas moving in the pipe center lowers the dynamic pressure here, so the droplets can be entrained under the pressure difference (see Figure 9f).This process reduces the film volume and makes the film discontinuous.The Kelvin-Helmholtz instability is mainly associated with flows that have tangential variation in the velocity field [38][39][40].This instability is caused by the hydrodynamic amplification of perturbations that arise at the gas-liquid interface with a discontinuity in the velocity field.The entrained droplet size is approximately equal to the height of the most unstable wavelength.
where ∆h w,K−H is the height of the wavelength.For this study, the calculated droplet diameter is 1.8 mm, which is close to the average droplet diameter in Figure 6.This indicates that the simulation results in this study fit the Kelvin-Helmholtz instability theory and are reasonable.
As the flow develops downstream, the entrained droplets travel to the wall and are deposited into the liquid film again.As the entrainment and deposition reach dynamic equilibrium, the flow is fully developed.The most important factor that affects droplet entrainment is the vortices near the gas-liquid interface.

Entrainment Changes in the Elbow
The entrainment changes when the fully developed flow enters the elbow.Figure 10 shows the phase distribution and velocity vectors of the flow in the elbow: the flow is restricted by the geometry, and the gas moves toward the elbow intrados as the flow first enters the elbow.Droplets are uniformly distributed in the pipe cross-section (elbow angle = 0 • ).When the flow reaches an elbow angle of 30 • , more droplets are deposited on the extrados, and the liquid film starts to take shape.The flow near the wall moves toward the intrados and forms the secondary flow.For an elbow angle of 60 • , most droplets are deposited on the extrados, and droplets in the gas become less frequent, while the film becomes thicker.Two counter-rotated vortices drive the film toward the extrados.For an elbow angle of 90 • , the film at the extrados thickens.Droplets occur again at the intrados under the strong effect of the vortices.After entering the downstream pipe, the film thickness at the extrados increases, and the vortices become weaker.
Unlike the roll-wave mechanism proposed by Ishii [26], which is reasonable for high gas velocity in a straight pipe, droplets re-entrained at the intrados seem to be caused by the mechanism of vortices sucking, since the gas velocity here is relatively low.As shown in Figure 4c, it is also found that the re-entrained droplet size at the intrados is smaller (1 mm-1.5 mm).This may be because vortices at the intrados provide more sucking and external stress force on the droplets for their breakup.

Conclusions
In this study, a verified CFD method is used to investigate entrainment behavior in a pipeline.Firstly, entrainment characteristics like droplet diameter and velocity distribution are analyzed in detail.Then, the entrainment mechanism is interpreted based on a vortex theory.Lastly, entrainment changes in the elbow are studied.The results are valuable and provide guidance for entrainment study in multiphase flow.The main conclusions are summarized below: (1) Under the annular inlet condition, the entrainment distribution along the upstream developing length from the entrance is close to a form of exponential relaxation.The droplet diameter in a developed annular flow follows a negative skewness distribution, with most droplet diameters ranging from 1.5 mm to 2 mm.As the radial distance (from the wall) increases, the fluctuation of the droplets becomes stronger, and the velocity difference between gas and droplets increases linearly.
(2) Turbulence bursts and vortices near the wall jointly deform the liquid film and crest so that they begin to take shape; then, the external stress force acting on the liquid crest overcomes the surface force, and droplets form.Finally, vortices sucking transforms the droplets into gas, and the entrained droplet size is approximately equal to the height of the most unstable wavelength.(3) As the annular flow enters the elbow, droplets hit the extrados and form a liquid film, the droplet number in gas decreases, and at an elbow angle of 60 • , there are no droplets in the pipe.Secondary flow promotes the film formation and expansion of the film to the upper and lower parts of the pipe.However, for an elbow angle of 90 • , droplets re-occur near the elbow intrados, and the re-entrained droplet size is much smaller than that in the straight upstream pipe.Vortices sucking at low gas velocity play an important role in this process.

Unique work description
In this study, a more accurate simulation using polygonal cells instead of hexahedral cells was carried out to investigate details of droplet distribution and entrainment changes in the elbow, which is unique and different from other simulations.

Recommendations and future work
Although polygonal cells can provide more accurate results, they require more time and more computational resources, especially for finer grids.Grid optimization is a promising topic.In the future, AI technology like cloud computing and large data can be integrated into CFD simulations to enhance efficiency and consume fewer resources.Combined with pipeline-monitoring technologies, CFD simulation can make the detection system's response faster and smarter.
Author Contributions: Z.X.: conceptualization, formal analysis, investigation, writing-original draft.X.C.: conceptualization, supervision.All authors have read and agreed to the published version of the manuscript.

Figure 1 .
Figure 1.Flow chart for the methodology.

Figure 2 .
Figure 2. Geometry configuration of the computational domain.

Figure 2 .
Figure 2. Geometry configuration of the computational domain.

Figure 4 .
Figure 4. Phase distribution and velocity vector of the pipeline: (a) overview; (b) enlarged view of the upstream; (c) enlarged view of elbow and downstream; (d) representative cross-sections.

Figure 5 .
Figure 5. Relationship between upstream length and area-weighted average void fraction of the pipe cross-sections.

Figure 6 .
Figure 6.Relationship between droplet diameter and droplet number distribution.

Figure 7 .
Figure 7. Relationship between radial distance (from the wall) and droplet diameter.

Figure 8 .
Figure 8. Velocity distribution according to radial distance (from the wall).

Figure 9 .
Figure 9. Mechanism of the entrainment process in annular flow.

Figure 10 .
Figure 10.Phases and velocity distribution for different elbow angles.

Table 2 .
Settings of the flow simulation.

Table 3 .
Boundary conditions of the flow simulation.

Table 4 .
Entrainment fraction obtained using different methods.

Table 5 .
Variance distribution of droplet diameters.