二维方腔流动是不可压缩黏性的典型流动,可以用来检验各种数值算法计算精度和可靠性,目前尚不能求得它的解析解。基于Matlab编程,采用SIMPLE算法求解二维方腔流动,得到流动达到稳定状态时各物理量的分布。结果表明:随着雷诺数增大,流函数最大值减小,最大流函数的位置也更趋于方腔的中心。SIMPLE算法可以有效求解方腔流动,具有较高的求解精度。