二相流计算的一种差分算法ξ.pdfVIP

( ) 　应用数学和力学 ,第 25 卷 第 5 期 2004 年 5 月 应用数学和力学编委会编 　　　　　　　　　　　 　 　　Applied Mathematics and Mechanics 重 庆 出 版 社 出 版 ( ) 文章编号 :1000-0887 2004 05-0488-09 二相流计算的一种差分算法 1 2 1 1 李 　强 , 　封建湖 , 　蔡体敏 , 　胡春波 ( 1. 西北工业大学 164 信箱 ,西安 710072 ; 2 . 长安大学 理学院 ,西安 710064) (李家春推荐) 摘要 : 　考虑两相流的力学行为 ,忽略相间的耗散作用 ,建立了 Euler 型的基本控制方程 　·状态方 程采用刚性状态方程 　·基于 Abgrall 提出的准则 ,在流动区域内,对可压两相流提出了一个精度较 高的 Euler 型数值方法 ,数值格式是 Godunov 型格式 ,对守恒型和非守恒型方程采用 HLLC 型和Lax- Friedrichs 型近似 Riemann 解算器 ,引入了速度驰豫和压强驰豫过程来代替两相间的相互作用 　·在 一维情形下给出数值算例 ,并且和 Saurel 的算例进行了比较 ,结果表明该算法不但精确而且稳定 , 且在间断处没有数值振荡 　· 关 　键 　词 : 　可压二相流 ; 　迟豫过程 ; 　近似 Riemann 解算器 中图分类号 : 　O359 　　　文献标识码 : 　A 引　　言 二相流问题的研究在现代工农业生产中具有很重要的地位 　·迄今 ,研究二相流有两种方 法 ,一种将流体相视为连续介质 ,将颗粒相视为离散体系 ,在拉格朗日坐标系下探讨颗粒运动 时物理量的变化 ,这种方法称为颗粒轨迹模型或拉格朗日法 　·另一种方法是把流体和颗粒两 相均视为同时充满流场而且相互作用的连续介质进行研究 ,该方法被称为双流体模型或欧拉 法 　·本文重点研究双流体模型　·对于双流体模型的数值模拟问题 , 目前由于人们对其流动机 理还不清楚 ,对其波结构还不甚了解 ,所以数值研究还处于初步阶段 ,现有的方法主要是将求 解守恒型问题的方法移植过来 　·在该方面已有一些研究 ,Toumi [ 1] 在 G. Dal Maso 、P. Le Floch[2 ] 等工作的基础之上 ,试图采用线性化方法求解 ,但如何选取合适的线性化路径这一问题尚未解 决 ;L . Sainaulieu[3 ] 在线性化方法的基础上采用有限体积法求解了二相流问题 ,但计算结果的误 差比较大 ,对间断的识别能力差 　·R . Saurel [4 ] 在 Abgrall [5 ] 工作基础上 ,利用 HLL 格式对守恒型 方程进行了离散 ,其数值模拟结果在间断附近有非物理振荡　·本文以D . A . Drew[6 ] 等提出的方 程为基本控制方程 , 以 HLLC 格式为基础 ,对方程进行数值求解 　· 1 　模 型 介 绍 本文中使用 D . A . Drew[6 ] 、Saurel 和 Gallouet [7 ] 中提出的可压二相流方程 : 收稿 日期 : 　2002-05- 17 ; 修订日期 : 　2003- 12-03 ( ) 基金项 目: 　航空科学基金资助项 目 03c53021 ( ) ( 作者简介 : 　李强