双重时间步预处理方法及其在空化流场求解中的应用研究.pdfVIP

中国．T．程热物理学会 叶轮机械气动热力学 学术会议论文 ． 编廿．：062027 双重时间步预处理方法及其在空化 流场求解中的应用 黄典贵 (上海大学力学研究所上海200072) Tel：021 dghuang@staff．shu．edu．cn 摘要：预处理技术可以统一高速与低速、可压与不可压流场的求解方法。空化流场中，不可压液 体与可压气体流场并存，流场中的相变导致了复杂的热物理过程。本文采用基于压力、速度与焓 为原始变量的预处理方法，用能量方程取代传统的非物理的空化模型，用Roe格式捕获间断面， 采用双重时问步方法来建立一套解决非定常空化流动问题的方法，算例表明，该方法可以捕获流 场中的非定常空泡。 关键词：非定常流场、空化、状态方程 1、前言 一 预处理技术是一种将高速流动、低速流动、可压流动与不可压流动统一起来的流体 力学计算方法，目前已经得到了广泛应用¨~51，对定常流动问题，这一方法通过引入伪 ．时间步而实现时间推进。目前的预处理技术大多以压力、速度与温度为原始变量，研究 表明，这一技术在求解无相变的流场时，效果是很好的。但是，当流场中发生相变时， 在相变区域，饱和压力与只是温度的函数，与密度无关，此时，将线性相关的压力与温 度作为独立的自变量进行空化流场的求解显然是不行的。本文作者提出了基于压力、速 度与焓为原始变量的预处理方法哺!，由压力与焓，可以构建液体区与饱和区各状态参数 之间的关系…。 对于非定常流动问题，采用双重时间策略，物理时间层用米跟踪流场的物理变化， 伪时间层用来进行每个物理时间步内的迭代哺吲。本文试图将双重时间策略引入以压力、 速度与焓为原始变量的预处理方法中，进行非定常空化流场的求解。 2、基本方程 (p，t／，y，w．h)为原始变量，构造出流体力学预处理方法‘¨。 I二+一——+一上+一二=一一+一一+一上+ at aQa芎aQaJ7 aQ8乞 aQ8专 aQ8n aQa专 677 其中，Q--(p，／．／，v，w’办)。r预处理矩阵，E、F和G分别是善，刁和f方向的无粘通鼙， E’、F’和G’分别是孝，刀和f方向上的粘性通量．O为源项，为了：1了省篇幅，r、E、F、G、 的详细表述参见文献‘61，E’、F‘和G’的详细表述参j^!‘相关论文或书籍。已经证明，方程 组(1)不管在高速还是低速条件下，都不会出现奇异，始终能够通过时间推进法获得 所求的解。 对于定常流场的求解，r为伪时间层，方程组(1)的对流项采用Roe格式离散， 扩散项采用中心差分离科们，在伪时间步上，采用一阶差分，在时间上向前推进，就可 以获得定常解。 r—AFQt+i+唾层‘+t+磊，t+-+吒G“-：如 (2) 这里，AQ‘?1=Q¨一Q‘，名E川=名E‘+彳‘唾(△QM)， 磊GM=磊F‘+B‘磊(△Q¨)，吒G川=露Gj+c‘毛(△Q川)， R跗=D‘+忙’一E”+F’一F’+G’一GV广 万表示空间差分，下标为差分的方向，A，B，C分别表示对应于原始变量Q的通量Jacobian 矩阵，k是伪时间步数。 对于非定常流场的求解，增加一个物理时间层，用二阶精度的后差格式离散物理时 间层陬引。 · ．(3)‘(3) r——+——+艿，¨1+万。“1+艿，“1=足‘，c … △f