Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite-Element”

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.

At least 19 records

SAVY-4000 Finite-Element Drop Test Analysis

PFE Auxiliary Systems conducted drop testing on SAVY-4000 containers to evaluate structural response under 12-foot drop conditions. In support of that effort, a finite-element modeling capability was developed to simulate drop response across multiple container sizes and impact orientations. The purpose of this work was to provide a consistent analysis framework that could support interpretation of testing, compare response trends across multiple configurations, and generate quantities of interest for later comparison with experimental data. More broadly, the analysis and testing were intended to assess whether the containers continued to perform their primary function after a 12-foot drop, namely maintaining structural integrity and containment of the contents. The modeling approach combined an implicit preload analysis with an explicit drop simulation so that each drop event began from a mechanically realistic assembled condition, including compression of the silicone O-ring. Separate models were developed for 2-quart, 5-quart, 12-quart, and 10-gallon containers. The results were evaluated in terms of strain-gauge response, collar-lid gap behavior, and accumulated plastic strain. In addition, parametric studies were performed on the 2-quart container to assess sensitivity to O-ring stiffness, friction, canister thickness, geometry tolerance, and mesh density. The simulations showed that predicted drop responses depended strongly on both container size and drop orientation. Gap metrics identified cases in which the predicted collar-lid opening exceeded the nominal O-ring cross-section threshold, while plastic strain metrics identified localized regions of elevated permanent deformation. Parametric studies showed that the predicted response was especially sensitive to the assumed O-ring stiffness and contact friction, while the geometry tolerance study produced smaller changes in the cases examined. The main value of this work was that it established a repeatable modeling and simulation workflow to support drop-test implementation, evaluate effects of future configuration changes, and understand modeling assumptions that most influenced predicted response. At the current stage, the results were viewed as preliminary model predictions rather than validated predictions. The next step would be to compare drop-test data to the model so that predictive values of the workflow could be refined and used with greater confidence to assess whether the containers maintained structural integrity and containment of the contents after a 12-foot drop.

42 ENGINEERING↗

Benchmarking of massively parallel phase-field codes for directional solidification

We present a detailed benchmark comparing two state-of-the-art phase-field implementations for simulating alloy solidification under experimentally relevant conditions. The study investigates the directional solidification of Al-3wt%Cu under high-velocity solidification conditions and SCN-0.46wt% camphor under microgravity conditions from National Aeronautics and Space Administration (NASA) DECLIC-DSI-R experiments. Both codes, one employing finite-difference discretization with uniform mesh and GPU-acceleration (GPU-PF) and the other one employing finite-element discretization with adaptive-mesh and CPU-parallelization (PRISMS-PF), solve the same quantitative phase-field formulation that incorporates an anti-trapping current for the solidification of dilute alloys. We evaluate the predictions of each code for dendritic morphology, primary spacing, and tip dynamics in both 2D and 3D, as well as their numerical convergence and computational performance. While existing benchmark problems have primarily focused on simplified or small-scale simulations, they do not reflect the computational and modeling challenges posed by employing experimentally relevant time and length scales. Our results provide a practical framework for assessing phase-field code performance as well as validating and facilitating their application in integrated computational materials engineering (ICME) workflows that require integration with realistic experimental data.

36 MATERIALS SCIENCE↗

Material-dependent photon ionizing radiation effects in Si and GaAs PIN diodes: A numerical investigation

We present a finite-element drift-diffusion-Poisson model in the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework to compare the radiation response of silicon (Si) and gallium arsenide (GaAs) PIN diodes under high-energy photon irradiation. The model solves coupled carrier continuity and Poisson’s equations with Shockley-Read-Hall recombination, and is verified against standard analytical J-V behavior. Using a simplified 1D geometry with ideal Ohmic contacts, we quantify device response under forward and reverse bias with a 100 MeV photon flux. Under forward bias, Si exhibits markedly greater radiation sensitivity than GaAs, including larger increases in current density, stronger local field and carrier-product perturbations, and higher recombination. Under reverse bias, GaAs shows larger radiation-induced photocurrent and broader current-density peaks near junctions, indicating an advantage for photodetection. Integrated steady-state recombination is consistently higher in Si across voltages. Under periodic photon pulses, GaAs produces higher-amplitude photoresponse and settles more rapidly than Si. These results highlight material-dependent trade-offs for radiation-tolerant, high-speed optoelectronics and provide guidance for selecting PIN architectures in aerospace, nuclear, and high-energy physics environments.

36 MATERIALS SCIENCE↗

Single-particle detection of enhanced polarizability in Au-decorated semiconducting nanorods via scanning dielectric microscopy

Hybrid nanostructures that combine semiconducting and metallic components offer great potential for photothermal therapy, optoelectronics, and sensing, by integrating tunable optical properties with enhanced light absorption and charge transport. Boosting the integrated performance of these hybrid systems demands techniques capable of probing local variations of the physical properties inaccessible to bulk analysis. Here, we report the single-particle dielectric characterization of hybrid, semiconducting bismuth sulfide (Bi 2 S 3 ) nanorods (NR) decorated with metallic Au nanoparticles (NP), employing scanning dielectric microscopy, which uses electrostatic force microscopy in combination with finite-element numerical simulations. We reveal a pronounced enhancement in the local dielectric response of Bi2S3 upon Au decoration, attributed to interfacial polarization and electron transfer from Au to the Bi 2 S 3 matrix, thus suggesting a enhanced metallic-like polarizability at the single-particle level. Numerical simulations show that the response is dominated by the vertical component of the permittivity and that the decorating metallic Au NP produce only moderate shielding of the semiconductor Bi 2 S 3 NR core, indicating that the large increase in the dielectric response originates primarily from intrinsic modifications within the NR. Overall, these findings provide direct insight into structure–property relationships at the single-particle level, supporting the rational design of advanced hybrid nanostructures with tailored electronic functionalities.

36 MATERIALS SCIENCE↗

Nonlinear electrical transport near the metal–insulator transition in V 4 O 7 thin films

Nonlinear electrical transport associated with correlated electronic states has been widely investigated in transition-metal oxides near metal–insulator transitions. Here, in this study, we investigate nonlinear transport and threshold switching in sputter-deposited V 4 O 7 thin films grown on fused silica substrates. Temperature-dependent transport measurements show a metal–insulator transition near 240 K with negligible thermal hysteresis, defining the temperature scale that governs the nonlinear electrical response. Current–voltage measurements reveal reproducible bipolar threshold switching over a broad temperature range, with the threshold voltage decreasing systematically as the transition temperature is approached. Electrothermal finite-element simulations reproduce the measured switching characteristics and show that switching occurs when localized Joule heating drives a confined region of the device toward the transition temperature. The simulated peak local temperature at threshold lies close to the intrinsic metal–insulator transition temperature of V 4 O 7 , establishing that the nonlinear electrical response originates from electrothermal feedback acting on the strongly temperature-dependent conductivity near the transition. These results identify V 4 O 7 thin films as a model correlated-oxide system in which volatile threshold switching emerges from proximity to an extended metal–insulator transition.

25 ENERGY STORAGE↗

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↗

Comparison between Alcator C-Mod ICRF experiments and 3D full wave simulations

Reliable modeling of ion cyclotron range of frequencies (ICRF) antenna performance is essential for interpreting present experiments and guiding the design of future reactors. In this work, a 3D model of the Alcator C-Mod field-aligned antenna is implemented in the Petra-M finite-element framework [S. Shiraiwa et al 2023 Nucl. Fusion 63 026024] and benchmarked against experimental results. Four experimental cases are examined. First, the simulated rectified sheath potentials on the antenna limiters are compared with measurements from a power tapering experiment. Second, proof-of-principle far-field sheath simulations are performed. In scenarios with low single-pass absorption, simulations predict enhanced sheath potentials on a distant poloidal limiter in the far field of the antenna, consistent with experimental observations. Third, the simulated antenna loading during edge localized modes follows experimental trends and appears to be dominated by the density gradient at the pedestal. Finally, Petra-M predicts the unintended excitation of high-$k$ ∥ modes during monopole phasing operation, in agreement with experimental evidence of poor wave coupling and accessibility to the plasma core. Overall, reasonable agreement is found between the simulations and experiments. At the same time, areas of imperfect agreement are identified. These provide important guidance on the limits of the current state-of-the-art modeling, which should be kept in mind when using it as a predictive tool for future reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Systematic Construction of Time-Dependent Hamiltonians for Microwave-Driven Josephson Circuits

Time-dependent electromagnetic drives are fundamental for controlling complex quantum systems, including superconducting Josephson circuits. In these devices, accurate time-dependent Hamiltonian models are imperative for predicting their dynamics and designing high-fidelity quantum operations. Existing numerical methods, such as black-box quantization (BBQ) and energy-participation ratio (EPR), excel at modeling the static Hamiltonians of Josephson circuits. However, these techniques do not fully capture the behavior of driven circuits stimulated by external microwave drives, nor do they include a generalized approach to account for the inevitable noise and dissipation that enter through microwave ports. Here, we introduce numerical techniques that leverage classical microwave simulations, efficiently executable in finite-element solvers, to obtain the time-dependent Hamiltonian of microwave-driven superconducting circuits with arbitrary geometries under charge, flux, or mixed electromagnetic modulation. Importantly, our techniques do not rely on a lumped-element description of the superconducting circuit, in contrast to previous approaches to tackling this problem. We demonstrate the versatility of our approach by characterizing the driven properties of realistic circuit devices in complex electromagnetic environments, including coherent dynamics due to charge and flux modulation, as well as drive-induced relaxation and dephasing. Our techniques offer a powerful toolbox for optimizing circuit designs and advancing practical applications in superconducting quantum computing.

Lu, Yao [Yale U.; Yale U. (main); Fermilab] (ORCID↗

Integrated Magnetics for Multiport Inductive and Conductive Power Transfer Architecture for Electric Vehicles

This paper proposes a shared magnetic arrangement for an electric-vehicle charger that supports both conductive and inductive charging while preserving independent power control at the two ports. In the proposed implementation, the high-frequency transformer used for the conductive path and the resonant inductor required by the inductive path are realized in an integrated form. By combining the magnetic structure and primary-side converter hardware, the charger can deliver improved hardware utilization and higher power density without a corresponding increase in cost or implementation complexity. The paper describes the magnetic realization and presents finite-element and circuit-level simulation results to confirm minimal interaction between the two charging paths. The study demonstrates that the proposed charger is a viable option for flexible, high-power EV charging systems.

Mukherjee, Subho [ORNL] (ORCID:0009000672297925)↗

A Unified Off-Board Charger for Three-Phase Inductive and Single-Phase Conductive EV Charging

This paper presents a unified off-board charger integrating a three-phase inductive power transfer (IPT) system and a single-phase conductive (plug-in) charger for electric vehicles (EVs). The proposed charger shares a high-frequency inverter and an integrated magnetic structure that enables three operating modes: conductive charging, inductive charging, and simultaneous charging. The transmitter coils for the three-phase IPT and the high-frequency transformer windings for the conductive charger are arranged on a common ferrite pad while maintaining magnetic decoupling between the two power transfer paths. The proposed unified architecture effectively increases charger utilization and power density without a proportional increase in cost, weight, or volume. Finite-element analysis using Ansys Maxwell confirms negligible coupling between the inductive and conductive coils, and the extracted parameters are employed in PLECS simulations to verify the independent operation of each charging mode.

Jo, Cheolhui [ORNL] (ORCID:0000000322692434)↗

A Unified Wireless Charger, On-Board Charger, and Auxiliary Power Module for Electric Vehicle Charging Systems

This paper proposes a unified electric vehicle (EV) charging architecture that integrates wireless power transfer (WPT), an on-board charger (OBC), and an auxiliary power module (APM) within a single architecture. By sharing a multi-functional magnetic structure and active switch bridges, the proposed topology eliminates additional transformers and converter stages, reducing hardware complexity and improving power density. A multipurpose magnetic design achieves magnetic decoupling among the WPT, OBC, and APM functions while maintaining the required coupling for each mode. Through electrical reconfiguration, the WPT operates as an LCC-S converter, whereas the OBC and APM operate as dual-active-bridge (DAB) converters. The system supports multiple operating modes, including simultaneous high-voltage and lowvoltage battery charging. Finite-element and circuit simulations verify the magnetic characteristics and system operation, demonstrating the feasibility of the proposed unified architecture for EV charging applications.

Jo, Cheolhui [ORNL] (ORCID:0000000322692434)↗

A Unified Design Theory for Multi-Port Polyphase Transformers Enabling Scalable Power-Multiplexed EV Fleet Charging Systems

This paper presents a unified analytical design theory for multi-port polyphase transformers, targeting scalable and isolated high-power Electric Vehicle (EV) fleet charging systems with power multiplexing capability. As fleet electrification accelerates, conventional one-to-one charger architectures face significant challenges in infrastructure cost, peak power demand, and low utilization of installed power electronics. Power-multiplexed charging architectures, which dynamically distribute power from a shared pool of converter modules across multiple vehicles, have emerged as a promising solution. However, such architectures require scalable, isolated multi-port power interfaces capable of routing energy among multiple inputs and outputs, whose design remains complex and dependent on iterative modeling. To address this gap, the proposed theory provides closed-form expressions for self-inductance, leakage inductance, and mutual coupling terms for arbitrary multi-phase, multi-port transformer structures. The formulation enables direct synthesis of isolated multi-input and multi-output resonant converter systems without reliance on geometry-specific finite-element analysis or extensive parameter extraction. This capability is particularly critical for power-multiplexed systems, where modular converter structures must interface with multiple vehicles while maintaining galvanic isolation and flexible power allocation. The effectiveness of the proposed framework is demonstrated through the design of a 360 kW multi-phase system operating over a 700–900 VDC input and 400–1250 VDC output range. PLECS simulation results confirm accurate prediction of system behavior and validate the applicability of the approach to multi-port, power-multiplexed charging scenarios. The proposed method significantly reduces design complexity while enabling scalable, cost-effective, and fully utilized EV fleet charging infrastructure.

Asa, Erdem [ORNL] (ORCID:0000000190884812)↗

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↗

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↗

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↗

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↗

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)↗

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↗