基于Lagrange方法封闭的两相湍流场方程模型-计算物理.PDF

第３３卷 第３期 计　 　 算　 　 物　 　 理 Ｖｏｌ．３３，Ｎｏ．３ ２０１６年５月 　 ＣＨＩＮＥＳＥ ＪＯＵＲＮＡＬ ＯＦ ＣＯＭＰＵＴＡＴＩＯＮＡＬ ＰＨＹＳＩＣＳ　 Ｍａｙ，２０１６ 文章编号：１００１⁃２４６Ｘ（２０１６）０３⁃０３０５⁃０６ 基于Ｌａｇｒａｎｇｅ方法封闭的两相湍流场方程模型 １ １，２，∗ ２ 王　 路 ，　 徐江荣 ，　 刘保银 （１． 南京理工大学能源与动力工程学院，南京２１００９４；２． 杭州电子科技大学能源研究所，杭州３１００１８） 摘　 要：两相湍流场方程模型采用基于Ｅｕｌｅｒ方法的一阶矩方程，而二阶矩方程由Ｌａｇｒａｎｇｅ方法得到，新模型的封 闭不需要附加其它假设．首先基于概率密度函数给出颗粒运动的连续方程和动量方程，其次由基于平均Ｌａｎｇｅｖｉｎ 方程的Ｌａｇｒａｎｇｅ模型推导得到颗粒二阶矩方程，最终获得封闭的二阶矩模型．将新模型用于气固两相壁面射流的 数值模拟，结果表明新模型合理有效． 关键词：场方程模型；Ｅｕｌｅｒ方法；Ｌａｇｒａｎｇｅ方法；数值模拟 中图分类号：Ｏ３５ 文献标志码：Ａ ０　 引言 雷诺方程在连续相湍流的数值模拟中被广泛应用，湍流数值模拟的零方程模型、单方程模型、双方程模 型、雷诺应力模型、代数应力模型等层出不穷，这些模型都是以平均Ｎａｖｉｅｒ⁃Ｓｔｏｋｅｓ方程为基础，而平均Ｎａｖｉｅｒ⁃ Ｓｔｏｋｅｓ方程本身是不封闭的，需要附加各种近似假设和经验公式加以封闭，由于其本身具有无法克服的缺 ［１］ 陷，他们的使用范围具有很大的局限性，比较具有代表性的标准ｋ⁃ε模型 是以Ｂｏｕｓｓｉｎｅｓｑ假设为基础，它 只对高雷诺数的湍流有效，而且对模拟旋流和绕流也有缺陷． 与连续相雷诺方程类似，离散相雷诺方程也在数值模拟中被广泛采用，其核心思想是在Ｅｕｌｅｒ体系中构 ［２］ 造方程使统计矩方程封闭，Ｒｅｅｋｓ基于概率密度函数（ＰＤＦ）框架得到双流体模型方程 ，利用 Ｃｈａｐｍａｎ⁃ ［３］ Ｅｎｓｋｏｇ近似法得到三阶矩方程使模型封闭．Ｓｉｍｏｎｉｎ 也采用概率密度函数方法推导颗粒的统计矩方程， ［４］ ［５］ Ｍｉｎｉｅｒ 和Ｐｅｒａｎｏ 以及Ｍｉｎｉｅｒ和Ｐｏｚｏｒｓｋｉ 由ＰＤＦ输运方程推导出类似雷诺展开得到的双流体模型，包括 ［６］ 两相雷诺应力和湍流动能输运方程，用代数方程对不封闭的方程组加以封闭．ｖａｎ Ｓｌｏｏｔｅｎ，Ｊａｙｅｓｈ 仿照Ｒｏｄｉ 方程对雷诺应力方程进行处理，得到了颗粒应力代数方程，与颗粒相的连续方程、动量方程组成封闭的方程 ［７］ ［８］ ［９］ 组．国内的周力行等人先后提出的ｋ⁃ε⁃ｋ 模型 、统一二阶矩模型 、代数应力模型 ，都采用经验公式或ｐ 者近似假设，与连续相雷诺方程一样，其适用范围自然受到种种限制． 两相湍流流动的另外一种处理思路是离散相在Ｌａｇｒａｎｇｅ体系下采用颗粒的轨道方法，连续相依然采用 ［１０］ ［１１］ 雷诺方程数值模拟，这是一种非常重要并被广泛采用的计算方法，Ｘｕ 和Ｐｏｐｅ ，Ｌｉ 和Ｍｏｄｅｓｔ 以及柳朝 ［１