Testing out new PFHub upload mechanism

U.S. national research laboratories and agencies play an integral role in advancing science and technology for the public good. The authors of this article, as research software engineers (RSEs) and allies from eight unique national R&D organizations, came together to explore RSE needs from the perspective of national institutions. We identified three key areas of improvement for future RSEs to pursue science in the national interest: community establishment, hiring and retention, and recognition. To retain and cultivate this essential talent, U.S. national institutions must evolve to support appropriate career pathways for RSEs, and to recognize and reward RSEs’ work.

Cobalt-based superalloys with two phase γ/γ′ microstructures offer great promise as candidates for next-generation high-temperature alloys for applications, such as turbine blades. It is essential to understand the thermodynamic and kinetic factors that influence the microstructural evolution of these alloys in order to optimize the alloy compositions and processing steps with a goal to improve their coarsening, creep and rafting behavior. We are using a continuum phase field approach to study the diffusion process and to predict the equilibrium shapes of Co-Al-W γ′ precipitates. In order to obtain quantitatively predictive capabilities, we extract chemical free energies for the γ/γ′ phases based on CALculation of PHAse Diagrams (CALPHAD) thermodynamic data and diffusion mobilities for Co alloys based on CALPHAD kinetic data. We also use experimental or first-principles data for other quantities, such as misfit strain and interface information, for the parameterization of our model. A particular focus of our study is to understand how different energy balances, misfit strain and kinetics affect the coarsening and rafting behavior of γ′ precipitates, and the sensitivity of the final precipitate shape to materials parameters. We find that the equilibrium shape of the precipitate results from a delicate competition between chemical, interfacial, and elastic energies, and it is very sensitive to changes in model parameters. We examine how modeling input parameters affect the equilibrium shape of precipitates and relate these parameters to experimentally available values.

We present nucleation phase field model benchmark problems, expanding on our previous benchmark problems on diffusion, precipitation, dendritic growth, linear elasticity, fluid flow and electrochemistry. Nucleation is the process in which either a new thermodynamic phase or a new structure is created, such as solidification from the melt, or self-assembly of particulates. Based on where the nucleation occurs, it can be divided into two main categories: homogeneous nucleation and heterogeneous nucleation. In the first nucleation benchmark problem, we focus on homogeneous nucleation for both single seed under different initial conditions and multiple seeds. The second nucleation benchmark problem focuses on athermal heterogeneous nucleation and nucleation behavior near the free growth limit with different undercooling driving forces.

The COVID-19 pandemic has by-and-large prevented in-person meetings since March 2020. While the increasing deployment of effective vaccines around the world is a very positive development, the timeline and pathway to "normality" is uncertain and the "new normal" we will settle into is anyone's guess. Particle physics, like many other scientific fields, has more than a year of experience in holding virtual meetings, workshops, and conferences. A great deal of experimentation and innovation to explore how to execute these meetings effectively has occurred. Therefore, it is an appropriate time to take stock of what we as a community learned from running virtual meetings and discuss possible strategies for the future. Continuing to develop effective strategies for meetings with a virtual component is likely to be important for reducing the carbon footprint of our research activities, while also enabling greater diversity and inclusion for participation. This report summarizes a virtual two-day workshop on Virtual Meetings held May 5-6, 2021 which brought together experts from both inside and outside of high-energy physics to share their experiences and practices with organizing and executing virtual workshops, and to develop possible strategies for future meetings as we begin to emerge from the COVID-19 pandemic. This report outlines some of the practices and tools that have worked well which we hope will serve as a valuable resource for future virtual meeting organizers in all scientific fields.

In this work, we continue our development of phase field model benchmark problems with the addition of a third set, complimenting our previously developed problems for diffusion, precipitation, dendritic growth and linear elasticity. These benchmark problems are being jointly developed by the Center for Hierarchical Materials Design (CHiMaD) and the National Institute of Standards and Technology (NIST) along with input from the phase field community. The first problem in this third set targets Stokes flow, with a particular emphasis on flow around an obstruction placed asymmetrically in the domain. While Stokes flow is not traditionally in the canon of phase field problems, it is a class of problems gaining importance in areas such as filtration and water purification. The second problem deals with coupled Cahn-Hilliard diffusion and electrostatic forces, which is an important area in energy storage and battery sciences. We present our own solutions and discuss sources of numerical errors for the Stokes problem as well as simple checks to avoid fundamental issues in the coupled diffusion-electrostatics problem. The latter problem contains some subtleties that we expand on in an Appendix.

Scientific communities struggle with the challenge of effectively and efficiently sharing content and data. An online portal provides a valuable space for scientific communities to discuss challenges and collate scientific results. Examples of such portals include the Micromagnetic Modeling Group (μMAG) [1], the Interatomic Potentials Repository (IPR) [2, 3] and on a larger scale the NIH Genetic Sequence Database (GenBank) [4]. In this work, we present a description of a generic web portal that leverages existing online services to provide a framework that may be adopted by other small scientific communities. The first deployment of the PFHub framework supports phase-field practitioners and code developers participating in an effort to improve quality assurance for phase-field codes. Funding statement: D.W. wishes to acknowledge the Materials Genome Initiative funding allocated to the National Institute of Standards and Technology. S.J.D wishes to acknowledge funding from the U.S. Department of Energy, Office of Basic Energy Sciences, Division of Materials Sciences and Engineering under Award #DE-SC0008637 as part of the Center for PRedictive Integrated Structural Materials Science (PRISMS Center) at University of Michigan. P.W.V. is grateful for the financial assistance under the award 70NANB14H012 from the National Institute of Standards and Technology as part of the Center for Hierarchical Materials Design (CHiMaD).

We present the second set of benchmark problems for phase field models that are being jointly developed by the Center for Hierarchical Materials Design (CHiMaD) and the National Institute of Standards and Technology (NIST) along with input from other members in the phase field community. As the integrated computational materials engineering (ICME) approach to materials design has gained traction, there is an increasing need for quantitative phase field results. New algorithms and numerical implementations increase computational capabilities, necessitating standard problems to evaluate their impact on simulated microstructure evolution as well as their computational performance. We propose one benchmark problem for solidification and dendritic growth in a single-component system, and one problem for linear elasticity via the shape evolution of an elastically constrained precipitate. We demonstrate the utility and sensitivity of the benchmark problems by comparing the results of (1) dendritic growth simulations performed with different time integrators and (2) elastically constrained precipitate simulations with different precipitate sizes, initial conditions, and elastic moduli. These numerical benchmark problems will provide a consistent basis for evaluating different algorithms, both existing and those to be developed in the future, for accuracy and computational efficiency when applied to simulate physics often incorporated in phase field models.

We study the evolution of prior columnar β phase, interface L phase, and α phase during directional solidification of a Ti-6Al-4V melt pool. Finite element simulations estimate the solidification temperature and velocity fields in the melt pool and analyze the stress field and thermal distortions in the solidified part during the laser powder bed fusion process. A phase-field model uses the temperature and velocity fields to predict the formation of columnar prior-β(Ti) phase. During the solidification of β phase from an undercooled liquid, the residual liquid below the solidus temperature within the β columns results in α phase. The finite element simulated stress and strain fields are correlated with the length scales and volume fractions of the microstructure fields. Finally, the coalescence behavior of the β(Ti) cells during solidification is illustrated. The above analyses are important as they can be used for proactive control of the subsequent modeling of the heat treatment processes.

Cellular or dendritic microstructures that result as a function of additive manufacturing solidification conditions in a Ni-based melt pool are simulated in the present work using three-dimensional phase-field simulations. A macroscopic thermal model is used to obtain the temperature gradient G and the solidification velocity V which are provided as inputs to the phase-field model. We extract the cell spacings, cell core compositions, and cell tip as well as mushy zone temperatures from the simulated microstructures as a function of V. Cell spacings are compared with different scaling laws that correlate to the solidification conditions and approximated by G-mV-n. Cell core compositions are compared with the analytical solutions of a dendrite growth theory and found to be in good agreement. Through analysis of the mushy zone, we extract a characteristic bridging plane, where the primary γ phase coalesces across the intercellular liquid channels at a γ fraction between 0.6 and 0.7. The temperature and the γ fraction in this plane are found to decrease with increasing V. The simulated microstructural features are significant as they can be used as inputs for the simulation of subsequent heat treatment processes.

Ni-rich γ cells and Nb-rich eutectic droplets that form during laser power bed fusion solidification of Ni-Nb alloys are studied using experiments and simulations. Finite element simulations estimate the cooling rates in the melt pool and phase-field simulations predict the resulting cellular microstructures. The cell and droplet spacings are determined as a function of cooling rate and fit to a power law. The formation of Laves phase is predicted for a critical composition of Nb in the liquid droplets. Finally, our simulations demonstrate that anisotropy in the γ orientation influences the Laves fraction significantly.

We use single-track laser melting experiments and simulations on Inconel 625 to estimate the dimensions and microstructure of the resulting melt pool. Our work is based on a design-of-experiments approach which uses multiple laser power and scan speed combinations. Single-track experiments generated melt pools of certain dimensions that showed reasonable agreement with our finite-element calculations. Phase-field simulations were used to predict the size and segregation of the cellular microstructure that formed along the melt-pool boundaries for the solidification conditions that changed as a function of melt-pool dimensions.

In this study, an alloy phase-field model is used to simulate solidification microstructures at different locations within a solidified molten pool. The temperature gradient G and the solidification velocity V are obtained from a macroscopic heat transfer finite element simulation and provided as input to the phase-field model. The effects of laser beam speed and the location within the melt pool on the primary arm spacing and on the extent of Nb partitioning at the cell tips are investigated. Simulated steady-state primary spacings are compared with power law and geometrical models. Cell tip compositions are compared to a dendrite growth model. The extent of non-equilibrium interface partitioning of the phase-field model is investigated. Although the phase-field model has an anti-trapping solute flux term meant to maintain local interface equilibrium, we have found that during simulations it was insufficient at maintaining equilibrium. This is due to the fact that the additive manufacturing solidification conditions fall well outside the allowed limits of this flux term.

Numerical simulations are used in this work to investigate aspects of microstructure and microseg-regation during rapid solidification of a Ni-based superalloy in a laser powder bed fusion additive manufacturing process. Thermal modeling by finite element analysis simulates the laser melt pool, with surface temperatures in agreement with in situ thermographic measurements on Inconel 625. Geometric and thermal features of the simulated melt pools are extracted and used in subsequent mesoscale simulations. Solidification in the melt pool is simulated on two length scales. For the multicomponent alloy Inconel 625, microsegregation between dendrite arms is calculated using the Scheil-Gulliver solidification model and DICTRA software. Phase-field simulations, using Ni-Nb as a binary analogue to Inconel 625, produced microstructures with primary cellular/dendritic arm spacings in agreement with measured spacings in experimentally observed microstructures and a lesser extent of microsegregation than predicted by DICTRA simulations. The composition profiles are used to compare thermodynamic driving forces for nucleation against experimentally observed precipitates identified by electron and X-ray diffraction analyses. Our analysis lists the precipitates that may form from FCC phase of enriched interdendritic compositions and compares these against experimentally observed phases from 1 h heat treatments at two temperatures: stress relief at 1143 K (870 °C) or homogenization at 1423 K (1150 °C).

This study presents windowless CdSe/CdTe thin film photovoltaic devices with in-plane patterning at a submicrometer length scale. The photovoltaic cells are fabricated upon two interdigitated comb electrodes prepatterned at micrometer length scale on an insulating substrate. CdSe is electrodeposited on one electrode, and CdTe is deposited by pulsed laser deposition over the entire surface of the resulting structure. Previous studies of symmetric devices are extended in this study. Specifically, device performance is explored with asymmetric devices having fixed CdTe contact width and a range of CdSe contact widths, and the devices are fabricated with improved dimensional tolerance. Scanning photocurrent microscopy (also known as laser beam induced current mapping) is used to examine local current collection efficiency, providing information on the spatial variation of performance that complements current-voltage and external quantum efficiency measurements of overall device performance. Modeling of carrier transport and recombination indicates consistency of experimental results for local and blanket illumination. Performance under simulated air mass 1.5 illumination exceeds 5% for all dimensions examined, and the best-performing device achieved 5.9% efficiency.

We develop a multicomponent phase-field model specially formulated to robustly simulate concentration variations from molar to atomic magnitudes across an interface, i.e., partition coefficients in excess of 10±23 such as may be the case with species which are predominant in one phase and insoluble in the other. Substitutional interdiffusion on a normal lattice and concurrent interstitial diffusion are included. The composition in the interface follows the approach of Kim, Kim, and Suzuki [Phys. Rev. E 60, 7186 (1999)] and is compared to that of Wheeler, Boettinger, and McFadden [Phys. Rev. A 45, 7424 (1992)] in the context of large partitioning. The model successfully reproduces analytical solutions for binary diffusion couples and solute trapping for the demonstrated cases of extremely large partitioning.