Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Magnetic and Mechanical Analysis of Bi-2212 Rutherford Cable in a Cos-Theta Sub-Scale Dipole Coil

The U.S. Magnet Development Program (US-MDP) explores high-field accelerator magnets compatible with operational conditions beyond the limits of Nb$_3$Sn technology. The ongoing R&D High-Temperature Superconductors (HTS) suggests using Bi$_2$Sr$_2$CaCu$_2$O$_{8-x}$ (Bi-2212) as superconducting element. Bi-2212 Rutherford cables maintain a high critical current (I$_C$) when exposed to a large external magnetic field. However, Bi-2212 exhibits an oversensitive stress-strain response when subject to large Lorentz forces. This paper reports on the magnetic and mechanical analysis of the Bi-2212 cosine-theta insert being developed at Fermilab for a hybrid magnet composed of two external layers of Nb$_3$Sn and two internal layers of Bi-2212. We performed a FEM analysis of the insert to estimate the HTS stress state in the coil's strands under magnetic and mechanical loads.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Exploring Structural Anisotropy in Amorphous Tb-Co via Changes in Medium-Range Ordering

Amorphous thin films grown by magnetron co-sputtering exhibit changes in atomic structure with varying growth and annealing temperatures. Structural variations influence the bulk properties of the films. Scanning nanodiffraction performed in a transmission electron microscope (TEM) is applied to amorphous Tb 17 Co 83 (a-Tb-Co) films deposited over a range of temperatures to measure relative changes in medium-range ordering (MRO). These measurements reveal an increase in MRO with higher growth temperatures and a decrease in MRO with higher annealing temperatures. The trend in MRO indicates a relationship between the growth conditions and local atomic ordering. By tilting select films, the TEM measures variations in the local atomic structure as a function of orientation within the films. The findings support claims that preferential ordering along the growth direction results from temperature-mediated adatom configurations during deposition, and that oriented MRO correlates with increased structural anisotropy, explaining the strong growth-induced perpendicular magnetic anisotropy found in rare earth–transition metal films. Beyond magnetic films, we propose the tilted FEM workflow as a method of extracting anisotropic structural information in a variety of amorphous materials with directionally dependent bulk properties, such as films with inherent bonding asymmetry grown by physical vapor deposition.

36 MATERIALS SCIENCE↗
DOE OSTI2026

The Use of Machine Learning Models for Predicting the Dielectric Strength of Gases

Technological advancements in high voltage systems have pushed sulfur hexafluoride (SF6) to its operational limits. Furthermore, this gas has other drawbacks including a high liquefaction temperature and a high global warming potential. Therefore, there has been an urgent need to find alternative gases with high dielectric strength (DS). In this work, density functional theory (DFT) is used to calculate molecular descriptors that are fed into an artificial neural network (ANN) and a random forest (RF). These machine learning (ML) models are then used to predict the DS for hundreds of molecules. A finite element model (FEM) is also used to calculate the electric field profile of multiple simple electrode geometries as the applied voltage to the system is increased. Results indicate that the random forest model has better generalization to unseen data than the neural network. The highest DS value predicted by the RF was 2.16 relative to the experimental DS of SF6. The results also demonstrate how choosing a gas with a higher DS and a geometry with minimal edges and corners can significantly increase the operating voltage of an electrical system. Due to its superior generalization, the RF represents the most promising path toward an accurate DS predictor once sufficient experimental data are available.

Mileski, Matthew [AFIT]↗
DOE OSTI2025

Structural Analysis of Multilayer Tapered Canted-Cosine-Theta Superconducting Magnet for the Interaction Region of Electron Ion Collider

The electron ion collider (EIC) will collide high energy and highly polarized hadron and electron beams with luminosities up to 10 34 cm -2 s -1 . The magnet designs at the interaction region (IR) are challenging due to the close proximity of the hadron and electron beams. A canted cosine theta (CCT) quadrupole magnet was designed due to the space restrictions at the IR. The CCT magnet has a tapered aperture, which allows higher gradients at the interaction point (IP) side of the magnet. The tapered geometry makes the magnet not symmetrical along the axial direction. In conclusion, to validate the design in terms of the mechanical strength and the allowed deformations, a detailed structural mechanical analysis for the tapered CCT quadrupole magnet has been carried out using the finite element method (FEM) considering both the shrinking due to cool down and the electromagnetic force on the conductors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Structural Analysis of a Serpentine Superconducting Magnet for the Interaction Region of the Electron Ion Collider

The electron ion collider (EIC) is under design to collide the high energy and highly polarized hadron with electron beams with luminosities up to 10 34 cm −2 s −1 . The superconducting magnet designs at the interaction region (IR) are quite challenging due to the close proximity of the hadron and electron beams. Several serpentine types of superconducting magnets have been designed due to the space restrictions at the IR. Here, to validate the designs in terms of the mechanical strength, structural integrity, as well as the allowed deformations, a detailed structural mechanical analysis for a 12-layer superconducting quadrupole serpentine magnet (Q1BpR) in the rear side hadron beams has been carried out using the finite element method (FEM) considering the pretension stress, shrinking due to cool down, and the electromagnetic force on the conductors. Anisotropy material properties and individual roving tension have been considered. The electro-magnetic (EM) simulations were performed using COMSOL, confirmed with RAT, the structural mechanical analysis were performed using COMSOL by considering the contact elements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Numerical Simulation of Irradiation Induced Swelling for STAR 4.1 Blanket

Here, this article presents the preliminary results of FEM implementation of irradiation-induced swelling, hardening, and creep effect using ANSYS user programmable features (UPFs) applied to material properties. The first liquid metal breeding blanket model for STAR 4.1 tokamak as a part of the virtual prototyping system is under development at Princeton Plasma Physics Laboratory (PPPL). Intense neutron irradiation produces significant changes in the physical and mechanical properties of Fe-(8%–9%)Cr-based reduced activation ferritic martensitic (RAFM) steels. An ANSYS model of the breeding blanket was built to involve these irradiation-induced material property changes with neutron fluence distribution mapped from the MCNP model. Simulation results show that significant structure deformation forms from nuclear swelling, stress, and plastic strain arising from DPA gradient through wall thickness. The material degradation effect is not negligible. This modeling feature can either assist in analyzing the structural behaviors with the influence of nuclear swelling or provide guidance to design the structure to withstand irradiation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The Shape Factor for Pits and Its Impact on Pit Stability

This study employs the finite element method (FEM) to predict the impact of pit shape on pit stability via shape factors of various pit geometries relevant to localized corrosion. Lower values of the shape factor indicate an increased ease in maintaining pit stability. Analyzed geometries include undercut pits, bispherical pit-within-pit structures, and covered pits with perforated (lacy) covers. The effect of the water layer thickness, transport properties of the electrolyte inside and outside the pit, and cathode location on the pit shape factor were also explored. Results show that occluded pits exhibit lower shape factors than open ones, with disk-shaped pits decreasing further as c/r ratio and occlusion angle increase. The findings also suggest minimal influence from a secondary pit if the primary remains active. An equation is presented that quantifies the impact of lacy covers, revealing significant shape factor value reduction. Additionally, high salt concentrations inside pits have a limited stabilizing effect compared to geometry, while thin water layers and adjacent cathodic/sinks reduce pit stability.

Shehi, A. (ORCID:0009000271548041)↗
DOE OSTI2025

Seismic Contingency Auto Generator

This code takes in premade earthquake scenario XML files from USGS, power grid data, and converts them into a contingency file (.con file) that can be used by power grid solvers. Within the .con file are a number (Specified by the user) of contingencies that have randomly failed power transformers based on their likelihood of failure and peak ground acceleration (PGA) value around the transformer. The transformers' likelihood of failure was calculated based on a variety of finite element modeling on various transformer designed for specific transformer voltage classes. Parameters from these FEM were used to create generic fragility curves for transformers within a specific voltage class, which correspond with earthquake PGA values to produced a probability of failure for a given earthquake scenario. More refined versions of this process, such as specifying specific transformer design categories within a voltage class, could also be applied in future iterations of the software.

Vaagensmith, Bjorn [Idaho National Laboratory (INL↗
DOE OSTI2025

Application of physics-informed neural networks (PINNs) solution to coupled thermal and hydraulic processes in silty sands

Abstract The accurate modeling of water and heat transport in soils is crucial for both geo-environmental and geothermal engineering. Traditional modeling methods are problematic because they require well-defined boundaries and initial conditions. Recently, physics-informed neural networks (PINNs), which incorporate partial differential equations (PDEs) to solve forward and inverse problems, have attracted increasing attention in machine learning research. In this study, we applied PINNs to tackle hydraulic and thermal transport coupling forward problems in silty sands. A fully connected deep neural network was utilized for training. This neural network model leverages automatic differentiation to apply the governing equations as constraints, based on the mathematical approximations established by the neural network itself. We conducted forward problems and compared the solutions derived from PINNs with those from Finite Element Method (FEM) simulations. The forward problem results demonstrate the PINNs model’s capability in predicting hydraulic transport, heat transport, and thermal–hydraulic coupling in silty sands under various boundary conditions. The PINNs exhibited great performance in simulating the thermal–hydraulic coupling problem. The accuracy of the PINNs solutions shows its potential for simulation in geotechnical engineering.

Feng, Yuan↗
DOE OSTI2024

Strain-concentration for fast, compact photonic modulation and non-volatile memory

A critical figure of merit (FoM) for electro-optic (EO) modulators is the transmission change per voltage, d T / d V . Conventional approaches in wave-guided modulators maximize d T / d V via a high EO coefficient or longer light-material interaction lengths but are ultimately limited by material losses and nonlinearities. Optical and RF resonances improve d T / d V at the cost of spectral non-uniformity, especially for high- Q optical cavity resonances. Here, we introduce an EO modulator based on piezo-strain-concentration of a photonic crystal cavity to address both trade-offs: (i) it eliminates the trade-off between d T / d V and waveguide loss—i.e., enhancement of the resonance tuning efficiency d v c / d V for the fixed EO coefficient, waveguide length, and cavity Q —and (ii) at high DC strains it exhibits a non-volatile (NV) cavity tuning Δ v c ,NV for passive memory and programming of multiple devices into resonance despite fabrication variations. The device is fabricated on a scalable silicon nitride-on-aluminum nitride platform. We measure d v c / d V =177±1MHz/V, corresponding to Δ v c =40±0.32GHz for a voltage spanning ±120V with an energy consumption of δ U /Δ v c =0.17nW/GHz. The modulation bandwidth is flat up to ω BW,3dB /2 π =3.2±0.07MHz for broadband DC-AC and 142±17MHz for resonant operation near a 2.8 GHz mechanical resonance. Optical extinction up to 25 dB is obtained via Fano-type interference. Strain-induced beam-buckling modes are programmable under a “read-write” protocol with a continuous, repeatable tuning range of 5±0.25GHz, allowing for storage and retrieval, which we quantify with mutual information of 2.4 bits and a maximum non-volatile excursion of 8 GHz. Using a full piezo-optical finite-element-model (FEM) we identify key design principles for optimizing strain-based modulators and chart a path towards achieving performance comparable to lithium niobate-based modulators and the study of high strain physics on-chip.

Wen, Y. Henry (ORCID:0009000685423628)↗
DOE OSTI2024

Modeling of Bi-2212 strand and Rutherford cable in a dipole coil insert

The U.S. Magnet Development Program (US-MDP) explores high-field accelerator magnets compatible with op- erational conditions beyond the limits of Nb3Sn technology. The ongoing R&D High-Temperature Superconductors (HTS) suggests using Bi2Sr2CaCu2O8−x (Bi-2212) as superconducting element. Bi-2212 Rutherford cables maintain a high critical current (IC ) when exposed to a large external magnetic field. However, Bi-2122 exhibits an oversensitive stress-strain response when subject to large Lorentz forces. This paper reports on the magnetic and mechanical analysis of the Bi-2212 cosine-theta insert being developed at Fermilab for a hybrid magnet composed of 2 external layers of Nb3Sn and 2 internal layers of Bi-2212. We performed a FEM analysis of the insert to estimate the HTS stress state in the coil s strands under magnetic and mechanical loads.

D'Agliano, A.↗
DOE OSTI2024

Near-net shaping of high-density, complex, spatially-tailored ceramics for extreme environments

The purpose of this subcontract is for LLNL to collaborate with the Subcontractor (UCSB) in the areas of binder jet additive manufacturing (BJAM), hot pressing (HP), finite element modeling (FEM), and multi-material sintering studies. LLNL wants to further develop its hot pressing, AM, and modeling capabilities through LLNL’s Early Career Faculty Program. This work will focus on developing a method to hot press complex-shaped ceramic parts through the co-deposition of a structural/functional ceramic and sacrificial mold material in one assembly via multi-material binder jet printing. The method aims to produce dense, near-net shaped parts with applications relevant to a wide variety of fields including aerospace, defense, and energy.

36 MATERIALS SCIENCE↗
DOE OSTI2024

R&D to Ensure a Scientific Basis for Qualification Tests and Standards (Final Report)

Project return on investment in a photovoltaic (PV) system depends increasingly on maintaining high energy yields, and the system lifetime is a major factor in levelized cost of electricity (LCOE). Thus, the rate of PV deployment and the success of these assets depends upon reliable long-term power generation. The overarching objective of this program is to improve photovoltaic (PV) module reliability via development of tests and standards. Where reliability problems or risk are discovered, we can design tests to ensure that these liabilities don't affect future generations of products. Customers can use these tests to understand which products are susceptible to certain degradation mechanisms, and manufacturers can use the tests to design unwanted characteristics out of their products. The work under this program identifies PV reliability needs, performs characterization that provides scientific understanding of targeted degradation mechanisms, and translates those data into practical and predictive test protocols and standards. Major accomplishments include: A model for polarization-type potential induced degradation (PID-p) was developed and validated against experimental data. NREL is currently leading a new edition of IEC 62804-1 for PID detection. PID-p can cause large losses in current and voltage for some module designs on cloudy days. Finite element modeling (FEM) and experiment was used to determine when cells crack in a module. It was shown that cells in landscape orientation are much more likely to crack than those on portrait orientation. Shortly thereafter, the first products with portrait-oriented cells were introduced. Studies of how to test for light and elevated temperature degradation (LeTID) culminated with the publication of IEC TS 63342. Software to predict the progression of LeTID was developed, validated, and made publicly available. Field validated tests and international standards for durability of PV module coatings abrasion, backsheets, and encapsulants were developed. Examples are IEC 62788-1-1, IEC 62788-2 ED2, IEC TS 62788-7-2, IEC 62788-7-3 ED1, IEC 63209-2. NREL led the development a high-temperature testing technical specification, and published guidelines that enable installers to determine whether higher-temperature testing is needed, simply based on location and mounting configuration. In a number of our case studies, variations in the bills of materials or workmanship have been associated with variations in reliability. These observations emphasize the importance of quality assurance to reliability. A framework for criticality (i.e. Pareto) analysis was developed and published. The framework helps us and other researchers determine what problems should be addressed for reliability research to have the biggest industry impact. NREL continues to participate actively in international standards development and stakeholder engagement activities, including organizing an annual PV Reliability Workshop. These activities are important for ensuring we address issues that are relevant and timely, and that we convey our results to those who may benefit.

14 SOLAR ENERGY↗
DOE OSTI2026

Multiphysics and Multiscale Simulation Methods for Electromagnetic Energy Assisted Fossil Fuel to Hydrogen Conversion (Final Scientific/Technical Report)

This report summarizes the technical accomplishments of the four-year research project “Multiphysics and Multiscale Simulation Methods for Electromagnetic Energy Assisted Fossil Fuel to Hydrogen Conversion” (Award No. DE-FE0032092), conducted at Howard University and the University of Houston (subawardee) from September 2021 to August 2025. The project successfully achieved all four major objectives: 1. 3D Structural Characterization – Developed 3D optical imaging and mechanical sectioning methods to characterize catalyst distribution and support morphology in nickel foam substrates. Successfully reconstructed 3D geometries and imported them into COMSOL Multiphysics for electromagnetic simulations. 2. EM Hotspot Simulation – Created all-frequency stable electromagnetic formulations and 3D nodal discontinuous Galerkin (NDG) methods for coupled electromagnetic-thermal-fluid problems in multiscale catalytic media. Demonstrated stable solutions from DC to microwave frequencies. 3. Multiphysics Coupling – Developed multiscale simulation methods coupling FEM electromagnetic solvers with thermal transport equations. Reactive molecular dynamics (ReaxFF MD) simulations were performed to investigate catalytic reaction mechanisms at the atomistic level. Demonstrated electromagnetic-thermal co-simulation capabilities for porous catalyst structures. 4. System Optimization – Designed and optimized EM-assisted catalytic systems using nickel foam and carbon foam structures, demonstrating significant temperature increases due to microwave heating. Observed and characterized plasma generation in carbon fiber catalysts. Investigated multiple reaction chamber geometries for improved microwave energy deposition. The project produced significant scientific contributions including 15+ peer-reviewed publications, trained multiple Ph.D. students and undergraduate researchers, and advanced the understanding of microwave-assisted hydrogen production from fossil fuels.

08 HYDROGEN↗
DOE OSTI2026

Numerical Investigation of Fluid Flow and Space Charge in Liquid Argon Time Projection Chamber (LArTPC) Detectors

Overview This project focused on developing a high-fidelity numerical framework to simulate the multiphysics environment within Liquid Argon Time Projection Chamber (LArTPC) detectors. The primary objective was to characterize the complex interplay between ion transport, background fluid dynamics, and electric field distortions—a critical factor for the calibration and sensitivity of next-generation High Energy Physics experiments, such as DUNE. Technical Achievements The research successfully yielded a hybrid numerical space-charge solver utilizing a Cell-Centered Finite Volume Method (FVM) for ion transport coupled with a Finite Element Method (FEM) for electric potential. Key accomplishments include: • Verification & Validation: The 3-D solver was rigorously verified against 1-D analytical solutions, demonstrating high numerical accuracy in predicting space-charge-induced field deviations. • Field Distortion Analysis: 3D simulations revealed that space charge effects introduce significant non-uniformities in the electric field. Critically, the research identified that background LAr flow velocities, when comparable to ion drift velocities, markedly exacerbate these distortions. • Technology Transfer: The resulting source code and comprehensive user manuals were successfully transferred to collaborators at Fermilab, providing a portable computational tool for the broader scientific community. Challenges and Future Directions While the space-charge solver achieved all performance metrics, the integrated fluid dynamics modeling encountered convergence challenges stemming from the extreme 200-fold disparity in length scales between the detector's 37 mm inlet pipes and the 8-meter global domain. To address this, the project has identified a clear technical pivot toward Hierarchical Geometric Adaptive Mesh Refinement (HG-AMR). By implementing an h-type refinement strategy with hanging nodes, future iterations of this solver will be capable of resolving localized high-gradient inlet flows without the prohibitive computational costs of regular grids. This advancement, combined with data-driven uncertainty quantification based on MicroBooNE-style calibration, will enable the precise modeling of detector responses in large-scale cryogenic environments where direct measurement remains difficult. Impact The computational tools developed under this award provide a foundation for enhancing the energy resolution and spatial reconstruction of noble liquid detectors. By bridging the gap between theoretical fluid dynamics and experimental field calibration, this work supports the DOE’s mission to advance the frontiers of neutrino physics and dark matter detection.

42 ENGINEERING↗
DOE OSTI2026

Positron Injector for Ce+BAF

A baseline concept for a continuous wave (CW) polarized positron injector was developed for the CEBAF at Jefferson Lab. This concept is based on the generation of CW longitudinally polarized positrons by a high-current, polarized electron beam 1 mA, 130-370 MeV, and 90% longitudinal polarization that passes through a rotating, water-cooled, tungsten target. The simulations for the Ce+BAF injector at the Low Energy Recirculator Facility (LERF) are performed, including positron beam generation, capture, energy selection, and acceleration to 123 MeV. The positron yield (or positron current) and longitudinal polarization are calculated considering the longitudinal and transverse CEBAF acceptances (<1% energy spread, 1.2 mm bunch length and normalized emittance of 100 mm·mrad). The impact of target thickness, electron beam energy, and transverse size on positron yield within the required emittance limit is evaluated. Target thermal and structural FEM analyses (ANSYS Fluent) are performed to determine the maximum electron beam current and minimum transverse beam size on the target.

Benesch, Jay [Thomas Jefferson National Accelerato↗
DOE OSTI2026

Final Technical Report - Center for Simulation of Fusion Relevant RF Actuators

We have developed a suite of 3D electromagnetic field solvers, both FEM and FDTD based, that account for the RF antenna and vacuum vessel geometries with unprecedented accuracy. Workflows were developed that make it possible to translate CAD models for the antenna and vacuum vessel to physics meshes for RF wave simulation. Nonlinear RF sheath formation has been incorporated self-consistently as a boundary condition in these solvers. We have also carried out extensive studies of the impact of RF sheaths on the ion energy angle distribution at plasma-material interfaces, using high fidelity particle-in-cell codes. Comprehensive simulation models were developed to assess the impact of blob-like edge turbulence on RF wave propagation and the impact of the RF ponderomotive force on the plasma scrape-off layer (SOL). A fluid transport solver for the far-SOL was also developed which accounts for the high parallel to perpendicular heat anisotropy on an unstructured mesh, thus making it possible to precisely represent an antenna structure in the presence of edge transport. Finally we have developed a hierarchy of core wave propagation and absorption models that self-consistently combine continuum Fokker Planck and Monte Carlo treatments of fast ion evolution with ICRF full-wave field solvers and continuum Fokker Planck treatments of fast electron evolution with both full-wave and ray tracing models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

SAM Finite Volume Method Development Status Update: GCR Application, Restart, and MultiApp

The System Analysis Module (SAM) is being developed as a modern system analysis code for advanced non-light-water-reactor safety analysis under the U.S. DOE NEAMS program. Previous feasibility studies have demonstrated that a staggered-grid finite volume method (SG-FVM), implemented under the MOOSE framework, can deliver more than an order of magnitude speedup over the existing continuous Galerkin finite element method (CG-FEM) solver for liquid-cooled, incompressible but thermally expandable flow systems. This work extends the previous effort to compressible, gas-cooled reactor applications, where pressure couples directly into the mass equation adding additional nonlinearity into the equation system. New code capabilities are implemented for pebble bed high-temperature gas-cooled reactor (PB-HTGR) analysis, including a pebble bed CoreChannel component, built-in pebble bed effective thermal conductivity model and channel-to-channel crossflow model. The capabilities are tested, benchmarked, and demonstrated for problems with increased level of model and physical complexities, including the HTTU effective thermal conductivity test, the SANA passive cooling test, and a demonstration case using the GPBR200 reactor design covering steady-state operation, DLOFC and PLOFC transients. Across all cases, the SG-FVM solver demonstrated strong robustness and efficiency, and the solutions agree well with reference results and data. The finding of this work proves that SG-FVM is a viable and efficient solver pathway for compressible, gas-cooled reactor system analysis in SAM. In addition, work has been done to successfully support SAM-FVM recover/restart code feature that is essential to reactor safety analysis applications, and MultiApp code feature that is essential to multi-scale and multi-physics simulations. In summary, this work continued from previous feasibility studies, and further demonstrated that the SG-FVM will serve as a strong foundation for SAM’s advanced solver algorithm for future deployment.

Zou, Ling↗