基于液体的统一对流扩散型流体力学方程组.pdfVIP

中国T程热物理学会 流体机械 学术会议论文 编号：067045 基于液体的统一’ 对流扩散型流体力学方程组 黄典贵 (上海人学力学研究所上海200072) Tel：021 dghuang@staff．shu．edu．cn 摘要：为了表达上的方便及求解格式的统一，通常采用统一的方程形式来表达连续方程-动量方 程、能量方程、湍动能方程和耗散方程等。除了连续方程外，其他方程都可以写成对流扩散方程 的形式，由于没有扩散项，连续方程比较特别，也相对小便处理。在液体区，通过合理的数学推 导，不作任何近似、假定与简化，本文得到一套全新的连续方程形式。该连续方程以压力为未知 变量，是对流扩散型的，使得所有的流体动力学方程组都具有完全统一的方程形式。 关键词：流动，对流扩散方程， 流体力学 1、前言 ’ ．压力场是驱动流场运动的一个非常重要的参数，然而压力项通常是流场求解过程中 的一个相当麻烦的项，早期二维流场求解时，通常应用流函数一涡函数方法，将压力项 消去，解得流函数与涡函数，然后得到速度场，爵由此间接推算出压力场，这是一种非 常成功的求解方法。可对于三维流动问题，由于三维流函数不便表达，一般只能采朋其 他的办法。 压力修正方法是求解不可压或低速流动的主要方法，该方法中，运动方程独立求解， 然后求出使连续方程得以满足所需要的压力修正量，但压力校正法的收敛速度与计算稳 定性远比不上流函数一涡函数法：对于可压和高马赫数流动问题，传统的计算流体力学 通常采用时间推进法来求解可压N—S方程组，在跨音或超音流场求解时，只要采用的计 算格式适当，时间推进的方法通常能够收敛，但对于低速或不可压流动，由丁流动速度与 音速之比很小，时间相关的N—S方程组中的对流项系数矩阵的刚度变得1卜常人，方程组趋 予奇异，时间推进法使得计算收敛很慢或雉以收敛¨一’。 事实上，绝对不可压的液体是不存在的，本文通过合理的数学推导，不需要作任何 简化与假定，在的液体区，’将连续方程转换为以压力为变量的对流扩散方程的形式，使 得所有的方程都具有完全统一的形式，可以应刚成熟的对流扩散问题解法来求解这套新 的流体力学方程组。 连续方程、动量方程、能量方程共五个方程13“’，经整理可以得到如下的对流扩散 统一形式： 134 号手+(旷一古Vr，]·V矽=与}△≯+{} 其中p为密度，旷：“7-+谚+w云，u、v和w分别为x，y。z方向上的速度，(7．，歹，云) 是单位欠餐。巾为通用变量，r．为通用扩散系数，S◆为源项。 对5个方程，通用变量巾，通用扩散系数r．和源项s◆的表达如下表： 表l 通用变量巾，遥用扩散系数r．和源项s． 方程 巾 ’ r． S◆ 连续 p +· E=o 砖=一p2V·矿 方程 X向动 U 巧=p+p。 量方程 $=一罢+譬(罢)·雩7塞·a瑟t：舐 Y向动 V e=p+p。+ 量方程 s．；：-考+丢(磅心(q豺昙(q刳 Z向动 C=P+P。 量方程 W s：=一鲁+未(。卦言(。务三(。鲁) 能量 h 方程 o2昔蚩 ．$=(号手+材罢+V考+w警)+(r了G+胪)， 其中：P为压力，h