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.展开更多
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.展开更多
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）.展开更多
Liquid film cooling as an advanced cooling technology is widely used in space vehicles.Stable operation of liquid film along the rocket combustion inner wall is crucial for thermal protection of rocket engines.The sta...Liquid film cooling as an advanced cooling technology is widely used in space vehicles.Stable operation of liquid film along the rocket combustion inner wall is crucial for thermal protection of rocket engines.The stability of liquid film is mainly determined by the characteristics of interfacial wave,which is rarely investigated right now.How to improve the stability of thin film has become a hot spot.In view of this,an advanced model based on the conventional Volume of Fluid(VOF)model is adopted to investigate the characteristics of interfacial wave in gas-liquid flow by using OpenFOAM,and the mechanism of formation and development of wave is revealed intuitively through numerical study.The effects from gas velocity,surface tension and dynamic viscosity of liquid(three factors)on the wave are studied respectively.It can be found that the gas velocity is critical to the formation and development of wave,and four modes of droplets generation are illustrated in this paper.Besides,a gas vortex near the gas-liquid interface can induce formation of wave easily,so changing the gas vortex state can regulate formation and development of wave.What’s more,the change rules of three factors influencing on the interfacial wave are obtained,and the surface tension has a negative effect on the formation and development of wave only when the surface tension coefficient is above the critical value,whereas the dynamic viscosity has a positive effect in this process.Lastly,the maximum height and maximum slope angle of wave will level off as the gas velocity increases.Meanwhile,the maximum slope angle of wave is usually no more than 38°,no matter what happens to the three factors.展开更多
A mathematical modei of two-dimensional turbulent gas-particle twophase flow based on the modified diffusion flux modei (DFM) and a numerical simulation method to analyze the gas-particle flow structures are developed...A mathematical modei of two-dimensional turbulent gas-particle twophase flow based on the modified diffusion flux modei (DFM) and a numerical simulation method to analyze the gas-particle flow structures are developed. The modified diffusion flux modei, in which the acceleration due to various forces is taken into account for the calculation of the diffusion velocity of particles, is applicable to the analysis of multi-dimensional gas-particle two-phase turbulent flow. In order to verify its accuracy and efficiency, the numerical simulation by DFM is compared with experimental studies and the prediction by k-ε-kp two-fluid modei, which shows a reasonable agreement. It is confirmed that the modified diffusion flux modei is suitable for simulating the multi-dimensional gas-particle two-phase flow.展开更多
Pressure drop and liquid hold-up are two very important fluid flow parameters in design and control of multiphase flow pipelines.Friction factors play an important role in the accurate calculation of pressure drop.Var...Pressure drop and liquid hold-up are two very important fluid flow parameters in design and control of multiphase flow pipelines.Friction factors play an important role in the accurate calculation of pressure drop.Various empirical and semi-empirical closure relations exist in the literature to calculate the liquid-wall,gas-wall and interfacial friction in two-phase pipe flow.However most of them are empirical correlations found under special experimental conditions.In this paper by modification of a friction model available in the literature,an improved semiempirical model is proposed.The proposed model is incorporated in the two-fluid correlations under equilibrium conditions and solved.Pressure gradient and velocity profiles are validated against experimental data.Using the improved model,the pressure gradient deviation from experiments diminishes by about 3%;the no-slip condition at the interface is satisfied and the velocity profile is predicted in better agreement with the experimental data.展开更多
Instead of the capillary plasma generator(CPG),a discharge rod plasma generator(DRPG)is used in the30 mm electrothermal-chemical(ETC)gun to improve the ignition uniformity of the solid propellant.An axisymmetric two-d...Instead of the capillary plasma generator(CPG),a discharge rod plasma generator(DRPG)is used in the30 mm electrothermal-chemical(ETC)gun to improve the ignition uniformity of the solid propellant.An axisymmetric two-dimensional interior ballistics model of the solid propellant ETC gun(2D-IB-SPETCG)is presented to describe the process of the ETC launch.Both calculated pressure and projectile muzzle velocity accord well with the experimental results.The feasibility of the 2D-IB-SPETCG model is proved.Depending on the experimental data and initial parameters,detailed distribution of the ballistics parameters can be simulated.With the distribution of pressure and temperature of the gas phase and the propellant,the influence of plasma during the ignition process can be analyzed.Because of the radial flowing plasma,the propellant in the area of the DRPG is ignited within 0.01 ms,while all propellant in the chamber is ignited within 0.09 ms.The radial ignition delay time is much less than the axial delay time.During the ignition process,the radial pressure difference is less than 5 MPa at the place 0.025 m away from the breech.The radial ignition uniformity is proved.The temperature of the gas increases from several thousand K(conventional ignition)to several ten thousand K(plasma ignition).Compare the distribution of the density and temperature of the gas,we know that low density and high temperature gas appears near the exits of the DRPG,while high density and low temperature gas appears at the wall near the breech.The simulation of the 2D-IB-SPETCG model is an effective way to investigate the interior ballistics process of the ETC launch.The 2D-IB-SPETC model can be used for prediction and improvement of experiments.展开更多
On August 30, 2020, a high-intensity storm that dropped 45.4 mm of rain in 5 hours hit the Heixiluo basin and triggered a landslide-generated debris flow event, causing fatalities and damage. The original source of th...On August 30, 2020, a high-intensity storm that dropped 45.4 mm of rain in 5 hours hit the Heixiluo basin and triggered a landslide-generated debris flow event, causing fatalities and damage. The original source of the debris flow was a large slope collapse on a steep hillside. The fallen debris mass was enlarged through sediment entrainment and slope collapse and ultimately buried a bridge at the gully entrance. Approximately 6.9×10;m;of material, including sediments and collapsed slope deposits in the gullies, was entrained, and the maximum erosion depth reached 17 m. A geomorphological analysis was initially performed based on a detailed field investigation to recognize the liquid and solid sources of the debris flow and the areas subjected to deposition and erosion. A map of the erosiondeposition distribution was obtained based on preand post-event DEMs. Using the rainfall estimated by the nearest rain gauge and the solid source estimated by the DEMs, runoff and debris flow propagation was simulated using a liquid-solid two-phase model that considers the effects of runoff and entrainment. The similarity between the estimated and simulated deposition-erosion volumes was satisfactory. The behaviour of debris flows captured in the simulation is broadly in line with the main features of the observed event.展开更多
Numerical simulations are performed both for the single airflow and air-PMtwo-phase flow in wall flow diesel participate filters (DPF) for the first time. The calculationdomain is divided into two regions. In. the inl...Numerical simulations are performed both for the single airflow and air-PMtwo-phase flow in wall flow diesel participate filters (DPF) for the first time. The calculationdomain is divided into two regions. In. the inlet and outlet flow channels, the simulations areperformed for the steady and laminar flow; In the porous filtration walls, the calculation model forflow in porous media is used. The Lagrange two-phase flow model is used to calculate the air-PMflow in DPF, for the dispersed phase (PM), its flow tracks are obtained by the integrating of theLagrange kinetic equation. The calculated velocity, pressure distribution and PM flow tracks in DPFare obtained, which exhibits the main flow characteristics in wall flow DPF and will be help for theoptimal design and performance prediction of wall flow DPF.展开更多
According to environmental and energy issues,renewable energy has been vigorously promoted.Now solar power is widely used in many areas but it is limited by the weather conditions and cannot work continuously.Heat sto...According to environmental and energy issues,renewable energy has been vigorously promoted.Now solar power is widely used in many areas but it is limited by the weather conditions and cannot work continuously.Heat storage is a considerable solution for this problem and thermochemical energy storage is the most promising way because of its great energy density and stability.However,this technology is not mature enough to be applied to the industry.The reactor is an important component in the thermochemical energy storage system where the charging and discharging process happens.In this paper,a spiral coil is proposed and used as a reactor in the thermochemical energy storage system.The advantages of the spiral coil include simple structure,small volume,and so on.To investigate the flow characteristics,the simulation was carried out based on energy-minimization multi-scale model(EMMS)and Eulerian two-phase model.CaCO_(3) particles were chosen as the reactants.Particle distribution was shown in the results.The gas initial velocity was set to 2 m·s^(-1),3 m·s^(-1),and 4 m·s^(-1).When the particles flowed in the coil,gravity,centrifugal force and drag force influenced their flow.With the Reynold numbers increasing,centrifugal and drag force got larger.Accumulation phenomenon existed in the coil and results showed with the gas velocity increasing,accumulation moved from the bottom to the outer wall of the coil.Besides,the accumulation phenomenon was stabilized whenφ>720°.Also due to the centrifugal force,a secondary flow formed,which means solid particles moved from the inside wall to the outside wall.This secondary flow could promote turbulence and mixing of particles and gas.In addition,when the particle volume fraction is reduced from 0.2 to 0.1,the accumulation at the bottom of the coil decreases,and the unevenness of the velocity distribution becomes larger.展开更多
A detached eddy simulation(DES) and a k-ε-based Reynolds-averaged Navier–Stokes(RANS) calculation on the co-current spray drying chamber is presented. The DES used here is based on the Spalart–Allmaras(SA) turbulen...A detached eddy simulation(DES) and a k-ε-based Reynolds-averaged Navier–Stokes(RANS) calculation on the co-current spray drying chamber is presented. The DES used here is based on the Spalart–Allmaras(SA) turbulence model, whereas the standard k-ε(SKE) was considered here for comparison purposes. Predictions of the mean axial velocity, temperature and humidity profile have been evaluated and compared with experimental measurements. The effects of the turbulence model on the predictions of the mean axial velocity, temperature and the humidity profile are most noticeable in the(highly anisotropic) spraying region. The findings suggest that DES provide a more accurate prediction(with error less than 5%) of the flow field in a spray drying chamber compared with RANS-based k-ε models. The DES simulation also confirmed the presence of anisotropic turbulent flow in the spray dryer from the analysis of the velocity component fluctuations and turbulent structure as illustrated by the Q-criterion.展开更多
Erosion in slurry pumps presents a persistent challenge in industrial applications.This study examines the erosion of the static components of a 150ZJ-C42 centrifugal slurry pump,currently in operation at a beneficiat...Erosion in slurry pumps presents a persistent challenge in industrial applications.This study examines the erosion of the static components of a 150ZJ-C42 centrifugal slurry pump,currently in operation at a beneficiation plant,under varying particle conditions.Utilizing high-precision three-dimensional reverse engineering,the pump’s flow passage geometry was reconstructed to facilitate detailed erosion analysis.Focusing on the front and rear baffles of the pump chamber,as well as the volute,erosion patterns were analyzed for different particle volume concentrations and sizes.The results reveal that the highest erosion damage consistently occurs near the volute tongue,with wear being most severe in regions adjacent to the partition plate near the rear cover.Erosion damage intensity in this area correlates positively with particle diameter.Notably,the average erosion rate in the volute surpasses that of the front and rear chamber liners,reaching a value as high as 6.03×10^(-7)kg·m^(-2)·s^(-1)at a particle concentration of 9%and diameter of 0.1 mm,adversely impacting pump stability.For the pump chamber baffles,increased erosion is observed at a particle diameter of 0.05 mm under constant volume concentration conditions,while higher particle concentrations exacerbate localized erosion.展开更多
In this paper, the stochastic particle-trajectory model is proposed for simulating the dynamic behavior of circulating fluidized bed (CFB). In our model, the motion of solid phase is obtained by calculating the indivi...In this paper, the stochastic particle-trajectory model is proposed for simulating the dynamic behavior of circulating fluidized bed (CFB). In our model, the motion of solid phase is obtained by calculating the individual particle trajectory while gas flow is obtained by solving the Navier-Stokes Equation including two-phase interaction. For the calculation of solid phase, the motion of each particle is decomposed into a collision process and a suspension process. In suspension process, the less important and/or unclear forces are described as a random force considering gravity, drag force and pressure gradient. As a result, the proposed model gives some numerical simulations of CFB. It indicates that the stochastic particle-trajectory model can be used to simulate qualitatively the annulus-core structure of CFB and the influences of stochastic factors cannot be ignored. In a CFB, the coupling of stochastic factors between two phases makes the radial voidage decreased. Moreover, the upward motion of particles is mitigated by both stochastic factors and turbulence between two phases.展开更多
Under the absolute coordinate system, the unsteady Reynolds averaged Navier-Stokes(URANS)equations and the k-ω SST turbulence model are solved using the finite volume method to simulate the aerodynamic characteristic...Under the absolute coordinate system, the unsteady Reynolds averaged Navier-Stokes(URANS)equations and the k-ω SST turbulence model are solved using the finite volume method to simulate the aerodynamic characteristics of large seaplane flying with the ground-effect above wavy surface. The velocity inlet wave-making method and the volume of fluid model are used to accurately simulate the linear regular waves and to precisely capture the free surface. This paper studies the influence of the sideslip angle on the aerodynamic characteristics of large seaplane when it is cruising above wavy water. The results show that the wave surface mainly affects the pressure distribution on the lower surface of the wing. When the sideslip angle varies from 0° to 8°,the varying of frequency of aerodynamic forces is consistent with the wave encounter frequency,and both periods are 0.6 s. With the increase of the sideslip angle,the lift coefficient and the pitching moment coefficient decrease. However,when the sideslip angle is smaller than 4°,the decrease amplitude is small,and the significant decrease occurs above 4° and during the whole process of the change of sideslip angle,the aerodynamic fluctuation amplitude is almost unchanged. As the drag coefficient increases with the increase of sideslip angle,significant increase also occurs when the value is greater than4°,and the fluctuation amplitude does not show any correlations. The rolling moment coefficient and yaw moment coefficient increase with the increase of the sideslip angle,and the fluctuation amplitudes of both increase linearly with the increase of the sideslip angle.展开更多
Flood discharge over top outlet of dam is simulated by 2 dimension water air two phase mathematical model. Distribution of dynamic pressure, turbulent kinetic energy (k), turbulent dissipation rate (ε), free wat...Flood discharge over top outlet of dam is simulated by 2 dimension water air two phase mathematical model. Distribution of dynamic pressure, turbulent kinetic energy (k), turbulent dissipation rate (ε), free water surface and velocity field have been obtained. The simulated results were tested by physical model, which shows that the computed results are identical with that of the physical model.展开更多
针对复杂管道系统内水流冲击滞留气团现象,采用VOF模型(Volume of Fluid Model)进行了数值模拟计算,并与一维模型进行了比较计算分析,结果表明:系统的最大压力并不总是气团的最大压力,有可能还会叠加水体对管壁的撞击而形成的突然升高...针对复杂管道系统内水流冲击滞留气团现象,采用VOF模型(Volume of Fluid Model)进行了数值模拟计算,并与一维模型进行了比较计算分析,结果表明:系统的最大压力并不总是气团的最大压力,有可能还会叠加水体对管壁的撞击而形成的突然升高压力。与实验实测结果的比较分析表明:采用VOF模型,能够较精细地仿真水流冲击滞留气团现象的气团形态、流场结构以及压力分布等的变化过程,其压力数值计算结果与实验实测基本吻合,其计算误差明显小于现有一维模型的计算误差,是深入研究该复杂瞬变流现象的有效方法。展开更多
基金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 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 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）.
文摘Liquid film cooling as an advanced cooling technology is widely used in space vehicles.Stable operation of liquid film along the rocket combustion inner wall is crucial for thermal protection of rocket engines.The stability of liquid film is mainly determined by the characteristics of interfacial wave,which is rarely investigated right now.How to improve the stability of thin film has become a hot spot.In view of this,an advanced model based on the conventional Volume of Fluid(VOF)model is adopted to investigate the characteristics of interfacial wave in gas-liquid flow by using OpenFOAM,and the mechanism of formation and development of wave is revealed intuitively through numerical study.The effects from gas velocity,surface tension and dynamic viscosity of liquid(three factors)on the wave are studied respectively.It can be found that the gas velocity is critical to the formation and development of wave,and four modes of droplets generation are illustrated in this paper.Besides,a gas vortex near the gas-liquid interface can induce formation of wave easily,so changing the gas vortex state can regulate formation and development of wave.What’s more,the change rules of three factors influencing on the interfacial wave are obtained,and the surface tension has a negative effect on the formation and development of wave only when the surface tension coefficient is above the critical value,whereas the dynamic viscosity has a positive effect in this process.Lastly,the maximum height and maximum slope angle of wave will level off as the gas velocity increases.Meanwhile,the maximum slope angle of wave is usually no more than 38°,no matter what happens to the three factors.
基金Special Funds for Major State Basic Research Projects of China(G1999022200)
文摘A mathematical modei of two-dimensional turbulent gas-particle twophase flow based on the modified diffusion flux modei (DFM) and a numerical simulation method to analyze the gas-particle flow structures are developed. The modified diffusion flux modei, in which the acceleration due to various forces is taken into account for the calculation of the diffusion velocity of particles, is applicable to the analysis of multi-dimensional gas-particle two-phase turbulent flow. In order to verify its accuracy and efficiency, the numerical simulation by DFM is compared with experimental studies and the prediction by k-ε-kp two-fluid modei, which shows a reasonable agreement. It is confirmed that the modified diffusion flux modei is suitable for simulating the multi-dimensional gas-particle two-phase flow.
基金supported by the Iran National Science Foundation(Grant 96006257)。
文摘Pressure drop and liquid hold-up are two very important fluid flow parameters in design and control of multiphase flow pipelines.Friction factors play an important role in the accurate calculation of pressure drop.Various empirical and semi-empirical closure relations exist in the literature to calculate the liquid-wall,gas-wall and interfacial friction in two-phase pipe flow.However most of them are empirical correlations found under special experimental conditions.In this paper by modification of a friction model available in the literature,an improved semiempirical model is proposed.The proposed model is incorporated in the two-fluid correlations under equilibrium conditions and solved.Pressure gradient and velocity profiles are validated against experimental data.Using the improved model,the pressure gradient deviation from experiments diminishes by about 3%;the no-slip condition at the interface is satisfied and the velocity profile is predicted in better agreement with the experimental data.
文摘Instead of the capillary plasma generator(CPG),a discharge rod plasma generator(DRPG)is used in the30 mm electrothermal-chemical(ETC)gun to improve the ignition uniformity of the solid propellant.An axisymmetric two-dimensional interior ballistics model of the solid propellant ETC gun(2D-IB-SPETCG)is presented to describe the process of the ETC launch.Both calculated pressure and projectile muzzle velocity accord well with the experimental results.The feasibility of the 2D-IB-SPETCG model is proved.Depending on the experimental data and initial parameters,detailed distribution of the ballistics parameters can be simulated.With the distribution of pressure and temperature of the gas phase and the propellant,the influence of plasma during the ignition process can be analyzed.Because of the radial flowing plasma,the propellant in the area of the DRPG is ignited within 0.01 ms,while all propellant in the chamber is ignited within 0.09 ms.The radial ignition delay time is much less than the axial delay time.During the ignition process,the radial pressure difference is less than 5 MPa at the place 0.025 m away from the breech.The radial ignition uniformity is proved.The temperature of the gas increases from several thousand K(conventional ignition)to several ten thousand K(plasma ignition).Compare the distribution of the density and temperature of the gas,we know that low density and high temperature gas appears near the exits of the DRPG,while high density and low temperature gas appears at the wall near the breech.The simulation of the 2D-IB-SPETCG model is an effective way to investigate the interior ballistics process of the ETC launch.The 2D-IB-SPETC model can be used for prediction and improvement of experiments.
基金supported by the CAS Original Innovation Program(Grant No.ZDBS-LY-DQC039)National Natural Science Foundation of China(Grant No.41907241)+1 种基金Foundation of Youth Innovation Promotion Association,Chinese Academy of Sciences(Grant No.2021373)CAS“Light of West China”Program and the Foundation for Young Scientists of the Institute of Mountain Hazards and Environment,CAS(Grant No.SDS-QN-1901)。
文摘On August 30, 2020, a high-intensity storm that dropped 45.4 mm of rain in 5 hours hit the Heixiluo basin and triggered a landslide-generated debris flow event, causing fatalities and damage. The original source of the debris flow was a large slope collapse on a steep hillside. The fallen debris mass was enlarged through sediment entrainment and slope collapse and ultimately buried a bridge at the gully entrance. Approximately 6.9×10;m;of material, including sediments and collapsed slope deposits in the gullies, was entrained, and the maximum erosion depth reached 17 m. A geomorphological analysis was initially performed based on a detailed field investigation to recognize the liquid and solid sources of the debris flow and the areas subjected to deposition and erosion. A map of the erosiondeposition distribution was obtained based on preand post-event DEMs. Using the rainfall estimated by the nearest rain gauge and the solid source estimated by the DEMs, runoff and debris flow propagation was simulated using a liquid-solid two-phase model that considers the effects of runoff and entrainment. The similarity between the estimated and simulated deposition-erosion volumes was satisfactory. The behaviour of debris flows captured in the simulation is broadly in line with the main features of the observed event.
文摘Numerical simulations are performed both for the single airflow and air-PMtwo-phase flow in wall flow diesel participate filters (DPF) for the first time. The calculationdomain is divided into two regions. In. the inlet and outlet flow channels, the simulations areperformed for the steady and laminar flow; In the porous filtration walls, the calculation model forflow in porous media is used. The Lagrange two-phase flow model is used to calculate the air-PMflow in DPF, for the dispersed phase (PM), its flow tracks are obtained by the integrating of theLagrange kinetic equation. The calculated velocity, pressure distribution and PM flow tracks in DPFare obtained, which exhibits the main flow characteristics in wall flow DPF and will be help for theoptimal design and performance prediction of wall flow DPF.
基金the financial support provided by Natural Science Foundation of Jiangsu Province (BK20180936)the Initial Funding of Scientific Research for the Introduction of Talents (YJ2021-41)
文摘According to environmental and energy issues,renewable energy has been vigorously promoted.Now solar power is widely used in many areas but it is limited by the weather conditions and cannot work continuously.Heat storage is a considerable solution for this problem and thermochemical energy storage is the most promising way because of its great energy density and stability.However,this technology is not mature enough to be applied to the industry.The reactor is an important component in the thermochemical energy storage system where the charging and discharging process happens.In this paper,a spiral coil is proposed and used as a reactor in the thermochemical energy storage system.The advantages of the spiral coil include simple structure,small volume,and so on.To investigate the flow characteristics,the simulation was carried out based on energy-minimization multi-scale model(EMMS)and Eulerian two-phase model.CaCO_(3) particles were chosen as the reactants.Particle distribution was shown in the results.The gas initial velocity was set to 2 m·s^(-1),3 m·s^(-1),and 4 m·s^(-1).When the particles flowed in the coil,gravity,centrifugal force and drag force influenced their flow.With the Reynold numbers increasing,centrifugal and drag force got larger.Accumulation phenomenon existed in the coil and results showed with the gas velocity increasing,accumulation moved from the bottom to the outer wall of the coil.Besides,the accumulation phenomenon was stabilized whenφ>720°.Also due to the centrifugal force,a secondary flow formed,which means solid particles moved from the inside wall to the outside wall.This secondary flow could promote turbulence and mixing of particles and gas.In addition,when the particle volume fraction is reduced from 0.2 to 0.1,the accumulation at the bottom of the coil decreases,and the unevenness of the velocity distribution becomes larger.
基金Supported by the Ministry of Education Malaysia through RACE(RDU121308)and FRGS(RDU130136)
文摘A detached eddy simulation(DES) and a k-ε-based Reynolds-averaged Navier–Stokes(RANS) calculation on the co-current spray drying chamber is presented. The DES used here is based on the Spalart–Allmaras(SA) turbulence model, whereas the standard k-ε(SKE) was considered here for comparison purposes. Predictions of the mean axial velocity, temperature and humidity profile have been evaluated and compared with experimental measurements. The effects of the turbulence model on the predictions of the mean axial velocity, temperature and the humidity profile are most noticeable in the(highly anisotropic) spraying region. The findings suggest that DES provide a more accurate prediction(with error less than 5%) of the flow field in a spray drying chamber compared with RANS-based k-ε models. The DES simulation also confirmed the presence of anisotropic turbulent flow in the spray dryer from the analysis of the velocity component fluctuations and turbulent structure as illustrated by the Q-criterion.
基金The authors gratefully acknowledge the filnancial support of the National Natural Science Foundation of China(Grant No.52369018)the Major Training Program of University Research and Innovation Platform of Gansu Provincial Department of Education(No.2024CXPT-09)+1 种基金the Administration of Central Funds Guiding the Local Science and Technology Development,China(Grant No.23ZYQA0320)the Double First-Class Key Program of Gansu Provincial Department of Education,Grant No.GCJ2022-38.
文摘Erosion in slurry pumps presents a persistent challenge in industrial applications.This study examines the erosion of the static components of a 150ZJ-C42 centrifugal slurry pump,currently in operation at a beneficiation plant,under varying particle conditions.Utilizing high-precision three-dimensional reverse engineering,the pump’s flow passage geometry was reconstructed to facilitate detailed erosion analysis.Focusing on the front and rear baffles of the pump chamber,as well as the volute,erosion patterns were analyzed for different particle volume concentrations and sizes.The results reveal that the highest erosion damage consistently occurs near the volute tongue,with wear being most severe in regions adjacent to the partition plate near the rear cover.Erosion damage intensity in this area correlates positively with particle diameter.Notably,the average erosion rate in the volute surpasses that of the front and rear chamber liners,reaching a value as high as 6.03×10^(-7)kg·m^(-2)·s^(-1)at a particle concentration of 9%and diameter of 0.1 mm,adversely impacting pump stability.For the pump chamber baffles,increased erosion is observed at a particle diameter of 0.05 mm under constant volume concentration conditions,while higher particle concentrations exacerbate localized erosion.
基金Supported by the National Natural Science Foundation of China (No. 29976024) Natural Science Foundation of Shaanxi Province (No. 2001SL12) Research Foundation of Multi-Phase Reaction Laboratory in Chinese Academy of Sciences.
文摘In this paper, the stochastic particle-trajectory model is proposed for simulating the dynamic behavior of circulating fluidized bed (CFB). In our model, the motion of solid phase is obtained by calculating the individual particle trajectory while gas flow is obtained by solving the Navier-Stokes Equation including two-phase interaction. For the calculation of solid phase, the motion of each particle is decomposed into a collision process and a suspension process. In suspension process, the less important and/or unclear forces are described as a random force considering gravity, drag force and pressure gradient. As a result, the proposed model gives some numerical simulations of CFB. It indicates that the stochastic particle-trajectory model can be used to simulate qualitatively the annulus-core structure of CFB and the influences of stochastic factors cannot be ignored. In a CFB, the coupling of stochastic factors between two phases makes the radial voidage decreased. Moreover, the upward motion of particles is mitigated by both stochastic factors and turbulence between two phases.
文摘Under the absolute coordinate system, the unsteady Reynolds averaged Navier-Stokes(URANS)equations and the k-ω SST turbulence model are solved using the finite volume method to simulate the aerodynamic characteristics of large seaplane flying with the ground-effect above wavy surface. The velocity inlet wave-making method and the volume of fluid model are used to accurately simulate the linear regular waves and to precisely capture the free surface. This paper studies the influence of the sideslip angle on the aerodynamic characteristics of large seaplane when it is cruising above wavy water. The results show that the wave surface mainly affects the pressure distribution on the lower surface of the wing. When the sideslip angle varies from 0° to 8°,the varying of frequency of aerodynamic forces is consistent with the wave encounter frequency,and both periods are 0.6 s. With the increase of the sideslip angle,the lift coefficient and the pitching moment coefficient decrease. However,when the sideslip angle is smaller than 4°,the decrease amplitude is small,and the significant decrease occurs above 4° and during the whole process of the change of sideslip angle,the aerodynamic fluctuation amplitude is almost unchanged. As the drag coefficient increases with the increase of sideslip angle,significant increase also occurs when the value is greater than4°,and the fluctuation amplitude does not show any correlations. The rolling moment coefficient and yaw moment coefficient increase with the increase of the sideslip angle,and the fluctuation amplitudes of both increase linearly with the increase of the sideslip angle.
文摘Flood discharge over top outlet of dam is simulated by 2 dimension water air two phase mathematical model. Distribution of dynamic pressure, turbulent kinetic energy (k), turbulent dissipation rate (ε), free water surface and velocity field have been obtained. The simulated results were tested by physical model, which shows that the computed results are identical with that of the physical model.
文摘针对复杂管道系统内水流冲击滞留气团现象,采用VOF模型(Volume of Fluid Model)进行了数值模拟计算,并与一维模型进行了比较计算分析,结果表明:系统的最大压力并不总是气团的最大压力,有可能还会叠加水体对管壁的撞击而形成的突然升高压力。与实验实测结果的比较分析表明:采用VOF模型,能够较精细地仿真水流冲击滞留气团现象的气团形态、流场结构以及压力分布等的变化过程,其压力数值计算结果与实验实测基本吻合,其计算误差明显小于现有一维模型的计算误差,是深入研究该复杂瞬变流现象的有效方法。
