Stable, High Order Accurate Adaptive Schemes for Long Time, Highly Intermittent Geophysics Problems
2013 (English)Report (Other academic)
Many geophysical phenomena are characterized by properties that evolve over a wide range of scales which introduce difficulties when attempting to model these features in one computational method. We have developed a high-order finite difference method for the elastic wave equation that is able to efficiently handle varying temporal scales in a single, stand-alone framework. We apply this method to earthquake cycle models characterized by extremely long interseismic periods interspersed with abrupt, short periods of dynamic rupture. Through the use of summation-by-parts operators and weak enforcement of boundary conditions we derive a provably stable discretization. Time stepping is achieved through the implicit θ-method which allows us to take large time steps during the intermittent period between earthquakes and adapts appropriately to fully resolve rupture.
Place, publisher, year, edition, pages
Linköping: Linköping University Electronic Press, 2013. , 26 p.
LiTH-MAT-R, ISSN 0348-2960 ; 10
high-order accuracy, stability, adaptive time-integration, summation-by-parts, weak boundary condition, earthquake cycle
IdentifiersURN: urn:nbn:se:liu:diva-98205ISRN: LiTH-MAT-R--2013/10--SEOAI: oai:DiVA.org:liu-98205DiVA: diva2:652989