EGU General Assembly 2021
© Author(s) 2021. This work is distributed under
the Creative Commons Attribution 4.0 License.

Poroelastoplastic modeling of borehole shear bands on high order curvilinear meshes using CUDA technology

Anatoly Vershinin1,2, Vladimir Levin1, and Yury Podladchikov1,3
Anatoly Vershinin et al.
  • 1Lomonosov Moscow State University, mechanical and mathematical, Moscow, Russian Federation (
  • 2IPE RAS, Moscow, Russia (
  • 3University of Lausanne, Lausanne, Switzerland (

The presentation describes an approach to solving problems of modeling the development of zones of localization of plastic deformations within the framework of a poroelastoplastic model generalizing Biot's model. A distinctive feature of this model is a two-way coupling between mechanical processes occurring in a porous elastoplastic matrix and a saturating viscous fluid.

  For the numerical solution of the problem, a variational formulation based on the Galerkin method and the isoparametric  spectral element method (SEM) is used to discretize the geometric model and PDEs on curvilinear unstructured SEM meshes. SEM orders up to the 15th were used for calculations.

  The software implementation of the developed algorithm based on SEM is performed using CUDA. A spectral element mesh is naturally mapped to a CUDA grid of SMs, and accordingly, each spectral element is mapped to a streaming block, within which individual nodes are processed by the corresponding threads within the block.

The research for this article is performed partially in Schmidt Institute of Physics of the Earth of the Russian Academy of Sciences and supported by the Russian Science Foundation under grant № 19-77-10062.

How to cite: Vershinin, A., Levin, V., and Podladchikov, Y.: Poroelastoplastic modeling of borehole shear bands on high order curvilinear meshes using CUDA technology, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-15600,, 2021.

Display materials

Display file