Для описания процесса течения лавовых потоков или роста лавовых куполов мы рассматриваем задачу численного моделирования растекания вязкой неоднородной несжимаемой жидкости под действием гравитационных сил по заданной поверхности. Математическая модель, описывающая такое движение, задается уравнением Навье-Стокса и уравнением неразрывности.

Программные коды для компьютерного моделирования разработаны на основе пакета SPlisHSPlasH, который представляет собой библиотеку с открытым исходным кодом для языка C++. Для моделирования динамики лавы был написан дополнительный код для переменной вязкости.

Подробнее с результатами можно ознакомиться в разделе Публикации. Так же можно скачать исходный код решателя для переменной вязкозти разделе загрузок.