Search

This work presents an open source, dislocation density based crystal plasticity modeling framework, $\rho$-CP. A Kocks-type thermally activated flow is used for accounting for the temperature and strain rate effects on the crystallographic shearing rate. Slip system-level mobile and immobile dislocation densities, as well slip system-level backstress, are used as internal state variables for representing the substructure evolution during plastic deformation. A fully implicit numerical integration scheme is presented for the time integration of the finite deformation plasticity model. The framework is implemented and integrated with the open source finite element solver, Multiphysics Object-Oriented Simulation Environment (MOOSE). Example applications of the model are demonstrated for predicting the anisotropic mechanical response of single and polycrystalline hcp magnesium, strain rate effects and cyclic deformation of polycrystalline fcc OFHC copper, and temperature and strain rate effects on the thermo-mechanical deformation of polycrystalline bcc tantanlum. Simulations of realistic Voronoi-tessellated microstructures as well as Electron Back Scatter Diffraction (EBSD) microstructures are demonstrated to highlight the model's ability to predict large deformation and misorientation development during plastic deformation., Comment: 30 pages, 19 figures, 5 tables, v2

This work presents recent the progress in the development of the Concurrent Atomistic-Continuum (CAC) method for coarse-grained space- and time-resolved atomistic simulations of phonon transport. Application examples, including heat pulses propagating across grain boundaries and phase interfaces, as well as the interactions between phonons and moving dislocations, are provided to demonstrate the capabilities of CAC. The simulation results provide visual evidence and reveal the underlying physics of a variety of phenomena including: phonon focusing, wave interference, dislocation drag, interfacial Kapitza resistance caused by quasi-ballistic phonon transport, etc. A new method to quantify fluxes in transient transport processes is also introduced.

Progress in the development of coupled atomistic-continuum methods for simulations of critical dynamic material behavior has been hampered by a spurious wave reflection problem at the atomistic-continuum interface. This problem is mainly caused by the difference in material descriptions between the atomistic and continuum models, which results in a mismatch in phonon dispersion relations. In this work, we introduce a new method based on atomistic dynamics of lattice coupled with a concurrent atomistic-continuum method to enable a full phonon representation in the continuum description. This then permits the passage of short-wavelength, high-frequency phonon waves from the atomistic to continuum regions. The benchmark examples presented in this work demonstrate that the new scheme enables the passage of all allowable phonons through the atomistic-continuum interface; it also preserves the wave coherency and energy conservation after phonons transport across multiple atomistic-continuum interfaces. This work is the first step towards developing a concurrent atomistic-continuum simulation tool for non-equilibrium phonon-mediated thermal transport in materials with microstructural complexity.

The propagation of a heat pulse in a single crystal and across grain boundaries (GBs) is simulated using a concurrent atomistic-continuum method furnished with a coherent phonon pulse model. With a heat pulse constructed based on a Bose-Einstein distribution of phonons, this work has reproduced the phenomenon of phonon focusing in single and polycrystalline materials. Simulation results provide visual evidence that the propagation of a heat pulse in crystalline solids with or without GBs is partially ballistic and partially diffusive, i.e., there is a co-existence of ballistic and diffusive thermal transport, with the long-wavelength phonons traveling ballistically while the short-wavelength phonons scatter with each other and travel diffusively. To gain a quantitative understanding of GB thermal resistance, the kinetic energy transmitted across GBs is monitored on the fly and the time-dependent energy transmission for each specimen is measured; the contributions of coherent and incoherent phonon transport to the energy transmission are estimated. Simulation results reveal that the presence of GBs modifies the nature of thermal transport, with the coherent long-wavelength phonons dominating the heat conduction in materials with GBs. In addition, it is found the phonon-GB interaction can result in the reconstruction of the GBs.

Hydrogen embrittlement of metals is widely observed, but its atomistic origins remain little understood and much debated. Combining a unique identification of interstitial sites through polyhedral tessellation and first-principles calculations, we study hydrogen adsorption at grain boundaries in a variety of face-centered cubic metals of Ni, Cu, gamma-Fe and Pd. We discover the chemomechanical origin of variation of adsorption energetics for interstitial hydrogen at grain boundaries. A general chemomechanical formula is established to provide accurate assessments of hydrogen trapping and segregation energetics at grain boundaries, and it also offers direct explanations for certain experimental observations. The present study deepens our mechanistic understanding of the role of grain boundaries in hydrogen embrittlement, and promises a viable path towards predictive microstructure engineering against hydrogen embrittlement in structural metals.