Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase field modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Ginzburg-Landau theory of the bcc-liquid interface kinetic coefficient

We extend the Ginzburg-Landau (GL) theory of atomically rough bcc-liquid interfaces outside of equilibrium. We use this extension to derive an analytical expression for the kinetic coefficient, which is the proportionality constant μ(nˆ) between the interface velocity along a direction nˆ normal to the interface and the interface undercooling. The kinetic coefficient is expressed as a spatial integral along the normal direction of a sum of gradient square terms corresponding to different nonlinear density wave profiles. Anisotropy arises naturally from the dependence of those profiles on the angles between the principal reciprocal lattice vectors K → i and nˆ. Values of the kinetic coefficient for the (100), (110), and (111) interfaces are compared quantitatively to the prediction of linear Mikheev-Chernov (MC) theory and previous molecular dynamics (MD) simulation studies of crystallization kinetics for a classical model of Fe. Additional MD simulations are carried out here to compute the relaxation time of density waves in the liquid in order to make this comparison free of fit parameters. The GL theory predicts an expression for μ similar to the MC theory but yields a better agreement with MD simulations for both its magnitude and anisotropy due to a fully nonlinear description of density wave profiles across the solid-liquid interface. In particular, the overall magnitude of μ predicted by GL theory is an order of magnitude larger than predicted by the MC theory. GL theory is also used to derive an inverse relation between μ and the solid-liquid interfacial free energy. Here, the general methodology used here to derive an expression for μ(nˆ) also applies to amplitude equations derived from the phase-field-crystal model, which only differ from GL theory by the choice of cubic and higher order nonlinearities in the free-energy density.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Microstructure prediction for Ti-22Al-25Nb in laser powder bed fusion

This work presents a physics-informed framework for predicting solidification morphology and defect susceptibility in additively manufactured Ti–22Al–25Nb across a broad processing space. The framework integrates solidification microstructure selection (SMS) analysis with a single-track defect-based printability map to establish a unified methodology linking processing parameters to both interfacial morphology and manufacturability. Thermal gradients G and solidification rates R are first computed using the Thermo-Calc Additive Manufacturing (TC-AM) module, a finite-interface-dissipation (FID) phase-field (PF) model coupled with CALPHAD method is then employed to systematically distinguish planar and dendritic regimes as functions of $G$ and $R$. By superimposing the printability map onto the morphology projections, a comprehensive process–structure framework is obtained. Across most processing conditions, the predicted microstructure is predominantly dendritic, while planar growth emerges only under selected laser power $P$ and scan speed $v$ combinations. In addition to morphology classification, the framework quantifies the dendritic area fraction and introduces a width-based morphology descriptor to characterize the spatial extent of planar/dendritic regions within the melt pool. It provides mechanistic insight into the interplay between solidification physics and defect formation, offering practical guidance for parameter selection and microstructural control in Ti–22Al–25Nb additive manufacturing (AM).

36 MATERIALS SCIENCE↗
DOE OSTI2025

Mesoscale Modeling of Dislocation Cell Structure Evolution and Radiation-Induced Segregation in Additively Manufactured Austenitic Stainless Steel

Structural alloys under irradiation develop radiation-induced segregation (RIS) at point defect sinks, which lead to undesired changes in the alloy's properties. Additively manufactured austenitic stainless steels are expected to yield a distinct response to irradiation damage owing to their unique as-printed dislocation cell structure. In this talk, we present the development of a mesoscale model for dislocation cell structure evolution coupled with RIS using the MOOSE framework. The evolution of network dislocations within spatially distinct cells and cell walls are modeled using climb-mediated edge annihilation and generation processes. Starting with initial microstructures of as-printed non-equilibrium segregations at dislocation cell walls, our 2D and 3D simulations explore the effects of irradiation temperature, dose rate, dislocation sink bias, and specimen thickness on microstructure evolution. The results will be compared against experimental characterizations of in-situ and ex-situ ion-irradiated samples, and the implications of evolving cell structure on irradiation damage response will be discussed.

additive manufacturing↗
DOE OSTI2025

In situ irradiation of spent nuclear fuels

To improve the economics of commercial nuclear reactors, nuclear vendors and utilities are seeking approval for increased burnup usage of the existing nuclear fleet. This is especially critical for meeting the clean energy initiative by the United States Government, Department of Energy (DOE). However, one of the key challenges the nuclear industry faces in this regard is that the regions exposed to high burnup and low temperatures exhibit a fine-grained microstructure with large bubbles known as high-burnup structure (HBS) [1]. The formation of HBS has been correlated to the diminished performance of the reactor, as well as fuel fragmentation and pulverization during transient and accidental conditions [2]. Therefore, it is paramount to understand the mechanisms for HBS formation along with its impact on the properties and performance of nuclear fuels. While existing programs, such as Nuclear Energy Advanced Modeling and Simulation (NEAMS) and Advanced Fuel Campaign (AFC) are focusing on evaluating the performance impact of HBS, the physical mechanisms contributing to HBS formation are still not fully understood. In addition, having predictive capabilities and sound understanding of the microstructural evolution of nuclear fuel is essential for accelerated development, qualification, and deployment of new nuclear materials and novel reactor designs for advanced nuclear reactors. There is a lack of consensus among researchers regarding the mechanisms leading to such restructuring observed in HBS. Grain subdivision due to polygonization versus recrystallization, continuous versus discrete recrystallization occurring in tandem or conjunction, etc., have been proposed and debated. In general, it is hypothesized that defect accumulation and dislocation interaction within the grains cause the realignment of dislocations into grain boundaries, leading to the new subgrain formation, which over time transforms into new grains. However, due to the lack of transient data, the importance of fission rate, irradiation, thermal, and stress history of the fuel on the restructuring could not be assessed. In situ microstructural evolution under various irradiation conditions is desired to bridge this gap. Alternatively, phase-field-based models have been developed to capture HBS formation via discrete recrystallization utilizing the classical nucleation approach [3–5]. However, in these models, the grain nucleation criteria are often defined based on empirical relations for burnup and fission gas density leading to dislocation density change. A mechanistic approach to capture the dislocation interaction with the microstructural features leading to grain subdivision is lacking.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

Density functional theory for carbon dioxide crystal

We present a density functional approach to describe the solid−liquid phase transition, interfacial and crystal structure, and properties of polyatomic CO{sub 2}. Unlike previous phase field crystal model or density functional theory, which are derived from the second order direct correlation function, the present density functional approach is based on the fundamental measure theory for hard-sphere repulsion in solid. More importantly, the contributions of enthalpic interactions due to the dispersive attractions and of entropic interactions arising from the molecular architecture are integrated in the density functional model. Using the theoretical model, the predicted liquid and solid densities of CO{sub 2} at equilibrium triple point are in good agreement with the experimental values. Based on the structure of crystal-liquid interfaces in different planes, the corresponding interfacial tensions are predicted. Their respective accuracies need to be tested.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Thermodynamic origin of nonvolatility in resistive memory

Electronic switches based on the migration of high-density point defects, or memristors, are poised to revolutionize post-digital electronics. Despite significant research, key mechanisms for filament formation and oxygen transport remain unresolved, hindering our ability to predict and design device properties. For example, experiments have achieved 10 orders of magnitude longer retention times than predicted by current models. Here, using electrical measurements, scanning probe microscopy, and first-principles calculations on tantalum oxide memristors, we reveal that the formation and stability of conductive filaments crucially depend on the thermodynamic stability of the amorphous oxygen-rich and oxygen-poor compounds, which undergo composition phase separation. Including the previously neglected effects of this amorphous phase separation reconciles unexplained discrepancies in retention and enables predictive design of key performance indicators such as retention stability. Furthermore, this result emphasizes non-ideal thermodynamic interactions as key design criteria in post-digital devices with defect densities substantially exceeding those of today’s covalent semiconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Theory of terahertz pulse transmission through ferroelectric nanomembranes

An analytical model is developed to predict the temporal evolution of the lattice polarization in ferroelectric nanomembranes upon the excitation by a terahertz (THz) electromagnetic pulse of an arbitrary waveform and the concurrent transmission of the THz pulse in both linear and nonlinear regimes. It involves the use of the perturbation method to solve the equation of motion for the lattice polarization in both unclamped and strained ferroelectric nanomembranes within the framework of Landau-Ginzburg-Devonshire theory. The model is applicable to perovskite oxides such as BaTiO 3 and SrTiO 3 , wurtzite Al 1−𝑥 ⁢Sc 𝑥 ⁢N, and trigonal LiNbO 3 . Our analytical model provides a theoretical basis for determining the thermodynamic and kinetic parameters of ferroelectric materials through a THz transmission experiment. The calculation results also suggest an approach to reversing the chirality of a circularly polarized THz pulse by harnessing the resonant polarization-photon coupling in ferroelectrics. This capability of chirality reversal, along with the high tunability from a strain applied along any arbitrarily oriented in-plane axis, provides new opportunities for THz wave modulation without relying on complex metasurface designs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Multiscale and Machine Learning Modeling for Additive Manufacturing

Additive manufacturing (AM) techniques provide the opportunity to simultaneously design new materials and components with complex structures in less time, enabling faster material developments. Even though compositionally similar, the texture of the materials produced by such techniques is significantly different from conventionally manufactured materials. Additively manufactured materials produces highly heterogeneous microstructure within a single build. Such variations in the microstructure make qualifying AM products challenging for extreme environment applications. Understanding the AM process and its influence on the materials’ microstructures/properties is paramount for evaluating the workability and performance of the manufactured materials. The performance of AM materials for advanced nuclear reactor applications is of interest to the Advanced Materials and Manufacturing Technologies (AMMT) program under the Department of Energy Office of Nuclear Energy. Hence, considering the microstructural variabilities in the AM products and their impact on the performance of the material, it is important to correlate the process conditions to the final product and establish a process-structure-property- performance (PSPP) correlation for AM materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Influence of chemical strains on the electrocaloric response, polarization morphology, tetragonality, and negative-capacitance effect of ferroelectric core-shell nanorods and nanowires

Using Landau-Ginzburg-Devonshire (LGD) approach, we proposed the analytical description of the influence of chemical strains on spontaneous polarization and the electrocaloric response in ferroelectric core-shell nanorods. We postulate that the nanorod core presents a defect-free single-crystalline ferroelectric material, and elastic defects are accumulated in the ultrathin shell, where they can induce tensile or compressive chemical strains. Finite-element modeling (FEM) based on the LGD approach reveals transitions of domain-structure morphology induced by chemical strains in the Ba Ti O 3 nanorods. Namely, tensile chemical strains induce and support the single-domain state in the central part of the nanorod, while the curled domain structures appear near the unscreened or partially screened ends of the rod. The vortexlike domains propagate toward the central part of the rod and fill it entirely, when the rod is covered by a shell with compressive chemical strains above some critical value. The critical value depends on the nanorod sizes, aspect ratio, and screening conditions at its ends. Both analytical theory and FEM predict that the tensile chemical strains in the shell increase the nanorod polarization, lattice tetragonality, and electrocaloric response well above the values corresponding to the bulk material. The physical reason for the increase is strong electrostriction coupling between the mismatch-type elastic strains induced in the core by chemical strains in the shell. Comparison with earlier XRD data confirmed an increase of the tetragonality ratio in tensile Ba Ti O 3 nanorods compared to the bulk material. Obtained analytical expressions, which are suitable for the description of strain-induced changes in a wide range of multiaxial ferroelectric core-shell nanorods and nanowires, can be useful for strain engineering of advanced ferroelectric nanomaterials for energy storage, harvesting, electrocaloric applications, and negative capacitance elements. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2024

In‐Operando Spatiotemporal Imaging of Coupled Film‐Substrate Elastodynamics During an Insulator‐to‐Metal Transition

Abstract The drive toward non‐von Neumann device architectures has led to an intense focus on insulator‐to‐metal (IMT) and the converse metal‐to‐insulator (MIT) transitions. Studies of electric field‐driven IMT in the prototypical VO 2 thin‐film channel devices are largely focused on the electrical and elastic responses of the films, but the response of the corresponding TiO 2 substrate is often overlooked, since it is nominally expected to be electrically passive and elastically rigid. Here, in‐operando spatiotemporal imaging of the coupled elastodynamics using X‐ray diffraction microscopy of a VO 2 film channel device on TiO 2 substrate reveals two new surprises. First, the film channel bulges during the IMT, the opposite of the expected shrinking in the film undergoing IMT. Second, a microns thick proximal layer in the substrate also coherently bulges accompanying the IMT in the film, which is completely unexpected. Phase‐field simulations of coupled IMT, oxygen vacancy electronic dynamics, and electronic carrier diffusion incorporating thermal and strain effects suggest that the observed elastodynamics can be explained by the known naturally occurring oxygen vacancies that rapidly ionize (and deionize) in concert with the IMT (MIT). Fast electrical‐triggering of the IMT via ionizing defects and an active “IMT‐like” substrate layer are critical aspects to consider in device applications.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Defects, Entropy, and the Stabilization of Alternative Phase Boundary Orientations in Battery Electrode Particles

Using a novel statistical approach that efficiently explores the space of possible defect configurations, our present study investigates the chemomechanical coupling between interfacial structural defects and phase boundary alignments within phase-separating electrode particles. Applied to the battery cathode material Li X FePO 4 as an example, the theoretical analysis reveals that small, defect-induced deviations from an ideal interface can lead to dramatic shifts in the orientations of phase boundaries between Li-rich and Li-lean phases, stabilizing otherwise unfavorable orientations. Significantly, this stabilization arises predominantly from configurational entropic factors associated with the presence of the interfacial defects rather than from absolute energetic considerations. The specific entropic factors pertain to the diversity of defect configurations and their contributions to rotational/orientational rigidity of phase boundaries. Comparison of the predictions with experimental observations indicates that the additional entropy contributions indeed play a dominant role under actual cycling conditions, leading to the conclusion that interfacial defects must be considered when analyzing the stability and evolution kinetics of the internal phase microstructure of strongly phase-separating systems. Possible implications for tuning the kinetics of (de)lithiation based on selective defect incorporation are discussed. Ultimately, this understanding can be generalized to the chemomechanics of other defective solid phase boundaries.

25 ENERGY STORAGE↗
DOE OSTI2018

Local Exact Controllability of Two-Phase Field Solidification Systems with Few Controls

We analyze a control problem for a phase field system modeling the solidification process of materials that allow two different types of crystallization coupled to a Navier–Stokes system and a nonlinear heat equation, with a reduced number of controls. We prove that this system is locally exactly controllable to suitable trajectories, with controls acting only on the motion and heat equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Intermediate Principal Stress Effects on the 3D Cracking Behavior of Flawed Rocks Under True Triaxial Compression

Crack initiation, growth, and coalescence in flawed rocks have been extensively studied for 2D (planar, penetrating) flaws under uniaxial/biaxial compression. However, little is known as to the mechanisms and processes of cracking from 3D flaws under true triaxial compression, where the intermediate principal stress (σ 2 ) is distinguished from the major and minor principal stresses. In this work, we systematically investigate the effects of σ 2 on the 3D cracking behavior of rock specimens with preexisting flaws, through the use of mechanistic simulations of mixed-mode fracture in rocks. We explore how two characteristics of σ 2 , namely, (i) its orientation with respect to the flaw and (ii) its magnitude, affect two aspects of the cracking behavior, namely, (i) the cracking pattern and (ii) the peak stress. Results show that the orientation of σ 2 exerts more control over the cracking pattern than the flaw inclination angle. The peak stress becomes highest when σ 1 is parallel to the flaw, whereas it becomes lowest when σ 2 is parallel to the flaw. Also, the effects of σ 2 magnitude are more significant when σ 2 becomes more oblique to the flaw plane. On the basis of our observations, we propose mechanisms underlying the cracking behavior of 3D flawed rocks under true triaxial compression.

58 GEOSCIENCES↗
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 OSTI2018

Phase field simulation of martensitic transformation in pre-strained nanocomposite shape memory alloys

We show in this paper how strain engineering alters the fundamental characteristic of a martensitic transformation (MT) and gives it a new set of properties including large quasi-linear elastic strain response with nearly vanishing hysteresis and low elastic modulus. The work is motivated and inspired by a recent experimental study on elastic and inelastic (transformation) strain matching in a pre-strained nano-composite with Nb nanowires embedded in a NiTi shape memory alloy matrix. In particular, we demonstrate by computer simulation that dislocations at Nb/NiTi interfaces produced by the pre-straining are responsible for the unprecedented properties. Microstructural evolution captured in the simulations reveals that local stress fields associated with the dislocations regulate the nucleation and growth of martensite, turning the otherwise sharp, strong first-order transition into a continuous, high-order like transition. The simulations predict that the stress-strain hysteresis and modulus of the composite decrease with increasing amount of pre-strain, which agrees well with the experimental measurement. This study suggests a design strategy by introducing non-uniform stress fields for enhanced properties of shape memory alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Morphologies of dealloying corrosion attack at grain boundaries

Dealloying corrosion at grain boundaries severely compromises the performance of polycrystalline materials across a wide variety of technological applications. The impact of this phenomenon depends upon the morphology and rate of intergranular dealloying, which can range from planar to wormhole-like patterns that rapidly advance into the alloy. Using 2D and 3D multi-phase field simulations, we reveal how diverse microstructures result from a fundamental interplay between alloy composition and a grain boundary migration mechanism that alters diffusional pathways of dealloying. Inside intergranular dealloying channels, corrosion product buildup can spawn new channels that branch into grain interiors, such that alloys can be degraded from the inside out. These processes further lead to an atypical coarsening mode assisted by diffusion in the dealloying agent. We summarize a unifying explanation for dealloying morphology selection in polycrystalline alloys, which provides an important step towards their optimization for dealloying corrosion environments.

Corrosion↗