基于输运方程的空化两相流算法设计.pdfVIP

第九届全国水动力学学术会议暨第二十二届全国水动力学研讨会文集 基于输运方程的空化两相流算法设计木 张凌新赵伟国邵雪明 浙江大学工程力学系流体工程研究所，杭州，310027 Email：圣b曼翌g!i翌g茎i望壁圣。j婪：皇堕坠：生旦：堡曼垒皇兰里壁圣j坚：金鱼坐：垒翌 摘要： 空化现象常见于各类液体流动中。本研究从数值计算的角度出发，对空化流计 算方法进行了研究。首先，以基于输运方程的空化两相流模型为基础，给出了空化流动的数 值计算方程。分析密度变化的来源，推导了带有相变源项的物质守恒方程；采用基于压力的 对该方程中的压力源项进行了隐式化处理。 关键词： 空化流；空化模型；基于压力的算法；数值计算 1 引言 水在常温下为液体，但当局部压力低于饱和蒸汽压时，水会汽化形成空化现象。空化现 象常见于各类水动力系统中，如螺旋桨、水翼、水下潜体和水利机械等。空化现象一旦发生， 流场的计算将变得比较困难，因为其中既有不可压的水。又有压缩性极强的水汽混合物。 求解流体运动有两类方法：一类是基于压力的，一类是基于密度的。基于密度的方法常 用于可压缩性流动中，然而对于水流场的计算不太有利，冈为此时水的密度变化比较微弱。 而基于压力的方法则能很好地处理不可压缩流动，著名的SIMPLE算法就是这一类方法。然而 区别于一般的不可压缩流动，空化流的密度场也是变化的，而且空化区的混合流体具有极强 的压缩性，所以基于压力的算法需要进一步的改进来适应这种特殊流动的计算。Senocak和 Shyy㈠’从可压缩性的角度考虑了密度与压力的耦合关系，类似于可压缩流体建立了空化流动 的压力修正关系。张凌新和符松㈠’也从这一角度将全速算法应用于空化流的算法设计。然而 上述方法在构建密度压力关系时还存在一定的任意性，需要进一步研究。 本研究采剧一种新的思路：即通过分析空化流动的相变来源，得到不包含密度随体导数 的物质守恒方程，在此基础上建立压力修上E方程。这一思路避免了密度与压力的直接耦合， 相变的物质守恒关系由空化模型确定，在理论上更为严谨。基于输运方程的两相流模型由于 1 其包含了更多的物理机制而被广泛的运用于空化流建模及计算，本研究将对Singhal等人13 提出的质量分数模型和Kunz等人…提出的体积分数模型进行分析，给出基于这两种输运方程 的空化流算法设计。 ‘感谢国家自然科学基金青年基金项目的支持 -317． 第九届全国水动力学学术会议暨第二十二届全国水动力学研讨会文集 2控制方程 假设水汽 两相拥有统一的流动变量，空化流的控制方程如-卜．所示 连续方程： 害+v．(pⅣ)_o’ 动量方程： (2) V 7。 掣再(胛)可一卯+v．f， 岔 湍流模型方程(七一s模式)： 掣帅妒丑…V·陋珈]， (3) 掣毗妒e，妻托：户知№珈]。 (4) 出于简洁，省略了各变量的湍流平均符号，这里p和∥均表示两相流场的混合物质属性。对 于空化流动来讲，还需要添加一个空化模型方程，我们将在后面给出描述物质分数的空化模 型方程。 为阐述方便，以二维坐标系为例。对中心点为尸的单元控制体积分，得到流动控制方程 的通用守恒形式： (5)