The low diffusion （LD） particle method, proposed by Burt and Boyd, is modified for the near-continuum two-phase flow simulations. The LD method has the advantages of easily coupling with the direct simulation Monte ...The low diffusion （LD） particle method, proposed by Burt and Boyd, is modified for the near-continuum two-phase flow simulations. The LD method has the advantages of easily coupling with the direct simulation Monte Carlo （DSMC） method for multi-scale flow simulations and dramatically reducing the numerical diffusion error and statistical scatter of the equilibrium particle methods. Liquidor solid-phase particles are introduced in the LD method. Their velocity and temperature updating are respectively, calculated from the motion equation and the temperature equation according to the local gas properties. Coupling effects from condensed phase to gas phase are modeled as momentum and energy sources, which are respectively, equal to the negative values of the total momentum and energy increase in liquid or solid phase. The modified method is compared with theoretical results for unsteady flows, and good agreements are obtained to indicate the reliability of the one-way gas-to-particle coupling models. Hybrid LD-DSMC algorithm is implemented and performed for nozzle discharging gas-liquid flow to show the prospect of the LD-DSMC scheme for multi-scale two-phase flow simulations.展开更多
Air entrapped in liquid metal during the mold filling process seriously affects the casting quality,thus it is important to track its behavior in the mold cavity.A liquid-gas two-phase flow model is developed to descr...Air entrapped in liquid metal during the mold filling process seriously affects the casting quality,thus it is important to track its behavior in the mold cavity.A liquid-gas two-phase flow model is developed to describe the mold filling process and predict the air entrapment defect.The model is based on the combination of SOLA and Level Set Method.The pressure and velocity fields are calculated by SOLA,and the interface movement is simulated by Level Set method as the most common interface tracking method in recent years.In order to validate the feasibility of the model,the liquid-gas two-phase simulation results were tested by the broken dam problem and the S-shaped experiment.Comparison between the experiments and simulation results show that Level Set method might be a very promising tool in two-phase flow simulation during the mold filling process.展开更多
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.展开更多
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.展开更多
This study investigates the droplet formation for the liquid–liquid two-phase flow within a square T-junction microchannel through numerical simulation using volume of fluid method and experimental visualization usin...This study investigates the droplet formation for the liquid–liquid two-phase flow within a square T-junction microchannel through numerical simulation using volume of fluid method and experimental visualization using high-speed camera imaging.The T-junction microchannel has a cross-sectional width of 0.6 mm and a total length of 27.3 mm.The solution of cyclohexane with 2%and 3%mass concentrations of sorbitan trioleate surfactant were used as the continuous phase,and water was used as the dispersed phase.Slug flow,characteristic of squeezing regime,were predominantly observed.The effects of liquid–liquid two-phase flow rate ratio,and dimensionless number on droplet size,and pressure drop were investigated.The squeezing regime was mapped for 0.0005≤Ca_(c)≤0.0052(capillary number)and 0.1≤q≤10(flow rate ratio).The pressure drops of slugs were in the range from 40 Pa to 200 Pa.The slug lengths were measured between 1 mm and 9 mm.A universal flow map dependent on Ca_(c)Re_(d)^(0.5) are projected to investigate the droplet formation behavior in T-junction microchannel.Correlation expressions are proposed to predict pressure drops and the slug lengths for liquid–liquid two-phase flow in a square T-junction microchannel,using dimensionless numbers such as flow rate ratio and capillary number.The result shows that large continuous phase flow rates facilitate smaller slugs,whereas higher dispersed phase flow rates result in longer shorts.展开更多
The Sichuan-Yunnan region,located at the southeastern margin of the Qinghai-Xizang Plateau,serves as a key channel for the southeastward extrusion of plateau material.The characteristics of crustal deformation and the...The Sichuan-Yunnan region,located at the southeastern margin of the Qinghai-Xizang Plateau,serves as a key channel for the southeastward extrusion of plateau material.The characteristics of crustal deformation and the mechanisms of deep material flow have been central topics of interest in geoscience research.In this work,a three-dimensional viscoelastic-plastic finite element model including the upper and mid-lower crust was established,constrained by GNSS horizontal crustal velocity observations and incorporating maj or active faults and geophysical survey data to explore the contribution of mid-lower crustal flow to surface deformation and its coupling with faults.Comparison of modeling experiments shows that relying solely on boundary loading or uniform layering assumptions fails to reproduce the GNSS observed velocities.We introduce a mid-lower crustal low-velocity weak zone,derived from the latest seismic velocity structure models.The new model improves the fit to GNSS observations.Tests of different viscosity coefficients in the low-velocity zone indicate an optimal viscosity range of 7.5×10^(19)-1×10^(20)Pa·s.Vertical profiles reveal that mid-lower crustal material motion is mainly concentrated at depths of 20-40 km,forming localized channelized flow in low-velocity zone with a typical Poiseuille velocity profile which indicates a ductile,fluid-like behavior with the lowvelocity zone serving as primary pathways for deep material transport.The results further show that under the geometric constraints of upper-crustal faults,the mid-lower crustal flow contributes approximately 1-3 mm/a to surface deformation,primarily concentrated along major faults.This indicates that faults play a key role in constraining and modulating the transmission of deep-seated dynamics to shallow surface deformation.However,the contribution of mid-lower crustal flow is also significant;neglecting its influence on surface deformation would lead to an incomplete understanding of the deformation pattern and bias the interpretation of block boundaries and crustal kinematic segmentation.展开更多
Although extensive research has been conducted on CO_(2)-enhanced coalbed methane(CO_(2)-ECBM)recovery,most prior studies have focused on the impact of gas adsorption-induced swelling on coal permeability under equili...Although extensive research has been conducted on CO_(2)-enhanced coalbed methane(CO_(2)-ECBM)recovery,most prior studies have focused on the impact of gas adsorption-induced swelling on coal permeability under equilibrium conditions.This paper presents a comprehensive thermo-hydro-mechanical-chemical(THMC)model that integrates thermal expansion and heat conduction(T),gas diffusion in the matrix and gas-water two-phase flow in the fractures(H),matrix and fracture deformation due to poroelasticity(M),and non-equilibrium binary gas adsorption-induced matrix swelling(C)during CO_(2)-ECBM recovery.The accuracy of the proposed model was verified through experimental data,and the model was simulated using finite element method(FEM)software.Simulation results indicate that the permeability evolution can be categorized into three stages.Ignoring the impact of water on gas adsorption properties would lead to an overestimation of the influence of adsorption-induced swelling,while disregarding non-equilibrium adsorption underestimates it.An examination of five designed cases identified critical factors influencing permeability.Parametric analysis shows that increases in the injection pressure,the injection temperature,and the initial permeability promote the competitive adsorption-induced swelling between CH_(4)and CO_(2),leading to increased CH_(4)production and CO_(2)storage.Conversely,an increase in initial water saturation reduces available gas flow space,decreasing both CH_(4)production and CO_(2)storage.Higher irreducible water saturation favors early gas recovery,while lower irreducible water saturation is more advantageous for long-term recovery.展开更多
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.展开更多
In order to speed underwater launch of minor-caliber weapons,a sealing device can be set in front of underwater muzzle to separate water,preventing the muzzle from water immersion.By establishing and simplifying the m...In order to speed underwater launch of minor-caliber weapons,a sealing device can be set in front of underwater muzzle to separate water,preventing the muzzle from water immersion.By establishing and simplifying the model of underwater weapon sealing device and unstructured mesh computing domain model based on computational fluid dynamics（CFD）,dynamic mesh and user defined function（UDF）,the N-S equation is solved and the numerical analysis and calculation of the complex two-phase flow inside the sealing device are carried out.The results show that the gas discharged from the sealing device is conducive to the formation of the projectile supercavity.When the projectile is launched at 5munder water,the shock wave before and after the projectile has impact on the box body up to 100 MPa,therefore the sealing device must be strong enough.The research results have the vital significance to the design of underwater weapon sealing device and the formation of the projectile supercavitation.展开更多
The two-phase flow in porous media is affected by multiple factors.In the present study,a two-dimensional numerical model of porous media was developed using the actual pore structure of the core sample.The phase fiel...The two-phase flow in porous media is affected by multiple factors.In the present study,a two-dimensional numerical model of porous media was developed using the actual pore structure of the core sample.The phase field method was utilized to simulate the impact of displacement velocity,the water-gas viscosity ratio,and the density ratio on the flow behavior of two-phase fluids in porous media.The effectiveness of displacement was evaluated by analyzing CO_(2)saturation levels.The results indicate that the saturation of CO_(2)in porous media increased as the displacement velocity increased.When the displacement velocity exceeded 0.01 m/s,there was a corresponding increase in CO_(2)saturation.Conversely,when the displacement velocity was below this threshold,the impact on CO_(2)saturation was minimal.An“inflection point,”M3,was present in the viscosity ratio.When the viscosity of CO_(2)is less than 8.937×10^(-5)Pa·s(viscosity ratio below M3),variations in the viscosity of CO_(2)had little impact on its saturation.Conversely,when the viscosity of CO_(2)exceeded 8.937×10^(-5)Pa·s(viscosity ratio greater than M3),saturation increased with an increase in the viscosity ratio.In terms of the density ratio,the saturation of CO_(2)increased monotonically with an increase in the density ratio.Similarly,increasing density ratios resulted in a monotonic increase in CO_(2)saturation,though this trend was less pronounced in numerical simulations.Analysis results of displacement within dead-end pores using pressure and velocity diagrams reveal eddy currents as contributing factors.Finally,the impact of pore throat structure on the formation of dominant channels was examined.展开更多
Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of...Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of unsteady flow characteristics on solid-liquid two-phase flow and pump performance were researched under design condition. The results show that clocking effect has a significant influence on the flow in pump, and the fluctuation of flow velocity and pressure is obvious, particularly near the volute tongue, at the position of small sections of volute and within diffuser. Clocking effect has a more influence on liquid-phase than on solid-phase, and the wake-jet structure of relative velocity of solid-phase is less obvious than liquid-phase near the volute tongue and the impeller passage outlet. The fluctuation of relative velocity of solid-phase flow is 7.6% smaller than liquid-phase flow at the impeller outlet on circular path. Head and radial forces of the impeller are 8.1% and 85.7% of fluctuation, respectively. The results provide a theoretical basis for further research for turbulence, improving efficient, reducing the hydraulic losses and wear. Finally, field tests were carried out to verify the operation and wear of slurry pump.展开更多
A mathematical model is set to evaluate the 3-D dense solid-liquid two-phaseturbulent flow in a non-clogging mud pump, the flow feature in the impeller channel is simulatedwith the tool of IPSA. Meanwhile, resort to T...A mathematical model is set to evaluate the 3-D dense solid-liquid two-phaseturbulent flow in a non-clogging mud pump, the flow feature in the impeller channel is simulatedwith the tool of IPSA. Meanwhile, resort to TECPLOT as the post-processor, the simulation results isvisualized. The results show the main flow characteristics: There exists backflow and aberrantvelocities at inlet area and a relative velocity slip between two phases; A jet-wake flow pattern isdiscerned around the shroud-suction side area; The relative velocity vector of solid phase iscloser to the pressure surface than that of liquid phase and the trend is more obvious with theincrease of diameter; The kinetic energy of turbulence k and the dissipation rate e reach theirpeaks at the corner of pressure and suction surface. The simulation results show a good agreementwith the experimental flow features in the impeller channel, which prove the turbulent model used isvalid and provide a theoretical design basis to non-clogging pumps.展开更多
An improved large eddy simulation using a dynamic second-order sub-grid-scale （SGS） stress model has been developed to model the governing equations of dense turbulent particle-liquid two-phase flows in a rotating c...An improved large eddy simulation using a dynamic second-order sub-grid-scale （SGS） stress model has been developed to model the governing equations of dense turbulent particle-liquid two-phase flows in a rotating coordinate system, and continuity is conserved by a mass-weighted method to solve the filtered governing equations. In the cur- rent second-order SGS model, the SGS stress is a function of both the resolved strain-rate and rotation-rate tensors, and the model parameters are obtained from the dimensional consistency and the invariants of the strain-rate and the rotation-rate tensors. In the numerical calculation, the finite volume method is used to discretize the governing equations with a staggered grid system. The SIMPLEC algorithm is applied for the solution of the discretized governing equations. Body- fitted coordinates are used to simulate the two-phase flows in complex geometries. Finally the second-order dynamic SGS model is successfully applied to simulate the dense turbu-lent particle-liquid two-phase flows in a centrifugal impeller. The predicted pressure and velocity distributions are in good agreement with experimental results.展开更多
The flow with solid-liquid two-phase media inside centrifugal pumps is very complicated and the relevant method for the hydraulic design is still immature so far. There exist two main problems in the operation of the ...The flow with solid-liquid two-phase media inside centrifugal pumps is very complicated and the relevant method for the hydraulic design is still immature so far. There exist two main problems in the operation of the two-phase flow pumps, i.e., low overall efficiency and severe abrasion. In this study, the three-dimensional, steady, incompressible, and turbulent solid-liquid two-phase flows in a low-specific-speed centrifugal pump are numerically simulated and analyzed by using a computational fluid dynamics （CFD） code based on the mixture model of the two-phase flow and the RNG k-~ two-equation turbulence model, in which the influences of rotation and curvature are fully taken into account. The coupling between impeller and volute is implemented by means of the frozen rotor method. The simulation results predicted indicate that the solid phase properties in two-phase flow, especially the concentration, the particle diameter and the density, have strong effects on the hydraulic performance of the pump. Both the pump head and the efficiency are reduced with increasing particle diameter or concentration. However, the effect of particle density on the performance is relatively minor. An obvious jet-wake flow structure is presented near the volute tongue and becomes more remarkable with increasing solid phase concentration. The suction side of the blade is subject to much more severe abrasion than the pressure side. The obtained results preliminarily reveal the characteristics of solid-liquid two-phase flow in the centrifugal pump, and are helpful for improvement and empirical correction in the hydraulic design of centrifugal pumps.展开更多
Streamline simulation is developed to simulate waterflooding in fractured reservoirs. Conventional reservoir simulation methods for fluid flow simulation in large and complex reservoirs are very costly and time consum...Streamline simulation is developed to simulate waterflooding in fractured reservoirs. Conventional reservoir simulation methods for fluid flow simulation in large and complex reservoirs are very costly and time consuming. In streamline method, transport equations are solved on one-dimensional streamlines to reduce the computation time with less memory for simulation. First, pressure equation is solved on an Eulerian grid and streamlines are traced. Defining the "time of flight", saturation equations are mapped and solved on streamlines. Finally, the results are mapped back on Eulerian grid and the process is repeated until the simulation end time. The waterflooding process is considered in a fractured reservoir using the dual porosity model. Afterwards, a computational code is developed to solve the same problem by the IMPES method and the results of streamline simulation are compared to those of the IMPES and a commercial software. Finally, the accuracy and efficiency of streamline simulator for simulation of two-phase flow in fractured reservoirs has been proved.展开更多
The flow field of gas and liquid in a φ150mm rotating-stream-tray (RST) scrubber is simulated by using computational fluid dynamic (CFD) method. The sismulation is based on the two-equation RNG κ-ε turbulence model...The flow field of gas and liquid in a φ150mm rotating-stream-tray (RST) scrubber is simulated by using computational fluid dynamic (CFD) method. The sismulation is based on the two-equation RNG κ-ε turbulence model, Eulerian multiphase model, and a real-shape 3D model with a huge number of meshes. The simulation results include detailed information about velocity, pressure, volume fraction and so on. Some features of the flow field are obtained: liquid is atomized in a thin annular zone; a high velocity air zone prevents water drops at the bottom from flying towards the wall; the pressure varies sharply at the end of blades and so on. The results will be helpful for structure optimization and engineering design.展开更多
Chokes are one of the most important components of downhole flow-control equipment. The particle erosion mathematical model, which considers particle-particle interaction, was established and used to simulate solid pa...Chokes are one of the most important components of downhole flow-control equipment. The particle erosion mathematical model, which considers particle-particle interaction, was established and used to simulate solid particle movement as well as particle erosion characteristics of the solid-liquid two-phase flow in a choke. The corresponding erosion reduction approach by setting ribs on the inner wall of the choke was advanced. This mathematical model includes three parts： the flow field simulation of the continuous carrier fluid by an Eulerian approach, the particle interaction simulation using the discrete particle hard sphere model by a Lagrangian approach and calculation of erosion rate using semiempirical correlations. The results show that particles accumulated in a narrow region from inlet to outlet of the choke and the dominating factor affecting particle motion is the fluid drag force. As a result, the optimization of rib geometrical parameters indicates that good anti-erosion performance can be achieved by four ribs, each of them with a height （H） of 3 mm and a width （B） of 5 mm equaling the interval between ribs （L）.展开更多
Debris flows are typical two-phase flows, which commonly accompany erosion in mountainous areas, and may destroy bridge engineering by scouring. In this study, a physically-based two-phase model is applied for the sim...Debris flows are typical two-phase flows, which commonly accompany erosion in mountainous areas, and may destroy bridge engineering by scouring. In this study, a physically-based two-phase model is applied for the simulation of debris flow scouring of bridge pier. In this model, the shear stress of debris flow on an erodible bed is considered to be a function of the solid shear stress, fluid shear stress, and volume fraction; accordingly, the erosion is incorporated into the two-phase model. Using a highaccuracy computational scheme based on the finite volume method, the model is employed for simulating a dynamic debris flow over an erodible bed. The numerical results are consistent with the experimental data, and verify the feasibility of the two-phase model. Moreover, a simple numerical test is performed to exhibit the fundamental behaviour of debris flow scouring of bridge pier, which shows that the degree of erosion on each side of the pier is higher compared to other areas. The scouring depth is influenced by the variations of solid volume fraction and velocity of debris flow and pier width.展开更多
文摘The low diffusion （LD） particle method, proposed by Burt and Boyd, is modified for the near-continuum two-phase flow simulations. The LD method has the advantages of easily coupling with the direct simulation Monte Carlo （DSMC） method for multi-scale flow simulations and dramatically reducing the numerical diffusion error and statistical scatter of the equilibrium particle methods. Liquidor solid-phase particles are introduced in the LD method. Their velocity and temperature updating are respectively, calculated from the motion equation and the temperature equation according to the local gas properties. Coupling effects from condensed phase to gas phase are modeled as momentum and energy sources, which are respectively, equal to the negative values of the total momentum and energy increase in liquid or solid phase. The modified method is compared with theoretical results for unsteady flows, and good agreements are obtained to indicate the reliability of the one-way gas-to-particle coupling models. Hybrid LD-DSMC algorithm is implemented and performed for nozzle discharging gas-liquid flow to show the prospect of the LD-DSMC scheme for multi-scale two-phase flow simulations.
基金National High Technology Research and Development Program of China(863program)(2006AA04Z140)National Natural Science Foundation of China(NSFC)(50605024)
文摘Air entrapped in liquid metal during the mold filling process seriously affects the casting quality,thus it is important to track its behavior in the mold cavity.A liquid-gas two-phase flow model is developed to describe the mold filling process and predict the air entrapment defect.The model is based on the combination of SOLA and Level Set Method.The pressure and velocity fields are calculated by SOLA,and the interface movement is simulated by Level Set method as the most common interface tracking method in recent years.In order to validate the feasibility of the model,the liquid-gas two-phase simulation results were tested by the broken dam problem and the S-shaped experiment.Comparison between the experiments and simulation results show that Level Set method might be a very promising tool in two-phase flow simulation during the mold filling process.
文摘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.
基金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.
基金supports for this project from the National Natural Science Foundation of China(22378295).
文摘This study investigates the droplet formation for the liquid–liquid two-phase flow within a square T-junction microchannel through numerical simulation using volume of fluid method and experimental visualization using high-speed camera imaging.The T-junction microchannel has a cross-sectional width of 0.6 mm and a total length of 27.3 mm.The solution of cyclohexane with 2%and 3%mass concentrations of sorbitan trioleate surfactant were used as the continuous phase,and water was used as the dispersed phase.Slug flow,characteristic of squeezing regime,were predominantly observed.The effects of liquid–liquid two-phase flow rate ratio,and dimensionless number on droplet size,and pressure drop were investigated.The squeezing regime was mapped for 0.0005≤Ca_(c)≤0.0052(capillary number)and 0.1≤q≤10(flow rate ratio).The pressure drops of slugs were in the range from 40 Pa to 200 Pa.The slug lengths were measured between 1 mm and 9 mm.A universal flow map dependent on Ca_(c)Re_(d)^(0.5) are projected to investigate the droplet formation behavior in T-junction microchannel.Correlation expressions are proposed to predict pressure drops and the slug lengths for liquid–liquid two-phase flow in a square T-junction microchannel,using dimensionless numbers such as flow rate ratio and capillary number.The result shows that large continuous phase flow rates facilitate smaller slugs,whereas higher dispersed phase flow rates result in longer shorts.
基金co-supported by the National Key Research and Development Program of China(2021YFC3000604)the Key Program of the National Natural Science Foundation of China(42130101)。
文摘The Sichuan-Yunnan region,located at the southeastern margin of the Qinghai-Xizang Plateau,serves as a key channel for the southeastward extrusion of plateau material.The characteristics of crustal deformation and the mechanisms of deep material flow have been central topics of interest in geoscience research.In this work,a three-dimensional viscoelastic-plastic finite element model including the upper and mid-lower crust was established,constrained by GNSS horizontal crustal velocity observations and incorporating maj or active faults and geophysical survey data to explore the contribution of mid-lower crustal flow to surface deformation and its coupling with faults.Comparison of modeling experiments shows that relying solely on boundary loading or uniform layering assumptions fails to reproduce the GNSS observed velocities.We introduce a mid-lower crustal low-velocity weak zone,derived from the latest seismic velocity structure models.The new model improves the fit to GNSS observations.Tests of different viscosity coefficients in the low-velocity zone indicate an optimal viscosity range of 7.5×10^(19)-1×10^(20)Pa·s.Vertical profiles reveal that mid-lower crustal material motion is mainly concentrated at depths of 20-40 km,forming localized channelized flow in low-velocity zone with a typical Poiseuille velocity profile which indicates a ductile,fluid-like behavior with the lowvelocity zone serving as primary pathways for deep material transport.The results further show that under the geometric constraints of upper-crustal faults,the mid-lower crustal flow contributes approximately 1-3 mm/a to surface deformation,primarily concentrated along major faults.This indicates that faults play a key role in constraining and modulating the transmission of deep-seated dynamics to shallow surface deformation.However,the contribution of mid-lower crustal flow is also significant;neglecting its influence on surface deformation would lead to an incomplete understanding of the deformation pattern and bias the interpretation of block boundaries and crustal kinematic segmentation.
基金the support from the National Natural Science Foundation of China(No.52079077)Natural Science Foundation of Hubei Provincial(2025AFB358).
文摘Although extensive research has been conducted on CO_(2)-enhanced coalbed methane(CO_(2)-ECBM)recovery,most prior studies have focused on the impact of gas adsorption-induced swelling on coal permeability under equilibrium conditions.This paper presents a comprehensive thermo-hydro-mechanical-chemical(THMC)model that integrates thermal expansion and heat conduction(T),gas diffusion in the matrix and gas-water two-phase flow in the fractures(H),matrix and fracture deformation due to poroelasticity(M),and non-equilibrium binary gas adsorption-induced matrix swelling(C)during CO_(2)-ECBM recovery.The accuracy of the proposed model was verified through experimental data,and the model was simulated using finite element method(FEM)software.Simulation results indicate that the permeability evolution can be categorized into three stages.Ignoring the impact of water on gas adsorption properties would lead to an overestimation of the influence of adsorption-induced swelling,while disregarding non-equilibrium adsorption underestimates it.An examination of five designed cases identified critical factors influencing permeability.Parametric analysis shows that increases in the injection pressure,the injection temperature,and the initial permeability promote the competitive adsorption-induced swelling between CH_(4)and CO_(2),leading to increased CH_(4)production and CO_(2)storage.Conversely,an increase in initial water saturation reduces available gas flow space,decreasing both CH_(4)production and CO_(2)storage.Higher irreducible water saturation favors early gas recovery,while lower irreducible water saturation is more advantageous for long-term recovery.
文摘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.
基金National Natural Science Foundation of China(No.51175481)
文摘In order to speed underwater launch of minor-caliber weapons,a sealing device can be set in front of underwater muzzle to separate water,preventing the muzzle from water immersion.By establishing and simplifying the model of underwater weapon sealing device and unstructured mesh computing domain model based on computational fluid dynamics（CFD）,dynamic mesh and user defined function（UDF）,the N-S equation is solved and the numerical analysis and calculation of the complex two-phase flow inside the sealing device are carried out.The results show that the gas discharged from the sealing device is conducive to the formation of the projectile supercavity.When the projectile is launched at 5munder water,the shock wave before and after the projectile has impact on the box body up to 100 MPa,therefore the sealing device must be strong enough.The research results have the vital significance to the design of underwater weapon sealing device and the formation of the projectile supercavitation.
基金National Science Foundation of China,Grant/Award Number:5227090113Shenzhen Science and Technology Program,Grant/Award Numbers:KCXFZ20230731093901003,KCXFZ20211020163816023The Project of Hetao Shenzhen-Hong Kong Science and Technology Innovation Cooperation Zone,Grant/Award Number:HZQB-KCZYB-2020083。
文摘The two-phase flow in porous media is affected by multiple factors.In the present study,a two-dimensional numerical model of porous media was developed using the actual pore structure of the core sample.The phase field method was utilized to simulate the impact of displacement velocity,the water-gas viscosity ratio,and the density ratio on the flow behavior of two-phase fluids in porous media.The effectiveness of displacement was evaluated by analyzing CO_(2)saturation levels.The results indicate that the saturation of CO_(2)in porous media increased as the displacement velocity increased.When the displacement velocity exceeded 0.01 m/s,there was a corresponding increase in CO_(2)saturation.Conversely,when the displacement velocity was below this threshold,the impact on CO_(2)saturation was minimal.An“inflection point,”M3,was present in the viscosity ratio.When the viscosity of CO_(2)is less than 8.937×10^(-5)Pa·s(viscosity ratio below M3),variations in the viscosity of CO_(2)had little impact on its saturation.Conversely,when the viscosity of CO_(2)exceeded 8.937×10^(-5)Pa·s(viscosity ratio greater than M3),saturation increased with an increase in the viscosity ratio.In terms of the density ratio,the saturation of CO_(2)increased monotonically with an increase in the density ratio.Similarly,increasing density ratios resulted in a monotonic increase in CO_(2)saturation,though this trend was less pronounced in numerical simulations.Analysis results of displacement within dead-end pores using pressure and velocity diagrams reveal eddy currents as contributing factors.Finally,the impact of pore throat structure on the formation of dominant channels was examined.
基金Project(51375498)supported by the National Natural Science Foundation of China
文摘Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of unsteady flow characteristics on solid-liquid two-phase flow and pump performance were researched under design condition. The results show that clocking effect has a significant influence on the flow in pump, and the fluctuation of flow velocity and pressure is obvious, particularly near the volute tongue, at the position of small sections of volute and within diffuser. Clocking effect has a more influence on liquid-phase than on solid-phase, and the wake-jet structure of relative velocity of solid-phase is less obvious than liquid-phase near the volute tongue and the impeller passage outlet. The fluctuation of relative velocity of solid-phase flow is 7.6% smaller than liquid-phase flow at the impeller outlet on circular path. Head and radial forces of the impeller are 8.1% and 85.7% of fluctuation, respectively. The results provide a theoretical basis for further research for turbulence, improving efficient, reducing the hydraulic losses and wear. Finally, field tests were carried out to verify the operation and wear of slurry pump.
文摘A mathematical model is set to evaluate the 3-D dense solid-liquid two-phaseturbulent flow in a non-clogging mud pump, the flow feature in the impeller channel is simulatedwith the tool of IPSA. Meanwhile, resort to TECPLOT as the post-processor, the simulation results isvisualized. The results show the main flow characteristics: There exists backflow and aberrantvelocities at inlet area and a relative velocity slip between two phases; A jet-wake flow pattern isdiscerned around the shroud-suction side area; The relative velocity vector of solid phase iscloser to the pressure surface than that of liquid phase and the trend is more obvious with theincrease of diameter; The kinetic energy of turbulence k and the dissipation rate e reach theirpeaks at the corner of pressure and suction surface. The simulation results show a good agreementwith the experimental flow features in the impeller channel, which prove the turbulent model used isvalid and provide a theoretical design basis to non-clogging pumps.
基金the National Natural Science Foundation of China(50779069 and 90510007)the Start-up Scientific Research Foundation of China Agricultural University(2006021)the Beijing Natural Science Foundation(3071002).
文摘An improved large eddy simulation using a dynamic second-order sub-grid-scale （SGS） stress model has been developed to model the governing equations of dense turbulent particle-liquid two-phase flows in a rotating coordinate system, and continuity is conserved by a mass-weighted method to solve the filtered governing equations. In the cur- rent second-order SGS model, the SGS stress is a function of both the resolved strain-rate and rotation-rate tensors, and the model parameters are obtained from the dimensional consistency and the invariants of the strain-rate and the rotation-rate tensors. In the numerical calculation, the finite volume method is used to discretize the governing equations with a staggered grid system. The SIMPLEC algorithm is applied for the solution of the discretized governing equations. Body- fitted coordinates are used to simulate the two-phase flows in complex geometries. Finally the second-order dynamic SGS model is successfully applied to simulate the dense turbu-lent particle-liquid two-phase flows in a centrifugal impeller. The predicted pressure and velocity distributions are in good agreement with experimental results.
基金supported by National Natural Science Foundation of China (Grant No. 21076198)Zhejiang Provincial Natural Science Foundation of China (Granted No. R1100530)National Basic Research Program of China (973 Program,Grant No. 2009CB724303)
文摘The flow with solid-liquid two-phase media inside centrifugal pumps is very complicated and the relevant method for the hydraulic design is still immature so far. There exist two main problems in the operation of the two-phase flow pumps, i.e., low overall efficiency and severe abrasion. In this study, the three-dimensional, steady, incompressible, and turbulent solid-liquid two-phase flows in a low-specific-speed centrifugal pump are numerically simulated and analyzed by using a computational fluid dynamics （CFD） code based on the mixture model of the two-phase flow and the RNG k-~ two-equation turbulence model, in which the influences of rotation and curvature are fully taken into account. The coupling between impeller and volute is implemented by means of the frozen rotor method. The simulation results predicted indicate that the solid phase properties in two-phase flow, especially the concentration, the particle diameter and the density, have strong effects on the hydraulic performance of the pump. Both the pump head and the efficiency are reduced with increasing particle diameter or concentration. However, the effect of particle density on the performance is relatively minor. An obvious jet-wake flow structure is presented near the volute tongue and becomes more remarkable with increasing solid phase concentration. The suction side of the blade is subject to much more severe abrasion than the pressure side. The obtained results preliminarily reveal the characteristics of solid-liquid two-phase flow in the centrifugal pump, and are helpful for improvement and empirical correction in the hydraulic design of centrifugal pumps.
文摘Streamline simulation is developed to simulate waterflooding in fractured reservoirs. Conventional reservoir simulation methods for fluid flow simulation in large and complex reservoirs are very costly and time consuming. In streamline method, transport equations are solved on one-dimensional streamlines to reduce the computation time with less memory for simulation. First, pressure equation is solved on an Eulerian grid and streamlines are traced. Defining the "time of flight", saturation equations are mapped and solved on streamlines. Finally, the results are mapped back on Eulerian grid and the process is repeated until the simulation end time. The waterflooding process is considered in a fractured reservoir using the dual porosity model. Afterwards, a computational code is developed to solve the same problem by the IMPES method and the results of streamline simulation are compared to those of the IMPES and a commercial software. Finally, the accuracy and efficiency of streamline simulator for simulation of two-phase flow in fractured reservoirs has been proved.
基金Supported by the National 863 Project (2001AA642030-1) and Zhejiang Provincial Key Research Project (010007037).
文摘The flow field of gas and liquid in a φ150mm rotating-stream-tray (RST) scrubber is simulated by using computational fluid dynamic (CFD) method. The sismulation is based on the two-equation RNG κ-ε turbulence model, Eulerian multiphase model, and a real-shape 3D model with a huge number of meshes. The simulation results include detailed information about velocity, pressure, volume fraction and so on. Some features of the flow field are obtained: liquid is atomized in a thin annular zone; a high velocity air zone prevents water drops at the bottom from flying towards the wall; the pressure varies sharply at the end of blades and so on. The results will be helpful for structure optimization and engineering design.
基金supported by the Fund of Innovation Research Group of National Natural Science Foundation of China (Grant NO.5052160450323001)Major Program of National Natural Science Foundation of China (Grant No.50536020)
文摘Chokes are one of the most important components of downhole flow-control equipment. The particle erosion mathematical model, which considers particle-particle interaction, was established and used to simulate solid particle movement as well as particle erosion characteristics of the solid-liquid two-phase flow in a choke. The corresponding erosion reduction approach by setting ribs on the inner wall of the choke was advanced. This mathematical model includes three parts： the flow field simulation of the continuous carrier fluid by an Eulerian approach, the particle interaction simulation using the discrete particle hard sphere model by a Lagrangian approach and calculation of erosion rate using semiempirical correlations. The results show that particles accumulated in a narrow region from inlet to outlet of the choke and the dominating factor affecting particle motion is the fluid drag force. As a result, the optimization of rib geometrical parameters indicates that good anti-erosion performance can be achieved by four ribs, each of them with a height （H） of 3 mm and a width （B） of 5 mm equaling the interval between ribs （L）.
基金Financial support from the NSFC-ICIMOD (41661144041)Key Research and Development Program (2017SZ0041)Sichuan Province Science and Technology Support Project (2016SZ0067)
文摘Debris flows are typical two-phase flows, which commonly accompany erosion in mountainous areas, and may destroy bridge engineering by scouring. In this study, a physically-based two-phase model is applied for the simulation of debris flow scouring of bridge pier. In this model, the shear stress of debris flow on an erodible bed is considered to be a function of the solid shear stress, fluid shear stress, and volume fraction; accordingly, the erosion is incorporated into the two-phase model. Using a highaccuracy computational scheme based on the finite volume method, the model is employed for simulating a dynamic debris flow over an erodible bed. The numerical results are consistent with the experimental data, and verify the feasibility of the two-phase model. Moreover, a simple numerical test is performed to exhibit the fundamental behaviour of debris flow scouring of bridge pier, which shows that the degree of erosion on each side of the pier is higher compared to other areas. The scouring depth is influenced by the variations of solid volume fraction and velocity of debris flow and pier width.
