A Barnes-Hut scheme for simulating fault slip
| Autor: | N. M. Beeler, T. E. Tullis |
|---|---|
| Jazyk: | angličtina |
| Rok vydání: | 2011 |
| Předmět: | |
| Zdroj: | Nonlinear Processes in Geophysics, Vol 18, Iss 2, Pp 133-146 (2011) |
| Druh dokumentu: | article |
| ISSN: | 1023-5809 1607-7946 |
| DOI: | 10.5194/npg-18-133-2011 |
| Popis: | To account for natural spatial and temporal complexity, large-scale, long-duration calculations are required for simulations of seismicity in fault zones that host large earthquakes. Without advances in computational methods, the rate of progress in "earthquake simulator" models and associated earthquake forecasts is limited by the rates at which computer speed and storage increase. To explore improvements in computational efficiency we develop the first implementation of the Barnes-Hut algorithm (Barnes and Hut, 1986) to calculate elastic interactions in a fault model. The Barnes-Hut method is an efficient, numerical scheme that treats local forces exactly and distant forces approximately. The approach is illustrated in example simulations of non-linear fault strength in plane strain. Rudimentary error analysis indicates that efficient calculations, where execution time scales with number of grid points (N) as N log N, can be conducted routinely with errors on the order of 0.1%. We expect the Barnes-Hut method to be well suited for conducting initial exploration of parameter space for fault simulations with non-linear constitutive equations, and for efficient calculations of stress interaction in complex fault systems. |
| Databáze: | Directory of Open Access Journals |
| Externí odkaz: |
|