Abstract:Fluid and melt transport in the solid mantle can be modeled as a two-phase flow in which the liquid flow is resisted by the compaction of the viscously deforming solid mantle. Given the wide impact of liquid transport on the geodynamical and geochemical evolution of the Earth, the so-called "compaction equations" are increasingly being incorporated into geodynamical modeling studies. When implementing these equations, it is common to use a regularization technique to handle the porosity singularity in the dry mantle. Moreover, it is also common to enforce a positive porosity (liquid fraction) to avoid unphysical negative values of porosity. However, the effects of this "capped" porosity on the liquid flow and mass conservation have not been quantitatively evaluated. Here, we investigate these effects using a series of 1- and 2-dimensional numerical models implemented using the commercial finite-element package COMSOL Multiphysics ® . The results of benchmarking experiments against a semi-analytical solution for 1- and 2-D solitary waves illustrate the successful implementation of the compaction equations. We show that the solutions are accurate when the element size is smaller than half of the compaction length. Furthermore, in time-evolving experiments where the solid is stationary (immobile), we show that the mass balance errors are similarly low for both the capped and uncapped (i.e., allowing negative porosity) experiments. When Couette flow, convective flow, or subduction corner flow of the solid mantle is assumed, the capped porosity leads to overestimations of the mass of liquid in the model domain and the mass flux of liquid across the model boundaries, resulting in intrinsic errors in mass conservation even if a high mesh resolution is used. Despite the errors in mass balance, however, the distributions of the positive porosity and peaks (largest positive liquid fractions) in both the uncapped and capped experiments are similar. Hence, the capping of porosity in the compaction equations can be reasonably used to assess the main pathways and first-order distribution of fluids and melts in the mantle.

Density-Driven Compactional Flow in Porous Media

Nonlinear Viscoelastic Compaction in Sedimentary Basins

A Mathematical Model for Voigt Poro-Visco-Plastic Deformation

Constitutive Model for Granular-Fluid Flows Based on Stress Decomposition

A Unified Approach to Mechanical Compaction, Pressure Solution, Mineral Reactions and the Temperature Distribution in Hydrocarbon Basins

Modeling flow in porous media with double porosity/permeability: Mathematical model, properties, and analytical solutions

Combinatorial Model of Solute Transport in Porous Media.

Elastic–viscoplastic constitutive theory of dense granular flow and its three-dimensional numerical realization

Theoretical modeling of squeezing flow in porous media under arbitrary boundary velocity

Influence of the compressibility of pore fluid on the elasto-plastic analysis of non-saturated porous media

Modeling pore fluid pressure and deformation in unsaturated soils under gravitational compaction and cyclic loading

Modeling flow in porous media with double porosity/permeability: A stabilized mixed formulation, error analysis, and numerical solutions

Fluid flow through anisotropic and deformable double porosity media with ultra-low matrix permeability: A continuum framework

Pressure Solution in Sedimentary Basins: Effect of Temperature Gradient

Modeling liquid transport in the Earth's mantle as two-phase flow: effect of an enforced positive porosity on liquid flow and mass conservation

Fluid flow channeling and mass transport with discontinuous porosity distribution

A Density-Dependent Constitutive Model of Rockfill as Well as a Dynamic Contact Technique for Simulation of Dynamic Compaction with MPM

A Unified Constitutive Model for Pressure Sensitive Shear Flow Transitions in Moderate Dense Granular Materials.

Flow Through Fissured Media with Dual Porosity and Nonlinear Elastic Behavior

Investigation of Inertial Effect in Simplified Porous Media Flow

A fully implicit and consistent finite element framework for modeling reservoir compaction with large deformation and nonlinear flow model. Part I: theory and formulation