Search NASA⌕ Search

SEARCH · Search NASA

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.

1,404 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Integrated Simulation of Weld Residual Stress Evolution and Crack Propagation Using XFEM

Nuclear power plant components operate in environments that promote multiple degradation mecha- nisms, several of which involve crack initiation and growth. An ongoing effort in the U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program is developing a general capability within the Multiphysics Object Oriented Simulation Environment (MOOSE) framework for simulating three-dimensional crack growth under a range of driving conditions, including fatigue, stress corrosion cracking (SCC), brittle fracture, and stress-relaxation cracking. This report demonstrates an end-to-end workflow that uses this capability to model weld-residual-stress-driven SCC in the J-groove weld of a pressurized-water reactor control rod drive mechanism penetration in the vessel head. The workflow consists of a thermomechanical welding simulation with temperature-dependent plasticity, followed by cooldown to ambient conditions, and a restart of the simulation using the MOOSE extended finite element method (XFEM) module to propagate a three-dimensional crack through the residual stress field. New welding capabilities were developed to properly initialize newly activated elements in the weld region, and robustness improvements were made to the mesh-based algorithm for defining cutting planes in the 3D XFEM algorithm, allowing it to handle complex crack fronts and stress fields. Together these advances allowed the simulated SCC crack to grow from an initial elliptical flaw in the weld, across the weld, through the tube wall, and almost to the triple point (where the weld, tube, and reactor pressure vessel head intersect) over roughly 36 years of simulated service. These results demonstrate a workflow that can be extended to fully three-dimensional welding simulations and more complex crack interaction problems.

42 - ENGINEERING↗
DOE OSTI2026

Field Validation of a Grid-Interactive Efficient Building Software Solution

The U.S. General Services Administration's (GSA's) Green Proving Ground (GPG) program, in partnership with the National Laboratory of the Rockies (NLR), completed a field study of a Grid-Interactive Efficient Buildings (GEB) software solution. The study focused on a single testbed facility to test the GEB functionality of the software solution, along with other features. The testbed facility - a courthouse - is a common building type in GSA's vast building portfolio, offering potentially impactful findings on a scalable level. The study evaluated Prescriptive Data's technology, Nantum OS, a connected building operating system ("GEB Solution") which aggregates multiple sources of previously siloed building data and combines that data with external sources, such as weather information or utility signals, into a single integrated platform. A GEB Solution is a type of Energy Management Information System (EMIS). EMIS is defined as a system of devices, data services, and software applications that communicates with any building system or third-party data source to aggregate and transform data into new capabilities to aid in the optimization of energy use at the building, campus, or agency level. This specific GEB Solution is an EMIS with ASO, automated system optimization, offering supervisory control of certain aspects of the Building Automation System (BAS). Multiple features were evaluated including, but not limited to, Continuous Demand Management to avoid setting new monthly kilowatt (kW) peaks, energy efficiency for reduction of kilowatt hours (kWh) and natural gas consumption, and automated demand response (ADR) for purposes of lowering demand during a utility called Demand Response (DR) event. The testbed facility was the Foley Federal Building and US Courthouse ("Foley Federal Building") located in Las Vegas, NV. This is a 209,496 sq. ft. building constructed in the 1960s with major renovations in 2004. The facility was a good candidate due to the large prevalence of office and courthouse spaces in the GSA portfolio of buildings. It also has many features which allow integration into and control of the building and a strong facilities team to assist with the study. Quantitative and qualitative performance objectives were developed using GSA's GPG GEB project template along with input from the vendor and building facility staff; these are outlined in Table 1. The quantitative performance objectives focused on continuous demand management, energy efficiency, and automated demand response. The qualitative performance objectives focused on the ease of installation and commissioning as well as the operability of the GEB solution. Other performance metrics that are reported on include carbon reduction, cost effectiveness, and occupant acceptance.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS2022

Optical properties of morphologically complex black carbon aerosols: Effects of coatings

Optical properties are computed for fractal-like aggregate black carbon (BC) aerosols coated with different substances. Two models are used for these aerosols: (i) the coated aggregate model (model I), where coating is added to the voids and surroundings of BC fractal-like aggregates; and (ii) the closed-cell aggregate model (model II), where coating is added concentrically to each monomer of BC fractal-like aggregates. Our results favor choosing the coated aggregate model (model I) to simulate scattering and absorption by coated BC aerosols because this model is morphologically more realistic, and because this model yields mass absorption cross section (MAC) and backscattering linear depolarization ratio (LDR) values that are consistent with field measurements. Moreover the corresponding computed degree of linear polarization (DoLP) and LDR values are very sensitive to changes in the coating volume fraction (fvol) and the coating refractive index (m). With the same absorbing BC core, the MAC value increases steadily with increasing fvol and m values. For example, using a density of 1.8 g/cm 3 and BC refractive index 1.95 + i0.79, the calculated MAC values for uncoated BC aggregates range from 6.2 to 6.8 m 2 /g at 0.55 µm. When coating material is applied to the BC aggregates, the calculated MAC values for model I particles increase to between 9.8 and 13.2 m 2 /g (depending on m) when fvol = 87.5%. The backscattering LDR values also tend to increase with the increasing m values for the shapes and sizes considered in our study. For model I particles, the backscattering LDR values span a wide range of 4.2–27.8% at a wavelength of 0.35 µm at fvol = 87.5% when m increases from 1.33 to 1.55. Our results are relevant to analyses of polarimetric and lidar observations of smoke particles, especially when these particles undergo hygroscopic growth.

Black carbon aerosols↗
NASA NTRS1993

Charge-exchange erosion studies of accelerator grids in ion thrusters

A particle simulation model is developed to study the charge-exchange grid erosion in ion thrusters for both ground-based and space-based operations. Because the neutral gas downstream from the accelerator grid is different for space and ground operation conditions, the charge-exchange erosion processes are also different. Based on an assumption of now electric potential hill downstream from the ion thruster, the calculations show that the accelerator grid erosion rate for space-based operating conditions should be significantly less than experimentally observed erosion rates from the ground-based tests conducted at NASA Lewis Research Center (LeRC) and NASA Jet Propulsion Laboratory (JPL). To resolve this erosion issue completely, we believe that it is necessary to accurately measure the entire electric potential field downstream from the thruster.

Xiaohang Peng↗
NASA NTRSDate not supplied

Helmholtz Motor: A Novel Electric Machine That Enables High Temperature Superconducting Armatures

Electrified aircraft are being developed to increase the efficiency and reduce the cost of operating subsonic transport aircraft. Achieving a substantial impact necessitates focusing on single- and twin-aisle aircraft which use propulsion systems with >20 MW ratings. For these aircraft, multi-MW superconducting electric machines are being developed due to their high specific power and high efficiency. Cryogenic electric machines are a promising technology for multi megawatt electric aircraft drivetrains. Second generation high temperature superconductors (HTS) are the ideal conductor for the field winding (typically the rotor) of superconducting machines, and they have attractive features for the armature winding (typically the stator) of these machines. However, their use in armature windings has been severely limited due to excessive AC losses in conventional electric machines. This presentation will present the working principle, design, and a performance trade study of a novel electric machine – termed a Helmholtz machine – that enables HTS armatures by significantly reducing AC losses. The working principle of the patent pending Helmholtz motor is to optimize the rotor-produced magnetic field (i.e., the field winding) so that it produces a magnetic field almost purely in the plane of the HTS. This is accomplished using 2 rotors that contain a set of matched Helmholtz coil pairs (or matched permanent magnet Halbach arrays). The armature (stator) HTS coils are positioned between the rotors and oriented to minimize the out of plane magnetic field. The armature HTS coils are also designed to minimize the out of plane component of their self field. The presentation will include the lessons learned in the sizing and design of this motor, including equations for how to match the magnet arrays. A performance trade study will be presented for a partially superconducting, radial flux version of the motor. The trade study involved an optimization of a motor design code based on analytical calculations (electromagnetic, AC loss, thermofluid, mechanical, and motor sizing) that determine the motor’s total efficiency and total specific power. The design code will be described in the presentation, with an emphasis on the calculation of magnetic fields, AC loss, and temperature distribution.

Superconducting electric machines↗
DOE OSTI2025

Strong anharmonicity dictates ultralow thermal conductivities of type-I clathrates

Type-I clathrate solids have attracted significant interest due to their ultralow thermal conductivities and sub- sequent promise for thermoelectric applications, yet the mechanisms underlying these properties are not well understood. Here, we extend the framework of vibrational dynamical mean-field theory (VDMFT) to calculate temperature-dependent thermal transport properties of solids using a many-body Green’s function approach. When applied to a coarse-grained description of 𝑋 8 Ga 16 Ge 30 , where 𝑋= Ba, Sr, we find that nonresonant scattering between cage acoustic modes and rattling modes leads to a reduction of acoustic phonon lifetimes and thus thermal conductivities. Moreover, we find that the moderate temperature dependence of conductivities above 300 K, which is consistent with experimental measurements, cannot be reproduced by textbook perturbation theory calculations, which predict a 𝑇 −1 dependence. Therefore, we suggest that nonperturbative anharmonic effects, including four- and higher-phonon scattering processes, are responsible for the ultralow thermal conductivities of type-I clathrates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2019

Rewet Temperature Correlations for Liquid-Nitrogen Boiling Pipe Flows Across Varying Flow Conditions and Orientations

In many convective liquid-vapor phase change heat transfer engineering applications, cryogenic fluids are widely used in industrial processes, spacecraft and cryosurgery systems, and so on. For example, cryogens are usually used as liquid fuels such as liquid hydrogen and oxygen in the rocket industry, liquid nitrogen (LN2) and helium are frequently used to cool superconducting magnetic device for medical applications. In these systems, proper transport, handling, and storage of cryogenic fluids are of extreme importance. Among all the cryogenic transport processes performed in room temperatures, quenching, also termed chilldown, is a unavoidable initial, transient phase-change heat transfer process that brings the system down to the cryogenic condition. The Leidenfrost temperature or rewet temperature that signals the end of film boiling is practically considered the completion point of a quenching process. Therefore, rewet temperature has been considered the most important parameter for the engineering design of cryogenic thermal management systems. As most of the previous correlations for predicting the Leidenfrost temperature and the rewet temperature have been basically developed for water, they are shown to disagree with recent liquid nitrogen pipe chilldown experiments in upward and downward flow directions over a wide range of flow rates, pressures, and degrees of inlet subcooling. In addition to a complete review of the literature, two new correlations are presented in this work, one based on bubble growth and another based on the theoretical maximum limit of superheat. Each correlation performs well over the entire data set.

liquid nitrogen↗
NASA NTRS2026

Thermochemical Stability of Ca2Yb8(SiO4)6O2 Apatite in Presence of Molten Calcium-Magnesium-Aluminosilicate (CMAS)

Thermochemical stability of ytterbium silicon oxyapatite Ca 2 Yb 8 (SiO 4 ) 6 O 2 (CYbS) in the presence of molten calcium-magnesium aluminosilicate (CMAS) has been investigated at elevated temperatures for consideration as a thermal and environmental barrier coating (T/EBC) material. CYbS apatite powder was synthesized from the constituent oxides via a solid-state reaction method. Hot-pressed apatite substrates were exposed to molten CMAS at 1200, 1300, and 1400 °C for 1, 10, and 50 h. Development of phases in the interaction region of the heat-treated specimens was monitored using scanning electron microscopy, transmission electron microscopy, high-angle annular dark-field imaging, selected area electron diffraction, and energy dispersive x-ray spectroscopy. Monoclinic cyclosilicate Ca 3 Yb 2 (Si 3 O 9 ) 2 formed from interaction of CYbS apatite with CaO in the CMAS melt at the apatite-CMAS reaction front and continued to nucleate and grow within the residual CMAS in diffusion couples annealed for 1 to 50 h at 1200 °C as well as in those heat treated at 1300 °C for 1 h. Residual CMAS was depleted of Ca when cyclosilicate was present. Dendritic wollastonite (CaSiO3) was observed within the residual CMAS in couples annealed at 1200 and 1300 °C. Ingress of molten CMAS, because of its exponential decrease in viscosity, occurred through open pores and along the grain boundaries of the apatite substrates without any detectable chemical reaction at 1300 and 1400 °C. Results of this study indicate that Ca 2 Yb 8 (SiO 4 ) 6 O 2 apatite has the potential to mitigate the CMAS corrosion up to about 1200 °C but not at higher temperatures.

Narottam P Bansal↗
DOE OSTI2025

Spin-strain interactions under hydrostatic pressure in α-RuCl 3

We investigate the effects of hydrostatic pressure on 𝛼−RuCl 3 , a prototypical material for the Kitaev spin model on a honeycomb lattice with a possible spin-liquid ground state. Using ultrasound measurements at pressures up to 1.16 GPa, we reveal significant modifications of the acoustic properties and the 𝐻−𝑇 phase diagram of this material. Hydrostatic pressure suppresses the three-dimensional magnetic order and induces a dimerization transition at higher pressures. At low pressures, the sound attenuation exhibits a linear temperature dependence, while above 0.28 GPa, it becomes nearly temperature independent, suggesting a shift in the phonon scattering regime dominated by Majorana fermions. These findings provide new insights into spin-strain interactions in Kitaev magnets and deliver a detailed characterization of the 𝐻−𝑇 phase diagram of 𝛼−RuCl 3 under hydrostatic pressure.

Focused ion beam↗
DOE OSTI2026

Resonant scattering at the center of the galaxy cluster PKS 0745-191 with XRISM

We report evidence of the resonant scattering effect at the center of the galaxy cluster PKS 0745-191 with XRISM. We analyzed XRISM/Resolve commissioning-phase observations of the distant cluster PKS 0745-191 ($z = 0.103$) with a $54$ ks exposure. The gain drift was corrected using the onboard modulated X-ray source (MXS), and spectra were extracted from all the pixels that were well illuminated by MXS, the core region (four central pixels, $\sim$100 kpc), and the surrounding region. A single-temperature collisional ionization equilibrium (CIE) model fits the full field-of-view spectrum with $kT \approx 6$ keV and a turbulent velocity of ${\approx}120$ km s$^{-1}$. From the core ($r < 50$ kpc) spectrum, we detected a ${\approx}22\%$ suppression of the Fe xxv He$\alpha$ resonance (w) line relative to the CIE prediction. We performed a Monte Carlo simulation to calculate the resonant scattering (RS) effect using radial profiles from Chandra data. The RS-inferred turbulence agrees with that by Resolve line-broadening, demonstrating that RS provides an independent and consistent constraint on ICM turbulence. These results highlight XRISM/Resolve’s potential for turbulence studies in galaxy clusters.

Astronomy and AstroPhysics↗
DOE OSTI2026

Role of Amorphous Chains in Nanoplastic Formation from Semicrystalline Polymers

Semicrystalline polymers, e.g., polyethylene terephthalate (PET), release micro- (100 nm to 1 mm) and nanoplastics (up to 100 nm) [MNPL] when they are degraded under quiescent conditions. However, the exact molecular mechanisms leading to material fragmentation into MNPL are unknown. Here, we monitor the evolution of chain molecular weight and the MNPL production kinetics during hydrolysis of PET pellets. We find that only ~0.6% of the amorphous phase ester bonds are hydrolyzed at the onset of MNPL release. Then, by combining a random scission model with measured amorphous spacings, we estimate that only ~15% of the stress transmitters in the amorphous phase, namely, bridges and entangled loops, have failed by this point. Thus, spontaneous fragmentation of the semicrystalline nanostructure occurs despite significant intercrystalline connectivity. We resolve this apparent contradiction by proposing that fragmentation can only occur when progressive tie-chain scission causes the material to undergo the ductile/brittle transition. Mechanistically, we speculate that the internal stresses responsible for material fragmentation are caused by processing (i.e., residual stresses) and/or sample densification induced by chemi-crystallization. We discuss additional factors that may affect our analysis and thus require further investigation, such as skin-core effects, preferential degradation of stress transmitters and/or recrystallization processes.

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

Thermal Bekenstein-Hawking entropy from the worldsheet

We define and compute the leading sphere diagram contribution to the entropy of the BTZ black hole supported by Kalb-Ramond flux in bosonic string theory. In a winding condensate description, integrating exactly over the constant mode for the radial direction of AdS 3 reduces the problem to one of the correlation functions of winding operators in the free theory. The volume of the residual PSL(2,ℂ) gauge group of the sphere is canceled by the action of conformal transformations on the winding interaction insertions. We formulate a precise version of the replica trick in terms of (infinitesimally) non-integer winding condensates to produce the entropy of the BTZ black hole. The resulting entropy can be calculated from the one-point function of a non-local operator on the worldsheet.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

The brighter-fatter effect in the JWST MIRI Si:As IBC detectors: I. Observations, impact on science, and modeling

Context.The Mid-Infrared Instrument (MIRI) on board theJames WebbSpace Telescope (JWST) uses three Si:As impurity band conduction (IBC) detector arrays. The output voltage level of each MIRI detector pixel is digitally recorded by sampling up the ramp. For uniform or low-contrast illumination, the pixel ramps become nonlinear in a predictable way, but in areas of high contrast, the nonlinearity curve becomes much more complex. The origin of the effect is poorly understood and currently not calibrated out of the data. Aims.We provide observational evidence of the brighter-fatter effect (BFE) in MIRI conventional and high-contrast coronagraphic imaging, low-resolution spectroscopy, and medium-resolution spectroscopy data, and we investigate the physical mechanism that gives rise to the effect on the MIRI detector pixel raw voltage integration ramps. Methods.We used public data from the JWST/MIRI commissioning and Cycle 1 phase. We also developed a numerical electrostatic model of the MIRI detectors using a modified version of the publicPoisson_CCDcode. Results.We find that the physical mechanism behind the BFE manifesting in MIRI data is fundamentally different to that of charge-coupled devices and photodiode arrays such as the Hawaii-XRG near-infrared detectors used by the NIRISS, NIRCam, and NIRSpec instruments on board JWST. Observationally, the BFE makes the JWST MIRI data yield 10–25% larger point sources and spectral line profiles as a function of the relative level of de-biasing of neighboring detector pixels. This broadening impacts the MIRI absolute flux calibration, time-series observations of faint companions, and point spread function modeling and subtraction. We also find that the intra-pixel 2D profile of the shrinking Si:As IBC detector depletion region directly impacts the accuracy of the pixel ramp nonlinearity calibration model.

Astronomy & Astrophysics↗
DOE OSTI2025

Ion Transport in Concentrated Crosslinked Solid Polymer Electrolytes

Crosslinking polymers is a common approach to create mechanically stable solid materials such as polymer electrolytes for lithium batteries. In conventional liquid electrolytes, the solvent molecules move freely to accommodate the field-induced motion of ions. However, in crosslinked polymer electrolytes, the rearrangement of polymer segments is constrained by the deformation limits of the network. Herein, we develop a new transport model that accounts for both the formation of concentration gradients and the elasticity of the electrolyte. The elasticity is incorporated by adding an additional term related to the entropy of crosslinked strands to the electrochemical potential of the salt. The resulting Crosslink Model contains two adjustable parameters: $\mathcal{N}$, the average number of monomers in a strand, and λ crit , the maximum strain the network can sustain. These solid-like constraints produce singularities in the governing transport equations, fundamentally altering the concentration profiles. Plateaus in salt concentrations emerge near the electrodes, and network elasticity introduces a strain overpotential. When compared to a Baseline Model ($\mathcal{N}$ → ∞, equivalent to concentrated solution theory), which predicts steepest gradients near the electrodes, both models yield similar current–voltage relationships. Model predictions are compared to electrochemical data for a poly(ethylene oxide)-based crosslinked polymer electrolyte.

Patel, Vivaan [University of California, Berkeley,↗
DOE OSTI2026

Superconducting phase diagram of multilayer square-planar nickelates

The discovery of superconductivity in square-planar nickelates has offered a rich materials platform to explore the origins of high-temperature superconductivity. However, experimental investigations have largely been limited to the infinite-layer R NiO 2 ( R , rare earth) nickelates. For this work we constructed a phase diagram of multilayer square-planar Nd n+1 Ni n O 2n+2 compounds and found signatures of superconductivity for dimensionality n = 4 to 8. Upon decreasing n , the superconducting anisotropy evolves owing to 4ƒ electron effects, and electronic structure characteristics approach cuprate-like behavior. Magnetic fluctuations persist from within the superconducting regime and into the overdoped, nonsuperconducting regime. The superconducting regime overlaps with that of chemically doped infinite-layer nickelates, demonstrating underlying commonalities as well as differences across varying structural realizations of square-planar nickelates. Our work establishes this layered template for creating new nickel-based superconductors.

36 MATERIALS SCIENCE↗
NASA NTRS1971

Battery simulation program

Computer program calculates spacecraft battery energy at specific times dictated by input sequence of recharge, discharge, and no activity phases. Recharge phases are assumed during times not specified, unless batteries are fully charged. Warnings are printed when energy falls below specified level. Program assumes two identical batteries.

Abreu, A.↗
DOE OSTI2026

Advanced blade-shaped thermal energy storage device: Development and application

Thermal energy storage (TES) using phase change materials (PCMs) is a promising approach for capturing and reusing excess thermal energy, yet widespread adoption is limited by low thermal conductivity, bulky configurations, and inadequate scalability. Here, this study presents a modular, blade-shaped TES prototype designed to address these challenges. The device integrates a lightweight aluminum shell, an embedded serpentine coil for active or passive heat exchange, and a cost-effective corrugated metal mesh for enhanced PCM thermal conductivity. With thickness-to-length and thickness-to-width ratios of 0.03 and 0.08, respectively, the blade-shaped TES achieves a compact, modular form factor suitable for space-constrained applications. Experimental testing demonstrated the efficient charge and discharge behavior of blade-shaped TES, capturing PCM superheating, phase-change transitions, and subcooling dynamics, with charging and discharging efficiencies of 94.9% and 94.6%, respectively. Also, the system can potentially achieve higher energy density than that of conventional TES designs. When integrated into a household refrigerator during the study, three blade-shaped TES modules successfully shifted 100% of peak-time compressor operation to off-peak hours, reducing energy consumption while maintaining more stable compartment temperatures. The blade-shaped TES's thin geometry, modularity, and enhanced thermal performance support scalable deployment across residential, commercial, and industrial applications, providing a versatile, cost-effective solution for high-efficiency, demand-flexible thermal energy management.

Blade-shaped↗
DOE OSTI2026

Modern insights into the mechanisms of neptunium oxalate decomposition

Neptunium oxalate (Np(C 2 O 4 ) 2 ·6H 2 O) is a historically relevant solid phase used in nuclear processing as a precursor for neptunium dioxide (NpO 2 ). Although Np oxalate has been synthesized and used for NpO 2 production for decades, the thermal decomposition mechanism of this phase remains poorly understood and has not been evaluated in over 30 years. Conflicting reports in historical literature suggest either a direct conversion from anhydrous oxalate to NpO 2 or a decomposition that includes the formation of Np carbonate or oxidized Np intermediate phases. In this work, we reexamine the decomposition pathway of Np(C 2 O 4 ) 2 ·6H 2 O using thermal analysis coupled with evolved gas analysis and temperature-dependent Raman spectroscopy to elucidate decomposition mechanisms and intermediate phases using modern analytical techniques. Thermal analysis revealed a three-stage decomposition process, including dehydration below 200 °C, oxalate breakdown between 170 and 370 °C, and NpO 2 formation by 500 °C. However, an unidentified plateau in the thermal data was observed during measurements. Raman spectroscopy confirmed the stages of decomposition, and in the analysis of potential intermediate phases, no carbonate phases or Np 2 O 5 were identified. Raman data suggest that residual oxalate or nonstoichiometric oxide are present during decomposition before pure NpO 2 is formed. These findings clarify aspects of the Np oxalate decomposition mechanism and address longstanding discrepancies in the literature, with a specific focus on Np-specific materials chemistry.

Lawson, Kathryn M. [Oak Ridge National Laboratory ↗