Численное моделирование процесса просачивания в трещиновато-пористый грунт в условиях Крайнего Севера

Приводится математическое моделирование комплексной мультифизичной задачи, актуальной для территорий Крайнего Севера и Арктики. Актуальность данной задачи характеризируется важностью процесса просачивания при формировании и оттаивании слоя вечной мерзлоты. Современные прикладные задачи по большей части требуют учета сложных геометрий, а также большого количества разных процессов и их взаимовлияния. Мультифизичная модель состоит из уравнения Ричардса для описания процесса просачивания, модели двойной пористости для описания естественной трещиноватости грунта, задачи Стефана для описания температурного режима грунта в условиях криолитозоны. Вычислительный алгоритм базируется на конечно-элементной аппроксимации по пространству на триангулированных сетках Делоне и использовании схемы расщепления по времени с применением линеаризации с предыдущего временного слоя.