1. Modeling Multi‐Material Structural Patterns in Tectonic Flow With a Discontinuous Galerkin Level Set Method.
- Author
-
Wú, Qǐháng, Lin, Shoufa, and Unger, André
- Subjects
- *
LEVEL set methods , *STRUCTURAL models , *STOKES flow , *GEOLOGICAL cross sections , *GEOLOGICAL mapping , *VISCOUS flow , *ELLIPSOIDS , *MATERIAL point method , *CREEP (Materials) - Abstract
We formulate a numerical framework, in both 2d and 3d, to model the structural patterns emerging from creeping viscous flow typically encountered in long‐term ductile lithospheric deformation by coupling the discontinuous Galerkin level set method with a finite element Stokes‐like flow solver. The level set formulation has the advantage of retaining information on the interface geometry, decreased memory requirement and improved computational efficiency from the two‐way particle‐mesh information transfer compared to particle‐in‐cell methods. Furthermore, our formulation fully exploits the advantages of the finite element method (e.g., the flexibility of mesh geometry and the ease of handling anisotropic materials) by using a unified finite element framework. The novelty of our formulation is the capability to offer a fully dynamic approach for modeling structural patterns resulting from a tectonic flow that is non‐steady and inhomogeneous. The material distribution and the finite deformation patterns predicted from the numerical model can be directly compared with geological map patterns (e.g., lithological distribution at specified depths and on cross‐sections) and field structural analyses (e.g., foliation, lineation and strain patterns), thus offering the possibility of ground‐truthing the modeling results by field evidence. As examples for potential applications of our method, we apply our method to the modeling of a competent inclusion in simple shear flow, as well as a Rayleigh‐Taylor type density overturn. Our models demonstrate good agreement with previous 2d benchmark results and produce 3d lithological and deformation patterns comparable to field observations. Plain Language Summary: A novel modeling method is proposed to simulate the ductile deformation patterns resulting from viscous tectonic flow. Domains of different material properties are efficiently captured by the scalar level set function, which describes the signed distance to the domain interface. The level set field, paired with passive tracer particles, enables the modeling of the spatial distribution of lithological units, regional fold patterns of stratigraphic layers as well as finite deformation patterns (in terms of foliation/lineation patterns, strain history and strain intensity) in 3d block model, horizontal slices at certain depth level and vertical cross‐sections. Simulated cross‐sections and geological maps can be constructed along any direction and at any depth, thus allowing direct comparison between numerical models and field data obtained from geological mapping. Key Points: We implement a method for simulating tectonic flow in both 2d and 3d capable of predicting ductile deformation patternsAreas of different material properties are efficiently captured with the level set method on a stationary meshTracer particles record simulated strain ellipsoids, thus facilitating comparisons between model results and field data [ABSTRACT FROM AUTHOR]
- Published
- 2023
- Full Text
- View/download PDF