Search NASASearch

SEARCH · Search NASA

Results for “finite element simulation”

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.

406 records · Page 10

Final Technical Report for DE-SC0022206

This project developed foundational genetic, genomic, and epigenetic tools for anaerobic fungi (Neocallimastigomycota), a group of microorganisms with exceptional natural abilities to deconstruct lignocellulosic biomass. Efficient biomass deconstruction remains a major barrier to economical production of renewable fuels, chemicals, and materials from agricultural and forestry residues. The project sought to enable mechanistic studies and future engineering of anaerobic fungi by improving genomic resources, establishing methods for gene expression, and investigating epigenetic regulation of biomass-degrading pathways. Major accomplishments included generation of the first chromosome-scale genome assemblies for multiple anaerobic fungal species, providing publicly available genomic resources that support both engineering and fundamental biological research. The project established the first reproducible system for heterologous gene expression in anaerobic fungi and identified genomic features and mobile genetic elements that may support future development of stable transformation technologies. In parallel, the project demonstrated direct conversion of untreated lignocellulosic biomass into fuels and specialty chemicals through a fungal-yeast bioprocess and identified anaerobic fungal enzymes with utility for metabolic engineering. The research also revealed that epigenetic regulation plays an important role in controlling fungal gene expression and enzyme production, identifying potential strategies for enhancing biomass degradation. Collectively, this work established anaerobic fungi as a tractable emerging platform for bioenergy and biomanufacturing research, generated valuable public resources, trained the next generation of researchers, and advanced DOE-BER goals related to predictive biology, sustainable bioprocessing, and the circular bioeconomy.

Solomon, Kevin [University of Delaware] (ORCID:000

Phase-field modeling of stored-energy-driven grain growth with intra-granular variation in dislocation density

Abstract We present a phase-field (PF) model to simulate the microstructure evolution occurring in polycrystalline materials with a variation in the intra-granular dislocation density. The model accounts for two mechanisms that lead to the grain boundary migration: the driving force due to capillarity and that due to the stored energy arising from a spatially varying dislocation density. In addition to the order parameters that distinguish regions occupied by different grains, we introduce dislocation density fields that describe spatial variation of the dislocation density. We assume that the dislocation density decays as a function of the distance the grain boundary has migrated. To demonstrate and parameterize the model, we simulate microstructure evolution in two dimensions, for which the initial microstructure is based on real-time experimental data. Additionally, we applied the model to study the effect of a cyclic heat treatment (CHT) on the microstructure evolution. Specifically, we simulated stored-energy-driven grain growth during three thermal cycles, as well as grain growth without stored energy that serves as a baseline for comparison. We showed that the microstructure evolution proceeded much faster when the stored energy was considered. A non-self-similar evolution was observed in this case, while a nearly self-similar evolution was found when the microstructure evolution is driven solely by capillarity. These results suggest a possible mechanism for the initiation of abnormal grain growth during CHT. Finally, we demonstrate an integrated experimental-computational workflow that utilizes the experimental measurements to inform the PF model and its parameterization, which provides a foundation for the development of future simulation tools capable of quantitative prediction of microstructure evolution during non-isothermal heat treatment.

Materials Science

SOLPS-ITER modeling of a dedicated divertor for negative triangularity operation on DIII-D

The design of a new dedicated divertor for negative triangularity (NT) operation on DIII-D with neutral baffles and pumping is informed by SOLPS-ITER transport modeling. This dedicated NT divertor is the latest step in a progression of NT shapes with various divertor characteristics explored on DIII-D, including NT shapes at reduced triangularity and a campaign with stronger shaping that included new armored components on the outboard side. SOLPS simulations played a key role in these divertor designs. Interpretive simulations, using cross-field diffusivities constrained by experimental data in the NT Shelf shape were used to inform the design of the 2023 armor campaign components. A similar procedure used armor campaign data to predict conditions for the dedicated NT divertor. The predictive simulations were used to assess the divertor fluxes, detachment threshold, pumped flux, and neutral leakage. For the dedicated NT divertor, SOLPS simulations and two-point-modeling were used to show the relative impact of magnetic topology (mainly longer connection length) and divertor closure on the divertor conditions relative to the armor campaign. It is predicted that the dedicated NT divertor reaches detachment (measured by target ion flux rollover) at a lower upstream density (≈(1.75−−2.4)×1⁢019m−3) as compared to the armor campaign shape. For the preliminary design geometry, divertor closure reduces the neutral leakage by ≈10%. Parametric optimization indicating further ≈20%–60% improvement in the leakage flux and recycled flux crossing the pump entrance is possible for relatively minor changes to the divertor and baffle layout.

Lore, Jeremy [ORNL] (ORCID:000000029192465X)

Estimating Return on Investment for Energy Technical Assistance Programs

The U.S. Department of Energy's Office of State and Community Energy Programs engaged the National Laboratory of the Rockies to assess the return on investment (ROI) of technical assistance (TA) programs that support state, local, and Tribal energy planning. Although TA delivers value through capacity building, stakeholder engagement, and knowledge transfer, these benefits are often intangible and challenging to monetize. This study reviews existing ROI frameworks and synthesizes the most relevant elements into a hybrid approach tailored to energy TA programs. The proposed framework integrates monetary and non-monetary outcomes through early logic model development, baseline data collection, and the use of proxies for intangible benefits. As a case study, this paper applies this approach to the Communities Local Energy Action Program (Communities LEAP), demonstrating how ROI can inform program design, data strategy, and performance assessment. Findings underscore that ROI should be applied selectively and planned from the outset to ensure data alignment and attribution accuracy. The framework offers TA practitioners a structured approach that can be leveraged for future programs to evaluate and communicate the multifaceted value of TA investments.

29 ENERGY PLANNING, POLICY, AND ECONOMY

Large language model-driven database for thermoelectric materials

Thermoelectric materials have the ability to convert waste heat into electricity, offering a valuable solution for energy harvesting. However, their widespread use is hindered by low conversion efficiency, the reliance on expensive rare earth elements, and the environmental and regulatory concerns associated with lead-based materials. A fast and cost-effective way to identify highly efficient thermoelectric materials is through data-driven methods. These approaches rely on robust and comprehensive datasets to train models. Although there are several databases on thermoelectric materials, there is still a need to collect and integrate experimental data from peer-reviewed research articles to capture diverse compositions and properties of materials. Here, in this work, we developed a comprehensive database of 7,123 thermoelectric compounds, containing key information such as chemical composition, structural detail, seebeck coefficient, electrical and thermal conductivity, power factor, and figure of merit (ZT). We used the GPTArticleExtractor workflow, powered by large language models (LLM), to extract and curate data automatically from the scientific literature published in Elsevier journals. This process enabled the creation of a structured database that addresses the challenges of manual data collection. The open access database could stimulate data-driven research and advance thermoelectric material analysis and discovery.

Database

Enter the AHU (36th Chamber of ASHRAE): A Multi-site Field Study of ASHRAE G36

Despite being recognized as the best practice for advanced building controls, ASHRAE Guideline 36 (G36) has seen slow adoption in retrofit cases. Decisionmakers lack credible field evidence to justify the time and person-power investment. Most prior analyses have relied on software simulations, which overlook implementation challenges and fail to persuade owners to move from models to real-world deployment. This paper presents a multi-site field study of G36 performance, drawing on measured results from 17 projects across diverse building types and climate zones. The analysis disaggregates outcomes by the most widely adopted air handling unit (AHU) based G36 strategies, including trim-and-respond approaches to supply air temperature (SAT) and duct static pressure (DSP) reset, and economizer controls. Results for controls re programming implementations are encouraging, with HVAC savings ranging from 2% - 49%, with a median of 18%. The range aligns with simulation study findings showing 1% - 46% savings through these three strategies. These findings show that even without capital investment in control infrastructure upgrades, existing building owners are reaping significant benefits from updating HVAC sequences of operation to industry best-practice solutions. The results give practitioners and decisionmakers a reference point on what to expect, which strategies deliver, and how field performance may compare to simulation.

Deshpande, Reva

Neutron detector response modeling in NOvA

Neutrons can present a significant challenge for neutrino experiments in which energy reconstruction is critical. With the ability to escape detection completely and with a weak correlation between their kinetic energy and any eventual energy deposition, it is difficult to fully account for neutrons produced in neutrino interactions. This in turn leads to significant model dependence when evaluating neutron-related systematic uncertainties. The NOvA experiment is a long-baseline neutrino oscillation experiment with a high-statistics sample of antineutrino data collected by its near detector. We report an excess relative to data of simulated neutron candidates with low energy depositions when using standard Geant4 physics lists. The simulation excess is traced to an overabundance of secondary photons produced from interactions of neutrons with kinetic energy greater than \SI{20}{\mega\eV}. Improved agreement with data is obtained by applying the data-driven neutron-on-carbon \menate model for neutrons between \SI{20}{\mega\eV} and ${\sim}$\SI{100}{\mega\eV}. With \menate, the residual oversimulation is more uniform across the calorimetric neutron energy spectrum, suggesting possible overproduction of primary neutrons by the GENIE neutrino interaction generator. These results motivate the adoption of \menate-supplemented Geant4 simulation as the nominal simulation in the production of future \nova simulation.

Abubakar, S.

Realization of two-dimensional discrete time crystals with anisotropic Heisenberg coupling

A discrete time crystal (DTC) is an out-of-equilibrium phase of matter that spontaneously breaks discrete time-translation symmetry. Previous studies have been limited to a set of models with Ising-like couplings - and mostly only in one dimension - thus precluding our understanding of the existence (or not) of DTCs in models with more realistic interactions. In this work, by combining the latest generation of IBM quantum processors with state-of-the-art tensor network methods, we demonstrate the existence of a DTC in a two-dimensional system governed by anisotropic Heisenberg interactions. We uncover a rich phase diagram encompassing spin-glass, ergodic, and time-crystalline phases, and identify the interplay of initialization, interaction anisotropy, and driving protocols in stabilizing the DTC phase. By extending the study of Floquet matter beyond simplified models, we lay the groundwork for exploring how driven systems bridge the gap between quantum coherence and emergent non-equilibrium thermodynamics.

Phase transitions and critical phenomena

Delineating the Impact of Transition‐Metal Crossover on Solid‐Electrolyte Interphase Formation with Ion Mass Spectrometry

Abstract Lithium‐metal batteries (LMB) employing cobalt‐free layered‐oxide cathodes are a sustainable path forward to achieving high energy densities, but these cathodes exhibit substantial transition‐metal dissolution during high‐voltage cycling. While transition‐metal crossover is recognized to disrupt solid‐electrolyte interphase (SEI) formation on graphite anodes, experimental evidence is necessary to demonstrate this for lithium‐metal anodes. In this work, advanced high‐resolution 3D chemical analysis is conducted with time‐of‐flight secondary‐ion mass spectrometry (TOF‐SIMS) to establish spatial correlations between the transition metals and electrolyte decomposition products found on cycled lithium‐metal anodes. Insights into the localization of various chemistries linked to crucial processes that define LMB performance, such as lithium deposition, SEI growth, and transition‐metal deposition are deduced from a precise elemental and spatial analysis of the SEI. Heterogenous transition‐metal deposition is found to perpetuate both heterogeneous SEI growth and lithium deposition on lithium‐metal anodes. These correlations are confirmed across various lithium‐metal anodes that are cycled with different cobalt‐free cathodes and electrolytes. An advanced electrolyte that is stable to higher voltages is shown to minimize transition‐metal crossover and its effects on lithium‐metal anodes. Overall, these results highlight the importance of maintaining uniform SEI coverage on lithium‐metal anodes, which is disrupted by transition‐metal crossover during operation at high voltages.

Chemistry

Machine Learning for Multipactor Susceptibility Prediction in Planar RF Gaps

Multipactor discharge is a nonlinear electron avalanche that limits the performance of high-power radio-frequency (RF) and vacuum electronic devices. Predicting multipactor susceptibility traditionally relies on Monte Carlo or particle-in-cell (PIC) simulations, which become computationally expensive for large parametric studies. In this work, we present a supervised machine-learning (ML) framework for prediction of multipactor susceptibility in a two-surface planar geometry. The models are trained using high-fidelity PIC simulation generated susceptibility data and learn the relationship between operational parameters, geometry, and material-dependent secondary electron emission properties. The proposed approach enables rapid reconstruction of susceptibility charts while preserving the physical structure of multipactor growth regions.

43 PARTICLE ACCELERATORS

Resonant Metasurfaces with Van Der Waals Hyperbolic Nanoantennas and Extreme Light Confinement

This work reports on a metasurface based on optical nanoantennas made of van der Waals material hexagonal boron nitride. The optical nanoantenna made of hyperbolic material was shown to support strong localized resonant modes stemming from the propagating high-k waves in the hyperbolic material. An analytical approach was used to determine the mode profile and type of cuboid nanoantenna resonances. An electric quadrupolar mode was demonstrated to be associated with a resonant magnetic response of the nanoantenna, which resembles the induction of resonant magnetic modes in high-refractive-index nanoantennas. The analytical model accurately predicts the modes of cuboid nanoantennas due to the strong boundary reflections of the high-k waves, a capability that does not extend to plasmonic or high-refractive-index nanoantennas, where the imperfect reflection and leakage of the mode from the cavity complicate the analysis. In the reported metasurface, excitations of the multipolar resonant modes are accompanied by directional scattering and a decrease in the metasurface reflectance to zero, which is manifested as the resonant Kerker effect. Van der Waals nanoantennas are envisioned to support localized resonances and can become an important functional element of metasurfaces and transdimensional photonic components. By designing efficient subwavelength scatterers with high-quality-factor resonances, this work demonstrates that this type of nanoantenna made of naturally occurring hyperbolic material is a viable substitute for plasmonic and all-dielectric nanoantennas in developing ultra-compact photonic components.

Chemistry

Cation-Diffusive Carbon Interlayers Stabilize Na Metal and Double the Current in Na-S Redox-Flow Batteries for Grid-Scale Energy Storage

The sodium-sulfur nonaqueous redox-flow batteries (Na-S NARFBs) using earth-abundant elements are highly attractive due to the low material cost and improved energy density for grid-scale energy storage. However, the low current performance, poor Na0/Na+ redox kinetics, and Na dendrite growth pose severe challenges. We introduce cation-diffusive layers (CDLs): thin and Na+ affinitive interlayers at the Na anode that direct Na+ transport and stabilize Na deposition. Benchmarking three archetypal materials—carbon paper (CP), glass microfiber paper (GF), and foam—across Na-Na and Na-Cu, and Na-S cells identifies CP as the optimum. CP reduces symmetric cell overpotential by more than 70%, achieves 98% Na plating-stripping efficiency, and doubles the Na-S cell current density from 0.5 to 1.0 mA cm−2 without sacrificing capacity or efficiency. Ex situ electrochemical and SEM/XPS analysis, combined with molecular dynamics (MD) studies, reveal that electron-rich carbon fibers disperse supporting salt aggregates, enrich near-surface Na+ density, and create ion transport pathways for fast Na0/Na+ exchange while mitigating membrane degradation. Because of the ion-centric mechanism, CDLs can be generalized to other metal-anode designs. Further, this work establishes CDL design rules—cationic affinity and appropriate micro/nanostructure—as a simple, scalable route to high-current, durable metal-anode flow batteries.

Wu, Wenda [ORNL] (ORCID:000900033307687X)

Spinel high-entropy oxides (FeNiCrMnZnX) 3 O 4 (X = Al, mg) as anode materials for high-performance lithium-ion batteries

To address the high cost, cobalt dependency, and resource constraints typical of conventional high-entropy oxide (HEO) anodes, this study reports the successful synthesis of two Co-free, six-component spinel-type HEOs(FeNiCrMnZnAl) 3 O 4 (HEO-Al) and (FeNiCrMnZnMg) 3 O 4 (HEO-Mg), via a sol-gel method. The distinct effects of Al 3+ and Mg 2+ incorporation on the electrochemical performance and lithium storage kinetics were systematically investigated. XRD, Raman, and TEM characterizations confirm that both materials possess a pure spinel phase, uniform particle size, and homogeneous elemental distribution. Notably, electrochemical evaluations reveal that HEO-Al delivers a superior reversible capacity of 480.7 mAh g −1 after 100 cycles at 0.1 A g −1 , and maintains 354.5 mAh g −1 after 1000 long-term cycles at 1 A g −1 , significantly outperforming HEO-Mg. Kinetic analysis indicates that HEO-Al exhibits lower charge transfer resistance, a higher Li + diffusion coefficient, and a pseudocapacitive contribution of up to 82%. Furthermore, these findings demonstrate that Al substitution effectively optimizes the structural stability and lithium storage kinetics of Co-free HEOs, providing a viable strategy for designing low-cost, highly stable HEO anode systems.

Anode materials

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)

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

String-Breaking Dynamics in Quantum Adiabatic and Diabatic Processes

Confinement prohibits isolation of color charges, e.g., quarks, in nature via a process called string breaking : the separation of two charges results in an increase in the energy of a color flux, visualized as a string, connecting those charges. Eventually, creating additional charges is energetically favored, hence breaking the string. Such a phenomenon can be probed in simpler models, including quantum spin chains, enabling enhanced understanding of string-breaking dynamics. A challenging task is to understand how string breaking occurs as time elapses, in an out-of-equilibrium setting. This work establishes the phenomenology of dynamical string breaking induced by a gradual increase of string tension over time. It, thus, goes beyond instantaneous quench processes and enables tracking the real-time evolution of strings in a more controlled setting. We focus on domain-wall confinement in a family of quantum Ising chains. Our results indicate that, for sufficiently short strings and slow evolution, string breaking can be described by the transition dynamics of a two-state quantum system akin to a Landau-Zener process. For longer strings, a more intricate spatiotemporal pattern emerges: the string breaks by forming a superposition of bubbles (domains of flipped spins of varying sizes), which involve highly excited states. We finally demonstrate that string breaking driven only by quantum fluctuations can be realized in the presence of sufficiently long-ranged interactions. This work holds immediate relevance for studying string breaking in quantum-simulation experiments.

Ising model

Bi‐Continuous W‐Rich Refractory High Entropy Alloy‐Cu Composite: Toward Material Innovation of Nuclear Reactor Coolant System

Abstract Refractory high‐entropy alloys (RHEAs) are considered promising candidate materials for next‐generation nuclear reactors due to their superior mechanical strength, irradiation resistance, and thermal stability at high temperatures. However, the significant positive heat of mixing between refractory alloying elements and Cu, commonly used in cooling systems, poses challenges in forming composite structures. This study addresses the issue using a liquid metal dealloying (LMD) process. A precursor alloy (WTaVTi) with a directional dendrite‐interdendrite structure is fabricated and reacted with molten Cu at 1200 °C for 96 h. This approach produced a RHEA‐Cu composite with a stable interface between RHEA (W 31.5 Ta 30.9 V 21.4 Ti 14.3 ) and Cu, featuring a spontaneously formed W‐rich interlayer that enhances interfacial bonding. The composite showed excellent irradiation resistance, with 30% less swelling under α‐ion irradiation than pure W. It also exhibited low thermal conductivity at room temperature, but reached ≈120 W m −1 ·K −1 at ≈650 °C, surpassing pure W. This temperature‐dependent rise in κ, with a positive gradient of +0.075 W m −1 ·K − 2 , is attributed to decreasing diffuse mismatch at elevated temperatures. The large‐scale reaction and stable microstructure achieved through LMD process highlight its industrial potential. This work offers a strategy for developing high‐performance materials by combining RHEA's radiation resistance with Cu's thermal conductivity for extreme environments.

Chemistry

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