By combining with an improved model on engraving process,a two-phase flow interior ballistic model has been proposed to accurately predict the flow and energy conversion behaviors of pyrotechnic actuators.Using comput...By combining with an improved model on engraving process,a two-phase flow interior ballistic model has been proposed to accurately predict the flow and energy conversion behaviors of pyrotechnic actuators.Using computational fluid dynamics(CFD),the two-phase flow and piston engraving characteristics of a pyrotechnic actuator are investigated.Initially,the current model was utilized to examine the intricate,multi-dimensional flow,and energy conversion characteristics of the propellant grains and combustion gas within the pyrotechnic actuator chamber.It was discovered that the combustion gas on the wall's constant transition from potential to kinetic energy,along with the combined effect of the propellant motion,are what create the pressure oscillation within the chamber.Additionally,a numerical analysis was conducted to determine the impact of various parameters on the pressure oscillation and piston motion,including pyrotechnic charge,pyrotechnic particle size,and chamber structural dimension.The findings show that decreasing the pyrotechnic charge will lower the terminal velocity,while increasing and decreasing the pyrotechnic particle size will reduce the pressure oscillation in the chamber.The pyrotechnic particle size has minimal bearing on the terminal velocity.The results of this investigation offer a trustworthy forecasting instrument for comprehending and creating pyrotechnic actuator designs.展开更多
The 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.展开更多
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 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.展开更多
The challenge of wide brine source and its additional problems come from the economy(energy consumption and other costs),security(re-dissolution of surrounding salt rocks),and environment(groundwater pollution by brin...The challenge of wide brine source and its additional problems come from the economy(energy consumption and other costs),security(re-dissolution of surrounding salt rocks),and environment(groundwater pollution by brine)of salt cavern oil storage are worth examining to improve the efficiency of oil storage.Against this background,this work presented an operating mode of salt cavern oil and gas co-storage and using natural gas displacement for petroleum recovery.A gas-oil two-phase flow model with gas dissolution and exsolution was proposed to evaluate the application prospects of the new method precisely.Numerical studies indicated that the gas void fraction at the wellhead under quasi-steady state conditions is approximately 0.153,which belongs to bubbly flow,and the pressure at the wellhead of the central tube increased from 5.54 to 6.12 MPa during the entire transient flow stage,with an increase of 10.47%.Compared to the traditional method of using brine as the working fluid,the pump pressure rises from 2.92 to 14.01 MPa.However,if the new mode can be linked with the salt cavern gas storage and when the initial wellhead gas pressure exceeds 13 MPa,the energy consumption of the new method will be lower than that of the traditional brine-based operational mode.A new empirical formula is proposed to determine the two-phase flow pattern under different operating parameters.A special focus was given to energy consumption for oil recovery,which grows roughly in accordance with the operating pressure and oil recovery rate.However,the energy cost per volume of crude oil remains almost unchanged.This work provided a new solution for the serious brine problem and is expected to achieve petroleum recovery through natural gas displacement.展开更多
Gas-liquid two-phase flow in fractal porous media is pivotal for engineering applications,yet it remains challenging to be accurately characterized due to complex microstructure-flow interactions.This study establishe...Gas-liquid two-phase flow in fractal porous media is pivotal for engineering applications,yet it remains challenging to be accurately characterized due to complex microstructure-flow interactions.This study establishes a pore-scale numerical framework integratingMonte Carlo-generated fractal porousmedia with Volume of Fluid(VOF)simulations to unravel the coupling among pore distribution characterized by fractal dimension(Df),flow dynamics,and displacement efficiency.A pore-scale model based on the computed tomography(CT)microstructure of Berea sandstone is established,and the simulation results are compared with experimental data.Good agreement is found in phase distribution,breakthrough behavior,and flow path morphology,confirming the reliability of the numerical simulation method.Ten fractal porous media models with Df ranging from 1.25~1.7 were constructed using a Monte-Carlo approach.The gas-liquid two-phase flow dynamics was characterized using the VOF solver across gas injection rates of 0.05-5m/s,inwhich the time-resolved two-phase distribution patternswere systematically recorded.The results reveal that smaller fractal dimensions(Df=1.25~1.45)accelerate fingering breakthrough(peak velocity is 1.73 m/s at Df=1.45)due to a bimodal pore size distribution dominated by narrow channels.Increasing Df amplifies vorticity generation by about 3 times(eddy viscosity is 0.033 Pa⋅s at Df=1.7)through reduced interfacial curvature,while tortuosity-driven pressure differentials transition from sharp increases(0.4~6.3 Pa at Df=1.25~1.3)to inertial plateaus(4.8 Pa at Df=1.7).A nonlinear increase in equilibrium gas volume fraction(fav=0.692 at Df=1.7)emerges from residual gas saturation and turbulence-enhanced dispersion.This behavior is further modulated by flow velocity,with fav peaking at 0.72 under capillary-dominated conditions(0.05 m/s),but decreasing to 0.65 in the inertial regime(0.5 m/s).The work quantitatively links fractal topology to multiphase flow regimes,demonstrating the critical role of Df in governing preferential pathways,energy dissipation,and phase distribution.展开更多
In multiphase pumps transporting gas-liquid two-phase flows,the high-speed rotation of the impeller induces complex deformations in bubble shapes within the flow domain,making the prediction of gasliquid two-phase dra...In multiphase pumps transporting gas-liquid two-phase flows,the high-speed rotation of the impeller induces complex deformations in bubble shapes within the flow domain,making the prediction of gasliquid two-phase drag forces highly challenging in numerical simulations.To achieve precise prediction of the drag forces on irregular bubbles within multiphase pumps,this study modifies the existing bubble drag force model and applies the revised model to the prediction of gas-liquid two-phase flow within multiphase pumps.The research findings indicate that the modified drag force model significantly enhances the accuracy of predicting flow characteristics within the pump,particularly under high gas volume fraction conditions.The simulation results for gas phase distribution and vorticity exhibit strong agreement with experimental data.The modified drag model better captures the accumulation of the gas phase at the suction side of the impeller outlet.It also accurately predicts the vortex characteristics induced by bubble backflow from the trailing edges of the diffuser.Additionally,the adjustment of the drag coefficient enhances the model’s ability to represent local flow field characteristics,thereby optimizing the performance simulation methods of multiphase pumps.Compared to traditional drag force models,the modified model reduces prediction errors in head and efficiency by 36.4%and 27.5%,respectively.These results provide important theoretical foundations and model support for improving the accuracy of gas-liquid two-phase flow simulations and optimizing the design of multiphase pumps under high gas volume fraction conditions.展开更多
The gas-droplet two-phase reacting flow in a model combustor with the V-gutter flame holder is studied by an Eulerian-Lagrangian large-eddy simulation (LES) approach. The k-equation subgrid-scale model is used to simu...The gas-droplet two-phase reacting flow in a model combustor with the V-gutter flame holder is studied by an Eulerian-Lagrangian large-eddy simulation (LES) approach. The k-equation subgrid-scale model is used to simulate the subgrid eddy viscosity, and the eddy-break-up (EBU) combustion subgrid-scale model is used to determine the chemical reaction rate. A two-step turbulent combustion subgrid-scale model is employed for calculating carbon monoxide CO concentration, and the NO subgrid-scale pollutant formation model for the evaluation of the rate of NO formation. The heat flux model is applied to the prediction of radiant heat transfer. The gas phase is solved with the SIMPLE algorithm and a hybrid scheme in the staggered grid system. The liquid phase equations are solved in a Lagrangian frame in reference of the particle-source-in-cell (PSIC) algorithm. From simulation results, the exchange of mass, moment and energy between gas and particle fields for the reacting flow in the afterburner with a V-gutter flame holder can be obtained. By the comparison of experimental and simulation results, profile temperature and pollutant of the outlet are quite in agreement with experimental data. Results show that the LES approach for predicting the two-phase instantaneous reacting flow and pollutant emissions in the afterburner is feasible.展开更多
The influence of the squeeze film between the tube and the support structure on flow-induced vibrations is a critical factor in tube bundles subjected to two-phase cross-flow.This aspect can significantly alter the th...The influence of the squeeze film between the tube and the support structure on flow-induced vibrations is a critical factor in tube bundles subjected to two-phase cross-flow.This aspect can significantly alter the threshold for fluidelastic instability and affect heat transfer efficiency.This paper presents a mathematical model incorporating the squeeze film force between the tube and the support structure.We aim to clarify the mechanisms underlying fluidelastic instability in tube bundle systems exposed to two-phase flow.Using a self-developed computer program,we performed numerical calculations to examine the influence of the squeeze film on the threshold of fluidelastic instability in the tube bundle system.Furthermore,we analyzed how the thickness and length of the squeeze film affect both the underlying mechanisms and the critical velocity of fluidelastic instability.展开更多
This work investigated the dynamic behavior of vertical pipes conveying gas-liquid two-phase flow when subjected to external excitations at both ends.Even with minimal excitation amplitude,resonance can occur when the...This work investigated the dynamic behavior of vertical pipes conveying gas-liquid two-phase flow when subjected to external excitations at both ends.Even with minimal excitation amplitude,resonance can occur when the excitation frequency aligns with the natural frequency of the pipe,significantly increasing the degree of operational risk.The governing equation of motion based on the Euler-Bernoulli beam is derived for the relative deflection with stationary simply supported ends,with the effects of the external excitations represented by source terms distributed along the pipe length.The fourth-order partial differential equation is solved via the generalized integral transform technique(GITT),with the solution successfully verified via comparison with results in the literature.A comprehensive analysis of the vibration phenomena and changes in the motion state of the pipe is conducted for three classes of external excitation conditions:same frequency and amplitude(SFSA),same frequency but different amplitudes(SFDA),and different frequencies and amplitudes(DFDA).The numerical results show that with increasing gas volume fraction,the position corresponding to the maximum vibration displacement shifts upward.Compared with conditions without external excitation,the vibration displacement of the pipe conveying two-phase flow under external excitation increases significantly.The frequency of external excitation has a significant effect on the dynamic behavior of a pipe conveying two-phase flow.展开更多
With the increasing miniaturization of systems and surging demand for power density,accurate prediction and control of two-phase flow pressure drop have become a core challenge restricting the performance of microchan...With the increasing miniaturization of systems and surging demand for power density,accurate prediction and control of two-phase flow pressure drop have become a core challenge restricting the performance of microchannel heat exchangers.Pressure drop,a critical hydraulic characteristic,serves as both a natural constraint for cooling systems and determines the power required to pump the working fluid through microchannels.This paper reviews the characteristics,prediction models,and optimization measures of two-phase flow pressure drop for low-boiling-point working fluids in microchannels.It systematically analyzes key influencing factors such as fluid physical properties,operating conditions,channel geometry,and flow patterns,and discusses the complex mechanisms of pressure drop under the coupling effect of multi-physical fields.Mainstream prediction models are reviewed:the homogeneous flow model simplifies calculations but shows large deviations at low quality;the separated flow model considers interphase interactions and can be applied to micro-scales after modification;the flow-pattern-based model performs zoned modeling but relies on subjective classification;machine learning improves prediction accuracy but faces the“black-box”problem.In terms of optimization,channel designs are improved through porous structures and micro-rib arrays,and flow rate distribution is optimized using splitters to balance pressure drop and heat transfer performance.This study provides theoretical support for microchannel thermal management in high-power-density devices.展开更多
Polymer flooding is an important means of improving oil recovery and is widely used in Daqing,Xinjiang,and Shengli oilfields,China.Different from conventional injection media such as water and gas,viscoelastic polymer...Polymer flooding is an important means of improving oil recovery and is widely used in Daqing,Xinjiang,and Shengli oilfields,China.Different from conventional injection media such as water and gas,viscoelastic polymer solutions exhibit non-Newtonian and nonlinear flow behavior including shear thinning and shear thickening,polymer convection,diffusion,adsorption,retention,inaccessible pore volume,and reduced effective permeability.However,available well test model of polymer flooding wells generally simplifies these characteristics on pressure transient response,which may lead to inaccurate results.This work proposes a novel two-phase numerical well test model to better describe the polymer viscoelasticity and nonlinear flow behavior.Different influence factors that related to near-well blockage during polymer flooding process,including the degree of blockage(inner zone permeability),the extent of blockage(composite radius),and polymer flooding front radius are explored to investigate these impacts on bottom hole pressure responses.Results show that polymer viscoelasticity has a significant impact on the transitional flow segment of type curves,and the effects of near-well formation blockage and polymer concentration distribution on well test curves are very similar.Thus,to accurately interpret the degree of near-well blockage in injection wells,it is essential to first eliminate the influence of polymer viscoelasticity.Finally,a field case is comprehensively analyzed and discussed to illustrate the applicability of the proposed model.展开更多
A second-order moment two-phase turbulence model for simulating dense gas-particle flows （USM-Θ model）, combining the unified second-order moment twophase turbulence model for dilute gas-particle flows with the kin...A second-order moment two-phase turbulence model for simulating dense gas-particle flows （USM-Θ model）, combining the unified second-order moment twophase turbulence model for dilute gas-particle flows with the kinetic theory of particle collision, is proposed. The interaction between gas and particle turbulence is simulated using the transport equation of two-phase velocity correlation with a two-time-scale dissipation closure. The proposed model is applied to simulate dense gas-particle flows in a horizontal channel and a downer. Simulation results and their comparison with experimental results show that the model accounting for both anisotropic particle turbulence and particle-particle collision is obviously better than models accounting for only particle turbulence or only particle-particle collision. The USM-Θ model is also better than the k-ε-kp-Θ model and the k-ε-kp-εp-Θ model in that the first model can simulate the redistribution of anisotropic particle Reynolds stress components due to inter-particle collision, whereas the second and third models cannot.展开更多
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.展开更多
基金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(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.
基金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.
基金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.
基金The financial support from the National Science and Technology Major Project,China(No.2024ZD1004107)the Natural Science Foundation of Wuhan(No.2024040701010062)is greatly appreciated。
文摘The challenge of wide brine source and its additional problems come from the economy(energy consumption and other costs),security(re-dissolution of surrounding salt rocks),and environment(groundwater pollution by brine)of salt cavern oil storage are worth examining to improve the efficiency of oil storage.Against this background,this work presented an operating mode of salt cavern oil and gas co-storage and using natural gas displacement for petroleum recovery.A gas-oil two-phase flow model with gas dissolution and exsolution was proposed to evaluate the application prospects of the new method precisely.Numerical studies indicated that the gas void fraction at the wellhead under quasi-steady state conditions is approximately 0.153,which belongs to bubbly flow,and the pressure at the wellhead of the central tube increased from 5.54 to 6.12 MPa during the entire transient flow stage,with an increase of 10.47%.Compared to the traditional method of using brine as the working fluid,the pump pressure rises from 2.92 to 14.01 MPa.However,if the new mode can be linked with the salt cavern gas storage and when the initial wellhead gas pressure exceeds 13 MPa,the energy consumption of the new method will be lower than that of the traditional brine-based operational mode.A new empirical formula is proposed to determine the two-phase flow pattern under different operating parameters.A special focus was given to energy consumption for oil recovery,which grows roughly in accordance with the operating pressure and oil recovery rate.However,the energy cost per volume of crude oil remains almost unchanged.This work provided a new solution for the serious brine problem and is expected to achieve petroleum recovery through natural gas displacement.
基金funded by the National Key R&D Program of China,China(Grant No.2023YFB4005500)National Natural Science Foundation of China,China(Grant Nos.52379113 and 52379114).
文摘Gas-liquid two-phase flow in fractal porous media is pivotal for engineering applications,yet it remains challenging to be accurately characterized due to complex microstructure-flow interactions.This study establishes a pore-scale numerical framework integratingMonte Carlo-generated fractal porousmedia with Volume of Fluid(VOF)simulations to unravel the coupling among pore distribution characterized by fractal dimension(Df),flow dynamics,and displacement efficiency.A pore-scale model based on the computed tomography(CT)microstructure of Berea sandstone is established,and the simulation results are compared with experimental data.Good agreement is found in phase distribution,breakthrough behavior,and flow path morphology,confirming the reliability of the numerical simulation method.Ten fractal porous media models with Df ranging from 1.25~1.7 were constructed using a Monte-Carlo approach.The gas-liquid two-phase flow dynamics was characterized using the VOF solver across gas injection rates of 0.05-5m/s,inwhich the time-resolved two-phase distribution patternswere systematically recorded.The results reveal that smaller fractal dimensions(Df=1.25~1.45)accelerate fingering breakthrough(peak velocity is 1.73 m/s at Df=1.45)due to a bimodal pore size distribution dominated by narrow channels.Increasing Df amplifies vorticity generation by about 3 times(eddy viscosity is 0.033 Pa⋅s at Df=1.7)through reduced interfacial curvature,while tortuosity-driven pressure differentials transition from sharp increases(0.4~6.3 Pa at Df=1.25~1.3)to inertial plateaus(4.8 Pa at Df=1.7).A nonlinear increase in equilibrium gas volume fraction(fav=0.692 at Df=1.7)emerges from residual gas saturation and turbulence-enhanced dispersion.This behavior is further modulated by flow velocity,with fav peaking at 0.72 under capillary-dominated conditions(0.05 m/s),but decreasing to 0.65 in the inertial regime(0.5 m/s).The work quantitatively links fractal topology to multiphase flow regimes,demonstrating the critical role of Df in governing preferential pathways,energy dissipation,and phase distribution.
基金funded by Sichuan Natural Science Foundation Outstanding Youth Science Foundation(No.2024NSFJQ0012)Key project of Regional Innovation and Development Joint Fund of National Natural Science Foundation(No.U23A20669)Sichuan Science and Technology Program(2022ZDZX0041).
文摘In multiphase pumps transporting gas-liquid two-phase flows,the high-speed rotation of the impeller induces complex deformations in bubble shapes within the flow domain,making the prediction of gasliquid two-phase drag forces highly challenging in numerical simulations.To achieve precise prediction of the drag forces on irregular bubbles within multiphase pumps,this study modifies the existing bubble drag force model and applies the revised model to the prediction of gas-liquid two-phase flow within multiphase pumps.The research findings indicate that the modified drag force model significantly enhances the accuracy of predicting flow characteristics within the pump,particularly under high gas volume fraction conditions.The simulation results for gas phase distribution and vorticity exhibit strong agreement with experimental data.The modified drag model better captures the accumulation of the gas phase at the suction side of the impeller outlet.It also accurately predicts the vortex characteristics induced by bubble backflow from the trailing edges of the diffuser.Additionally,the adjustment of the drag coefficient enhances the model’s ability to represent local flow field characteristics,thereby optimizing the performance simulation methods of multiphase pumps.Compared to traditional drag force models,the modified model reduces prediction errors in head and efficiency by 36.4%and 27.5%,respectively.These results provide important theoretical foundations and model support for improving the accuracy of gas-liquid two-phase flow simulations and optimizing the design of multiphase pumps under high gas volume fraction conditions.
文摘The gas-droplet two-phase reacting flow in a model combustor with the V-gutter flame holder is studied by an Eulerian-Lagrangian large-eddy simulation (LES) approach. The k-equation subgrid-scale model is used to simulate the subgrid eddy viscosity, and the eddy-break-up (EBU) combustion subgrid-scale model is used to determine the chemical reaction rate. A two-step turbulent combustion subgrid-scale model is employed for calculating carbon monoxide CO concentration, and the NO subgrid-scale pollutant formation model for the evaluation of the rate of NO formation. The heat flux model is applied to the prediction of radiant heat transfer. The gas phase is solved with the SIMPLE algorithm and a hybrid scheme in the staggered grid system. The liquid phase equations are solved in a Lagrangian frame in reference of the particle-source-in-cell (PSIC) algorithm. From simulation results, the exchange of mass, moment and energy between gas and particle fields for the reacting flow in the afterburner with a V-gutter flame holder can be obtained. By the comparison of experimental and simulation results, profile temperature and pollutant of the outlet are quite in agreement with experimental data. Results show that the LES approach for predicting the two-phase instantaneous reacting flow and pollutant emissions in the afterburner is feasible.
基金financially supported by the National Natural Science Foundation of China(Grant No.12072336).
文摘The influence of the squeeze film between the tube and the support structure on flow-induced vibrations is a critical factor in tube bundles subjected to two-phase cross-flow.This aspect can significantly alter the threshold for fluidelastic instability and affect heat transfer efficiency.This paper presents a mathematical model incorporating the squeeze film force between the tube and the support structure.We aim to clarify the mechanisms underlying fluidelastic instability in tube bundle systems exposed to two-phase flow.Using a self-developed computer program,we performed numerical calculations to examine the influence of the squeeze film on the threshold of fluidelastic instability in the tube bundle system.Furthermore,we analyzed how the thickness and length of the squeeze film affect both the underlying mechanisms and the critical velocity of fluidelastic instability.
基金financially supported by the Key Research and Development Program of Shandong Province(Grant Nos.2022CXGC020405,2023CXGC010415 and 2025TSGCCZZB0238)the National Natural Science Foundation of China(Grant No.52171288)the financial support from CNPq,FAPERJ,ANP,Embrapii,and China National Petroleum Corporation(CNPC).
文摘This work investigated the dynamic behavior of vertical pipes conveying gas-liquid two-phase flow when subjected to external excitations at both ends.Even with minimal excitation amplitude,resonance can occur when the excitation frequency aligns with the natural frequency of the pipe,significantly increasing the degree of operational risk.The governing equation of motion based on the Euler-Bernoulli beam is derived for the relative deflection with stationary simply supported ends,with the effects of the external excitations represented by source terms distributed along the pipe length.The fourth-order partial differential equation is solved via the generalized integral transform technique(GITT),with the solution successfully verified via comparison with results in the literature.A comprehensive analysis of the vibration phenomena and changes in the motion state of the pipe is conducted for three classes of external excitation conditions:same frequency and amplitude(SFSA),same frequency but different amplitudes(SFDA),and different frequencies and amplitudes(DFDA).The numerical results show that with increasing gas volume fraction,the position corresponding to the maximum vibration displacement shifts upward.Compared with conditions without external excitation,the vibration displacement of the pipe conveying two-phase flow under external excitation increases significantly.The frequency of external excitation has a significant effect on the dynamic behavior of a pipe conveying two-phase flow.
基金supported by the Beijing Municipal Science&Technology Commission(Z231100006123010).
文摘With the increasing miniaturization of systems and surging demand for power density,accurate prediction and control of two-phase flow pressure drop have become a core challenge restricting the performance of microchannel heat exchangers.Pressure drop,a critical hydraulic characteristic,serves as both a natural constraint for cooling systems and determines the power required to pump the working fluid through microchannels.This paper reviews the characteristics,prediction models,and optimization measures of two-phase flow pressure drop for low-boiling-point working fluids in microchannels.It systematically analyzes key influencing factors such as fluid physical properties,operating conditions,channel geometry,and flow patterns,and discusses the complex mechanisms of pressure drop under the coupling effect of multi-physical fields.Mainstream prediction models are reviewed:the homogeneous flow model simplifies calculations but shows large deviations at low quality;the separated flow model considers interphase interactions and can be applied to micro-scales after modification;the flow-pattern-based model performs zoned modeling but relies on subjective classification;machine learning improves prediction accuracy but faces the“black-box”problem.In terms of optimization,channel designs are improved through porous structures and micro-rib arrays,and flow rate distribution is optimized using splitters to balance pressure drop and heat transfer performance.This study provides theoretical support for microchannel thermal management in high-power-density devices.
基金supported by the National Natural Science Foundation of China(52104049)the Young Elite Scientist Sponsorship Program by Beijing Association for Science and Technology(BYESS2023262)。
文摘Polymer flooding is an important means of improving oil recovery and is widely used in Daqing,Xinjiang,and Shengli oilfields,China.Different from conventional injection media such as water and gas,viscoelastic polymer solutions exhibit non-Newtonian and nonlinear flow behavior including shear thinning and shear thickening,polymer convection,diffusion,adsorption,retention,inaccessible pore volume,and reduced effective permeability.However,available well test model of polymer flooding wells generally simplifies these characteristics on pressure transient response,which may lead to inaccurate results.This work proposes a novel two-phase numerical well test model to better describe the polymer viscoelasticity and nonlinear flow behavior.Different influence factors that related to near-well blockage during polymer flooding process,including the degree of blockage(inner zone permeability),the extent of blockage(composite radius),and polymer flooding front radius are explored to investigate these impacts on bottom hole pressure responses.Results show that polymer viscoelasticity has a significant impact on the transitional flow segment of type curves,and the effects of near-well formation blockage and polymer concentration distribution on well test curves are very similar.Thus,to accurately interpret the degree of near-well blockage in injection wells,it is essential to first eliminate the influence of polymer viscoelasticity.Finally,a field case is comprehensively analyzed and discussed to illustrate the applicability of the proposed model.
基金the Special Funds for Major State Basic Research of China(G-1999-0222-08)the National Natural Science Foundation of China(50376004)Ph.D.Program Foundation,Ministry of Education of China(20030007028)
文摘A second-order moment two-phase turbulence model for simulating dense gas-particle flows （USM-Θ model）, combining the unified second-order moment twophase turbulence model for dilute gas-particle flows with the kinetic theory of particle collision, is proposed. The interaction between gas and particle turbulence is simulated using the transport equation of two-phase velocity correlation with a two-time-scale dissipation closure. The proposed model is applied to simulate dense gas-particle flows in a horizontal channel and a downer. Simulation results and their comparison with experimental results show that the model accounting for both anisotropic particle turbulence and particle-particle collision is obviously better than models accounting for only particle turbulence or only particle-particle collision. The USM-Θ model is also better than the k-ε-kp-Θ model and the k-ε-kp-εp-Θ model in that the first model can simulate the redistribution of anisotropic particle Reynolds stress components due to inter-particle collision, whereas the second and third models cannot.
基金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.
