A numerical code for the simulation of the dynamics of multicomponent fluids based on the highly stable and accurate finite element algorithm by Hauke and Hughes (1998) has been implemented. This numerical method allows a unified approach forcompressible and incompressible transient flows. Therefore, it is particularly suitable for the simulation of the dynamics within magma chambers and along volcanic conduits, where a wide range of Mach and Reynolds numbers occurs. The balance equations of mass, momentum, energy, and composition are solved for the unknowns pressure, velocity, temperature, and composition of a homogeneous mixture with properties dependent on the local conditions. The equations are discretized in time and space with Galerkin least-squares and discontinuity-capturing stabilizing techniques. The conservation equations for chemical components have been added to the original Hauke and Hughes (1998) formulation, along with the corresponding stabilization terms. The linear non-symmetric system of discretized equations is solved with a preconditioned GMRES. The code is written in C++, picking up FE and mathematical tools from the open source OFELI, Diffpack, and MTL libraries. The computational results have been validated on classical test cases in a wide range of flow conditions from compressible to incompressible. Applications to magma chamber and conduit flow dynamic problems show several features of the multidimensional transient dynamics before, during, and after volcanic eruptions.
A numerical code for the simulation of transient multicomponent magma dynamics
BARSANTI, MICHELE;
2004-01-01
Abstract
A numerical code for the simulation of the dynamics of multicomponent fluids based on the highly stable and accurate finite element algorithm by Hauke and Hughes (1998) has been implemented. This numerical method allows a unified approach forcompressible and incompressible transient flows. Therefore, it is particularly suitable for the simulation of the dynamics within magma chambers and along volcanic conduits, where a wide range of Mach and Reynolds numbers occurs. The balance equations of mass, momentum, energy, and composition are solved for the unknowns pressure, velocity, temperature, and composition of a homogeneous mixture with properties dependent on the local conditions. The equations are discretized in time and space with Galerkin least-squares and discontinuity-capturing stabilizing techniques. The conservation equations for chemical components have been added to the original Hauke and Hughes (1998) formulation, along with the corresponding stabilization terms. The linear non-symmetric system of discretized equations is solved with a preconditioned GMRES. The code is written in C++, picking up FE and mathematical tools from the open source OFELI, Diffpack, and MTL libraries. The computational results have been validated on classical test cases in a wide range of flow conditions from compressible to incompressible. Applications to magma chamber and conduit flow dynamic problems show several features of the multidimensional transient dynamics before, during, and after volcanic eruptions.I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.
