可压缩流体力学问题地有限分析格式.pdfVIP

中国T程热物理学会 叶轮机械气动热力学 第f。一届年会 编号：052047 可压缩流体力学问题的有限分析格式 黄典贵 (上海大学力学研究所上海200072) Tel：021 dghuang@staff．shu．edu．cn 摘要：求解对流扩散方程的有限分析格式具有自动迎风性质、计算稳定性好、收敛快等特点。对 于可压缩问题，通过合理的数学推导．不作任何近似与简化．就可以将原始的连续方程转化为以 压力表示的对流扩散方程的形式，使得所有的流体动力学方程都具有完全统一的对流扩散形式， 并培出了这种三维对流扩散方程组的有限精确差分计算格式。 关键词：可艇缩流动，对流扩散方程，流体力学，有限分析格式 1 前言 对于经典的流体力学偏微分方程组，有限差分格式用数值微商替代微分，导出相应 的线性代数方程组：有限元法从定解问题的弱解形式(变分形式)出发，将古典的 是在局部单元上线性化微分方程，求局部单元上的精确解，从而构成整体的线性代数方 程组。研究表明，他具有自动迎风性质，计算稳定性好，收敛快等特点…。 压力项通常是粘性流体力学求解中比较麻烦的一项。二维流动问题中的流函数一涡 函数方法，通过消去压力项，解得流函数与涡函数分布，然后得到速度场及压力场。对 于三维流动问题，由于三维流函数不存在，流函数一涡函数方法不再适用。 对不可压或低马赫数流动问题，传统的计算流体力学大多是采用压力修正方法来求 解的，在这～方法中，运动方程独立求解，然后求出使连续方程得以满足所需要的压力修 正鲑，但压力校正法的收敛速度与计算稳定性远比不上流函数一涡函数法：另一方面，对 于可压和高马赫数流动问题，传统的计算流体力学通常采用时间推进法来求解可压N～s 方程组。当求解跨音或超音流动问题时，只要采用的计算格式适当，时间推进的方法通常 能够能够收敛，但对丁低速或不可压流动，由于流动速度与音速之比很小，时间相关的 N—s方程组中的对流项系数矩阵的刚度变得非常大，方程组趋于奇异，时间推进法使得计 算收敛很慢或难以收敛“’…。 2 控制方程 连续方程、动龄方科、能鼙方程及k一￡两方稗湍流模型方程。、”共七个方程，经整 理可以得到如F的对流扩散统一形式： 水动力学囝防科技重点实验空摧金资助项目，卜海市教委科技发展基会资助项目。 167 卟獬谢倒悟 E2= [(窆)2+(鲁]2+(舅)2]告 E3 ㈣博+剀詈 Ed b愧蒜埘∽；蓑棚祁s斜鲁 设定计算域({，n)为～个单位体积的正方体，对式(趴计算域的二个方向上的网 格间距分别为(6{，6。，6‘)，令算子： 陪以：寒蝎筹 卜一B；嚣慢筹 (4) 卜一c；是地等 式(3)可以写成： 尝也眦卅¨^忡。+鲁 (5) L；巾?“=j百：：；}i了k^巾；1二’～2ch(A．)巾?+1+e一^咖篙1] L．十j“=ji：；最吾了}8’十茹1—2ch(B．Ⅻ?“+e—BJm篙1] (6) c。十：：：J Le十：“2j吾晶[e。‘十：：}～2ch(c。如：“+e 迎^m删憷卦列吣m呲彻梧H戊”等，B，=警 169