一维气液两相漂移模型的ausmv算法研究.pdf

　 ， ３５ １２ Ａｐｐｌｉｅｄ Ｍａｔｈｅｍａｔｉｃｓ ａｎｄ Ｍｅｃｈａｎｉｃｓ　 　 ２０１４年１２月１５ 日出版 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 Ｖｏｌ．３５，Ｎｏ．１２，Ｄｅｃ．１５，２０１４ 文章编号：１０００⁃０８８７（２０１４）１２⁃１３７３⁃１０ ⓒ 应用数学和力学编委会，ＩＳＳＮ １０００⁃０８８７ 一维气液两相漂移模型的ＡＵＳＭＶ算法研究∗ 徐朝阳，　 孟英峰，　 魏　 纳，　 李　 皋，　 万里平 （油气藏地质及开发工程国家重点实验室（西南石油大学），成都 ６１０５００） 摘要：　 将ＡＵＳＭＶ（ａｄｖｅｃｔｉｏｎ ｕｐｓｔｒｅａｍ ｓｐｌｉｔｔｉｎｇ ｍｅｔｈｏｄ Ｖ）格式从计算气体动力学问题扩展至一维 等温瞬态气液两相管流．阐述了采用ＡＵＳＭＶ格式构建气液两相漂移模型数值通量的方法及边界 单元的处理方法．采用Ｒｕｎｇｅ⁃Ｋｕｔｔａ方法与经典的保单调ＭＵＳＣＬ（ｍｏｎｏｔｏｎｅｕｐｓｔｒｅａｍ⁃ｃｅｎｔｒｅｄｓｃｈｅｍｅｓ ｆｏｒ ｃｏｎｓｅｒｖａｔｉｏｎ ｌａｗｓ）方法结合Ｖａｎ Ｌｅｅｒ 限制器，构建具有二阶时间和空间精度的数值计算方法．计 算经典Ｚｕｂｅｒ⁃Ｆｉｎｄｌａｙ激波管问题和复杂漂移关系变质量流动问题并与可靠的参考结果进行了对 比．分析表明：ＡＵＳＭＶ格式应用于气液两相流动漂移模型时计算效率高、精度高、耗散效应和色散 效应小，低流速条件下能够精确地描述间断． 关　 键　 词：　 一维等温气液两相流；　 漂移模型；　 ＡＵＳＭＶ；　 数值计算 ＋ 中图分类号：　 Ｏ３５９ ．１；Ｏ２４１．８２　 　 　 文献标志码：　 Ａ ｄｏｉ：１０．３８７９／ ｊ．ｉｓｓｎ．１０００⁃０８８７．２０１４．１２．００９ 引　 　 言 气液两相管流问题在石油、天然气、化工等领域是重要的研究课题，高效、高精度的数值计 算方法是保障理论研究与工程应用的重要基础．现在广泛应用的气液两相流动模型为漂移模 型（ｄｒｉｆｔ⁃ｆｌｕｘ ｍｏｄｅｌ，ＤＦＭ）和分相流模型（ｔｗｏ⁃ｆｌｕｉｄ ｍｏｄｅｌ，ＴＦＭ）．分相流模型ＴＦＭ 由Ｄｅｌｈａｙｅ于 １９６８年提出，模型考虑各相间作用力，计算精度高，但是形式复杂，计算量较大．漂移模型由 Ｚｕｂｅｒ等提出，相比分相流模型具有形式简单和计算量小的优点，且计算精度满足工程需 ［１⁃４］ 求 ．漂移模型相比单相流体的Ｎａｖｉｅｒ⁃Ｓｔｏｋｅｓ方程组更为复杂，Ｂｅｎｚｏｎｉ⁃Ｇａｖａｇｅ，Ｔｈéｒｏｎ和Ｂａｒ⁃ ｎｅａ分析了漂移模型控制方程组的数学性质，他们的研究表明：在现实工程中的各相流速区域 ［５⁃７］ 内（０～１００ ｍ／ ｓ），漂移模型属于双曲型方程组 ．因为，漂移模型的流通量无法完全以守恒参 数表示，其Ｊａｃｏｂｉ矩阵没有解析特征值．所以，研究不依赖（近似）Ｒｉｅｍａｎｎ求解器或Ｊａｃｏｂｉ矩 阵构造数值通量的计算方法成为重要的研究方向． ＡＵＳＭ（ａｄｖｅｃｔｉｏｎ ｕｐｓｔｒｅａｍ ｓｐｌｉｔｔｉｎｇ ｍｅｔｈｏｄ）类格式是成功不依赖（近似）Ｒｉｅｍａｎｎ求解器或 Ｊａｃｏｂｉ矩阵计算Ｎ⁃Ｓ方程组的著名方法．１９９３年，Ｌｉｏｕ 和Ｓｔｅｆｆｅｎ结合通量差分分裂格式（ｆｌｕｘ ｄｉｆｆｅｒｅｎｃｅ ｓｐｌｉｔｔｉｎｇ，ＦＤＳ）和矢通量分裂格式（ｆｌｕｘ ｖｅｃｔｏｒ ｓｐｌｉｔｔｉｎｇ，ＦＶＳ）的优点，提出了著名的 ［８］ ＡＵＳＭ格式，其基本思路是将无黏通量分为对流通量项和压力通量项分别处理 ．在ＡＵＳＭ格 ∗ 收稿日期：　 ２０１４⁃０６⁃２３；修订日期：　 ２０１４⁃０９⁃２２ 基金项目：　 国家科技重大专项（２０１１ＺＸ０５０２１⁃００３）；国家自然科学基金 作者简介：　 徐朝阳（１９８５—），男，四川仪陇人，博士生（Ｅ⁃ｍａｉｌ：０４０１１ｘｚｙ＠ｓｉｎａ．ｃｏｍ）； 孟英峰（１９５４—），男，