传热与流体流动的数值计算(69章).pptVIP

传热与流体流动的数值计算 [美] S.V. 帕坦卡 著 同济大学机械工程学院 朱 彤 第六章 流场的计算 6-1 制订一个特殊程序的必要性 主要的困难 计算速度场真正的困难在于未知压力场。 压力梯度构成动量方程中源项的一部分，然而还没有一个可以用来求得压力场的明显方程。 以涡量为基础的方法 二维问题中，通过交叉微分从两个动量方程中消去压力，就可以导出一个涡量输运方程。 把涡量的定义与二维稳态情况下的流函数定义相结合就是“流函数/涡量法”。 优点：压力消失了，代替处理连续性方程和两个动量方程问题，只需解两个方程求得流函数和涡量。 缺点：壁面上的涡流量值很难给定；由涡量求出压力；此方法不能外推到三维情况。 6-2 某些有关的困难 压力梯度项的表达 二维情况 6-3 一种解决困难的妙法——交错的网格 把速度分量放在与所有其它变量不同的网格上 采用交错网格时，速度分量是对位于控制容积表面上的点进行计算的。 交错网格的直接结果是：无需对有关的速度分量进行任何内插运算，就可以计算出通过控制容积表面的质量流量。 优点： 对于典型的控制容积，离散化的连续性方程将含有相邻速度分量的差，避免了图6.4波形速度场会满足连续性方程的情况，只有“合理”的速度场才有可能满足连续性方程； 两相邻的网格点之间的压力差现在成了位于这两个网格点之间的速度分量的自然驱动力。 代价 计算机程序必须提供有关速度分量位置的全部指令值和几何信息； 进行相当冗长的内插。 6-4 动量方程 如果压力场已知，就可以应用对通用变量所推得的公式求得动量方程的解。采用交错网格使离散化动量方程与主网格点处计算的其他变量φ略有不同。 x方向动量方程的一个交错的控制容积如图。 差值pP-pE可以用来计 算作用在速度u的控 制容积上的压力合力。 X方向动量方程离散化 Y方向动量方程离散化 对于图6.7中二维情况，u的四个相邻点示于控制容积的外面。 对于三维情况，应当包括u的六个相邻点。 下图表示对y方向动量方程的控制容积，它在y方向是错位的： 6-5 压力和速度的修正 为了改进压力估计值p*，假定正确的压力p=p*+p’。其中p’称为压力修正。类似有，速度修正u’，v’，w’。 略去右式第一项，得到速度修正公式： 6-6 压力修正方程 连续性方程：（假设密度ρ与压力没有直接关系） 在图示阴影部分的控制容积内积分 应用速度修正公式来代替所有的速度分量。整理得到对p’的离散化方程： 6-7 “SIMPLE”算法 Semi-Implicit Method for Pressure-Linked Equations（解压力耦合方程的半隐式法） 计算进行的顺序 估计压力场p*。 求解动量方程一得到u*、v*、w*。 解p’方程。 以p’加p*计算p。 利用速度修正公式由带星号的速度值计算u、v、w。 求解那些通过源项、流体特性等影响流场的其他一些物理量φ（如温度、浓度以及紊度等物理量）的离散化方程。 把经过修正的压力p处理成一个新的估计的压力p*，返回到第二步，重复全部过程，直至求得收敛的解为止。 压力修正方程的讨论 讨论在推导速度修正公式（6.17）中去掉项Σanbunb’的目的，并证实这样做不会对计算带来明显危害。 略去它能够把p’方程写成通用的φ方程相同的形式，并采用一个依次的，每次求解一个变量的解法。 该项代表压力修正对速度的间接或隐式的影响，这里采用的仅仅是部分的、而不是完全的隐式方案。 由SIMPLE所给定的收敛解并不包含任何忽略此项所产生的误差。 “最后”一次迭代得到的收敛的解不会受推导p’方程时所作的任何近似的影响。 质量源b可以用作流体流动解的收敛性的一个有用指标，迭代应当进行到计算域内各处b都变得足够小为止。 可以把压力修正方程看成为仅仅把我们引向正确压力场的一个中间算法，对最终的解没有直接影响。 但过程的收敛速度将与p’方程的特定公式有关。略去项过多可能导致发散。 除非采用某种欠松弛，否则6.6中所推导出来的压力修正方程也会趋于发散。 每次迭代期间速度并非保持在它们的带星号条件下，而是要用速度修正公式不断进行修正。 推导p’方程中，认为密度ρ已知，并且没有考虑压力对密度的影响。 P’方程非常象是对流传导问题的离散化方程。 P’方程的导热形性质意味着对于任何一种空间坐标该方程都不具备单向的特性。 压力修正方程的边界条件 已知边界上的压力（速度未知） 此时边界上的p’值为0，类似于热传导问题中给定温度边界条件的情况。 规定垂直于边界的速度分量 压力的相对特性 讨论一个常密度稳态流场，其中在所有边界位置上的法向速度都是已知的。则由p’方程来确定p’的绝对值是没有任何意义的。只有压力的差值不会因p’场上加上一个任意常数而有所变化，