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 des...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.展开更多
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.展开更多
Air-bubble generator is the key part of the self-inspiration type swirl flotation machines,whose flow field structure has a great effect on flotation.The multiphase volume of fluid(VOF),standard k-ε turbulent model a...Air-bubble generator is the key part of the self-inspiration type swirl flotation machines,whose flow field structure has a great effect on flotation.The multiphase volume of fluid(VOF),standard k-ε turbulent model and the SIMPLE method were chosen to simulate the present model;the first order upwind difference scheme was utilized to perform a discrete solution for momentum equation.The distributing law of the velocity,pressure,turbulent kinetic energy of every section along the flow direction of air-bubble generator was analyzed.The results indicate that the bubbles are heavily broken up in the middle cross section of throat sect and the entrance of diffuser sect along the flow direction,and the turbulent kinetic energy of diffuser sect is larger than the entrance of throat sect and mixing chamber.展开更多
Dedicated experiments and numerical simulations have been conducted to investigate the splitting characteristics of a gas-liquid two phase flow at a T junction.The experiments were carried out for different gas-liquid...Dedicated experiments and numerical simulations have been conducted to investigate the splitting characteristics of a gas-liquid two phase flow at a T junction.The experiments were carried out for different gas-liquid velocities.The flow rates in the two branches were measured accurately to determine how the two considered phases distribute in the two outlets.The experimental results have shown that when the two outlet pressures are asymmetric,the two-phase flow always tends to flow into the outlet which has a lower pressure.As the inlet liquid velocity increases,however,the two-phase flow gradually tends to split evenly.Compared with the experiment results,the pressure difference between the two outlets can be determined more accurately by means of numerical simulation.The trends of experimental results and simulations are in very good agreement.展开更多
Erosion-corrosion of liquid-solid two-phase flow occurring in a pipe with sudden expansion in cross-section is numerically simulated in this paper. The global model for erosion-corrosion process includes three main co...Erosion-corrosion of liquid-solid two-phase flow occurring in a pipe with sudden expansion in cross-section is numerically simulated in this paper. The global model for erosion-corrosion process includes three main components: the liquid-solid two-phase flow model, erosion model and corrosion model. The Eulerian-Lagrangian approach is used to simulate liquid-solid two-phase flow, while the stochastic trajectory model was adopted to obtain properties of particle phase. Two-way coupling effect between the fluid and the particle phase is considered in the model. The accuracy of the models is tested by the data in the reference. The comparison shows that the model is basically correct and feasible.展开更多
The flow focusing nozzle is a new type of nozzle that performs effective atomization of the discrete phase by means of high-speed motion of the continuous phase.The flow pattern and its morphological changes have a si...The flow focusing nozzle is a new type of nozzle that performs effective atomization of the discrete phase by means of high-speed motion of the continuous phase.The flow pattern and its morphological changes have a significant effect on the atomization, but the influence of different parameters on the morphological change of the flow pattern remains unclear.The flow focusing pattern and morphological changes in the two-phase flow inside the nozzle were simulated numerically, based on the volume of fluid method.The results demonstrate that the ratio of the nozzle-to-capillary distance and capillary diameter, the gas–liquid velocity ratio, and capillary diameter have significant effects on the flow pattern.When the ratio of the nozzle-to-capillary distance H and capillary diameter D increases, or the capillary diameter D increases, the flow pattern tends to transform into a laminar form; however, when the gas–liquid velocity ratio V increases, the flow pattern tends to transform into a turbulence form.Furthermore, we define the cone-shaped expansion rate, cone-shaped focusing rate,and cone angle in order to study the morphological changes in the cone shape inside the nozzle.The results indicate that the morphological change of the cone shape and flow pattern transformation is interrelated.When the cone shape tends to be unstable, the flow pattern changes towards flow blurring, whereas, a stable cone indicates that the flow tends to exhibit a droplet pattern.展开更多
To investigate the movement mechanism of debris flow,a two-dimensional,two-phase,depthintegrated model is introduced.The model uses Mohr-Coulomb plasticity for the solid rheology,and the fluid stress is modeled as a N...To investigate the movement mechanism of debris flow,a two-dimensional,two-phase,depthintegrated model is introduced.The model uses Mohr-Coulomb plasticity for the solid rheology,and the fluid stress is modeled as a Newtonian fluid.The interaction between solid and liquid phases,which plays a major role in debris flow movement,is assumed to consist of drag and buoyancy forces.The applicability of drag force formulas is discussed.Considering the complex interaction between debris flow and the bed surface,a combined friction boundary condition is imposed on the bottom,and this is also discussed.To solve the complex model equations,a numerical method with second-order accuracy based on the finite volume method is proposed.Several numerical experiments are performed to verify the feasibilities of model and numerical schemes.Numerical results demonstrate that different solid volume fractions substantially affect debris flow movement.展开更多
文摘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.
文摘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.
基金Project supported by the Scientific Research Foundation of Hebei University of Technology of China
文摘Air-bubble generator is the key part of the self-inspiration type swirl flotation machines,whose flow field structure has a great effect on flotation.The multiphase volume of fluid(VOF),standard k-ε turbulent model and the SIMPLE method were chosen to simulate the present model;the first order upwind difference scheme was utilized to perform a discrete solution for momentum equation.The distributing law of the velocity,pressure,turbulent kinetic energy of every section along the flow direction of air-bubble generator was analyzed.The results indicate that the bubbles are heavily broken up in the middle cross section of throat sect and the entrance of diffuser sect along the flow direction,and the turbulent kinetic energy of diffuser sect is larger than the entrance of throat sect and mixing chamber.
基金the National Science and Technology Major Project of China(No.2016ZX05028-004-003).
文摘Dedicated experiments and numerical simulations have been conducted to investigate the splitting characteristics of a gas-liquid two phase flow at a T junction.The experiments were carried out for different gas-liquid velocities.The flow rates in the two branches were measured accurately to determine how the two considered phases distribute in the two outlets.The experimental results have shown that when the two outlet pressures are asymmetric,the two-phase flow always tends to flow into the outlet which has a lower pressure.As the inlet liquid velocity increases,however,the two-phase flow gradually tends to split evenly.Compared with the experiment results,the pressure difference between the two outlets can be determined more accurately by means of numerical simulation.The trends of experimental results and simulations are in very good agreement.
基金Supported by the National Natural Science Foundation of China(N.59831030).
文摘Erosion-corrosion of liquid-solid two-phase flow occurring in a pipe with sudden expansion in cross-section is numerically simulated in this paper. The global model for erosion-corrosion process includes three main components: the liquid-solid two-phase flow model, erosion model and corrosion model. The Eulerian-Lagrangian approach is used to simulate liquid-solid two-phase flow, while the stochastic trajectory model was adopted to obtain properties of particle phase. Two-way coupling effect between the fluid and the particle phase is considered in the model. The accuracy of the models is tested by the data in the reference. The comparison shows that the model is basically correct and feasible.
基金Supported by the National Natural Science Foundation of China(51776016,51606006)Beijing Natural Science Foundation(3172025,3182030)+2 种基金National Key Research and Development Program(2017YFB0103401)National Engineering Laboratory for Mobile Source Emission Control Technology(NELMS2017A10)the Talents Foundation of Beijing Jiaotong University(2018RC017)
文摘The flow focusing nozzle is a new type of nozzle that performs effective atomization of the discrete phase by means of high-speed motion of the continuous phase.The flow pattern and its morphological changes have a significant effect on the atomization, but the influence of different parameters on the morphological change of the flow pattern remains unclear.The flow focusing pattern and morphological changes in the two-phase flow inside the nozzle were simulated numerically, based on the volume of fluid method.The results demonstrate that the ratio of the nozzle-to-capillary distance and capillary diameter, the gas–liquid velocity ratio, and capillary diameter have significant effects on the flow pattern.When the ratio of the nozzle-to-capillary distance H and capillary diameter D increases, or the capillary diameter D increases, the flow pattern tends to transform into a laminar form; however, when the gas–liquid velocity ratio V increases, the flow pattern tends to transform into a turbulence form.Furthermore, we define the cone-shaped expansion rate, cone-shaped focusing rate,and cone angle in order to study the morphological changes in the cone shape inside the nozzle.The results indicate that the morphological change of the cone shape and flow pattern transformation is interrelated.When the cone shape tends to be unstable, the flow pattern changes towards flow blurring, whereas, a stable cone indicates that the flow tends to exhibit a droplet pattern.
基金Financial support from the NSFC-ICIMOD(41661144041)NSFC(Grant No.41772312)+1 种基金Key Research and Development Program(2017SZ0041)Sichuan Province Science and Technology Support Project(2016SZ0067)
文摘To investigate the movement mechanism of debris flow,a two-dimensional,two-phase,depthintegrated model is introduced.The model uses Mohr-Coulomb plasticity for the solid rheology,and the fluid stress is modeled as a Newtonian fluid.The interaction between solid and liquid phases,which plays a major role in debris flow movement,is assumed to consist of drag and buoyancy forces.The applicability of drag force formulas is discussed.Considering the complex interaction between debris flow and the bed surface,a combined friction boundary condition is imposed on the bottom,and this is also discussed.To solve the complex model equations,a numerical method with second-order accuracy based on the finite volume method is proposed.Several numerical experiments are performed to verify the feasibilities of model and numerical schemes.Numerical results demonstrate that different solid volume fractions substantially affect debris flow movement.
