哈尔滨工业大学计算传热学第八章求解椭圆型问题的涡量流函数法分析报告.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个方程 *