哈尔滨工业大学计算传热学第八章求解椭圆型问题的涡量流函数法.ppt

§8-1 涡量-流函数方程 * 第八章 求解椭圆型问题的涡量-流函数法 问题：原始变量法中压力项处理的困难。 ① 为避免锯齿、棋盘压力场采用交错网格。 ② 交错网格引起u、v与P的耦合关系。 在很多情况中压力不是重要参数，尤其是换热计算并不关心压力分布。 5个方程—5个未知变量 §8-2 涡量-流函数方程的离散 界面流速的插值 解涡量方程和能量方程时要用到控制容积界面的速度 §8-3 涡量-流函数离散方程的求解 离散方程求解步骤 ----求物性参数 8-4 B2 B3 B1 B6 B5 B4 §8-4　ω、ψ边界条件的确定 Ψ 可根据直观的物理意义确定边界，而ω的物理意义并不直观。 一般意义边界条件举例 一、壁面边界B1 B6 １、ψ 如果壁面　u＝v＝０ 则ψ＝常数 当只有一个壁面（相联）时， 取ψ＝０ 当有多个不相联壁面时，一个面取ψ＝０，其他的根据 总流量值来确定 ２、ω 一种办法：台劳级数展开 而 在壁面上 故 即 （１８） 将涡量用流函数表示，公式（１８）具有一阶精度。 ２、移动边界 实际无限远，人为确定一个有限的计算域B3，仍应距绕流物体足够远。 y ∞ B3 二、上边界 B3 有两种情况 １、静止壁面，则其与B1相同。 ωB3的计算与B1相同 于是 即ψB3＝常数 0 B3 因为对称 故υB5≡0 故 三、对称边界B５ １、无旋 因为对称线是流量，故ψB5＝常数 或 ψB5＝０ 且 所以 四、进口边界B2 假定 1 2 3 j+1 j j+2 2、有旋 给定ψ（y）B2 但　　　　很难表达 假设 则 所以 L3 L2 L1 j 五、出口边界B4 １、充分发展 所以 即 即 ２、 ３、 规定函数的导数，限制较少，导数阶数越高，限制越少． 六、尖角 １、凹尖 　　根据物理意义 ２、凸角 　　Ψc＝常量，根据上游条件确定 （１）双值法 　C点即可作为B1面的点，又可作为B6面的点，可以分别按壁面涡公式写 固C是奇点，可以不是单值，上游离散用　　，下游用 a (j,j+1) (i,j) (i+1,j) c B1 B6 （２）平均法 （３）分离点法 　　从上游看C点没有分离 从下游看C点没有分离处，流动已分离，不是壁面涡 §8-5 三维流动问题的涡量-流函数方程 矢量：包括3个分量，即3个方程 *