University of Crete

The energy of a deformed crystal is calculated in the context of a central force lattice model in two dimensions. When the crystal shape is a lattice polygon, it is shown that the energy equals the bulk elastic energy, plus the boundary integral of a surface energy density, plus the sum over the vertices of a corner energy function. This is an exact result when the interatomic potential has finite range; for an infinite-range potential it is asymptotically valid as the lattice parameter tends to zero. The surface energy density is obtained explicitly as a function of the deformation gradient and boundary normal. The corner energy is found as an explicit function of the deformation gradient and the normals of the two facets meeting at the corner. A new bond counting approach is used, which reduces the problem to certain lattice point problems of number theory. The approach is then extended to more general convex regions with possibly curved boundary. The resulting surface energy density depends on the unit normal in a striking way. It is continuous at irrational directions, discontinuous at rational ones and nowhere differentiable. The method also yields an explicit interfacial energy for twin and phase boundaries.

Self adjusted meshes have important benefits approximating PDEs with solutions that exhibit nontrivial characteristics. When appropriately chosen, they lead to efficient, accurate and robust algorithms. Error control is also important, since appropriate analysis can provide guarantees on how accurate the approximate solution is through a posteriori estimates. Error control may lead to appropriate adaptive algorithms by identifying areas of large errors and adjusting the mesh accordingly. Error control and associated adaptive algorithms for important equations in Mathematical Physics is an open problem.

In this talk we consider the main structure of an algorithm which permits mesh redistribution with time and the nontrivial characteristics associated with it. We present improved algorithms and we discuss successful approaches towards error control for model problems (linear and nonlinear) of parabolic or hyperbolic type.