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展开更多
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.展开更多
低氧稀释(moderate and intense low-oxygen dilution,MILD)燃烧具有传热均匀、NO_x污染物排放低的特点。新一代MILD煤粉燃烧技术主要通过高速射流引起强烈的湍流混合来实现。其中,煤颗粒在高速射流中的扩散行为非常关键。目前对球形颗...低氧稀释(moderate and intense low-oxygen dilution,MILD)燃烧具有传热均匀、NO_x污染物排放低的特点。新一代MILD煤粉燃烧技术主要通过高速射流引起强烈的湍流混合来实现。其中,煤颗粒在高速射流中的扩散行为非常关键。目前对球形颗粒在气固射流中扩散行为的研究已非常深入,然而化石燃料属于典型的非球形颗粒,其在射流中的扩散行为与球形颗粒具有一定的差异,该类非球形颗粒在高速射流下的扩散特性值得进一步研究。为此,该文采用玻璃珠、玻璃渣和煤粉等颗粒开展了宽Re范围下的高速两相圆射流实验,通过激光Doppler相位分析技术(phase-Doppler anemometry,PDA)获取并分析了颗粒的质量浓度、速度及湍动能分布随球形度、粒径以及射流速度的变化规律。结果表明:非球形颗粒在射流中的质量浓度、速度、湍动能分布与粒径较小的球形颗粒具有一定相似性,但其扩散行为不能仅通过Stokes数进行定量表征,除曳力之外,升力对非球形颗粒扩散也具有一定影响;与球形颗粒相比,非球形颗粒的扩散更为显著,其主要原因是其径向湍动能显著增强所致;射流速度的增加促进了颗粒的剪切层集聚和径向扩散,对非球形颗粒的促进作用更强。展开更多
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 展开更多
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.展开更多
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 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 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.展开更多
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 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.展开更多
Particle-resolved direct numerical flow solvers predominantly use a projection method to decouple the non-linear mass and momentum conservation equations.The computing performance of such solvers often decays beyond O...Particle-resolved direct numerical flow solvers predominantly use a projection method to decouple the non-linear mass and momentum conservation equations.The computing performance of such solvers often decays beyond O(1000)cores due to the cost of solving at least one large three-dimensional pressure Poisson problem per time step.The parallelization may perform moderately well only or even poorly sometimes despite using an efficient algebraic multigrid preconditioner[38].We present an accurate and scalable solver using a direction splitting algorithm[12]to transform all three-dimensional parabolic/elliptic problems(and in particular the elliptic pressure Poisson problem)into a sequence of three one-dimensional parabolic sub-problems,thus improving its scalability up to multiple thousands of cores.We employ this algorithm to solve mass and momentum conservation equations in flows laden with fixed non-spherical rigid bodies.We consider the presence of rigid bodies on the(uniform or non-uniform)fixed Cartesian fluid grid by modifying the diffusion and divergence stencils on the impacted grid node near the rigid body boundary.Compared to[12],we use a higher-order interpolation scheme for the velocity field to maintain a secondorder stress estimation on the particle boundary,resulting in more accurate dimensionless coefficients such as drag C_(d)and lift C_(l).We also correct the interpolation scheme due to the presence of any nearby particle to maintain an acceptable accuracy,making the solver robust even when particles are densely packed in a sub-region of the computational domain.We present classical validation tests involving a single or multiple(up to O(1000))rigid bodies and assess the robustness,accuracy and computing speed of the solver.We further show that the Direction Splitting solver is∼5 times faster on 5120 cores than our solver[38]based on a classical projection method[5].展开更多
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.展开更多
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.展开更多
基金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
基金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.
基金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
基金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.
文摘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.
基金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.
基金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).
基金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.
基金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.
基金support of the University of British Columbia via its Four Year Doctoral Fellowship programThe authors greatly appreciate the financial support of the Natural Sciences and Engineering Research Council of Canada(NSERC)via Anthony Wachs’s Discovery Grant RGPIN-2016-06572+1 种基金This research was enabled by the support provided by Compute Canada(http://www.computecanada.ca)through Anthony Wachs’s 2020,2021,and 2022 Resources for Research Groups allocation qpf-764-abThis research was also supported in part through computational resources and services provided by Advanced Research Computing at the University of British Columbia.
文摘Particle-resolved direct numerical flow solvers predominantly use a projection method to decouple the non-linear mass and momentum conservation equations.The computing performance of such solvers often decays beyond O(1000)cores due to the cost of solving at least one large three-dimensional pressure Poisson problem per time step.The parallelization may perform moderately well only or even poorly sometimes despite using an efficient algebraic multigrid preconditioner[38].We present an accurate and scalable solver using a direction splitting algorithm[12]to transform all three-dimensional parabolic/elliptic problems(and in particular the elliptic pressure Poisson problem)into a sequence of three one-dimensional parabolic sub-problems,thus improving its scalability up to multiple thousands of cores.We employ this algorithm to solve mass and momentum conservation equations in flows laden with fixed non-spherical rigid bodies.We consider the presence of rigid bodies on the(uniform or non-uniform)fixed Cartesian fluid grid by modifying the diffusion and divergence stencils on the impacted grid node near the rigid body boundary.Compared to[12],we use a higher-order interpolation scheme for the velocity field to maintain a secondorder stress estimation on the particle boundary,resulting in more accurate dimensionless coefficients such as drag C_(d)and lift C_(l).We also correct the interpolation scheme due to the presence of any nearby particle to maintain an acceptable accuracy,making the solver robust even when particles are densely packed in a sub-region of the computational domain.We present classical validation tests involving a single or multiple(up to O(1000))rigid bodies and assess the robustness,accuracy and computing speed of the solver.We further show that the Direction Splitting solver is∼5 times faster on 5120 cores than our solver[38]based on a classical projection method[5].
基金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.
基金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.
