Hostname: page-component-586b7cd67f-gb8f7 Total loading time: 0 Render date: 2024-11-23T10:42:45.468Z Has data issue: false hasContentIssue false

About an analytical verification of quasi-continuum methods with Γ-convergence techniques

Published online by Cambridge University Press:  05 April 2013

Mathias Schäffner
Affiliation:
Chair for Mathematics in the Sciences, Department of Mathematics, University of Würzburg, Emil-Fischer-Str. 40, D-97074 Würzburg, Germany
Anja Schlömerkemper
Affiliation:
Chair for Mathematics in the Sciences, Department of Mathematics, University of Würzburg, Emil-Fischer-Str. 40, D-97074 Würzburg, Germany
Get access

Abstract

Quasi-continuum (QC) methods are computational techniques, which reduce the complexity of atomistic simulations in a static setting while keeping information on small-scale structures and effects. The main idea is to couple atomistic and continuum models and thus to obtain quite detailed but still not too expensive numerical simulations.

We aim at a mathematically rigorous verification of QC methods by means of discrete to continuum limits. In this article we present our first results for the so-called quasi-nonlocal QC method in the context of fracture mechanics. To this end we start from a one-dimensional chain of atoms with nearest and next-to-nearest neighbour interactions of Lennard-Jones type. This is considered as a fully atomistic model of which the Γ-limits (of zeroth and first order) for an infinite number of atoms are known [7].

The QC models we construct are equal to this fully atomistic model in the atomistic region; in the continuum regime we approximate the next-to-nearest neighbour interactions by some nearest neighbour potential which is related to the so-called Cauchy-Born rule. Further we choose certain representative atoms in order to coarsen the mesh in the continuum region. It turns out that the selection of the representative atoms is crucial and influences the Γ-limits.

We regard a QC model as good if the Γ-limits of zeroth and first order or at least their minimal values and minimizers are the same as those of the fully atomistic model. Our analysis shows that, while in an elastic regime only the size of the atomistic region matters, in the case of fracture a proper choice of the representative atoms is an essential ingredient.

Type
Articles
Copyright
Copyright © Materials Research Society 2013 

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

REFERENCES

Blanc, X., Le Bris, C. and Legoll, F., Analysis of a prototypical multiscale method coupling atomistic and continuum mechanics, M2AN Math. Model. Numer. Anal. 39 (2005).CrossRefGoogle Scholar
Braides, A., Γ-convergence for beginners, Oxford Univ. Press, 2002.CrossRefGoogle Scholar
Braides, A. and Cicalese, M., Surface energies in nonconvex discrete systems, Math. Models Methods Appl. Sci. 17 (2007), 9851037.CrossRefGoogle Scholar
W. E, Principles of multiscale modelling, Cambridge University Press, Cambridge, 2011.Google Scholar
Español, M., Kochmann, D., Conti, S. and Ortiz, M., A Γ-convergence analysis of the quasicontinuum method, Submitted to SIAM Multiscale Modeling and Simulation.Google Scholar
Ortner, C., A priori and a posteriori analysis of the quasi-nonlocal quasi-continuum method in 1D, Math. Comp. 80 (2011), 12651285.CrossRefGoogle Scholar
Scardia, L., Schlömerkemper, A. and Zanini, C., Boundary layer energies for nonconvex discrete systems, Math. Models Methods Appl. Sci., 21 (2011), 777817.Google Scholar
Schäffner, M., Schlömerkemper, A., A Γ-convergence approach to quasi-continuum methods in 1D, in preparation.Google Scholar
Tadmor, E.B. and Miller, R.E., Modeling Materials: Continuum, Atomistic and Multiscale Techniques, Cambridge University Press, 2011.CrossRefGoogle Scholar