Gas-liquid two-phase flow in fractal porous media is pivotal for engineering applications,yet it remains challenging to be accurately characterized due to complex microstructure-flow interactions.This study establishe...Gas-liquid two-phase flow in fractal porous media is pivotal for engineering applications,yet it remains challenging to be accurately characterized due to complex microstructure-flow interactions.This study establishes a pore-scale numerical framework integratingMonte Carlo-generated fractal porousmedia with Volume of Fluid(VOF)simulations to unravel the coupling among pore distribution characterized by fractal dimension(Df),flow dynamics,and displacement efficiency.A pore-scale model based on the computed tomography(CT)microstructure of Berea sandstone is established,and the simulation results are compared with experimental data.Good agreement is found in phase distribution,breakthrough behavior,and flow path morphology,confirming the reliability of the numerical simulation method.Ten fractal porous media models with Df ranging from 1.25~1.7 were constructed using a Monte-Carlo approach.The gas-liquid two-phase flow dynamics was characterized using the VOF solver across gas injection rates of 0.05-5m/s,inwhich the time-resolved two-phase distribution patternswere systematically recorded.The results reveal that smaller fractal dimensions(Df=1.25~1.45)accelerate fingering breakthrough(peak velocity is 1.73 m/s at Df=1.45)due to a bimodal pore size distribution dominated by narrow channels.Increasing Df amplifies vorticity generation by about 3 times(eddy viscosity is 0.033 Pa⋅s at Df=1.7)through reduced interfacial curvature,while tortuosity-driven pressure differentials transition from sharp increases(0.4~6.3 Pa at Df=1.25~1.3)to inertial plateaus(4.8 Pa at Df=1.7).A nonlinear increase in equilibrium gas volume fraction(fav=0.692 at Df=1.7)emerges from residual gas saturation and turbulence-enhanced dispersion.This behavior is further modulated by flow velocity,with fav peaking at 0.72 under capillary-dominated conditions(0.05 m/s),but decreasing to 0.65 in the inertial regime(0.5 m/s).The work quantitatively links fractal topology to multiphase flow regimes,demonstrating the critical role of Df in governing preferential pathways,energy dissipation,and phase distribution.展开更多
In multiphase pumps transporting gas-liquid two-phase flows,the high-speed rotation of the impeller induces complex deformations in bubble shapes within the flow domain,making the prediction of gasliquid two-phase dra...In multiphase pumps transporting gas-liquid two-phase flows,the high-speed rotation of the impeller induces complex deformations in bubble shapes within the flow domain,making the prediction of gasliquid two-phase drag forces highly challenging in numerical simulations.To achieve precise prediction of the drag forces on irregular bubbles within multiphase pumps,this study modifies the existing bubble drag force model and applies the revised model to the prediction of gas-liquid two-phase flow within multiphase pumps.The research findings indicate that the modified drag force model significantly enhances the accuracy of predicting flow characteristics within the pump,particularly under high gas volume fraction conditions.The simulation results for gas phase distribution and vorticity exhibit strong agreement with experimental data.The modified drag model better captures the accumulation of the gas phase at the suction side of the impeller outlet.It also accurately predicts the vortex characteristics induced by bubble backflow from the trailing edges of the diffuser.Additionally,the adjustment of the drag coefficient enhances the model’s ability to represent local flow field characteristics,thereby optimizing the performance simulation methods of multiphase pumps.Compared to traditional drag force models,the modified model reduces prediction errors in head and efficiency by 36.4%and 27.5%,respectively.These results provide important theoretical foundations and model support for improving the accuracy of gas-liquid two-phase flow simulations and optimizing the design of multiphase pumps under high gas volume fraction conditions.展开更多
By combining with an improved model on engraving process,a two-phase flow interior ballistic model has been proposed to accurately predict the flow and energy conversion behaviors of pyrotechnic actuators.Using comput...By combining with an improved model on engraving process,a two-phase flow interior ballistic model has been proposed to accurately predict the flow and energy conversion behaviors of pyrotechnic actuators.Using computational fluid dynamics(CFD),the two-phase flow and piston engraving characteristics of a pyrotechnic actuator are investigated.Initially,the current model was utilized to examine the intricate,multi-dimensional flow,and energy conversion characteristics of the propellant grains and combustion gas within the pyrotechnic actuator chamber.It was discovered that the combustion gas on the wall's constant transition from potential to kinetic energy,along with the combined effect of the propellant motion,are what create the pressure oscillation within the chamber.Additionally,a numerical analysis was conducted to determine the impact of various parameters on the pressure oscillation and piston motion,including pyrotechnic charge,pyrotechnic particle size,and chamber structural dimension.The findings show that decreasing the pyrotechnic charge will lower the terminal velocity,while increasing and decreasing the pyrotechnic particle size will reduce the pressure oscillation in the chamber.The pyrotechnic particle size has minimal bearing on the terminal velocity.The results of this investigation offer a trustworthy forecasting instrument for comprehending and creating pyrotechnic actuator designs.展开更多
The gas-droplet two-phase reacting flow in a model combustor with the V-gutter flame holder is studied by an Eulerian-Lagrangian large-eddy simulation (LES) approach. The k-equation subgrid-scale model is used to simu...The gas-droplet two-phase reacting flow in a model combustor with the V-gutter flame holder is studied by an Eulerian-Lagrangian large-eddy simulation (LES) approach. The k-equation subgrid-scale model is used to simulate the subgrid eddy viscosity, and the eddy-break-up (EBU) combustion subgrid-scale model is used to determine the chemical reaction rate. A two-step turbulent combustion subgrid-scale model is employed for calculating carbon monoxide CO concentration, and the NO subgrid-scale pollutant formation model for the evaluation of the rate of NO formation. The heat flux model is applied to the prediction of radiant heat transfer. The gas phase is solved with the SIMPLE algorithm and a hybrid scheme in the staggered grid system. The liquid phase equations are solved in a Lagrangian frame in reference of the particle-source-in-cell (PSIC) algorithm. From simulation results, the exchange of mass, moment and energy between gas and particle fields for the reacting flow in the afterburner with a V-gutter flame holder can be obtained. By the comparison of experimental and simulation results, profile temperature and pollutant of the outlet are quite in agreement with experimental data. Results show that the LES approach for predicting the two-phase instantaneous reacting flow and pollutant emissions in the afterburner is feasible.展开更多
The influence of the squeeze film between the tube and the support structure on flow-induced vibrations is a critical factor in tube bundles subjected to two-phase cross-flow.This aspect can significantly alter the th...The influence of the squeeze film between the tube and the support structure on flow-induced vibrations is a critical factor in tube bundles subjected to two-phase cross-flow.This aspect can significantly alter the threshold for fluidelastic instability and affect heat transfer efficiency.This paper presents a mathematical model incorporating the squeeze film force between the tube and the support structure.We aim to clarify the mechanisms underlying fluidelastic instability in tube bundle systems exposed to two-phase flow.Using a self-developed computer program,we performed numerical calculations to examine the influence of the squeeze film on the threshold of fluidelastic instability in the tube bundle system.Furthermore,we analyzed how the thickness and length of the squeeze film affect both the underlying mechanisms and the critical velocity of fluidelastic instability.展开更多
This work investigated the dynamic behavior of vertical pipes conveying gas-liquid two-phase flow when subjected to external excitations at both ends.Even with minimal excitation amplitude,resonance can occur when the...This work investigated the dynamic behavior of vertical pipes conveying gas-liquid two-phase flow when subjected to external excitations at both ends.Even with minimal excitation amplitude,resonance can occur when the excitation frequency aligns with the natural frequency of the pipe,significantly increasing the degree of operational risk.The governing equation of motion based on the Euler-Bernoulli beam is derived for the relative deflection with stationary simply supported ends,with the effects of the external excitations represented by source terms distributed along the pipe length.The fourth-order partial differential equation is solved via the generalized integral transform technique(GITT),with the solution successfully verified via comparison with results in the literature.A comprehensive analysis of the vibration phenomena and changes in the motion state of the pipe is conducted for three classes of external excitation conditions:same frequency and amplitude(SFSA),same frequency but different amplitudes(SFDA),and different frequencies and amplitudes(DFDA).The numerical results show that with increasing gas volume fraction,the position corresponding to the maximum vibration displacement shifts upward.Compared with conditions without external excitation,the vibration displacement of the pipe conveying two-phase flow under external excitation increases significantly.The frequency of external excitation has a significant effect on the dynamic behavior of a pipe conveying two-phase flow.展开更多
With the increasing miniaturization of systems and surging demand for power density,accurate prediction and control of two-phase flow pressure drop have become a core challenge restricting the performance of microchan...With the increasing miniaturization of systems and surging demand for power density,accurate prediction and control of two-phase flow pressure drop have become a core challenge restricting the performance of microchannel heat exchangers.Pressure drop,a critical hydraulic characteristic,serves as both a natural constraint for cooling systems and determines the power required to pump the working fluid through microchannels.This paper reviews the characteristics,prediction models,and optimization measures of two-phase flow pressure drop for low-boiling-point working fluids in microchannels.It systematically analyzes key influencing factors such as fluid physical properties,operating conditions,channel geometry,and flow patterns,and discusses the complex mechanisms of pressure drop under the coupling effect of multi-physical fields.Mainstream prediction models are reviewed:the homogeneous flow model simplifies calculations but shows large deviations at low quality;the separated flow model considers interphase interactions and can be applied to micro-scales after modification;the flow-pattern-based model performs zoned modeling but relies on subjective classification;machine learning improves prediction accuracy but faces the“black-box”problem.In terms of optimization,channel designs are improved through porous structures and micro-rib arrays,and flow rate distribution is optimized using splitters to balance pressure drop and heat transfer performance.This study provides theoretical support for microchannel thermal management in high-power-density devices.展开更多
Polymer flooding is an important means of improving oil recovery and is widely used in Daqing,Xinjiang,and Shengli oilfields,China.Different from conventional injection media such as water and gas,viscoelastic polymer...Polymer flooding is an important means of improving oil recovery and is widely used in Daqing,Xinjiang,and Shengli oilfields,China.Different from conventional injection media such as water and gas,viscoelastic polymer solutions exhibit non-Newtonian and nonlinear flow behavior including shear thinning and shear thickening,polymer convection,diffusion,adsorption,retention,inaccessible pore volume,and reduced effective permeability.However,available well test model of polymer flooding wells generally simplifies these characteristics on pressure transient response,which may lead to inaccurate results.This work proposes a novel two-phase numerical well test model to better describe the polymer viscoelasticity and nonlinear flow behavior.Different influence factors that related to near-well blockage during polymer flooding process,including the degree of blockage(inner zone permeability),the extent of blockage(composite radius),and polymer flooding front radius are explored to investigate these impacts on bottom hole pressure responses.Results show that polymer viscoelasticity has a significant impact on the transitional flow segment of type curves,and the effects of near-well formation blockage and polymer concentration distribution on well test curves are very similar.Thus,to accurately interpret the degree of near-well blockage in injection wells,it is essential to first eliminate the influence of polymer viscoelasticity.Finally,a field case is comprehensively analyzed and discussed to illustrate the applicability of the proposed model.展开更多
Clayey-silt natural gas hydrate reservoirs in the South China Sea exhibit loose and unconsolidated structures, heterogeneous pore structures, high clay mineral contents, and strong hydrophilicity. These characteristic...Clayey-silt natural gas hydrate reservoirs in the South China Sea exhibit loose and unconsolidated structures, heterogeneous pore structures, high clay mineral contents, and strong hydrophilicity. These characteristics complicate the gas-water two-phase flow process in porous media following hydrate decomposition, posing challenges for efficient development. This study examines the transport response of clayey-silt reservoir samples from the Shenhu area using gas-water two-phase flow experiments and CT scanning to explore changes in pore structure, gas-water distribution, and relative permeability under varying flow conditions. The results indicate that pore heterogeneity significantly influences flow characteristics. Gas preferentially displaces water in larger pores, forming fracture-like pores, which serve as preferential flow channels for gas migration. The preferential flow channels enhance gas-phase permeability up to 19 times that of the water phase when fluid pressures exceed total stresses. However,small pores retain liquid, leading to a high residual water saturation of 0.561. CT imaging reveals that these hydro-fractures improve gas permeability but also confine gas flow to specific channels. Pore network analysis shows that gas injection expands the pore-throat network, enhancing connectivity and forming fracture-like pores. Residual water remains trapped in smaller pores and throats, while structural changes, including new fractures, improve gas flow pathways and overall connectivity. Relative permeability curves demonstrate a narrow gas-water cocurrent-flow zone, a right-shifted iso-permeability point and high reservoir capillary pressure, indicating a strong "water-blocking" effect. The findings suggest that optimizing reservoir stimulation techniques to enhance fracture formation, reduce residual water saturation, and improve gas flow capacity is critical for efficient hydrate reservoir development.展开更多
A second-order moment two-phase turbulence model for simulating dense gas-particle flows （USM-Θ model）, combining the unified second-order moment twophase turbulence model for dilute gas-particle flows with the kin...A second-order moment two-phase turbulence model for simulating dense gas-particle flows （USM-Θ model）, combining the unified second-order moment twophase turbulence model for dilute gas-particle flows with the kinetic theory of particle collision, is proposed. The interaction between gas and particle turbulence is simulated using the transport equation of two-phase velocity correlation with a two-time-scale dissipation closure. The proposed model is applied to simulate dense gas-particle flows in a horizontal channel and a downer. Simulation results and their comparison with experimental results show that the model accounting for both anisotropic particle turbulence and particle-particle collision is obviously better than models accounting for only particle turbulence or only particle-particle collision. The USM-Θ model is also better than the k-ε-kp-Θ model and the k-ε-kp-εp-Θ model in that the first model can simulate the redistribution of anisotropic particle Reynolds stress components due to inter-particle collision, whereas the second and third models cannot.展开更多
Shale gas production involves complex gas-water two-phase flow,with flow patterns in proppant-filled fractures playing a critical role in determining production efficiency.In this study,3D geometric models of 40/70 me...Shale gas production involves complex gas-water two-phase flow,with flow patterns in proppant-filled fractures playing a critical role in determining production efficiency.In this study,3D geometric models of 40/70 mesh ceramic particles and quartz sand proppant clusters were elaborated using computed tomography(CT)scanning.These models were used to develop a numerical simulation framework based on the lattice Boltzmann method(LBM),enabling the investigation of gas-water flow behavior within proppant-filled fractures under varying driving forces and surface tensions.Simulation results at a closure pressure of 15 MPa have revealed that ceramic particles exhibit a simpler and more porous internal structure than quartz sand of the same size.Under identical flow conditions,ceramic proppants demonstrate higher fluid replacement efficiency.Replacement efficiency increases with higher porosity,greater driving force,and lower surface tension.Furthermore,fluid displacement is strongly influenced by pore geometry:flow is faster in straighter and wider channels,with preferential movement through larger pores forming dominant flow paths.The replacement velocity exhibits a characteristic time evolution,initially rapid,then gradually decreasing,correlating positively with the development of these dominant channels.展开更多
Cleat serves as the primary flow pathway for coalbed methane(CBM)and water.However,few studies consider the impact of local contact on two-phase flow within cleats.A visual generalized model of endogenous cleats was c...Cleat serves as the primary flow pathway for coalbed methane(CBM)and water.However,few studies consider the impact of local contact on two-phase flow within cleats.A visual generalized model of endogenous cleats was constructed based on microfluidics.A microscopic and mesoscopic observation technique was proposed to simultaneously capture gas-liquid interface morphology of pores and throat and the two-phase flow characteristics in entire cleat system.The local contact characteristics of cleats reduced absolute permeability,which resulted in a sharp increase in the starting pressure.The reduced gas flow capacity narrowed the co-infiltration area and decreased water saturation at the isotonic point in a hydrophilic environment.The increased local contact area of cleats weakened gas phase flow capacity and narrowed the co-infiltration area.Jumping events occurred in methane-water flow due to altered porosity caused by local contact in cleats.The distribution of residual phases changed the jumping direction on the micro-scale as well as the dominant channel on the mesoscale.Besides,jumping events caused additional energy dissipation,which was ignored in traditional two-phase flow models.This might contribute to the overestimation of relative permeability.The work provides new methods and insights for investigating unsaturated flow in complex porous media.展开更多
Vibration of flexible pipelines in the marine environment affects the flow characteristics of the transported materials inside the pipelines,which is related to transportation efficiency and energy consumption,thereby...Vibration of flexible pipelines in the marine environment affects the flow characteristics of the transported materials inside the pipelines,which is related to transportation efficiency and energy consumption,thereby necessitating further investigation.In this study,the flow characteristics of particle-liquid two-phase flow transported upward in flexible pipelines are investigated based on the computational fluid dynamics-discrete element method(CFD-DEM).Typical forms of vibration including standing wave vibration and traveling wave vibration are employed and compared with a stationary pipeline.Results reveal that particles in the upward-traveling-wave vibrating pipeline still mainly distribute in the middle of the pipeline,while particles in the standing-wave vibrating pipeline exhibit periodic transverse aggregation near the pipe wall,and the fluctuations of particle concentration and particle z-direction velocity over time in each cross section of the pipeline are more obviously suppressed.When the propagation direction of the vibration wave changes from the same direction as the particle transport to static and then to the opposite direction,its hindering and regulating effect on the particles gradually increases,and the pipeline pressure drop gradually decreases.展开更多
Deep-sea mineral resource transportation predominantly utilizes hydraulic pipeline methodology.Environmental factors induce vibrations in flexible pipelines,thereby affecting the internal flow characteristics.Therefor...Deep-sea mineral resource transportation predominantly utilizes hydraulic pipeline methodology.Environmental factors induce vibrations in flexible pipelines,thereby affecting the internal flow characteristics.Therefore,real-time monitoring of solid–liquid two-phase flow in pipelines is crucial for system maintenance.This study develops an autoencoder-based deep learning framework to reconstruct three-dimensional solid–liquid two-phase flow within flexible vibrating pipelines utilizing sparse wall information from sensors.Within this framework,separate X-model and F-model with distinct hidden-layer structures are established to reconstruct the coordinates and flow field information on the computational domain grid of the pipeline under traveling wave vibration.Following hyperparameter optimization,the models achieved high reconstruction accuracy,demonstrating R^(2)values of 0.990 and 0.945,respectively.The models’robustness is evaluated across three aspects:vibration parameters,physical fields,and vibration modes,demonstrating good reconstruction performance.Results concerning sensors show that 20 sensors(0.06%of total grids)achieve a balance between accuracy and cost,with superior accuracy obtained when arranged along the full length of the pipe compared to a dense arrangement at the front end.The models exhibited a signal-to-noise ratio tolerance of approximately 27 dB,with reconstruction accuracy being more affected by sensor failures at both ends of the pipeline.展开更多
The volume fraction of the solid and liquid phase of debris flows, which evolves simultaneously across terrains, largely determines the dynamic property of debris flows. The entrainment process significantly influence...The volume fraction of the solid and liquid phase of debris flows, which evolves simultaneously across terrains, largely determines the dynamic property of debris flows. The entrainment process significantly influences the amplitude of the volume fraction. In this paper, we present a depth-averaged two-phase debris-flow model describing the simultaneous evolution of the phase velocity and depth, the solid and fluid volume fractions and the bed morphological evolution. The model employs the Mohr–Coulomb plasticity for the solid stress, and the fluid stress is modeled as a Newtonian viscous stress. The interfacial momentum transfer includes viscous drag and buoyancy. A new extended entrainment rate formula that satisfies the boundary momentum jump condition （Iverson and Ouyang, 2015） is presented. In this formula, the basal traction stress is a function of the solid volume fraction and can take advantage of both the Coulomb and velocity-dependent friction models. A finite volume method using Roe’s Riemann approximation is suggested to solve the equations. Three computational cases are conducted and compared with experiments or previous results. The results show that the current computational model and framework are robust and suitable for capturing the characteristics of debris flows.展开更多
In photothermal power(solar energy)generation systems,purging residual molten salt from pipelines using highpressure gas poses a significant challenge,particularly in clearing the bottom of regulating valves.Ineffecti...In photothermal power(solar energy)generation systems,purging residual molten salt from pipelines using highpressure gas poses a significant challenge,particularly in clearing the bottom of regulating valves.Ineffective purging can lead to crystallization of the molten salt,resulting in blockages.To address this issue,understanding the gas-liquid two-phase flow dynamics during high-pressure gas purging is crucial.This study utilizes the Volume of Fluid(VOF)model and adaptive dynamic grids to simulate the gas-liquid two-phase flow during the purging process in a DN50 PN50 conventional molten salt regulating valve.Initially,the reliability of the CFD simulations is validated through comparisons with experimental data and findings from the literature.Subsequently,simulation experiments are conducted to analyze the effects of various factors,including purge flow rates,initial liquid accumulation masses,purge durations,and the profiles of the valve bottom flow channels.The results indicate that the purging process comprises four distinct stages:Initial violent surge stage,liquid discharge stage,liquid partial fallback stage,liquid dissipation stage.For an initial liquid height of 17 mm at the bottom of the valve,the critical purge flow rate lies between 3 and 5 m/s.Notably,the critical purge flow rate is independent of the initial liquid accumulation mass.As the purge gas flow rate increases,the volume of liquid discharged also increases.Beyond the critical purge flow rate,higher purge gas velocities lead to shorter purge durations.Interestingly,the residual liquid mass after purging remains unaffected by the initial liquid accumulation.Additionally,the flow channel profile at the bottom of the valve significantly influences both the critical purge speed and the efficiency of the purging process.展开更多
Gas-solid two-phase turbulent flows,mass transfer,heat transfer and catalytic cracking reactions areknown to exert interrelated influences in commercial fluid catalytic cracking(FCC)riser reactors.In the presentpaper,...Gas-solid two-phase turbulent flows,mass transfer,heat transfer and catalytic cracking reactions areknown to exert interrelated influences in commercial fluid catalytic cracking(FCC)riser reactors.In the presentpaper,a three-dimensional turbulent gas-solid two-phase flow-reaction model for FCC riser reactors was devel-oped.The model took into account the gas-solid two-phase turbulent flows,inter-phase heat transfer,masstransfer,catalytic cracking reactions and their interrelated influence.The k-V-k_P two-phase turbulence modelwas employed and modified for the two-phase turbulent flow patterns with relatively high particle concentration.Boundary conditions for the flow-reaction model were given.Related numerical algorithm was formed and a nu-merical code was drawn up.Numerical modeling for commercial FCC riser reactors could be carried out with thepresented model.展开更多
Merits of the Fernandes model(Fernandes et al.1983)for two-phase slug flow in verticaltubes are reviewed in this paper.While predicting many macroscopic parameters of slug flow in verti-cal tubes,it fails to present c...Merits of the Fernandes model(Fernandes et al.1983)for two-phase slug flow in verticaltubes are reviewed in this paper.While predicting many macroscopic parameters of slug flow in verti-cal tubes,it fails to present correctly the trend that the average voidage in liquid slugs increases asthe rising velocity of Taylor bubbles is increased.It is also desirable to extend its application toelectrolyte systems, and to churn flow conditions.Based on the diagnostic analysis,the model equa-tion for gas entrainment by falling liquid film is reformulated and the influence of surface tension isalso accounted for.Development of the falling liquid film is recognized in the revised model in or-der to suit the case of short Taylor bubbles as well.The modified model predicts the variation of av-erage voidage in liquid slugs in good agreement with available experimental data.展开更多
The two-phase flow models are commonly used in industrial applications, such as nuclear, power, chemical-process, oil-and-gas, cryogenics, bio-medical, micro-technology and so on. This is a survey paper on the study o...The two-phase flow models are commonly used in industrial applications, such as nuclear, power, chemical-process, oil-and-gas, cryogenics, bio-medical, micro-technology and so on. This is a survey paper on the study of compressible nonconservative two-fluid model, drift-flux model and viscous liquid-gas two-phase flow model. We give the research developments of these three two-phase flow models, respectively. In the last part, we give some open problems about the above models.展开更多
A full second-order moment (FSM) model and an algebraic stress (ASM) two-phase turbulence modelare proposed and applied to predict turbulent bubble-liquid flows in a 2D rectangular bubble column. Predictiongives the b...A full second-order moment (FSM) model and an algebraic stress (ASM) two-phase turbulence modelare proposed and applied to predict turbulent bubble-liquid flows in a 2D rectangular bubble column. Predictiongives the bubble and liquid velocities, bubble volume fraction, bubble and liquid Reynolds stresses and bubble-liquidvelocity correlation. For predicted two-phase velocities and bubble volume fraction there is only slight differencebetween these two models, and the simulation results using both two models are in good agreement with the particleimage velocimetry (PIV) measurements. Although the predicted two-phase Reynolds stresses using the FSM are insomewhat better agreement with the PIV measurements than those predicted using the ASM, the Reynolds stressespredicted using both two models are in general agreement with the experiments. Therefore, it is suggested to usethe ASM two-phase turbulence model in engineering application for saving the computation time.展开更多
基金funded by the National Key R&D Program of China,China(Grant No.2023YFB4005500)National Natural Science Foundation of China,China(Grant Nos.52379113 and 52379114).
文摘Gas-liquid two-phase flow in fractal porous media is pivotal for engineering applications,yet it remains challenging to be accurately characterized due to complex microstructure-flow interactions.This study establishes a pore-scale numerical framework integratingMonte Carlo-generated fractal porousmedia with Volume of Fluid(VOF)simulations to unravel the coupling among pore distribution characterized by fractal dimension(Df),flow dynamics,and displacement efficiency.A pore-scale model based on the computed tomography(CT)microstructure of Berea sandstone is established,and the simulation results are compared with experimental data.Good agreement is found in phase distribution,breakthrough behavior,and flow path morphology,confirming the reliability of the numerical simulation method.Ten fractal porous media models with Df ranging from 1.25~1.7 were constructed using a Monte-Carlo approach.The gas-liquid two-phase flow dynamics was characterized using the VOF solver across gas injection rates of 0.05-5m/s,inwhich the time-resolved two-phase distribution patternswere systematically recorded.The results reveal that smaller fractal dimensions(Df=1.25~1.45)accelerate fingering breakthrough(peak velocity is 1.73 m/s at Df=1.45)due to a bimodal pore size distribution dominated by narrow channels.Increasing Df amplifies vorticity generation by about 3 times(eddy viscosity is 0.033 Pa⋅s at Df=1.7)through reduced interfacial curvature,while tortuosity-driven pressure differentials transition from sharp increases(0.4~6.3 Pa at Df=1.25~1.3)to inertial plateaus(4.8 Pa at Df=1.7).A nonlinear increase in equilibrium gas volume fraction(fav=0.692 at Df=1.7)emerges from residual gas saturation and turbulence-enhanced dispersion.This behavior is further modulated by flow velocity,with fav peaking at 0.72 under capillary-dominated conditions(0.05 m/s),but decreasing to 0.65 in the inertial regime(0.5 m/s).The work quantitatively links fractal topology to multiphase flow regimes,demonstrating the critical role of Df in governing preferential pathways,energy dissipation,and phase distribution.
基金funded by Sichuan Natural Science Foundation Outstanding Youth Science Foundation(No.2024NSFJQ0012)Key project of Regional Innovation and Development Joint Fund of National Natural Science Foundation(No.U23A20669)Sichuan Science and Technology Program(2022ZDZX0041).
文摘In multiphase pumps transporting gas-liquid two-phase flows,the high-speed rotation of the impeller induces complex deformations in bubble shapes within the flow domain,making the prediction of gasliquid two-phase drag forces highly challenging in numerical simulations.To achieve precise prediction of the drag forces on irregular bubbles within multiphase pumps,this study modifies the existing bubble drag force model and applies the revised model to the prediction of gas-liquid two-phase flow within multiphase pumps.The research findings indicate that the modified drag force model significantly enhances the accuracy of predicting flow characteristics within the pump,particularly under high gas volume fraction conditions.The simulation results for gas phase distribution and vorticity exhibit strong agreement with experimental data.The modified drag model better captures the accumulation of the gas phase at the suction side of the impeller outlet.It also accurately predicts the vortex characteristics induced by bubble backflow from the trailing edges of the diffuser.Additionally,the adjustment of the drag coefficient enhances the model’s ability to represent local flow field characteristics,thereby optimizing the performance simulation methods of multiphase pumps.Compared to traditional drag force models,the modified model reduces prediction errors in head and efficiency by 36.4%and 27.5%,respectively.These results provide important theoretical foundations and model support for improving the accuracy of gas-liquid two-phase flow simulations and optimizing the design of multiphase pumps under high gas volume fraction conditions.
基金supported by the National Natural Science Foundation of China(Grant No.11972194).
文摘By combining with an improved model on engraving process,a two-phase flow interior ballistic model has been proposed to accurately predict the flow and energy conversion behaviors of pyrotechnic actuators.Using computational fluid dynamics(CFD),the two-phase flow and piston engraving characteristics of a pyrotechnic actuator are investigated.Initially,the current model was utilized to examine the intricate,multi-dimensional flow,and energy conversion characteristics of the propellant grains and combustion gas within the pyrotechnic actuator chamber.It was discovered that the combustion gas on the wall's constant transition from potential to kinetic energy,along with the combined effect of the propellant motion,are what create the pressure oscillation within the chamber.Additionally,a numerical analysis was conducted to determine the impact of various parameters on the pressure oscillation and piston motion,including pyrotechnic charge,pyrotechnic particle size,and chamber structural dimension.The findings show that decreasing the pyrotechnic charge will lower the terminal velocity,while increasing and decreasing the pyrotechnic particle size will reduce the pressure oscillation in the chamber.The pyrotechnic particle size has minimal bearing on the terminal velocity.The results of this investigation offer a trustworthy forecasting instrument for comprehending and creating pyrotechnic actuator designs.
文摘The gas-droplet two-phase reacting flow in a model combustor with the V-gutter flame holder is studied by an Eulerian-Lagrangian large-eddy simulation (LES) approach. The k-equation subgrid-scale model is used to simulate the subgrid eddy viscosity, and the eddy-break-up (EBU) combustion subgrid-scale model is used to determine the chemical reaction rate. A two-step turbulent combustion subgrid-scale model is employed for calculating carbon monoxide CO concentration, and the NO subgrid-scale pollutant formation model for the evaluation of the rate of NO formation. The heat flux model is applied to the prediction of radiant heat transfer. The gas phase is solved with the SIMPLE algorithm and a hybrid scheme in the staggered grid system. The liquid phase equations are solved in a Lagrangian frame in reference of the particle-source-in-cell (PSIC) algorithm. From simulation results, the exchange of mass, moment and energy between gas and particle fields for the reacting flow in the afterburner with a V-gutter flame holder can be obtained. By the comparison of experimental and simulation results, profile temperature and pollutant of the outlet are quite in agreement with experimental data. Results show that the LES approach for predicting the two-phase instantaneous reacting flow and pollutant emissions in the afterburner is feasible.
基金financially supported by the National Natural Science Foundation of China(Grant No.12072336).
文摘The influence of the squeeze film between the tube and the support structure on flow-induced vibrations is a critical factor in tube bundles subjected to two-phase cross-flow.This aspect can significantly alter the threshold for fluidelastic instability and affect heat transfer efficiency.This paper presents a mathematical model incorporating the squeeze film force between the tube and the support structure.We aim to clarify the mechanisms underlying fluidelastic instability in tube bundle systems exposed to two-phase flow.Using a self-developed computer program,we performed numerical calculations to examine the influence of the squeeze film on the threshold of fluidelastic instability in the tube bundle system.Furthermore,we analyzed how the thickness and length of the squeeze film affect both the underlying mechanisms and the critical velocity of fluidelastic instability.
基金financially supported by the Key Research and Development Program of Shandong Province(Grant Nos.2022CXGC020405,2023CXGC010415 and 2025TSGCCZZB0238)the National Natural Science Foundation of China(Grant No.52171288)the financial support from CNPq,FAPERJ,ANP,Embrapii,and China National Petroleum Corporation(CNPC).
文摘This work investigated the dynamic behavior of vertical pipes conveying gas-liquid two-phase flow when subjected to external excitations at both ends.Even with minimal excitation amplitude,resonance can occur when the excitation frequency aligns with the natural frequency of the pipe,significantly increasing the degree of operational risk.The governing equation of motion based on the Euler-Bernoulli beam is derived for the relative deflection with stationary simply supported ends,with the effects of the external excitations represented by source terms distributed along the pipe length.The fourth-order partial differential equation is solved via the generalized integral transform technique(GITT),with the solution successfully verified via comparison with results in the literature.A comprehensive analysis of the vibration phenomena and changes in the motion state of the pipe is conducted for three classes of external excitation conditions:same frequency and amplitude(SFSA),same frequency but different amplitudes(SFDA),and different frequencies and amplitudes(DFDA).The numerical results show that with increasing gas volume fraction,the position corresponding to the maximum vibration displacement shifts upward.Compared with conditions without external excitation,the vibration displacement of the pipe conveying two-phase flow under external excitation increases significantly.The frequency of external excitation has a significant effect on the dynamic behavior of a pipe conveying two-phase flow.
基金supported by the Beijing Municipal Science&Technology Commission(Z231100006123010).
文摘With the increasing miniaturization of systems and surging demand for power density,accurate prediction and control of two-phase flow pressure drop have become a core challenge restricting the performance of microchannel heat exchangers.Pressure drop,a critical hydraulic characteristic,serves as both a natural constraint for cooling systems and determines the power required to pump the working fluid through microchannels.This paper reviews the characteristics,prediction models,and optimization measures of two-phase flow pressure drop for low-boiling-point working fluids in microchannels.It systematically analyzes key influencing factors such as fluid physical properties,operating conditions,channel geometry,and flow patterns,and discusses the complex mechanisms of pressure drop under the coupling effect of multi-physical fields.Mainstream prediction models are reviewed:the homogeneous flow model simplifies calculations but shows large deviations at low quality;the separated flow model considers interphase interactions and can be applied to micro-scales after modification;the flow-pattern-based model performs zoned modeling but relies on subjective classification;machine learning improves prediction accuracy but faces the“black-box”problem.In terms of optimization,channel designs are improved through porous structures and micro-rib arrays,and flow rate distribution is optimized using splitters to balance pressure drop and heat transfer performance.This study provides theoretical support for microchannel thermal management in high-power-density devices.
基金supported by the National Natural Science Foundation of China(52104049)the Young Elite Scientist Sponsorship Program by Beijing Association for Science and Technology(BYESS2023262)。
文摘Polymer flooding is an important means of improving oil recovery and is widely used in Daqing,Xinjiang,and Shengli oilfields,China.Different from conventional injection media such as water and gas,viscoelastic polymer solutions exhibit non-Newtonian and nonlinear flow behavior including shear thinning and shear thickening,polymer convection,diffusion,adsorption,retention,inaccessible pore volume,and reduced effective permeability.However,available well test model of polymer flooding wells generally simplifies these characteristics on pressure transient response,which may lead to inaccurate results.This work proposes a novel two-phase numerical well test model to better describe the polymer viscoelasticity and nonlinear flow behavior.Different influence factors that related to near-well blockage during polymer flooding process,including the degree of blockage(inner zone permeability),the extent of blockage(composite radius),and polymer flooding front radius are explored to investigate these impacts on bottom hole pressure responses.Results show that polymer viscoelasticity has a significant impact on the transitional flow segment of type curves,and the effects of near-well formation blockage and polymer concentration distribution on well test curves are very similar.Thus,to accurately interpret the degree of near-well blockage in injection wells,it is essential to first eliminate the influence of polymer viscoelasticity.Finally,a field case is comprehensively analyzed and discussed to illustrate the applicability of the proposed model.
基金the National Natural Science Foundation of China (Nos. 42302143, 42172159)China Geological Survey Project (No. DD20211350)support from the G. Albert Shoemaker endowment
文摘Clayey-silt natural gas hydrate reservoirs in the South China Sea exhibit loose and unconsolidated structures, heterogeneous pore structures, high clay mineral contents, and strong hydrophilicity. These characteristics complicate the gas-water two-phase flow process in porous media following hydrate decomposition, posing challenges for efficient development. This study examines the transport response of clayey-silt reservoir samples from the Shenhu area using gas-water two-phase flow experiments and CT scanning to explore changes in pore structure, gas-water distribution, and relative permeability under varying flow conditions. The results indicate that pore heterogeneity significantly influences flow characteristics. Gas preferentially displaces water in larger pores, forming fracture-like pores, which serve as preferential flow channels for gas migration. The preferential flow channels enhance gas-phase permeability up to 19 times that of the water phase when fluid pressures exceed total stresses. However,small pores retain liquid, leading to a high residual water saturation of 0.561. CT imaging reveals that these hydro-fractures improve gas permeability but also confine gas flow to specific channels. Pore network analysis shows that gas injection expands the pore-throat network, enhancing connectivity and forming fracture-like pores. Residual water remains trapped in smaller pores and throats, while structural changes, including new fractures, improve gas flow pathways and overall connectivity. Relative permeability curves demonstrate a narrow gas-water cocurrent-flow zone, a right-shifted iso-permeability point and high reservoir capillary pressure, indicating a strong "water-blocking" effect. The findings suggest that optimizing reservoir stimulation techniques to enhance fracture formation, reduce residual water saturation, and improve gas flow capacity is critical for efficient hydrate reservoir development.
基金the Special Funds for Major State Basic Research of China(G-1999-0222-08)the National Natural Science Foundation of China(50376004)Ph.D.Program Foundation,Ministry of Education of China(20030007028)
文摘A second-order moment two-phase turbulence model for simulating dense gas-particle flows （USM-Θ model）, combining the unified second-order moment twophase turbulence model for dilute gas-particle flows with the kinetic theory of particle collision, is proposed. The interaction between gas and particle turbulence is simulated using the transport equation of two-phase velocity correlation with a two-time-scale dissipation closure. The proposed model is applied to simulate dense gas-particle flows in a horizontal channel and a downer. Simulation results and their comparison with experimental results show that the model accounting for both anisotropic particle turbulence and particle-particle collision is obviously better than models accounting for only particle turbulence or only particle-particle collision. The USM-Θ model is also better than the k-ε-kp-Θ model and the k-ε-kp-εp-Θ model in that the first model can simulate the redistribution of anisotropic particle Reynolds stress components due to inter-particle collision, whereas the second and third models cannot.
文摘Shale gas production involves complex gas-water two-phase flow,with flow patterns in proppant-filled fractures playing a critical role in determining production efficiency.In this study,3D geometric models of 40/70 mesh ceramic particles and quartz sand proppant clusters were elaborated using computed tomography(CT)scanning.These models were used to develop a numerical simulation framework based on the lattice Boltzmann method(LBM),enabling the investigation of gas-water flow behavior within proppant-filled fractures under varying driving forces and surface tensions.Simulation results at a closure pressure of 15 MPa have revealed that ceramic particles exhibit a simpler and more porous internal structure than quartz sand of the same size.Under identical flow conditions,ceramic proppants demonstrate higher fluid replacement efficiency.Replacement efficiency increases with higher porosity,greater driving force,and lower surface tension.Furthermore,fluid displacement is strongly influenced by pore geometry:flow is faster in straighter and wider channels,with preferential movement through larger pores forming dominant flow paths.The replacement velocity exhibits a characteristic time evolution,initially rapid,then gradually decreasing,correlating positively with the development of these dominant channels.
基金the financial support from the National Natural Science Foundation of China (No.42102127)the Postdoctoral Research Foundation of China (No.2024 M751860)。
文摘Cleat serves as the primary flow pathway for coalbed methane(CBM)and water.However,few studies consider the impact of local contact on two-phase flow within cleats.A visual generalized model of endogenous cleats was constructed based on microfluidics.A microscopic and mesoscopic observation technique was proposed to simultaneously capture gas-liquid interface morphology of pores and throat and the two-phase flow characteristics in entire cleat system.The local contact characteristics of cleats reduced absolute permeability,which resulted in a sharp increase in the starting pressure.The reduced gas flow capacity narrowed the co-infiltration area and decreased water saturation at the isotonic point in a hydrophilic environment.The increased local contact area of cleats weakened gas phase flow capacity and narrowed the co-infiltration area.Jumping events occurred in methane-water flow due to altered porosity caused by local contact in cleats.The distribution of residual phases changed the jumping direction on the micro-scale as well as the dominant channel on the mesoscale.Besides,jumping events caused additional energy dissipation,which was ignored in traditional two-phase flow models.This might contribute to the overestimation of relative permeability.The work provides new methods and insights for investigating unsaturated flow in complex porous media.
基金supported by the National Natural Science Foundation of China(Grant Nos.52471293,12372270)the National Youth Science Foundation of China(Grant Nos.52101322,52108375)+3 种基金the National Key Research and Development Program of China(Grant No.2023YFC2811600)supported by the Program for Intergovernmental International S&T Cooperation Projects of Shanghai Municipality(Grant No.22160710200)the Oceanic Interdisciplinary Program of Shanghai Jiao Tong University(Grant No.SL2022PT101)the Open Fund of the State Key Laboratory of Coastal and Offshore Engineering,Dalian University of Technology(Grant No.LP2415).
文摘Vibration of flexible pipelines in the marine environment affects the flow characteristics of the transported materials inside the pipelines,which is related to transportation efficiency and energy consumption,thereby necessitating further investigation.In this study,the flow characteristics of particle-liquid two-phase flow transported upward in flexible pipelines are investigated based on the computational fluid dynamics-discrete element method(CFD-DEM).Typical forms of vibration including standing wave vibration and traveling wave vibration are employed and compared with a stationary pipeline.Results reveal that particles in the upward-traveling-wave vibrating pipeline still mainly distribute in the middle of the pipeline,while particles in the standing-wave vibrating pipeline exhibit periodic transverse aggregation near the pipe wall,and the fluctuations of particle concentration and particle z-direction velocity over time in each cross section of the pipeline are more obviously suppressed.When the propagation direction of the vibration wave changes from the same direction as the particle transport to static and then to the opposite direction,its hindering and regulating effect on the particles gradually increases,and the pipeline pressure drop gradually decreases.
基金financial support by the National Natural Science Foundation of China(Nos.52471293 and 12372270)the National Youth Science Foundation of China(Nos.52101322 and 52108375)+3 种基金the Program for Intergovernmental International S&T Cooperation Projects of Shanghai Municipality,China(Nos.24510711100 and 22160710200)The Oceanic Interdisciplinary Program of Shanghai Jiao Tong University(No.SL2022PT101)is also gratefully acknowledgedfunded by the Open Fund of the State Key Laboratory of Coastal and Offshore Engineering of Dalian University of Technology(No.LP2415)the National Key R&D Program of China(No.2023YFC2811600)is gratefully acknowledged.
文摘Deep-sea mineral resource transportation predominantly utilizes hydraulic pipeline methodology.Environmental factors induce vibrations in flexible pipelines,thereby affecting the internal flow characteristics.Therefore,real-time monitoring of solid–liquid two-phase flow in pipelines is crucial for system maintenance.This study develops an autoencoder-based deep learning framework to reconstruct three-dimensional solid–liquid two-phase flow within flexible vibrating pipelines utilizing sparse wall information from sensors.Within this framework,separate X-model and F-model with distinct hidden-layer structures are established to reconstruct the coordinates and flow field information on the computational domain grid of the pipeline under traveling wave vibration.Following hyperparameter optimization,the models achieved high reconstruction accuracy,demonstrating R^(2)values of 0.990 and 0.945,respectively.The models’robustness is evaluated across three aspects:vibration parameters,physical fields,and vibration modes,demonstrating good reconstruction performance.Results concerning sensors show that 20 sensors(0.06%of total grids)achieve a balance between accuracy and cost,with superior accuracy obtained when arranged along the full length of the pipe compared to a dense arrangement at the front end.The models exhibited a signal-to-noise ratio tolerance of approximately 27 dB,with reconstruction accuracy being more affected by sensor failures at both ends of the pipeline.
基金Financial support from NSFC(Grant No.41572303,4151001059,41101008)Key Projects in the National Science & Technology Pillar Program(2014BAL05B01)CAS "Light of West China" Program
文摘The volume fraction of the solid and liquid phase of debris flows, which evolves simultaneously across terrains, largely determines the dynamic property of debris flows. The entrainment process significantly influences the amplitude of the volume fraction. In this paper, we present a depth-averaged two-phase debris-flow model describing the simultaneous evolution of the phase velocity and depth, the solid and fluid volume fractions and the bed morphological evolution. The model employs the Mohr–Coulomb plasticity for the solid stress, and the fluid stress is modeled as a Newtonian viscous stress. The interfacial momentum transfer includes viscous drag and buoyancy. A new extended entrainment rate formula that satisfies the boundary momentum jump condition （Iverson and Ouyang, 2015） is presented. In this formula, the basal traction stress is a function of the solid volume fraction and can take advantage of both the Coulomb and velocity-dependent friction models. A finite volume method using Roe’s Riemann approximation is suggested to solve the equations. Three computational cases are conducted and compared with experiments or previous results. The results show that the current computational model and framework are robust and suitable for capturing the characteristics of debris flows.
文摘In photothermal power(solar energy)generation systems,purging residual molten salt from pipelines using highpressure gas poses a significant challenge,particularly in clearing the bottom of regulating valves.Ineffective purging can lead to crystallization of the molten salt,resulting in blockages.To address this issue,understanding the gas-liquid two-phase flow dynamics during high-pressure gas purging is crucial.This study utilizes the Volume of Fluid(VOF)model and adaptive dynamic grids to simulate the gas-liquid two-phase flow during the purging process in a DN50 PN50 conventional molten salt regulating valve.Initially,the reliability of the CFD simulations is validated through comparisons with experimental data and findings from the literature.Subsequently,simulation experiments are conducted to analyze the effects of various factors,including purge flow rates,initial liquid accumulation masses,purge durations,and the profiles of the valve bottom flow channels.The results indicate that the purging process comprises four distinct stages:Initial violent surge stage,liquid discharge stage,liquid partial fallback stage,liquid dissipation stage.For an initial liquid height of 17 mm at the bottom of the valve,the critical purge flow rate lies between 3 and 5 m/s.Notably,the critical purge flow rate is independent of the initial liquid accumulation mass.As the purge gas flow rate increases,the volume of liquid discharged also increases.Beyond the critical purge flow rate,higher purge gas velocities lead to shorter purge durations.Interestingly,the residual liquid mass after purging remains unaffected by the initial liquid accumulation.Additionally,the flow channel profile at the bottom of the valve significantly influences both the critical purge speed and the efficiency of the purging process.
文摘Gas-solid two-phase turbulent flows,mass transfer,heat transfer and catalytic cracking reactions areknown to exert interrelated influences in commercial fluid catalytic cracking(FCC)riser reactors.In the presentpaper,a three-dimensional turbulent gas-solid two-phase flow-reaction model for FCC riser reactors was devel-oped.The model took into account the gas-solid two-phase turbulent flows,inter-phase heat transfer,masstransfer,catalytic cracking reactions and their interrelated influence.The k-V-k_P two-phase turbulence modelwas employed and modified for the two-phase turbulent flow patterns with relatively high particle concentration.Boundary conditions for the flow-reaction model were given.Related numerical algorithm was formed and a nu-merical code was drawn up.Numerical modeling for commercial FCC riser reactors could be carried out with thepresented model.
文摘Merits of the Fernandes model(Fernandes et al.1983)for two-phase slug flow in verticaltubes are reviewed in this paper.While predicting many macroscopic parameters of slug flow in verti-cal tubes,it fails to present correctly the trend that the average voidage in liquid slugs increases asthe rising velocity of Taylor bubbles is increased.It is also desirable to extend its application toelectrolyte systems, and to churn flow conditions.Based on the diagnostic analysis,the model equa-tion for gas entrainment by falling liquid film is reformulated and the influence of surface tension isalso accounted for.Development of the falling liquid film is recognized in the revised model in or-der to suit the case of short Taylor bubbles as well.The modified model predicts the variation of av-erage voidage in liquid slugs in good agreement with available experimental data.
基金supported by the National Natural Science Foundation of China(11722104,11671150)supported by the National Natural Science Foundation of China(11571280,11331005)+3 种基金supported by the National Natural Science Foundation of China(11331005,11771150)by GDUPS(2016)the Fundamental Research Funds for the Central Universities of China(D2172260)FANEDD No.201315
文摘The two-phase flow models are commonly used in industrial applications, such as nuclear, power, chemical-process, oil-and-gas, cryogenics, bio-medical, micro-technology and so on. This is a survey paper on the study of compressible nonconservative two-fluid model, drift-flux model and viscous liquid-gas two-phase flow model. We give the research developments of these three two-phase flow models, respectively. In the last part, we give some open problems about the above models.
基金Supported by the Special Funds for Major State Basic Research Projects, PRC(G1999-0222-08) and the National Natural Science Foundation of China(No. 19872039).
文摘A full second-order moment (FSM) model and an algebraic stress (ASM) two-phase turbulence modelare proposed and applied to predict turbulent bubble-liquid flows in a 2D rectangular bubble column. Predictiongives the bubble and liquid velocities, bubble volume fraction, bubble and liquid Reynolds stresses and bubble-liquidvelocity correlation. For predicted two-phase velocities and bubble volume fraction there is only slight differencebetween these two models, and the simulation results using both two models are in good agreement with the particleimage velocimetry (PIV) measurements. Although the predicted two-phase Reynolds stresses using the FSM are insomewhat better agreement with the PIV measurements than those predicted using the ASM, the Reynolds stressespredicted using both two models are in general agreement with the experiments. Therefore, it is suggested to usethe ASM two-phase turbulence model in engineering application for saving the computation time.
