可压缩颗粒圆孔射流并行直接数值模拟算法研究.pdfVIP

高性能计算应用 59 可压缩颗粒圆孔射流并行直接数值模拟算法研究  李德波 广东电网公司电力科学研究院 广州510060 ldbyx@126.com  樊建人 罗坤 岑可法 浙江大学 能源清洁利用国家重点实验室 杭州310027 摘要： 本文对空间模式发展的气固两相圆孔射流中颗粒与流体双相耦合作用进行了并行环境下 的直接数值模拟算法研究。气相流场采用可压缩的N-S方程直接求解。计算颗粒场时，采用 Lagrangian方法跟踪实际的颗粒运动。利用并行求解算法，实现了颗粒穿越边界面的模拟。为 了模拟颗粒对流体的作用，考虑了颗粒和流体的双相耦合。在本文的计算条件下，颗粒的直径 远小于网格的间距，平均的Kolmogorov尺度和网格的间距在一个量级，保证了DNS的要求。气 相和颗粒相的应力与实验的对比研究表明，本文的颗粒并行程序是可信的。 关键词：气固两相，直接数值模拟，并行算法 引言 为了实现大规模的并行计算，整个流场和颗粒 气固两相流动在自然界是一种最广泛的流动形 场的计算采用并行求解。气相场采用区域分解算法 [5] 式。电厂煤粉颗粒的燃烧，喷雾燃烧等等。在实际 进行并行计算。本文中提到的区域分解算法，与文 [5] 的颗粒运动中，颗粒的扩散成为影响系统效率和稳 献中的算法有较大的差异 。空间离散采用显式8阶 [6] 定性的重要因素。气固两相流动的理论研究也是多 格式 ，时间项采用5步4阶Runge-kutta离散。采用 相流研究的一个热点。DNS和LES方法在多相流研 显式格式的最大优势是能克服并行求解的困难。目 究中，得到了广泛的应用。在颗粒圆孔射流研究方 前，在追求高精度离散Navier-Stokes方程文献中，有 [1] 许多研究者采用高阶精致格式离散。精致格式离散 面，DNS方法研究还比较少。Yuu et all首先进行了 圆孔射流直接数值模拟和实验研究。实验中流体的 后的方程是三对角方程组，需要采用追赶法求解。 Reynolds数为1700，数值模拟和实验取得了基本一致 通过我们反复的数值实验表明，紧致格式的并行效 [2] 率要远低于显式格式，尤其是大规模的数值模拟计 的结果。随后，Yuu与其合作者进行了圆孔射流的 LES数值模拟研究和实验研究。模拟的Reynolds数为 算中，这种差异还会更加明显。本文的区域分解只 10000。通过LDV实际测量，得到的平均速度，湍流 是并行求解过程中，人为的将求解的区域分解成若 强度和Reynolds stress与实验是吻合比较好。 干个小的子区域，区域的个数取决于模拟的CPU的 在颗粒两相研究中，根据已有的文献来看，用 数目。但是在CPU与CPU之间的边界上的数据，完 可压缩Navier-Stokes方程求解流场，还比较少。本文 全是靠传值的方法实现，没有采用任何近似求解的 [3,4] 算法，保证我们得到的并行求解的结果理论上与单 采用最近发展起来的无反射边界条件 ，求解了可 压缩的N-S方程。克服了边界上反射波对内场的影 CPU求解的结果是完全一致的。每个CPU的网格数受 响，成功实现了射流在边界上运动的过程。可压缩