传热与流体流动数值计算(6~9课)-精品文档.pptVIP

传热与流体流动的数值计算;第六章 流场的计算;6-2 某些有关的困难;;压力梯度项的表达 二维情况;积分连续性方程;离散化的连续性方程要求相间点而不是相邻点上的速度相等。 其结果是完全不合乎实际速度场却满足离散化连续方程。在二、三维情况下，对所有速度分量都可以构成类似分布图形，速度分量满足连续性方程但却不能认为是合理的或有意义的解。 ;把速度分量放在与所有其它变量不同的网格上 采用交错网格时，速度分量是对位于控制容积表面上的点进行计算的。;交错网格的直接结果是：无需对有关的速度分量进行任何内插运算，就可以计算出通过控制容积表面的质量流量。 优点： 对于典型的控制容积，离散化的连续性方程将含有相邻速度分量的差，避免了图6.4波形速度场会满足连续性方程的情况，只有“合理”的速度场才有可能满足连续性方程； 两相邻的网格点之间的压力差现在成了位于这两个网格点之间的速度分量的自然驱动力。 代价 计算机程序必须提供有关速度分量位置的全部指令值和几何信息； 进行相当冗长的内插。;如果压力场已知，就可以应用对通用变量所推得的公式求得动量方程的解。采用交错网格使离散化动量方程与主网格点处计算的其他变量φ略有不同。 x方向动量方程的一个交错的控制容积如图。 差值pP-pE可以用来计 算作用在速度u的控 制容积上的压力合力。;X方向动量方程离散化;Y方向动量方程离散化;对于图6.7中二维情况，u的四个相邻点示于控制容积的外面。 对于三维情况，应当包括u的六个相邻点。 下图表示对y方向动量方程的控制容积，它在y方向是错位的：;为了改进压力估计值p*，假定正确的压力p=p*+p’。其中p’称为压力修正。类似有，速度修正u’，v’，w’。 略去右式第一项，得到速度修正公式：;连续性方程：（假设密度ρ与压力没有直接关系） 在图示阴影部分的控制容积内积分;应用速度修正公式来代替所有的速度分量。整理得到对p’的离散化方程：;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’场上加上一个任意常数而有所变化，因而压力是一个相对变量，而非绝对变量。 如果p’的绝对值不是单值的，接待数方程的迭代法收敛于一个解，它的绝对值取决于初始的估计值。但直接解法会遇到一个退化矩阵，不能给出一个解。解决办法是任意规定一个控制容积内的p’值，而求解对其余控制容积的p’值。 观察不确定p’场令一个方法是注意：对所有控制容积的连续性方程并不是一组线性独立的方程组。 在许多实际问题中，绝对压力值要比所遇到的局部压差值大得多。最好是把在一个适宜的网格点处的压力p取为0，并把它作为该参考值的相对压力。;问题的提出 大多