高等流体力学第九章..ppt

第九章 流场的计算 9.1 压力梯度项的表达 上式可看出动量方程包含两个相间而不是相邻网格点之间的压力差，在效果上压力取的是较实际所采用的网格为粗的网格上的值，导致精度下降 9.2 锯齿形压力场 在任一网格点P处，对应的pW－pE＝0。结果是错误的 9.3 棋盘形压力场 9.4 连续性方程 在如图所示的控制容积内积分 9.5 波形速度场 这些速度满足连续性方程但却不能认为是合理的或有意义的 9.6 交错网格 如果不一定要在同样的网格点上计算所有的变量，而对每一个因变量采用不同的网格，成为一种交错网格。 采用交错网格时，速度分量是对位于控制容积表面上的点进行计算 9.7 速度分量交错网格 9.7 速度分量交错网格 9.8 SIMPLE算法 SIMPLE——用以进行流场计算的程序，是解压力耦合方程的半隐式方法。 9.8 SIMPLE算法 SIMPLE算法的基本思想： 对于给定的压力场，求解离散形式的动量方程，得出速度场，对压力场加以修正，据此求得新的速度场，检查是否收敛，若不收敛，用修正后的压力值作为给定的压力场，开始下一层次的计算，如此反复，直到获得收敛解。 SIMPLE算法的两个关键问题： ①如何构造压力修正方程？ ②如何构造速度修正方程？ 9.8.1 交错网格关于速度的 和 的动量方程的离散形式： （1） （2） 速度修正方程 9.8.1 速度修正方程 设有初始的猜测压力场 ，根据前两式有： （3） （4） 定义压力修正值 ， （5） 9.8.1 速度修正方程 定义速度修正值 和 ，则有： （6） （7） 将正确的压力场 代入动量离散方程（1）（2），得到正确的速度场（ ）。 9.8.1 速度修正方程 假定源项 不变，方程（1）和（2）减去（3）和（4），得： （8） （9） 9.8.1 速度修正方程 引入压力修正值与速度修正值的表达式（5）（6）和（7），方程（8）和（9）可写成： （10） （11） 可以看出，由压力修正值可可求出速度修正值（ ）。 9.8.1 速度修正方程 为简化上式的求解过程，可近似处理：略去方程中与速度修正值相关的 和 ，有： 其中， （12） （13） （14） ， 9.8.1 速度修正方程 将（12）和（13）代入式（6）和（7），有： （15） （16） 对于 和 ，也存在类似的表达式； 上式表明，如果已知压力修整值 ，便可对猜测的速度场 作出相应的速度修正，得到正确的速度场 。 9.9 压力修正方程的讨论 在推导速度修正式 时摒弃 项 这样做不会对计算带来明显的危害。 1.略去 项可以把p’方程写成与通用的Φ方 程相同的形式，并采用一个一次的，每次求解一个 变量的解法。 2.在SIMPLE算法中，忽略 是一种认可。 3.由SIMPLE算法所给定的收敛解并不包含任何由忽 略 所产生的误差。 4.最后一次迭代之后，所有的解就是收敛的，对于所有的控制容积，压力修正方程的质量源b将实际上变为0。。所以收敛解不会受推导p’ 方程所作的任何近似的影响。 5.当计算域内各处质量源b的值变得足够小时就可以停止迭代。 6.压力方程可看成是引向正确的压力场的一个中间算法，对最终解没有直接影响。 7.过程的收敛速度与p’方程的特定公式有关，如果略去的项太多，可能会导致发散。 9.9 压力修正方程的讨论 开始 假设一个速度分布，用与计算首次迭代 时的动量离散方程中的系数和常数项 假设一个压力场，即给定压力猜测值 根据当前速度场及压力场，计算动 量离散方程等方程中的系数和常数项 步骤1：求解动量离散方程，式 与 步骤2：根据速度 ，求解 压力修正方程， 步骤3：修正压力与速度，式（5），（15），和（16） 步骤4：求解所有其他的离散化输运方程 （视需要进行） 收敛否？ 结束 令 否 是 (3) (4) (5) (15) (16) 输运方程 SIMPLE算法流程图 压力修正方程 流场数值计算的主要方法 流场计算的基本过程是在空间上用有限体积法或其他类似方法将计算域离散成许多小的体积单元，在每个体积单元上对离散后的控制方程组进行求解。 本质：对离散后的控制方程组的求解。 流场数值计算方法分类图 流场