第六章——流场的计算.pptVIP

第六章流场的计算流场计算的问题：（1）运动方程中的对流项包括非线性量(速度)（2）每个速度分量既出现在运动方程中，又出现在连续方程中（3）压力项出现在运动方程中，但却没有可以直接求解压力的方程（不可压缩流体流动问题）解决的办法是迭代法解决办法是非原始变量法和原始变量法迭代法是从一个估计的速度场开始，通过迭代逐步逼近速度的收敛值非原始变量法是从控制方程中消去压力的方法原始变量法是直接以原始变量作为因变量进行流场求解，该类方法也称基本变量法原始变量法主要有SIMPLE算法,SIMPLER算法，SIMPLEC算法和PIOS算法基本变量法求解的有关困难（1）运动方程中压力梯度离散所遇到的困难（2）连续方程离散所遇到的困难WPEWPEWPE(a)压力梯度(b)控制容积(c)波状压力！！！锯齿波问题通常由一阶导数项产生解决方案——交错网格标量主控制容积交错网格示意图WEPSNWEPSN标量控制容积矢量控制容积交错网格对边界的标量进行计算时，恰好是矢量控制容积的节点，无需进行线性差值例如：在运动方程的离散过程中,控制容积为速度控制容积，压力导数的离散为：！！！此时相邻两点间的压力差构成了压力导数，从根本上解决了一般网格系统时遇到的锯齿波问题采用交错网格的代价是，所以存储于主节点的物性值在求解动量方程的速度时必须通过插值才能得到所需位置上的值运动方程的离散相比一般标量的离散，运动方程离散出现新特性：（1）对运动方程积分所用的控制容积不是主控制容积，而是各自的控制容积（2）运动方程中的压力梯度项从源项中分离出来WEPSN矢量控制容积SIMPLE算法压力与速度的修正（连续方程用于求解压力分布）对于离散的运动方程，只有压力场已知，或是按照某种方法估计出来才能求解。除非采用正确的压力场，所得的速度场将不会满足连续性方程基于估计的压力场不满足连续方程的速度场用表示。来自于求解下列方程组：对估计的压力场进行改进，其修正值为速度修正值为正确的速度为：正确的压力为：压力修正方程：压力修正应满足的条件是：根据所改进的速度场能满足连续性方程WEPSNNWNESWSE把连续方程在时间间隔内对主控制容积积分将修正后的速度表达式代入连续方程得SIMPLE算法的基本思想SIMPLE(Semi-ImplicitMethodforPressure-LinkedEquations)算法是求解压力耦合方程的半隐式法。在得到速度修正方程式的过程中，略去了速度修正项，而只保留压力修正项，去掉速度修正项就成为“半隐”，而保留这一部分时，就成为一个全隐的代数方程SIMPLE算法的计算步骤：（1）假定一个压力场p*（2）求解运动方程，得（4）利用速度修正方程得（3）求解压力修正方程得，进而算出（5）利用改进后的速度场通过源项，物性等与速度耦合的标量（6）把作为一个新的估计压力返回到第二步，重复整个过程第六章流场的计算**