Geophysical Finite-Element Simulation Tool (GeoFEST): Algorithms and Validation for Quasistatic Regional Faulted Crust Problems
Geophysical Finite-Element Simulation Tool (GeoFEST): Algorithms and Validation for Quasistatic Regional Faulted Crust Problems
复制标题
地球物理有限元模拟工具 (GeoFEST)：准静态区域断层地壳问题的算法和验证
DOI：
--
复制
发表时间：
2008
期刊：
影响因子：
--
通讯作者：
A. Donnellan
中科院分区：
文献类型：
--
作者：
J. Parker;G. Lyzenga;C. Norton;C. Zuffada;M. Glasscoe;J. Lou;A. Donnellan
GeoFEST（地球物理有限元模拟工具）是一个二维和三维有限元软件包，用于地球物理和其他连续域应用中的固体应力和应变建模。它是NASA的HockeSim项目的特色高性能应用程序之一。该程序的目标是在UNIX系统上编译和运行，并在不同的系统上运行，包括顺序和消息传递并行系统。椭圆型偏微分方程的解是通过有限元基采样得到的，导致稀疏线性系统主要通过共轭梯度迭代求解到公差水平;在序列系统上，还支持线性系统的直接反演的Crout分解。支持的物理模型包括各向同性线性弹性和牛顿和幂律粘弹性，通过隐式准静态时间步进。除了三角形，四边形，四面体和六面体连续体元素，GeoFEST支持分裂节点故障，体力和表面牵引力。该软件和相关的网格细化策略已在各种测试用例上进行了验证，并与解析解进行了严格的比较。这些包括一个盒子形的域与强加的运动在一个表面上，一对走滑断层的跨越安排，和两个社区商定的基准情况下：走滑断层在封闭的盒子，和四分之一域圆形故障问题。该代码的科学应用包括建模的静态和瞬态的同震和震后地球变形，地球对冰川，大气和水文负荷的响应，以及其他涉及地质介质的整体变形的情况。
GeoFEST (Geophysical Finite Element Simulation Tool) is a two- and three-dimensional finite element software package for the modeling of solid stress and strain in geophysical and other continuum domain applications. It is one of the featured high-performance applications of the NASA QuakeSim project. The program is targeted to be compiled and run on UNIX systems, and is running on diverse systems including sequential and message-passing parallel systems. Solution to the elliptical partial differential equations is obtained by finite element basis sampling, resulting in a sparse linear system primarily solved by conjugate gradient iteration to a tolerance level; on sequential systems a Crout factorization for the direct inversion of the linear system is also supported. The physics models supported include isotropic linear elasticity and both Newtonian and power-law viscoelasticity, via implicit quasi-static time stepping. In addition to triangular, quadrilateral, tetrahedral and hexahedral continuum elements, GeoFEST supports split-node faulting, body forces, and surface tractions. This software and related mesh refinement strategies have been validated on a variety of test cases with rigorous comparison to analytical solutions. These include a box-shaped domain with imposed motion on one surface, a pair of strike slip faults in stepover arrangement, and two community-agreed benchmark cases: a strike slip fault in an enclosing box, and a quarter-domain circular fault problem. Scientific applications of the code include the modeling of static and transient co- and post-seismic earth deformation, Earth response to glacial, atmospheric and hydrological loading, and other scenarios involving the bulk deformation of geologic media.