Search NASASearch

SEARCH · Search NASA

Results for “Diffuse interface methods”

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.

486 records · Page 11

First Wall Design of a Tokamak Pilot Plant Using a Monte Carlo Model for 3-D Heat Flux Deposition

We present a method for calculating the heat fluxes deposited on nonaxisymmetric tokamak first wall components, allowing for a first-of-its-kind model for power handling in the tokamak far scrape-off layer (SOL). The DIV3D Monte Carlo model features strict global power conservation and can calculate the finite cross-field plasma transport into magnetically-shadowed regions, which is significant when dealing with meter-scale shadows introduced by components such as poloidal limiters or antennas. As a case study, we apply the DIV3D model to inform the distribution of first wall poloidal limiters in an ARC-class reactor device. We demonstrate that discrete protection limiters can efficiently reduce peak heat fluxes on recessed breeder wall components in the presence of significant far-SOL plasma fluxes. By varying the toroidal periodicity and radial standoff depth of the limiters, we demonstrate one of the tradeoffs that must be considered in first wall design: more limiters provide greater protection, but at the cost of reduced breeding performance. We also present the impact that radial misalignments between limiters would have on first wall power loading.

Monte Carlo methods

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

Mu2e Production Target Thermal Test Project

The Mu2e experiment requires a production target that is capable of operating under extreme thermal conditions caused by an 8 GeV proton beam. This project’s objective supports the development of the Mu2e Pro- duction Target by testing the Stickman model’s thermal behavior. Angel Flores Luviano has assisted Jonathan Williams in progressing this project by contributing to the development of a Radiative Cooling Test Fixture (RCTF) that will be used to evaluate the thermal behavior of the new production target model. Engineering calculations were performed to an- alyze thermal performance and pressure drop within the cold well and wa- ter cooling circuit. These calculations also determined the optimal sizing for key components of the water system. An engineering note was made to document these calculations. CAD models of the water circuit piping, thermocouple mounting bars, and radiator chimney were developed, and prototype test rig components were fabricated using 3D printing. Future work will focus on continued RCTF development which includes control system integration, heater hardware design, and interfaces that can be scaled up to increase thermal capacity.

Flores Luviano, Angel Augusto [Unlisted, US]

Classical-Quantum Algorithm for Solving Stochastic Programs

Stochastic programming provides a rigorous mathematical framework for making decisions under uncertainty in a risk-aware manner. Two-stage stochastic programming is, perhaps, the simplest form of this framework. Here the first-stage variables represent decisions that must be made "here and now" in the face of uncertainty, while the second-stage variables are decisions made after uncertain events. However, the broad adoption of stochastic programming has been hindered by computational challenges caused by the two-stage stochastic programming formulation which requires solving an ensemble of optimization problems. Using quantum amplitude estimation (QAE), quantum computers have shown the theoretic ability to compute expectations with Monte-Carlo methods with quadratically fewer samples than classical methods. In this work, we present a quantum algorithm for computing the expectation term using QAE for given first-stage decisions. Further, we detail methods of computing gradient information from the quantum calculation enabling the application of classical gradient-based optimization techniques. The result is a classical-quantum hybrid method of solving two-stage stochastic programs. These techniques are demonstrated with computational experiments based an engineering optimization problem.

97 MATHEMATICS AND COMPUTING

Robust and High Throughput Tools for Characterization of Wet Cellulosic Biogas Feedstocks and Co-Products: Cooperative Research and Development Final Report, CRADA Number CRD-21-21472

This project will develop standardized analytical procedures to support understanding of the chemical composition of wet manure feedstocks, and the digestate formed following anaerobic digestion of organic waste to produce renewable natural gas (RNG). The developed methods will be validated to improve component mass closure, eliminate ambiguities around key constituents, and document the precision and accuracy of the methods. Following validation of the analytical laboratory methods for characterization of wet organic waste feedstocks and digestate, NLR will use the same samples to build and validate near-infrared (NIR) models for high throughput analysis, potentially deployable at commercial facilities, enabling increased production of low carbon intensity RNG from wet organic wastes.

09 BIOMASS FUELS

Complex orders and chirality in the classical Kitaev-Γ model

It is well recognized that the low-energy physics of many Kitaev materials is governed by two dominant energy scales, the Ising-type Kitaev coupling 𝐾 and the symmetric off-diagonal Γ coupling. An understanding of the interplay between these two scales is therefore the natural starting point toward a quantitative description that includes subdominant perturbations that are inevitably present in real materials. This study focuses on the classical 𝐾−Γ model on the honeycomb lattice, with a specific emphasis on the region 𝐾< 0 and Γ > 0 , which is the most relevant for the available materials and which remains enigmatic in both quantum and classical limits, despite much effort. We employ large-scale Monte Carlo simulations on specially designed finite-size clusters and unravel the presence of a complex multisublattice magnetic order in a wide region of the phase diagram, whose structure is characterized in detail. We show that this order can be quantified in terms of a coarse-grained scalar-chirality order, featuring a counterrotating modulation on the two spin sublattices. Here, we also provide a comparison to previous studies and discuss the impact of quantum fluctuations on the phase diagram.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Evaluation of high-temperature properties of binary, eutectic Ni-X alloys (X = B, Ca, Ce, La, Y, Zr)

We investigated the high-temperature behavior of binary eutectic nickel-base alloys with low-solubility elements, towards developing alloys with new combinations of structural and functional (e.g. thermal or corrosion) properties. Based on thermodynamic screening, we cast binary near-eutectic Ni-B, Ni-Ca, Ni-Ce, Ni-La, Ni-Y, and Ni-Zr alloys and characterized their microstructure and high-temperature properties. All alloys (except Ni-Ca) contain a fine eutectic microstructure, containing submicron alternating lamellae of γ-Ni and continuous intermetallic phases. The eutectic microstructure did not significantly coarsen at 700°C (except for Ni-B), but measurably coarsened in all alloys at 900°C, most quickly in Ni-B, followed by Ni-Ce, Ni-La, Ni-Y, and Ni-Zr. Coarsening in these alloys occurs by fault migration; diffusion within the intermetallic, and its interfacial energy, likely influence coarsening kinetics. At room-temperature, all alloys possess negligible ductility due to the continuous, brittle intermetallic. At 700 and 900°C, the Ni-Zr alloy showed the best combination of strength and ductility. Oxidation experiments at 900°C on Ni-Zr and Ni-Ce revealed heavy internal oxidation of the continuous intermetallic phase. Overall, the Ni-Zr alloy showed good high-temperature coarsening resistance and strength, though ductility and oxidation resistance require improvement for use in structural applications. Possible routes for improving alloy properties via additive manufacturing are discussed.

Ekaputra, Clement [ORNL] (ORCID:0000000307166479)

Ferroelectric Fractals: Switching Mechanism of Wurtzite AlN

The advent of wurtzite ferroelectrics is enabling new ferroelectric devices for computer memory that have the potential to bypass the von Neumann bottleneck due to their robust polarization and silicon compatibility. However, the atomistic switching mechanism of wurtzites is still undetermined due to the limitations of density functional theory simulation size and experimental temporal and spatial resolution. Thus, physics-informed materials engineering to reduce coercive field and breakdown in these devices has been limited. In this work, the atomistic mechanism of domain wall migration and domain growth in aluminum nitride-based wurtzites is uncovered using molecular dynamics and Monte Carlo simulations. We reveal the anomalous switching mechanism of fast 1D single columns of atoms propagating from a slow-moving 2D fractallike domain wall. We find that the critical nucleus is a single aluminum ion that breaks its bond with one nitrogen and bonds to another nitrogen; this creates a cascade that flips atoms directly only in the same column, due to the extreme locality (sharpness) of the domain walls in wurtzites. We further show how the fractallike shape of the domain wall in the 2D plane breaks assumptions in the Kolmogorov, Avrami, and Ishibashi (KAI) model and leads to the anomalously fast switching in wurtzite structured ferroelectrics.

36 MATERIALS SCIENCE

Capturing Secondary Kinetic Instabilities in Three‐Dimensional Dayside Reconnection Using an Improved Gradient‐Based Closure

Magnetic reconnection is a highly dynamic process that excites a wide variety of kinetic waves and instabilities. Transverse current sheet instabilities such as the lower-hybrid drift and secondary drift-kink instabilities in particular have been shown by kinetic simulations to modify the reconnection and introduce significant turbulence and mixing to the reconnection layer. Past studies using the ten-moment fluid model to capture important kinetic physics such as the electron inertia and full representation of the pressure tensor proved advantageous to a two-fluid representation of reconnection, but the model struggled when using a local relaxation closure for the heat flux to replicate the current sheet instabilities and subsequent mixing seen in kinetic simulations. This work uses the Gkeyll software framework to perform simulations of asymmetric reconnection based on the 16 October 2015 MMS crossing of a diffusion region, the Burch event. An improved gradient-based heat flux closure is implemented, showing significant improvement in secondary kinetic instabilities that grow in the current sheet. These instabilities generate turbulence which leads to growth of secondary magnetic islands and flux ropes.

Bradshaw, K. [Princeton University, NJ (United Sta

Role of Fluid and Temperature in Fracture Mechanics and Coupled THMC Processes for Enhanced Geothermal Systems

The ability to sustain high energy extraction efficiency from Enhanced Geothermal Systems (EGS) is affected by the ability to measure, characterize, and predict the effect of different perturbations that arise from fluid injection rates, fluid temperature, and shut-in conditions on existing and stimulated fracture systems throughout a subsurface reservoir’s lifecycle. The difficulty arises from limited knowledge of the interaction of fluid and temperature driven fractures with frictional interfaces that govern local deformation and frictional behavior in the surrounding rock. Thermo-poro-mechanical coupled processes tend to dominate this interaction and strongly affect the permeability and local stress distributions that impact efficiency and prevent optimal stimulation conditions over multi-decadal time frames.

Pyrak-Nolte, Laura [Purdue Univ., West Lafayette,

Microstructure, mechanical, and thermal properties of compositionally complex (Hf,Zr,Nb,Ti)B 2 ‒LaB 6 ceramics

Novel compositionally complex borides, (Hf,Zr,Nb,Ti)B 2 and (Hf,Zr,Nb,Ti)B 2 ‒LaB 6 , were fabricated using spark plasma sintering process. (Hf,Zr,Nb,Ti)B 2 ‒LaB 6 exhibits a dual-phase microstructure, in which (Hf,Zr,Nb,Ti)B 2 is a primary phase with the hexagonal structure and LaB 6 is a secondary phase with a cubic structure. The mechanical properties of both (Hf,Zr,Nb,Ti)B 2 and (Hf,Zr,Nb,Ti)B 2 ‒LaB 6 are comparable, with a combination of high hardness and moderate fracture toughness. Thermal diffusivity and conductivity of (Hf,Zr,Nb,Ti)B 2 are much lower than the individual transition metal borides but are significantly increased by the addition of LaB 6 . Herein, it is implied that the thermal properties of boride ceramics can be controlled through the appropriate design of principal metal element compositions.

borides

Influence of Powder Characteristics and Processing Methods on Creep Performance of Powder Metallurgy Hot Isostatic Pressed SS316

The U.S. nuclear energy expansion goals are driving the demand for manufacturing routes that can rapidly produce large, complex, near net shape components. Powder metallurgy hot isostatic pressing (PM HIP) is an advanced manufacturing technique that can be economically scaled-up, while alleviating the supply chain challenges that forging and casting face in terms of cost and lead time constraints. This makes PM-HIP a viable technology to aid and accelerate large-scale part manufacturing for nuclear applications. However, large-scale qualification and deployment of this technology require a thorough understanding of the influence of powder feedstock quality, powder handling history, hot isostatic pressing (HIP) parameters, and subsequent heat treatment on microstructural evolution and elevated temperature mechanical performance. The present work focusses on 316 austenitic stainless steel (SS316) which is most commonly used in high temperature environments for nuclear applications Results from this study show that PM HIPed SS316 meets ASME tensile requirements at room temperature and at elevated temperature. However, creep performance of PM-HIPed SS316 remains inferior to its wrought counterpart, demonstrating that tensile performance alone is not a reliable metric for long duration high temperature integrity. Further, this report delineates powder derived microstructural features that govern creep damage, with key evidences pointing to “microstructural inheritance” from gas atomized powder feedstocks. Commercial SS316 powders of varying chemical compositions and recycling histories were studies, and the results showed large differences in elemental segregation, oxide surface layers and secondary phase distributions. Multi-scale characterization revealed segregation of chromium, molybdenum, manganese and silicon at the boundaries and the precipitation of manganese-, silicon-, and molybdenum-oxides. During HIP consolidation, these surface oxides transform into decorated prior particle boundaries (PPBs) and grain boundary inclusions that persist through conventional post-HIP solution annealing treatment. The retained oxides in post-HIP microstructures were found to influence grain growth, precipitation behavior, and ultimately creep cavitation and fracture. Such post-HIP heat treatments are therefore limited by a complex trade-off between grain growth, and oxide coarsening which aggravate creep damage by acting as nucleation sites for cavities. The objective of this work is to establish an integrated processing–structure–property framework for PM-HIP 316 stainless steel by investigating the influence of powder feedstock characteristics in pre- and post-HIP processing as well as to understand the significance of post-HIP heat treatment on microstructural evolution and creep properties. The results from this report emphasize the significance of powder feedstock integrity in improving creep performance of PM-HIPed SS316, by highlight the effect of rapid solidification, elemental segregation, oxide formation and powder recycling on microstructural defect inheritance following HIP consolidation. Rather than considering HIP processing, solution annealing, and mechanical performance independently, this report treats powder production, HIP consolidation, post-HIP thermal processing, and creep deformation as interconnected stages within a continuous metallurgical process. The resulting framework will provide a scientific basis for developing feedstock engineering strategies capable of improving long-term reliability of PM-HIP stainless steels and accelerating their qualification for advanced nuclear applications.

Ajjarapu, Pavan [Oak Ridge National Laboratory (OR

Mesh-based multiphysics coupling acceleration for fusion neutronics through clustering for fusion blanket applications

Accurate modeling of particle transport within fusion blankets is essential for predicting performance metrics such as heat deposition and the tritium breeding ratio (TBR). However, high-fidelity coupling of thermal fluids from computational fluid dynamics (CFD) to neutronics simulations often incurs significant computational costs due to the complexity of surface intersection calculations in Monte Carlo codes. This paper presents an accelerated multiphysics coupling method for neutronics that utilizes hierarchical agglomerative clustering to map complex material property distributions to a neutronics model. Implemented within the fusion reactor design and assessment (FREDA) framework, the method leverages existing Python packages to automate the creation of clustered geometries for OpenMC. The approach is demonstrated on a sector model of an ARC-class tokamak with an immersion molten salt blanket, and an simple geometry with varying isotopic concentrations. Results show that the clustering method significantly reduces computational burden without compromising fidelity, providing a foundation for agile iteration of neutronics simulations involving multiple coupled material properties.

Bae, Jin Whan [ORNL] (ORCID:0000000326548907)

Glancing Angle Deposition in Gas Sensing: Bridging Morphological Innovations and Sensor Performances

Glancing Angle Deposition (GLAD) has emerged as a versatile and powerful nanofabrication technique for developing next-generation gas sensors by enabling precise control over nanostructure geometry, porosity, and material composition. Through dynamic substrate tilting and rotation, GLAD facilitates the fabrication of highly porous, anisotropic nanostructures, such as aligned, tilted, zigzag, helical, and multilayered nanorods, with tunable surface area and diffusion pathways optimized for gas detection. This review provides a comprehensive synthesis of recent advances in GLAD-based gas sensor design, focusing on how structural engineering and material integration converge to enhance sensor performance. Key materials strategies include the construction of heterojunctions and core–shell architectures, controlled doping, and nanoparticle decoration using noble metals or metal oxides to amplify charge transfer, catalytic activity, and redox responsiveness. GLAD-fabricated nanostructures have been effectively deployed across multiple gas sensing modalities, including resistive, capacitive, piezoelectric, and optical platforms, where their high aspect ratios, tailored porosity, and defect-rich surfaces facilitate enhanced gas adsorption kinetics and efficient signal transduction. These devices exhibit high sensitivity and selectivity toward a range of analytes, including NO2, CO, H2S, and volatile organic compounds (VOCs), with detection limits often reaching the parts-per-billion level. Emerging innovations, such as photo-assisted sensing and integration with artificial intelligence for data analysis and pattern recognition, further extend the capabilities of GLAD-based systems for multifunctional, real-time, and adaptive sensing. Finally, current challenges and future research directions are discussed, emphasizing the promise of GLAD as a scalable platform for next-generation gas sensing technologies.

Chemistry

Mitigation of polysulfide shuttle effect in Li-S batteries through catalytic disproportionation reaction

Polysulfides are poorly retained within porous cathodes and readily diffuse into the electrolyte over time, leading to the well-known shuttle effect that undermines the reversibility of Li-S batteries. Here, in this study, we demonstrate that catalytic disproportionation of polysulfides provides an effective pathway to suppress this process by rapidly converting dissolved species into solid sulfur and sulfides, thereby preventing their migration into the electrolyte. Fundamentally, the sluggish kinetics of sulfur redox reactions are responsible for the accumulation and redistribution of soluble polysulfides in the bulk electrolyte. By accelerating these kinetics, catalyzed disproportionation not only confines sulfur within the conductive cathode matrix but also promotes the homogeneous precipitation of Li₂S₂/Li₂S, which enhances electrochemical reversibility and cycling stability. Using nitrogen-doped carbon (NC800) as a model catalyst, we reveal its ability to drive a pseudo-16-electron reduction pathway, leading to a single dominant Li₂S product and uniform deposition within the porous framework. In contrast, a non-catalytic carbon (KB) yields multiple polysulfide intermediates and heterogeneous deposition. The mechanistic insights provided here highlight the pivotal role of catalytic disproportionation in reshaping sulfur redox pathways and offer a rational strategy for mitigating polysulfide shuttling in practical Li-S pouch cells.

25 ENERGY STORAGE

Impact of the Annealing Temperature on the Local Se Distribution and Optical Properties in CdSexTe1-x Solar Cells by Multimodal X-ray Microscopy

The introduction of Se in CdTe solar-cell absorbers has led to multiple beneficial effects in the devices, including an increased short-circuit current and an extended charge-carrier lifetime. A critical manufacturing step of CdSexTe1-x solar cells is CdCl2 annealing at high temperatures, during which Se diffuses through the absorber layer. Understanding the local Se distribution in the absorber and its impact on optical properties of the material is key for the further development of this thin-film solar cell technology. In a multimodal x-ray scanning microscopy approach, we investigated the impact of the annealing temperature of CdSexTe1-x solar cells on the Se distribution and on optical characteristics of the material. For this purpose, we exploited an upgraded x-ray excited optical luminescence (XEOL) detection unit that enables simultaneous spectrally and temporally resolved XEOL mapping at synchrotron facilities that are ideal for studying thin-film solar-cell materials with hard x-rays featuring high penetration depth and lateral resolution.

14 SOLAR ENERGY

Integration of Online Cross-Section Generation Capability with Depletion and Transient Solvers in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE)-based reactor multiphysics analysis application jointly developed by Argonne and Idaho National Laboratories under the DOENE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In FY25, an online crosssection generation capability based on the Self-Shielding Application Programming Interface (SSAPI) was demonstrated for TRISO-fueled reactor problems under steady-state conditions. This fiscal year, that capability was extended to support depletion and transient multiphysics calculations, enabling high-fidelity analyses that generate self-shielded cross sections on the fly from the actual evolving composition and temperature states rather than from pre-tabulated libraries. For depletion, a two-way coupling was established in which SSAPI computes compact-averaged self-shielded cross sections that the depletion solver then uses to advance the Bateman equations, with the updated compositions returned to SSAPI at each step; the depletion module was refactored to support both library-based and SSAPI-based cross sections, and additional logic was added to track daughter isotopes and to exclude minor isotopes for efficiency. For transient analysis, the SSAPI multigroup library was extended with the kinetics data required for time-dependent calculations, the Improved Quasi-Static (IQS) scheme was coupled with SSAPI, and several supporting capabilities were implemented, including a self-shielding treatment that lets control rods and drums move within a self-shielded model, which had previously been impossible and had ruled out rod- and drum-movement transients with on-the-fly cross sections altogether, a new mixing scheme for delayed-neutron precursor decay constants, a checkpoint-based restart workflow, and performance improvements such as pointwise cross-section interpolation and the bypassing of unnecessary Dancoff factor calculations. The implemented capabilities were verified against Serpent Monte Carlo solutions. For depletion, a prismatic pin-cell problem based on a Next Generation Nuclear Plant (NGNP) Very High Temperature Reactor benchmark showed excellent agreement, with eigenvalue differences within 200 pcm over the entire burnup range (up to 140 MWD/kgU) and fission-product and actinide inventories agreeing to within 0.8% and 2.5%, respectively; a heat-pipe microreactor assembly problem with a much higher fuel loading confirmed the same behavior and quantified the bias introduced when the multigroup equivalence effect is neglected. For transient analysis, a pin-cell problem with a step reactivity insertion and temperature feedback reproduced the analytically expected asymptotic power and showed close agreement between the direct and IQS solutions, and a two-dimensional microreactor core problem with control-drum rotation exercised the new moving-drum self-shielding treatment and demonstrated successful coupling of the online crosssection generation with both the direct and IQS transient methods. The capability was further exercised on a full-core pebble-bed problem, in which Griffin was coupled with the System Analysis Module (SAM) to simulate load-following operation of the gPBR with the Doppler feedback resolved at the TRISO fuel kernel temperature. These developments in Griffin provide a convenient, high-fidelity approach to cross-section generation for advanced thermal reactors with geometrically complex and highly heterogeneous configurations, including TRISO-fueled prismatic and pebble-bed systems, and support steady-state, depletion, and transient multiphysics calculations. They also enable self-shielded cross sections to be evaluated directly at the actual coupled state of the system, thereby establishing a foundation for high-fidelity, fully coupled multiphysics analysis of advanced reactors

Park, H.

Robust Spectral Anomaly Detection in EELS Spectral Images via 3D Convolutional Variational Autoencoders

Abstract A 3D Convolutional Variational Autoencoder (3D‐CVAE) is introduced for automated anomaly detection in electron energy‐loss spectroscopy spectrum imaging (EELS‐SI) data. This approach leverages the full 3D structure of EELS‐SI data to detect subtle spectral anomalies while preserving both spatial and spectral correlations across the datacube. By employing cross‐entropy loss and training on bulk spectra, the model learns to reconstruct bulk features characteristic of the defect‐free material. In exploring methods for anomaly detection, both the 3D‐CVAE approach and principal component analysis (PCA) are evaluated, testing their performance using FeL‐edge ΔEpeak shifts designed to simulate material defects. These results show that 3D‐CVAE achieves superior anomaly detection and maintains consistent performance across various shift magnitudes. The method demonstrates clear bimodal separation between bulk and anomalous spectra, enabling reliable classification. Further analysis verifies that lower‐dimensional representations are robust to anomalies in the data. While performance advantages over PCA diminish with decreasing anomaly concentration, our method maintains high reconstruction quality even in challenging, noise‐dominated spectral regions. This approach provides a robust framework for unsupervised automated detection of spectral anomalies in EELS‐SI data, particularly valuable for analyzing complex material systems.

Chemistry