The gas-liquid countercurrent flow pattern is complex and the bubble migration velocity is difficult to predict in the process of bullheading well killing.The experiment on bubble migration in gas-liquid countercurren...The gas-liquid countercurrent flow pattern is complex and the bubble migration velocity is difficult to predict in the process of bullheading well killing.The experiment on bubble migration in gas-liquid countercurrent flow in annulus is carried out under different working conditions to reveal how the wellbore inclination angle,liquid phase property and countercurrent liquid velocity affect the bubble deformation and bubble migration trajectory/velocity,and to establish a bubble migration velocity prediction model.The bubbles in the countercurrent flow mainly migrate in two modes:free rising of isolated bubbles,and interactive rising of multiple bubbles.The bubbles migrate by an S-shaped trajectory in the countercurrent flow.With the increase of countercurrent liquid velocity,the lateral oscillation of bubbles is intensified.The increases of wellbore inclination angle,liquid density and liquid viscosity make the bubble migration trajectory gradually to be linear.The bubble is generally ellipsoidal during its rising.The wellbore inclination angle has little effect on the degree of bubble deformation.The bubbles are ellipsoidal during rising,with little influence of wellbore inclination angle on bubble deformation.With the increase of liquid viscosity and density,the aspect ratio of the bubble decreases.As the wellbore inclination angle increases,the bubble migration velocity gradually decreases.As the liquid viscosity increases,the bubble migration velocity decreases.As the liquid density increases,the bubble migration velocity increases slightly.The established bubble migration velocity prediction model yields errors within±15%,and demonstrates broad applicability across a wide range of operating conditions.展开更多
A mathematic model of two-phase flow and a physical model of two-dimensional (2D) vertical section for the plate-type structured packing Mellapak 250.Y were set up and verified. The models were used to study the influ...A mathematic model of two-phase flow and a physical model of two-dimensional (2D) vertical section for the plate-type structured packing Mellapak 250.Y were set up and verified. The models were used to study the influence of packing’s surface microstructure on the continuity of liquid film and the amount of liquid holdup. Simulation results show that the round corner shape and micro wavy structure are favorable in remaining the continuity of liquid film and increasing the amount of liquid holdup. The appropriate liquid flow rate was determined by investigating different liquid loadings to obtain an unbroken liquid film on the packing surface. The pressure difference between inlet and outlet for gas phase allowed gas and liquid to flow countercurrently in a 2D computational domain. The direction change of gas flow occurred near the phase interface area.展开更多
In order to evaluate CCFL （countercurrent flow limitation） characteristics in a PWR （pressurized water reactor） hot leg under reflux condensation, numerical simulations have been conducted using a 2F （two-fluid）...In order to evaluate CCFL （countercurrent flow limitation） characteristics in a PWR （pressurized water reactor） hot leg under reflux condensation, numerical simulations have been conducted using a 2F （two-fluid） model and a VOF （volume of fluid） method implemented in the CFD （computational fluid dynamics） software, FLUENT6.3.26. The 2F model gave good agreement with CCFL data in low pressure conditions but did not give good results for high pressure steam-water conditions. In the previous study, the computational grid and schemes were improved in the VOF method to improve calculations in circular tubes, and the calculated CCFL characteristics agreed well with the UPTF （Upper Plenum Test Facility） data at 1.5 MPa. In this study, therefore, using the 2F model and the computational grid previously improved for the VOF calculations, numerical simulations were conducted for steam-water flows at 1.5 MPa under PWR full-scale conditions. In the range of medium gas volumetric fluxes, the calculated CCFL characteristics agreed well with the values calculated by the VOF method and the UPTF data at 1.5 MPa. This indicated that the reference set of the interfacial drag correlations employed in this study could be applied not only to low pressures but also to high pressures.展开更多
To clarify the countercurrent flow in a PWR hot leg under reflux condensation, numerical simulations of countercurrent air-water flow for a 1/15th scale model of the PWR hot leg were conducted using the two-fluid mode...To clarify the countercurrent flow in a PWR hot leg under reflux condensation, numerical simulations of countercurrent air-water flow for a 1/15th scale model of the PWR hot leg were conducted using the two-fluid model implemented in CFD software. In this paper, the effect of expansion of the inclined pipe, which is the actual plant geometry, was evaluated. When increasing the air velocity, CCFL characteristics and the mechanism of flow pattern transition had significant differences between the case with and without expansion of the inclined pipe. CCFL characteristics were mitigated in the case with expansion. The effect of computational grid size was also discussed. When the supplied water velocity was small, the predicted flow pattern transition point agreed well with the measured data by increasing the number of cells. On the other hand, when the air velocity was decreasing, there were no significant differences in each case.展开更多
The coherent structures and the chaotic phenomena in the transition of the axisymmetric countercurrent mixing shear flow were investigated experimentally. Two kinds of self-excited oscillation modes could exist in the...The coherent structures and the chaotic phenomena in the transition of the axisymmetric countercurrent mixing shear flow were investigated experimentally. Two kinds of self-excited oscillation modes could exist in the axisymmetric countercurrent mixing shear flow. One is the shear layer self-excited oscillation mode corresponding to the high Reynolds number regime and the other is the jet column self-excited oscillation mode corresponding to the low Reynolds number regime in the case of the velocity ratio ranging from I to 1.5. Analyzing the auto-power spectrum, self-correlation-function and three dimensional reconstructed phase trajectory, the route to chaos through three Hopf bifurcations intercepted by an intermittence of the dynamical system corresponding to the axisymmetric countercurrent mixing shear flow was discovered when the velocity ratio is equal to 1.32.展开更多
The complex liquid film behaviors at flooding in an inclined pipe were investigated with computational fluid dynamic(CFD) approaches. The liquid film behaviors included the dynamic wave characteristics before flooding...The complex liquid film behaviors at flooding in an inclined pipe were investigated with computational fluid dynamic(CFD) approaches. The liquid film behaviors included the dynamic wave characteristics before flooding and the transition of flow pattern when flooding happened. The influences of the surface tension and liquid viscosity were specially analyzed. Comparisons of the calculated velocity at the onset of flooding with the available experimental results showed a good agreement. The calculations verify that the fluctuation frequency and the liquid film thickness are almost unaffected by the superficial gas velocity until the flooding is triggered due to the Kelvin–Helmholtz instability. When flooding triggered at the superficial liquid velocity larger than0.15 m·s-1, the interfacial wave developed to slug flow, while it developed to entrainment flow when it was smaller than 0.08 m·s-1. The interfacial waves were more easily torn into tiny droplets with smaller surface tension, eventually evolving into the mist flow. When the liquid viscosity increases, the liquid film has a thicker holdup with more intensive fluctuations, and more likely developed to the slug flow.展开更多
The development of a simplified 2-D numerical model was described forwind-driven circulation in reservoir u-sing standard k-ε turbulence model to specify eddy viscositydistribution. The governing equations are transf...The development of a simplified 2-D numerical model was described forwind-driven circulation in reservoir u-sing standard k-ε turbulence model to specify eddy viscositydistribution. The governing equations are transformed and solved on variable vertical grids, whichallows refinement at the surface and bottom boundaries. The results of the model simulation for floware compared with analytical solutions for laminar and turbulent flows, experimental data in awind-flume and wind wave tank. The sensitivity analysis results show that use of large number ofdepth layers increases the accuracy for the bottom counter-current flow. Prediction of surface driftwas not very sensitive to surface grid refinement. The model was also applied to Baisha reservoirfor an assumed wind condition and showed to be able to simulate the general features of surfacedrift and return flow under variable flow depth. The model can serve as alternative means ofstudying wind-driven flow beside experiments. It also reduced the problem complexity associated with3-D circulation models while faithfully reproducing the drift and near bottom return currents.展开更多
基金Supported by the National Natural Science Foundation of China(U21B2069,52274020,52288101,52274022)National Key Research and Development Program of China(2022YFC2806504)。
文摘The gas-liquid countercurrent flow pattern is complex and the bubble migration velocity is difficult to predict in the process of bullheading well killing.The experiment on bubble migration in gas-liquid countercurrent flow in annulus is carried out under different working conditions to reveal how the wellbore inclination angle,liquid phase property and countercurrent liquid velocity affect the bubble deformation and bubble migration trajectory/velocity,and to establish a bubble migration velocity prediction model.The bubbles in the countercurrent flow mainly migrate in two modes:free rising of isolated bubbles,and interactive rising of multiple bubbles.The bubbles migrate by an S-shaped trajectory in the countercurrent flow.With the increase of countercurrent liquid velocity,the lateral oscillation of bubbles is intensified.The increases of wellbore inclination angle,liquid density and liquid viscosity make the bubble migration trajectory gradually to be linear.The bubble is generally ellipsoidal during its rising.The wellbore inclination angle has little effect on the degree of bubble deformation.The bubbles are ellipsoidal during rising,with little influence of wellbore inclination angle on bubble deformation.With the increase of liquid viscosity and density,the aspect ratio of the bubble decreases.As the wellbore inclination angle increases,the bubble migration velocity gradually decreases.As the liquid viscosity increases,the bubble migration velocity decreases.As the liquid density increases,the bubble migration velocity increases slightly.The established bubble migration velocity prediction model yields errors within±15%,and demonstrates broad applicability across a wide range of operating conditions.
基金Supported by the National Key Basic R&D Program ("973" Program, No. 2009CB219905 and 2009CB219907)the Program for Changjiang Scholars and Innovative Research Teams in Universities (No. IRT0936)
文摘A mathematic model of two-phase flow and a physical model of two-dimensional (2D) vertical section for the plate-type structured packing Mellapak 250.Y were set up and verified. The models were used to study the influence of packing’s surface microstructure on the continuity of liquid film and the amount of liquid holdup. Simulation results show that the round corner shape and micro wavy structure are favorable in remaining the continuity of liquid film and increasing the amount of liquid holdup. The appropriate liquid flow rate was determined by investigating different liquid loadings to obtain an unbroken liquid film on the packing surface. The pressure difference between inlet and outlet for gas phase allowed gas and liquid to flow countercurrently in a 2D computational domain. The direction change of gas flow occurred near the phase interface area.
文摘In order to evaluate CCFL （countercurrent flow limitation） characteristics in a PWR （pressurized water reactor） hot leg under reflux condensation, numerical simulations have been conducted using a 2F （two-fluid） model and a VOF （volume of fluid） method implemented in the CFD （computational fluid dynamics） software, FLUENT6.3.26. The 2F model gave good agreement with CCFL data in low pressure conditions but did not give good results for high pressure steam-water conditions. In the previous study, the computational grid and schemes were improved in the VOF method to improve calculations in circular tubes, and the calculated CCFL characteristics agreed well with the UPTF （Upper Plenum Test Facility） data at 1.5 MPa. In this study, therefore, using the 2F model and the computational grid previously improved for the VOF calculations, numerical simulations were conducted for steam-water flows at 1.5 MPa under PWR full-scale conditions. In the range of medium gas volumetric fluxes, the calculated CCFL characteristics agreed well with the values calculated by the VOF method and the UPTF data at 1.5 MPa. This indicated that the reference set of the interfacial drag correlations employed in this study could be applied not only to low pressures but also to high pressures.
文摘To clarify the countercurrent flow in a PWR hot leg under reflux condensation, numerical simulations of countercurrent air-water flow for a 1/15th scale model of the PWR hot leg were conducted using the two-fluid model implemented in CFD software. In this paper, the effect of expansion of the inclined pipe, which is the actual plant geometry, was evaluated. When increasing the air velocity, CCFL characteristics and the mechanism of flow pattern transition had significant differences between the case with and without expansion of the inclined pipe. CCFL characteristics were mitigated in the case with expansion. The effect of computational grid size was also discussed. When the supplied water velocity was small, the predicted flow pattern transition point agreed well with the measured data by increasing the number of cells. On the other hand, when the air velocity was decreasing, there were no significant differences in each case.
文摘The coherent structures and the chaotic phenomena in the transition of the axisymmetric countercurrent mixing shear flow were investigated experimentally. Two kinds of self-excited oscillation modes could exist in the axisymmetric countercurrent mixing shear flow. One is the shear layer self-excited oscillation mode corresponding to the high Reynolds number regime and the other is the jet column self-excited oscillation mode corresponding to the low Reynolds number regime in the case of the velocity ratio ranging from I to 1.5. Analyzing the auto-power spectrum, self-correlation-function and three dimensional reconstructed phase trajectory, the route to chaos through three Hopf bifurcations intercepted by an intermittence of the dynamical system corresponding to the axisymmetric countercurrent mixing shear flow was discovered when the velocity ratio is equal to 1.32.
基金Supported by the Major State Basic Research Development Program of China(2011CB706501)the National Natural Science Foundation of China(51276157)
文摘The complex liquid film behaviors at flooding in an inclined pipe were investigated with computational fluid dynamic(CFD) approaches. The liquid film behaviors included the dynamic wave characteristics before flooding and the transition of flow pattern when flooding happened. The influences of the surface tension and liquid viscosity were specially analyzed. Comparisons of the calculated velocity at the onset of flooding with the available experimental results showed a good agreement. The calculations verify that the fluctuation frequency and the liquid film thickness are almost unaffected by the superficial gas velocity until the flooding is triggered due to the Kelvin–Helmholtz instability. When flooding triggered at the superficial liquid velocity larger than0.15 m·s-1, the interfacial wave developed to slug flow, while it developed to entrainment flow when it was smaller than 0.08 m·s-1. The interfacial waves were more easily torn into tiny droplets with smaller surface tension, eventually evolving into the mist flow. When the liquid viscosity increases, the liquid film has a thicker holdup with more intensive fluctuations, and more likely developed to the slug flow.
文摘The development of a simplified 2-D numerical model was described forwind-driven circulation in reservoir u-sing standard k-ε turbulence model to specify eddy viscositydistribution. The governing equations are transformed and solved on variable vertical grids, whichallows refinement at the surface and bottom boundaries. The results of the model simulation for floware compared with analytical solutions for laminar and turbulent flows, experimental data in awind-flume and wind wave tank. The sensitivity analysis results show that use of large number ofdepth layers increases the accuracy for the bottom counter-current flow. Prediction of surface driftwas not very sensitive to surface grid refinement. The model was also applied to Baisha reservoirfor an assumed wind condition and showed to be able to simulate the general features of surfacedrift and return flow under variable flow depth. The model can serve as alternative means ofstudying wind-driven flow beside experiments. It also reduced the problem complexity associated with3-D circulation models while faithfully reproducing the drift and near bottom return currents.
