GPU-based implementation of discrete element method for simulation of the geological fault geometry and position. / Lisitsa, Vadim V.; Tcheverda, Vladimir A.; Volianskaia, Victoria V.
In: Supercomputing Frontiers and Innovations, Vol. 5, No. 3, 2018, p. 46-50.Research output: Contribution to journal › Article › peer-review
}
TY - JOUR
T1 - GPU-based implementation of discrete element method for simulation of the geological fault geometry and position
AU - Lisitsa, Vadim V.
AU - Tcheverda, Vladimir A.
AU - Volianskaia, Victoria V.
PY - 2018
Y1 - 2018
N2 - We present an algorithm for numerical simulation of the geological fault formation. The approach is based on the discrete elements method, which allows modeling of the deformations and structural discontinuity of the Upper part of the Earth crust. In the discrete elements method, the medium is represented as an combination of discrete particles which interact as elastic or viscoelastic bodies. Additionally, external potential forces, for example gravitational forces, may be introduced. At each time step the full set of forces acting at each particle is computed, after that the position of the particle is evaluated on the base of Newtonian mechanics. We implement the algorithm using CUDA technology to simulate single statistical realization of the model, whereas MPI is used to parallelize with respect to different statistical realizations. Obtained numerical results show that for low dip angles of the tectonic displacements relatively narrow faults form, whereas high dip angles of the tectonic displacements lead to a wide V-shaped deformation zones.
AB - We present an algorithm for numerical simulation of the geological fault formation. The approach is based on the discrete elements method, which allows modeling of the deformations and structural discontinuity of the Upper part of the Earth crust. In the discrete elements method, the medium is represented as an combination of discrete particles which interact as elastic or viscoelastic bodies. Additionally, external potential forces, for example gravitational forces, may be introduced. At each time step the full set of forces acting at each particle is computed, after that the position of the particle is evaluated on the base of Newtonian mechanics. We implement the algorithm using CUDA technology to simulate single statistical realization of the model, whereas MPI is used to parallelize with respect to different statistical realizations. Obtained numerical results show that for low dip angles of the tectonic displacements relatively narrow faults form, whereas high dip angles of the tectonic displacements lead to a wide V-shaped deformation zones.
KW - CUDA
KW - Discrete element method
KW - Geological faults
KW - Statistical simulation
UR - http://www.scopus.com/inward/record.url?scp=85057583282&partnerID=8YFLogxK
U2 - 10.14529/jsfi180307
DO - 10.14529/jsfi180307
M3 - Article
AN - SCOPUS:85057583282
VL - 5
SP - 46
EP - 50
JO - Supercomputing Frontiers and Innovations
JF - Supercomputing Frontiers and Innovations
SN - 2409-6008
IS - 3
ER -
ID: 25773929