Stochastic simulation algorithms for iterative solution of the LamГЁ equation
I.A. Aksyuk, A.E. Kireeva, K.K. Sabelfeld, D.D. Smirnov
Institute of Computational Mathematics and Mathematical Geophysics of Siberian Branch of Russian Academy of Sciences, Novosibirsk, Russia
Keywords: meshless stochastic algorithm, random walk on spheres, global random walk algorithm, randomized algorithm for solving linear equations
Abstract
In this paper, iterative stochastic simulation algorithms for the Lamè equation describing the displacements of an isotropic elastic body are constructed. Three different stochastic methods are proposed: the first one is based on a global algorithm of random walk on spheres to compute the solution and its derivatives for an anisotropic diffusion equation. It does not use grids and does not require large amounts of RAM. The second method is based on a randomized algorithm for solving large systems of linear equations and requires the introduction of a grid. The third method is also grid-based and uses a random walk algorithm. All three methods implement an iterative process, at each step of which anisotropic diffusion equations are solved. The paper provides a comparative analysis of the proposed methods and discusses the limits of applicability of each of them.
|