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