: Supercomputing: 5th Russian Supercomputing Days, RuSCDays 2019 (Moscow, Russia, September 23-24, 2019), Revised Selected Papers. Communications in Computer and Information Science, Vol. 1129
In this paper, we present an algorithm for numerical simulation of the tectonic movements leading to the formation of geological faults. We use the discrete element method, so that the media are presented as an agglomeration of elastic, visco-elastic, or elasto-plastic interacting particles. This approach can naturally handle finite deformations and can account for the structural discontinuities is the Earth crust. 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.