Abstract:
We consider an application of the Message Passing Interface (MPI) technology for parallelization of the program code which solves equation of the linear elasticity theory. The solution of this equation describes the propagation of elastic waves in demormable rigid bodies. The solution of such direct problem of seismic wave propagation is of interest in seismics and geophysics. Our implementation of solver uses grid-characteristic method to make simulations. We consider technique to reduce time of communication between MPI processes during the simulation. This is important when it is necessary to conduct modeling in complex problem formulations, and still maintain the high level of parallelism effectiveness, even when thousands of processes are used. A solution of the problem of effective communication is extremely important when several computational grids with arbirtrary geometry of contacts between them are used in the calculation. The complexity of this task increases if an independent distribution of the grid nodes between processes is allowed. In this paper, a generalized approach is developed for processing contact conditions in terms of nodes reinterpolation from a given section of one grid to a certain area of the second grid. An efficient way of parallelization and establishing effective interprocess communications is proposed. For provided example problems we provide wave fileds and seismograms for both 2D and 3D formulations. It is shown that the algorithm can be realized both on Cartesian and on structured (curvilinear) computational grids. The considered statements demonstrate the possibility of carrying out calculations taking into account the surface topographies and curvilinear geometry of curvilinear contacts between the geological layers. Application of curvilinear grids allows to obtain more accurate results than when calculating only using Cartesian grids. The resulting parallelization efficiency is almost 100% up to 4096 processes (we used 128 processes as a basis to find efficiency). With number of processes larger than 4096, an expected gradual decrease in efficiency is observed. The rate of decline is not great, so at 16384 processes the parallelization efficiency remains at 80%.
The reported study was funded by RFBR according to the research project No. 18-07-00914 A. This work has been carried
out using computing resources of the federal collective usage center Complex for Simulation and Data Processing for Megascience Facilities at NRC “Kurchatov Institute”, http://ckp.nrcki.ru/.
Citation:
A. M. Ivanov, N. I. Khokhlov, “Parallel implementation of the grid-characteristic method in the case of explicit contact boundaries”, Computer Research and Modeling, 10:5 (2018), 667–678
\Bibitem{IvaKho18}
\by A.~M.~Ivanov, N.~I.~Khokhlov
\paper Parallel implementation of the grid-characteristic method in the case of explicit contact boundaries
\jour Computer Research and Modeling
\yr 2018
\vol 10
\issue 5
\pages 667--678
\mathnet{http://mi.mathnet.ru/crm678}
\crossref{https://doi.org/10.20537/2076-7633-2018-10-5-667-678}
Linking options:
https://www.mathnet.ru/eng/crm678
https://www.mathnet.ru/eng/crm/v10/i5/p667
This publication is cited in the following 13 articles:
V. Stetsyuk, N. Khokhlov, “Modeling the Wave Propagation near the Earth's Surface with Taking Relief into Account”, Phys. Part. Nuclei, 55:3 (2024), 525
V. S. Surov, “Mnogomernyi uzlovoi metod kharakteristik dlya giperbolicheskikh sistem”, Kompyuternye issledovaniya i modelirovanie, 13:1 (2021), 19–32
I. B. Petrov, P. V. Stognii, N. I. Khokhlov, “Mathematical modeling of 3D dynamic processes near a fracture using the Schoenberg fracture model”, Dokl. Math., 104:2 (2021), 254–257
Igor Petrov, Vasily Golubev, Alexey Shevchenko, 2021 Ivannikov Memorial Workshop (IVMEM), 2021, 61
N. Khokhlov, “A software package for modeling the propagation of dynamic wave disturbances in heterogeneous multi-scale media”, J. Phys.: Conf. Ser., 2056:1 (2021), 012012
Polina V. Stognii, Nikolay I. Khokhlov, Smart Innovation, Systems and Technologies, 214, Smart Modelling For Engineering Systems, 2021, 125
M. V. Muratov, V. V. Ryazanov, V. A. Biryukov, D. I. Petrov, I. B. Petrov, “Inverse Problems of Heterogeneous Geological Layers Exploration Seismology Solution by Methods of Machine Learning”, Lobachevskii J Math, 42:7 (2021), 1728
P. V. Stognii, I. B. Petrov, K. A. Beklemysheva, V. A. Miryaha, “Computer Exploration of the Ice Samples Strength Using Different Numerical Methods”, Lobachevskii J Math, 41:12 (2020), 2714
P. V. Stognii, N. I. Khokhlov, I. B. Petrov, “The numerical modeling of the elastic waves
propagation in the geological media with gas cavities using the grid-characteristic method”, Num. Anal. Appl., 13:3 (2020), 271–281
M. V. Muratov, V. A. Biryukov, I. B. Petrov, “Solution of the fracture detection problem by machine learning methods”, Dokl. Math., 101:2 (2020), 169–171
M. V. Muratov, V. A. Biryukov, I. B. Petrov, “Geological Fractures Detection by Methods of Machine Learning”, Lobachevskii J Math, 41:4 (2020), 533
N. I. Khokhlov, I. B. Petrov, “Primenenie setochno-kharakteristicheskogo metoda dlya resheniya zadach rasprostraneniya dinamicheskikh volnovykh vozmuschenii na vysokoproizvoditelnykh vychislitelnykh sistemakh”, Trudy ISP RAN, 31:6 (2019), 237–252
Katerina Beklemysheva, Alexey Ermakov, Alena Favorskaya, “Numerical simulation of cone object destruction under a short high-energy pulse”, Procedia Computer Science, 159 (2019), 1095