Parallel implementation of the grid-characteristic method in the case of explicit contact boundaries
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...
| Main Authors: | , |
|---|---|
| Format: | Article |
| Language: | Russian |
| Published: |
Institute of Computer Science
2018-10-01
|
| Series: | Компьютерные исследования и моделирование |
| Subjects: | |
| Online Access: | http://crm.ics.org.ru/uploads/crmissues/crm_2018_5/2018_05_07.pdf |
| _version_ | 1828750074771931136 |
|---|---|
| author | Andrey M. Ivanov Nikolai Igorevich Khokhlov |
| author_facet | Andrey M. Ivanov Nikolai Igorevich Khokhlov |
| author_sort | Andrey M. Ivanov |
| collection | DOAJ |
| description | 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%. |
| first_indexed | 2024-12-10T20:32:38Z |
| format | Article |
| id | doaj.art-997bc171d8694549b3f27728533ecfd6 |
| institution | Directory Open Access Journal |
| issn | 2076-7633 2077-6853 |
| language | Russian |
| last_indexed | 2024-12-10T20:32:38Z |
| publishDate | 2018-10-01 |
| publisher | Institute of Computer Science |
| record_format | Article |
| series | Компьютерные исследования и моделирование |
| spelling | doaj.art-997bc171d8694549b3f27728533ecfd62022-12-22T01:34:37ZrusInstitute of Computer ScienceКомпьютерные исследования и моделирование2076-76332077-68532018-10-0110566767810.20537/2076-7633-2018-10-5-667-6782735Parallel implementation of the grid-characteristic method in the case of explicit contact boundariesAndrey M. IvanovNikolai Igorevich KhokhlovWe 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%.http://crm.ics.org.ru/uploads/crmissues/crm_2018_5/2018_05_07.pdfparallel programminggrid-characteristic methodMPIstructured grids |
| spellingShingle | Andrey M. Ivanov Nikolai Igorevich Khokhlov Parallel implementation of the grid-characteristic method in the case of explicit contact boundaries Компьютерные исследования и моделирование parallel programming grid-characteristic method MPI structured grids |
| title | Parallel implementation of the grid-characteristic method in the case of explicit contact boundaries |
| title_full | Parallel implementation of the grid-characteristic method in the case of explicit contact boundaries |
| title_fullStr | Parallel implementation of the grid-characteristic method in the case of explicit contact boundaries |
| title_full_unstemmed | Parallel implementation of the grid-characteristic method in the case of explicit contact boundaries |
| title_short | Parallel implementation of the grid-characteristic method in the case of explicit contact boundaries |
| title_sort | parallel implementation of the grid characteristic method in the case of explicit contact boundaries |
| topic | parallel programming grid-characteristic method MPI structured grids |
| url | http://crm.ics.org.ru/uploads/crmissues/crm_2018_5/2018_05_07.pdf |
| work_keys_str_mv | AT andreymivanov parallelimplementationofthegridcharacteristicmethodinthecaseofexplicitcontactboundaries AT nikolaiigorevichkhokhlov parallelimplementationofthegridcharacteristicmethodinthecaseofexplicitcontactboundaries |