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.展开更多
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.展开更多
The two-phase flow models are commonly used in industrial applications, such as nuclear, power, chemical-process, oil-and-gas, cryogenics, bio-medical, micro-technology and so on. This is a survey paper on the study o...The two-phase flow models are commonly used in industrial applications, such as nuclear, power, chemical-process, oil-and-gas, cryogenics, bio-medical, micro-technology and so on. This is a survey paper on the study of compressible nonconservative two-fluid model, drift-flux model and viscous liquid-gas two-phase flow model. We give the research developments of these three two-phase flow models, respectively. In the last part, we give some open problems about the above models.展开更多
Sediment deposition in the pumping station has a huge negative impact on unit operation.The three-dimensional CFD method has been used to simulate inlet structure flow in pumping station based on the Eulerian solid- l...Sediment deposition in the pumping station has a huge negative impact on unit operation.The three-dimensional CFD method has been used to simulate inlet structure flow in pumping station based on the Eulerian solid- liquid two-phase flow model. The numerical results of the preliminary scheme show that sediment deposition occurs in the forebay of pumping station because of poor flow pattern therein. In order to improve hydraulic configuration in the forebay,one modified measure reconstructs water diversion weir shape,and another measure sets a water retaining sill in the approach channel. The simulation results of the modified scheme prove that back flow in the forebay has been eliminated and the sediment deposition region has also been reduced.展开更多
Rainfall infiltration on a soil slope is usually an unsaturated seepage process that can be described by a water-air two-phase flow model.The effect of pore air pressure on rainfall infiltration has been widely recogn...Rainfall infiltration on a soil slope is usually an unsaturated seepage process that can be described by a water-air two-phase flow model.The effect of pore air pressure on rainfall infiltration has been widely recognized and validated by means of numerical simulations and laboratory experiments.However,whether a slope can actually seal pore air continues to be debated by researchers.In this study,a water-air two-phase flow model is used to simulate the rainfall infiltration process on a soil slope,and a field experiment is conducted to realistically test the sealing conditions of a slope.According to the numerical simulation,the areas of water and air flow in and out on the slope surface are relatively stable and can be classified as the“inhalation zone”and“overflow zone”,respectively.Intermittent rainfall on the soil slope has an amplifying effect on pore air pressure because rainfall intensity is usually at the millimeter level,and it causes pore air pressure to reach the cm level.A field experiment was performed to determine whether a slope can realistically seal pore air and subsequently verify the regularity of rainfall infiltration.Air pressure sensors were buried in the slope to monitor the pore air pressures during the rainfall process.The monitoring results show that the pore air pressure in the slope changed,which indicates that the slope can seal air.Moreover,the amplification effects of intermittent rainfall on pore air pressure were observed for natural rainfall,which agrees well with the numerical simulation results.展开更多
Movement and growth of dendrites are common phenomena during solidification.To numerically investigate these phenomena,two-phase flow model is employed to formulate the FSI(fluid-structure interaction)problem during d...Movement and growth of dendrites are common phenomena during solidification.To numerically investigate these phenomena,two-phase flow model is employed to formulate the FSI(fluid-structure interaction)problem during dendritic solidification.In this model,solid is assumed to have huge viscosity to maintain its own shape and an exponential expression is constructed to describe variable viscosity across s-l(solid-liquid)interface.With an effective preconditioner for saddle point structure,we build a N-S(Navier-Stokes)solver robust to tremendous viscosity ratio(as large as 10^(10))between solid and liquid.Polycrystalline solidification is computed by vector-valued phase field model,which is computationally convenient to handle contact between dendrites.Locations of dendrites are updated by solving advection equations.Orientation change due to dendrite's rotation has been considered as well.Calculation is accelerated by two-level time stepping scheme,adaptive mesh refinement,and parallel computation.Settlement and growth of a single dendrite and multiple dendrites in Al-Cu alloy were simulated,showing the availability of the provided model to handle anisotropic growth,motion and impingement of dendrites.This study lays foundation to simulate solidification coupled with deformation in the future.展开更多
In the industrial process of producing the strong phosphoric acid(SPA),clarification of the solution is essential to the ultimate product.However,the large viscosity of sediment and the induced interface interaction r...In the industrial process of producing the strong phosphoric acid(SPA),clarification of the solution is essential to the ultimate product.However,the large viscosity of sediment and the induced interface interaction result in difficulties when the SPA is clarified.CFD numerical methodology was applied to simulate internal flow field and performance of the low speed scraper based on Mixture solidliquid two-phase flow model.Sediment deposition was generated by loading solid particles at the bottom of clarifying vessel.The moving mesh and RNG k-εmodel were used to simulate the rotational turbulent flow in clarifying tank.Variables studied,amongst others,were the scraper rotation speed and the mounting height,which could affect the solid suspension height.Features of flow field and solid volume fraction distribution in computational domain were presented and analyzed.The numerical reports of the scraper torque and velocities of inlet and outlet filed were obtained.It seems the torque value of rotatio-nal axis and particle suspending height augment with an increasing rotating speed.Meanwhile,a high revolving speed is good for the deposition discharge.The particle fraction distribution in meridional surface and horizontal surface at fixed rotation speed were analyzed to determine the corresponding optimal installation height.The simulating results reflect the flow field is marginally stirred by the scraper and proper working parameters are obtained,in which case the comprehensive properties of the scraper and the clarifying tank are superior.展开更多
The outflow problem for the viscous two-phase flow model in a half line is investigated in the present paper.The existence and uniqueness of the stationary solution is shown for both supersonic state and sonic state a...The outflow problem for the viscous two-phase flow model in a half line is investigated in the present paper.The existence and uniqueness of the stationary solution is shown for both supersonic state and sonic state at spatial far field,and the nonlinear time stability of the stationary solution is also established in the weighted Sobolev space with either the exponential time decay rate for supersonic flow or the algebraic time decay rate for sonic flow.展开更多
Mathematical modeling of pulsating heat pipes through ‘first’ principles is a contemporary problem which remains quite elusive. Simplifications and assumptions made in all the modeling approaches developed so far re...Mathematical modeling of pulsating heat pipes through ‘first’ principles is a contemporary problem which remains quite elusive. Simplifications and assumptions made in all the modeling approaches developed so far render them unsuitable for engineering design. In this paper, a more realistic modeling scheme is presented which provides considerable try for thought toward the next progressive step. At high enough heat flux level, closed loop pulsating heat pipes experience a bulk internal unidirectional fluid circulation. Under such a condition, conventional two-phase flow modeling in capillary tubes may be applied. This has been attempted for single-loop PHPs. A homogeneous model and a separated two-fluid flow model based on simultaneous conservation of mass, momentum and energy, have been developed for an equivalent ‘open flow’ system. The model allows prediction of two-phase flow parameters in each sub-section of the device thereby providing important insights into its operation. The concept of ‘void fraction constraint’ in pulsating heat pipe operation is introduced and its relevance to future modeling attempts is outlined.展开更多
In this paper, we study the global existence and uniqueness of strong solutions for the Baer-Nunziato two-phase flow model in a bounded domain with a no-slip boundary. The global existence and uniqueness of strong sol...In this paper, we study the global existence and uniqueness of strong solutions for the Baer-Nunziato two-phase flow model in a bounded domain with a no-slip boundary. The global existence and uniqueness of strong solutions are obtained when the initial value is near the equilibrium state in H<sup>2</sup> (Ω). Furthermore, the exponential convergence rates of the pressure and velocity are also proved by delicate energy methods.展开更多
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.展开更多
Polymer flooding is an important means of improving oil recovery and is widely used in Daqing,Xinjiang,and Shengli oilfields,China.Different from conventional injection media such as water and gas,viscoelastic polymer...Polymer flooding is an important means of improving oil recovery and is widely used in Daqing,Xinjiang,and Shengli oilfields,China.Different from conventional injection media such as water and gas,viscoelastic polymer solutions exhibit non-Newtonian and nonlinear flow behavior including shear thinning and shear thickening,polymer convection,diffusion,adsorption,retention,inaccessible pore volume,and reduced effective permeability.However,available well test model of polymer flooding wells generally simplifies these characteristics on pressure transient response,which may lead to inaccurate results.This work proposes a novel two-phase numerical well test model to better describe the polymer viscoelasticity and nonlinear flow behavior.Different influence factors that related to near-well blockage during polymer flooding process,including the degree of blockage(inner zone permeability),the extent of blockage(composite radius),and polymer flooding front radius are explored to investigate these impacts on bottom hole pressure responses.Results show that polymer viscoelasticity has a significant impact on the transitional flow segment of type curves,and the effects of near-well formation blockage and polymer concentration distribution on well test curves are very similar.Thus,to accurately interpret the degree of near-well blockage in injection wells,it is essential to first eliminate the influence of polymer viscoelasticity.Finally,a field case is comprehensively analyzed and discussed to illustrate the applicability of the proposed model.展开更多
Clayey-silt natural gas hydrate reservoirs in the South China Sea exhibit loose and unconsolidated structures, heterogeneous pore structures, high clay mineral contents, and strong hydrophilicity. These characteristic...Clayey-silt natural gas hydrate reservoirs in the South China Sea exhibit loose and unconsolidated structures, heterogeneous pore structures, high clay mineral contents, and strong hydrophilicity. These characteristics complicate the gas-water two-phase flow process in porous media following hydrate decomposition, posing challenges for efficient development. This study examines the transport response of clayey-silt reservoir samples from the Shenhu area using gas-water two-phase flow experiments and CT scanning to explore changes in pore structure, gas-water distribution, and relative permeability under varying flow conditions. The results indicate that pore heterogeneity significantly influences flow characteristics. Gas preferentially displaces water in larger pores, forming fracture-like pores, which serve as preferential flow channels for gas migration. The preferential flow channels enhance gas-phase permeability up to 19 times that of the water phase when fluid pressures exceed total stresses. However,small pores retain liquid, leading to a high residual water saturation of 0.561. CT imaging reveals that these hydro-fractures improve gas permeability but also confine gas flow to specific channels. Pore network analysis shows that gas injection expands the pore-throat network, enhancing connectivity and forming fracture-like pores. Residual water remains trapped in smaller pores and throats, while structural changes, including new fractures, improve gas flow pathways and overall connectivity. Relative permeability curves demonstrate a narrow gas-water cocurrent-flow zone, a right-shifted iso-permeability point and high reservoir capillary pressure, indicating a strong "water-blocking" effect. The findings suggest that optimizing reservoir stimulation techniques to enhance fracture formation, reduce residual water saturation, and improve gas flow capacity is critical for efficient hydrate reservoir development.展开更多
Shale gas production involves complex gas-water two-phase flow,with flow patterns in proppant-filled fractures playing a critical role in determining production efficiency.In this study,3D geometric models of 40/70 me...Shale gas production involves complex gas-water two-phase flow,with flow patterns in proppant-filled fractures playing a critical role in determining production efficiency.In this study,3D geometric models of 40/70 mesh ceramic particles and quartz sand proppant clusters were elaborated using computed tomography(CT)scanning.These models were used to develop a numerical simulation framework based on the lattice Boltzmann method(LBM),enabling the investigation of gas-water flow behavior within proppant-filled fractures under varying driving forces and surface tensions.Simulation results at a closure pressure of 15 MPa have revealed that ceramic particles exhibit a simpler and more porous internal structure than quartz sand of the same size.Under identical flow conditions,ceramic proppants demonstrate higher fluid replacement efficiency.Replacement efficiency increases with higher porosity,greater driving force,and lower surface tension.Furthermore,fluid displacement is strongly influenced by pore geometry:flow is faster in straighter and wider channels,with preferential movement through larger pores forming dominant flow paths.The replacement velocity exhibits a characteristic time evolution,initially rapid,then gradually decreasing,correlating positively with the development of these dominant channels.展开更多
Cleat serves as the primary flow pathway for coalbed methane(CBM)and water.However,few studies consider the impact of local contact on two-phase flow within cleats.A visual generalized model of endogenous cleats was c...Cleat serves as the primary flow pathway for coalbed methane(CBM)and water.However,few studies consider the impact of local contact on two-phase flow within cleats.A visual generalized model of endogenous cleats was constructed based on microfluidics.A microscopic and mesoscopic observation technique was proposed to simultaneously capture gas-liquid interface morphology of pores and throat and the two-phase flow characteristics in entire cleat system.The local contact characteristics of cleats reduced absolute permeability,which resulted in a sharp increase in the starting pressure.The reduced gas flow capacity narrowed the co-infiltration area and decreased water saturation at the isotonic point in a hydrophilic environment.The increased local contact area of cleats weakened gas phase flow capacity and narrowed the co-infiltration area.Jumping events occurred in methane-water flow due to altered porosity caused by local contact in cleats.The distribution of residual phases changed the jumping direction on the micro-scale as well as the dominant channel on the mesoscale.Besides,jumping events caused additional energy dissipation,which was ignored in traditional two-phase flow models.This might contribute to the overestimation of relative permeability.The work provides new methods and insights for investigating unsaturated flow in complex porous media.展开更多
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.展开更多
A two-phase flow model accelerated by graphical processing unit(GPU)is developed to solve fluid-solid interaction(FSI)using the sharp-interface immersed boundary method(IBM).This model solves the incompressible Navier...A two-phase flow model accelerated by graphical processing unit(GPU)is developed to solve fluid-solid interaction(FSI)using the sharp-interface immersed boundary method(IBM).This model solves the incompressible Navier-Stokes equations using the projection-based fractional step method in a fixed staggered Cartesian grid system.A volume of fluid(VOF)method with second-order accuracy is employed to trace the free surface.To represent the intricate surface geometry,the structure is discretized using the unstructured triangle mesh.Additionally,a ray tracing method is employed to classify fluid and solid points.A high-order stable scheme has been introduced to reconstruct the local velocity at interface points.Three FSI problems,including wave evolution around a breakwater,interaction between a periodic wave train and a moving float,and a 3-D moving object interacting with the free surface,were investigated to validate the accuracy and stability of the proposed model.The numerical results are in good agreement with the experimental data.Additionally,we evaluated the computational performance of the proposed GPU-based model.The GPU-based model achieved a 42.29 times speedup compared with the single-core CPU-based model in the three-dimension test.Additionally,the results regarding the time cost of each code section indicate that achieving more significant acceleration is associated with solving the turbulence,advection,and diffusion terms,while solving the pressure Poisson equation(PPE)saves the most time.Furthermore,the impact of grid number on computational efficiency indicates that as Fluid-solid interaction(FSI)immersed boundary method(IBM)graphical processing unit(GPU)two-phase flow moving rigid bodythe number of grids increases,the GPU-based model outperforms the multi-core CPU-based model.展开更多
Due to the coupling between the hydrodynamic equation and the phase-field equation in two-phase incompressible flows,it is desirable to develop efficient and high-order accurate numerical schemes that can decouple the...Due to the coupling between the hydrodynamic equation and the phase-field equation in two-phase incompressible flows,it is desirable to develop efficient and high-order accurate numerical schemes that can decouple these two equations.One popular and efficient strategy is to add an explicit stabilizing term to the convective velocity in the phase-field equation to decouple them.The resulting schemes are only first-order accurate in time,and it seems extremely difficult to generalize the idea of stabilization to the second-order or higher version.In this paper,we employ the spectral deferred correction method to improve the temporal accuracy,based on the first-order decoupled and energy-stable scheme constructed by the stabilization idea.The novelty lies in how the decoupling and linear implicit properties are maintained to improve the efficiency.Within the framework of the spatially discretized local discontinuous Galerkin method,the resulting numerical schemes are fully decoupled,efficient,and high-order accurate in both time and space.Numerical experiments are performed to validate the high-order accuracy and efficiency of the methods for solving phase-field models of two-phase incompressible flows.展开更多
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.展开更多
基金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(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.
基金supported by the National Natural Science Foundation of China(11722104,11671150)supported by the National Natural Science Foundation of China(11571280,11331005)+3 种基金supported by the National Natural Science Foundation of China(11331005,11771150)by GDUPS(2016)the Fundamental Research Funds for the Central Universities of China(D2172260)FANEDD No.201315
文摘The two-phase flow models are commonly used in industrial applications, such as nuclear, power, chemical-process, oil-and-gas, cryogenics, bio-medical, micro-technology and so on. This is a survey paper on the study of compressible nonconservative two-fluid model, drift-flux model and viscous liquid-gas two-phase flow model. We give the research developments of these three two-phase flow models, respectively. In the last part, we give some open problems about the above models.
基金Chinese National Foundation of Natural Science-Key Projects(51339005)
文摘Sediment deposition in the pumping station has a huge negative impact on unit operation.The three-dimensional CFD method has been used to simulate inlet structure flow in pumping station based on the Eulerian solid- liquid two-phase flow model. The numerical results of the preliminary scheme show that sediment deposition occurs in the forebay of pumping station because of poor flow pattern therein. In order to improve hydraulic configuration in the forebay,one modified measure reconstructs water diversion weir shape,and another measure sets a water retaining sill in the approach channel. The simulation results of the modified scheme prove that back flow in the forebay has been eliminated and the sediment deposition region has also been reduced.
基金sponsored by The National Natural Science Foundation of China(Grant Nos.51939004 and 51279090)The National Key Research and Development Program of China(2017YFC1501100)the Hubei Key Laboratory of Construction and Management in Hydropower Engineering(2020KSD11).
文摘Rainfall infiltration on a soil slope is usually an unsaturated seepage process that can be described by a water-air two-phase flow model.The effect of pore air pressure on rainfall infiltration has been widely recognized and validated by means of numerical simulations and laboratory experiments.However,whether a slope can actually seal pore air continues to be debated by researchers.In this study,a water-air two-phase flow model is used to simulate the rainfall infiltration process on a soil slope,and a field experiment is conducted to realistically test the sealing conditions of a slope.According to the numerical simulation,the areas of water and air flow in and out on the slope surface are relatively stable and can be classified as the“inhalation zone”and“overflow zone”,respectively.Intermittent rainfall on the soil slope has an amplifying effect on pore air pressure because rainfall intensity is usually at the millimeter level,and it causes pore air pressure to reach the cm level.A field experiment was performed to determine whether a slope can realistically seal pore air and subsequently verify the regularity of rainfall infiltration.Air pressure sensors were buried in the slope to monitor the pore air pressures during the rainfall process.The monitoring results show that the pore air pressure in the slope changed,which indicates that the slope can seal air.Moreover,the amplification effects of intermittent rainfall on pore air pressure were observed for natural rainfall,which agrees well with the numerical simulation results.
基金financially supported by the National Key Research and Development Program(No.2018YFA0702900)the National Natural Science Foundation of China(Nos.51774265 ,51701225)+3 种基金the National Science and Technology Major Project of China(No.2019ZX06004010)the Key Program of the Chinese Academy of Sciences(No.ZDRW-CN-2017-1)the Program of CAS Interdisciplinary Innovation Teamand Youth Innovation Promotion Association,CAS。
文摘Movement and growth of dendrites are common phenomena during solidification.To numerically investigate these phenomena,two-phase flow model is employed to formulate the FSI(fluid-structure interaction)problem during dendritic solidification.In this model,solid is assumed to have huge viscosity to maintain its own shape and an exponential expression is constructed to describe variable viscosity across s-l(solid-liquid)interface.With an effective preconditioner for saddle point structure,we build a N-S(Navier-Stokes)solver robust to tremendous viscosity ratio(as large as 10^(10))between solid and liquid.Polycrystalline solidification is computed by vector-valued phase field model,which is computationally convenient to handle contact between dendrites.Locations of dendrites are updated by solving advection equations.Orientation change due to dendrite's rotation has been considered as well.Calculation is accelerated by two-level time stepping scheme,adaptive mesh refinement,and parallel computation.Settlement and growth of a single dendrite and multiple dendrites in Al-Cu alloy were simulated,showing the availability of the provided model to handle anisotropic growth,motion and impingement of dendrites.This study lays foundation to simulate solidification coupled with deformation in the future.
基金Graduate Research and Innovation Program in Jiangsu Province(KYZZ16_0286)
文摘In the industrial process of producing the strong phosphoric acid(SPA),clarification of the solution is essential to the ultimate product.However,the large viscosity of sediment and the induced interface interaction result in difficulties when the SPA is clarified.CFD numerical methodology was applied to simulate internal flow field and performance of the low speed scraper based on Mixture solidliquid two-phase flow model.Sediment deposition was generated by loading solid particles at the bottom of clarifying vessel.The moving mesh and RNG k-εmodel were used to simulate the rotational turbulent flow in clarifying tank.Variables studied,amongst others,were the scraper rotation speed and the mounting height,which could affect the solid suspension height.Features of flow field and solid volume fraction distribution in computational domain were presented and analyzed.The numerical reports of the scraper torque and velocities of inlet and outlet filed were obtained.It seems the torque value of rotatio-nal axis and particle suspending height augment with an increasing rotating speed.Meanwhile,a high revolving speed is good for the deposition discharge.The particle fraction distribution in meridional surface and horizontal surface at fixed rotation speed were analyzed to determine the corresponding optimal installation height.The simulating results reflect the flow field is marginally stirred by the scraper and proper working parameters are obtained,in which case the comprehensive properties of the scraper and the clarifying tank are superior.
基金the paper is supported by the National Natural Science Foundation of China(Nos.11871047,11671384,11931010)the key research project of Academy for Multidisciplinary Studies,Capital Normal University,and by the Capacity Building for Sci-Tech Innovation-Fundamental Scientific Research Funds(No.007/20530290068).
文摘The outflow problem for the viscous two-phase flow model in a half line is investigated in the present paper.The existence and uniqueness of the stationary solution is shown for both supersonic state and sonic state at spatial far field,and the nonlinear time stability of the stationary solution is also established in the weighted Sobolev space with either the exponential time decay rate for supersonic flow or the algebraic time decay rate for sonic flow.
基金German National Science Foundation (DFG)(No. GR-412/22)
文摘Mathematical modeling of pulsating heat pipes through ‘first’ principles is a contemporary problem which remains quite elusive. Simplifications and assumptions made in all the modeling approaches developed so far render them unsuitable for engineering design. In this paper, a more realistic modeling scheme is presented which provides considerable try for thought toward the next progressive step. At high enough heat flux level, closed loop pulsating heat pipes experience a bulk internal unidirectional fluid circulation. Under such a condition, conventional two-phase flow modeling in capillary tubes may be applied. This has been attempted for single-loop PHPs. A homogeneous model and a separated two-fluid flow model based on simultaneous conservation of mass, momentum and energy, have been developed for an equivalent ‘open flow’ system. The model allows prediction of two-phase flow parameters in each sub-section of the device thereby providing important insights into its operation. The concept of ‘void fraction constraint’ in pulsating heat pipe operation is introduced and its relevance to future modeling attempts is outlined.
文摘In this paper, we study the global existence and uniqueness of strong solutions for the Baer-Nunziato two-phase flow model in a bounded domain with a no-slip boundary. The global existence and uniqueness of strong solutions are obtained when the initial value is near the equilibrium state in H<sup>2</sup> (Ω). Furthermore, the exponential convergence rates of the pressure and velocity are also proved by delicate energy methods.
基金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.
基金supported by the National Natural Science Foundation of China(52104049)the Young Elite Scientist Sponsorship Program by Beijing Association for Science and Technology(BYESS2023262)。
文摘Polymer flooding is an important means of improving oil recovery and is widely used in Daqing,Xinjiang,and Shengli oilfields,China.Different from conventional injection media such as water and gas,viscoelastic polymer solutions exhibit non-Newtonian and nonlinear flow behavior including shear thinning and shear thickening,polymer convection,diffusion,adsorption,retention,inaccessible pore volume,and reduced effective permeability.However,available well test model of polymer flooding wells generally simplifies these characteristics on pressure transient response,which may lead to inaccurate results.This work proposes a novel two-phase numerical well test model to better describe the polymer viscoelasticity and nonlinear flow behavior.Different influence factors that related to near-well blockage during polymer flooding process,including the degree of blockage(inner zone permeability),the extent of blockage(composite radius),and polymer flooding front radius are explored to investigate these impacts on bottom hole pressure responses.Results show that polymer viscoelasticity has a significant impact on the transitional flow segment of type curves,and the effects of near-well formation blockage and polymer concentration distribution on well test curves are very similar.Thus,to accurately interpret the degree of near-well blockage in injection wells,it is essential to first eliminate the influence of polymer viscoelasticity.Finally,a field case is comprehensively analyzed and discussed to illustrate the applicability of the proposed model.
基金the National Natural Science Foundation of China (Nos. 42302143, 42172159)China Geological Survey Project (No. DD20211350)support from the G. Albert Shoemaker endowment
文摘Clayey-silt natural gas hydrate reservoirs in the South China Sea exhibit loose and unconsolidated structures, heterogeneous pore structures, high clay mineral contents, and strong hydrophilicity. These characteristics complicate the gas-water two-phase flow process in porous media following hydrate decomposition, posing challenges for efficient development. This study examines the transport response of clayey-silt reservoir samples from the Shenhu area using gas-water two-phase flow experiments and CT scanning to explore changes in pore structure, gas-water distribution, and relative permeability under varying flow conditions. The results indicate that pore heterogeneity significantly influences flow characteristics. Gas preferentially displaces water in larger pores, forming fracture-like pores, which serve as preferential flow channels for gas migration. The preferential flow channels enhance gas-phase permeability up to 19 times that of the water phase when fluid pressures exceed total stresses. However,small pores retain liquid, leading to a high residual water saturation of 0.561. CT imaging reveals that these hydro-fractures improve gas permeability but also confine gas flow to specific channels. Pore network analysis shows that gas injection expands the pore-throat network, enhancing connectivity and forming fracture-like pores. Residual water remains trapped in smaller pores and throats, while structural changes, including new fractures, improve gas flow pathways and overall connectivity. Relative permeability curves demonstrate a narrow gas-water cocurrent-flow zone, a right-shifted iso-permeability point and high reservoir capillary pressure, indicating a strong "water-blocking" effect. The findings suggest that optimizing reservoir stimulation techniques to enhance fracture formation, reduce residual water saturation, and improve gas flow capacity is critical for efficient hydrate reservoir development.
文摘Shale gas production involves complex gas-water two-phase flow,with flow patterns in proppant-filled fractures playing a critical role in determining production efficiency.In this study,3D geometric models of 40/70 mesh ceramic particles and quartz sand proppant clusters were elaborated using computed tomography(CT)scanning.These models were used to develop a numerical simulation framework based on the lattice Boltzmann method(LBM),enabling the investigation of gas-water flow behavior within proppant-filled fractures under varying driving forces and surface tensions.Simulation results at a closure pressure of 15 MPa have revealed that ceramic particles exhibit a simpler and more porous internal structure than quartz sand of the same size.Under identical flow conditions,ceramic proppants demonstrate higher fluid replacement efficiency.Replacement efficiency increases with higher porosity,greater driving force,and lower surface tension.Furthermore,fluid displacement is strongly influenced by pore geometry:flow is faster in straighter and wider channels,with preferential movement through larger pores forming dominant flow paths.The replacement velocity exhibits a characteristic time evolution,initially rapid,then gradually decreasing,correlating positively with the development of these dominant channels.
基金the financial support from the National Natural Science Foundation of China (No.42102127)the Postdoctoral Research Foundation of China (No.2024 M751860)。
文摘Cleat serves as the primary flow pathway for coalbed methane(CBM)and water.However,few studies consider the impact of local contact on two-phase flow within cleats.A visual generalized model of endogenous cleats was constructed based on microfluidics.A microscopic and mesoscopic observation technique was proposed to simultaneously capture gas-liquid interface morphology of pores and throat and the two-phase flow characteristics in entire cleat system.The local contact characteristics of cleats reduced absolute permeability,which resulted in a sharp increase in the starting pressure.The reduced gas flow capacity narrowed the co-infiltration area and decreased water saturation at the isotonic point in a hydrophilic environment.The increased local contact area of cleats weakened gas phase flow capacity and narrowed the co-infiltration area.Jumping events occurred in methane-water flow due to altered porosity caused by local contact in cleats.The distribution of residual phases changed the jumping direction on the micro-scale as well as the dominant channel on the mesoscale.Besides,jumping events caused additional energy dissipation,which was ignored in traditional two-phase flow models.This might contribute to the overestimation of relative permeability.The work provides new methods and insights for investigating unsaturated flow in complex porous media.
文摘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 Key Research and Development Program of Yunnan Province(Grant No.202203AA080009)the Open Fund of State Key Laboratory of Hydraulics and Mountain River Engineering,Sichuan University(Grant No.SKHL2208).
文摘A two-phase flow model accelerated by graphical processing unit(GPU)is developed to solve fluid-solid interaction(FSI)using the sharp-interface immersed boundary method(IBM).This model solves the incompressible Navier-Stokes equations using the projection-based fractional step method in a fixed staggered Cartesian grid system.A volume of fluid(VOF)method with second-order accuracy is employed to trace the free surface.To represent the intricate surface geometry,the structure is discretized using the unstructured triangle mesh.Additionally,a ray tracing method is employed to classify fluid and solid points.A high-order stable scheme has been introduced to reconstruct the local velocity at interface points.Three FSI problems,including wave evolution around a breakwater,interaction between a periodic wave train and a moving float,and a 3-D moving object interacting with the free surface,were investigated to validate the accuracy and stability of the proposed model.The numerical results are in good agreement with the experimental data.Additionally,we evaluated the computational performance of the proposed GPU-based model.The GPU-based model achieved a 42.29 times speedup compared with the single-core CPU-based model in the three-dimension test.Additionally,the results regarding the time cost of each code section indicate that achieving more significant acceleration is associated with solving the turbulence,advection,and diffusion terms,while solving the pressure Poisson equation(PPE)saves the most time.Furthermore,the impact of grid number on computational efficiency indicates that as Fluid-solid interaction(FSI)immersed boundary method(IBM)graphical processing unit(GPU)two-phase flow moving rigid bodythe number of grids increases,the GPU-based model outperforms the multi-core CPU-based model.
基金supported by the NSFC Grant no.12271492the Natural Science Foundation of Henan Province of China Grant no.222300420550+1 种基金supported by the NSFC Grant no.12271498the National Key R&D Program of China Grant no.2022YFA1005202/2022YFA1005200.
文摘Due to the coupling between the hydrodynamic equation and the phase-field equation in two-phase incompressible flows,it is desirable to develop efficient and high-order accurate numerical schemes that can decouple these two equations.One popular and efficient strategy is to add an explicit stabilizing term to the convective velocity in the phase-field equation to decouple them.The resulting schemes are only first-order accurate in time,and it seems extremely difficult to generalize the idea of stabilization to the second-order or higher version.In this paper,we employ the spectral deferred correction method to improve the temporal accuracy,based on the first-order decoupled and energy-stable scheme constructed by the stabilization idea.The novelty lies in how the decoupling and linear implicit properties are maintained to improve the efficiency.Within the framework of the spatially discretized local discontinuous Galerkin method,the resulting numerical schemes are fully decoupled,efficient,and high-order accurate in both time and space.Numerical experiments are performed to validate the high-order accuracy and efficiency of the methods for solving phase-field models of two-phase incompressible flows.
文摘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.
