节点文献

均匀湍流内有限惯性颗粒碰撞率的直接模拟

Direct Numerical Simulation the Collision Rate of Finite Inertia Particles Suspended in Homogeneous Isotropic Turbulence

【作者】 陈胤密

【导师】 郑楚光; 柳朝晖;

【作者基本信息】 华中科技大学 , 热能工程, 2004, 硕士

【摘要】 燃烧源可吸入颗粒物对于人体健康和自然环境有严重的危害作用,了解其在燃烧过程中生成和长大过程的机理是在燃烧过程中控制颗粒物生成的前提,其核心问题之一是颗粒-颗粒的碰撞和凝并过程。而由于湍流对细微颗粒具有复杂的输运作用和局部聚集效应,细颗粒的碰撞率将不同于常规的几何碰撞率。研究具有有限惯性、并有局部富集条件下颗粒的碰撞率是目前两相流体力学和气溶胶动力学的研究热点之一。本文在系统总结有关颗粒平均碰撞率模型的研究进展和相关数学模型、算法的基础上,基于后溯算法,构建了能够模拟三维空间内具105~106个颗粒碰撞的直接模拟程序,其中颗粒动量方程采用二阶预估-校正格式,颗粒位移采用二阶Adams积分,气相湍流场采用拟谱方法求解,颗粒位置上流体的速度由气相场用四阶拉格朗日插值得到,通过采用相邻网格、网格排除和颗粒排序等方法对程序进行优化,程序计算量较常规方法可低2个数量级,满足了在微机上模拟巨量颗粒碰撞需要。使用该程序,首先对均匀剪切流内St<<1的颗粒(轻颗粒流)和St>1且速度呈高斯分布颗粒(重颗粒流)两种工况进行数值验证,其结果和文献中的解析解一致,证明了算法和程序的可靠性,并讨论了计算参数的选取原则。进一步的,本文对各向均匀同性湍流内有限惯性颗粒的聚集现象及其对碰撞率的影响进行了模拟计算,采用径向分布函数、归一化的条件颗粒浓度等参数定量表征了不同St数颗粒的局部聚集程度;对St=1的颗粒,受局部聚集的影响,其颗粒碰撞率可为零惯性颗粒的20~30倍。

【Abstract】 Inhaled particulate matter originating from combustion does serious harm to human health and environment. Mechanism of formation and growth of inhaled particulate matter in combustion is the precondition of its control during combustion process, where the controlling process is inter-particle collision and coagulation. Because turbulence imposes complex transport and accumulation effects on fine particle, the collision rate of fine particles there may different from those predicted by the custom model. Research on collision rate of finite inertial particle with preferential concentration is one of hotspots in gas-particle flow mechanism and aerosol dynamics.After a brief summary of the research developments in particle average collision rate model, correlated mathematic models and numerical simulation algorithms, a retroactive algorithm-based program was built for direct numerical simulation (DNS) the collisions between particles in 3-dimenional space, with particle number up to 105~106. The second-order forecast-correction method was applied in solution of particle’s momentum equation, where the fluid velocity at the particle location is determined using a fourth-order Lagrangian interpolation. The particle position was advanced by a second-order Adams integral method. The fluid velocity field was generated by a pseudo spectral method. By optimizing algorithms such as cell index, exclusive search and particle sequence et al., the program can fit the requirements of simulating collisions between huge numbers of particles on personal computer, whose computation load is two orders of magnitude less than that of custom algorithm. Verifications of the code were done by two cases: St<<1 particle in uniform shear flow and St>1 particle whose velocity distribution was Gaussian. The results were consistent with analytical solutions in literatures,which proved reliability of retroactive algorithm and the program. The principle of computational parameter selection was also discussed. Moreover, DNS study of finite inertial particle with preferential concentration in isotropic turbulence and its influence on collision rate were performed. The accumulation level of various St number particle was quantified by parameters such as radial distribution function, and conditional expectation of normalized particle number density et al.; It showed that collision rate of St =1 particle was 20~30 times larger than that of zero-inertial particle because of local preferential concentration.

  • 【分类号】TK16
  • 【被引频次】3
  • 【下载频次】284
节点文献中: 

本文链接的文献网络图示:

本文的引文网络