Please use this identifier to cite or link to this item:
|Title:||Gamra : simple meshing for complex earthquakes||Authors:||Landry, Walter
Immersed Interface Method
|Issue Date:||2016||Source:||Landry, W., & Barbot, S. (2016). Gamra: Simple meshing for complex earthquakes. Computers & Geosciences, 90, 49-63. doi:10.1016/j.cageo.2016.02.014||Series/Report no.:||Computers & Geosciences||Abstract:||The static offsets caused by earthquakes are well described by elastostatic models with a discontinuity in the displacement along the fault. A traditional approach to model this discontinuity is to align the numerical mesh with the fault and solve the equations using finite elements. However, this distorted mesh can be difficult to generate and update. We present a new numerical method, inspired by the Immersed Interface Method (Leveque and Li, 1994), for solving the elastostatic equations with embedded discontinuities. This method has been carefully designed so that it can be used on parallel machines on an adapted finite difference grid. We have implemented this method in Gamra, a new code for earth modeling. We demonstrate the correctness of the method with analytic tests, and we demonstrate its practical performance by solving a realistic earthquake model to extremely high precision.||URI:||https://hdl.handle.net/10356/105359
|ISSN:||0098-3004||DOI:||10.1016/j.cageo.2016.02.014||Rights:||© 2016 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).||Fulltext Permission:||open||Fulltext Availability:||With Fulltext|
|Appears in Collections:||EOS Journal Articles|
Items in DR-NTU are protected by copyright, with all rights reserved, unless otherwise indicated.