第6章求解椭圆型流动与换热问题的原始变量法.pptVIP

马良栋2008年3月24日 求解流动与换热问题的方法 分离式求解：u, v, p 各类变量独立地进行求解，即在一组给定的代数方程系数下，先用迭代法求解一类变量而保持其他变量为常数，如此依次求解各类变量。 §6.1 流场求解中的两个关键问题 一、压力梯度的离散 一、压力梯度的离散（续） 不包含pi, j 动量离散方程无法检测出附加 在其上面的棋盘形压力场。 采用2-d 压差构造的离散格式 无法检测出不合理压力场。 压力梯度离散必须用中心差分！ 且只能用1-d 压差！ 交叉网格！ 同位网格！ 二、压力的求解 压力没有控制方程。 压力与速度的关系隐含在连续性方程中，如果压力场是正确的，则据此压力场解出的速度场必满足连续性方程。 压力与速度耦合问题： 1、如何构造求解压力场的方程？ 2、在假定初始压力分布后，如何构造计算压力修正值的方程？ §6.2 交叉网格及动量方程的离散 一、交叉网格上速度分量位置的安排 u,v, p(其他标量场及物性参数)分别存储于三套网格上 主控制容积 u 控制容积 v 控制容积 二、交叉网格上动量方程的离散 特点：(1)积分用的控制容积是u, v各自的控制容积。 (2)压力梯度项从源项中分离出来。 三、交叉网格上的插值 1、界面上的流量 ue北界面上的流量 2、界面上的密度 3、界面上的扩导（ue北界面上的扩导） 四、采用交叉网格的注意事项 1、三类变量的节点编号方法 速度矢量箭头所指向的 主节点的编号为该速度 的编号。 2、与边界相邻的速度控制容积与内部速度控制容积不同 §6.3 求解Navier-Stokes方程的压力修正方法 一、压力修正方法的基本思想 1、假设一个压力场，记为p*。 2、利用p*求解动量离散方程，得出相应的速度u*, v*。 3、利用质量守恒方程来改进压力场，要求与改进后的压力场相对应的速度场能满足连续性方程。 4、以p*+p′ ，u*+u′ ，v*+v′作为本层次的解，并据此开始下一层次的迭代计算。 关键问题： (1)如何获得压力修正值p′，使与p*+p′相对应的u*+u′，v*+v′能满足连续性方程？ (2)获得了p′后，如何确定u′ ，v′ ？ 二、速度修正值的计算公式 三、求解压力修正值的代数方程 四、压力修正值方程的边界条件 b为一个控制容积不满足连续性的剩余质量的大小，可作为收敛判据。 RMAX=max(b(i, j)) RSUM= Σ(b(i, j)) §6.4 SIMPLE算法的计算步骤及举例 SIMPLE：Semi-Implicit Method for Pressure Linked Equations CLEAR ：Coupled and Linked Equations Algorithm Revised 一、SIMPLE算法的计算步骤 (1) 假定一个速度分布（u0，v0），以此计算动量离散方程中的系数及常数项； (2) 假定一个压力场p* ； (3) 依次求解两个动量方程，得u*，v* ； (4) 求解压力修正值方程，得p′ ； (5) 根据 p′改进速度值； (6) 利用改进后的速度场求解与速度场耦合的变量； (7) 利用改进后的速度场重新计算动量离散方程的系数, 用改进后的压力场作为下一层次迭代计算的初值, 重复(3)~(7)直到收敛。 二、SIMPLE算法应用举例 用SIMPLE算法确定 假定： §6.5 SIMPLE算法的讨论及收敛判据 一、SIMPLE算法的讨论 1、SIMPLE算法的简化假定 (1)速度场u0，v0与压力p*的假定各自独立地进行。 说明：随着迭代的进行两者间会逐渐趋于协调。 (2)速度修正值计算式略去了邻点速度修正值的影响。 说明：迭代趋于收敛时，邻点速度修正值已无影响。 (3)采用线性化的动量方程，离散方程系数及b为定值。 说明：迭代趋于收敛时，两迭代层次之间的量趋于对应相等。 简化不影响收敛的解，但会影响收敛速度及健壮性。 2、p′边界条件的假定 (1)边界法向速度已知。 (2)边界压力已知。 令压力修正值方程与边界相应的系数为零，相当于“绝热型” 边界条件。 要求：计算区域满足总体质量守恒。 4、速度与压力修正值的亚松弛 p′修正速度是合适的，但修正压力则被夸大。 亚松弛： 为限制相邻两层次之间的变化，以利于非线性问题迭代收敛（避免发散）。 二、流场迭代求解收敛的判据 1、内迭代与外迭代 内迭代：在一组确定的系数及源项