We present the first set of benchmark problems for phase field models that are being developed by the Center for Hierarchical Materials Design (CHiMaD) and the National Institute of Standards and Technology (NIST). While many scientific research areas use a limited set of well-established software, the growing phase field community continues to develop a wide variety of codes and lacks benchmark problems to consistently evaluate the numerical performance of new implementations. Phase field modeling has become significantly more popular as computational power has increased and is now becoming mainstream, driving the need for benchmark problems to validate and verify new implementations. We follow the example set by the micromagnetics community to develop an evolving set of benchmark problems that test the usability, computational resources, numerical capabilities and physical scope of phase field simulation codes. In this paper, we propose two benchmark problems that cover the physics of solute diffusion and growth and coarsening of a second phase via a simple spinodal decomposition model and a more complex Ostwald ripening model. We demonstrate the utility of benchmark problems by comparing the results of simulations performed with two different adaptive time stepping techniques, and we discuss the needs of future benchmark problems. The development of benchmark problems will enable the results of quantitative phase field models to be confidently incorporated into integrated computational materials science and engineering (ICME), an important goal of the Materials Genome Initiative.

Composite materials comprised of ferroelectric nanoparticles in a dielectric matrix are being actively investigated for a variety of functional properties attractive for a wide range of novel electronic and energy harvesting devices. However, the dependence of these functionalities on shapes, sizes, orientation and mutual arrangement of ferroelectric particles is currently not fully understood. In this study, we utilize a time-dependent Ginzburg-Landau approach combined with coupled-physics finite-element-method based simulations to elucidate the behavior of polarization in isolated spherical PbTiO or BaTiO nanoparticles embedded in a dielectric medium, including air. The equilibrium polarization topology is strongly affected by particle diameter, as well as the choice of inclusion and matrix materials, with monodomain, vortex-like and multidomain patterns emerging for various combinations of size and materials parameters. This leads to radically different polarization vs. electric field responses, resulting in highly tunable size-dependent dielectric properties that should be possible to observe experimentally. Our calculations show that there is a critical particle size below which ferroelectricity vanishes. For the PbTiO particle, this size is 2 and 3.4 nm, respectively, for high- and low-permittivity media. For the BaTiO particle, it is ∼3.6 nm regardless of the medium dielectric strength.

Cobalt-based alloys with γ/γ microstructures have the potential to become the next generation of superalloys, but alloy compositions and processing steps must be optimized to improve coarsening, creep, and rafting behavior. While these behaviors are different than in nickel-based superalloys, alloy development can be accelerated by understanding the thermodynamic factors influencing microstructure evolution. In this work, we develop a phase field model informed by first-principles density functional theory and experimental data to predict the equilibrium shapes of Co-Al-W γ precipitates. Three-dimensional simulations of single and multiple precipitates are performed to understand the effect of elastic and interfacial energy on coarsened and rafted microstructures; the elastic energy is dependent on the elastic stiffnesses, misfit strain, precipitate size, applied stress, and precipitate spatial distribution. We observe characteristic microstructures dependent on the type of applied stress that have the same γ morphology and orientation seen in experiments, indicating that the elastic stresses arising from coherent γ/γ interfaces are important for morphological evolution during creep. The results also indicate that the narrow γ channels between γ precipitates are energetically favored, and provide an explanation for the experimentally observed directional coarsening that occurs without any applied stress.