Звучит как типичный кошмар вычислительной гидродинамики – диффузия тепла, конвекция, фазовый переход, турбулентность… Начнут с уравнения теплопроводности, добавят источник для магмы, потом сведут это с решателем уравнений Навье-Стокса для потока жидкости. Скорее всего, понадобится нерегулярная сетка или конечно-элементная модель, а для давления – метод прогонки. В коде – тонна матричных операций, куча временных шагов и масса моментов «почему это нестабильно?». Настоящее удовольствие – когда температура взрывается, градиент давления запускает извержение, и ты можешь увидеть, как симуляция взорвётся в реальном времени – только не забудь подготовить профилировщик, а то будет «вылет».