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.展开更多
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.展开更多
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.展开更多
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 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.展开更多
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.展开更多
Two actual rocks drilled from a typical ultra-deep hydrocarbon reservoir in the Tarim Basin are selected to conduct in-situ stress-loading micro-focus CT scanning experiments.The gray images of rock microstructure at ...Two actual rocks drilled from a typical ultra-deep hydrocarbon reservoir in the Tarim Basin are selected to conduct in-situ stress-loading micro-focus CT scanning experiments.The gray images of rock microstructure at different stress loading stages are obtained.The U-Net fully convolutional neural network is utilized to achieve fine semantic segmentation of rock skeleton,pore space,and microfractures based on CT slice images of deep rocks.The three-dimensional digital rock models of deformed multiscale fractured-porous media at different stress loading stages are thereafter reconstructed,and the equivalent fracture-pore network models are finally extracted to explore the underlying mechanisms of gas-water two-phase flow at the pore-scale.Results indicate that,in the process of insitu stress loading,both the deep rocks have experienced three stages:linear elastic deformation,nonlinear plastic deformation,and shear failure.The micro-mechanical behavior greatly affects the dynamic deformation of rock microstructure and gas-water two-phase flow.In the linear elastic deformation stage,with the increase in in-situ stress,both the deep rocks are gradually compacted,leading to decreases in average pore radius,pore throat ratio,tortuosity,and water-phase relative permeability,while the coordination number nearly remains unchanged.In the plastic deformation stage,the synergistic influence of rock compaction and existence of micro-fractures typically exert a great effect on pore-throat topological properties and gas-water relative permeability.In the shear failure stage,due to the generation and propagation of micro-fractures inside the deep rock,the topological connectivity becomes better,fluid flow paths increase,and flow conductivity is promoted,thus leading to sharp increases in average pore radius and coordination number,rapid decreases in pore throat ratio and tortuosity,as well as remarkable improvement in relative permeability of gas phase and waterphase.展开更多
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.展开更多
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.展开更多
Energetic Semiconductor bridge(ESCB)based on reactive multilayered films(RMFs)has a promising application in the miniature and intelligence of initiator and pyrotechnics device.Understanding the ignition enhancement m...Energetic Semiconductor bridge(ESCB)based on reactive multilayered films(RMFs)has a promising application in the miniature and intelligence of initiator and pyrotechnics device.Understanding the ignition enhancement mechanism of RMFs on semiconductor bridge(SCB)during the ignition process is crucial for the engineering and practical application of advanced initiator and pyrotechnics devices.In this study,a one-dimensional(1D)gas-solid two-phase flow ignition model was established to study the ignition process of ESCB to charge particles based on the reactivity of Al/MoO_(3) RMFs.In order to fully consider the coupled exothermic between the RMFs and the SCB plasma during the ignition process,the heat release of chemical reaction in RMFs was used as an internal heat source in this model.It is found that the exothermal reaction in RMFs improved the ignition performance of SCB.In the process of plasma rapid condensation with heat release,the product of RMFs enhanced the heat transfer process between the gas phase and the solid charge particle,which accelerated the expansion of hot plasma,and heated the solid charge particle as well as gas phase region with low temperature.In addition,it made up for pressure loss in the gas phase.During the plasma dissipation process,the exothermal chemical reaction in RMFs acted as the main heating source to heat the charge particle,making the surface temperature of the charge particle,gas pressure,and gas temperature rise continuously.This result may yield significant advantages in providing a universal ignition model for miniaturized ignition devices.展开更多
In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering comp...In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering complex boundary shapes.Utilizing radial basis function point interpolation,the method approximates shape functions for unknown functions within the nodal influence domain.The shape functions constructed by the aforementioned meshless interpolation method haveδ-function properties,which facilitate the handling of essential aspects like the controlled bottom-hole flow pressure in horizontal wells.Moreover,the meshless method offers greater flexibility and freedom compared to grid cell discretization,making it simpler to discretize complex geometries.A variational principle for the flow control equation group is introduced using a weighted least squares meshless method,and the pressure distribution is solved implicitly.Example results demonstrate that the computational outcomes of the meshless point cloud model,which has a relatively small degree of freedom,are in close agreement with those of the Discrete Fracture Model(DFM)employing refined grid partitioning,with pressure calculation accuracy exceeding 98.2%.Compared to high-resolution grid-based computational methods,the meshless method can achieve a better balance between computational efficiency and accuracy.Additionally,the impact of fracture half-length on the productivity of horizontal wells is discussed.The results indicate that increasing the fracture half-length is an effective strategy for enhancing production from the perspective of cumulative oil production.展开更多
In response to escalating challenges in energy conservation and emission reduction,this study delves into the complexities of heat transfer in two-phase flows and adjustments to combustion processes within coal-fired ...In response to escalating challenges in energy conservation and emission reduction,this study delves into the complexities of heat transfer in two-phase flows and adjustments to combustion processes within coal-fired boilers.Utilizing a fusion of hybrid modeling and automation technologies,we develop soft measurement models for key combustion parameters,such as the net calorific value of coal,flue gas oxygen content,and fly ash carbon content,within theDistributedControl System(DCS).Validated with performance test data,thesemodels exhibit controlled root mean square error(RMSE)and maximum absolute error(MAXE)values,both within the range of 0.203.Integrated into their respective automatic control systems,thesemodels optimize two-phase flow heat transfer,finetune combustion conditions,and mitigate incomplete combustion.Furthermore,this paper conducts an in-depth exploration of the generationmechanismof nitrogen oxides(NO_(x))and low oxygen emission reduction technology in coal-fired boilers,demonstrating a substantial reduction in furnace exit NO_(x) generation by 30%to 40%and the power supply coal consumption decreased by 1.62 g/(kW h).The research outcomes highlight the model’s rapid responsiveness,enabling prompt reflection of transient variations in various economic indicator parameters.This provides a more effective means for real-time monitoring of crucial variables in coal-fired boilers and facilitates timely combustion adjustments,underscoring notable achievements in boiler combustion.The research not only provides valuable and practical insights into the intricacies of two-phase flow heat transfer and heat exchange but also establishes a pioneering methodology for tackling industry challenges.展开更多
According to a mathematical model for dense two-phase flows presented in the previous pape[1],a dense two-phase flow in a vertical pipeline is analytically solved, and the analytic expressions of velocity of each cont...According to a mathematical model for dense two-phase flows presented in the previous pape[1],a dense two-phase flow in a vertical pipeline is analytically solved, and the analytic expressions of velocity of each continuous phase and dispersed phase are respectively derived. The results show that when the drag force between twophasesdepends linearly on their relative velocity, the relative velocity profile in the pipeline coincides with Darcy's law except for the thin layer region near the pipeline wall, and that the theoretical assumptions in the dense two-phase flow theory mentioned are reasonable.展开更多
This paper presents a mass and momentum conservative semi-implicit finite volume(FV)scheme for complex non-hydrostatic free surface flows,interacting with moving solid obstacles.A simplified incompressible Baer-Nunzia...This paper presents a mass and momentum conservative semi-implicit finite volume(FV)scheme for complex non-hydrostatic free surface flows,interacting with moving solid obstacles.A simplified incompressible Baer-Nunziato type model is considered for two-phase flows containing a liquid phase,a solid phase,and the surrounding void.According to the so-called diffuse interface approach,the different phases and consequently the void are described by means of a scalar volume fraction function for each phase.In our numerical scheme,the dynamics of the liquid phase and the motion of the solid are decoupled.The solid is assumed to be a moving rigid body,whose motion is prescribed.Only after the advection of the solid volume fraction,the dynamics of the liquid phase is considered.As usual in semi-implicit schemes,we employ staggered Cartesian control volumes and treat the nonlinear convective terms explicitly,while the pressure terms are treated implicitly.The non-conservative products arising in the transport equation for the solid volume fraction are treated by a path-conservative approach.The resulting semi-implicit FV discretization of the mass and momentum equations leads to a mildly nonlinear system for the pressure which can be efficiently solved with a nested Newton-type technique.The time step size is only limited by the velocities of the two phases contained in the domain,and not by the gravity wave speed nor by the stiff algebraic relaxation source term,which requires an implicit discretization.The resulting semi-implicit algorithm is first validated on a set of classical incompressible Navier-Stokes test problems and later also adds a fixed and moving solid phase.展开更多
Oil reservoirs with low permeability and porosity that are in the middle and late exploitation periods in China＇s onshore oil fields are mostly in the high-water-cut production stage.This stage is associated with sev...Oil reservoirs with low permeability and porosity that are in the middle and late exploitation periods in China＇s onshore oil fields are mostly in the high-water-cut production stage.This stage is associated with severely non-uniform local-velocity flow profiles and dispersed-phase concentration（of oil droplets） in oil-water two-phase flow,which makes it difficult to measure water holdup in oil wells.In this study,we use an ultrasonic method based on a transmission-type sensor in oil-water two-phase flow to measure water holdup in lowvelocity and high water-cut conditions.First,we optimize the excitation frequency of the ultrasonic sensor by calculating the sensitivity of the ultrasonic field using the finite element method for multiphysics coupling.Then we calculate the change trend of sound pressure level attenuation ratio with the increase in oil holdup to verify the feasibility of the employed diameter for the ultrasonic sensor.Based on the results,we then investigate the effects of oildroplet diameter and distribution on the ultrasonic field.To further understand the measurement characteristics of the ultrasonic sensor,we perform a flow loop test on vertical upward oilwater two-phase flow and measure the responses of the optimized ultrasonic sensor.The results show that the ultrasonic sensor yields poor resolution for a dispersed oil slug in water flow（D OS/W flow）,but the resolution is favorable for dispersed oil in water flow（D O/W flow） and very fine dispersed oil in water flow（VFD O/W flow）.This research demonstrates the potential application of a pulsed-transmission ultrasonic method for measuring the fraction of individual components in oil-water two-phase flow with a low mixture velocity and high water cut.展开更多
基金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.
基金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.
基金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.
文摘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.
文摘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 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(No.52174043)the Beijing Natural Science Foundation(No.3242019)+1 种基金the CNPC Innovation Foundation(No.2022DQ02-0208)the State Key Laboratory of Deep Oil and Gas(No.SKLD0G2024-KFZD-06).
文摘Two actual rocks drilled from a typical ultra-deep hydrocarbon reservoir in the Tarim Basin are selected to conduct in-situ stress-loading micro-focus CT scanning experiments.The gray images of rock microstructure at different stress loading stages are obtained.The U-Net fully convolutional neural network is utilized to achieve fine semantic segmentation of rock skeleton,pore space,and microfractures based on CT slice images of deep rocks.The three-dimensional digital rock models of deformed multiscale fractured-porous media at different stress loading stages are thereafter reconstructed,and the equivalent fracture-pore network models are finally extracted to explore the underlying mechanisms of gas-water two-phase flow at the pore-scale.Results indicate that,in the process of insitu stress loading,both the deep rocks have experienced three stages:linear elastic deformation,nonlinear plastic deformation,and shear failure.The micro-mechanical behavior greatly affects the dynamic deformation of rock microstructure and gas-water two-phase flow.In the linear elastic deformation stage,with the increase in in-situ stress,both the deep rocks are gradually compacted,leading to decreases in average pore radius,pore throat ratio,tortuosity,and water-phase relative permeability,while the coordination number nearly remains unchanged.In the plastic deformation stage,the synergistic influence of rock compaction and existence of micro-fractures typically exert a great effect on pore-throat topological properties and gas-water relative permeability.In the shear failure stage,due to the generation and propagation of micro-fractures inside the deep rock,the topological connectivity becomes better,fluid flow paths increase,and flow conductivity is promoted,thus leading to sharp increases in average pore radius and coordination number,rapid decreases in pore throat ratio and tortuosity,as well as remarkable improvement in relative permeability of gas phase and waterphase.
文摘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 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.
基金supported by the National Natural Science Foundation of China(Grant Nos.22275092,52102107 and 52372084)the Fundamental Research Funds for the Central Universities(Grant No.30923010920)。
文摘Energetic Semiconductor bridge(ESCB)based on reactive multilayered films(RMFs)has a promising application in the miniature and intelligence of initiator and pyrotechnics device.Understanding the ignition enhancement mechanism of RMFs on semiconductor bridge(SCB)during the ignition process is crucial for the engineering and practical application of advanced initiator and pyrotechnics devices.In this study,a one-dimensional(1D)gas-solid two-phase flow ignition model was established to study the ignition process of ESCB to charge particles based on the reactivity of Al/MoO_(3) RMFs.In order to fully consider the coupled exothermic between the RMFs and the SCB plasma during the ignition process,the heat release of chemical reaction in RMFs was used as an internal heat source in this model.It is found that the exothermal reaction in RMFs improved the ignition performance of SCB.In the process of plasma rapid condensation with heat release,the product of RMFs enhanced the heat transfer process between the gas phase and the solid charge particle,which accelerated the expansion of hot plasma,and heated the solid charge particle as well as gas phase region with low temperature.In addition,it made up for pressure loss in the gas phase.During the plasma dissipation process,the exothermal chemical reaction in RMFs acted as the main heating source to heat the charge particle,making the surface temperature of the charge particle,gas pressure,and gas temperature rise continuously.This result may yield significant advantages in providing a universal ignition model for miniaturized ignition devices.
文摘In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering complex boundary shapes.Utilizing radial basis function point interpolation,the method approximates shape functions for unknown functions within the nodal influence domain.The shape functions constructed by the aforementioned meshless interpolation method haveδ-function properties,which facilitate the handling of essential aspects like the controlled bottom-hole flow pressure in horizontal wells.Moreover,the meshless method offers greater flexibility and freedom compared to grid cell discretization,making it simpler to discretize complex geometries.A variational principle for the flow control equation group is introduced using a weighted least squares meshless method,and the pressure distribution is solved implicitly.Example results demonstrate that the computational outcomes of the meshless point cloud model,which has a relatively small degree of freedom,are in close agreement with those of the Discrete Fracture Model(DFM)employing refined grid partitioning,with pressure calculation accuracy exceeding 98.2%.Compared to high-resolution grid-based computational methods,the meshless method can achieve a better balance between computational efficiency and accuracy.Additionally,the impact of fracture half-length on the productivity of horizontal wells is discussed.The results indicate that increasing the fracture half-length is an effective strategy for enhancing production from the perspective of cumulative oil production.
文摘In response to escalating challenges in energy conservation and emission reduction,this study delves into the complexities of heat transfer in two-phase flows and adjustments to combustion processes within coal-fired boilers.Utilizing a fusion of hybrid modeling and automation technologies,we develop soft measurement models for key combustion parameters,such as the net calorific value of coal,flue gas oxygen content,and fly ash carbon content,within theDistributedControl System(DCS).Validated with performance test data,thesemodels exhibit controlled root mean square error(RMSE)and maximum absolute error(MAXE)values,both within the range of 0.203.Integrated into their respective automatic control systems,thesemodels optimize two-phase flow heat transfer,finetune combustion conditions,and mitigate incomplete combustion.Furthermore,this paper conducts an in-depth exploration of the generationmechanismof nitrogen oxides(NO_(x))and low oxygen emission reduction technology in coal-fired boilers,demonstrating a substantial reduction in furnace exit NO_(x) generation by 30%to 40%and the power supply coal consumption decreased by 1.62 g/(kW h).The research outcomes highlight the model’s rapid responsiveness,enabling prompt reflection of transient variations in various economic indicator parameters.This provides a more effective means for real-time monitoring of crucial variables in coal-fired boilers and facilitates timely combustion adjustments,underscoring notable achievements in boiler combustion.The research not only provides valuable and practical insights into the intricacies of two-phase flow heat transfer and heat exchange but also establishes a pioneering methodology for tackling industry challenges.
文摘According to a mathematical model for dense two-phase flows presented in the previous pape[1],a dense two-phase flow in a vertical pipeline is analytically solved, and the analytic expressions of velocity of each continuous phase and dispersed phase are respectively derived. The results show that when the drag force between twophasesdepends linearly on their relative velocity, the relative velocity profile in the pipeline coincides with Darcy's law except for the thin layer region near the pipeline wall, and that the theoretical assumptions in the dense two-phase flow theory mentioned are reasonable.
基金funded by the Italian Ministry of Education,University and Research(MIUR)in the frame of the Departments of Excellence Initiative 2018-2027 attributed to DICAM of the University of Trento(grant L.232/2016)in the frame of the PRIN 2017 project Innovative numerical methods for evolutionary partial differential equations and applications,the PRIN 2022 project High order structure-preserving semi-implicit schemes for hyperbolic equations.D.is member of INdAM GNCS and was also co-funded by the European Union NextGenerationEU(PNRR,Spoke 7 CN HPC).Views and opinions expressed are however those of the author(s)only and do not necessarily reflect those of the European Union or the European Research Council.Neither the European Union nor the granting authority can be held responsible for them.
文摘This paper presents a mass and momentum conservative semi-implicit finite volume(FV)scheme for complex non-hydrostatic free surface flows,interacting with moving solid obstacles.A simplified incompressible Baer-Nunziato type model is considered for two-phase flows containing a liquid phase,a solid phase,and the surrounding void.According to the so-called diffuse interface approach,the different phases and consequently the void are described by means of a scalar volume fraction function for each phase.In our numerical scheme,the dynamics of the liquid phase and the motion of the solid are decoupled.The solid is assumed to be a moving rigid body,whose motion is prescribed.Only after the advection of the solid volume fraction,the dynamics of the liquid phase is considered.As usual in semi-implicit schemes,we employ staggered Cartesian control volumes and treat the nonlinear convective terms explicitly,while the pressure terms are treated implicitly.The non-conservative products arising in the transport equation for the solid volume fraction are treated by a path-conservative approach.The resulting semi-implicit FV discretization of the mass and momentum equations leads to a mildly nonlinear system for the pressure which can be efficiently solved with a nested Newton-type technique.The time step size is only limited by the velocities of the two phases contained in the domain,and not by the gravity wave speed nor by the stiff algebraic relaxation source term,which requires an implicit discretization.The resulting semi-implicit algorithm is first validated on a set of classical incompressible Navier-Stokes test problems and later also adds a fixed and moving solid phase.
基金supported by the National Natural Science Foundation of China(Nos.51527805,11572220 and 41174109)
文摘Oil reservoirs with low permeability and porosity that are in the middle and late exploitation periods in China＇s onshore oil fields are mostly in the high-water-cut production stage.This stage is associated with severely non-uniform local-velocity flow profiles and dispersed-phase concentration（of oil droplets） in oil-water two-phase flow,which makes it difficult to measure water holdup in oil wells.In this study,we use an ultrasonic method based on a transmission-type sensor in oil-water two-phase flow to measure water holdup in lowvelocity and high water-cut conditions.First,we optimize the excitation frequency of the ultrasonic sensor by calculating the sensitivity of the ultrasonic field using the finite element method for multiphysics coupling.Then we calculate the change trend of sound pressure level attenuation ratio with the increase in oil holdup to verify the feasibility of the employed diameter for the ultrasonic sensor.Based on the results,we then investigate the effects of oildroplet diameter and distribution on the ultrasonic field.To further understand the measurement characteristics of the ultrasonic sensor,we perform a flow loop test on vertical upward oilwater two-phase flow and measure the responses of the optimized ultrasonic sensor.The results show that the ultrasonic sensor yields poor resolution for a dispersed oil slug in water flow（D OS/W flow）,but the resolution is favorable for dispersed oil in water flow（D O/W flow） and very fine dispersed oil in water flow（VFD O/W flow）.This research demonstrates the potential application of a pulsed-transmission ultrasonic method for measuring the fraction of individual components in oil-water two-phase flow with a low mixture velocity and high water cut.
