A benchmark for two-dimensional numerical subduction modeling using COMSOL Multiphysics®

유수환,; 이창열,

doi:10.14770/jgsk.2018.54.6.683

Cited by 3 publications

(1 citation statement)

References 10 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…In the present study, we investigate the effect of the regularization with the capped porosity on the liquid transport and mass conservation in the two-phase flow model for the Earth's mantle, using the commercial finite element package COMSOL Multiphysics® (COMSOL, hereafter). COMSOL has been used previously in the context of mantle convection and successfully benchmarked (e.g., Lee, 2013;Yu and Lee, 2018;Trim et al, 2021). COMSOL was used to study the liquid transport in the mantle wedge of subduction zones, considering a simplified porous flow model without the incorporation of the effect of matrix compaction (e.g., Wada and Behn, 2015;Lee et al, 2021;Kim et al, 2022).…”

Section: Introductionmentioning

confidence: 99%

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

Lee,

Cerpa,

Han

et al. 2023

Preprint

View full text Add to dashboard Cite

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 the liquid transport on the geodynamical and geochemical evolution of the Earth, the so-called “compaction equations” are more and more incorporated in geodynamical modeling studies. The implementation of these equations requires a regularization technique to handle the porosity singularity in 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 transport and mass conservation have not been quantitatively evaluated. Here, we investigate these effects using a series of 1- and 2-dimensional numerical models 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 capped and uncapped experiments (i.e., allowing negative porosity). 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 mass flux of liquid across the model boundaries, resulting in intrinsic errors in mass conservation even if high mesh resolution is used. Despite the errors in mass balance, however, the general trends of porosity evolution in the capped experiments are similar to those in the uncapped experiments. Hence, the use of the regularization of the compaction equations with the enforced positive porosity is reasonable for modeling fluid and melt transport in a deforming mantle.

show abstract

Section: Introductionmentioning

confidence: 99%

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

Lee,

Cerpa,

Han

et al. 2023

Preprint

View full text Add to dashboard Cite

show abstract

Numerical or analog modeling study

Lee¹

2018

JGSK

View full text Add to dashboard Cite

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

Lee,

Cerpa,

Han

et al. 2024

Solid Earth

View full text Add to dashboard Cite

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.

show abstract

A benchmark for two-dimensional numerical subduction modeling using COMSOL Multiphysics®

Cited by 3 publications

References 10 publications

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

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

Numerical or analog modeling study

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

Contact Info

Product

Resources

About