LaCoDe : A Lagrangian two-dimensional thermo-mechanical code for large-strain compressible visco-elastic geodynamical modeling. / de Montserrat, Albert; Morgan, Jason P.; Hasenclever, Jörg.

In: Tectonophysics, Vol. 767, 228173, 20.09.2019, p. 1-19.

Research output: Contribution to journalArticle

Published

Documents

  • Accepted Manuscript

    Accepted author manuscript, 8 MB, PDF document

    Embargo ends: 8/08/20

    Licence: CC BY-NC-ND Show licence

Abstract

We present LaCoDe (Lagrangian Compressible Deformation), a MATLAB solver for the Stokes equations for compressible non-Newtonian visco-elastic flow in two dimensions, based on a Lagrangian formulation of the Finite Element Method. The incompressible Boussinesq approximation is a widespread assumption in numerical models of lithospheric deformation, thus potentially masking a significant contribution of mechanisms linked to volumetric changes that occur in the asthenospheric mantle and the lithosphere. LaCoDe employs a compressible formulation of the Stokes equations designed to address such volume-changing processes. First, we provide a description of the equations governing the deformation of Earth rocks and detailed overview of the algorithm, its numerical implementation, treatment of the non-linearities rising from the compressible formulation and ineleastic deformation, and the remeshing algorithm that tracks and transfers the physical fields from a highly-distorted mesh to a high-quality one. LaCoDe is then benchmarked by comparing numerical results to analytical solutions for the bending of a thin elastic beam under a constant uniform load, flow around a rigid inclusion, Rayleigh-Taylor instability, stress build-up in a visco-elastic Maxwell body, and Couette flow with viscous heating. The Rayleigh-Taylor instability test is further used to demonstrate the accuracy of the remeshing algorithm. The importance of including volumetric strain for geodynamic processes is illustrated by two numerical experiments: i) volumetric-strain inducing phase changes in amagmatic slow-spreading ridges, and ii) subducting slabs.

Original languageEnglish
Article number228173
Pages (from-to)1-19
Number of pages19
JournalTectonophysics
Volume767
Early online date8 Aug 2019
DOIs
Publication statusPublished - 20 Sep 2019
This open access research output is licenced under a Creative Commons Attribution-NonCommercial-NoDerivs 3.0 Unported License.

ID: 34470760