Вестник МГТУ. 2016, №1.1.

Мельников Н. Н. и др. Оценка области оттаивания криолитозоны… 30 Анализ результатов расчетов Приведем несколько примеров результатов численных экспериментов. Изображение динамики (с шагом 4 года) пространственного распределения изоповерхностей температуры в сечении, проходящем практически через середину реакторного зала, представлено на рис. 2. В данном расчете коэффициент теплопроводности обделки равен 2 Вт/(м·К), что позволяет получить наиболее презентабельное изображение. Хорошо виден постепенный разогрев ММГП и переход льда порового пространства породы в жидкую фазу вокруг реакторного модуля с максимальной скоростью разогрева на начальной стадии эксплуатации АСММ и постепенным снижением скорости разогрева в последующие годы. При меньших значениях коэффициента теплопроводности обделки (0,5 и 0,05 Вт/(м·К)) области разогрева, а значит, и оттаивания, заметно уменьшаются, графическая картина становится менее презентабельной. Более наглядными являются другие виды графиков, посредством которых можно получить интересные количественные показатели. Например, на рис. 3 (с шагом 2 года) представлена динамика пространственного распределения температуры вдоль оси Y от внешней границы области моделирования до внешней стенки модуля при вариации коэффициента теплопроводности обделки. При этом рис. 3, а–в соответствуют указанному выше базовому значению пористости ММГП, т. е. 0,05, а рис. 3, г – увеличенному вдвое, т. е. 0,10. Естественно, что пространственное распределение температуры вдоль оси Y с противоположной стороны модуля будет идентичным. Аналогичные распределения температуры вдоль оси Х представлены на рис. 4 (с шагом 2 года). Отличие в поведении кривых рис. 3 и 4 имеется и объясняется различием в размерах сторон модулей вдоль осей Х и Y . Как видно на рис. 1, максимальная протяженность подземной выработки отвечает оси Y . Поэтому из физики процесса следует ожидать максимальные глубины оттаивания вдоль оси Х . 4 года 8 лет 12 лет 16 лет Рис. 2. Динамика пространственного распределения температуры в сечении модели через середину реакторного модуля (коэффициент теплопроводности обделки 2 Вт/(м ⋅ К))