Material science, investigating the relationship between the structure of materials at atomic or molecular scales and their macroscopic properties, has driven the development of modern science and technology. The most important features of the microstructure of an engineering material are the crystalline defects that are manipulated to control its behavior. The following picture demonstrates different defects in crystal.
The primary step to simulation of crystal defects is to model the physical process in mathematics. It is impossible, however, to use only atomistic models to simulate all material properties because of the sheer number of atoms that would be required for such a simulation. Fortunately, since defects occupy only a small proportion of bulk crystals, it is more advisable to model the elastic field using more efficient continuum models wherever they are reasonable to apply. This project is devoted to propose an optimized coupling mechanism between not only atomistic and continuum, but rather between the (now distinct) atomistic and nonlinear continumm, or atomistic and linear elastic continuum for point crystal defects.
In physical language, the atomistic problem objects to find the atomistic configuration that minimise the total energy of the physical system. In mathematics, with a reference deformation on lattice , if we define the atomistic energy difference functional
then the atomistic problem is formulated as the following variational problem: find
However, due to the large number of atoms, it is unrealistic to perform defects simulation with the atomistic model.
Since the defects only occupy a small part of the lattice, the deformation is relatively 'smooth' outside the defect region, where we can approximate the atomistic potential in terms of nonlinear elasticity theory. Assume that the computational domain is decomposed as , and that is a regular mesh of and is the associated linear interpolation operator.
where and is the Cauchy-Born strain energy function
In order to couple the atomistic and continuum model in a consistent way, we aslo need to modify the potentials and iteractions on the interface. Here, we take a simple Geometric Reconstruction Atomistic Continuum (GRAC) method to achive the ghost force free coupling. In particular, we define GRAC energy
and the A/N coupling energy
then the A/N coupling problem seeks to find
The nonlinear elasticity model can be further approximated by its linearisation wherever the deformation is small. In this case, the whole space , the Cauchy Born energy function in is linearised as follows
If we define the linearisation energy
and also the A/N/L energy functional
then the A/N/L coupling problem seeks to find
We run a simulation of the Atomistic-Nonlinear Elasticity model to demonstrate the method of the steepest descent model. The setup defect is given on the first image. The animation contains 2 figures. The left figure shows the simulated defect and right shows the modification proposed for the next step to update the system on the left (although this simulation samples 1:10). We exit the simulation when the gradient
We acknowledge and thank the help of our supervisors Dr Christoph Ortner and Dr Andreas Dedner. We would also like to acknowledge the funding body EPSRC and the support from MASDOC CDT.
Contact: O.Dunbar at warwick.ac.uk, Yulong.Lu at warwick.ac.uk, Luke.D.Williams at warwick.ac.uk