Abstract. Numerical simulations of volcanic processes play a fundamental role in understanding the dynamics of magma storage, ascent, and eruption. The recent extraordinary progress in computer performance and improvements in numerical modeling techniques allow simulating multiphase systems in mechanical and thermodynamical disequilibrium. Nonetheless, the growing complexity of these simulations requires the development of flexible computational tools that can easily switch between sub-models and solution techniques. In this work we present MagmaFOAM, a library based on the open-source computational fluid dynamics software OpenFOAM that incorporates models for solving the dynamics of multiphase, multicomponent magmatic systems. Retaining the modular structure of OpenFOAM, MagmaFOAM allows runtime selection of the solution technique depending on the physics of the specific process and sets a solid framework for in-house and community model development, testing, and comparison. MagmaFOAM models thermomechanical nonequilibrium phase coupling and phase change, and it implements state-of-the-art multiple volatile saturation models and constitutive equations with composition-dependent and space–time local computation of thermodynamic and transport properties. Code testing is performed using different multiphase modeling approaches for processes relevant to magmatic systems: Rayleigh–Taylor instability for buoyancy-driven magmatic processes, multiphase shock tube simulations propaedeutical to conduit dynamics studies, and bubble growth and breakage in basaltic melts. Benchmark simulations illustrate the capabilities and potential of MagmaFOAM to account for the variety of nonlinear physical and thermodynamical processes characterizing the dynamics of volcanic systems.

No abstract

The paper deals with the mathematical modelling of a gas slug rising in a cylindrical duct filled with an incompressible liquid. This research is motivated by a phenomenon commonly observed during Strombolian eruptions at basaltic volcanoes, that is, mildly explosive events driven by a large bubble of magmatic gas (a slug) rising up the conduit and bursting at the surface. The model is compared with the laboratory experiments described in [15] and we prove that the constancy of the slug ascent velocity observed in these experiments is fully justifiable. The model is developed both in the case of a Newtonian fluid and in the case of a non-Newtonian power law fluid, more suitable for magmas, and applied to a volcanic scenario.

Abstract. Numerical simulations of volcanic processes play a fundamental role in understanding the dynamics of magma storage, ascent and eruption. The recent extraordinary progress in computer performance and improvements in numerical modeling techniques allow simulating multiphase systems in mechanical and thermodynamical disequilibrium. Nonetheless, the growing complexity of these simulations requires the development of flexible computational tools that can easily switch between sub-models and solution techniques. In this work we present MagmaFOAM, a library based on the open source computational fluid dynamics software OpenFOAM, that incorporates models for solving the dynamics of multiphase, multicomponent magmatic systems. Retaining the modular structure of OpenFOAM, MagmaFOAM allows run-time selection of the solution technique depending on the physics of the specific process, and sets a solid framework for in-house and community model development, testing and comparison. MagmaFOAM models thermo-mechanical non-equilibrium phase coupling and phase change, and implements state-of-the-art multiple volatile saturation models and constitutive equations with composition-dependent and space-time local computation of thermodynamic and transport properties. Code testing is performed using different multiphase modeling approches for processes relevant to magmatic systems: Rayleigh-Taylor instability, for buyoancy-driven magmatic processes; multiphase shock tube simulations, propedeutical to conduit dynamics studies; bubble growth and breakage in basaltic melts. Benchmark simulations illustrate the capabilities and potential of MagmaFOAM to account for the variety of non-linear physical and thermodynamical processes characterizing the dynamics of volcanic systems.