基于能量转换的空化流动计算方法.pdfVIP

中固T程热物理学会 叶轮机械气动热力学 学术会议论文 编号：062026 基于能量转换的空化流动计算方法 黄典贵1蔡瑞忠2 1(上海大学力学研究所上海200072) Tel：021 dghuang@staff．shu．edu．cn 2(清华大学热能系) 摘要：窄化流场中，不町骶液体与可压气体并存；流场中的相变导致了复杂的热物理过程；汽泡 壁附近密度变化大，流场不连续。现有的计算流体力学理论与方法无法对如此复的物理现象准确 模拟。本文采用基于压力、速度与焓为原始变量的预处理方法，用能量方程取代传统的非物理的 空化模型，用Roe格式捕获问断面，构建既适用于饱和区、也适用于液体及气体区的饱和参数求 解方法，建立一套解决空化流动问题的新方法，算例表明，新方法可以捕获流场中的空泡。 关键词：流场、空化、状态方程 ，， …、 1、前言 空化问题研究，包含了不少前沿的基础理论问题。1)空泡与其周罔的水之间的流 场结构复杂，空泡内是可压缩的气体，空泡外是几乎不可压缩的液体；2)随着水下运动 体的高速前进，空泡不断更新，是一个典型的非定常过程；3)水下潜体在高速前进时， 相变导致了复杂的热物理过程：4)在空泡壁面上，密度变化人，流场存在很强的间断面。 现有的流体力学理论与计算方法无法从根本上解决这类复杂问题。 与紊流问题相当，空化流动也是一个非常复杂的研究领域，对空化问题的理论研 究，目前有很多的方法，但大多都需要很多的假定与简化，与实际的物理问题相差很大。 基于均相平衡的空化流动模型，求解非定常的雷诺平均的N-S方程，获得水下潜体周围 空化流场中的非定常两相流动结构n吲是当今国内外空化问题研究的主要方法。但这一 方法有其局限性，即它放弃了流动过程中，特别是空化过程中必须遵守的物理规律一 能量守恒定律，而引入相平衡方程及非物理的空化模型。为什么说该空化模型是非物理 的呢?因为在饱和区，饱和压力是饱和温度的函数，饱和区内工质的干度或密度是焓的 函数，与压力无关，但该空化模型人为强制建立了相变区域内热物性上本来就不存在的 干度或密度与压力的关系，导致这一方法在理论上不完整。当然，紊流与空化有时是密 切相关的，也有不少学者研究紊流模式对空化的影响隅_’。 事实上，空化过程是汽液两相之间通过质量与能量传递来实现的，．放弃能量守恒 定律来研究空化问题显然是不完整的。本文采用基于压力、速度与焓为原始变量的预处 理方法，用能量方程取代的空化模型，用Roe格式捕获间断面，构建既适用于饱和区、 也适用于液体及气体区的状态参数求解方法，建立一套解决空化流动问题的新方法。． 123 2、基本方程及离散格式 基于三维曲线坐标系下的守恒型雷诺平均Navier—Stokes方程，以压力、速度与焓 (P．U，V，W。h)为原始变量，构造出新的流体力学预处理方法…。由连续方程与能苗方 程求得的压力与焓，反映了空化流动过程中的质量与能鼙转换规律，由压力与焓，无论 在饱和区(两相区)，还是在单相区，都可以唯一确定工质的比容、温度笆其他状态参 数，不需要引入非物理的空化模型。 (1) Ot aQa孳jaQ8q8Q8 aQa毛 aQ8q aQa毛 l二+一二+～二‘一二=一～+一一+一U．， 其中，Q=(p，甜，Kw，磊)，r预处理矩阵，E、F和G分别是善，7和f方向的无粘通量， F、Fv和G’分别是f，，7和f方向上的粘性通量，D为源项，为了节省篇幅，r、E、F、G、 的详细表述参见文献㈣，E”、F。和G’的详细表述参见相关论文或书籍。已经证明，方程 组(1)不管在高速还是低速条件下，都不会出现奇异．始终能够通过时间推进法获得 所求的解。 在计算域中，考虑一个有限体积倒梢，其中心点位于(i，j，k)，体积为d