求解椭圆型流动与换热问题的原始变量法-热流问题的数值计算-课件-06.ppt

求解回流问题流场的基本方法： 求解有回流问题的流场时，可以用速度、压力（或密度）作为基本变量，成为原始变量法。分为以密度为基本变量（即以u, v, ?为变量）及以压力为基本变量，即以u, v, p为变量。以密度为基本变量时，连续方程是求解密度的控制方程，解出密度后用状态方程去求压力。 6.1动量方程的源项及流场求解中的关键问题 6.1.1 动量方程中的源项 6.1.2不可压缩流动求解中的关键问题 压力梯度 项的离散为动量方程数值求解的关键问题之一。 若采用均分网格，将各项取中心差分，得： ⑵压力的一阶导数以源项的形式出现在动量方程中。 采用分离式求解各变量的离散方程时，由于压力没有独立的方程，需要设计一种专门的方法，以使在迭代求解过程中压力的值能不断的得到改进。 6.2 交叉网格及动量方程的离散 为解决流场中不合理的压力场的检测应当使相邻两点间的压差出现在动量离散方程。 为保证计算的准确度及对压力特性的模拟，相邻两点间的压差应当是压力梯度中心差分的表达式的组成部分。 6.2.1交错网格(staggered-grid)上速度分量位置的安排 6.2.2交叉网格上动量方程的离散 交错网格系统中，关于u, v的离散方程可通过对u, v各自的控制容积作积分得出。此时压力梯度的离散形式对ue为 ，ve为 ，相邻两点间的压力差构成了 ， 。 对ue控制容积： 6.6.2 SIMPLEC算法（SIMPLE Consistent） SIMPLE算法中，速度修正完全归结为由压力项的直接作用，忽略速度修正的作用，导致了不协调。 SIMPLEC算法的收敛性远优于SIMPLE。 SIMPLEC的计算步骤与SIMPLE基本相同，区别在于： 1.以 代替SIMPLE算法中的 ； 2.SIMPLEC算法中，不再亚松弛 Date修正方案的基本思想： 1.每一层次上动量离散方程的求解分成两个阶段进行，预测阶段和校正阶段。预测阶段采用SIMPLE算法 得出压力修正值 及速度 ，从而获得压力及速度的预测值： 本次迭代计算所获得的压力与速度值为：  把上两式代入连续性方程的离散形式可得有关的代数方程：  6.6.5 SIMPLEST算法 SIMPLEST的计算步骤与SIMPLE相同，只是在SIMPLEST中对于对流-扩散项的格式作了明确的规定： 6.6.6 不同算法的比较 在SIMPLE，SIMPLER，SIMPLEC中，SIMPLE能获得收敛解的松弛药捻子变化范围较小。 在收敛速度上，SIMPLER，SIMPLEC优于SIMPLE，SIMPLEX应用不广泛。 各种算法的特点见下表： 6.7加速SIMPLE系列算法收敛的一些方法  Patankar认为：速度亚松弛因子au, av可取0.5，而压力亚松弛因子ap可取0.8。 一般认为：常数c取1或1.1，au值尽可能地取大，一般取0.7~0.8。 6.7.2显式修正法（explicit correction step） 获得一个层次的 后，把它们代入动量离散方程的右端。考虑到实际计算中把亚松弛组织到代数方程的求解过程中，因而此时动量离散方程为： Eu表示u方程的时步倍率。另一方面由p*解出u*的方程为：  略去邻点速度修正值项，得速度修正值计算公式为：  将式（6）~（7）代入连续性方程得离散形式，可得 （2） 附加源项法 附加源项法中把法向一阶导数为零视为热流密度给定的特例。采用附加源项法迭代收敛速度要比边界值更新法快得多。 以上两种方法仅适用于出口边界无回流的情形 6.9封闭腔内自然对流换热的控制方程 1. Boussinesq假设 流体中的粘性耗散略而不计； 除密度外其它物性为常数； 对密度仅考虑动量方程中与体积有关的项，其余各项中的密度作为常数。 2. ?封闭腔内自然对流换热的控制方程 采用Boussinesq假设，因 ，整理得 6.7.1选择合适的松弛因子之值 引入参数E，使得 相应于?由0变到1，E由0变到，E的变化范围比?大得多，E称为