Вестник МГТУ. 2017, том 20, № 1/1.

Козырев С. А. и др. Методический подход к оценке времени проветривания… 6 Рис. 1. Геометрия моделируемой области (основная выработка 30×5×5 м; тупиковая выработка 5×15×5 м) Fig. 1. Geometry of the modelled area (the main excavation 30×5×5 м; the blind excavation 5×15×5 м) Численная реализация Как и в предыдущих работах авторов и их коллег [1–4], задача численного моделирования решалась в два этапа. На первом этапе выполнялся расчет аэродинамических характеристик (пространственные распределения скоростного поля и коэффициентов турбулентной вязкости), а на втором – решается нестационарное конвективно-диффузионное уравнение, описывающее процесс выноса пассивной примеси из загазованной части тупиковой выработки до достижения уровня ПДК. Для вычисления аэродинамических характеристик применяется стандартная ( k – ε)-модель турбулентности. В качестве примера результатов расчетов на первом этапе на рис. 2 изображено скоростное поле в форме векторов, а также по отдельным сечениям приведено распределение турбулентной динамической вязкости. Численные значения (см. легенду) свидетельствуют о том, что в тупиковой выработке прогнозируется именно турбулентный режим (расчетные значения на несколько порядков больше динамической вязкости воздуха). Рис. 2. Расчетное поле скорости и распределение (по отдельным сечениям) турбулентной динамической вязкости Fig. 2. Calculated velocity field and distribution (for individual sections) of turbulent dynamic viscosity На рис. 3 приведено распределение турбулентной динамической вязкости в сечении на полувысоте выработок в плоской постановке. Можно отметить существенную неоднородность поля скорости в тупиковой