1. Introduction
Gas-liquid two-phase flow is a very common type of flow widely used for energy in the chemical and power industry, as well as other fields [
1,
2,
3,
4,
5]. Pneumatic conveying is a classic case in which gas carries discrete particles or another fluid downstream of the pipeline using interfacial forces. From the perspective of energy transport, this process is the gas phase that acts as the power source, working directly on another phase in order to realize energy conversion and transfer. Taking ejector pressurization as an example, the high-speed air flow formed by high-pressure steam is the power source in the gas-water two-phase flow injector [
6]. The gas phase directly expands the water in a variable section pipeline where both are accelerated to form a supersonic saturated gas-water mixture. After a condensation shock, a sudden pressure increment occurs, leading to considerable pressure gain in comparison to the original steam. The low water-head hydropower system proposed by Sun et al. [
4] is another case involving interphase mechanical energy transfer where the dropping water directly pressurizes the natural air in the gas-water energy conversion equipment and then the high-pressure air pumps the water upward through the interfacial forces in the high-pressure pump.
When a certain mass flow of high-pressure air is introduced into the water nozzle during water jet propulsion, the liquid phase is also accelerated, while the air flow expands, gaining additional kinetic energy increments and introducing thrust augmentation, or efficiency enhancement [
5,
7,
8,
9,
10,
11,
12]. In contrast to turbomachinery, such as gas turbines, water pumps, etc., the gas in the two-phase flow directly pushes water to form a powerful jet. Today, improving efficiency and reducing energy loss are some of the most important issues in engineering. To achieve the optimal use of energy, much work has been performed on the optimization of energy systems, such as optimizing the gas microturbine cycle by bee algorithm [
13], adopting phase change materials in freezers [
14] and optimizing indirect flat evaporative coolers using exergy analysis [
15]. Therefore, before we optimize the use of gas energy in two-phase flows, it is important to understand the approaches through which the gas transfers its potential pressure energy to the liquid’s kinetic energy in variable cross-sectional tubes.
Gas-liquid two-phase flow patterns in horizontal or vertical pipelines can be divided into stratified, wavy, annular, bubbly and mist flows, etc. [
2,
16]. The gas-liquid interfacial area and spatial distribution vary considerably for different flow patterns; the coupling mechanism and coupling strength between two phases are also different. Since the work of gas to liquid is accomplished through momentum coupling along the gas-liquid interface, the delivery capability of gas power to liquid must change with the flow pattern [
17].
In separated flows, such as stratified, wavy and annular flows, the two phases are segregated by a continuous interface. The physical properties and velocities of the two phases are different, so the mathematical model consists of hydrodynamic equations for every single phase, as well as the appropriate kinematic and dynamic relations at the interface. Free streamline theory is one of the successful examples of this processing strategy, although its interface relation is very simple [
18]. The one-dimensional two fluid model is widely used in engineering applications to reduce the computational difficulties caused by multi-dimensional effects and to simultaneously capture the main characteristics of the flow accurately [
19]. One-dimensional models focus on the average quantities in which the time, volume and ensemble average are adopted, disregarding details of the flow field, and making it necessary to introduce additional relationships to close the equations [
20,
21,
22].
Fontalvo et al. [
17] studied the influence of different interface friction formulations, and the formation and propagation of interfacial waves by a two-fluid model for vertical annular flow, where the interfacial fraction factor derived by Wallis et al. was introduced as the closure condition. Castello-Branco et al. [
19] analyzed the stability of vertical annular flow using a transient one-dimensional two fluid model where the interfacial shear stress was determined by the correlation of Whalley and Hewitt. They found that the closure relationship notably influences the calculated wave frequencies and growth rates. Emamzadeh [
1] accurately predicted the transition from the stratified regime to the annular regime in a horizontal pipe flow using his two-fluid model in which the interphase relation derived by Ottens is selected. By considering droplet entrainment and gas condensation, Stevanovic and Studovic [
23] developed a simple one dimensional model to describe steady vertical annular flow, or the horizontal stratified flow of gas, liquid and liquid droplets, using Wallis’ correlation to determine the interphase shear stress. Utilizing a similar model, Sugawara explored the deposition and entrainment of droplets in annular flow [
24].
Using bubbly flow as an example in a dispersed flow, a number of bubbles with different physical properties and slip velocities are scattered in liquid, which can be presented by a trajectory model or a two fluid model. Wang and Chen [
25] achieved an accurate prediction of the relative velocities between phases and good agreement with the subsonic flow experiments by means of their bubble model for a convergent-divergent pipe flow in which drag, virtual mass and pressure gradient forces were introduced and the modified Rayleigh-Plesset equation is adopted to describe the dynamic property of the bubble. Gowing et al. [
8] applied a similar model to simulate a steady one-dimensional bubbly flow in three different convergent nozzles; their predicted thrust and efficiency were nearly identical to the experimental results. Fu et al. [
26] analyzed bubbly flow in convergent-divergent nozzles using a two fluid model and discussed the influence of mass flow rate, nozzle area ratio and bubble radius on the jet thrust. Wu et al. [
9] developed a trajectory model for bubbly flow in which the drag, virtual mass, pressure gradient, gravity and Saffman forces for bubbles are taken into account. A modified Rayleigh-Plesset equation and the Keller-Herring equation are used to express the bubble dynamics. Three-dimensional numerical simulation of the flow in a convergent-divergent nozzle showed good agreement with the experiments, validating the concept of thrust augmentation by air injection.
For some specific flow regimes, such as cavitation flow, flash flow [
27] and bubbly flow with very small bubbles, a homogeneous two-phase flow model with acceptable accuracy and greater convenience is suitable where the relative velocity is negligible, which is considered to be the limit of gas-liquid flow with infinite velocity relaxation rates and maximum momentum exchange between the two phases. LeMartelot et al. [
27] developed three homogeneous flow models, representing three limit states: interphase mechanical equilibrium, mechanical and thermal equilibrium, and thermodynamic equilibrium. Städtke [
21] also derived a homogeneous flow model from the two-fluid model.
Variable cross-section flow is attracting more interest, and new flow methods have been developed. Neil and Stuart [
28] measured the jet thrust of an annular flow nozzle with variable cross-section and compared it with a two-fluid model where the interfacial friction factor was developed by Ambrosini et al. [
29]. It was demonstrated that the predicted thrusts compared reasonably well with experiments. Noticeably, Pandey and Singh [
30] investigated the peristaltic transport of Herschel–Bulkley fluids in variable cross-section tubes by using a model where non-linear governing equations are linearized by a low Reynolds number and long wavelength approximations. For the pipe networks with complex geometric structures, Baranovskii [
31] proposed a novel network model to describe steady-state 3D flows by rejecting averaging of the velocity field and applying conjugation conditions that provide the mass balance for interior joints of the network. Based on this model, a feedback optimal control problem was also studied by Baranovskii [
32]. Zahedi and Babaee Rad [
33] studied the effects of 90-degree horizontal elbows with different curvature radii on air–water slug flow behaviors experimentally and numerically.
For the gas-liquid interaction of a two-phase flow, much attention has been directed to interphase forces and their mechanisms, or model expressions. Less attention has been given to the characteristics of momentum transfer and mechanical energy conversion between gas and liquid in the tube—especially in variable cross-section nozzles—which is important to the study and application of water jet propulsion and to making full use of gas energy. In this work, which focuses on flow in variable cross-section nozzles that generate thrust, the mechanical work performed on the liquid and the energy conversion efficiency during gas expansion and flow acceleration were analyzed using a two-fluid model for annular and bubbly flows, and a homogeneous model with mechanical equilibrium, selected to represent the growing interphase coupling in different flow patterns. The interphase energy and work transfer were calculated quantitatively and compared systematically for three flow patterns, which were rarely involved before. These are crucial for the selection of flow patterns in nozzle design and can provide direction for optimization to achieve maximum efficiency gains.
2. Parameters of Energy Transfer and Conversion Capabilities
In the gas-liquid flow of a variable cross-section pipe, the gas phase acts as the energy source, or the fuel [
8], and it performs expansion directly on the fluid, resulting in an increase in fluid kinetic energy. To obtain a systematic description of the work transfer and energy conversion, four coefficients were introduced, as follows.
(a) Work factor of gas
In gas-liquid two-phase flows, the work factor of gas is introduced to measure the proportion of gas output performed by interfacial forces to the theoretically available work of the gas.
Taking the gas-liquid flow in the pipe as a steady flow system, the energy equation of gas is obtained by neglecting the gravitational potential energy and wall friction, which excludes deviations from different experiential wall friction formulas for various flow patterns.
where
is the mass flow rate of gas,
is the specific enthalpy of gas and
is the gas velocity. The gas produced shaft work to liquid through the gas-liquid interfaces at a rate of
. Equation (1) shows that
is just a definite part of the enthalpy drop of the gas and can be obtained directly.
However, Equation (1) does not provide an explanation of how the power is generated, making it necessary to calculate the work--or the power--by forces. Because it is a one-dimensional steady flow, the streamlines coincide with the path lines. Then the flow parameters can be specifically described by a spatial coordinate, or equivalently expressed by time coordinates. Thus, the overall work performed by the gas on the liquid can be obtained by
, where
is the interfacial force and
is the flowing time of a particular fluid particle in the system. The work factor of gas is defined as:
In Equation (2), the denominator is the effort of the gas to maintain this work output. This coefficient represents the capability of the gas to output work through interfacial forces.
(b) Interface work transfer coefficient
The interfacial force is a pair of interaction forces for gas and liquid. The velocities of gas and liquid are usually different, thus the mechanical power of
on the liquid can be expressed as
, where
is liquid velocity. The interface work transfer coefficient
is a ratio of the work received by the liquid to the output work of the gas performed by the gas-liquid interfacial force.
It is apparent that greater slip velocity results in a smaller interface work transfer coefficient, which means greater work loss.
(c) Interphase work-to-energy conversion coefficient
The interphase work-to-energy conversion coefficient is one of the most popular metrics in gas-liquid two-phase flow; it is defined as the ratio between the mechanical work received by the liquid and the theoretical work available to the gas, which demonstrates how much of the gas energy is used to accelerate the liquid. The expression for this coefficient is:
Apparently, .
(d) Interphase mechanical efficiency
Interphase mechanical efficiency is defined as the ratio of the total mechanical energy of the gas and liquid at the outlet of the pipe to the energy at the inlet considered not only from the perspective of the gas but from the overall perspective of the gas-liquid two-phase nozzle flow that generates the thrust.
where
is the specific internal energy of gas. This efficiency represents the total energy loss caused by the transfer of mechanical energy between phases.
3. Annular Flow in Variable Cross-Section Tubes
Based on the two fluid model proposed by Wallis [
18], a one-dimensional model was established for annular flow in a variable cross-section pipe. The current model disregarded the exchange of mass and heat between phases; only momentum exchange was included to study the mechanism and law of work and energy transfer between the phases.
3.1. Flow Model
Regardless of the mixing process in front of the pipeline, the annular flow was regarded as a steady flow with an unchanged flow pattern during the entire process.
The mass equations for gas and liquid are given by:
where
and
are the volume fraction of gas and liquid, respectively;
and
are the density of gas and liquid, respectively;
and
are the mass transfer to gas and liquid, respectively; and
is the cross-sectional area of the pipeline. Assuming that there is no phase transition, gas concentration or liquid evaporation in the annular flow, the source term in Equations (6) and (7) is zero:
.
The momentum equations for gas and liquid are given by:
where
is the shared static pressure of gas, liquid and interface, and
and
are the momentum transfer to gas and liquid, respectively. Under ideal wall assumptions,
is the momentum transfer from gas to liquid through interfacial forces.
is the velocity associated with the interfacial mass transfer.
The energy equations of gas and liquid are given by:
where
is the enthalpy associated with mass transfer. As a result of the fast flow velocity and short staying time in the tube, the heat transfer in the flow was insufficient, so heat exchanges between phases and external heat exchanges were all ignored, resulting in an adiabatic flow.
Because the pressure was not high, the liquid was considered to be an incompressible flow with a constant density. Gases, however, are compressible, and the ideal gas equation of state was chosen. To ensure the symmetry of the gas-liquid governing equations, the stiffened gas equation of state (SG EOS) was used to describe the state of the gas and the liquid:
where
and
are the specific heat ratio of gas and liquid, respectively, and
,
,
and
are constants determined by the properties of gas and liquid, respectively. In Equations (12) and (13),
indicates how stiffened it is compared to the ideal gas and
represents the internal energy offset. The SG EOS takes into account the compressibility of liquid and is widely used in compressible two phase flow models. For more details, see Yeom and Choi [
34].
The void fraction for gas and liquid reads:
All the variables in Equations (6)–(14)are related only to the position x of the tube axis, where the tube profile A = A(x) is a geometric parameter that is known when the tube profile is given. There are nine unknowns: , and nine ordinary differential equations. After the momentum transfer model and were established, the equations were closed and could be solved using the Runge-Kutta method.
3.2. Closure Condition
The momentum exchange between gas and liquid in annular flow relies on interfacial shear force, which is influenced mainly by interfacial wave properties, such as wave type, length, velocity and amplitude, as well as other characteristic parameters [
1,
35,
36,
37] and entrainment [
3]. However, due to the complex mechanisms of interfacial waves and entrainment, an analytical model for interfacial shear force is difficult to obtain; thus, empirical, or semi-empirical, relations--fitted by experiments--are widely used. In these relationships, the interfacial shear stress is obtained from the measured pressure drop based on the momentum balance equation of the gas core; its expression refers to the knowledge of wall friction in a single-phase tube flow, which can be expressed as [
3,
20,
37,
38,
39]:
where
is the interfacial friction factor obtained through experiment, which can be expressed by dimensionless parameters such as the Reynolds number of the gas and liquid, dimensionless liquid thickness, etc., after dimensional analysis. Wallis proposed one of the earliest models for vertical annular flow [
20]. Based on the Wallis model, many improvements were made to develop different
models suitable for flow conditions in a specific range [
3,
24,
35,
36,
40,
41,
42,
43].
To obtain a correlation applicable to a wider flow range, Aliyu et al. [
3] fitted the interfacial friction factor with 332 experimental points and obtained a correlation applicable to the ranges: gas fraction, 0.1–1; pressure, 0.09–0.6 MPa; and pipe diameter, 5–127 mm. The formula is:
in which
t is the liquid thickness,
D is the pipe diameter,
,
, the gas Froude number
, and gas superficial velocity
. According to Equation (16), the momentum exchange term can be expressed as:
The interfacial shear force is a pair of interaction forces; thus, is obtained to close Equations (6)–(14).
3.3. Validation of the Annular Flow Model
This work emphasizes accurately describing the interphase work and energy transfer in variable cross-sectional tubes; thus, the interfacial shear stress model is crucial. To verify the current model, a marine two-phase nozzle was tested, and the pressure drop and nozzle thrust were measured.
(1) Set-up
Figure 1 presents the schematic of the testing system of this nozzle. The pipeline system is divided into two parts: water supply and gas supply. The water supplied by a water tank is driven by an XWL42.15 high-pressure plunger pump produced by the AR Company of Italy with a maximum pressure of 15 MPa. After that, the water enters an NXQ hydraulic bladder accumulator, which can suppress the pressure pulsation to generate a steady water flow. The water flow rate can be adjusted in the range of 0–65 L/min by controlling the motor speed with different frequency inputs. The nitrogen supplied by a high-pressure nitrogen cylinder passes through a pressure-reducing valve to obtain a stable pressure source. Then, a sonic nozzle is used to control the gas flow rate. A check valve at the end is used to prevent backflow. Finally, the water and gas are introduced into the mixing chamber and then the test nozzle.
(2) The test nozzle
The structure of the experimental nozzle is shown in
Figure 2a. A mixing chamber, in which the water flows through a porous rectifier in the center and the radially injected gas flows through the annular slit, is fixed on the support frame. The testing nozzle is attached to the mixing chamber, and it has an inlet diameter of 11.5 mm, a nozzle length of 25 mm and an outlet diameter of 5.5 mm. The nozzle profile is a Wittosinski curve. Herein, the flow pattern in the nozzle is a unique annular flow with a central water jet surrounded by annular gas flow. In practical application of the underwater nozzle, this configuration can easily produce an annular flow without a long gas-liquid mixer.
(3) Measurement devices
The thrust, water flow rate, nitrogen gas flow rate and pressure were measured in this experiment. A turbine flow meter (LWGB-6) with a range of 0–6
and an accuracy of
F.S. was used to measure the water flow rate. The gas flow rate was indirectly calculated by the pressure before the sonic nozzle. There were six pressure measuring points distributed on the nozzle wall axially, as shown in
Figure 2a,b, which were measured by six piezoresistive pressure transmitters (Cyyz11) from Star sensor manufacturing Co., Ltd., Langfang, China. The measuring range was 0–2.5 MPa, and the accuracy class was ±
F.S. The thrust was measured by a target disk with a diameter of 72 mm, as shown in
Figure 2c. The force sensor adopted the HZC-01 sensor made by Cheng Yin company, with a range of 0–50 kg and an accuracy of ±
F.S. The pressure, flow rates and thrust signals were collected to the computer via a data acquisition card for further analysis.
(4) Testing results
As shown in
Table 1, five tests were conducted with different gas and water flow rates, and the nozzle inlet pressures were also measured. The nozzle outlet pressure in these tests was 1 atm. Based on these conditions, the nozzle inlet conditions for the two fluid model could be obtained. To predict the interfacial shear force in this unique annular flow nozzle by using the present model, a simple transformation was adopted. Its principle is shown in
Figure 3, where the two nozzles have the same diameters, and the liquid film thickness in the left nozzle (a) is identical to that of the gas in the right nozzle (b). If the flow parameters, such as velocity and density, fulfill
,
,
and
, their interfacial shear stress is considered to be identical at this cross-section. Then, the pressure variation and thrust were calculated and compared with experimental results.
Figure 4a shows the calculated and measured pressure distributions along the nozzle axis for the five tests. The results indicate that the theoretical model can predict the first five pressures well, with only the last value being larger than the experiments. The difference was mainly caused by the enhanced interphase momentum exchange near the nozzle exit, where the two phases are mixed to a larger extent due to the instability of this flow pattern. According to
Figure 4b, the predicted thrusts agreed reasonably with experiments, and the most significant deviation was about −10%. It is demonstrated that the predicting accuracy of this model was not sensitive to the boundary conditions in the current tests, and it can accurately describe the interphase work and energy transfer.
3.4. Results and Discussion for Annular Flow
The annular flow was simulated using the variable cross-section pipeline shown in
Figure 5a under the working conditions of
0.1 MPa,
0.5 kg/s and
0.042 kg/s. The inlet diameter
was 50 mm, the outlet diameter
was 20 mm, and the length diameter ratio
was 5.
In the subsonic case (shown in
Figure 5a–e), while the cross-section converged linearly as the gas and liquid flowed downstream (
Figure 5a), the relative thickness of the liquid layer increased as the gas velocity became faster than the liquid velocity (
Figure 5b), and the dimensionless interfacial stresses increased rapidly after x/L > 0.6 (
Figure 5b) due to significant interphase velocity slip. The resulting interfacial forces climbed more slowly than the stresses because of the decreasing interfacial area. The total pressure of the gas gradually decreased and that of the liquid increased at the same time (
Figure 5c), indicating that the mechanical energy of gas was transferred to liquid by the work of the interphase force.
The work factor of gas
, defined in Equation (2), was used to measure the energy conversion capability in this process. The interphase work transfer coefficient
, defined in Equation (3), represented work loss caused by interfacial slip.
Figure 5d shows that all three coefficients decreased along the coordinate direction, despite enhanced energy conversion and transfer (
Figure 5e). As shown in
Figure 5a–e, the interfacial shear force and interphase energy transfer were continuously intensified by decreasing the cross-section area of the pipeline to control expansion and accelerate the gas process, which converted more gas work into the mechanical energy of the liquid.
Figure 6a–e illustrates the influence of the inlet conditions and geometry parameters. Overall, the interphase mechanical efficiency
was the least affected of the four parameters defined in
Section 2, with an absolute variation of only 0.5%. A greater mass flow rate of gas resulted in a higher
and
but led to a fall in
, as shown in
Figure 6a, which means that more gas made it possible for the liquid to receive the work of the gas. This benefit came at the cost of reducing the capability of the gas work output. Although greater liquid mass flow can also bring a higher
,
was nearly independent of this, and
showed a slight trend to decrease (
Figure 6b).
Higher gas total pressure had a greater ability to work on the liquid, resulting in a larger
; however, the resulting greater velocity slip was not conducive to energy transfer, leading to a smaller
. Thus, the effect of the interphase work-to-energy conversion coefficient
from the total gas pressure at the inlet was relatively small, implying that a reasonably lower gas total pressure supply is sufficient in engineering design (
Figure 6c).
In the convergent nozzle, a smaller area ratio
resulted in greater convergence and higher exit velocity. According to
Figure 6d, when the area ratio
was becoming smaller,
,
,
and
were all decreasing, indicating that a high jet speed can pay the price for efficient loss of energy transfer and conversion. If the length-diameter ratio increased, the residence time of the gas became longer, resulting in a higher interphase work-to-energy conversion coefficient
. However, the interphase mechanical efficiency
decreased, as shown in
Figure 6e.
From the perspective of correlation sensitivity, the gas flow rate is the most important inlet condition that can significantly affect the work factor of the gas and the interphase work transfer coefficient . In general, the interphase mechanical energy transfer is significantly more affected by structural parameters than by inlet conditions. By changing the cross-section area of the pipeline, more gas work is transferred to the mechanical energy of liquid, resulting in greater momentum of the jet; however, the efficiency of the energy transfer is decreased. For two-phase jet propulsion, the main goal of the nozzle design is to generate a larger thrust with higher propulsion efficiency, indicating the importance of conducting fine design and optimization of two-phase flow nozzles.
4. Bubbly Flow in Variable Cross-Section Tubes
Based on the bubbly flow model derived by Wang and Chen [
25], a one-dimensional steady bubbly flow model for the variable cross-section pipeline was established, disregarding the mass and heat exchange between phases. The gas was assumed to be an adiabatic flow, and the effects of boundary layer and turbulence were ignored.
4.1. Flow Model
The respective mass equations for the gas and liquid are:
The void fraction equation reads:
For simplicity, the bubbles whose coalescence and fission are ignored are assumed to be spherical and dispersed on any cross-section of the pipe. With no mass transfer at the interface, the relationship between the bubble volume,
, and gas density,
, is obtained as
, where the subscript
indicates the pipe inlet. Thus, the mass equation for a single bubble can be expressed as a differential form with the bubble radius
R:
As for the annular flow in
Section 3, the effect of gravity on the pipeline flow was not considered, and wall friction was also ignored. The combined momentum equation for the two-phase mixture is:
For a non-deformable bubble moving in the fluid, interfacial forces can be decomposed into skin and form drag force, virtual mass, Basset, wall lubrication, lift and turbulent dispersion force [
22,
44]. The last three forces are relatively small and suitable to the multi-dimensional case; they were not considered in the current one-dimensional model. The Basset force is suitable for unsteady flow and was also excluded. The bubble momentum equation is:
where
is the substantial derivative based on the liquid velocity, and
is the drag coefficient. The first term on the left side of Equation (23) represents the virtual mass force
, the second term is the additional virtual mass force caused by radial deformation of bubbles
and the fourth term is the drag force
; the inertial force on the right side can be ignored due to the negligible mass of a bubble.
Considering the inertial effects related to the growth and collapse of the bubbles, a modified Rayleigh-Plesset equation was introduced to relate the bubble radius to the local pressure:
where
,
is the surface tension coefficient, and
is the saturation vapor pressure. (Assume that the gas is adiabatic, and its specific heat ratio
k is a constant; here,
k = 1.4.).
According to Equation (24), four interfacial stresses were induced by the expansion of the bubbles and lead to the pressure difference between gases and liquids:
Formula (25) denotes the interfacial stress due to the radial deformation of bubbles and the inertia effect of liquid, resulting in a pressure difference between liquid at the bubble surface and at infinite distance. Formula (26) is the normal stress at the bubble surface due to liquid viscosity. Formula (27) denotes interfacial stress due to surface tension. Formula (28) is the interfacial stress resulting from the average pressure of surrounding fluids caused by the motion of the bubbles.
Among these interfacial forces, the directions of virtual mass force , drag force and deformation virtual mass force are opposite to the bubble velocity; thus, they are the driving forces for liquids.
In Equation (1), the power of interphase forces on gas can be disassembled into:
Though the duration of
, the work of these forces on gas can be expressed as:
,
,
, where
is the number of bubbles over the time of
,
is the number flow rate of bubbles and
represents the work of interfacial stress (25)~(28) on gases. However, the expression of
cannot be expressed analytically, but can be obtained indirectly through Equations (1) and (29). The total pressure increments of liquid result from the work of interfacial forces, and the power can be decomposed into:
The work of each force on liquid can be expressed as: , and .
For steady flow, in Equations (22)–(24) is considered. Equations (18)–(24) Are a set of second-order, non-linear ordinary differential equations with seven unknowns (). They are closed and, given the initial conditions, can be solved by the fourth-order Runge-Kutta method with variable step size.
4.2. Validation of the Bubbly Flow Model
This bubbly flow model used in variable cross-section tubes is well-established and has been verified extensively. Wang and Chen [
25] compared this bubbly flow model with experimental data, showing that it can reasonably predict the subsonic bubbly flow and the phase relative velocity. After that, Zhang et al. [
11] utilized the same model to simulate the bubbly nozzle flow of a two-phase underwater ramjet. Compared with the experiment, the predicted thrust was consistent with the measured data, within ±20% deviations. This model, therefore, is sufficient to simulate the interphase work and energy transfer.
4.3. Results and Discussion for Bubbly Flow
The convergent pipeline shown in
Figure 5a was adopted with flow conditions of
0.1 MPa,
1.5 kg/s,
0.001 kg/s and
1 mm, and an identical inlet velocity for gas and liquid.
As illustrated in
Figure 7a, by decreasing the cross-section linearly, the bubble radius increased continuously throughout the flow, indicating that the bubble was expanding and working on the liquid. The gas and liquid continued to accelerate in the pipeline, but the bubbles moved faster than the liquid.
Figure 7b shows that the forces on a single bubble increased sharply in the second half of the pipeline, indicating that intense momentum exchange occurred. Among the interfacial forces, drag force and virtual mass force were dominant, and deformation virtual mass force was relatively small.
Figure 7c indicates that the mechanical energy of gas was transferred to liquid.
Figure 7d shows the distribution of the three coefficients along the pipeline. That the gas work factor
was always close to one shows excellent energy conversion capability for bubbly flow, and the interphase work-to-energy conversion coefficient
was affected mainly by
(i.e., the interfacial slip). As shown in
Figure 7e, almost all of the enthalpy drop of gas (black line) was converted to the work
(red line), resulting in the work factor of gas
being up to 99.75%, and only 0.25% of the gas enthalpy drop was to increase its kinetic energy (
Figure 8). About 80.32% of
was transferred to liquid (dark blue line in
Figure 7e), which mainly depended on drag force (66.25%) and virtual mass force (29.78%), and the energy transferred by the deformation virtual mass force and other forces accounted for less than 4% (
Figure 8). In this process, the main energy loss resulted from the interphase velocity slip, implying that reducing velocity slip is the key to improving work and energy transfer between phases. Noticeably, it should start by improving the drag force and virtual mass force.
For bubbly flow in the pipeline, sketched in
Figure 5a, with different inlet conditions and different geometries, the steady bubbly flow model described above was adopted to simulate the flow characteristics, as shown in
Figure 9a–d. In general, the work factor of gas
was always close to one, which means that the gas output nearly all of its available energy to liquid; thus, the variation and value of interphase work-to-energy conversion coefficient
were always similar to the interface work transfer coefficient
. Meanwhile, the interphase mechanical efficiency
also showed a slight trend of variation, with a value of about 0.99~1.
As shown in
Figure 9a, a greater mass flow rate of gas injected into the bubbly flow helped to accelerate the liquid and decrease the velocity slip, leading to a larger
, but the cost was the increment of total mechanical energy loss, which resulted in a smaller
. Though greater mass flow rate of liquid also caused a fall of
,
and
were almost unchanged, indicating that increasing the liquid flow rate apparently did not change the relative movement of the bubbles inside it, so the interphase work transfer capability remained unchanged (
Figure 9b).
For bubbly flow in a convergent nozzle, the smaller area ratio contributed to smaller
,
and
(
Figure 9c), which also showed that the greater generated thrust resulted not only from the greater interphase work loss, but also from the higher loss of total mechanical energy of the two phases. According to
Figure 9d, a larger length-diameter ratio meant a longer interacting time between the bubbles and liquid, and, thus, more work could be transferred to the liquid through interfacial forces, resulting in larger
and
, which suggests an available way to improve gas energy usage in the future.
From the perspective of correlation sensitivity, the gas flow rate was the only inlet condition that affected the interface work transfer coefficient , and, thus, the interphase work-to-energy conversion coefficient . For bubbly flow, the inlet conditions and geometry parameters mainly affected the interphase work transfer, while the energy conversion of gas was always maintained at a high level. Compared with the inlet conditions, the geometry parameters of the pipelines could impact the interphase work transfer to a greater extent, implying that the optimization design of the pipe geometry is crucial.
5. Homogeneous Flow in Mechanical Equilibrium
In the homogeneous flow model, two phases evolved in mechanical equilibrium, so the gas, liquid and interface—having the same pressure and velocity at any cross section—represented limiting momentum exchange coupled by interfacial forces. In the absence of shock waves, each phase was assumed to be an isentropic flow, which greatly simplified the model. This method of the corresponding model, derived by Kapila et al. [
45], was given by LeMartelot et al. [
27].
The flow in the pipeline was also assumed to be a one-dimensional steady flow. According to LeMartelot et al. [
27], at the left and right sides of each cross-section, the following relations read:
where the subscript
k represents the gas (
k =
g) and liquid (
k =
l), respectively. The mass fraction of gas and liquid can be expressed as
. According to SG EOS (Equations (12) and (13)) and the isentropic condition, the density of gas and liquid reads
and
, respectively. Thus, the mixture density can be expressed as a function of pressure:
.
The total enthalpy of the mixture can be given by
, and the static enthalpy of gas and liquid can be expressed by the SG EOS
and
, respectively. According to the conservation of total enthalpy, the velocity at any cross-section can be expressed as a function of pressure:
For subsonic or transonic flows, and given inlet conditions
and tube profile
A(
x), the mass flow rate can be determined from the outlet section, which depends only upon the outlet pressure:
The velocity at any cross section
of the tube can be obtained from the mass flow rate formula
, where the pressure
must be determined. According to the conservation of total enthalpy:
only the pressure
is unknown, and it can be solved using the Newton-Raphson method. Once the local pressure was obtained, other quantities were determined by using the above corresponding relations.
No pressure difference and no velocity slip mean that the momentum coupling between phases reached the maximum, and there were infinite interfacial forces to ensure an identical velocity at all times. Thus, the interface work transfer coefficient satisfied , and, theoretically, the interphase work-to-energy conversion coefficient reached its highest value in homogeneous flow.
The pipe geometry shown in
Figure 5a was utilized, and the inlet conditions were as follows: liquid flow rate, 1.43 kg/s; gas flow rate, 0.001 kg/s; and outlet pressure, 0.1 MPa. The calculated interphase work-to-energy conversion coefficient
was 0.9981, indicating that only a small part of the gas enthalpy drop was converted into its kinetic energy, and most of that was transferred to liquid. Since there was no relative slip between phases, energy loss was zero, and the interphase mechanical efficiency
was one.
In homogeneous flow, the interphase work-to-energy conversion coefficient was close to one, and the variation of flow conditions had very little effect on it; therefore, sensitivity analysis and detailed discussion are not included here.
6. Comparison among Three Flow Regimes
From annular flow and bubbly flow to homogeneous flow, the momentum coupling between phases gradually strengthened and finally achieved a maximum level. These four coefficients (defined in
Section 2 for the three flow regimes) are shown in
Table 2 and were obtained under the same pipe geometry (
Figure 5a) and flow conditions (the liquid mass flow rate was 1.3679 kg/s, the gas mass flow rate was 0.0012 kg/s, and the outlet pressure was 0.1 MPa).
Under the same conditions, the interphase work-to-energy conversion coefficient and the interphase mechanical efficiency gradually increased from annular flow to bubbly flow to homogeneous flow, indicating that the liquid in the homogeneous flow could receive more transferred energy by interphase force, paying the lowest price to the energy conversion. Thus, from the perspective of interphase mechanical energy transfer, homogeneous flow was the limit of the gas-liquid two phase flow. Additionally, although the interphase energy transfer ability of bubbly flow was obviously weaker than that of the homogeneous flow, it was much stronger than that of the annular flow, and the interphase work-to-energy conversion coefficient was nearly ten times higher than that of the annular flow.
Among the three flow patterns, the difference in the work transfer coefficient was more apparent than in the gas work factor, implying that the velocity slip is the main influencing factor on interphase work-to-energy conversion coefficient . The velocity slip between the phases in the annular flow was greatest, so its interphase work transfer coefficient was only 0.1058. A higher degree of gas dispersion in liquid contributed to a greater transfer of mechanical energy, implying that enhancing interphase mixing is the best way to reduce energy loss.
7. Conclusions
This work established a two fluid model of annular flow, a two fluid model of bubbly flow and a homogeneous flow model in which only interphase momentum exchanges were considered in all three one-dimensional models. To capture the interphase energy conversion and work transfer, an interfacial shear stress model, proposed by Aliyu et al. [
3], was adopted to close the current annular flow model. Drag force, virtual mass force and additional virtual mass force were included in the bubbly flow model; the interfacial forces due to the radial deformation of bubbles involved in the modified Rayleigh-Plesset equation were also analyzed. Based on the three flow regimes, work efficiency and energy transfer characteristics were discussed and compared systematically by introducing four new coefficients: gas work factor
, interface work transfer coefficient
, interphase work-to-energy conversion coefficient
and interphase mechanical efficiency
. The main conclusions are as follows:
(1) For both the annular and bubbly flows, the gas flow rate was the most important inlet condition to significantly promote the interphase work transfer. Geometry parameters appeared to affect interphase mechanical energy transfer similarly for these two flow patterns and—to a larger extent—in comparison with inlet conditions;
(2) By changing the cross-sectional area of the pipeline, the interphase energy transfer was intensified to generate a greater jet momentum, but and were lower compared to the constant section tube, indicating the importance of carrying out the fine design and optimization of two-phase flow nozzles. In general, a higher mass flow rate of gas, a larger length-diameter ratio or a greater area ratio can produce higher . A smaller mass flow rate of liquid, or larger area ratio, contributes to a higher ;
(3) From annular flow to bubbly flow to homogeneous flow, the interphase momentum coupling mechanisms were distinguished, the result being that the influence of the inlet conditions and geometry parameters on interphase energy conversion and transfer was not identical. The higher degree of gas dispersion in the liquid was conducive to the conversion of gas work into the mechanical energy of the liquid. Among these three flow patterns, the annular flow had the lowest interphase work-to-energy conversion coefficient and mechanical efficiency (, ); the bubbly flow was second (, ); and the homogeneous flow was the limit for interphase mechanical energy conversion because the gas-liquid momentum coupling reached the maximum, resulting in the highest conversion coefficient and mechanical efficiency (, ).
The results of this study provide insight into the application and understanding of interphase momentum transfer and mechanical energy conversion in gas-liquid variable cross-section tube flow. The work potential of gas in other flow patterns still requires further study. The influences of mass and heat transfer between phases on mechanical energy conversion are also a topic for future research.