The phase-field model of crystallographic slip of Ortiz and Koslowski et al was extended here to slip processes that required the activation of multiple slip systems, and the resulting model was applied to the investigation of finite twist boundary arrays. The distribution of slip over a slip plane was described by means of multiple integer-valued phase fields. It was shown how all the terms in the total energy of the crystal, including the long-range elastic energy and the Peierls interplanar energy, could be written explicitly in terms of the multi-phase field. The model was used to ascertain stable dislocation structures arising in an array of finite twist boundaries. These structures were found to consist of regular square or hexagonal dislocation networks separated by complex dislocation pile-ups over the intervening transition layers.

A Multi-Phase Field Model of Planar Dislocation Networks. M.Koslowski, M.Ortiz: Modelling and Simulation in Materials Science and Engineering, 2004, 12[6], 1087-97