Numerical simulation of magma-rock interaction at Krafla volcano using OpenFOAM software and a simplified thermal model

Gerbault, Muriel; Melnik, Oleg; Borisova, Anastassia

doi:10.5194/se-16-1307-2025

Articles | Volume 16, issue 11

https://doi.org/10.5194/se-16-1307-2025

© Author(s) 2025. This work is distributed under
the Creative Commons Attribution 4.0 License.

https://doi.org/10.5194/se-16-1307-2025

© Author(s) 2025. This work is distributed under
the Creative Commons Attribution 4.0 License.

Articles | Volume 16, issue 11

Research article

|

12 Nov 2025

Research article |

| 12 Nov 2025

Numerical simulation of magma-rock interaction at Krafla volcano using OpenFOAM software and a simplified thermal model

Muriel Gerbault, Oleg Melnik, and Anastassia Borisova

Data sets

Solver and run file for the paper "Numerical simulation of magma-rock interaction at Krafla volcano using OpenFOAM software and a simplified thermal model" Muriel Gerboult and Oleg Melnik https://doi.org/10.5281/zenodo.17426563

Model code and software

Solver and run file for the paper "Numerical simulation of magma-rock interaction at Krafla volcano using OpenFOAM software and a simplified thermal model" Muriel Gerboult and Oleg Melnik https://doi.org/10.5281/zenodo.17426563

Video supplement

Numerical simulation of magma-rock interaction at Krafla volcano using OpenFOAM software and a simplified thermal model Oleg Melnik and Muriel Gerboult https://doi.org/10.5446/71779

Numerical simulation of magma-rock interaction at Krafla volcano using OpenFOAM software and a simplified thermal model Oleg Melnik and Muriel Gerboult https://doi.org/10.5446/71778

Short summary

We present a 2D numerical modelling of magma-rock interaction after a sill injection into shallow crust. A solver based on OpenFOAM package simulates heat and momentum transfer between magmas of different compositions, with contrasting properties. Two scenarios are considered to explain the presence of fresh rhyolitic fragments excavated during drilling at Krafla caldera in 2009: partial melting of felsic crust triggered by either a 300 m thick rhyolite intrusion or a 100 m thick basalt sill.