惯性颗粒所见被动标量统计特性的直接数值模拟研讨.pdf

中国工程热物理学会 多相流 学术会议论文 编号：094258 惯性颗粒所见被动标量统计特性的直接 数值模拟研究 刘亚明，柳朝晖，郑楚光 （华中科技大学，武汉，湖北，430074） （Freebirdlym2002@163.com, 0278409） 摘要：本文采用直接数值模拟方法，对惯性颗粒所见均匀各向同性湍流中具有平均标量梯度的被动标 量场统计特性进行了研究，探讨了颗粒惯性的影响。模拟结果表明：颗粒所见标量的自相关特性随颗 粒惯性的增加而单调减少；颗粒所见标量脉动能随颗粒stokes 数的增加先减少再增大，在st~0.1 附近 达到最小值，而颗粒所见标量耗散率随颗粒stokes 数的变化行为与颗粒所见标量的脉动能恰恰相反； 在颗粒路径上，应变率和涡量与标量耗散率之间的关联系数也进行了统计，在St~1.0 附近分别达到最 大和最小，其值为0.3 和0.25 左右，；此外，反映颗粒路径上局部混合状态的颗粒所见力和热时间尺度 比值 rp ，比Euler 方法的统计结果要偏低15~20%，我们推荐在相关模型中采用1.77 作为模型常数 关键词：被动标量，惯性颗粒，直接数值模拟，Lagrangian 方法 0. 前言： 在自然界和工程流动中，湍流是流体运动中最普遍的现象，它一个非常有用的特性 就是能迅速加快流场中的标量混合，这些标量可以是温度场，也可以是反应物的浓度场。 被动标量，只单向考虑湍流对标量的输运作用，而忽略标量对湍流的反向影响。对被动 标量的研究一直以来是湍流研究领域中一个热点，一方面它有助于加深对湍流自身的理 解，另一方面，对于改进和发展非等温反应多相流等工程应用模型也具有重要意义[1] 。 近些年来，随着计算机资源的迅速发展，数值模拟研究的方法已经同理论分析和实 验测量一样成为标量场研究的重要手段，其中不采用任何假设求解流体NS 方程和标量 对流扩散方程的直接数值模拟方法已经成为分析简单流动中标量场各种细观机理的关 键途径[2-8] 。在非等温/反应的气固两相流中，影响标量场混合的因素更为复杂，一方面 颗粒与标量场（温度场）之间存在着直接的相间交换，另一方面颗粒可以改变流场的动 力学结构来间接影响标量场的混合，迄今为止，只有少数该方面的直接数值模拟研究报 道。Sato 等人[9]通过对具有平均标量梯度的均匀各向同性湍流中标量场的直接数值模拟 研究，发现在平均标量梯度方向上速度和温度脉动具有很强的关联；Jaberi[10]在大尺度 上进行强制得到稳态的速度场和标量场，其数值模拟研究表明颗粒的温度脉动强度是颗 粒的stokes 数，标量Pr 数和流场雷诺数的递减函数；Jaberi 和Mashayek[11]随后对Jaberi[10] 的工作进行了扩展，在自由衰减的均匀流场和分别考虑颗粒与流场之间单向和双向耦合 的条件下进行了类似的讨论。不过，这些文献都是基于标量场和流场全场Euler 平均的 结果，而从改进和发展非等温/反应气固两相流模型的角度，尤其是对于物理上更为合理 的衔接拉氏动力学模型与相应的欧拉模型的 PDF 方法，颗粒路径上所见的流体微团标 量特性更为关键。Yeung[12] 曾经采用lagrangian 的方法，统计了零惯性的示踪粒子所见 标量的相关信息，发现所见标量的脉动强度并非是马尔科夫过程，并且具有比流体自身 更长的关联时间尺度；He 等[13]也曾部分讨论过惯性颗粒所见标量的相关统计量与颗粒 参数之间的关系，考虑了颗粒的响应时间对颗粒所见标量自相关和颗粒温度脉动强度的 影响。本文是对He 等人[13]工作的延伸，对其它重要的统计量，尤其是所见标量耗散率 的有关特性，进行了讨论，例如颗粒所见标量耗散率的强度，颗粒路径上标量耗散与流 场动力学结构等物理机理进行了讨论。 在本文中，主要讨论颗粒惯性的影响，为避免其它因素的干扰把颗粒视为惰性颗粒，即 湍流的变动和与标量场之间的交换不予考虑。 1. 数值模拟方法 1.1 控制方程