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.展开更多
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.展开更多
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.展开更多
As space technology advances,thermal control systems must effectively collect and dissipate heat from distributed,multi-source environments.Loop heat pipe is a highly reliable two-phase heat transfer component,but it ...As space technology advances,thermal control systems must effectively collect and dissipate heat from distributed,multi-source environments.Loop heat pipe is a highly reliable two-phase heat transfer component,but it has several limitations when addressing multi-source heat dissipation.Inspired by the transport and heat dissipation system of plants,large trees achieve stable and efficient liquid supply under the influence of two driving forces:capillary force during transpiration in the leaves(pull)and root pressure generated by osmotic pressure in the roots(push).The root pressure provides an effective liquid supply with a driving force exceeding 2 MPa,far greater than the driving force in conventional capillary-pumped two-phase loops.Research has shown that osmotic heat pipes offer a powerful driving force,and combining osmotic pressure with capillary force has significant advantages.Therefore,this paper designs a multi-evaporator,dual-drive two-phase loop,using both osmotic pressure and capillary force to solve the multi-source heat dissipation challenge.First,a transmembrane water flux model for the osmotic pressure-driven device was established to predict the maximum heat transfer capacity of the dual-drive two-phase loop.Then,an experimental setup for a multi-evaporator“osmotic pressure+capillary force”dual-drive two-phase loop was constructed,capable of transferring at least 235 W of power under a reverse gravity condition of 20 m.The study also analyzed the effects of reverse gravity height,heat load distribution among the three evaporators,startup sequence,and varying branch resistances on the performance of the dual-drive two-phase loop.展开更多
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.展开更多
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.展开更多
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.展开更多
In this study, the three-dimensional non-premixed two-phase kerosene/air rotating detonation engines with different isolator configurations and throat area ratios are simulated by the Eulerian-Lagrangian method. The e...In this study, the three-dimensional non-premixed two-phase kerosene/air rotating detonation engines with different isolator configurations and throat area ratios are simulated by the Eulerian-Lagrangian method. The effects of the divergence, straight, and convergence isolators on the rotating detonation wave dynamics and the upstream oblique shock wave propagation mechanism are analyzed. The differences in the rotating detonation wave behaviors between ground and flight operations are clarified.The results indicate that the propagation regimes of the upstream oblique shock wave depend on the isolator configurations and operation conditions. With a divergence isolator, the airflow is accelerated throughout the isolator and divergence section, leading to a maximum Mach number(~1.8) before the normal shock. The total pressure loss reaches the largest, and the detonation pressure drops. The upstream oblique shock wave can be suppressed within the divergence section with the divergence isolator.However, for the straight and convergence isolators, the airflow in the isolator with a larger ψ_(1)(0.3 and0.4) can suffer from the disturbance of the upstream oblique shock wave. The critical incident angle is around 39° at ground operation conditions. The upstream oblique shock wave tends to be suppressed when the engine operates under flight operation conditions. The critical pressure ratio β_(cr0) is found to be able to help in distinguishing the propagation regimes of the upstream oblique shock wave. Slightly below or above the β_(cr0) can obtain different marginal propagation results. The high-speed airflow in the divergence section affects the fuel droplet penetration distance, which deteriorates the reactant mixing and the detonation area. Significant detonation velocity deficits are observed and the maximum velocity deficit reaches 26%. The results indicate the engine channel design should adopt different isolator configurations based on the purpose of total pressure loss or disturbance suppression. This study can provide useful guidance for the channel design of a more complete two-phase rotating detonation engine.展开更多
Two-phase partitioning bioreactors(TPPBs)have been widely used because they overcome the mass-transfer limitation of hydrophobic volatile organic compounds(VOCs)in waste gas biological treatments.Understanding the mec...Two-phase partitioning bioreactors(TPPBs)have been widely used because they overcome the mass-transfer limitation of hydrophobic volatile organic compounds(VOCs)in waste gas biological treatments.Understanding the mechanisms of mass-transfer enhancement in TPPBs would enable efficient predictions for further industrial applications.In this study,influences of gradually increasing silicone oil ratio on the TPPB was explored,and a 94.35%reduction of the n-hexane partition coefficient was observed with 0.1 vol.%silicone,which increased to 80.7%along with a 40-fold removal efficiency enhancement in the stabilised removal period.The elimination capacity increased from 1.47 to 148.35 g/(m^(3)·h),i.e.a 101-fold increase compared with that of the single-phase reactors,when 10 vol.%(3 Critical Micelle Concentration)silicone oil was added.The significantly promoted partition coefficient was the main reason for the mass transfer enhancement,which covered the negative influences of the decreased total mass-transfer coefficient with increasing silicone oil volume ratio.The gradually rising stirring rate was benefit to the n-hexane removal,which became negative when the dominant resistance shifted from mass transfer to biodegradation.Moreover,a mass-transfer-reaction kinetic model of the TPPB was constructed based on the balance of n-hexane concentration,dissolved oxygen and biomass.Similar to the mechanism,the partition factor was predicted sensitive to the removal performance,and another five sensitive parameters were found simultaneously.This forecasting method enables the optimisation of TPPB performance and provides theoretical support for hydrophobic VOCs degradation.展开更多
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.展开更多
Stratified flow is a common phenomenon in horizontal tubes of two-phase flow systems. However, the existing methods for calculating the wetted angle of the flat interface model and the central angle of the two-circle ...Stratified flow is a common phenomenon in horizontal tubes of two-phase flow systems. However, the existing methods for calculating the wetted angle of the flat interface model and the central angle of the two-circle model rely on solving implicit transcendental equations, which require iterative numerical root-finding methods,thereby introducing computational complexity and inefficiency. This paper proposes the high-precision explicit approximate solutions for the two models, directly correlating the geometric parameters with the flow parameters, thus significantly enhancing the efficiency and accuracy of two-phase flow analysis.展开更多
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 exosomes hold significant potential in disease diagnosis and therapeutic interventions.The objective of this study was to investigate the potential of aqueous two-phase systems(ATPSs)for the separation of bovine m...The exosomes hold significant potential in disease diagnosis and therapeutic interventions.The objective of this study was to investigate the potential of aqueous two-phase systems(ATPSs)for the separation of bovine milk exosomes.The milk exosome partition behaviors and bovine milk separation were investigated,and the ATPSs and bovine milk whey addition was optimized.The optimal separation conditions were identified as 16%(mass)polyethylene glycol 4000,10%(mass)dipotassium phosphate,and 1%(mass)enzymatic hydrolysis bovine milk whey.During the separation process,bovine milk exosomes were predominantly enriched in the interphase,while protein impurities were primarily found in the bottom phase.The process yielded bovine milk exosomes of 2.0×10^(11)particles per ml whey with high purity(staining rate>90%,7.01×10^(10)particles per mg protein)and high uniformity(polydispersity index<0.03).The isolated exosomes were characterized and identified by transmission electron microscopy,zeta potential and size distribution.The results demonstrated aqueous two-phase extraction possesses a robust capability for the enrichment and separation of exosomes directly from bovine milk whey,presenting a novel approach for the large-scale isolation of exosomes.展开更多
Liquid film cooling as an advanced cooling technology is widely used in space vehicles.Stable operation of liquid film along the rocket combustion inner wall is crucial for thermal protection of rocket engines.The sta...Liquid film cooling as an advanced cooling technology is widely used in space vehicles.Stable operation of liquid film along the rocket combustion inner wall is crucial for thermal protection of rocket engines.The stability of liquid film is mainly determined by the characteristics of interfacial wave,which is rarely investigated right now.How to improve the stability of thin film has become a hot spot.In view of this,an advanced model based on the conventional Volume of Fluid(VOF)model is adopted to investigate the characteristics of interfacial wave in gas-liquid flow by using OpenFOAM,and the mechanism of formation and development of wave is revealed intuitively through numerical study.The effects from gas velocity,surface tension and dynamic viscosity of liquid(three factors)on the wave are studied respectively.It can be found that the gas velocity is critical to the formation and development of wave,and four modes of droplets generation are illustrated in this paper.Besides,a gas vortex near the gas-liquid interface can induce formation of wave easily,so changing the gas vortex state can regulate formation and development of wave.What’s more,the change rules of three factors influencing on the interfacial wave are obtained,and the surface tension has a negative effect on the formation and development of wave only when the surface tension coefficient is above the critical value,whereas the dynamic viscosity has a positive effect in this process.Lastly,the maximum height and maximum slope angle of wave will level off as the gas velocity increases.Meanwhile,the maximum slope angle of wave is usually no more than 38°,no matter what happens to the three factors.展开更多
Fluid-structure interaction(FSI)of gas-liquid two-phase fow in the horizontal pipe is investigated numerically in the present study.The volume of fluid model and standard k-e turbulence model are integrated to simulat...Fluid-structure interaction(FSI)of gas-liquid two-phase fow in the horizontal pipe is investigated numerically in the present study.The volume of fluid model and standard k-e turbulence model are integrated to simulate the typical gas-liquid two-phase fow patterns.First,validation of the numerical model is conducted and the typical fow patterns are consistent with the Baker chart.Then,the FSI framework is established to investigate the dynamic responses of the interaction between the horizontal pipe and gas-liquid two-phase fow.The results show that the dynamic response under stratified fow condition is relatively flat and the maximum pipe deformation and equivalent stress are 1.8 mm and 7.5 MPa respectively.Meanwhile,the dynamic responses induced by slug fow,wave fow and annular fow show obvious periodic fuctuations.Furthermore,the dynamic response characteristics under slug flow condition are maximum;the maximum pipe deformation and equivalent stress can reach 4mm and 17.5 MPa,respectively.The principal direction of total deformation is different under various flow patterns.Therefore,the periodic equivalent stress will form the cyclic impact on the pipe wall and affect the fatigue life of the horizontal pipe.The present study may serve as a reference for FSI simulation under gas-liquid two-phase transport conditions.展开更多
Accurate prediction of the frictional pressure drop is important for the design and operation of subsea oil and gas transporting system considering the length of the pipeline. The applicability of the correlations to ...Accurate prediction of the frictional pressure drop is important for the design and operation of subsea oil and gas transporting system considering the length of the pipeline. The applicability of the correlations to pipeline-riser flow needs evaluation since the flow condition in pipeline-riser is quite different from the original data where they were derived from. In the present study, a comprehensive evaluation of 24prevailing correlation in predicting frictional pressure drop is carried out based on experimentally measured data of air-water and air-oil two-phase flows in pipeline-riser. Experiments are performed in a system having different configuration of pipeline-riser with the inclination of the downcomer varied from-2°to-5°to investigated the effect of the elbow on the frictional pressure drop in the riser. The inlet gas velocity ranges from 0.03 to 6.2 m/s, and liquid velocity varies from 0.02 to 1.3 m/s. A total of885 experimental data points including 782 on air-water flows and 103 on air-oil flows are obtained and used to access the prediction ability of the correlations. Comparison of the predicted results with the measured data indicate that a majority of the investigated correlations under-predict the pressure drop on severe slugging. The result of this study highlights the requirement of new method considering the effect of pipe layout on the frictional pressure drop.展开更多
Understanding fingering, as a challenge to stable displacement during the immiscible flow, has become a crucial phenomenon for geological carbon sequestration, enhanced oil recovery, and groundwater protection. Typica...Understanding fingering, as a challenge to stable displacement during the immiscible flow, has become a crucial phenomenon for geological carbon sequestration, enhanced oil recovery, and groundwater protection. Typically governed by gravity, viscous and capillary forces, these factors lead invasive fluids to occupy pore space irregularly and incompletely. Previous studies have demonstrated capillary numbers,describing the viscous and capillary forces, to quantificationally induce evolution of invasion patterns.While the evolution mechanisms of invasive patterns have not been deeply elucidated under the constant capillary number and three variable parameters including velocity, viscosity, and interfacial tension.Our research employs two horizontal visualization systems and a two-phase laminar flow simulation to investigate the tendency of invasive pattern transition by various parameters at the pore scale. We showed that increasing invasive viscosity or reducing interfacial tension in a homogeneous pore space significantly enhanced sweep efficiency, under constant capillary number. Additionally, in the fingering crossover pattern, the region near the inlet was prone to capillary fingering with multi-directional invasion, while the viscous fingering with unidirectional invasion was more susceptible occurred in the region near the outlet. Furthermore, increasing invasive viscosity or decreasing invasive velocity and interfacial tension promoted the extension of viscous fingering from the outlet to the inlet, presenting that the subsequent invasive fluid flows toward the outlet. In the case of invasive trunk along a unidirectional path, the invasive flow increased exponentially closer to the outlet, resulting in a significant decrease in the width of the invasive interface. Our work holds promising applications for optimizing invasive patterns in heterogeneous porous media.展开更多
Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinni...Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinning and shear thickening,polymer convection,diffusion,adsorption retention,inaccessible pore volume and reduced effective permeability.Meanwhile,the flux density and fracture conductivity along the hydraulic fracture are generally non-uniform due to the effects of pressure distribution,formation damage,and proppant breakage.In this paper,we present an oil-water two-phase flow model that captures these complex non-Newtonian and nonlinear behavior,and non-uniform fracture characteristics in fractured polymer flooding.The hydraulic fracture is firstly divided into two parts:high-conductivity fracture near the wellbore and low-conductivity fracture in the far-wellbore section.A hybrid grid system,including perpendicular bisection(PEBI)and Cartesian grid,is applied to discrete the partial differential flow equations,and the local grid refinement method is applied in the near-wellbore region to accurately calculate the pressure distribution and shear rate of polymer solution.The combination of polymer behavior characterizations and numerical flow simulations are applied,resulting in the calculation for the distribution of water saturation,polymer concentration and reservoir pressure.Compared with the polymer flooding well with uniform fracture conductivity,this non-uniform fracture conductivity model exhibits the larger pressure difference,and the shorter bilinear flow period due to the decrease of fracture flow ability in the far-wellbore section.The field case of the fall-off test demonstrates that the proposed method characterizes fracture characteristics more accurately,and yields fracture half-lengths that better match engineering reality,enabling a quantitative segmented characterization of the near-wellbore section with high fracture conductivity and the far-wellbore section with low fracture conductivity.The novelty of this paper is the analysis of pressure performances caused by the fracture dynamics and polymer rheology,as well as an analysis method that derives formation and fracture parameters based on the pressure and its derivative curves.展开更多
A numerical study based on a two-dimensional two-phase SPH(Smoothed Particle Hydrodynamics)model to analyze the action of water waves on open-type sea access roads is presented.The study is a continuation of the analy...A numerical study based on a two-dimensional two-phase SPH(Smoothed Particle Hydrodynamics)model to analyze the action of water waves on open-type sea access roads is presented.The study is a continuation of the analyses presented by Chen et al.(2022),in which the sea access roads are semi-immersed.In this new configuration,the sea access roads are placed above the still water level,therefore the presence of the air phase becomes a relevant issue in the determination of the wave forces acting on the structures.Indeed,the comparison of wave forces on the open-type sea access roads obtained from the single and two-phase SPH models with the experimental results shows that the latter are in much better agreement.So in the numerical simulations,a two-phaseδ-SPH model is adopted to investigate the dynamical problems.Based on the numerical results,the maximum horizontal and uplifting wave forces acting on the sea access roads are analyzed by considering different wave conditions and geometries of the structures.In particular,the presence of the girder is analyzed and the differences in the wave forces due to the air cushion effects which are created below the structure are highlighted.展开更多
Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale...Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale gas well considering the two-phase flow of gas and water.The model accounts for the influence of natural fractures and matrix properties on the fracturing process and directly applies post-fracturing formation pressure and water saturation distribution to subsequent well shut-in and production simulation,allowing for a more accurate fracturing-production integrated simulation.The results show that the reservoir physical properties have great impacts on fracture propagation,and the reasonable prediction of formation pressure and reservoir fluid distribution after the fracturing is critical to accurately predict the gas and fluid production of the shale gas wells.Compared with the conventional method,the proposed model can more accurately simulate the water and gas production by considering the impact of fracturing on both matrix pressure and water saturation.The established model is applied to the integrated fracturing-production simulation of practical horizontal shale gas wells.The simulation results are in good agreement with the practical production data,thus verifying the accuracy of the model.展开更多
基金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 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.
基金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.
基金Science Foundation for Distinguished Young Scholars 2020-JCJQ-ZQ-042 Intelligent and Bionic Spacecraft Thermal Control Technology Inspired by Tree Sap Transport Principle.
文摘As space technology advances,thermal control systems must effectively collect and dissipate heat from distributed,multi-source environments.Loop heat pipe is a highly reliable two-phase heat transfer component,but it has several limitations when addressing multi-source heat dissipation.Inspired by the transport and heat dissipation system of plants,large trees achieve stable and efficient liquid supply under the influence of two driving forces:capillary force during transpiration in the leaves(pull)and root pressure generated by osmotic pressure in the roots(push).The root pressure provides an effective liquid supply with a driving force exceeding 2 MPa,far greater than the driving force in conventional capillary-pumped two-phase loops.Research has shown that osmotic heat pipes offer a powerful driving force,and combining osmotic pressure with capillary force has significant advantages.Therefore,this paper designs a multi-evaporator,dual-drive two-phase loop,using both osmotic pressure and capillary force to solve the multi-source heat dissipation challenge.First,a transmembrane water flux model for the osmotic pressure-driven device was established to predict the maximum heat transfer capacity of the dual-drive two-phase loop.Then,an experimental setup for a multi-evaporator“osmotic pressure+capillary force”dual-drive two-phase loop was constructed,capable of transferring at least 235 W of power under a reverse gravity condition of 20 m.The study also analyzed the effects of reverse gravity height,heat load distribution among the three evaporators,startup sequence,and varying branch resistances on the performance of the dual-drive two-phase loop.
基金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.
文摘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 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.
基金supported by the National Natural Science Foundation of China (Grant No. 12202204)the Natural Science Foundation of Jiangsu Province (Grant No. BK20220953)+1 种基金the Fundamental Research Funds for the Central Universitiesthe Science and Technology Association's Young Talent Nurturing Program of Jiangsu Province (Grant No. JSTJ-2024-004)
文摘In this study, the three-dimensional non-premixed two-phase kerosene/air rotating detonation engines with different isolator configurations and throat area ratios are simulated by the Eulerian-Lagrangian method. The effects of the divergence, straight, and convergence isolators on the rotating detonation wave dynamics and the upstream oblique shock wave propagation mechanism are analyzed. The differences in the rotating detonation wave behaviors between ground and flight operations are clarified.The results indicate that the propagation regimes of the upstream oblique shock wave depend on the isolator configurations and operation conditions. With a divergence isolator, the airflow is accelerated throughout the isolator and divergence section, leading to a maximum Mach number(~1.8) before the normal shock. The total pressure loss reaches the largest, and the detonation pressure drops. The upstream oblique shock wave can be suppressed within the divergence section with the divergence isolator.However, for the straight and convergence isolators, the airflow in the isolator with a larger ψ_(1)(0.3 and0.4) can suffer from the disturbance of the upstream oblique shock wave. The critical incident angle is around 39° at ground operation conditions. The upstream oblique shock wave tends to be suppressed when the engine operates under flight operation conditions. The critical pressure ratio β_(cr0) is found to be able to help in distinguishing the propagation regimes of the upstream oblique shock wave. Slightly below or above the β_(cr0) can obtain different marginal propagation results. The high-speed airflow in the divergence section affects the fuel droplet penetration distance, which deteriorates the reactant mixing and the detonation area. Significant detonation velocity deficits are observed and the maximum velocity deficit reaches 26%. The results indicate the engine channel design should adopt different isolator configurations based on the purpose of total pressure loss or disturbance suppression. This study can provide useful guidance for the channel design of a more complete two-phase rotating detonation engine.
基金supported by the National Key Research and Development Program of China(No.2022YFC3702000)the National Natural Science Foundation of China(No.52070169)the Project of Bureau of Science and Technology of Zhoushan,China(No.2022C41013).
文摘Two-phase partitioning bioreactors(TPPBs)have been widely used because they overcome the mass-transfer limitation of hydrophobic volatile organic compounds(VOCs)in waste gas biological treatments.Understanding the mechanisms of mass-transfer enhancement in TPPBs would enable efficient predictions for further industrial applications.In this study,influences of gradually increasing silicone oil ratio on the TPPB was explored,and a 94.35%reduction of the n-hexane partition coefficient was observed with 0.1 vol.%silicone,which increased to 80.7%along with a 40-fold removal efficiency enhancement in the stabilised removal period.The elimination capacity increased from 1.47 to 148.35 g/(m^(3)·h),i.e.a 101-fold increase compared with that of the single-phase reactors,when 10 vol.%(3 Critical Micelle Concentration)silicone oil was added.The significantly promoted partition coefficient was the main reason for the mass transfer enhancement,which covered the negative influences of the decreased total mass-transfer coefficient with increasing silicone oil volume ratio.The gradually rising stirring rate was benefit to the n-hexane removal,which became negative when the dominant resistance shifted from mass transfer to biodegradation.Moreover,a mass-transfer-reaction kinetic model of the TPPB was constructed based on the balance of n-hexane concentration,dissolved oxygen and biomass.Similar to the mechanism,the partition factor was predicted sensitive to the removal performance,and another five sensitive parameters were found simultaneously.This forecasting method enables the optimisation of TPPB performance and provides theoretical support for hydrophobic VOCs degradation.
文摘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.
基金supported by the General Research Fund from the Research Grants Council of the Hong Kong Special Administrative Region of China (No. PolyU 15210624)。
文摘Stratified flow is a common phenomenon in horizontal tubes of two-phase flow systems. However, the existing methods for calculating the wetted angle of the flat interface model and the central angle of the two-circle model rely on solving implicit transcendental equations, which require iterative numerical root-finding methods,thereby introducing computational complexity and inefficiency. This paper proposes the high-precision explicit approximate solutions for the two models, directly correlating the geometric parameters with the flow parameters, thus significantly enhancing the efficiency and accuracy of two-phase flow analysis.
基金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.
基金supported by the National Natural Science Foundation of China(22378350).
文摘The exosomes hold significant potential in disease diagnosis and therapeutic interventions.The objective of this study was to investigate the potential of aqueous two-phase systems(ATPSs)for the separation of bovine milk exosomes.The milk exosome partition behaviors and bovine milk separation were investigated,and the ATPSs and bovine milk whey addition was optimized.The optimal separation conditions were identified as 16%(mass)polyethylene glycol 4000,10%(mass)dipotassium phosphate,and 1%(mass)enzymatic hydrolysis bovine milk whey.During the separation process,bovine milk exosomes were predominantly enriched in the interphase,while protein impurities were primarily found in the bottom phase.The process yielded bovine milk exosomes of 2.0×10^(11)particles per ml whey with high purity(staining rate>90%,7.01×10^(10)particles per mg protein)and high uniformity(polydispersity index<0.03).The isolated exosomes were characterized and identified by transmission electron microscopy,zeta potential and size distribution.The results demonstrated aqueous two-phase extraction possesses a robust capability for the enrichment and separation of exosomes directly from bovine milk whey,presenting a novel approach for the large-scale isolation of exosomes.
文摘Liquid film cooling as an advanced cooling technology is widely used in space vehicles.Stable operation of liquid film along the rocket combustion inner wall is crucial for thermal protection of rocket engines.The stability of liquid film is mainly determined by the characteristics of interfacial wave,which is rarely investigated right now.How to improve the stability of thin film has become a hot spot.In view of this,an advanced model based on the conventional Volume of Fluid(VOF)model is adopted to investigate the characteristics of interfacial wave in gas-liquid flow by using OpenFOAM,and the mechanism of formation and development of wave is revealed intuitively through numerical study.The effects from gas velocity,surface tension and dynamic viscosity of liquid(three factors)on the wave are studied respectively.It can be found that the gas velocity is critical to the formation and development of wave,and four modes of droplets generation are illustrated in this paper.Besides,a gas vortex near the gas-liquid interface can induce formation of wave easily,so changing the gas vortex state can regulate formation and development of wave.What’s more,the change rules of three factors influencing on the interfacial wave are obtained,and the surface tension has a negative effect on the formation and development of wave only when the surface tension coefficient is above the critical value,whereas the dynamic viscosity has a positive effect in this process.Lastly,the maximum height and maximum slope angle of wave will level off as the gas velocity increases.Meanwhile,the maximum slope angle of wave is usually no more than 38°,no matter what happens to the three factors.
基金the National Natural Science Foundation of China(No.51779143)the Oceanic Interdisciplinary Program of Shanghai Jiao Tong University(No.SL2020ZD101)the Cultivation of Scientific Research Ability of Young Talents of Shanghai Jiao Tong University(No.19X100040072)。
文摘Fluid-structure interaction(FSI)of gas-liquid two-phase fow in the horizontal pipe is investigated numerically in the present study.The volume of fluid model and standard k-e turbulence model are integrated to simulate the typical gas-liquid two-phase fow patterns.First,validation of the numerical model is conducted and the typical fow patterns are consistent with the Baker chart.Then,the FSI framework is established to investigate the dynamic responses of the interaction between the horizontal pipe and gas-liquid two-phase fow.The results show that the dynamic response under stratified fow condition is relatively flat and the maximum pipe deformation and equivalent stress are 1.8 mm and 7.5 MPa respectively.Meanwhile,the dynamic responses induced by slug fow,wave fow and annular fow show obvious periodic fuctuations.Furthermore,the dynamic response characteristics under slug flow condition are maximum;the maximum pipe deformation and equivalent stress can reach 4mm and 17.5 MPa,respectively.The principal direction of total deformation is different under various flow patterns.Therefore,the periodic equivalent stress will form the cyclic impact on the pipe wall and affect the fatigue life of the horizontal pipe.The present study may serve as a reference for FSI simulation under gas-liquid two-phase transport conditions.
基金the support of the Opening Fund of State Key Laboratory of Multiphase Flow in Power Engineering(SKLMF-KF-2102)。
文摘Accurate prediction of the frictional pressure drop is important for the design and operation of subsea oil and gas transporting system considering the length of the pipeline. The applicability of the correlations to pipeline-riser flow needs evaluation since the flow condition in pipeline-riser is quite different from the original data where they were derived from. In the present study, a comprehensive evaluation of 24prevailing correlation in predicting frictional pressure drop is carried out based on experimentally measured data of air-water and air-oil two-phase flows in pipeline-riser. Experiments are performed in a system having different configuration of pipeline-riser with the inclination of the downcomer varied from-2°to-5°to investigated the effect of the elbow on the frictional pressure drop in the riser. The inlet gas velocity ranges from 0.03 to 6.2 m/s, and liquid velocity varies from 0.02 to 1.3 m/s. A total of885 experimental data points including 782 on air-water flows and 103 on air-oil flows are obtained and used to access the prediction ability of the correlations. Comparison of the predicted results with the measured data indicate that a majority of the investigated correlations under-predict the pressure drop on severe slugging. The result of this study highlights the requirement of new method considering the effect of pipe layout on the frictional pressure drop.
基金supported by the National Natural Science Foundation of China Joint Fund Project (Grant/Award Number: U20B6003)National Natural Science Foundation of China (Grant/Award Number: 52304054)。
文摘Understanding fingering, as a challenge to stable displacement during the immiscible flow, has become a crucial phenomenon for geological carbon sequestration, enhanced oil recovery, and groundwater protection. Typically governed by gravity, viscous and capillary forces, these factors lead invasive fluids to occupy pore space irregularly and incompletely. Previous studies have demonstrated capillary numbers,describing the viscous and capillary forces, to quantificationally induce evolution of invasion patterns.While the evolution mechanisms of invasive patterns have not been deeply elucidated under the constant capillary number and three variable parameters including velocity, viscosity, and interfacial tension.Our research employs two horizontal visualization systems and a two-phase laminar flow simulation to investigate the tendency of invasive pattern transition by various parameters at the pore scale. We showed that increasing invasive viscosity or reducing interfacial tension in a homogeneous pore space significantly enhanced sweep efficiency, under constant capillary number. Additionally, in the fingering crossover pattern, the region near the inlet was prone to capillary fingering with multi-directional invasion, while the viscous fingering with unidirectional invasion was more susceptible occurred in the region near the outlet. Furthermore, increasing invasive viscosity or decreasing invasive velocity and interfacial tension promoted the extension of viscous fingering from the outlet to the inlet, presenting that the subsequent invasive fluid flows toward the outlet. In the case of invasive trunk along a unidirectional path, the invasive flow increased exponentially closer to the outlet, resulting in a significant decrease in the width of the invasive interface. Our work holds promising applications for optimizing invasive patterns in heterogeneous porous media.
基金This work is supported by the National Natural Science Foundation of China(No.52104049)the Young Elite Scientist Sponsorship Program by Beijing Association for Science and Technology(No.BYESS2023262)Science Foundation of China University of Petroleum,Beijing(No.2462022BJRC004).
文摘Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinning and shear thickening,polymer convection,diffusion,adsorption retention,inaccessible pore volume and reduced effective permeability.Meanwhile,the flux density and fracture conductivity along the hydraulic fracture are generally non-uniform due to the effects of pressure distribution,formation damage,and proppant breakage.In this paper,we present an oil-water two-phase flow model that captures these complex non-Newtonian and nonlinear behavior,and non-uniform fracture characteristics in fractured polymer flooding.The hydraulic fracture is firstly divided into two parts:high-conductivity fracture near the wellbore and low-conductivity fracture in the far-wellbore section.A hybrid grid system,including perpendicular bisection(PEBI)and Cartesian grid,is applied to discrete the partial differential flow equations,and the local grid refinement method is applied in the near-wellbore region to accurately calculate the pressure distribution and shear rate of polymer solution.The combination of polymer behavior characterizations and numerical flow simulations are applied,resulting in the calculation for the distribution of water saturation,polymer concentration and reservoir pressure.Compared with the polymer flooding well with uniform fracture conductivity,this non-uniform fracture conductivity model exhibits the larger pressure difference,and the shorter bilinear flow period due to the decrease of fracture flow ability in the far-wellbore section.The field case of the fall-off test demonstrates that the proposed method characterizes fracture characteristics more accurately,and yields fracture half-lengths that better match engineering reality,enabling a quantitative segmented characterization of the near-wellbore section with high fracture conductivity and the far-wellbore section with low fracture conductivity.The novelty of this paper is the analysis of pressure performances caused by the fracture dynamics and polymer rheology,as well as an analysis method that derives formation and fracture parameters based on the pressure and its derivative curves.
基金supported by the New Cornerstone Science Foundation through the XPLORER PRIZE and the National Natural Science Foundation of China(Grant No.52088102).
文摘A numerical study based on a two-dimensional two-phase SPH(Smoothed Particle Hydrodynamics)model to analyze the action of water waves on open-type sea access roads is presented.The study is a continuation of the analyses presented by Chen et al.(2022),in which the sea access roads are semi-immersed.In this new configuration,the sea access roads are placed above the still water level,therefore the presence of the air phase becomes a relevant issue in the determination of the wave forces acting on the structures.Indeed,the comparison of wave forces on the open-type sea access roads obtained from the single and two-phase SPH models with the experimental results shows that the latter are in much better agreement.So in the numerical simulations,a two-phaseδ-SPH model is adopted to investigate the dynamical problems.Based on the numerical results,the maximum horizontal and uplifting wave forces acting on the sea access roads are analyzed by considering different wave conditions and geometries of the structures.In particular,the presence of the girder is analyzed and the differences in the wave forces due to the air cushion effects which are created below the structure are highlighted.
基金Supported by the National Natural Science Foundation of China(52374043)Key Program of the National Natural Science Foundation of China(52234003).
文摘Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale gas well considering the two-phase flow of gas and water.The model accounts for the influence of natural fractures and matrix properties on the fracturing process and directly applies post-fracturing formation pressure and water saturation distribution to subsequent well shut-in and production simulation,allowing for a more accurate fracturing-production integrated simulation.The results show that the reservoir physical properties have great impacts on fracture propagation,and the reasonable prediction of formation pressure and reservoir fluid distribution after the fracturing is critical to accurately predict the gas and fluid production of the shale gas wells.Compared with the conventional method,the proposed model can more accurately simulate the water and gas production by considering the impact of fracturing on both matrix pressure and water saturation.The established model is applied to the integrated fracturing-production simulation of practical horizontal shale gas wells.The simulation results are in good agreement with the practical production data,thus verifying the accuracy of the model.
