SIMPLE算法求解方腔内粘性不可压流动VIP免费

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