Numerical Methods for Multilattices

Among the efficient numerical methods based on atomistic models, the quasicontinuum (QC) method has attracted growing interest in recent years. The QC method was first developed for crystalline materials with Bravais lattice and was later extended to multilattices (Tadmor et al, 1999). Another existing numerical approach to modeling multilattices is homogenization. In the present paper we review the existing numerical methods for multilattices and propose another concurrent macro-to-micro method in the numerical homogenization framework. We give a unified mathematical formulation of the new and the existing methods and show their equivalence. We then consider extensions of the proposed method to time-dependent problems and to random materials.Comment: 31 page

Non-local energetics of random heterogeneous lattices

In this paper, we study the mechanics of statistically non-uniform two-phase elastic discrete structures. In particular, following the methodology proposed in (Luciano and Willis, Journal of the Mechanics and Physics of Solids 53, 1505-1522, 2005), energetic bounds and estimates of the Hashin-Shtrikman-Willis type are developed for discrete systems with a heterogeneity distribution quantified by second-order spatial statistics. As illustrated by three numerical case studies, the resulting expressions for the ensemble average of the potential energy are fully explicit, computationally feasible and free of adjustable parameters. Moreover, the comparison with reference Monte-Carlo simulations confirms a notable improvement in accuracy with respect to approaches based solely on the first-order statistics.Comment: 32 pages, 8 figure

Computational Multiscale Methods

Computational Multiscale Methods play an important role in many modern computer simulations in material sciences with different time scales and different scales in space. Besides various computational challenges, the meeting brought together various applications from many disciplines and scientists from various scientific communities

Mechanics of Materials: Mechanics of Interfaces and Evolving Microstructure

Emphasis in modern day efforts in mechanics of materials is increasingly directed towards integration with computational materials science, which itself rests on solid physical and mathematical foundations in thermodynamics and kinetics of processes. Practical applications demand attention to length and time scales which are sufficiently large to preclude direct application of quantum mechanics approaches; accordingly, there are numerous pathways to mathematical modelling of the complexity of material structure during processing and in service. The conventional mathematical machinery of energy minimization provides guidance but has limited direct applicability to material systems evolving away from equilibrium. Material response depends on driving forces, whether arising from mechanical, electromagnetic, or thermal fields. When microstructures evolve, as during plastic deformation, progressive damage and fracture, corrosion, stress-assisted diffusion, migration or chemical/thermal aging, the associated classical mathematical frameworks are often ad hoc and heuristic. Advancing new and improved methods is a major focus of 21st century mechanics of materials of interfaces and evolving microstructure

Large scale ab-initio simulations of dislocations

We present a novel methodology to compute relaxed dislocations core configurations, and their energies in crystalline metallic materials using large-scale ab-intio simulations. The approach is based on MacroDFT, a coarse-grained density functional theory method that accurately computes the electronic structure with sub-linear scaling resulting in a tremendous reduction in cost. Due to its implementation in real-space, MacroDFT has the ability to harness petascale resources to study materials and alloys through accurate ab-initio calculations. Thus, the proposed methodology can be used to investigate dislocation cores and other defects where long range elastic effects play an important role, such as in dislocation cores, grain boundaries and near precipitates in crystalline materials. We demonstrate the method by computing the relaxed dislocation cores in prismatic dislocation loops and dislocation segments in magnesium (Mg). We also study the interaction energy with a line of Aluminum (Al) solutes. Our simulations elucidate the essential coupling between the quantum mechanical aspects of the dislocation core and the long range elastic fields that they generate. In particular, our quantum mechanical simulations are able to describe the logarithmic divergence of the energy in the far field as is known from classical elastic theory. In order to reach such scaling, the number of atoms in the simulation cell has to be exceedingly large, and cannot be achieved with the state-of-the-art density functional theory implementations

Atomistic Models of Materials: Mathematical Challenges

[no abstract available

Fracture of solar-grade anisotropic polycrystalline Silicon: A combined phase field–cohesive zone model approach

Artículo Open Access en el sitio web del editor. Pago por publicar en abierto.This work presents a novel computational framework to simulate fracture events in brittle anisotropic polycrystalline materials at the microscopical level, with application to solar-grade polycrystalline Silicon. Quasi-static failure is modeled by combining the phase field approach of brittle fracture (for transgranular fracture) with the cohesive zone model for the grain boundaries (for intergranular fracture) through the generalization of the recent FE-based technique published in [M. Paggi, J. Reinoso, Comput. Methods Appl. Mech. Engrg., 31 (2017) 145–172] to deal with anisotropic polycrystalline microstructures. The proposed model, which accounts for any anisotropic constitutive tensor for the grains depending on their preferential orientation, as well as an orientation-dependent fracture toughness, allows to simulate intergranular and transgranular crack growths in an efficient manner, with or without initial defects. One of the advantages of the current variational method is the fact that complex crack patterns in such materials are triggered without any user-intervention, being possible to account for the competition between both dissipative phenomena. In addition, further aspects with regard to the model parameters identification are discussed in reference to solar cells images obtained from transmitted light source. A series of representative numerical simulations is carried out to highlight the interplay between the different types of fracture occurring in solar-grade polycrystalline Silicon, and to assess the role of anisotropy on the crack path and on the apparent tensile strength of the material.Unión Europea FP/2007–2013/ERC 306622Ministerio de Economía y Competitividad MAT2015–71036-P y MAT2015–71309-PJunta de Andalucía P11-TEP-7093 y P12-TEP- 105