SIMPLE算法求解方腔内粘性不可压流动学习资料(共32页).doc

SIMPLE算法求解方腔内粘性不可压流动目录一、问题描述1二、离散格式2交错网格3方程离散4三、SIMPLE算法基本思想7边界条件处理8虚拟网格处理9方程求解10输出变量处理12SIMPLE算法流程图14四、程序中主要变量的意义14五、计算结果与讨论16函数最大值16变量等值线图17主要结论21六、源程序21一、问题描述 假设 的方腔内充满粘性不可压缩流体，左、右、下壁固定，上壁以 运动，试求 时的定常解，方腔如图1所示。图1 方腔内流动示意图二、离散格式 本算例采用求解不可压缩流动的经典算法，即SIMPLE算法，求解方腔内粘性不可压缩流体运动的定常解。SIMPLE算法的全称为Semi-Implicit Method for Pressure-Linked Equations，即求解压力关联方程的半隐式算法。采用SIMPLE算法时，为了避免中心差分格式将“棋盘”型