The fluid temperature statistics along particle trajectories is crucial to under-stand the mechanisms of turbulent non-isothermal or reactive fluid-particle flow,especially for the Lagrangian model of non-isothermal p...The fluid temperature statistics along particle trajectories is crucial to under-stand the mechanisms of turbulent non-isothermal or reactive fluid-particle flow,especially for the Lagrangian model of non-isothermal particle-laden turbulent flow.In the present study,direct numerical simulations were utilized to generate temperature field statistics in particle-laden incompressible stationary homogeneous isotropic turbulent flows,which is focused on the effect of particle response time on the Lagrangian statistics of the particle and the fluid temperature seen by particles.It shows that,for the particles withτp/τk<1,the ratio of the fluid intensity seen by particle to fluid temperature intensities deceased asτp/τk increased;while for larger particles(τp/τk>1),the trend is inversed.For small parti-cles(τp/τk<5),the Lagrangian autocorrelation coefficient of the particle temperature R_(p)^(T)decreases as the particle inertia(τp/τk)increases.The trend is reversed for larger particles.The autocorrelation of fluid temperature along the particle path,R_(pf)^(T),decreased as the particle inertia increased.And as the particle inertia increased,the autocorrelation coeffi-cient of the fluid temperature seen by particle decreased more rapidly than that of the particle temperature.The mean temperature gradient contributes to the correlation be-tween the particles velocity component and temperature fluctuations in the direction of the gradient.For the particles withτp/τk<1,the magnitude of the correlation coefficient in-creases as the particle inertia increases,while this value is independent of the particle time constant for larger particles.展开更多
Particle-laden flows in a horizontal channel were investigated by means of a two-phase particle image velocimetry （PIV） technique. Experiments were performed at a Reynolds number of 6826 and the flow is seeded with ...Particle-laden flows in a horizontal channel were investigated by means of a two-phase particle image velocimetry （PIV） technique. Experiments were performed at a Reynolds number of 6826 and the flow is seeded with polythene beads of two sizes, 60μm and 110μm. One was slightly smaller than and the other was larger than the Kolmogorov length scale. The particle loadings were relatively low, with mass loading ratio ranging from 5 ×10^-4 to 4 × 10^-2 and volume fractions from 6×10×-7 to 4.8×10^-5, respectively. The results show that the presence of particles can dramatically modify the turbulence even under the lowest mass loading ratio of 5 × 10^-4. The mean flow is attenuated and de- creased with increasing particle size and mass loading. The turbulence intensities are enhanced in all the cases concerned. With the increase of the mass loading, the intensities vary in a complicated manner in the case of small particles, indicating complicated particle-turbulence interactions; whereas they increase monotonously in the case of large particles. The particle velocities and concentrations are also given. The particles lag behind the fluid in the center region but lead in the wall region, and this trend is more prominent for the large particles. The streamwise particle fluctuations are larger than the gas fluctuations for both sizes of particles, however their varying trend with the mass loadings is not so clear. The wallnormal fluctuations increase with increasing mass loadings. They are smaller in the 60μm particle case but larger in the 110μm particle case than those of the gas phase. It seems that the small particles follow the fluid motion to certain extent while the larger particles are more likely dominated by their own inertia. Finally, remarkable non-uniform distributions of particle concentration are observed, especially for the large particles. The inertia of particles is proved to be very important for the turbulence modification and particles behaviors and thus should be considered in horizonta展开更多
Fully developed vertical turbulent channel flow with particle transport wasinvestigated by use of Large Eddy Simulation (LES) approach coupled with dynamic the Sub-Grid Scale(SGS) model. It was assumed that the motion...Fully developed vertical turbulent channel flow with particle transport wasinvestigated by use of Large Eddy Simulation (LES) approach coupled with dynamic the Sub-Grid Scale(SGS) model. It was assumed that the motion of each particle is followed in a Lagrangian frame ofreference driven by the forces exerted by fluid motion and gravity under the condition of one-waycoupling. The goal of this study is to investigate the effectiveness of the LES technique forpredicting particle transport in turbulent flow and the behavior of particle-laden turbulent channelflow for three kinds of particles at different Stokes numbers. To depict the behavior ofparticle-laden turbulent channel flow, statistical quantities including particle fluctuation andfluid-particle velocity correlation, and visualization of the particle number density field wereanalyzed.展开更多
The dynamic and thermal performance of particle-laden turbulent flow is investigated via direction numerical simulation combined with the Lagrangian point-particle tracking under the condition of two-way coupling, wit...The dynamic and thermal performance of particle-laden turbulent flow is investigated via direction numerical simulation combined with the Lagrangian point-particle tracking under the condition of two-way coupling, with a focus on the contributions of particle feedback effect to momentum and heat transfer of turbulence. We take into account the effects of particles on flow drag and Nusselt number and explore the possibility of drag reduction in conjunction with heat transfer enhancement in particle-laden turbulent flows. The effects of particles on momentum and heat transfer are analyzed, and the possibility of drag reduction in conjunction with heat transfer enhancement for the prototypical case of particle-laden turbulent channel flows is addressed. We present results of turbulence modification and heat transfer in turbulent particle-laden channel flow, which shows the heat transfer reduction when large inertial particles with low specific heat capacity are added to the flow. However, we also found an enhancement of the heat transfer and a small reduction of the flow drag when particles with high specific heat capacity are involved. The present results show that particles, which are active agents, interact not only with the velocity field, but also the temperature field and can cause a dissimilarity in momentum and heat transport. This demonstrates that the possibility to increase heat transfer and suppress friction drag can be achieved with addition of particles with different thermal properties.展开更多
This paper presents the development and validation of a fully coupled computational fluid dynamics—discrete element method—volume of fluid(CFD-DEM-VOF)model to simulate the complex behavior of particle-laden flows w...This paper presents the development and validation of a fully coupled computational fluid dynamics—discrete element method—volume of fluid(CFD-DEM-VOF)model to simulate the complex behavior of particle-laden flows with free surfaces.The coupling between the fluid and particle phases is established through the implemented continuity,momentum,and alpha transport equation.The coupled particle forces such as drag,pressure gradient,dense virtual mass,viscous,and interface forces are also integrated,with drag and dense virtual mass forces being dependent on local porosity.The integrated conservative alpha transport equation ensures phase volume conservation during interactions between particles and water.Additionally,we have implemented a trilinear interpolation method designed to operate on unstructured hexahedral meshes.This method has been tested for its ability to properly resolve the coupling effects in the numerical simulations,particularly in cases with a relatively low cell-size ratio.The model is validated through three distinct test cases:single particle water entry,dam break with particles,and water entry of a group of particles case.The experimental setup is built to study the dynamics of the water entry of a group of particles,where three key flow features are analyzed:the evolution of average particle velocity,cavity shape,and particle dispersion cloud profiles in water.The tests involve four different scenarios,including two different water levels(16.1 and 20.1 cm)and two different particle densities(2650 and 4000 kg/m3).High-speed videometry and particle tracking velocimetry(using ImageJ/TrackMate)methods are employed for experimental data acquisition.It is demonstrated that numerical results are in excellent agreement with theoretical predictions and experimental data.The study highlights the significance of vortices in cavity shaping and particle dispersion.The validated CFD-DEM-VOF model constitutes a robust tool for simulating particle-laden flows,contributing valuable insights into the complex 展开更多
Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boun...Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boundary(IB)method developed in our previous work.For the moving structure modeled using the sharp interface IB method,a recursive box method is developed for efficiently classifying the background grid nodes.For the particles modeled using the diffuse interface IB method,a‘master-slave’approach is adopted.For the particle-particle interaction(PPI)and particle-structure interaction(PSI),a fast algorithm for classifying the active and inactive Lagrangian points,which discretize the particle surface,is developed for the‘dry’contact approach.The results show that the proposed recursive box method can reduce the classifying time from 52seconds to 0.3 seconds.Acceptable parallel efficiency is obtained for cases with different particle concentrations.Furthermore,the lubrication model is utilized when a particle approaches a wall,enabling an accurate simulation of the rebounding phenomena in the benchmark particle-wall collision problem.At last,the capability of the proposed computational framework is demonstrated by simulating particle-laden turbulent channel flows with rough walls.展开更多
The study of nano- and submicron Brownian particle-laden turbulent flow has wide industrial applicability and hence has received much attention. The purpose of the present paper is to provide and review some researche...The study of nano- and submicron Brownian particle-laden turbulent flow has wide industrial applicability and hence has received much attention. The purpose of the present paper is to provide and review some researches in this field. The topics are related to the universality, particularity, complexity and importance of nano- and submicron Brownian particle-laden turbulent flow, the models of particle general dynamical equation, the collision behavior of particles. Finally, several open research issues are identified.展开更多
The particle modulations to turbulence in round jets were experimentally studied by means of two-phase velocity measurements with Phase Doppler Anemometer （PDA）. Laden with very large particles, no significant atten...The particle modulations to turbulence in round jets were experimentally studied by means of two-phase velocity measurements with Phase Doppler Anemometer （PDA）. Laden with very large particles, no significant attenuations of turbulence intensities were measured in the farfields, due to small two-phase slip velocities and particle Reynolds number. The gas-phase turbulence is enhanced by particles in the near-fields, but it is significantly attenuated by the small particles in the far-fields. The smaller particles have a more profound effect on the attenuation of turbulence intensities. The enhancements or attenuations of turbulence intensities in the far-fields depends on the energy production, transport and dissipation mechanisms between the two phases, which are determined by the particle prop- erties and two-phase velocity slips. The non-dimensional parameter CTI is introduced to represent the change of turbulence intensity.展开更多
Axis-symmetric spheroids, such as rod-like and disk-like particles, have been found to orient preferentially in near-wall turbulence by both experiment and numerical simulation. In current work we examined the orienta...Axis-symmetric spheroids, such as rod-like and disk-like particles, have been found to orient preferentially in near-wall turbulence by both experiment and numerical simulation. In current work we examined the orientation of inertialess spheroids in a turbulent channel flow at medium friction Reynolds number Reτ=100 given based on the half of channel height. Both elongated prolate spheroid and flat oblate spheroid are considered and further compared with the reference case of spherical particle. The statistical results show that in near wall region the prolate spheroids tend to align in the streamwise direction while the oblate spheroids prefer to orient in the wallnormal direction, which are consistent with earlier observation in low Reynolds number (Reτ=180)wall turbulence. Around the channel center we found that the orientation of spheroids is not fully isotropic, even though the fluid vorticity are almost isotropic. The mechanism that gives rise to such particle orientations in wall-turbulence has been found to be related to fluid Lagrangian stretching and compression (Zhao and Andersson 2016). Therefore, we computed the left Cauchy-Green strain tensor along Lagrangian trajectories of tracer spheroids in current flow field and analyzed the fluid Lagrangian stretching and compression. The results indicated that, similar to the earlier observations, the directions of the Lagrangian stretching and compression in near-wall region are in the streamwise and wall-normal directions, respectively. Furthermore, cross over the channel the prolate spheroids aligned with the direction of Lagrangian stretching but oblate spheroids oriented with the direction of Lagrangian compression. The weak anisotropy of orientations of fluid Lagrangian stretching and compression observed at the channel center could be the reason for the aforementioned modest anisotropic orientation of spheroids in channel central region.展开更多
The conventional point-particle approach for treating the dispersed phase in a continuous flowfield is extended by taking into account the effect of finite particle size, using a Gaussian interpolation from Lagrangian...The conventional point-particle approach for treating the dispersed phase in a continuous flowfield is extended by taking into account the effect of finite particle size, using a Gaussian interpolation from Lagrangian points to the Eulerian field.The inter-phase exchange terms in the conservation equations are distributed over the volume encompassing the particle size, as opposed to the Dirac delta function generally used in the point-particle approach.The proposed approach is benchmarked against three different flow configurations in a numerical framework based on large eddy simulation（LES） turbulence closure.First, the flow over a circular cylinder is simulated for a Reynolds number of 3900 at 1 atm pressure.Results show good agreement with experimental data for the mean streamwise velocity and the vortex shedding frequency in the wake region.The calculated flowfield exhibits correct physics, which the conventional point-particle approach fails to capture.The second case deals with diesel jet injection in quiescent environment over a pressure range of 1.1–5.0 MPa.The calculated jet penetration depth closely matches measurements.It decreases with increasing chamber pressure, due to enhanced drag force in a denser fluid environment.Finally, water and acetone jet injection normal to air crossflow is studied at1 atm.The calculated jet penetration and Sauter mean diameter of liquid droplets compare very well with measurements.展开更多
This article investigates the near-field dynamics in a particle-laden round turbulent jet in a large-eddy simulation （LES）. A point-force two-way coupling model is adopted in the simulation to reveal the particle mo...This article investigates the near-field dynamics in a particle-laden round turbulent jet in a large-eddy simulation （LES）. A point-force two-way coupling model is adopted in the simulation to reveal the particle modulation of turbulence. The particles mainly excite the initial instability of the jet and bring about the earlier breakup of vortex rings in the near-field. The flow fluc- tuating intensity either in the axial or in the radial directions is hence increased by particles. The article also describes the mean velocity modulated by particles. The changing statistical velocity induced by particle modulation implies the effects of modulation of the local flow structures. This study is expected to be useful to the control of two-phase turbulent jets.展开更多
The alignment of elongated fibers and thin disks is known to be significantly influenced by the presence of fluid coherent structures in near-wall turbulence(Cui et al.2021).However,this earlier study is confined to t...The alignment of elongated fibers and thin disks is known to be significantly influenced by the presence of fluid coherent structures in near-wall turbulence(Cui et al.2021).However,this earlier study is confined to the spheroids with infinitely large or small aspect ratio,and the shape effect of finite aspect ratio on the alignment is not considered.The current study investigates the shape-dependent alignment of inertialess spheroids in structure-dominated regions of channel flow.With utilizing an ensemble-averaged approach for identifying the structure-dominated regions,we analyze the eigensystem of the linear term matrix in the Jeffery equation,which is governed by both particle shape and local fluid velocity gradients.In contrast to earlier conventional analysis based on local vorticity and strain rate,our findings demonstrate that the eigensystem of the Jeffery equation offers a convenient,effective,and universal framework for predicting the alignment behavior of inertialess spheroids in turbulent flows.By leveraging the eigensystem of the Jeffery equation,we uncover a diverse effect of fluid coherent structures on spheroid alignment with different particle shapes.Furthermore,we provide explanations for both shape-independent alignments observed in vortical-core regions and shape-dependent alignments around near-wall streamwise vortices.展开更多
A robust and efficient solver coupling computational fluid dynamics(CFD)with discrete element method(DEM)is developed to simulate particle-laden flows in various physical settings.An interpolation algorithm suitable f...A robust and efficient solver coupling computational fluid dynamics(CFD)with discrete element method(DEM)is developed to simulate particle-laden flows in various physical settings.An interpolation algorithm suitable for unstructured meshes is proposed to translate between mesh-based Eulerian fields and particle-based Lagrangian quantities.The interpolation scheme reduces the mesh-dependence of the averaging and interpolation procedures.In addition,the fluid-particle interaction terms are treated semi-implicitly in this algorithm to improve stability and to maintain accuracy.Finally,it is demonstrated that sub-stepping is desirable for fluid-particle systems with small Stokes numbers.A momentum-conserving sub-stepping technique is introduced into the fluid-particle coupling procedure,so that problems with a wide range of time scales can be solved without resorting to excessively small time steps in the CFD solver.Several numerical examples are presented to demonstrate the capabilities of the solver and the merits of the algorithm.展开更多
The absence of sub-grid scale(SGS) motions leads to severe errors in particle pair dynamics, which represents a great challenge to the large eddy simulation of particle-laden turbulent flow. In order to address this i...The absence of sub-grid scale(SGS) motions leads to severe errors in particle pair dynamics, which represents a great challenge to the large eddy simulation of particle-laden turbulent flow. In order to address this issue,data from direct numerical simulation(DNS) of homogenous isotropic turbulence coupled with Lagrangian particle tracking are used as a benchmark to evaluate the corresponding results of filtered DNS(FDNS). It is found that the filtering process in FDNS will lead to a non-monotonic variation of the particle collision statistics, including radial distribution function, radial relative velocity, and the collision kernel. The peak of radial distribution function shifts to the large-inertia region due to the lack of SGS motions, and the analysis of the local flowstructure characteristic variable at particle position indicates that the most effective interaction scale between particles and fluid eddies is increased in FDNS. Moreover,this scale shifting has an obvious effect on the odd-order moments of the probability density function of radial relative velocity, i.e. the skewness, which exhibits a strong correlation to the variance of radial distribution function in FDNS.As a whole, the radial distribution function, together with radial relative velocity, can compensate the SGS effects for the collision kernel in FDNS when the Stokes number based on the Kolmogorov time scale is greater than 3.0. However,it still leaves considerable errors for St< 3.0.展开更多
This work focuses on numerically studying hydrodynamic interaction between a passive particle and a self-propelled particle,termed a squirmer,by using a two-dimensional lattice Boltzmann method(LBM).It is found that t...This work focuses on numerically studying hydrodynamic interaction between a passive particle and a self-propelled particle,termed a squirmer,by using a two-dimensional lattice Boltzmann method(LBM).It is found that the squirmer can capture a passive particle and propel it simultaneously,provided the passive particle is situated within the squirmer's wake.Our research shows that the critical capture distance,which determines whether the particle is captured,primarily depends on the intensity of the squirmer's dipolarity.The stronger dipolarity of squirmer results in an increased critical capture distance.Conversely,the Reynolds number is found to have minimal influence on this interaction.Interestingly,the passive particle,when driven by the squirmer's wake,contributes to a reduction in the squirmer's drag.This results in a mutual acceleration for both particles.Our findings can provide valuable perspectives for formulating the principles of reducing the drag of micro-swimmers and help to achieve the goal of using micro-swimmers to transport goods without physical tethers.展开更多
This study investigates turbulent particle-laden channel flows using direct numerical simulations employing the Eulerian-Lagrangian method.A two-way coupling approach is adopted to explore the mutual interaction betwe...This study investigates turbulent particle-laden channel flows using direct numerical simulations employing the Eulerian-Lagrangian method.A two-way coupling approach is adopted to explore the mutual interaction between particles and fluid flow.The considered cases include flow with particle Stokes number varying from St=2 up to St=100 while maintaining a constant Reynolds number of Reτ=180 across all cases.A novel vortex identification method,Liutex(Rortex),is employed to assess its efficacy in capturing near-wall turbulent coherent structures and their interactions with particles.The Liutex method provides valuable information on vortex strength and vectors at each location,enabling a detailed examination of the complex interaction between fluid and particulate phases.As widely acknowledged,the interplay between clockwise and counterclockwise vortices in the near-wall region gives rise to low-speed streaks along the wall.These low-speed streaks serve as preferential zones for particle concentration,depending upon the particle Stokes number.It is shown that the Liutex method can capture these vortices and identify the location of low-speed streaks.Additionally,it is observed that the particle Stokes number(size)significantly affects both the strength of these vortices and the streaky structure exhibited by particles.Furthermore,a quantitative analysis of particle behavior in the near-wall region and the formation of elongated particle lines was carried out.This involved examining the average fluid streamwise velocity fluctuations at particle locations,average particle concentration,and the normal velocity of particles for each set of particle Stokes numbers.The investigation reveals the intricate interplay between particles and near-wall structures and the significant influence of particles Stokes number.This study contributes to a deeper understanding of turbulent particle-laden channel flow dynamics.展开更多
The instantaneous degradation of erosion surface of ethylene propylene diene monomer(EPDM)insulation subjected to the particle-laden flow in two operating conditions was measured by using a real-time X-ray radiography...The instantaneous degradation of erosion surface of ethylene propylene diene monomer(EPDM)insulation subjected to the particle-laden flow in two operating conditions was measured by using a real-time X-ray radiography system.The images of its erosion state and dynamic ablation rate were obtained.And the charring-layer was analyzed by using SEM and energy spectrum.The experimental results indicate that the erosion rate of EPDM insulation layer impacted by low speed and low concentration particle flow is relatively small in the 1st second since the motor starting,but increases rapidly in 1 to 2.5 s,while the erosion rate of EPDM insulation layer impacted by high speed and high concentration particle flow increases rapidly in the 1st second;the ablation rate at the section eroded intensively by particle flow increases at first,then decreases,and goes to stabilization after 4.5 s;the higher speed and concentration particle flow are,the deeper particles get into charring layer,which lead to more thermal increment and thinner charring layer.展开更多
Non-spherical particles exist widely in natural and industrial fluid systems and the motions of nonspherical particles are significantly different from that of spherical particles.In this paper,a simplified model of n...Non-spherical particles exist widely in natural and industrial fluid systems and the motions of nonspherical particles are significantly different from that of spherical particles.In this paper,a simplified model of non-spherical particles considering particle drag correction,lift,and rotation was established.Based on the Eulerian-Lagrangian simulation,the dispersion characteristics of spherical and nonspherical particles with different Stokes numbers in a high-speed turbulent jet were analyzed and compared considering the effect of particle rotation.The results show that,the differences in particle dispersion and radial velocity fluctuation between non-spherical particles and spherical particles in the jet are significant,especially when Stokes number is large.Moreover,the effects of different type of forces on the dispersion of non-spherical particles and spherical particles were compared in detail,which revealed that the change of the Magnus force caused by the increase in the angular velocity of non-spherical particles plays a dominant role in the differences of particle dispersion.展开更多
The micro- and macro-time scales in two-phase turbulent channel flows are investigated using the direct nu- merical simulation and the Lagrangian particle trajectory methods for the fluid- and the particle-phases, res...The micro- and macro-time scales in two-phase turbulent channel flows are investigated using the direct nu- merical simulation and the Lagrangian particle trajectory methods for the fluid- and the particle-phases, respectively. Lagrangian and Eulerian time scales of both phases are cal- culated using velocity correlation functions. Due to flow anisotropy, micro-time scales are not the same with the theo- retical estimations in large Reynolds number （isotropic） tur- bulence. Lagrangian macro-time scales of particle-phase and of fluid-phase seen by particles are both dependent on particle Stokes number. The fluid-phase Lagrangian inte- gral time scales increase with distance from the wall, longer than those time scales seen by particles. The Eulerian inte- gral macro-time scales increase in near-wall regions but de- crease in out-layer regions. The moving Eulerian time scales are also investigated and compared with Lagrangian integral time scales, and in good agreement with previous measure- ments and numerical predictions. For the fluid particles the micro Eulerian time scales are longer than the Lagrangian ones in the near wall regions, while away from the walls the micro Lagrangian time scales are longer. The Lagrangian integral time scales are longer than the Eulerian ones. The results are useful for further understanding two-phase flow physics and especially for constructing accurate prediction models of inertial particle dispersion.展开更多
基金This work was partially supported by the State Key Project of Fundamental Research,Ministry of Science and Technology,China(Grant Nos.G1999022207,2002CB211600)the National Natural Science Foundation of China(Grant Nos.50276021,50576027)Program for New Century Excellent Talents in University,Ministry of Education,China(Grant No.NCET-04-0708).
文摘The fluid temperature statistics along particle trajectories is crucial to under-stand the mechanisms of turbulent non-isothermal or reactive fluid-particle flow,especially for the Lagrangian model of non-isothermal particle-laden turbulent flow.In the present study,direct numerical simulations were utilized to generate temperature field statistics in particle-laden incompressible stationary homogeneous isotropic turbulent flows,which is focused on the effect of particle response time on the Lagrangian statistics of the particle and the fluid temperature seen by particles.It shows that,for the particles withτp/τk<1,the ratio of the fluid intensity seen by particle to fluid temperature intensities deceased asτp/τk increased;while for larger particles(τp/τk>1),the trend is inversed.For small parti-cles(τp/τk<5),the Lagrangian autocorrelation coefficient of the particle temperature R_(p)^(T)decreases as the particle inertia(τp/τk)increases.The trend is reversed for larger particles.The autocorrelation of fluid temperature along the particle path,R_(pf)^(T),decreased as the particle inertia increased.And as the particle inertia increased,the autocorrelation coeffi-cient of the fluid temperature seen by particle decreased more rapidly than that of the particle temperature.The mean temperature gradient contributes to the correlation be-tween the particles velocity component and temperature fluctuations in the direction of the gradient.For the particles withτp/τk<1,the magnitude of the correlation coefficient in-creases as the particle inertia increases,while this value is independent of the particle time constant for larger particles.
基金The project supported by the National Natural Science Foundation of China （50276021）, and Program for New Century Excellent Talents in University, Ministry of Education （NCET-04-0708） The English text was polished by Yunming Chen.
文摘Particle-laden flows in a horizontal channel were investigated by means of a two-phase particle image velocimetry （PIV） technique. Experiments were performed at a Reynolds number of 6826 and the flow is seeded with polythene beads of two sizes, 60μm and 110μm. One was slightly smaller than and the other was larger than the Kolmogorov length scale. The particle loadings were relatively low, with mass loading ratio ranging from 5 ×10^-4 to 4 × 10^-2 and volume fractions from 6×10×-7 to 4.8×10^-5, respectively. The results show that the presence of particles can dramatically modify the turbulence even under the lowest mass loading ratio of 5 × 10^-4. The mean flow is attenuated and de- creased with increasing particle size and mass loading. The turbulence intensities are enhanced in all the cases concerned. With the increase of the mass loading, the intensities vary in a complicated manner in the case of small particles, indicating complicated particle-turbulence interactions; whereas they increase monotonously in the case of large particles. The particle velocities and concentrations are also given. The particles lag behind the fluid in the center region but lead in the wall region, and this trend is more prominent for the large particles. The streamwise particle fluctuations are larger than the gas fluctuations for both sizes of particles, however their varying trend with the mass loadings is not so clear. The wallnormal fluctuations increase with increasing mass loadings. They are smaller in the 60μm particle case but larger in the 110μm particle case than those of the gas phase. It seems that the small particles follow the fluid motion to certain extent while the larger particles are more likely dominated by their own inertia. Finally, remarkable non-uniform distributions of particle concentration are observed, especially for the large particles. The inertia of particles is proved to be very important for the turbulence modification and particles behaviors and thus should be considered in horizonta
文摘Fully developed vertical turbulent channel flow with particle transport wasinvestigated by use of Large Eddy Simulation (LES) approach coupled with dynamic the Sub-Grid Scale(SGS) model. It was assumed that the motion of each particle is followed in a Lagrangian frame ofreference driven by the forces exerted by fluid motion and gravity under the condition of one-waycoupling. The goal of this study is to investigate the effectiveness of the LES technique forpredicting particle transport in turbulent flow and the behavior of particle-laden turbulent channelflow for three kinds of particles at different Stokes numbers. To depict the behavior ofparticle-laden turbulent channel flow, statistical quantities including particle fluctuation andfluid-particle velocity correlation, and visualization of the particle number density field wereanalyzed.
基金supported by the National Natural Science Foundation of China (Grants 11272198, 11572183)
文摘The dynamic and thermal performance of particle-laden turbulent flow is investigated via direction numerical simulation combined with the Lagrangian point-particle tracking under the condition of two-way coupling, with a focus on the contributions of particle feedback effect to momentum and heat transfer of turbulence. We take into account the effects of particles on flow drag and Nusselt number and explore the possibility of drag reduction in conjunction with heat transfer enhancement in particle-laden turbulent flows. The effects of particles on momentum and heat transfer are analyzed, and the possibility of drag reduction in conjunction with heat transfer enhancement for the prototypical case of particle-laden turbulent channel flows is addressed. We present results of turbulence modification and heat transfer in turbulent particle-laden channel flow, which shows the heat transfer reduction when large inertial particles with low specific heat capacity are added to the flow. However, we also found an enhancement of the heat transfer and a small reduction of the flow drag when particles with high specific heat capacity are involved. The present results show that particles, which are active agents, interact not only with the velocity field, but also the temperature field and can cause a dissimilarity in momentum and heat transport. This demonstrates that the possibility to increase heat transfer and suppress friction drag can be achieved with addition of particles with different thermal properties.
基金Sepro Mineral Systems through the Mitacs Accelerate project No:IT12396.We would also like to express our appreciation to Compute Canada and UBC ARC for supporting PIME Lab by granting us access to their high-performance computing platforms.The simulations presented in this work were conducted using the Compute Canada Beluga and UBC ARC Sockeye clusters,both of which contributed equally to the computational resources utilized.
文摘This paper presents the development and validation of a fully coupled computational fluid dynamics—discrete element method—volume of fluid(CFD-DEM-VOF)model to simulate the complex behavior of particle-laden flows with free surfaces.The coupling between the fluid and particle phases is established through the implemented continuity,momentum,and alpha transport equation.The coupled particle forces such as drag,pressure gradient,dense virtual mass,viscous,and interface forces are also integrated,with drag and dense virtual mass forces being dependent on local porosity.The integrated conservative alpha transport equation ensures phase volume conservation during interactions between particles and water.Additionally,we have implemented a trilinear interpolation method designed to operate on unstructured hexahedral meshes.This method has been tested for its ability to properly resolve the coupling effects in the numerical simulations,particularly in cases with a relatively low cell-size ratio.The model is validated through three distinct test cases:single particle water entry,dam break with particles,and water entry of a group of particles case.The experimental setup is built to study the dynamics of the water entry of a group of particles,where three key flow features are analyzed:the evolution of average particle velocity,cavity shape,and particle dispersion cloud profiles in water.The tests involve four different scenarios,including two different water levels(16.1 and 20.1 cm)and two different particle densities(2650 and 4000 kg/m3).High-speed videometry and particle tracking velocimetry(using ImageJ/TrackMate)methods are employed for experimental data acquisition.It is demonstrated that numerical results are in excellent agreement with theoretical predictions and experimental data.The study highlights the significance of vortices in cavity shaping and particle dispersion.The validated CFD-DEM-VOF model constitutes a robust tool for simulating particle-laden flows,contributing valuable insights into the complex
基金Project supported by the National Natural Science Foundation of China(Nos.12202456 and12172360)the Basic Science Center Program for“Multiscale Problems in Nonlinear Mechanics”of the National Natural Science Foundation of China(No.11988102)the China Postdoctoral Science Foundation(No.2021M693241)。
文摘Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boundary(IB)method developed in our previous work.For the moving structure modeled using the sharp interface IB method,a recursive box method is developed for efficiently classifying the background grid nodes.For the particles modeled using the diffuse interface IB method,a‘master-slave’approach is adopted.For the particle-particle interaction(PPI)and particle-structure interaction(PSI),a fast algorithm for classifying the active and inactive Lagrangian points,which discretize the particle surface,is developed for the‘dry’contact approach.The results show that the proposed recursive box method can reduce the classifying time from 52seconds to 0.3 seconds.Acceptable parallel efficiency is obtained for cases with different particle concentrations.Furthermore,the lubrication model is utilized when a particle approaches a wall,enabling an accurate simulation of the rebounding phenomena in the benchmark particle-wall collision problem.At last,the capability of the proposed computational framework is demonstrated by simulating particle-laden turbulent channel flows with rough walls.
基金Project supported by the Major Program of the National Natural Science Foundation of China(Grant No.11132008)
文摘The study of nano- and submicron Brownian particle-laden turbulent flow has wide industrial applicability and hence has received much attention. The purpose of the present paper is to provide and review some researches in this field. The topics are related to the universality, particularity, complexity and importance of nano- and submicron Brownian particle-laden turbulent flow, the models of particle general dynamical equation, the collision behavior of particles. Finally, several open research issues are identified.
基金supported by the National Natural Science Foundation of China (50876053 and 50706021).
文摘The particle modulations to turbulence in round jets were experimentally studied by means of two-phase velocity measurements with Phase Doppler Anemometer （PDA）. Laden with very large particles, no significant attenuations of turbulence intensities were measured in the farfields, due to small two-phase slip velocities and particle Reynolds number. The gas-phase turbulence is enhanced by particles in the near-fields, but it is significantly attenuated by the small particles in the far-fields. The smaller particles have a more profound effect on the attenuation of turbulence intensities. The enhancements or attenuations of turbulence intensities in the far-fields depends on the energy production, transport and dissipation mechanisms between the two phases, which are determined by the particle prop- erties and two-phase velocity slips. The non-dimensional parameter CTI is introduced to represent the change of turbulence intensity.
基金the financial support from the National Natural Science Foundation of China(91752205,11702158 and 11490551)granted by the Programme for Supercomputing(NN2649K)
文摘Axis-symmetric spheroids, such as rod-like and disk-like particles, have been found to orient preferentially in near-wall turbulence by both experiment and numerical simulation. In current work we examined the orientation of inertialess spheroids in a turbulent channel flow at medium friction Reynolds number Reτ=100 given based on the half of channel height. Both elongated prolate spheroid and flat oblate spheroid are considered and further compared with the reference case of spherical particle. The statistical results show that in near wall region the prolate spheroids tend to align in the streamwise direction while the oblate spheroids prefer to orient in the wallnormal direction, which are consistent with earlier observation in low Reynolds number (Reτ=180)wall turbulence. Around the channel center we found that the orientation of spheroids is not fully isotropic, even though the fluid vorticity are almost isotropic. The mechanism that gives rise to such particle orientations in wall-turbulence has been found to be related to fluid Lagrangian stretching and compression (Zhao and Andersson 2016). Therefore, we computed the left Cauchy-Green strain tensor along Lagrangian trajectories of tracer spheroids in current flow field and analyzed the fluid Lagrangian stretching and compression. The results indicated that, similar to the earlier observations, the directions of the Lagrangian stretching and compression in near-wall region are in the streamwise and wall-normal directions, respectively. Furthermore, cross over the channel the prolate spheroids aligned with the direction of Lagrangian stretching but oblate spheroids oriented with the direction of Lagrangian compression. The weak anisotropy of orientations of fluid Lagrangian stretching and compression observed at the channel center could be the reason for the aforementioned modest anisotropic orientation of spheroids in channel central region.
基金sponsored by the William R.T.Oakes Endowment of the Georgia Institute of Technology
文摘The conventional point-particle approach for treating the dispersed phase in a continuous flowfield is extended by taking into account the effect of finite particle size, using a Gaussian interpolation from Lagrangian points to the Eulerian field.The inter-phase exchange terms in the conservation equations are distributed over the volume encompassing the particle size, as opposed to the Dirac delta function generally used in the point-particle approach.The proposed approach is benchmarked against three different flow configurations in a numerical framework based on large eddy simulation（LES） turbulence closure.First, the flow over a circular cylinder is simulated for a Reynolds number of 3900 at 1 atm pressure.Results show good agreement with experimental data for the mean streamwise velocity and the vortex shedding frequency in the wake region.The calculated flowfield exhibits correct physics, which the conventional point-particle approach fails to capture.The second case deals with diesel jet injection in quiescent environment over a pressure range of 1.1–5.0 MPa.The calculated jet penetration depth closely matches measurements.It decreases with increasing chamber pressure, due to enhanced drag force in a denser fluid environment.Finally, water and acetone jet injection normal to air crossflow is studied at1 atm.The calculated jet penetration and Sauter mean diameter of liquid droplets compare very well with measurements.
基金National Natural Science Foundation of China (50706021)Ph.D.Programs Foundation of Ministry of Education of China (20070003018)TNList Cross-discipline Foundation
文摘This article investigates the near-field dynamics in a particle-laden round turbulent jet in a large-eddy simulation （LES）. A point-force two-way coupling model is adopted in the simulation to reveal the particle modulation of turbulence. The particles mainly excite the initial instability of the jet and bring about the earlier breakup of vortex rings in the near-field. The flow fluc- tuating intensity either in the axial or in the radial directions is hence increased by particles. The article also describes the mean velocity modulated by particles. The changing statistical velocity induced by particle modulation implies the effects of modulation of the local flow structures. This study is expected to be useful to the control of two-phase turbulent jets.
基金supported by the Natural Science Foundation of China(Grant Nos.92252104,92252204,12388101,and 12302285)the China Postdoctoral Science Foundation(Grant No.2022M721849).
文摘The alignment of elongated fibers and thin disks is known to be significantly influenced by the presence of fluid coherent structures in near-wall turbulence(Cui et al.2021).However,this earlier study is confined to the spheroids with infinitely large or small aspect ratio,and the shape effect of finite aspect ratio on the alignment is not considered.The current study investigates the shape-dependent alignment of inertialess spheroids in structure-dominated regions of channel flow.With utilizing an ensemble-averaged approach for identifying the structure-dominated regions,we analyze the eigensystem of the linear term matrix in the Jeffery equation,which is governed by both particle shape and local fluid velocity gradients.In contrast to earlier conventional analysis based on local vorticity and strain rate,our findings demonstrate that the eigensystem of the Jeffery equation offers a convenient,effective,and universal framework for predicting the alignment behavior of inertialess spheroids in turbulent flows.By leveraging the eigensystem of the Jeffery equation,we uncover a diverse effect of fluid coherent structures on spheroid alignment with different particle shapes.Furthermore,we provide explanations for both shape-independent alignments observed in vortical-core regions and shape-dependent alignments around near-wall streamwise vortices.
文摘A robust and efficient solver coupling computational fluid dynamics(CFD)with discrete element method(DEM)is developed to simulate particle-laden flows in various physical settings.An interpolation algorithm suitable for unstructured meshes is proposed to translate between mesh-based Eulerian fields and particle-based Lagrangian quantities.The interpolation scheme reduces the mesh-dependence of the averaging and interpolation procedures.In addition,the fluid-particle interaction terms are treated semi-implicitly in this algorithm to improve stability and to maintain accuracy.Finally,it is demonstrated that sub-stepping is desirable for fluid-particle systems with small Stokes numbers.A momentum-conserving sub-stepping technique is introduced into the fluid-particle coupling procedure,so that problems with a wide range of time scales can be solved without resorting to excessively small time steps in the CFD solver.Several numerical examples are presented to demonstrate the capabilities of the solver and the merits of the algorithm.
基金supported by the National Natural Science Foundation of China (Grants 51390494, 51306065, and 51276076)the Foundation of State Key Laboratory of Coal Combustion (Grant FSKLCCB1702)
文摘The absence of sub-grid scale(SGS) motions leads to severe errors in particle pair dynamics, which represents a great challenge to the large eddy simulation of particle-laden turbulent flow. In order to address this issue,data from direct numerical simulation(DNS) of homogenous isotropic turbulence coupled with Lagrangian particle tracking are used as a benchmark to evaluate the corresponding results of filtered DNS(FDNS). It is found that the filtering process in FDNS will lead to a non-monotonic variation of the particle collision statistics, including radial distribution function, radial relative velocity, and the collision kernel. The peak of radial distribution function shifts to the large-inertia region due to the lack of SGS motions, and the analysis of the local flowstructure characteristic variable at particle position indicates that the most effective interaction scale between particles and fluid eddies is increased in FDNS. Moreover,this scale shifting has an obvious effect on the odd-order moments of the probability density function of radial relative velocity, i.e. the skewness, which exhibits a strong correlation to the variance of radial distribution function in FDNS.As a whole, the radial distribution function, together with radial relative velocity, can compensate the SGS effects for the collision kernel in FDNS when the Stokes number based on the Kolmogorov time scale is greater than 3.0. However,it still leaves considerable errors for St< 3.0.
基金Project supported by the National Natural Science Foundation of China (Grant Nos.12132015 and 11972336)。
文摘This work focuses on numerically studying hydrodynamic interaction between a passive particle and a self-propelled particle,termed a squirmer,by using a two-dimensional lattice Boltzmann method(LBM).It is found that the squirmer can capture a passive particle and propel it simultaneously,provided the passive particle is situated within the squirmer's wake.Our research shows that the critical capture distance,which determines whether the particle is captured,primarily depends on the intensity of the squirmer's dipolarity.The stronger dipolarity of squirmer results in an increased critical capture distance.Conversely,the Reynolds number is found to have minimal influence on this interaction.Interestingly,the passive particle,when driven by the squirmer's wake,contributes to a reduction in the squirmer's drag.This results in a mutual acceleration for both particles.Our findings can provide valuable perspectives for formulating the principles of reducing the drag of micro-swimmers and help to achieve the goal of using micro-swimmers to transport goods without physical tethers.
基金supported by the National Natural Science Foundation of China (Grant No.92052202)the Fundamental Research Funds for the Central Universities (Grant No.lzujbky-2021-ey19).
文摘This study investigates turbulent particle-laden channel flows using direct numerical simulations employing the Eulerian-Lagrangian method.A two-way coupling approach is adopted to explore the mutual interaction between particles and fluid flow.The considered cases include flow with particle Stokes number varying from St=2 up to St=100 while maintaining a constant Reynolds number of Reτ=180 across all cases.A novel vortex identification method,Liutex(Rortex),is employed to assess its efficacy in capturing near-wall turbulent coherent structures and their interactions with particles.The Liutex method provides valuable information on vortex strength and vectors at each location,enabling a detailed examination of the complex interaction between fluid and particulate phases.As widely acknowledged,the interplay between clockwise and counterclockwise vortices in the near-wall region gives rise to low-speed streaks along the wall.These low-speed streaks serve as preferential zones for particle concentration,depending upon the particle Stokes number.It is shown that the Liutex method can capture these vortices and identify the location of low-speed streaks.Additionally,it is observed that the particle Stokes number(size)significantly affects both the strength of these vortices and the streaky structure exhibited by particles.Furthermore,a quantitative analysis of particle behavior in the near-wall region and the formation of elongated particle lines was carried out.This involved examining the average fluid streamwise velocity fluctuations at particle locations,average particle concentration,and the normal velocity of particles for each set of particle Stokes numbers.The investigation reveals the intricate interplay between particles and near-wall structures and the significant influence of particles Stokes number.This study contributes to a deeper understanding of turbulent particle-laden channel flow dynamics.
基金Sponsored by the National Nature Science Foundation of China(50976095)
文摘The instantaneous degradation of erosion surface of ethylene propylene diene monomer(EPDM)insulation subjected to the particle-laden flow in two operating conditions was measured by using a real-time X-ray radiography system.The images of its erosion state and dynamic ablation rate were obtained.And the charring-layer was analyzed by using SEM and energy spectrum.The experimental results indicate that the erosion rate of EPDM insulation layer impacted by low speed and low concentration particle flow is relatively small in the 1st second since the motor starting,but increases rapidly in 1 to 2.5 s,while the erosion rate of EPDM insulation layer impacted by high speed and high concentration particle flow increases rapidly in the 1st second;the ablation rate at the section eroded intensively by particle flow increases at first,then decreases,and goes to stabilization after 4.5 s;the higher speed and concentration particle flow are,the deeper particles get into charring layer,which lead to more thermal increment and thinner charring layer.
基金supported by National Science and Technology Major Project of China(2019-I-0022-0021)Seed fund of Shanxi Research Institute for Clean Energy,Tsinghua Universitythe National Natural Science Foundation of China(51761125011)。
文摘Non-spherical particles exist widely in natural and industrial fluid systems and the motions of nonspherical particles are significantly different from that of spherical particles.In this paper,a simplified model of non-spherical particles considering particle drag correction,lift,and rotation was established.Based on the Eulerian-Lagrangian simulation,the dispersion characteristics of spherical and nonspherical particles with different Stokes numbers in a high-speed turbulent jet were analyzed and compared considering the effect of particle rotation.The results show that,the differences in particle dispersion and radial velocity fluctuation between non-spherical particles and spherical particles in the jet are significant,especially when Stokes number is large.Moreover,the effects of different type of forces on the dispersion of non-spherical particles and spherical particles were compared in detail,which revealed that the change of the Magnus force caused by the increase in the angular velocity of non-spherical particles plays a dominant role in the differences of particle dispersion.
基金supported by the National Natural Science Foundation of China (11132005 and 50706021)
文摘The micro- and macro-time scales in two-phase turbulent channel flows are investigated using the direct nu- merical simulation and the Lagrangian particle trajectory methods for the fluid- and the particle-phases, respectively. Lagrangian and Eulerian time scales of both phases are cal- culated using velocity correlation functions. Due to flow anisotropy, micro-time scales are not the same with the theo- retical estimations in large Reynolds number （isotropic） tur- bulence. Lagrangian macro-time scales of particle-phase and of fluid-phase seen by particles are both dependent on particle Stokes number. The fluid-phase Lagrangian inte- gral time scales increase with distance from the wall, longer than those time scales seen by particles. The Eulerian inte- gral macro-time scales increase in near-wall regions but de- crease in out-layer regions. The moving Eulerian time scales are also investigated and compared with Lagrangian integral time scales, and in good agreement with previous measure- ments and numerical predictions. For the fluid particles the micro Eulerian time scales are longer than the Lagrangian ones in the near wall regions, while away from the walls the micro Lagrangian time scales are longer. The Lagrangian integral time scales are longer than the Eulerian ones. The results are useful for further understanding two-phase flow physics and especially for constructing accurate prediction models of inertial particle dispersion.
