Results for “Finite-Element” — Page 3

DOE OSTI2025

Addendum to SAND2023-09604 Xyce lumped-element transmission line model verification to support Empire-Cable cable SGEMP analyses

This report supplements the Verification of Empire-Cable SAND report by expanding on the use of Xyce to simulate the coupling to a transmission line cable model. While Empire-Cable solves its governing equations on a high-order, finite-element mesh with an an implicit-in-time formulation, Xyce must use a first order graph for the circuit and explicit-in-time approach to be compatible with non-linear electrical device models. Thus, given the different solution methodologies in Xyce as compared to Empire-Cable, the convergence rates are expected to be different but the overall quality of the solution should be the same. The original four canonical problems studied in the Empire-Cable verification report are replicated here running in Xyce using transmission line modeling parameters from the verification report. Overall, agreement between the codes is excellent with Xyce’s convergence rates being limited mostly to first order due to the circuit network approximation of a transmission line being a first order approximation.

42 ENGINEERING↗

DOE OSTI2026

Equivalent equilibrium ion distributions in ion sources and beams

Ion sources and ion beams often contain ions of more than one charge-to-mass ratio, due to the presence of multiple ionization states and/or multiple atomic and molecular ion species in the source. We show that, under certain conditions, analysis of these systems can be simplified using the notion of “equivalent equilibrium ion distributions” (EEIDs). These distributions have identical ion trajectories, space-charge distributions, and electric-field distributions, but can be composed of ion species with arbitrary charge-to-mass ratios. We use this result to obtain multi-species generalizations for the perveance of an ion diode, which relates the diode voltage to the extracted beam current, and for the perveance of a drifting ion beam, which determines the axial length-scale for beam expansion due to the self-electric field. Applied to a quasi-neutral plasma ion source containing cold ions and hot electrons, we find that the location of the beam-plasma interface (meniscus), where ions are extracted from the plasma, depends on the electron density and ion energy per unit charge at the interface, independent of the ion composition. For certain types of numerical simulations, EEIDs allow one to calculate the beam optics of an intense multi-species ion source by doing a simulation for a single ion species with an artificial charge-to-mass ratio chosen to minimize computation time. We illustrate these results using the 3D particle-in-cell, finite-element code Aleph. The simulations show how ions in a plasma source self-consistently accelerate to the supersonic speed needed for steady, radially expanding flow. Multi-species equilibria are subject to streaming instabilities because, while the different species have the same kinetic energy per unit charge, their velocities scale with the square root of the charge-to-mass ratio. The simulations show that a strong ion streaming instability can occur in a multi-species plasma ion source.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2026

Dynamic Multipoles in the RCS

The Rapid Cycling Synchrotron (RCS) of the Electron--Ion Collider (EIC) will accelerate polarized electrons up to 18 GeV at a repetition rate of 1 Hz. Dynamic multipoles driven by eddy currents in the vacuum chamber can affect beam dynamics and polarization preservation during the energy ramp. Building on the analytical treatment of S.Y. Lee for a ramped normal-conducting dipole magnet, we use conformal transformations to derive compact expressions for dynamic multipoles in higher-order ramped magnets with circular beam pipes. We further derive analytical frequency-domain transfer functions for both self-response and cross-response multipoles, providing compact dynamical models for their evolution during arbitrary ramp waveforms. The resulting analytical predictions are benchmarked against finite-element simulations, including the effects of practical pole geometry and beam-pipe misalignment. Finally, the expected dynamic multipoles in the RCS magnets are estimated and their implications for beam dynamics and polarization preservation are discussed.

43 PARTICLE ACCELERATORS↗

DOE OSTI2026

Estimation of Fission Product Transport Parameters for Cesium in the AGR-3/4 TRISO Fuel Experiment

A one-dimensional (1D) finite-element model of fission product transport in the AGR-3/4 experiment has been developed using the Multiphysics Object Oriented Simulation Environment (MOOSE) framework and implemented in the fuel performance code, BISON. The model resolves capsule-specific geometries, materials, and temperature histories and simulates radial migration of fission products from the fuel compact through the inner ring, outer ring, and into the sink ring. Model parameters governing diffusion and sorption were estimated for key fission products – cesium (Cs), and europium (Eu) – by simultaneously fitting modeled isotopic concentration profiles and total ring inventories to a post-irradiation experimental measurement. These data include gamma scanning, liquid scintillation for Sr-90, radial deconsolidation leach-burn-leach analysis, tomographic reconstructions, and destructive physical sampling. A mortar-based interfacial sorption framework was implemented to enforce physically consistent mass transfer and flux conservation across gas gaps. Two classes of parameter sets were derived: a least-squares best-fit, and a safety-oriented conservative-fit, what applies strong penalties for underprediction of sink inventories. Across all twelve capsules, the model successfully reproduces the dominant radial transport trends for Cs, Sr, with decreasing concentrations from the compact outward through successive rings. Cs behavior is captured most consistently, while strontium predictions reveal systematic trade-offs between compact accuracy and conservative sink-ring bounding. The results demonstrate that sink ring weighted calibration provides conservative, safety relevant bounds on low temperature fission product transport, but at the cost of underpredicting compact inventories for Sr isotopes. These discrepancies highlight the need for additional physics, including fast-slow diffusion model, incorporating trapping mechanism in the transport behavior. Overall, this work establishes a robust, capsule-specific modeling framework for AGR-3/4 fission product transport and provides a defensible basis for parameter selection in source-term and fuel performance analyses for high temperature gas-cooled reactors.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

Determining Optimal Magnetometer Configuration on MAGIS-100

Long-baseline atom interferometers such as the Matter-wave Atomic Gradiometer Interferometric Sensor (MAGIS-100) require stringent control and continuous characterization of background magnetic fields and spatial gradients to prevent systemic phase shifts that mimic ultralight dark matter or gravitational wave signatures. Because direct sensor placement within the ultra-high vacuum beam pipe is infeasible, in-situ magnetic field monitoring relies on external sensor arrays situated in the surrounding annular region. This work demonstrates a field reconstruction framework for a 5.3-meter MAGIS-100 modular section using finite-element Opera simulations. Transverse magnetic fields are expanded using a cylindrical multipole framework as informed by Fermilab’s Muon g-2 experiment, with magnetometer array configurations optimized via Fisher information matrix D-optimality. Inverting external sensor readings through a Gauss-Newton scheme recovers interior tube fields across distinct axial positions. In the discontinuity-averse uniform region (slice pair P4), the model achieves sub-noise-floor performance with a cross-validated root-mean-square error (RMSE) of $6.7227 \times 10^{-4}\text{ A/m}$ ($0.845\times$ sensor noise floor) and an interior field coefficient of variation of $1.71\%$. An elbow criterion in the Fisher bounds establishes $n_{\text{max}} = 2$ as the optimal multipole truncation order to prevent noise amplification from over-parameterization, with $n_{\text{max}} = 3$ (sextupole) order chosen for analysis to demonstrate further complexity and cross-pair comparison. Furthermore, analytical differentiation of the fitted multipole coefficients yields dense spatial maps of the transverse Jacobian gradient matrix $\nabla \mathbf{H}$ along with propagated $1\sigma$ uncertainty bounds across the beam region ($r \le 2.75\text{ in}$). This operational framework confirms that external magnetometer arrays can reliably monitor magnetic field uniformity and spatial gradients along the 100-meter flight path given appropriate sampling for any complexity order.

Appleby, Darwin [William Rainey Harper Coll.] (ORC↗

DOE OSTI2026

Determining Optimal Magnetometer Configuration on MAGIS-100

Long-baseline atom interferometers such as the Matter-wave Atomic Gradiometer Interferometric Sensor (MAGIS-100) require stringent control and continuous characterization of background magnetic fields and spatial gradients to prevent systemic phase shifts that mimic ultralight dark matter or gravitational wave signatures. Because direct sensor placement within the ultra-high vacuum beam pipe is infeasible, in-situ magnetic field monitoring relies on external sensor arrays situated in the surrounding annular region. This work demonstrates a field reconstruction framework for a 5.3-meter MAGIS-100 modular section using finite-element Opera simulations. Transverse magnetic fields are expanded using a cylindrical multipole framework as informed by Fermilab’s Muon g-2 experiment, with magnetometer array configurations optimized via Fisher information matrix D-optimality. Inverting external sensor readings through a Gauss-Newton scheme recovers interior tube fields across distinct axial positions. In the discontinuity-averse uniform region (slice pair P4), the model achieves sub-noise-floor performance with a cross-validated root-mean-square error (RMSE) of $6.7227 \times 10^{-4}\text{ A/m}$ ($0.845\times$ sensor noise floor) and an interior field coefficient of variation of $1.71\%$. An elbow criterion in the Fisher bounds establishes $n_{\text{max}} = 2$ as the optimal multipole truncation order to prevent noise amplification from over-parameterization, with $n_{\text{max}} = 3$ (sextupole) order chosen for analysis to demonstrate further complexity and cross-pair comparison. Furthermore, analytical differentiation of the fitted multipole coefficients yields dense spatial maps of the transverse Jacobian gradient matrix $\nabla \mathbf{H}$ along with propagated $1\sigma$ uncertainty bounds across the beam region ($r \le 2.75\text{ in}$). This operational framework confirms that external magnetometer arrays can reliably monitor magnetic field uniformity and spatial gradients along the 100-meter flight path given appropriate sampling for any complexity order.

Appleby, Darwin [William Rainey Harper Coll.] (ORC↗

DOE OSTI2026

Specimen sizing and remaining life sample calculations for high temperature reactor material surveillance

Advanced high-temperature nuclear reactors rely on structural components that will operate for decades under combined mechanical, thermal, and environmental loading. Materials surveillance programs are a promising strategy for managing the resulting uncertainty in long-term structural integrity by monitoring degradation in service using passively actuated mechanical test articles. Previous reports have developed a simplified, spreadsheetimplementable framework for sizing these test articles and for inferring accumulated creep damage and remaining life from ex-situ test data. This report advances that work toward practical deployment by providing a sample-problem book: a collection of worked, end-toend examples in which the ASME Section III, Division 5 design analysis of a representative high-temperature reactor component is carried through to a sized surveillance article, verified with detailed finite-element analysis, and concluded with a remaining-life assessment based on an assumed ex-situ creep-rate measurement on the retrieved specimen. The report also summarizes ongoing ANL engagement with ASTM Committee E10 on Nuclear Technology and Applications toward drafting a standard covering surveillance procedures for advanced reactors.

Barua, Bipul (ORCID:0000000247184113)↗

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 OSTI2025

Hybrid Basis and Multi-Center Grid Method for Strong-Field Processes

We present a time-dependent framework that combines a hybrid basis, consisting of Gaussian-type orbitals (GTOs) and finite-element discrete-variable representation (FEDVR) functions, with a multicenter grid to simulate strong-field and attosecond dynamics in atoms and molecules. The method incorporates the construction of the orthonormal hybrid basis, the evaluation of electronic integrals, a unitary time-propagation scheme, and the extraction of optical and photoelectron observables. Its accuracy and robustness are benchmarked on one-electron systems such as atomic hydrogen and the dihydrogen cation (H$^+_2$) through comparisons with essentially-exact reference results for bound-state energies, high-harmonic generation spectra, photoionization cross sections, and photoelectron momentum distributions. This work establishes the groundwork for its integration with quantum-chemistry methods, which is already operational but will be detailed in future work, thereby enabling ab initio simulations of correlated polyatomic systems in intense ultrafast laser fields.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2025

Effect of sintering temperature on adhesion of spray-on piezoelectric transducers

Conventionally sol-gel spray-on transducers require a high-temperature (> 700 ◦C) sintering process; however, this process can affect the microstructure of the substrate material. For mechanical elbows and valves utilized for fluid transport in the energy sector, the components are designed to have a specific microstructure, and deviations from these specifications can create weak points in the system. For this reason it is important to investigate how the temperature of the deposition process affects the substrate. This paper investigates the effect of high-temperature and low-temperature (< 150 ◦C) processing conditions on the surface composition of the substrate. Furthermore, the resultant transducers from high- and low-temperature fabrication processes are compared to determine if a low-temperature processing method is feasible. For these studies a sol-gel spray-on process is employed to deposit piezoelectric ceramics onto a stainless-steel 316L substrate. Energy-dispersive X-ray spectroscopy is utilized to determine the composition of the substrate surface before and after transducer deposition. Results indicate that the high-temperature processing conditions may alter the surface composition of the metal due to a diffusion of the metal into the ceramic, which results in a metal surface that is bonded to the ceramic. Furthermore, it is shown that low-temperature processing of spray-on transducers is a viable method for transducer fabrication where the resultant transducers meet the industry minimum requirement of 30 dB signalto-noise ratio. In parallel simulation calculations, finite-element method (FEM) studies were performed to model the adhesive strength of the low-temperature processed transducer to the substrate surface. Comparisons between the simulations and experiments suggest that the bond strength is much greater than the commercial gel bonds and closer to hardened epoxy glue bonds. These results indicate that spray-on transducers fabricated under lowtemperature processing conditions are a viable solution for leave-in-place monitoring of structures.

M. Sinding, Kyle↗

DOE OSTI2025

Microwave Debinding of Ceramics Produced by Additive Manufacturing

Microwave-assisted debinding offers an innovative approach to processing 3D-printed engineering ceramics. Ceramic parts produced with photopolymer-based additive manufacturing techniques require thermal debinding to remove polymers that bind ceramic particles during the printing process. Microwave-assisted debinding significantly reduces processing times and offers energy savings over conventional thermal treatments due to higher heating rates and more uniform heating. Finite-element modeling of the debinding process allows for the optimization and prediction of temperature distribution, stress development, and potential defects, leading to improved process control and material quality.

36 - MATERIALS SCIENCE↗

DOE OSTI2026

Broad frequency tuning of a Nb$_{3}$Sn superconducting microwave cavity for dark matter searches

We demonstrate a novel broad-frequency tuning mechanism for superconducting microwave cavities designed for dark matter searches. Using a Nb$_3$Sn-coated cigar-shaped cavity operating at approximately 9 GHz, we achieve continuous frequency tuning exceeding 1 GHz by mechanically separating the two cavity halves: a "tuning-by-opening" technique. Finite-element method simulations predict that radiative losses do not degrade the quality factor even for large openings, as a closed cavity with an intrinsic quality factor of $10^7$ maintains this value for apertures up to 9 mm, corresponding to a tuning range from 9.0 to 7.5 GHz. Experimental validation using both copper ring spacers and a continuous sliding mechanism confirms $Q_0$ values exceeding the dark matter quality factor across the entire explored frequency range, despite mechanical imperfections and film non-uniformities. This tuning approach avoids inserting elements into the resonant volume, making it particularly suitable for high-Q superconducting cavities in axion haloscope experiments and readily applicable to REBCO-based implementations capable of operating in multi-tesla magnetic fields.

Maiello, D. [Padua U.; INFN, Padua] (ORCID:0009000↗

DOE OSTI2026

Optimal Hairpin Winding Configuration for EV Traction Motors to Enhance Active CMV Cancellation Capability of NPL.X Inverter

Active common-mode voltage (CMV) cancellation in neutral-point-less (NPL.X) inverter is an effective approach for reducing common-mode electromagnetic interference and alleviating voltage stress associated with high-speed switching in EV traction drives. However, the existing hairpin winding configuration often exhibits impedance imbalance due to slot leakage inductance and manufacturing constraints, limiting the effectiveness of active CMV cancellation. This paper investigates symmetric hairpin winding configurations for a dual-three-phase EV traction machine to improve impedance symmetry. Considering the manufacturing constraints of six-layer hairpin windings, four feasible symmetric winding configurations are analyzed. Two-dimensional finite-element analysis is performed to extract the impedance parameters of each configuration. The extracted impedance parameters are used to quantify the impedance mismatch between complementary phase pairs and evaluate its influence on active CMV cancellation. The results show that the fully interleaved winding configuration achieves perfectly balanced phase impedances with 0 % mismatch, whereas the mismatch in the other configurations ranges from approximately 2.4 % to 0.8 %. Consequently, the resulting total CMV overshoot is significantly suppressed to near zero voltage under high-voltage and high switching frequency operating conditions, demonstrating the importance of winding symmetry for maximizing the active CMV cancellation capability of the NPL.X inverter.

Lee, Kangbeen [Purdue Univ., West Lafayette, IN (U↗

DOE OSTI2026

In-Situ Magnetic Field Reconstruction in the MAGIS-100 Experiment

Long-baseline atom interferometers such as the Matter-wave Atomic Gradiometer Interferometric Sensor (MAGIS-100) require stringent control and continuous characterization of background magnetic fields and spatial gradients to prevent systemic phase shifts that mimic ultralight dark matter or gravitational wave signatures. Because direct sensor placement within the ultra-high vacuum beam pipe is infeasible, in-situ magnetic field monitoring relies on external sensor arrays situated in the surrounding annular region. This work demonstrates a field reconstruction framework for a 5.3-meter MAGIS-100 modular section using finite-element Opera simulations. Transverse magnetic fields are expanded using a cylindrical multipole framework as informed by Fermilab’s Muon g-2 experiment, with magnetometer array configurations optimized via Fisher information matrix D-optimality. Inverting external sensor readings through a Gauss-Newton scheme recovers interior tube fields across distinct axial positions. In the discontinuity-averse uniform region (slice pair P4), the model achieves sub-noise-floor performance with a cross-validated root-mean-square error (RMSE) of $6.7227 \times 10^{-4}\text{ A/m}$ ($0.845\times$ sensor noise floor) and an interior field coefficient of variation of $1.71\%$. An elbow criterion in the Fisher bounds establishes $n_{\text{max}} = 2$ as the optimal multipole truncation order to prevent noise amplification from over-parameterization, with $n_{\text{max}} = 3$ (sextupole) order chosen for analysis to demonstrate further complexity and cross-pair comparison. Furthermore, analytical differentiation of the fitted multipole coefficients yields dense spatial maps of the transverse Jacobian gradient matrix $\nabla \mathbf{H}$ along with propagated $1\sigma$ uncertainty bounds across the beam region ($r \le 2.75\text{ in}$). This operational framework confirms that external magnetometer arrays can reliably monitor magnetic field uniformity and spatial gradients along the 100-meter flight path given appropriate sampling for any complexity order.

Appleby, Darwin [William Rainey Harper Coll.; Ferm↗

DOE OSTI2025

Physics-based stabilized finite element approximations of the Poisson–Nernst–Planck equations

We present and analyze two stabilized finite element methods for solving numerically the Poisson–Nernst–Planck equations. The stabilization we consider is carried out by using a shock detector and a discrete graph Laplacian operator for the ion equations, whereas the discrete equation for the electric potential need not be stabilized. Discrete solutions stemmed from the first algorithm preserve both maximum and minimum discrete principles. For the second algorithm, its discrete solutions are conceived so that they hold discrete principles and obey an entropy law provided that an acuteness condition is imposed for meshes. Remarkably the latter is found to be unconditionally stable. We validate our methodology through transient numerical experiments that show convergence toward steady-state solutions.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Polynomial range estimation as a troubled-cell indicator for high-order methods

Two troubled-cell indicators based on polynomial range estimation methods are used to flag cells that may violate positivity constraints. One method uses interval extension, and the second uses the range enclosure property of the Bernstein polynomial basis. Furthermore, both methods reduce compute time for the positivity preserver by limiting its application to a subset of cells. The Bernstein polynomial method remains effective as the problem dimensionality increases. Interval extension applied to the internal energy equation permits the use of the troubled-cell indicators for rational functions, though performance suffers compared to directly applying the indicators to polynomial functions.

42 ENGINEERING↗

DOE OSTI2025

A three-dimensional laser ray-tracing methodology for radiation-hydrodynamics simulations

We report on a methodology for performing laser ray-tracing in three spatial dimensions for radiation-hydrodynamics simulation codes. Our method, which is an extension of that developed in Haines et al., Comput. Fluids 201, 104478 (2020), utilizes an automatically generated separate mesh for the laser ray-tracing from the radiation-hydrodynamics mesh. This enables the laser mesh to be tailored to minimize ray noise with significantly fewer rays than would be required when the ray-tracing is performed on the radiation-hydrodynamics mesh, primarily by allowing the use of high-aspect-ratio cells that are not suitable for hydrodynamics solvers. For a planar target, we show that our method provides a ≈ 100× reduction in computational expense to achieve a fixed level of ray noise relative to ray-tracing directly on the radiation-hydrodynamics mesh. The relatively low ray requirement also enables efficient computation of cross-beam energy transfer. Each cell in the logically cubic laser mesh is a non-convex dodecahedron with triangular sides, and numerical integration of the ray trajectories and inverse bremsstrahlung is performed by mapping each cell to the unit cube. We will describe our methodology in detail as well as its implementation in the xRAGE radiation-hydrodynamics code, discuss performance, and present the results from applying the methodology to test problems with analytic solutions for laser ray-tracing through a quadratic density gradient with an analytic solution as well as for a laser-driven heat front. In 3D radiation-hydrodynamics simulations of laser-driven experiments performed on the National Ignition Facility, laser ray-tracing with our methodology uses less than 1% of total computational time while introducing acceptably low levels of ray noise.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Assessment of errors in analytic modeling of permanent magnet electron spectrometers for laser-plasma accelerators

Here, we assess the error in three treatments of a critical electron beam diagnostic for laser-plasma acceleration (LPA) experiments: a permanent magnet electron spectrometer. Since LPA electron spectrometers are often difficult to calibrate due to the scarcity of well-characterized, tunable electron beam sources in the appropriate energy range and the mechanical complexity of electron spectrometers, the standard of calibration is Hall probe measurements. We first compare the electron spectrometer performance between SIMION calculations and the Hall probe measurements. We find up to a 60% (<8%) error in determining the absolute energy for electrons below (above) 80 MeV when using the SIMION-modeled magnetic field vs the measured magnetic field due to SIMION overpredicting fringe field strength. The difference in spectrometer energy resolution is within ∼10% between the modeled and measured fields. We then assess a simple block model, which is commonly used. We find that the simple block model presented here sufficiently captures the predicted as-built spectrometer performance for any application provided that the user can tolerate uncertainties on absolute energy determination up to 5% and errors on energy resolution up to 1%.

Electromagnetic optics↗