Felsic melt migration via porous flow – a numerical modeling approach

2021 ◽  
Author(s):  
Petra Maierová ◽  
Pavlína Hasalová ◽  
Karel Schulmann
Keyword(s):  
Numerical Modeling ◽  
Solid Rock ◽  
Driving Forces ◽  
Solid Matrix ◽  
Unusual Behavior ◽  
Two Phase ◽  
Porous Flow ◽  
Melt Migration ◽  
Metasedimentary Rocks ◽  
Melt Transport

<p>Melting of the continental crust and subsequent melt transport has been most thoroughly described in the case of metasedimentary rocks. In these rocks segregation and migration of melt occur either through an interconnected network of veins and melt-rich layers (leucosome) or in form of diapirs. For these rocks, porous flow of melt at grain scale is mostly regarded only as a transient stage of separation of melt from the solid rock.</p><p>An entirely different style of melting and melt transport occurs in the case of felsic metaigneous rocks. We use the example from the Bohemian Massif, the eastern European Variscan belt, where metaigneous migmatites were studied in large detail. Here, melt did not segregate from the solid rock but migrated pervasively along most of the grain boundaries and equilibrated with the host rock. This equilibration resulted in formation of a continuous sequence of texturally, geochemically and compositionally different migmatites.</p><p>The question arises, what are the conditions and driving forces for this unusual behavior. We attempt to address this question by means of numerical modeling of two-phase flow (i.e. flow of porous solid matrix and melt), using the open-source finite-element ASPECT code (aspect.geodynamics.org). Most previous numerical studies of this process were either purely generic or focused on the melting of the mantle. In order to study this process in crustal conditions, we set up a 2D crustal-scale thermo-mechanical model that includes melting and freezing. We investigate the role of material properties (viscosity, solidus and liquidus temperatures, solid matrix permeability, melt composition) and thermal and velocity boundary conditions, as well as the effect of grid resolution. The results are discussed in terms of realistic parameter values and possible styles of melt migration and deformation of the matrix.</p>

The Numerical Modeling of Spherical Explosion in the Foam

2016 ◽  
Vol 11 (1) ◽  
pp. 60-65 ◽  
Author(s):  
R.Kh. Bolotnova ◽  
E.F. Gainullina
Keyword(s):  
Shock Wave ◽  
Numerical Modeling ◽  
Initial Pressure ◽  
Water Volume ◽  
Symmetric Model ◽  
Two Phase ◽  
Experimental Conditions ◽  
Initial Water ◽  
Aqueous Foam ◽  
Spherical Explosion

The spherical explosion propagation process in aqueous foam with the initial water volume content α10=0.0083 corresponding to the experimental conditions is analyzed numerically. The solution method is based on the one-dimensional two-temperature spherically symmetric model for two-phase gas-liquid mixture. The numerical simulation is built by the shock capturing method and movable Lagrangian grids. The amplitude and the width of the initial pressure pulse are found from the amount of experimental explosive energy. The numerical modeling results are compared to the real experiment. It’s shown, that the foam compression in the shock wave leads to the significant decrease in velocity and in amplitude of the shock wave.

Numerical Modeling of Evolution of Boundary Between Incompressible Gas and Liquid in Two-Dimensional Region

2006 ◽  
Vol 4 ◽  
pp. 224-236
Author(s):  
A.S. Topolnikov
Keyword(s):  
Numerical Modeling ◽  
Interphase Boundary ◽  
Incompressible Liquid ◽  
Stokes Equations ◽  
Numerical Code ◽  
Navier Stokes ◽  
Two Phase ◽  
Navier Stokes Equations ◽  
Incompressible Media ◽  
Good Agreement

The paper is devoted to numerical modeling of Navier–Stokes equations for incompressible media in the case, when there exist gas and liquid inside the rectangular calculation region, which are separated by interphase boundary. The set of equations for incompressible liquid accounting for viscous, gravitational and surface (capillary) forces is solved by finite-difference scheme on the spaced grid, for description of interphase boundary the ideology of Level Set Method is used. By developed numerical code the set of hydrodynamic problems is solved, which describe the motion of two-phase incompressible media with interphase boundary. As a result of numerical simulation the solutions are obtained, which are in good agreement with existing analytical and experimental solutions.

Numerical modeling of two-phase binary fluid mixing using mixed finite elements

2012 ◽  
Vol 16 (4) ◽  
pp. 1101-1124 ◽  
Author(s):  
Shuyu Sun ◽  
Abbas Firoozabadi ◽  
Jisheng Kou
Keyword(s):  
Numerical Modeling ◽  
Finite Elements ◽  
Mixed Finite Elements ◽  
Fluid Mixing ◽  
Two Phase ◽  
Binary Fluid

Melting and melt transport mechanisms in the dynamic earth: from melt segregation, extraction to the formation of crustal magmatic systems

2021 ◽  
Author(s):  
Harro Schmeling
Keyword(s):  
Continental Crust ◽  
Dynamic Models ◽  
Elevated Temperatures ◽  
Two Phase Flow ◽  
Solidus Temperature ◽  
Phase Flow ◽  
Transport Mechanisms ◽  
Two Phase ◽  
Melt Transport ◽  
Open Question

<p><strong>Introduction</strong></p><p>At various regions within the dynamic earth melts are generated due to decompressional melting, reduction of the solidus temperature due to volatiles or due to elevated temperatures. They segregate from these partially molten regions, rise by various transport mechanisms and may form crustal magmatic systems where they are emplaced or erupt. The physics of various aspects of this magmatic cycle will be addressed.</p><p><strong>Melt transport mechanisms</strong></p><p>Starting from a partially molten region by which mechanism(s) does the melt segregate out of the melt source region and rise through the mantle or crust? The basic mechanism is two-phase flow, i.e. a liquid phase percolates through a solid, viscously deforming matrix. The corresponding equations and related issues such as compaction or effective matrix rheology are addressed. Beside simple Darcy flow, special solutions of the equations are addressed such as solitary porosity waves. Depending on the bulk to shear viscosity ratio of the matrix and the non-dimensional size of these waves, they show a variety of features: they may transport melt over large distances, or they show transitions from rising porosity waves to diapiric rise or to fingering. Other solutions of the equations lead to channeling, either mechanically or chemically driven. One open question is how do such channels transform into dykes which have the potential of rising through sub-solidus overburden. A recent hypothesis addresses the possibility that rapid melt percolation may reach the thermal non-equilibrium regime, i.e. the local temperature of matrix and melt may evolve differently.  Once dykes have been formed they may propagate upwards driven by melt buoyancy and controlled by the ambient stress field. Often in dynamic models the complexities of melt transport are simplified by parameterized melt extraction. The limitations of such simplifications will be addressed.</p><p><strong>Modelling magmatic systems in thickened continental crust </strong></p><p>Once basaltic melts rise from the mantle, they may underplate continental crust and generate silicic melts. Early dynamic models (Bittner and Schmeling, 1995, Geophys. J. Int.) showed that such silicic magma bodies may rise to mid-crustal depth by diapirism. More recent approaches (e.g. Blundy and Annan, 2016, Elements) emplace sill intrusions into the crust at various levels and calculate the thermal and melting effects responsible for the formation of mush zones. Recently Schmeling et al. (2019, Geophys. J. Int.) self-consistently modelled the formation of crustal magmatic systems, mush zones and magma bodies by including two-phase flow, melting/solidification and effective power-law rheology. In these models melt is found to rise to mid-crustal depths by a combination of compaction/decompaction assisted two-phase flow, sometimes including solitary porosity waves, diapirism or fingering. An open question in these models is whether or how dykes may self-consistently form to transport the melts to shallower depth. First models which combine elastic dyke-propagation (Maccaferri et al., 2019, G-cubed) with the two-phase flow crustal models are promising.</p>

Experimental and Numerical Modeling of Gaseous Carbon Dioxide Injecting Into Aqueous Water in Bottle Filler System

2001 ◽  
Author(s):  
Y. H. Zheng ◽  
R. S. Amano
Keyword(s):  
Carbon Dioxide ◽  
Numerical Modeling ◽  
Tap Water ◽  
Co2 Injection ◽  
Operating Characteristics ◽  
Two Phase ◽  
Carbonation Process ◽  
Filling System ◽  
Gas System ◽  
Gaseous Carbon Dioxide

Abstract An efficient enhancement of the carbonation rate in the bottle filling stage can substantially increase the production in beverage industries. The bottle filling system currently used in most of the manufacturers can still be improved for a better performance of carbonation by designing the injection tube system. This paper reports on an experimental and numerical mass transfer modeling that can simulate the dissolution process of gaseous carbon dioxide into aqueous water in the bottle filler system. In order to establish the operating characteristics of the bottle filler system, an ordinary tap water and pure carbon dioxide were used as the liquid-gas system. The two-phase numerical modeling was developed that can serve as a framework for the continuous improvement of the design of the carbonation process in the bottle filler system. For an optimal design of CO2 injection tube and flow conditions, a computational fluid dynamics (CFD) approach is one of the most power tools. However, since only limited experimental data are available in the open literature to verify the computational results, an experiment study was performed to obtain measurements of CO2 level, temperature, and pressure during the carbonation process in the bottle filled with liquid. Both experimental and numerical studies of various flow condition and different sizes of injection tube are presented in this paper.

Sign in / Sign up

Export Citation Format

Share Document