Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase field”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Phase II Field Demonstration at Lansing Smith Generating Plant, Southport, Florida (Final Report)

The Final Technical Report: Field Test Design and Pressure Management Strategies for Phase II Field Demonstration of Optimal Pressure Control, Plume Management, and Produced Water Strategies presents the culmination of multi-year efforts under the U.S. Department of Energy’s Brine Extraction and Storage Test (BEST) program, led by EPRI in partnership with Gulf Power at the Plant Smith site near Panama City, Florida. The project was designed to evaluate and demonstrate the technical feasibility of managing subsurface pressures and fluid movement associated with large-scale CO₂ injection, using low-salinity water as a proxy fluid. Through a combination of field injection testing, reservoir modeling, and optimization studies, the research team developed and refined active and passive brine extraction strategies aimed at controlling injection-induced pressure buildup, mitigating risks of fault activation, and managing plume migration. The field demonstration incorporated a new injection well (TIW-2), a new monitoring/extraction well (TEMW-A), and an existing well (TIW-1) repurposed for passive pressure relief. Complementary geophysical monitoring was designed to track plume development and evaluate the effectiveness of pressure management strategies. The report details the integrated workflow encompassing wellfield development, regulatory permitting, model calibration, and survey design. It includes a comprehensive electromagnetic (EM) modeling and inversion study used to develop a cost-effective, time-lapse geophysical monitoring plan capable of imaging the evolving low-salinity plume within the highly saline Lower Tuscaloosa formation. Reservoir simulation results guided the design of the 17-month injection program and the timing of active extraction to maintain formation pressures below a hypothetical fault reactivation threshold. Supporting analyses evaluated potential injectivity risks related to clay fines migration and geochemical incompatibility, identifying practical mitigation measures such as salinity and pH control. The final design recommends a series of crosswell EM surveys—conducted before, during, and after injection—as the most effective and economical approach for plume imaging, supplemented by continuous downhole pressure and flow monitoring. Collectively, the study provides a field-ready framework for cost-effective pressure management and monitoring in support of future CO₂ storage projects.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2020

Magnetic field induced phase transition in spinel GeNi 2 O 4

Cubic spinel GeNi 2 O 4 exhibits intriguing magnetic properties with two successive antiferromagnetic phase transitions (T N1 ≈ 12.1 and T N2 ≈ 11.4 K) with the absence of any structural transition. We have performed detailed heat capacity and magnetic measurements in different crystallographic orientations. A new magnetic phase in presence of magnetic field (H ≥ 4 T) along the [111] direction is revealed, which is not observed when the magnetic field is applied along the [100] and [110] directions. High field neutron powder diffraction measurements confirm such a change in magnetic phase, which could be ascribed to a spin reorientation in the presence of magnetic field. In conclusion, a strong magnetic anisotropy and competing magnetic interactions play a crucial role on the complex magnetic behavior in this cubic system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Multiscale modeling for high burnup structure formation and associated pulverization

This report summarizes the lower length scale modeling work performed in the fiscal year 2023 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to capture the microstrutural evolution and associated pulverization criteria in the high burnup regions. This year, the resolution mechanisms within the cluster dynamics code has been updated and it influence on the bubble growth in HBS regions at the mesoscale is studied. To improve the accuracy of phase-field models of fission gas bubble growth, a new Helmholtz free energy for high-density Xe gas is derived from a virial equation of state. Two previously used strategies for representing net vacancy production in phase-field models of fission gas bubble growth are compared with each other and with analytical models of bubble growth. The vacancy source only model is found to be more convenient to parameterize realistically compared with the vacancy source+sink model. Furthermore, the phase-field-fracture simulation have been performed using MD- informed failure stress values and realistic HBS structure obtained from the phase-field simulations. We also present the uncertainty bands on prediction of the critical stress to account for the effect of the lower length scale variabilities on the failure criteria at the mesoscale. It is observed that pulverization may occur in partially restructured regions with restructuring fraction as low as 17%. In light of this, an update to the BISON’s pulverization criteria is recommended.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Mesoscale Modeling of the Effects of Accelerated Burnup on UO2 Microstructural Evolution

Accelerating the nuclear fuel qualification process will rely on some combination of advanced modeling and simulation techniques with accelerated irradiation testing and separate effects experiments to enable the development of new fuel concepts in a shorter time frame. One of the key challenges to successfully leveraging accelerated irradiation tests will be understanding the artifacts that may be introduced with accelerated accumulation of dose and/or burnup. This work presents phase field (MARMOT) simulations of the evolution of representative 2D UO2 microstructures up to 40 MWd/kgU. Simulations were performed under both commercial light water reactor fuel conditions as well as those that would be expected for highly accelerated (~10x) burnup conditions similar to those used in the MiniFuel irradiations in Oak Ridge National Laboratory’s High Flux Isotope Reactor. The phase field model was coupled with a discrete nucleation algorithm to model re- structuring at high burnup. The effect of the different fission rates in both microstructures was investigated at two temperatures: 650?C and 800?C. The lower temperature simulations both showed an onset of restructuring at nearly 60 MWd/kgU. More extensive restructuring was obtained in the MiniFuel microstructure compared with that of the PWR fuel. At 800?C, no restructuring was obtained as a result of the thermally activated diffusion of Xe atoms and U vacancies to fission gas bubbles, which reduces the nucleation driving force. These results highlight the importance of using modeling and simulation tools to inform the environmental conditions during targeted accelerated irradiation tests to extract the most useful fuel performance data.

accelerated fuel qualification, Phase Field, Restr↗
DOE OSTI2017

Simulations of irradiated-enhanced segregation and phase separation in Fe–Cu–Mn alloys

For reactor pressure vessel steels, the addition of Cu, Mn, and Ni has a positive effect on mechanical, corrosion and radiation resistance properties. However, experiments show that radiation-enhanced segregation and/or phase separation is one of important material property degradation processes. In this work, we developed a model integrating rate theory and phase-field approaches to investigate the effect of irradiation on solute segregation and phase separation. The rate theory is used to describe the accumulation and clustering of radiation defects while the phase-field approach describes the effect of radiation defects on phase stability and microstructure evolution. The Fe-Cu-Mn ternary alloy is taken as a model system. The free energies used in the phase-field model are from CALPHAD. Spatial dependent radiation damage from atomistic simulations is introduced into the simulation cell for a given radiation dose rate. The radiation effect on segregation and phase separation is taken into account through the defect concentration dependence of solute mobility. With the model the effect of temperatures and radiation rates on Cu and Mn segregation and Cu-rich phase nucleation are systematically investigated. The segregation and nucleation mechanisms are analyzed. The simulations demonstrated that the nucleus of Cu precipitates has a core-shell composition profile, i.e., Cu rich at center and Mn rich at the interface, in good agreement with the theoretical calculation as well as experimental observations.

Li, Boyan↗
DOE OSTI2014

Field-current phase diagrams of in-plane spin transfer torque memory cells with low effective magnetization storage layers

Field-current phase diagrams were measured on in-plane anisotropy Co{sub 60}Fe{sub 20}B{sub 20} magnetic tunnel junctions to obtain the spin transfer torque (STT) field-current switching window. These measurements were used to characterise junctions with varying free layer thicknesses from 2.5 down to 1.1 nm having a reduced effective demagnetizing field due to the perpendicular magnetic anisotropy at CoFeB/MgO interface. Diagrams were obtained with 100 ns current pulses, of either same or alternating polarity. When consecutive pulses have the same polarity, it is possible to realize the STT switching even for conditions having a low switching probability. This was evidenced in diagrams with consecutive pulses of alternating polarity, with 100% switching obtained at 4.7 MA/cm{sup 2}, compared to the lower 3.4 MA/cm{sup 2} value for same polarity pulses. Although the low level of the current density window is higher in alternating polarity diagrams, the field window in both diagrams is the same and therefore independent of the pulse polarity sequence.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

A chain stretch-based gradient-enhanced model for damage and fracture in elastomers

Similar to quasi-brittle materials, it has been recently shown that elastomers can exhibit a macroscopically diffuse damage zone that accompanies the fracture process. In this study, we introduce a stretch-based gradient-enhanced damage (GED) model that allows the fracture to localize and also captures the development of a physically diffuse damage zone. This capability contrasts with the paradigm of the phase field method for fracture, where a sharp crack is numerically approximated in a diffuse manner. Capturing fracture localization and diffuse damage in our approach is achieved by considering nonlocal effects that encompass network topology, heterogeneity, and imperfections. These considerations motivate the use of a statistical damage function dependent upon the nonlocal deformation state. From this model, fracture toughness is realized as an output. While GED models have been classically utilized for damage modeling of structural engineering materials (e.g., concrete), they face challenges when trying to capture the cascade from damage to fracture, often leading to damage zone broadening (de Borst and Verhoosel, 2016). This deficiency contributed to the popularity of the phase-field method over the GED model for elastomers and other quasi-brittle materials. Other groups have proceeded with damage-based GED formulations that prove identical to the phase-field method (Lorentz et al., 2012), but these inherit the aforementioned limitations. To address this issue in a thermodynamically consistent framework, we implement two modeling features (a nonlocal driving force bound and a simple relaxation function) specifically designed to capture the evolution of a physically meaningful damage field and the simultaneous localization of fracture, thereby overcoming a longstanding obstacle in the development of these nonlocal strain- or stretch-based approaches. Here, we discuss several numerical examples to understand the features of the approach at the limit of incompressibility, and compare them to the phase-field method as a benchmark for the macroscopic response and fracture energy predictions.

Elastomers↗
DOE OSTI2023

Discovery of Nanoscale Electric Field‐Induced Phase Transitions in ZrO 2

Abstract The emergence of ferroelectric and antiferroelectric properties in the semiconductor industry's most prominent high‐k dielectrics, HfO 2 and ZrO 2 , is leading to technology developments unanticipated a decade ago. Yet the failure to clearly distinguish ferroelectric from antiferroelectric behavior is impeding progress. Band‐excitation piezoresponse force microscopy and molecular dynamics are used to elucidate the nanoscale electric field‐induced phase transitions present in ZrO 2 ‐based antiferroelectrics. Antiferroelectric ZrO 2 is clearly distinguished from a closely resembling pinched La‐doped HfO 2 ferroelectric. Crystalline grains in the range of 3 – 20 nm are imaged independently undergoing reversible electric field induced phase transitions. The electrically accessible nanoscale phase transitions discovered in this study open up an unprecedented paradigm for the development of new nanoelectronic devices.

36 MATERIALS SCIENCE↗
NASA NTRS1983

Sintering, properties and fabrication of Si3N4 + Y2O3 based ceramics

Pure silicon nitride shows a remarkable resistance to sintering without the use of densification additives. The present investigation is concerned with results which show the effect of chemical content on sinterability, taking into account the composition, raw material impurities, and processing contaminants. Aspects of sintering are discussed along with strength characteristics, and oxidation relations. Attention is given to phase field I and II materials, phase field III and IV materials, tungsten carbide and oxidation at 600 C, and studies involving shape fabrication by injection molding. It was found that in sintering Si3N4 + Y2O3 an increase in the amount of Y2O3 and, in particular, the addition of Al2O3 enhances the fluidity of the liquid phase.

Quackenbush, C. L.↗
DOE OSTI2017

Formation of gapless Z{sub 2} spin liquid phase manganites in the (Sm{sub 1–y}Gd{sub y}){sub 0.55}Sr{sub 0.45}MnO{sub 3} system in zero magnetic field: Topological phase transitions to states with low and high density of 2D-vortex pairs induced by the magnetic field

The evolution of the ground state of the manganese spin ensemble in the (Sm{sub 1–y}Gd{sub y}){sub 0.55}Sr{sub 0.45}MnO{sub 3} in the case of isovalent substitution of rare-earth samarium ions with large radii with gadolinium ions with significantly smaller radii is studied. The measured temperature dependences of the ac magnetic susceptibility and the field dependences of the dc magnetizations are analyzed using the Heisenberg–Kitaev model describing the transition from the ordered spin state with classical isotropic AFM exchange to the frustrated spin state with quantum highly anisotropic FM exchange. A continuous transition from the 3D ferromagnetic state of manganese spins in the initial sample with y = 0 to zigzag AFM ordering of CE-type spins in ab planes for y = 0.5, coexisting in samples with y = 0.5, 0.6, and 0.7 at temperatures below T{sub N} ≅ 48.5 K with a disordered phase such as a quantum Griffiths phase is identified. As the gadolinium concentration further increases, the CE-type zigzag AFM structure is molten, which leads to the appearance of an unusual phase in Gd{sub 0.55}Sr{sub 0.45}MnO{sub 3} in the temperature range close to the absolute zero. This phase has characteristic features of a gapless Z{sub 2} quantum spin liquid in zero external magnetic field. The step changes in the magnetization isotherms measured at 4.2 K in the field range of ±75 kOe are explained by quantum phase transitions of the Z{sub 2} spin liquid to a phase with topological order in weak magnetic fields and a polarized phase in strong fields. The significant difference between critical fields and magnetization jumps in isotherms indicates the existence of hysteretic phenomena in quantum spin liquid magnetization–demagnetization processes caused by the difference between localization–delocalization of 2D vortex pairs induced by a magnetic field in a quantum spin liquid with disorder.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Meso-scale modeling of UO2 nuclear fuel to high burnup

To improve the economics of light water reactors for commercial nuclear energy generation, utility operators are seeking to obtain regulatory approval to run UO2 fuel to higher levels of burnup. One potential impediment to obtaining this approval is the phenomenon of fuel fragmentation, relocation, and dispersal (FFRD). FFRD can result when fuel experiences a rapid temperature transient, such as that occurring during a Loss Of Coolant Accident (LOCA). FFRD has historically been most associated with the rim region in UO2 fuel pellets, where the phenomenon of fragmentation is also referred to as pulverization due to the small size of the fragments. More recent evidence suggests that the so-called “dark zone” (due to its appearance in micrographs) that can be observed in the mid-radial regions of high burnup fuel is also susceptible to FFRD. Although empirical fuel performance models have been developed that can adequately predict pulverization in the rim region under typical LWR conditions, a scientific understanding of what underlies fuel restructuring and subsequent FFRD is lacking even in the rim region, and no models are currently available for the behavior the dark zone. To address these challenges, the U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has employed a multi-scale modeling approach to improve scientific understanding and develop new fuel performance models. In this talk, I will focus on meso-scale efforts, which form a crucial link between atomic-scale and engineering-scale models. Phase-field modeling combined with cluster dynamics is used to predict the restructuring process in the rim region. Phase-field fracture modeling, informed by atomistic simulations, is used to predict the onset of pulverization in the rim region. Combining these techniques together allows the extent of rim pulverization to be predicted. The formation and evolution of the dark zone has also been simulated with the phase-field method, using an improved approach to vacancy source term parameterization. The work shows the important impact of microstructure on fuel performance.

fracture↗
DOE OSTI2020

Dendrite Growth Morphology Modeling in Liquid and Solid Electrolytes

The main goal of this project is to develop a multi-scale modeling approach that connects micron-scale phase-field models and atomic-scale density functional theory (DFT)-based simulations via parameter- and relationship-passing in order to predict Li-metal dendrite morphology evolution, in both liquid and solid electrolytes. The key hypothesis of the DFT-informed phase-field multiscale modeling approach is that it can capture the electrochemical-mechanical driving forces and incorporate the roles of nano-meter-thin solid electrolyte interphase (SEI) in liquid electrolytes as well as of the microstructures of micro-meter-thick solid electrolytes (SEs) for all-solid-state batteries. In this project, we have formulated and implemented phase-field models to incorporate the electrochemical driving forces in liquid electrolytes and then incorporate mechanical driving forces to simulate dendrite growth in solid electrolytes with resolved microstructures. We have implemented two treatments for the SEI: an explicit model to include the microstructure of the SE or SEI in the phase field model and an implicit model to simulate the impact of nano-meter thick SEI in liquid electrolytes by varying the electrode/electrolyte interfacial properties. The key interfacial properties, including the electronic and ionic transport properties, the charge transfer reaction kinetics, and mechanical properties, were computed by DFT-based calculations. At the DFT-based model, one key advancement is to directly predict the charge transfer reaction kinetics at a complex Li/SEI/electrolyte interface by linking DFT with density functional tight binding (DFTB) calculations. As the main accomplishments, we have demonstrated two successful predictions in both solid electrolyte and liquid electrolyte based on this multiscale approach. The predicted intergranular Li dendrite growth in LLZO revealed the importance of trapped electrons at internal interfaces in the microstructure of LLZO. The predicted electroplating morphology of mossy Li and faceted Mg agreed well with experiments. The insights provided by the multiscale model and the model enabled electrolyte and SEI design will accelerate the development of Li-metal electrode for high energy density batteries, that meet DOE’s target on cell density (>350 Wh/kg) and cost below $100/kWhuse for EV applications.

25 ENERGY STORAGE↗
DOE OSTI2020

Computational thermomechanics for crystalline rock. Part II: Chemo-damage-plasticity and healing in strongly anisotropic polycrystals

We present a thermal–mechanical–chemical-phase field model that captures the multi-physical coupling effects of precipitation creeping, crystal plasticity, anisotropic fracture, and crack healing in polycrystalline rock at various temperature and strain-rate regimes. This model is solved via a fast Fourier transfer solver with an operator-split algorithm to update displacement, temperature and phase field, and chemical concentration incrementally. In nuclear waste disposal in salt formation, brine inside the crystal salt may migrate along the grain boundary and cracks due to the gradient of interfacial energy and pressure. This migration has a significant implication on the permeability evolution, creep deformation, and crack healing within rock salt but is difficult to incorporate implicitly via effective medium theories compared with computational homogenization. As such, we introduce a thermodynamic framework and a corresponding computational implementation that explicitly captures the brine diffusion along the grain boundary and crack at the grain scale. Meanwhile, the anisotropic fracture and healing are captured via a high-order phase field that represents the regularized crack region in which a newly derived non-monotonic driving force is used to capture the fracture and healing due to the solution–precipitation. Numerical examples are presented to demonstrate the capacity of the thermodynamic framework to capture the multiphysics material behaviors of rock salt.

42 ENGINEERING↗
DOE OSTI2022

GrainNN: A neighbor-aware long short-term memory network for predicting microstructure evolution during polycrystalline grain formation

High fidelity simulations of grain formation in alloys are an indispensable tool for process-to-mechanical-properties characterization. Such simulations, however, can be computationally expensive as they require fine spatial and temporal discretizations. Their cost becomes an obstacle to parametric studies and ensemble runs and ultimately makes downstream tasks like optimal control and uncertainty quantification challenging. To enable such downstream tasks, we introduce GrainNN, an efficient and accurate reduced-order model for epitaxial grain growth in additive manufacturing conditions. GrainNN is a sequence-to-sequence long-short-term-memory (LSTM) deep neural network that evolves the dynamics of manually crafted features. Its innovations are (1) an attention mechanism with grain-microstructure-specific transformer architecture; and (2) an overlapping combination of several clones of the network to generalize to grain configurations that are different from those used for training. This design enables GrainNN to predict grain formation for unseen physical parameters, grain number, domain size and geometry. Furthermore, GrainNN not only reconstructs the quantities of interest but also can be pointwise accurate. In our numerical experiments, we use a polycrystalline phase field method to both generate the training data and assess GrainNN. For multiparametric, ensemble simulations with many grains, GrainNN can be orders of magnitude faster than phase field simulations, while delivering 5%–15% pointwise error. Additionally, this speedup includes the cost of the phase field simulations for generating training data.

36 MATERIALS SCIENCE↗
NASA NTRS1982

Superalloy composition modeling

Development of a predictive method for determination of the gamma/gamma prime phase fields, i.e., gamma prime volume fraction as a function of the multicomponent composition, is described. The cluster variation method used for binary alloys in which the precipitated phase is coherent with the matrix phase is extended for application to the multicomponent coherent gamma/gamma prime nickel-based superalloys. It is shown that the cluster variation method can accurately describe the equilibrium (incoherent) gamma/gamma prime phase fields in the binary Ni-Al phase diagram. The gamma/gamma prime phase field for the Ni-Cr-Al ternary phase diagram is computed as a function of temperature. A reasonable fit results between the calculated and the experimental diagrams. The modeling of the six-component Ni-Cr-Al-Co-Mo-Ti base superalloy and the effect of Ni substitution of Co are discussed.

Barefoot, J.↗
DOE OSTI2014

The scaling of weak field phase-only control in Markovian dynamics

We consider population transfer in open quantum systems, which are described by quantum dynamical semigroups (QDS). Using second order perturbation theory of the Lindblad equation, we show that it depends on a weak external field only through the field's autocorrelation function, which is phase independent. Therefore, for leading order in perturbation, QDS cannot support dependence of the population transfer on the phase properties of weak fields. We examine an example of weak-field phase-dependent population transfer, and show that the phase-dependence comes from the next order in the perturbation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗