• central scheme
• compressible
• magma ascent
• Riemann
• two phase

In this thesis a multidimensional compressible two-phase flow model with two pressures, two velocities and a single temperature is studied. The system of partial differential equations, derived using the theory of thermodynamically compatible systems, is written in conservative form and the hyperbolicity of the system is shown. Furthermore, for a more complete analysis, the characteristic polynomial is computed. The single temperature model is also compared with the classical Baer-Nunziato model, showing that, although they can be written in similar way, they present important differences.

In order to solve the governing equations, several numerical schemes for the discretization in space, based on the finite volume methods, are considered. Several numerical fluxes, the Lax-Friedrichs, the Generalized FORCE, the Kurganov and Tadmor, the Kurganov, Noelle and Petrova and finally a modification of the Lax-Friedrichs numerical fluxes are introduced. For the discretization in time, two approaches for solving the system of partial differential equations, the splitting and the Runge-Kutta approach, are presented. All the numerical schemes are implemented using the open source CFD (Computational Fluid Dynamics) toolbox "OpenFOAM" (Open source Field Operation And Manipulation, developed by OpenCFD Ltd).

After the definition of the physical model and of the numerical methods, the robustness, the correctness and the accuracy of the proposed schemes are investigated performing several monodimensional and multidimensional well-known numerical tests in presence of rarefaction and shock waves. Two shock-bubble interaction experiments have been reproduced numerically, comparing the results with laboratory observations. The results show a really good agreement between the simulation and the laboratory experiments. The physical model appears to be accurate and the numerical methods robust, being able to properly resolve detailed flow features as shock-wave refractions, reflections and diffractions.

Finally, an application of the single temperature model to the volcanological field is presented, proposing a new model for the initial stages of magma ascent in a conduit during explosive eruptions. The magma ascent model is derived from the single temperature model adding a new transport equation for the gas dissolved in the liquid phase and showing through the calculation of the characteristic polynomial that the new equation does not alter the hyperbolicity of the system. Using one of the numerical schemes proposed, the initial phases of an explosive eruption at Soufriére Hills Volcano are simulated, focusing the attention on the effect of disequilibrium processes.

To conclude, in the appendix, the derivation of the single temperature model using the thermodynamically compatible systems is presented and some details regarding the implementation of the numerical schemes proposed in this thesis using the OpenFOAM framework are given.