第二讲_流动基本方程.docVIP

计算流体力学（研究生）参考材料——能源与动力学院.黄国平编 第 PAGE 14页/共 NUMPAGES 14页 第二章 基本流动方程 本项目研究需要计算进气道内外流耦合的复杂三维粘性流动，采用的流场控制方程为三维的Navier-Stokes方程组，它是连续方程、动量方程、能量方程的联立方程组。本章给出了积分形式和微分形式的N-S方程组，以及在有限体积法离散中需使用的坐标转换后N-S方程组，另外还给出了湍流计算使用的平均化后的湍流N-S方程。 §2.1 积分形式N-S方程组 在直角坐标系下，忽略重力做功和辐射传热的积分形式N-S方程组可写为如下的矢量形式： 质量守恒方程 （2-1） 动量守恒方程 （2-2） 能量守恒方程 （2-3） 这是一个关于时间的双曲型方程组，其中E代表总比能（内能和动能之和）；为能量流矢量（energy flux vector），这里假定能量流矢量仅仅表达分子能量的运输，它可以用Fourier定律来描述，即 （2-4） K为热传导系数，根据Prandtl数：为常数，及Sutherland公式： 来确定，其中和 是与流体有关的常数，详见表2-1。 对于理想气体，总比能可表述为下述张量缩并的形式： 并且总焓为： 这里ｅ代表内能，Ｈ为总焓，γ为比热比。 方程（2-2）中是粘性应力张量（viscous stress tensor），对于牛顿流体，假定应力张量随着变形率张量连续变化（Stokes假设） （2-5） 变形率张量表为 为粘性系数，层流的粘性系数与压力的关系可以忽略不计，于是可采用Sutherland公式 ，其中和 是与流体有关的常数，详见表2-1。 表2-1 气体的粘性系数和热传导系数 气体种类 （N·S/㎡） （K） （W/m·K） （K）R （J/㎏·K）空气1.458700×10-6110.5552.5010830×10-3194.4287.07氮气N21.399312×10-6106.6672.3611054×10-3166.7296.80氧气O21.751720×10-6138.8892.6954795×10-3222.2259.90 §2.2 微分形式N-S方程组 在直角坐标系下，忽略重力做功和辐射传热的微分形式N-S方程组可以按守恒型变量写为如下的矢量形式： （2-6） 其中，u是守恒变量，f、g、h 是对流通量，a、b、c 粘性通量，它们的表达式为： （2-7） 粘性通量中，应力和热流量的表达式为： （2-8） 上述公式中，t是时间，x， y， z是空间坐标，是密度，是三个坐标方向的速度分量，E 是单位质量气体的总能，p是压强，T是温度，为粘性系数，K为热传导系数。对于量热完全气体，可给出下述形式的气体状态方程： （2-9） §2.3 方程坐标变换 N-S方程组可以从笛卡尔坐标系转换到任意曲线坐标系下，其中 （2-10） 这里介绍的坐标变换是由Viviand[1]和Vinokur[2]发展的，变换的选取是如此进行，以便使曲线空间网格间隔是均匀的，且为单位长度。这将导致计算域ξ、η 和ζ为矩形域且为均匀的规则网格。所以在数值计算式中标准的非加权差分格式可以应用。原笛卡尔空间将被称为物理域，在物理域的点与计算域中的点将形成一一对应的关系（奇异点例外）。在奇异点区域，通常是一个物理点对应许多个计算点（这通常出现在计算边界上）。 按式（2-10）进行变换，对于某变量，有： 令上式的变换矩阵为J，则其行列式称为雅可比行列式，记作： 则有 变换后得守恒型N-S方程：