Implicit Euler simulation of one-dimensional Burridge-Knopoff model of earthquakes with set-valued friction laws

Xiaogang Xiong, Ryo Kikuuwe, Motoji Yamamoto

Research output: Contribution to journalArticle

3 Citations (Scopus)

Abstract

In the simulations of Burridge-Knopoff (BK) model of earthquakes, the friction force laws are important to produce earthquake-like stick-slip behaviors. Some friction force laws are set-valued and the BK model with them can produce consistent results with observed data of earthquakes in some aspects. However, it is cumbersome to simulate the BK model with set-valued laws by conventional explicit integration methods. In the presence of set-valued laws, the explicit integration methods can easily lead to the numerical chattering, violated constraints on the velocity of force laws, and the difficulty of identifying the states of blocks of the BK model. This paper employs an implicit Euler integration method to simulate the BK model with symmetric and asymmetric set-valued laws. This method removes the numerical chattering in the BK model, even in the cases of large time step sizes. It can easily detect the stuck or slipping state of a block element. Comparing to previous results integrated by explicit integration methods in the literature, the results integrated by this implicit method show smoother curves and lower irregularities in the magnitude distribution of events.

Original languageEnglish
Pages (from-to)1039-1057
Number of pages19
JournalAdvances in Computational Mathematics
Volume41
Issue number6
DOIs
Publication statusPublished - Dec 1 2015

Fingerprint

Earthquake
Euler
Friction
Earthquakes
Simulation
Chattering
Model
Stick-slip
Implicit Method
Irregularity
Curve

All Science Journal Classification (ASJC) codes

  • Computational Mathematics
  • Applied Mathematics

Cite this

Implicit Euler simulation of one-dimensional Burridge-Knopoff model of earthquakes with set-valued friction laws. / Xiong, Xiaogang; Kikuuwe, Ryo; Yamamoto, Motoji.

In: Advances in Computational Mathematics, Vol. 41, No. 6, 01.12.2015, p. 1039-1057.

Research output: Contribution to journalArticle

@article{caf89163464947d5a302e37c85ca4eeb,
title = "Implicit Euler simulation of one-dimensional Burridge-Knopoff model of earthquakes with set-valued friction laws",
abstract = "In the simulations of Burridge-Knopoff (BK) model of earthquakes, the friction force laws are important to produce earthquake-like stick-slip behaviors. Some friction force laws are set-valued and the BK model with them can produce consistent results with observed data of earthquakes in some aspects. However, it is cumbersome to simulate the BK model with set-valued laws by conventional explicit integration methods. In the presence of set-valued laws, the explicit integration methods can easily lead to the numerical chattering, violated constraints on the velocity of force laws, and the difficulty of identifying the states of blocks of the BK model. This paper employs an implicit Euler integration method to simulate the BK model with symmetric and asymmetric set-valued laws. This method removes the numerical chattering in the BK model, even in the cases of large time step sizes. It can easily detect the stuck or slipping state of a block element. Comparing to previous results integrated by explicit integration methods in the literature, the results integrated by this implicit method show smoother curves and lower irregularities in the magnitude distribution of events.",
author = "Xiaogang Xiong and Ryo Kikuuwe and Motoji Yamamoto",
year = "2015",
month = "12",
day = "1",
doi = "10.1007/s10444-014-9398-4",
language = "English",
volume = "41",
pages = "1039--1057",
journal = "Advances in Computational Mathematics",
issn = "1019-7168",
publisher = "Springer Netherlands",
number = "6",

}

TY - JOUR

T1 - Implicit Euler simulation of one-dimensional Burridge-Knopoff model of earthquakes with set-valued friction laws

AU - Xiong, Xiaogang

AU - Kikuuwe, Ryo

AU - Yamamoto, Motoji

PY - 2015/12/1

Y1 - 2015/12/1

N2 - In the simulations of Burridge-Knopoff (BK) model of earthquakes, the friction force laws are important to produce earthquake-like stick-slip behaviors. Some friction force laws are set-valued and the BK model with them can produce consistent results with observed data of earthquakes in some aspects. However, it is cumbersome to simulate the BK model with set-valued laws by conventional explicit integration methods. In the presence of set-valued laws, the explicit integration methods can easily lead to the numerical chattering, violated constraints on the velocity of force laws, and the difficulty of identifying the states of blocks of the BK model. This paper employs an implicit Euler integration method to simulate the BK model with symmetric and asymmetric set-valued laws. This method removes the numerical chattering in the BK model, even in the cases of large time step sizes. It can easily detect the stuck or slipping state of a block element. Comparing to previous results integrated by explicit integration methods in the literature, the results integrated by this implicit method show smoother curves and lower irregularities in the magnitude distribution of events.

AB - In the simulations of Burridge-Knopoff (BK) model of earthquakes, the friction force laws are important to produce earthquake-like stick-slip behaviors. Some friction force laws are set-valued and the BK model with them can produce consistent results with observed data of earthquakes in some aspects. However, it is cumbersome to simulate the BK model with set-valued laws by conventional explicit integration methods. In the presence of set-valued laws, the explicit integration methods can easily lead to the numerical chattering, violated constraints on the velocity of force laws, and the difficulty of identifying the states of blocks of the BK model. This paper employs an implicit Euler integration method to simulate the BK model with symmetric and asymmetric set-valued laws. This method removes the numerical chattering in the BK model, even in the cases of large time step sizes. It can easily detect the stuck or slipping state of a block element. Comparing to previous results integrated by explicit integration methods in the literature, the results integrated by this implicit method show smoother curves and lower irregularities in the magnitude distribution of events.

UR - http://www.scopus.com/inward/record.url?scp=84949101670&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84949101670&partnerID=8YFLogxK

U2 - 10.1007/s10444-014-9398-4

DO - 10.1007/s10444-014-9398-4

M3 - Article

AN - SCOPUS:84949101670

VL - 41

SP - 1039

EP - 1057

JO - Advances in Computational Mathematics

JF - Advances in Computational Mathematics

SN - 1019-7168

IS - 6

ER -