Engineering Papers⌕ Search

Engineering topics

Voorhees, P. W.

Publications and source records attributed to Voorhees, P. W..

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Phase field benchmark problems for nucleation

In this work, 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.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Phase field benchmark problems targeting fluid flow and electrochemistry

Here, 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.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Phase field benchmark problems for dendritic growth and linear elasticity

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 solidifiication 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. As a result, 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.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Predicting the morphologies of γ' precipitates in cobalt-based superalloys

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. Here, 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.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Coarsening of complex microstructures following spinodal decomposition

Coarsening plays a pivotal role in materials engineering, but our understanding of the dynamics of coarsening in morphologically complex systems is still limited. In this paper, we examine the correlations between the interfacial velocity and interfacial morphologies, and then predict the evolution of mean curvature based on the correlations. Three simulated structures with varying volume fractions, two bicontinuous and one nonbicontinuous, are generated using the Cahn-Hilliard equation. We find general correlations between interfacial velocity and mean curvature, as well as between interfacial velocity and the surface Laplacian of the mean curvature. Furthermore, we find that the probability of finding a patch of interface with a given normal velocity and the same local principal curvatures is described well by a Gaussian distribution, independent of the principal curvature values and the volume fractions of the structures. We also find that average interfacial velocity is described by a polynomial of the mean curvature and the net curvature. Based on this finding, we develop a semi-analytical approach to predicting the rate of change of the mean curvature, which determines the morphological evolution of complex microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Benchmark problems for numerical implementations of phase field models

Here, 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.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Ostwald ripening of faceted Si particles in an Al-Si-Cu melt

The microstructural evolution of an Al-Si-Cu alloy during Ostwald ripening is imaged via synchrotron-based, four-dimensional (i.e., space and time resolved) X-ray tomography. Samples of composition Al-32 wt%Si-15 wt%Cu were annealed isothermally at 650 °C, in the two-phase solid-liquid regime, while tomographic projections were collected in situ over the course of five hours. Advances in experimental methods and computational approaches enable us to characterize the local interfacial curvatures and velocities during ripening. The sequence of three-dimensional reconstructions and interfacial shape distributions shows highly faceted Si particles in a copper-enriched liquid, that become increasingly isotropic or rounded over time. In addition, we find that the coarsening rate constant is approximately the same in the binary and ternary systems. By coupling these experimental measurements with CALPHAD modeling and ab initio molecular dynamics simulation, we assess the influence of Cu on the coarsening process. Lastly, we find the unusual “pinning” of microstructure at the junction between rough and smooth interfaces and suggest a mechanism for this behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2015

The three-dimensional morphology of growing dendrites

The processes controlling the morphology of dendrites have been of great interest to a wide range of communities, since they are examples of an out-of-equilibrium pattern forming system, there is a clear connection with battery failure processes, and their morphology sets the properties of many metallic alloys. We determine the three-dimensional morphology of free growing metallic dendrites using a novel X-ray tomographic technique that improves the temporal resolution by more than an order of magnitude compared to conventional techniques. These measurements show that the growth morphology of metallic dendrites is surprisingly different from that seen in model systems, the morphology is not self-similar with distance back from the tip, and that this morphology can have an unexpectedly strong influence on solute segregation in castings. These experiments also provide benchmark data that can be used to validate simulations of free dendritic growth.

36 MATERIALS SCIENCE↗
NASA NTRS2013

Coarsening in Solid-liquid Mixtures: Overview of Experiments on Shuttle and ISS

The microgravity environment on the Shuttle and the International Space Station (ISS) provides the ideal condition to perform experiments on Coarsening in Solid-Liquid Mixtures (CSLM) as deleterious effects such as particle sedimentation and buoyancy-induced convection are suppressed. For an ideal system such as Lead-Tin in which all the thermophysical properties are known, the initial condition in microgravity of randomly dispersed particles with local clustering of solid Tin in eutectic liquid Lead-Tin matrix, permitted kinetic studies of competitive particle growth for a range of volume fractions. Verification that the quenching phase of the experiment had negligible effect of the spatial distribution of particles is shown through the computational solution of the dynamical equations of motion, thus insuring quench-free effects from the coarsened microstructure measurements. The low volume fraction experiments conducted on the Shuttle showed agreement with transient Ostwald ripening theory, and the steady-state requirement of LSW theory was not achieved. More recent experiments conducted on ISS with higher volume fractions have achieved steady-state condition and show that the kinetics follows the classical diffusion limited particle coarsening prediction and the measured 3D particle size distribution becomes broader as predicted from theory.

Ostwald ripening↗
NASA NTRS1993

Morphological instability in epitaxially strained dislocation-free solid films - Linear stability theory

The morphological instability of a growing epitaxially strained dislocation-free solid film is analyzed. An evolution equation for the film surface is derived in the dilute limit of vacancies based on surface diffusion driven by a stress-dependent chemical potential. From the time-dependent linear stability problem the conditions for which a growing film is unstable are determined. It is found that the instability is driven by the lattice mismatch between the film and the substrate; however, low temperatures as well as elastically stiff substrates are stabilizing influences. The results also reveal that the critical film thickness for instability depends on the growth rate of the film itself. Detailed comparison with experimental observations indicates that the instability described exhibits many of the observed features of the onset of the 'island instability'.

Spencer, B. J.↗
NASA NTRS1992

On the morphological development of second-phase particles in elastically-stressed solids

We examine the equilibrium shape and the dynamics of the morphological evolution of a coherent misfitting particle in an elastically anisotropic medium. No a priori assumptions are made on the possible particle morphologies; the particle evolves in a manner consistent with both the diffusion and elastic fields surrounding the particle and the thermodynamics of interfaces in stressed solids. Through these calculations we find the thermodynamically stable equilibrium shape of a misfitting particle. By thermodynamically stable we mean that all the equilibrium conditions, elastic, chemical and interfacial are satisfied simultaneously, in contrast to previous treatments which determined the equilibrium shape by minimizing the sum of the elastic and interfacial energies for certain classes of particle shapes. We find that the equilibrium particle morphologies are not simple geometric shapes, have fourfold symmetry and a continuously varying interfacial curvature with position along the interface. Furthermore, the results show that the times required for a particle to evolve to its equilibrium morphology are likely to be within those which are accessible experimentally. In addition, we develop a general approach for determining the equilibrium shape of a particle in an elastically stressed solid.

Voorhees, P. W.↗
NASA NTRS1992

The effect of compositionally-generated elastic stresses on morphological instability during directional solidification

The effect of compositionally-generated elastic stresses on the conditions for morphological instability during directional solidification of a dilute, ideal binary alloy is investigated using the Gibbs-Thomson equations for an elastically-stressed solid with zero surface stress in equilibrium with a liquid. It is found that these stresses lead to a small stabilization of the Mullins and Sekerka cellular mode of instability. The steady mode is stabilized by a stress-induced modification to the interfacial concentration of the solid, which in turn alters the amount of solute rejected to inhibit the growth of perturbations. The presence of elastic stresses could generate a new oscillatory instability which is most likely to be found in experiments near absolute stability for materials with segregation coefficients near unity and large solute expansion coefficients.

Spencer, B. J.↗
NASA NTRS1991

Measurement and analysis of grain boundary grooving by volume diffusion

Experimental measurements of isothermal grain boundary grooving by volume diffusion are carried out for Sn bicrystals in the Sn-Pb system near the eutectic temperature. The dimensions of the groove increase with a temporal exponent of 1/3, and measurement of the associated rate constant allows the determination of the product of the liquid diffusion coefficient D and the capillarity length Gamma associated with the interfacial free energy of the crystal-melt interface. The small-slope theory of Mullins is generalized to the entire range of dihedral angles by using a boundary integral formulation of the associated free boundary problem, and excellent agreement with experimental groove shapes is obtained. By using the diffusivity measured by Jordon and Hunt, the present measured values of Gamma are found to agree to within 5 percent with the values obtained from experiments by Gunduz and Hunt on grain boundary grooving in a temperature gradient.

Hardy, S. C.↗
NASA NTRS1990

Coarsening in binary solid-liquid mixtures

A theory of Ostwald ripening has been developed for a solid-liquid mixture cosisting of a low volume fraction array of spherical solid particles in a liquid wherein the coarsening process proceeds via the transport of both heat and mass. It is found that the simultaneous transport of heat and mass during ripening does not alter the exponents of the temporal power laws governing the ripening process from their classical values but does alter the amplitudes of these power laws. The growth rate of the cube of the average particle radius, the rate constant, is found to depend both on the alloy solute concentration and the ratio of the thermal to solutal diffusivities. In most metallic systems, a large decrease in the rate constant can be expected with small additions of solute to a pure metal. Possible extensions of this theory to the analogous problem of ripening in isothermal ternary alloys are also discussed.

Voorhees, P. W.↗