气固两相流动大规模并行直接数值模拟算法研究进展与展望.PDFVIP

高性能计算应用 51 气固两相流动大规模并行直接数值模拟算法 研究进展与展望  李德波 广东电网公司电力科学研究院 广州510060 ldbyx@126.com  樊建人 罗坤 岑可法 浙江大学 能源清洁利用国家重点实验室 杭州310027 摘要： 对欧拉-拉格朗日下气固两相流动并行直接数值模拟算法的最新研究成果进行了系统的总 结。主要涉及到气固两相流大规模并行直接数值模拟空间导数离散格式的选取，流体相和颗粒 相并行区域划分的策略，颗粒相并行数据传递的逻辑控制方法，颗粒相碰撞并行的实现方法， 三维计算区域二个方向分解并行计算的若干关键问题，以及并行加速比和并行效率的计算结 果。指明了解决欧拉-拉格朗日下气固两相流动并行直接数值模拟算法若干关键问题，以及今 后需要解决的若干关键问题。 关键词：欧拉-拉格朗日，气固两相流，直接数值模拟，颗粒碰撞算法，大规模并行计算 引言 计算方法的研究，深入研究两相流动的难题，比如 气固两相流动的数值模拟主要有两大类方法， 湍流调制，优先聚集，颗粒与流体的耦合之间的质 一种是流体相在欧拉坐标系下求解，颗粒相在拉格 量、动量耦合效应等等，都需要依赖于准确的试验 朗日下进行颗粒轨道计算，这种方法称为欧拉-拉格 数据，这里包括直接的试验研究和数值模拟的数据 [4-5] 朗日方法；另外一种方法不仅流体相在欧拉坐标系 。因此，开展两相流动并行计算方法的研究显得 下进行计算，将颗粒看作一种拟流体，采用与流体 非常必要。可以预计在不久的将来，随着并行计算 相类似的控制方程，在欧拉坐标系下求解，这种方 技术的不断发展，必然会产生更加精细的两相流动 [1-5] 的计算方法，得到的颗粒与流体作用的内在的信息 法称为欧拉-欧拉方法 。本文主要介绍在欧拉-拉 格朗日框架下，进行气固两相流体大规模并行直接 更加细致，也必将大量的依赖于大规模的并行计算 [10] 数值模拟的最新研究进展。 机，因此这样大规模颗粒并行算法显得尤为重要 。 [10-12,16-25] 随着大规模并行计算机硬件和软件技术的发 李德波等 系统研究了气固两相流动大 展，为两相流动的直接数值模拟研究提供了新的平 规模并行直接数值模拟算法。研究了颗粒两相流并 台。以前两相流动的研究，往往受到计算资源的限 行直接数值模拟空间离散选取的方法，流体相和颗 制，只能考虑颗粒与流体之间的单相耦合，过去研 粒相并行区域划分的策略，颗粒相并行数据传递的 究者气固两相流数值模拟中也采用计算颗粒的概 逻辑控制方法，颗粒相碰撞并行的实现方法，三维 念，其主要目的是为了节省计算量，减少实际跟踪 计算区域二个方向区域分解并行计算的若干关键问 的颗粒的数目。同时颗粒与颗粒之间的碰撞，颗粒 题，以及并行加速比和并行效率的计算结果。颗粒 [6-9] 两相流并行数值模拟的结果与实验进行了严格的对 与流体之间的双向耦合作用往往不考虑 。在并 行计算技术快速发展的今天，开展两相流动并