欧拉-拉格朗日下颗粒相大规模数值计算关键问题研究.pdfVIP

高性能计算应用 43 欧拉-拉格朗日下颗粒相大规模数值计算 关键问题研究  李德波 徐齐胜 曾庭华 廖永进沈跃良 广东电网公司电力科学研究院 广州510060 ldbyx@126.com 摘要： 本文对气固两相流数值模拟中，离散相数值计算方法在商业软件(Ansys Fluent 14.0)和直接 数值模拟中差异性进行了对比研究。主要对离散相数值算法具体实现方法进行了深入分析。本 文分析对比结果对于今后进行气固两相流动工程问题的商业软件模拟，提供重要的理论上的参 考，以及不同气固两相工程问题数值模拟中商业软件模拟需要主要的关键问题，提供有价值的 参考。 关键词：气固两相，商业软件，Ansys Fluent 14.0，湍流扩散，直接数值模拟 前言 问题的研究，实际工程问题的数值模拟还需要较长 气固两相流动是自然界最具有代表性的流动形 的时间；第二种方法是雷诺平均数值模拟(RANS)。 式。在火力发电厂煤粉燃烧过程，化工领域中气固 这种方法是将瞬态的Navier-Stokes方程进行雷诺平 相分选等都涉及到气固两相流动。目前国际上研究 均，忽略到湍流脉动的信息，得到时均场的控制方 气固两相流动的方法主要分为两类：一种是将流体 程组。这种方法是目前计算代价最小的数值模拟技 作为连续介质，在欧拉坐标系内考虑流体，而颗粒 术，但是只能得到时均湍流场的信息，无法反应湍 相作为离散相，在拉格朗日坐标系下来描述；第二 流脉动的细节信息；第三种数值模拟方法是大涡数 种方法将流体作为连续介质，颗粒群当成拟连续介 值模拟方法（LES）。这种方法的特点是将湍流不同 质或者拟流体来考虑，流体和颗粒相在空间共存和 尺度的涡结构进行空间过滤，大尺度的涡结构采用 互相渗透，统一在欧拉坐标系加以描述。随着商业 直接数值模拟方法求解，过滤尺度以下的湍流场采 用模型求解。这种方法是上述的DNS和RANS的一种 软件的不断发展(Ansys Fluent 14.0, Ansys CFX, Start- [1-3] 折中方案，避免了湍流所有尺度上的精确的求解， CD)等 ，针对气固两相流动的数值模拟的商业软件 (Ansys Fluent 14.0, Ansys CFX, Start-CD)[1-3]也广泛应 仅仅在过滤尺度以下采用模型进行求解。大涡数值 用到工业界实际问题的研究中，研究结果为实际工 模拟方法是目前学术界和工业界发展比较快的一种 业过程提供了重要的参考信息，提高了工业过程设 方法，已经有研究者采用LES方法进行实际工程问 [4-19] 计的效率。 题的数值模拟 ，在目前的计算机硬件技术下， 进行湍流气固两相流动的数值模拟，目前主要 LES方法可以解决部分中等网格规模的工程问题。可 [4-19] 以预计在不久的将来，随着大规模并行计算机硬件 有三种尺度上的数值模拟 。一种是直接数值模拟 方法(Direct numerical simulation)。这种方法是精确地 技术和发展，LES方法将是工业界进行实际问题数值 求解描述流体运动规律的Navier-Stokes方程，得到空 模拟的首选方法，因为LES数值模拟的精度要明显 间和时间上湍流所有尺度的信息，是目前为止数值 高于RANS，同时计算量的代价要小于DNS。正是在 模拟精度最高的一种模拟方法，同时也是计算量最 上述三种不同湍流尺度上的数值模拟，也导致了目 高的数值