Analysis of a Sinclair-type domain decomposition solver for atomistic/continuum coupling
The "flexible boundary condition" method, introduced by Sinclair and coworkers in the 1970s, remains among the most popular methods for simulating isolated two-dimensional crystalline defects, embedded in an effectively infinite atomistic domain. In essence, the method can be characterized as a doma...
Saved in:
| Main Author: | |
|---|---|
| Format: | Article |
| Published: |
2019
|
| Subjects: | |
| Online Access: | https://arxiv.org/abs/1912.10530 http://arxiv.org/abs/1912.10530 |
| Tags: |
No Tags, Be the first to tag this record!
|
| _version_ | 1879173296882188288 |
|---|---|
| author | Hodapp, M. |
| author_browse | Hodapp, M. |
| author_facet | Hodapp, M. |
| author_sort | Hodapp, M. |
| building | arXiv |
| collection | Acceso abierto |
| description | The "flexible boundary condition" method, introduced by Sinclair and coworkers in the 1970s, remains among the most popular methods for simulating isolated two-dimensional crystalline defects, embedded in an effectively infinite atomistic domain. In essence, the method can be characterized as a domain decomposition method which iterates between a local anharmonic and a global harmonic problem, where the latter is solved by means of the lattice Green function of the ideal crystal. This local/global splitting gives rise to tremendously improved convergence rates over related alternating Schwarz methods. In a previous publication (Hodapp et al., 2019, Comput. Methods in Appl. Mech. Eng. 348), we have shown that this method also applies to large-scale three-dimensional problems, possibly involving hundreds of thousands of atoms, using fast summation techniques exploiting the low-rank nature of the asymptotic lattice Green function.
Here, we generalize the Sinclair method to bounded domains and develop an implementation using a discrete boundary element method to correct the infinite solution with respect to a prescribed far-field condition, thus preserving the advantage of the original method of not requiring a global spatial discretization. Moreover, we present a detailed convergence analysis and show for a one-dimensional problem that the method is unconditionally stable under physically motivated assumptions. To further improve the convergence behavior, we develop an acceleration technique based on a relaxation of the transmission conditions between the two subproblems. Numerical examples for linear and nonlinear problems are presented to validate the proposed methodology. |
| doi_str_mv | 10.1137/19M130861X |
| first_indexed | 2026-10-15T23:55:05Z |
| format | Article |
| fullrecord | <oai_dc:dc xmlns:oai_dc="http://www.openarchives.org/OAI/2.0/oai_dc/" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://www.openarchives.org/OAI/2.0/oai_dc/ https://www.openarchives.org/OAI/2.0/oai_dc.xsd"><identifier>arxiv-1912.10530</identifier><datestamp>2021-10-15</datestamp>
<dc:title>Analysis of a Sinclair-type domain decomposition solver for atomistic/continuum coupling</dc:title>
<dc:creator>Hodapp, M.</dc:creator>
<dc:subject>Numerical Analysis</dc:subject>
<dc:subject>Materials Science</dc:subject>
<dc:subject>65N55, 70C20, 65N80</dc:subject>
<dc:description>The "flexible boundary condition" method, introduced by Sinclair and coworkers in the 1970s, remains among the most popular methods for simulating isolated two-dimensional crystalline defects, embedded in an effectively infinite atomistic domain. In essence, the method can be characterized as a domain decomposition method which iterates between a local anharmonic and a global harmonic problem, where the latter is solved by means of the lattice Green function of the ideal crystal. This local/global splitting gives rise to tremendously improved convergence rates over related alternating Schwarz methods. In a previous publication (Hodapp et al., 2019, Comput. Methods in Appl. Mech. Eng. 348), we have shown that this method also applies to large-scale three-dimensional problems, possibly involving hundreds of thousands of atoms, using fast summation techniques exploiting the low-rank nature of the asymptotic lattice Green function.
Here, we generalize the Sinclair method to bounded domains and develop an implementation using a discrete boundary element method to correct the infinite solution with respect to a prescribed far-field condition, thus preserving the advantage of the original method of not requiring a global spatial discretization. Moreover, we present a detailed convergence analysis and show for a one-dimensional problem that the method is unconditionally stable under physically motivated assumptions. To further improve the convergence behavior, we develop an acceleration technique based on a relaxation of the transmission conditions between the two subproblems. Numerical examples for linear and nonlinear problems are presented to validate the proposed methodology. </dc:description>
<dc:description>accepted for publication in SIAM Multiscale Modeling and Simulation</dc:description>
<dc:date>2019-12-22</dc:date>
<dc:date>2021-03-12</dc:date>
<dc:type>text</dc:type>
<dc:identifier>https://arxiv.org/abs/1912.10530</dc:identifier>
<dc:identifier>http://arxiv.org/abs/1912.10530</dc:identifier>
<dc:identifier>Multiscale Model. Simul., 19(4), 1499-1537 (2021)</dc:identifier>
<dc:identifier>doi:10.1137/19M130861X</dc:identifier>
</oai_dc:dc> |
| id | arxiv-1912.10530 |
| institution | arXiv |
| last_indexed | 2026-10-15T23:55:05Z |
| publishDate | 2019 |
| publishDateRange | 2019 |
| publishDateSort | 2019 |
| record_format | oai_dc |
| spelling | arxiv-1912.105302021-10-15 Analysis of a Sinclair-type domain decomposition solver for atomistic/continuum coupling Hodapp, M. Numerical Analysis Materials Science 65N55, 70C20, 65N80 The "flexible boundary condition" method, introduced by Sinclair and coworkers in the 1970s, remains among the most popular methods for simulating isolated two-dimensional crystalline defects, embedded in an effectively infinite atomistic domain. In essence, the method can be characterized as a domain decomposition method which iterates between a local anharmonic and a global harmonic problem, where the latter is solved by means of the lattice Green function of the ideal crystal. This local/global splitting gives rise to tremendously improved convergence rates over related alternating Schwarz methods. In a previous publication (Hodapp et al., 2019, Comput. Methods in Appl. Mech. Eng. 348), we have shown that this method also applies to large-scale three-dimensional problems, possibly involving hundreds of thousands of atoms, using fast summation techniques exploiting the low-rank nature of the asymptotic lattice Green function. Here, we generalize the Sinclair method to bounded domains and develop an implementation using a discrete boundary element method to correct the infinite solution with respect to a prescribed far-field condition, thus preserving the advantage of the original method of not requiring a global spatial discretization. Moreover, we present a detailed convergence analysis and show for a one-dimensional problem that the method is unconditionally stable under physically motivated assumptions. To further improve the convergence behavior, we develop an acceleration technique based on a relaxation of the transmission conditions between the two subproblems. Numerical examples for linear and nonlinear problems are presented to validate the proposed methodology. accepted for publication in SIAM Multiscale Modeling and Simulation 2019-12-22 2021-03-12 text https://arxiv.org/abs/1912.10530 http://arxiv.org/abs/1912.10530 Multiscale Model. Simul., 19(4), 1499-1537 (2021) doi:10.1137/19M130861X |
| spellingShingle | Numerical Analysis Materials Science 65N55, 70C20, 65N80 Hodapp, M. Analysis of a Sinclair-type domain decomposition solver for atomistic/continuum coupling |
| title | Analysis of a Sinclair-type domain decomposition solver for atomistic/continuum coupling |
| title_full | Analysis of a Sinclair-type domain decomposition solver for atomistic/continuum coupling |
| title_fullStr | Analysis of a Sinclair-type domain decomposition solver for atomistic/continuum coupling |
| title_full_unstemmed | Analysis of a Sinclair-type domain decomposition solver for atomistic/continuum coupling |
| title_short | Analysis of a Sinclair-type domain decomposition solver for atomistic/continuum coupling |
| title_sort | analysis of a sinclair type domain decomposition solver for atomistic continuum coupling |
| topic | Numerical Analysis Materials Science 65N55, 70C20, 65N80 |
| url | https://arxiv.org/abs/1912.10530 http://arxiv.org/abs/1912.10530 |
| work_keys_str_mv | AT hodappm analysisofasinclairtypedomaindecompositionsolverforatomisticcontinuumcoupling |