低马赫数完全气体流动的计算方法研讨.pdfVIP

第 i-届全国计算流体力学会议论文集 低马赫数完全气体流动的计算方法研究 向大平 邓，1、刚 毛枚良 中口空气动力研宪与疚爪中心篇空气动力李研宪所 ，闷仲阳621000 摘要 ;文献川提出了一种计算低速流动的新方法，该方法去除了不可压方法要求速度场散度为零 的限制，同时也省却了可压方法中能it方程的计算 ‘在这种方法的基础上，作者进一步优化了算法，确定 了该方法马赫数的适用范围，讨论了参数口对汁算的影响，计算了大t的典型算例，获得了与其它数值计 算和实验结果都符合得很好的结果，从而证明了这种新的计算方法的合理性和有效性 关键词:低马赫数流杯计算方法数值模拟 1 引言 通常。两种策略被用来求解低马赫数流动的问题:一种方法是求解可压缩N-S方程，另 一种是求解不可压缩N-S方程。用可压缩N-S方程的困难主要在于能量方程与连续性方程、 动量方程辊合使得数值求解的代价过大，而且，在马赫数很低的情况下这种方法的收敛极为 缓慢。求解不可压N-S方程要求速度场散度处处为零，而这也不是一件十分容易的事。 由于控制方程的刚性严重，用原始变量不可压缩N-S方程求解低马赫数流动问题的收 敛十分顷难 从数学本质上分析，不可压流场中的压力场是椭圆性质而速度场是双曲性质的， 流场中某一点压力的变化向四周的传播是不需要时间的，速度影响的传播却是需要时间的。 因此，求解不可压N-S方程的关键在于在时间推进求速度场的同时求出压力场。一类通常 的求解方法就是首先用时间推进方法求解原始变量的动量方程求出速度场，然后采用迭代法 求解压力Poisson方程求出压力场。但是大量的数值试验表明Poisson方程的收敛速度极慢， 而且这种方法也引起较大的速度散度误差。不可压流动的模型建立的墓础在于很大的压力改 变只引起密度或者说速度场散度极小的改变，因此，如果计算过程中速度散度的误差较大， 必然引起很大的压力误差。 通过分析完全气体的特征，文献P]提出了一个新的计算低马赫数流动的方法。新方法 去除了不可压方法要求速度散度场为零的限制，同时也省却了可压方法中能量方程的计算。 2 低马赫数流动一个新方法的控制方程 事实上自然界的一切流动都是可压的，以往的不可压方法假设密度不变，从而要求在 计算过程中强制流场的速度散度处处为零，而这一要求恰恰就是采用不可压方法计算的最大 难点。因此不管从物理本质还是计算效率来看，不可压方法都是不尽合理的。以空气为介质， 邓小刚、庄逢甘P1分析了低马赫数流动的特点，发现温度的变化与密度的变化有相同的量级， 都正比于M二o因此提出简化低马赫数流动控制方程的另一方法是忽略状态方程中温度的变 化，从而得到了一种计算正压流体低马赫数流动的新方法， pDDpt_a`v (1) p，_DD1t;二一8p十ra,i; (2) 257 p=constant.p (3) 其中口是一常数，建议其取值范围是 15r5:1.56 (4) 3 数值方法 3.1控制方程的无量纲化及坐标变换 以“.， P.,P。和特征长度D作为无量纲化参数将方程(()一3〔)无量纲化，取(3) 中constant为 11 器=--a,v; (5) 。会一，十佘ajej, 6 3 ( P=pa ) 7 ( 气代表应变率张量· 1 式中雷诺数Re= Mm二u =‘YM a 在正交的曲线坐标系