FLUENT中应用DPM模型时双R分布的详细说明.docx

F L U E N T 中 应 用 D P M 模 型 时 双 R 分 布 的 详 细 说 明 使用动网格的模型在应用 DPM 模型进行计算时， InjectionType 不能使用 surface 。 关于 rosin-rammler 分布 举例说明，有一组颗粒服从这样一种粒径分布，见下表： Diameter MassFraction Range( μm) inRange 0-70 0.05 70-100 0.10 100-120 0.35 120-150 0.30 150-180 0.15 180-200 0.05 定义一个变量 Yd，其定义为：比指定粒径 d 大的颗粒的质量分数。那么上面所说的颗粒的粒径分布所对 应的 Yd 就是： MassFractionwith Diameter, d( μm) DiameterGreaterthan d ,Yd 70 0.95 100 0.85 120 0.50 150 0.20 180 0.05 200 (0.00) Rosin-Rammler 分布函数假定粒径 d 和 Yd 只见存在这样一种指数关系： ????????????? Yd=(e -( d/ dm ) ) n ????????????? （1 ） 其中 d[size=10.5pt]m 为平均粒径（ MeanDiameter ）； n 为传播系数（ SpreadParameter ）。为了获得 上述两种数值，需要找到 d 和 Yd 的关系。 MassFractionwith Diameter, d ( μm) DiameterGreaterthan d ,Yd 70 0.95 100 0.85 120 0.50 150 0.20 180 0.05 200 (0.00) Yd=e -1 ≈ 0.368 所对应的 d 值即为 d[size=10.5pt]m ，由于上表中没有 0.368 ，所以需要根据已有数值进行 拟合，得到曲线如下： 根据上图找到 Yd=0.368 所对应的 d 值，在这里 d[size=10.5pt]m [size=10.5pt]= d≈ 131 μm。 得到 d[size=10.5pt]m 后，根据式（ 1 ）可以得到式（ 2 ） ?????????? n=ln(-ln Yd )/ln( d/ dm)???? （ 2 ） 用式（ 2 ）来计算 n 值。共有 6 种粒径，但是只能得到 5 个 n 值，因为最后一个 Yd=0 。将 5 个 n 值进 行平均，最终便得到 n 的大小。在这个例子中 n=4.52 。 关于 n 的计算，可以使用 Excel 来做。个人感觉还是很好使的。 掌握了 rosin-rammler 分布规则后，便可在使用 DPM 模型时根据实际情况设定不同的粒径了。 但需要注意的是 rosin-rammler 分布是一种连续分布 然而，在 fluent 计算的时候，需要设定 [size=10.5pt]Numberofdiameters ，并且在计算结果中可以 看到，粒径的分布并不是连续的。如果粒径分成 10 级，那么只能看见 10 种粒径。而且，均从同一个颗 粒源发射出。 关于 InjectionType 说明 ??? 常用 InjectionType 是 group 和 surface ，采用 group 形式时需要设定第一个颗粒源和最后一个颗 粒源的位置，而且这种形式可以设定颗粒源的个数；但是 surface 就不能设定颗粒源的个数，颗粒源的 位置默认分布在发射面的每一个网格上。所以在网格数目很大时，第二种形式会有较多的颗粒源，这要 消耗巨大的计算时间。而第一种形式由于可控治颗粒源数量，故可以设定相对较少的颗粒源。但麻烦的 是在三维模型中需要设置不止一个颗粒源。