笔者应用三维k -ε紊流模型描述压入式局部通风工作面风流的流动过程 ,并采用控制容积法导出了描述流体流动方程的离散化方程式 ;用计算流体力学的方法求得了三维k -ε紊流模型的数值解。采用SIMPLE(压力耦合方程式的半阴解法 )算法解算流场 ;TDMA(三对角线算法 )和Gauss Seidel法结合通过线顺法求解离散方程。利用弱松弛法防止非线性方程组迭代求解过程中的发散现象 ,获得了掘进工作面空间的速度分布 ;并将模拟解算结果与在具有实际尺寸的拱形巷道模型中测得的实验结果进行了对比。数值计算结果与模型实验测定结果非常一致 ,验证了笔者的数学模型和数值模拟方法的正确性。