Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite element method (FEM)”

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.

424 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A Full-Induction Magnetohydrodynamics Solver for Liquid Metal Fusion Blankets in Vertex-CFD

Multiphysics modeling of liquid metal fusion blankets, which produce tritium and convert energy of neutrons created via fusion reactions into heat, is crucial for predicting performance, ensuring structural integrity, and optimizing energy production. While traditional blanket modeling of liquid metal flows during normal steady operating conditions commonly employs the inductionless approximation of the magnetohydrodynamics (MHD) equations, transient scenarios, when the plasma-confining magnetic field varies on millisecond time scales, require a full-induction MHD approach that dynamically evolves the magnetic field via the time-dependent induction equation. This paper presents the formulation, implementation, and initial verification of a full-induction MHD solver integrated within the open-source Vertex-CFD framework, which aims to achieve tight multiphysics coupling, a flexible software design enabling easy extension and addition of physics models, and performance portability across computing platforms. The solver utilizes finite element spatial discretization, implicit Runge–Kutta time integration, and an inexact Newton method to solve the resulting discrete nonlinear system, leveraging Trilinos packages for efficient computation. Verification against selected benchmark problems demonstrates accuracy and robustness of the solver. Furthermore, when the solver is applied to an idealized blanket model in 2.5D and full 3D, results obtained with Vertex-CFD are in good agreement with recently published quasi-2D simulations. These findings establish a computational foundation for future simulations of transient MHD phenomena in liquid metal blankets with Vertex-CFD, and open avenues for future extensions and performance optimizations.

Endeve, Eirik [ORNL] (ORCID:0000000312519507)↗
DOE OSTI2026

Bacterial Bioleaching and Biorecovery for Biomining Unconventional Rare Earth Element Feedstocks

Bacterial metabolic interactions with rare earth elements (REEs) can be harnessed for biomining unconventional feedstocks like abandoned coal-mine drainage (AMD). Pennsylvania has ~500 AMD passive remediation systems that can precipitate REE rich solids. REEs include yttrium and the lanthanide series that are used in modern energy and technology. Bacteria that metabolically interact with REEs can be used for biomining in an affordable efficient process that does not require hazardous chemical additives. Currently, the microbial metal mechanisms that contribute to REE biorelease and biorecovery are poorly understood. Our work shows acidogenic bacterial isolates (Bacillus mycoides JR07 and Bacillus pseudomycoides KB7) successfully bioleach a mixed REE solution from AMD solids by their organic acid production and biofilm formation. Further, our work shows the potential for bacterial lanthanide-dependent enzymes to recover lanthanides from a mixed REE solution; here we have bacterial isolate Methylobacterium sp. B3 that can recover soluble lanthanum. Whole genome sequencing of Methylobacterium sp. B3 predict lanthanide-dependent methanol dehydrogenase XoxF. Understanding the microbial metabolism and genes involved in the REE release and recovery is crucial to optimize the biomining of AMD solids. Our work addresses the growing need to develop novel REE mining methods from unconventional feedstocks.

biogeochemistry↗
DOE OSTI2024

Phonon screening and dissociation of excitons at finite temperatures from first principles

The properties of excitons, or correlated electron–hole pairs, are of paramount importance to optoelectronic applications of materials. A central component of exciton physics is the electron–hole interaction, which is commonly treated as screened solely by electrons within a material. However, nuclear motion can screen this Coulomb interaction as well, with several recent studies developing model approaches for approximating the phonon screening of excitonic properties. While these model approaches tend to improve agreement with experiment, they rely on several approximations that restrict their applicability to a wide range of materials, and thus far they have neglected the effect of finite temperatures. Here, we develop a fully first-principles, parameter-free approach to compute the temperature-dependent effects of phonon screening within the ab initio GW -Bethe–Salpeter equation framework. We recover previously proposed models of phonon screening as well-defined limits of our general framework, and discuss their validity by comparing them against our first-principles results. We develop an efficient computational workflow and apply it to a diverse set of semiconductors, specifically AlN, CdS, GaN, MgO, and SrTiO 3 . We demonstrate under different physical scenarios how excitons may be screened by multiple polar optical or acoustic phonons, how their binding energies can exhibit strong temperature dependence, and the ultrafast timescales on which they dissociate into free electron–hole pairs.

Science & Technology - Other Topics↗
DOE OSTI2024

Determining Stellar Elemental Abundances from DESI Spectra with the Data-driven Payne

Abstract Stellar abundances for a large number of stars provide key information for the study of Galactic formation history. Large spectroscopic surveys such as the Dark Energy Spectroscopic Instrument (DESI) and LAMOST take median-to-low-resolution (R≲ 5000) spectra in the full optical wavelength range for millions of stars. However, the line-blending effect in these spectra causes great challenges for elemental abundance determination. Here we employDD-Payne, a data-driven method regularized by differential spectra from stellar physical models, to the DESI early data release spectra for stellar abundance determination. Our implementation delivers 15 labels, including effective temperatureT eff , surface gravity log g , microturbulence velocityv mic , and the abundances for 12 individual elements, namely C, N, O, Mg, Al, Si, Ca, Ti, Cr, Mn, Fe, and Ni. Given a spectral signal-to-noise ratio of 100 per pixel, the internal precisions of the label estimates are about 20 K forT eff , 0.05 dex for log g , and 0.05 dex for most elemental abundances. These results agree with the theoretical limits from the Crámer–Rao bound calculation within a factor of 2. The majority of the accreted halo stars contributed by the Gaia–Enceladus–Sausage are discernible from the disk and in situ halo populations in the resultant [Mg/Fe]–[Fe/H] and [Al/Fe]–[Fe/H] abundance spaces. We also provide distance and orbital parameters for the sample stars, which spread over a distance out to ∼100 kpc. The DESI sample has a significantly higher fraction of distant (or metal-poor) stars than the other existing spectroscopic surveys, making it a powerful data set for studying the Galactic outskirts. The catalog is publicly available.

Astronomy & Astrophysics↗
DOE OSTI2026

Design of monolithic piezoelectric bimorph mirrors made from lithium niobate

Recent advances in deformable mirrors based on monolithic piezoelectric substrates, such as lithium niobate, have the potential to improve the image quality of X-ray optical systems on synchrotron and free-electron laser beamlines and microscopes. However, the quantitative relationship between design parameters and the deformed shape has not been readily available in the literature. We present an analytical model, validated through finite element analysis, enabling calculation of tangential and sagittal curvatures based on mirror dimensions, crystallographic orientation, and applied voltage. We demonstrate that through the selection of material orientation, it is possible to achieve different deformed shapes (sphere, cylinder, or hyperbolic paraboloid). This methodology can be generalized to other piezoelectric materials and substrate-integrated actuator systems.

Marzari, Francesco [University of Trento (Italy); ↗
DOE OSTI2026

The critical role of intrinsic defects and many-body interactions on the stability of MnBi2Te4

Intrinsic antisite defects pose a major challenge to understanding and predicting the exotic properties of the layered topological magnetic insulator MnBi2Te4 (MBT). In this work, we study the origin of the abundance of intrinsic defects in MBT, including many-body defect–defect interactions and many-body electronic correlations. Until now, ab initio methods have struggled to explain thermodynamic stability and properties influenced by defect behavior in MBT. We model native Mn–Bi antisite defects in MBT at finite temperatures using a cluster expansion that includes defect–defect interactions. To overcome the limitations of conventional density functional theory (DFT), we introduce a hybrid approach that incorporates high-accuracy quantum Monte Carlo (QMC) calculations, introducing missing correlations. This strategy allows for accurate estimation of defect energetics and finite-temperature properties. We compute the configurational free energy, defect concentration, and configurational heat capacity, revealing a second-order order–disorder phase transition near the experimental synthesis temperature. Our study provides the first theoretical insight into the thermodynamics of intrinsic defects in MBT. The negative free energy relative to pristine MBT at synthesis temperatures indicates that Mn–Bi antisite formation is thermodynamically spontaneous. We also present a broadly applicable general framework for correcting low-level theoretical theories using highly accurate many-body corrections from QMC.

Ghaffar, Abdul [ORNL] (ORCID:0000000241190168)↗
DOE OSTI2025

The phase diagram of quantum chromodynamics in one dimension on a quantum computer

The quantum chromodynamics (QCD) phase diagram, which reveals the state of strongly interacting matter at different temperatures and densities, is key to answering open questions in physics, ranging from the behaviour of particles in neutron stars to the conditions of the early universe. However, classical simulations of QCD face significant computational barriers, such as the sign problem at finite matter densities. Quantum computing offers a promising solution to overcome these challenges. Here, we take an important step toward exploring the QCD phase diagram with quantum devices by preparing thermal states in one-dimensional non-Abelian gauge theories. We experimentally simulate the thermal states of SU(2) and SU(3) gauge theories at finite densities on a trapped-ion quantum computer using a variational method. This is achieved by introducing two features: Firstly, we add motional ancillae to the existing qubit register to efficiently prepare thermal probability distributions. Secondly, we introduce charge-singlet measurements to enforce colour-neutrality constraints. This work pioneers the quantum simulation of QCD at finite density and temperature for two and three colours, laying the foundation to explore QCD phenomena on quantum platforms.

Quantum information↗
DOE OSTI2026

A Novel Low-Profile High-Efficiency Three-Phase Matrix Transformer

High step-down isolated DC-DC conversion from an 800 V DC bus to low-voltage, high-current outputs is required in automotive auxiliary converters and data center power supplies. In such applications, conventional transformer-based converters require large turns ratios, which increase winding resistance, leakage inductance, and magnetic height. This paper proposes a novel low-profile three-phase matrix transformer that realizes a large effective voltage ratio through flux division among multiple secondary legs, without increasing the physical turns count of each winding. As a result, the proposed structure reduces copper usage and transformer height while preserving the voltage conversion capability of a conventional three-phase transformer. Finite element analysis shows that the proposed design reduces magnetic height by 27%, ferrite volume by 34%, and copper volume by 28%. Circuit-level simulations of an 800 V/12 V,3 kW CLLLC dual-active-bridge converter further show that the lower winding resistance reduces total system loss by 91% and increases DC-DC efficiency from 82.6% to 97.6% at 3 kW output.

Inoue, Shuntaro [ORNL] (ORCID:0000000262637627)↗
DOE OSTI1969

Application of Nuclear Technology to Art Identification Problems: First Annual Report

Throughout the centuries, as man has become increasingly affluent, his interest in antiquities and objects of art, and the price he has been willing to pay for them, has increased. Inevitably, one result has been to attract the unscrupulous to the production and sale of counterfeit paintings, sculptures, and other works of art. The skill of forgers varies but often is highly sophisticated. However, dealers, museum curators, and other experts have also become more skillful in the detection of forgeries. Scientific tools of increasing sensitivity and sophistication have gradually been applied to the examination of the materials of art and archaeology. Such tools frequently support and render more reliable the judgment of the experts. The quality of forgeries may improve further as their makers become acquainted with the new methods of detection and in turn learn to circumvent or confound the methods. Therefore, the development of still more advanced methods of examination has great utility in the art world. This is particularly true if the ultimate effect is to make circumvention of the methods of examination so costly that the economic incentives for producing forgeries may be substantially reduced, if not eliminated.

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

Predicting microstructurally sensitive fatigue‐crack path in WE43 magnesium using high‐fidelity numerical modeling and three‐dimensional experimental characterization

Abstract Microstructurally small fatigue‐crack growth in polycrystalline materials is highly three‐dimensional due to sensitivity to local microstructural features (e.g., grains). One requirement for modeling microstructurally sensitive crack propagation is establishing the criteria that govern crack evolution, including crack deflection. Here, a high‐fidelity finite‐element modeling framework is used to assess the performance and validity of various crack‐growth criteria, including slip‐based metrics (e.g., fatigue‐indicator parameters), as potential criteria for predicting three‐dimensional crack paths in polycrystalline materials. The modeling framework represents cracks as geometrically explicit discontinuities and involves voxel‐based remeshing, mesh‐gradation control, and a crystal‐plasticity constitutive model. The predictions are compared to experimental measurements of WE43 magnesium samples subject to fatigue loading, for which three‐dimensional grain structures and fatigue‐crack surfaces were measured post‐mortem using near‐field high‐energy x‐ray diffraction microscopy and x‐ray computed tomography. Findings from this work are expected to improve the predictive capabilities of simulations involving microstructurally small fatigue‐crack growth in polycrystalline materials.

Engineering↗
DOE OSTI2024

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↗
DOE OSTI2026

Numerical-heating effects in atmospheric pressure streamer discharges simulated with a PIC code

Artificial heating in plasma simulations is a well-known phenomenon which occurs when, among other things, the Debye length is poorly resolved by the simulation mesh. Here, in this work, the degree to which numerical-heating occurs during a simulation of a nanosecond atmospheric pressure streamer discharge is examined. The streamer is simulated using a two-dimensional finite-element, particle-in-cell code Empire, which uses direct simulation Monte Carlo for binary particle interactions. Initially, an estimate of the numerical-heating rate applied to Empire is performed using a simple plasma model. Second, a positive atmospheric pressure streamer discharge simulation is performed to study the effects of numerical heating on plasma density, electron temperature, and streamer velocity. The nominal Debye length is approximately 1 μm and the amount of numerical heating introduced in the simulation is varied by using mesh sizes ranging from 2 μm to 20 μm. A measurable numerical heating quantity is proposed that can be used to estimate the appropriate element size and quantify the numerical-heating that can be expected over the simulation time for an atmospheric pressure streamer. In conclusion while Δx/λ D violations can be an issue it is not likely to be an issue with streamer discharges that are temporally short and occur in environments where collision frequencies are high. This result validates the rationale of grid size choices for a large amount of previously published works where Δx/λ D violation was not clearly addressed. Primary finding of this work is that numerical heating is of minor concern for plasma simulations where electron–neutral collisions are numerous such that multiple collisions can occur within a single plasma period.

Nikic, Dejan [University of New Mexico, Albuquerqu↗
DOE OSTI2025

AL-CRADA-2025-01 Final Report: American-Made Wind Turbine Materials Recycling Prize Accelerate! Phase 2

This report summarizes the Phase 2 achievements of the Wind Turbine Materials Recycling Prize project led by Critical Materials Recycling (CMR) in partnership with Ames National Laboratory. The team demonstrated the technical feasibility and economic viability of recovering and reusing rare-earth Nd-Fe-B magnets from end-of-life (EOL) wind turbine generators. Through mechanical disassembly, demagnetization, and precision processing, magnets were harvested and characterized to assess structural, compositional, and magnetic properties. Results showed that recovered magnets performed on par with commercial-grade counterparts, enabling their integration into actuator and motor prototypes without redesign. System-level modeling and finite-element simulations validated operational performance, while Life Cycle Cost (LCC) analysis revealed up to 20% material cost savings and 2–3% reductions in overall motor system costs. Life Cycle Assessment (LCA) further confirmed substantial environmental benefits, including a ~94% reduction in global warming potential relative to virgin magnet production. The project’s success was enabled by a tightly integrated partnership between industry and national lab experts, illustrating a viable path for domestic magnet-to-magnet recycling. These findings support a scalable, circular solution for rare-earth material recovery, reducing dependence on imported critical minerals and advancing sustainable clean energy technologies.

Lograsso, Thomas [Ames Laboratory (AMES), Ames, IA↗
DOE OSTI2026

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↗
DOE OSTI2026

Evaluating Fracture Behavior of Bioinspired Alumina-YSZ Composites through Static and Dynamic Mechanical Testing

Ceramic materials are known for their high hardness and strength but are limited by low toughness and sudden failure. Inspired by the microstructure of dental enamel, which features undulating rods that promote crack deflection and energy absorption, architected specimens were designed and manufactured. A bespoke ball-on-ring (BoR) testing apparatus measured the static biaxial rupture strength of these bioinspired materials. Additionally, impact testing with spherical steel projectiles assessed their behavior under dynamic conditions. Monolithic alumina specimens were compared to architected alumina reinforced with yttria-stabilized zirconia (YSZ) rods produced via direct ink write 3D printing. BoR tests revealed that monolithic alumina had a Weibull modulus of 10.53 and a characteristic strength of 372.3 MPa, while the composite specimens showed a Weibull modulus of 5.46 and a characteristic strength of 213.3 MPa. Although the composite failed at lower stresses and projectile velocities, it exhibited notable crack deflection and fracture resistance, with cracks being effectively interrupted by the rod inclusions. The composite specimens also resulted in fewer and larger fragments upon impact. Finite element simulations confirmed the effectiveness of the rod structures in enhancing fracture resistance.

Bioinspired materials↗
DOE OSTI2026

SQUID Readout of a High-$Q$ Superconducting $LC$ Resonator

We demonstrate the readout of a superconducting $LC$ resonator, with $Q \approx 2 \times 10^{6}$ at approximately 250 kHz, using a dc SQUID followed by a SQUID series-array amplifier readout chain. We find that $Q$ depends on the SQUID flux-bias point, increasing near the shallow-slope point of the SQUID modulation curve and decreasing near the steep-slope point, consistent with the previously observed SQUID damping effects. From this variation, we infer the SQUID effective input impedance. We further infer the effective temperature of the resonator circuit from the resonance peak in the noise spectrum and show that it also depends on the SQUID flux-bias point, suggesting a possible contribution from SQUID back-action noise. This system provides a working prototype for future experiments based on lumped-element resonators with SQUID readout, in particular low-mass axion searches within the DMRadio program.

FOS: Physical sciences↗
DOE OSTI2024

Cosmic neutrino decoupling and its observable imprints: insights from entropic-dual transport

Abstract Very different processes characterize the decoupling of neutrinos to form the cosmic neutrino background (CνB) and the much later decoupling of photons from thermal equilibrium to form the cosmic microwave background (CMB). The CνB emerges from the fuzzy, energy-dependent neutrinosphere and encodes the physics operating in the early universe in the temperature rangeT∼ 10 MeV toT∼ 10 keV. This is the epoch where beyond Standard Model (BSM) physics, especially in the neutrino sector, may be influential in setting the light element abundances, the necessarily distorted fossil neutrino energy spectra, and other light particle energy density contributions. Here we use techniques honed in extensive CMB studies to analyze the CνB as calculated in detailed neutrino energy transport and nuclear reaction simulations of the protracted weak decoupling and primordial nucleosynthesis epochs. Our moment method, relative entropy, and differential visibility approach can leverage future high precision CMB and light element primordial abundance measurements to provide new insights into the CνB and any BSM physics it encodes. We demonstrate that the evolution of the energy spectrum of the CνB throughout the weak decoupling epoch is accurately captured in the Standard Model by only three parameters per species, a non-trivial conclusion given the deviation from thermal equilibrium and the impact of the decrease of electron-positron pairs. Furthermore, we can interpret each of the three parameters as physical characteristics of a non-equilibrium system. Though the treatment presented here makes some simplifying assumptions including ignoring neutrino flavor oscillations, the success of our compact description within the Standard Model motivates its use also in BSM scenarios. We further demonstrate how observations of primordial light element abundances can be used to place constraints on the CνB energy spectrum, deriving response functions that can be applied for general deviations from a thermal spectrum. Combined with the description of those deviations that we develop here, our methods provide a convenient and powerful framework to constrain the impact of BSM physics on the CνB.

Astronomy & Astrophysics↗