High order accurate and explicit time-stable solvers are well suited for hyperbolic wave propagation problems. As a result of the complexities of real geometries, internal interfaces and nonlinear boundary and interface conditions, discontinuities and sharp wave fronts may become fundamental features of the solution. Thus, geometrically flexible and adaptive numerical algorithms are critical for high fidelity and efficient simulations of wave phenomena in many applications. Adaptive curvilinear meshes hold promise to minimise the effort to represent complicated geometries or heterogeneous material data avoiding the bottleneck of feature-preserving meshing. To enable the design of stable DG methods on three space dimensional (3D) curvilinear elements we construct a structure preserving anti-symmetric coordinate transformation motivated by the underlying physics. Using a physics-based numerical penalty-flux, we develop a 3D provably energy-stable discontinuous Galerkin finite element approximation of the elastic wave equation in geometrically complex and heterogenous media. By construction, our numerical flux is upwind and yields a discrete energy estimate analogous to the continuous energy estimate. The ability to treat conforming and non-conforming curvilinear elements allows for flexible adaptive mesh refinement strategies. The numerical scheme has been implemented in ExaHyPE, a simulation engine for parallel dynamically adaptive simulations of wave problems on adaptive Cartesian meshes. We present 3D numerical experiments of wave propagation in heterogeneous isotropic and anisotropic elastic solids demonstrating stability and high order accuracy. We demonstrate the potential of our approach for computational seismology in a regional wave propagation scenario in a geologically constrained 3D model including the geometrically complex free-surface topography of Mount Zugspitze, Germany.
翻译:由于真实的地貌、内部界面和非线性边界及界面条件的复杂性,不连贯和尖锐的波面可能会成为解决方案的基本特征。因此,几何灵活和适应性的数字算法对于在许多应用中高忠诚度和对波现象的有效模拟至关重要。适应性的卷尾线网舍有望最大限度地减少代表复杂的地貌或各种材料数据的努力,避免地貌保存网状的瓶颈。为了能够在三种空间的深度(3D)上设计稳定的DG方法,我们建造了一个结构,以保持由基础物理学驱动的反对称协调转变。我们使用基于物理的数字惩罚峰值,开发出一个3D可辨别不固定的温度加勒金定值元素,在对地貌复杂和变异性媒体中的弹性波状模型。通过构建,我们的数字流是上方,并产生一个类似于连续的DGD方法,在三个空间的深度(3D)上层线线系中,我们构建了一个稳定度的轨道结构,在不断的轨迹上,我们可以将一个不精确的轨道上变动的变动的变动的变动的机图。我们用来处理一个不动的模型的变动的变动的变动的模型。