DynEarthSol2D: An efficient unstructured finite element method to study long-term tectonic deformation
DynEarthSol2D: An efficient unstructured finite element method to study long-term tectonic deformation
复制标题
DOI：
10.1002/jgrb.50148
复制
发表时间：
2013-05-01
影响因子：
3.9
通讯作者：
Calo, V. M.
中科院分区：
文献类型：
--
作者：
Choi, E.;Tan, E.;Calo, V. M.
许多构造问题需要将岩石圈视为可压缩的弹性材料，根据施加的应力水平和温度条件，该材料也可以粘性流动或以脆性方式破裂。我们提出了一种灵活的方法来解决由此产生的复杂材料响应，这对所使用的离散化和流变模型提出了严峻的挑战。这种鲁棒、自适应、二维有限元方法使用非结构化网格求解拉格朗日形式的动量平衡和热方程。该方法的实现随本文的发布而向公众发布，并命名为 DynEarthSol2D（可在 http://bitbucket.org/tan2/dynearthsol2 获取）。求解器在剪切应变集中（局部化）的地方使用条件网格自适应性，并在标记粒子的辅助下使用保守映射，以在重新网格化期间保留相和相边界。我们详细介绍了求解器，并根据文献中的解析和数值解决方案在许多基准问题中对其进行了验证。这些结果使我们能够验证和验证我们的软件框架，并显示与康体佳算法的快速拉格朗日分析的早期实现相比，其性能提高了一个数量级。
Many tectonic problems require to treat the lithosphere as a compressible elastic material, which can also flow viscously or break in a brittle fashion depending on the stress level applied and the temperature conditions. We present a flexible methodology to address the resulting complex material response, which imposes severe challenges on the discretization and rheological models used. This robust, adaptive, two-dimensional, finite element method solves the momentum balance and the heat equation in Lagrangian form using unstructured meshes. An implementation of this methodology is released to the public with the publication of this paper and is named DynEarthSol2D (available at http://bitbucket.org/tan2/dynearthsol2). The solver uses contingent mesh adaptivity in places where shear strain is focused (localization) and a conservative mapping assisted by marker particles to preserve phase and facies boundaries during remeshing. We detail the solver and verify it in a number of benchmark problems against analytic and numerical solutions from the literature. These results allow us to verify and validate our software framework and show its improved performance by an order of magnitude compared against an earlier implementation of the Fast Lagrangian Analysis of Continua algorithm.