全速解法在湍流跨音速流动的应用.pdfVIP

全速解法在湍流跨音速流动的应用 郭阳张景明符松 (清华大学工程力学系北京100084) 摘要本文对不可压流动的常用算法SIMPLE算法进行了推广．使其能计算从亚音速到 超音速一定马赫数范围内的流动。这里，可压缩流动和不可压缩流动的数值算法实现了统 一，称为全速解法。本文对全速解法在二维流动计算中的应用性进行了初步的研究。我们 采用了非交错网格的有限体积方法对控制方程进行离敞，并用动量插值法来求得连续方程 中单元边界上的变量值。 本文对全速解法在=维层流的计算效果进行了考核，而后又将此算法在湍流跨音速流 动中应用。计算表明本方法是成功的，能够很好的反映各种马赫数下的流场特性。 1引言 流体的运动可分为可压缩流动和不可压缩流动两种。 传统上．不可压缩流动的计算通 常采用SIMPLE算法”1，这种方法通过将连续方程转化为压力修正方程来实现压力和速度 ． 的祸合。另一方面．可压缩流动的计算则往往采用以罾度作为基本变量的通量分裂法，再 从状态方程中求得压力。 遗隧的是通量分裂法不便用于不可压，镭或低马赫数沉动的计算． 因为在这些情况下密度的变化很小．压力密度之间的耦合也变得很弱。为了避免这个困难， 人们做了很多努力，包括在连续方程中加入人工粘性作为时间微分项。但是天有在达到稳 定状态时解才是有效的。如果采用压力作为基本变量将避免基子密度的万法中可能遇到的 问题。这是因为无论流动马赫数是多少。压力的变化总是存在的。因此以压力作为基本变 量可以使汁算适用于整个马赫数范围。这样就可以用统一的算法采解决可压缩和不可压缩 流动的计算．使算法的适用范围大大拓宽。 SIMPLE算法在不可压缩流动中的应用已经比较成熟。而在可压缩流动中这种算法的研 究还不够完善。本文在现有研究成果【=l【“的基础上，加入了一些新的算法洛式．包括使用 Cartesian坐标系中的速度分量作为基本变量。在输运方程的离散化中提出多种用于求解单 。 ’ 元边界上变量值的对流插值格式p】。为了避免复杂的数据存储，我们采用了非交错网格的 动量插值法来求得连续方程中单元边界上的变量值。 这些方法被成功地运用于二维流场的计算，计算结果表明了改进后的SIMPLE算法在可 压缩流场中的灵活性和高效性。 2 数值方法 21控制方程的离散化 设流体为常比热的完全气体，其本构关系满足STOKES假设，则控制方程组为连续方 程，动量方程(N-s方程)，能量方程及完全气体的状态方程。我们将其推广到非正交曲 线坐标系中．并在定常状况下把它们改写成标准的对流扩散形式： 畦 ÷(C中一D∞)=dS中 (5) (“t 其中中可以是l，U，v，T等。因篇幅原因，方程的具体形式略去。 96 在中心为P的单元控制体上对(5)式进行积分，得到方程的通用守恒形式，然后将其中的 通量分解为对流分量和扩散分量再进行离散，可得最终离散方程的形式如下： n5W，E，S，N (6) Apffp=2A加m肋+勋 H6 2．2连续方程中压力、速度、密度的耦合 SIMPLE算法的关键是如何通过连续方程建立起合理的压力修正方程。可压缩流动引入 了密度变化，它的连续方程在离散化时与不可压缩流动很不相同。对此．我们将仔细讨论。 制体为例来加以解释； 连续方程控制体积分后得÷(∥)。一(矽)。+(尸功。一(尸叼。=0 (7) 由于我们把压力作为基本变量，所以连续方程必须转化成压力或其相关量的方程。根据 得到的通量肼’并不一定满足连续方程，所以可能会出现质量残量：