针对冒落采空区风流移动同时存在紊流、过渡流、层流的非线性渗流问题,结合有限元方法,提出以变渗透系数的达西(Darcy)渗流反复迭代求得非线性渗流。计算中根据迭代方法的特点,为改进迭代算法振荡性收敛效果,引入动态收敛因子(即初始收敛因子和阻尼收敛因子),讨论了收敛因子对收敛过程的影响。通过合理联合选取初始收敛因子和阻尼收敛因子,使迭代的收敛速度更快。与Darcy渗流相比,采空区非线性渗流速度场趋向于平缓,计算结果与实际流态更接近。得到采空区风压等值线与流函数线正交的流网、漏风强度分布(速度场),给出了采空区流态变化与过渡的分布结果。迭代方法满足工程要求。
In allusion to the nonlinear seepage problem that turbulent flow, transition flow laminar flow coexist with the movement of airflow in a goaf, combined with finite element method, an iterative algorithm to solve nonlinear seepage equation based on Darcy seepage with variable permeability is presented. According to characters of iterative method, a dynamic convergence factor in terms of initial convergence factor and damping convergence factor is introduced to improve the oscillation convergence effect of iterative algorithm. The influence of convergence factor on convergence process is discussed, and the convergence velocity can be accelerated by rationally choosing the initial convergence factor and the damping convergence factor. In comparison to the Darcy seepage, the goaf nonlinear seepage field tends to mild, and the computing result is closer to real flow pattern. The orthogonal flow net formed by air pressure isoline and flow function isoline of goaf, the distribution of air leakage intensity (velocity field) are obtained. At the same time, the distribution result of the variation and transition of flow pattern in goaf is achieved.