HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS		QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

This Article Figures Only Full Text Full Text (PDF) Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in Web of Science Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Web of Science (6) Citing Articles via Google Scholar Google Scholar Articles by Dalguer, L. A. Articles by Day, S. M. Search for Related Content GeoRef GeoRef Citation

Comparison of Fault Representation Methods in Finite Difference Simulations of Dynamic Rupture

¹ Department of Geological Sciences
San Diego State University
San Diego, California 92182

Assessing accuracy of numerical methods for spontaneous rupture simulation is challenging because we lack analytical solutions for reference. Previous comparison of a boundary integral method (BI) and finite-difference method (called DFM) that explicitly incorporates the fault discontinuity at velocity nodes (traction- at-split-node scheme) shows that both converge to a common, grid-independent solution and exhibit nearly identical power-law convergence rates with respect to grid spacing x. We use this solution as a reference for assessing two other proposed finite-difference methods, the thick fault (TF) and stress glut (SG) methods, both of which approximate the fault-jump conditions through inelastic increments to the stress components (inelastic-zone schemes). The TF solution fails to match the qualitative rupture behavior of the reference solution and has quantitative misfits in root- mean-square rupture time of 30% for the smallest computationally feasible x (with 9 grid-point resolution of cohesive zone, denoted _c = 9). For sufficiently small values of x, the SG method reproduces the qualitative features of the reference solution, but rupture velocity remains systematically low for SG relative to the reference solution, and SG lacks the well-defined power-law convergence seen for BI and DFM. The rupture-time error for SG, with _c 9, remains well above uncertainty in the reference solution, and the split-node method attains comparable accuracy with _c 1/4 as large (and computation timescales as (_c)⁴). Thus, accuracy is highly sensitive to the formulation of the fault-jump conditions: The split-node method attains power-law convergence. The SG inelastic-zone method achieves solutions that are qualitatively meaningful and quantitatively reliable to within a few percent, but convergence is uncertain, and SG is computationally inefficient relative to the split-node approach. The TF inelastic-zone method does not achieve qualitatively meaningful solutions to the 3D test problem and is sufficiently computationally inefficient that it is not feasible to explore convergence quantitatively.

This article has been cited by other articles:

				L. A. Dalguer, H. Miyake, S. M. Day, and K. Irikura Surface Rupturing and Buried Dynamic-Rupture Models Calibrated with Statistical Observations of Past Earthquakes Bulletin of the Seismological Society of America, June 1, 2008; 98(3): 1147 - 1161. [Abstract] [Full Text] [PDF]

				K. B. Olsen, S. M. Day, J. B. Minster, Y. Cui, A. Chourasia, D. Okaya, P. Maechling, and T. Jordan TeraShake2: Spontaneous Rupture Simulations of Mw 7.7 Earthquakes on the Southern San Andreas Fault Bulletin of the Seismological Society of America, June 1, 2008; 98(3): 1162 - 1185. [Abstract] [Full Text] [PDF]

				S. M. Corish, C. R. Bradley, and K. B. Olsen Assessment of a Nonlinear Dynamic Rupture Inversion Technique Applied to a Synthetic Earthquake Bulletin of the Seismological Society of America, June 1, 2007; 97(3): 901 - 914. [Abstract] [Full Text] [PDF]