Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element methods”

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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Bayesian optimization for inverse calibration of expensive computer models: A case study for Johnson-Cook model in machining

Inverse model calibration for identifying the constitutive model parameters can be computationally demanding for expensive-to-evaluate simulation models. Here, this paper presents a modified Bayesian optimization (BO) method, denoted as BO-bound, that incorporates theoretical bounds on the quantity of interest. A case study for the inverse calibration of the Johnson Cook (J-C) flow stress model parameters is presented using machining (cutting) force data. The results show fast calibration of the five J-C parameters within 25 simulations. In general, the BO-bound method is applicable for inverse calibration of any expensive simulation models as well as optimization problems with known bounds.

Bayesian optimization↗
DOE OSTI2025

Orbit-averaging and deposition accuracy for runaway electron beams in hybrid kinetic-MHD simulations of the runaway plateau

We develop a new procedure that combines the kinetic orbit runaway electron code (KORC) and the NIMROD extended-magnetohydrodynamic code to simulate runaway electrons (REs) in the post-disruption plateau. KORC integrates guiding-center orbits, with a barycentric-based binary search strategy providing initial guesses for the Newton–Raphson logical-to-physical coordinate inversion, ensuring reliable particle-to-mesh mapping in NIMROD, whose fields remain static for the present study. Samples are drawn in accord with experimental parallel current profiles of RE beams during the plateau phase. Deposition in NIMROD is verified through comparison with a Python-based finite-element code that ensures periodicity in the poloidal direction and continuity at the magnetic axis. Accurate representation of near-axis fields requires finer mesh resolution to prevent under- and overshoots in current density from orbit inaccuracies. Yet, at a fixed particle count, increasing mesh resolution amplifies statistical noise in the deposited fields. An orbit-averaging method accumulates partial current deposits over multiple kinetic steps and reduces the statistical noise with little added computational cost. By coupling kinetic routines from KORC directly into the NIMROD codebase, these developments lay essential groundwork for future self-consistent KORC–NIMROD coupling.

Algorithms and data structure↗
DOE OSTI2025

Numerical Simulation Studies of Ultrasonic De-Icing for Heating, Ventilation, Air Conditioning, and Refrigeration Structures

Ice accumulation on heating, ventilation, air conditioning, and refrigeration (HVACR) structures presents significant operational challenges. These challenges include reduced efficiency, increased energy consumption, and potential damage to equipment. Traditional de-icing methods, such as chemical treatments, mechanical scraping, or heating-based techniques, are often labor-intensive, costly, and environmentally harmful. Here, this study uniquely investigates ultrasonic de-icing as an energy-efficient alternative for HVACR applications, focusing on the specific structural geometries found in these systems. A comprehensive numerical simulation framework was developed using finite element analysis to explore ultrasonic wave propagation across four distinct HVACR structures. Key parameters such as ultrasonic frequency, power levels, and the number and placement of actuators were examined for their impact on ice detachment efficiency. Results from simulations on a plate structure reveal that ultrasonic excitation can propagate effectively across large areas (at least 150 × 150 mm), enhancing the de-icing coverage. Lower frequency (e.g., 30 to 45 kHz) excitation results in greater displacement, improving de-icing performance, while increased actuator numbers with the same total power input also enhance effectiveness. Two actuators seem sufficient for the de-icing of a 300 × 300 mm plate. For tube-and-fin structures, specific high-power ultrasonic frequencies selectively excite the fin plates, demonstrating efficient ice removal when actuated on the tube. However, optimal performance requires careful design of actuator placement and vibration modes to accommodate the irregular shapes of these structures.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

MolSym : A Python package for handling symmetry in molecular quantum chemistry

A consideration of the point group symmetry of molecules is often advantageous from a computational efficiency standpoint and sometimes necessary for the correct treatment of chemical physics problems. Many modern electronic structure software packages include a treatment of symmetry, but these are sometimes incomplete or unusable outside of that program’s environment. Therefore, we have developed the MolSym package for handling molecular symmetry and its associated functionalities to provide a platform for including symmetry in the implementation and development of other methods. Features include point group detection, molecule symmetrization, arbitrary generation of symmetry element sets and character tables, and symmetry adapted linear combinations of real spherical harmonic basis functions, Cartesian displacement coordinates, and internal coordinates. We present some of the advantages of using molecular symmetry as achieved by MolSym, particularly with respect to Hartree–Fock theory, and the reduction of finite difference displacements in gradient/Hessian computations. Furthermore, this package is designed to be easily integrated into other software development efforts and may be extended to further symmetry applications.

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

The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas

JOREK is a massively parallel fully implicit non-linear extended magneto-hydrodynamic (MHD) code for realistic tokamak X-point plasmas. It has become a widely used versatile simulation code for studying large-scale plasma instabilities and their control and is continuously developed in an international community with strong involvements in the European fusion research programme and ITER organization. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and physics studies performed with the code. A dedicated section highlights some of the verification work done for the code. A hierarchy of different physics models is available including a free boundary and resistive wall extension and hybrid kinetic-fluid models. The code allows for flux-surface aligned iso-parametric finite element grids in single and double X-point plasmas which can be extended to the true physical walls and uses a robust fully implicit time stepping. Particular focus is laid on plasma edge and scrape-off layer (SOL) physics as well as disruption related phenomena. Among the key results obtained with JOREK regarding plasma edge and SOL, are deep insights into the dynamics of edge localized modes (ELMs), ELM cycles, and ELM control by resonant magnetic perturbations, pellet injection, as well as by vertical magnetic kicks. Also ELM free regimes, detachment physics, the generation and transport of impurities during an ELM, and electrostatic turbulence in the pedestal region are investigated. Regarding disruptions, the focus is on the dynamics of the thermal quench (TQ) and current quench triggered by massive gas injection and shattered pellet injection, runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events. Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring TQs triggered by locked modes, and radiative collapses are being studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Specimen Size Artifacts Associated with a Glovebox Deployable Laser Flash Diffusivity Measurement System

To accelerate nuclear fuel qualification and deployment efforts, Oak Ridge National Laboratory has led the world in developing vehicles for accelerating burnup accumulation in nuclear fuel by testing miniature samples in the High Flux Isotope Reactor. However, a challenge remains to establish parity between post irradiation characterization methods for miniature fuel specimens in comparison to conventional fuel elements. This work aims to identify artifacts associated with the thermal analysis of miniature disk specimens in comparison to samples of standardized dimensions through a combination of experimental and computational tools. Using a newly procured Netzsch 427 laser flash diffusivity system, preliminary data on mild steel specimens were collected, and finite element techniques were used to model heat transfer using 1D conditions to assess issues associated with pulse width and data collection frequency effects. These preliminary investigations found that significant variability in measured thermal diffusivity exist within datasets. These variations are not necessarily attributed to limitations with respect to data collection rate; nor are they associated with limitations on pulse width. Instead, these artifacts may be a result of limitations in detector intensity. Future work aims to design a focusing lens to better measure the smaller surface area of miniature specimens analyzed using the laser flash analysis system.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

An assessment of the utility of multirate time integration for the modeling of laser powder bed fusion

Finite element simulation of the laser powder bed fusion process is made challenging by the inherently multiscale nature of the process. When using the typical global time stepping techniques, slowly-evolving regions of the domain receive the same numerical treatment as the regions with the highest temperature rates. The current work details the implementation and evaluation of an implicit multirate method which is able to advance different regions of the domain with distinct time steps, depending on their current solution rate. Previous work indicates that this representation of the temporal scales of the problem can lead to significant reductions in wall clock run time, and it is shown herein that multirate time integration, when used with a uniform process-scale mesh, can result in speedups of approximately 19, 42, and 87 times, for domains with edge lengths of approximately 1 mm, 2 mm, and 5 mm, respectively. When used in conjunction with $\ h$-refinement (limited to two levels of refinement), the resulting speedups (taken relative to the uniform mesh) are around 28, 85, and 665 times. The method is demonstrated to converge as indicated by the literature, and its use with an AM-Bench domain is demonstrated.

42 ENGINEERING↗
DOE OSTI2021

A 3D Full Stress Tensor Model for Oklahoma

Abstract The stress tensor is an important property for upper crustal studies such as those that involve pore fluids and earthquake hazards. At tectonic plate scale, plate boundary forces and mantle convection are the primary drivers of the stress field. In many local settings (10–100 s of km and <10 km depth) in tectonic plate interiors, we can simplify by assuming a constant background stress field that is perturbed by local heterogeneity in density and elasticity. Local stress orientation and sometimes magnitude can be estimated from earthquake and borehole‐based observations when available. Modeling of the local stress field often involves interpolating sparse observations. We present a new method to estimate the 3D stress field in the upper crust and demonstrate it for Oklahoma. We created a 3D material model by inverting multiple types of geophysical observations simultaneously. Integrating surface‐wave dispersion, local travel times and gravity observations produces a model of P‐wave velocity, S‐wave velocity, and density. The stress field can then be modeled using finite element simulations. The simulations are performed using our simplified view of the local stress field as the sum of a constant background stress field that is perturbed by local density and elasticity heterogeneity and gravitational body forces. An orientation of N82°E, for the maximum compressive tectonic force, best agrees with previously observed stress orientations and faulting types in Oklahoma. The gravitational contribution of the horizontal stress field has a magnitude comparable to the tectonic contribution for the upper 5 km of the subsurface.

58 GEOSCIENCES↗
DOE OSTI2026

Determining Stability Margins in Adiabatic Superconducting Magnets with 3-D Finite Element Analysis

Superconducting magnets play a key role in the development of experiments at Fermilab; understanding the operating stability of these can allow us to utilize more potent magnets for future experiments (like the proposed Muon Collider), optimize the design of magnets in more immediate experiments (like Mu2e), and research the future use of more exotic materials (like high-temperature superconductors). This summer, I developed a 3-D parametric FEA program in ANSYS Mechanical APDL that simulates quench in superconducting magnets, and I also developed a parametric MATLAB program that predicts thermal behavior in magnet quench using the MIITS method. These programs can provide useful quenching parameters (like minimum quench energy and normal zone propagation velocity) for different cases of quench, leading to the previously mentioned objective of magnet design optimization. To test the programs, preliminary cases were run and the data produced was compared and analyzed. The results of these analyses, as well as the program operating methods, are discussed in this project.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

The black hole weak gravity conjecture with multiple charges

We study the effect of higher-derivative corrections on asymptotically flat, fourdimensional, dyonic black holes in low-energy models of gravity coupled to N U(1) gauge fields. For large extremal black holes, the leading O (1/Q 2 ) correction to the extremality bound is calculated from the most general low-energy effective action containing operators with up to four derivatives. Motivated by the multi-charge generalization of the Weak Gravity Conjecture, we analyze the necessary kinematic conditions for an asymptotically large extremal black hole to decay into a multi-particle state of extremal black holes. In the large black hole regime, we show that the convex hull condition degenerates to the requirement that a certain quartic form constructed from the Wilson coefficients of the fourderivative effective operators, is everywhere positive. Using on-shell unitarity methods, we show that higher-derivative operators are renormalized at one-loop only if they generate local, on-shell matrix elements that are invariant tensors of the electromagnetic duality group U(N). The one-loop logarithmic running of the four-derivative Wilson coefficients is calculated and shown to imply the positivity of the extremality form at some finite value of Q 2 . This result generalizes an argument recently given by Charles [1], and shows that under the given assumptions the multi-charge Weak Gravity Conjecture is not a Swampland criterion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Adaptive finite differencing in high accuracy electronic structure calculations

Abstract A multi-order Adaptive Finite Differencing (AFD) method is developed for the kinetic energy operator in real-space, grid-based electronic structure codes. It uses atomic pseudo orbitals produced by the corresponding pseudopotential codes to optimize the standard finite difference (SFD) operators for improved precision. Results are presented for a variety of test systems and Bravais lattice types, including the well-known Δ test for 71 elements in the periodic table, the Mott insulator NiO, and borax decahydrate, which contains covalent, ionic, and hydrogen bonds. The tests show that an 8th-order AFD operator leads to the same average Δ value as that achieved by plane-wave codes and is typically far more accurate and has a much lower computational cost than a 12th-order SFD operator. The scalability of real-space electronic calculations is demonstrated for a 2016-atom NiO cell, for which the computational time decreases nearly linearly when scaled from 18 to 144 CPU-GPU nodes.

Briggs, E. L. (ORCID:0000000343983492)↗
DOE OSTI2024

A novel digital lifecycle for Material‐Process‐Microstructure‐Performance relationships of thermoplastic olefins foams manufactured via supercritical fluid assisted foam injection molding

Abstract This research significantly enhances the applicability of thermoplastic olefins (TPOs) in the automotive industry using supercritical N 2 as a physical foaming agent, effectively addressing the limitations of traditional chemical agents. It merges experimental results with simulations to establish detailed material‐process‐microstructure‐performance (MP2) relationships, targeting 5–20% weight reductions. This innovative approach labeled digital lifecycle (DLC) helps accurately predict tensile, flexural, and impact properties based on the foam microstructure, along with experimentally demonstrating improved paintability. The study combines process simulations with finite element models to develop a comprehensive digital model for accurately predicting mechanical properties. Our findings demonstrate a strong correlation between simulated and experimental data, with about a 5% error across various weight reduction targets, marking significant improvements over existing analytical models. This research highlights the efficacy of physical foaming agents in TPO enhancement and emphasizes the importance of integrating experimental and simulation methods to capture the underlying foaming mechanism to establish material‐process‐microstructure‐performance (MP2) relationships. Highlights Establishes a material‐process‐microstructure‐performance (MP2) for TPO foams Sustainably produces TPO foams using supercritical (ScF) N 2 with 20% lightweighting Shows enhanced paintability for TPO foam improved surface aesthetics Digital lifecycle (DLC) that predicts both foam microstructure and properties DLC maps process effects & microstructure onto FEA mesh for precise prediction

Engineering↗
DOE OSTI2022

Vibrational levels of a generalized Morse potential

A Generalized Morse Potential (GMP) is an extension of the Morse Potential (MP) with an additional exponential term and an additional parameter that compensate for MP’s erroneous behavior in the long range part of the interaction potential. Because of the additional term and parameter, the vibrational levels of the GMP cannot be solved analytically, unlike the case for the MP. We present several numerical approaches for solving the vibrational problem of the GMP based on Galerkin methods, namely, the Laguerre Polynomial Method (LPM), the Symmetrized LPM, and the Polynomial Expansion Method (PEM), and apply them to the vibrational levels of the homonuclear diatomic molecules B 2 , O 2 , and F 2 , for which high level theoretical near full configuration interaction (CI) electronic ground state potential energy surfaces and experimentally measured vibrational levels have been reported. Overall, the LPM produces vibrational states for the GMP that are converged to within spectroscopic accuracy of 0.01 cm –1 in between 1 and 2 orders of magnitude faster and with much fewer basis functions/grid points than the Colbert–Miller Discrete Variable Representation (CN-DVR) method for the three homonuclear diatomic molecules examined in this study.

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

Multiscale investigation of thermomechanical and compositional developments in Ni alloy 718 under laser processing

Laser processing has been widely employed in various applications due to its exceptional spatial resolution. However, the rapid temperature gradients generated in localized areas present significant challenges for experimental characterization using conventional instruments. To characterize Ni alloy 718 during laser processing, we employed in-situ synchrotron X-ray diffraction with a high-speed detector, a method particularly well-suited for probing processes with high temporal and spatial resolution. Through a series of in-situ experiments, we investigated the local variations in the evolution of microstructures and thermomechanical behaviors within a keyhole mode melt pool. The in situ macroscopic thermomechanical behaviors were quantified using an empirical model derived from diffraction patterns, with experimental results showing reasonable agreement with finite element analysis. Various laser parameters were tested to assess their influences on the residual strains in the melt pools. The results revealed that the residual strain in the keyhole mode melt pool is relatively insensitive to variations in the parameters and is smaller than that in the melt pool created under conduction mode laser scanning. Additionally, we analyzed the shapes of individual diffraction spots, providing insights into the plastic behaviors and compositional developments in the resolidified alloy. The analysis confirmed that compositional variations in a dendritic microstructure manifest as asymmetric broadening of the diffraction spots.

Ni alloy 718↗
DOE OSTI2020

High precision magnetic field measurement and mapping of the LEReC 180 degree bending magnet using very low field NMR with hall combined probe (140 - 350 gauss)

The Relativistic Heavy Ion Collider (RHIC) at BNL are using the Low Energy RHIC Electron Cooling (LEReC) to conduct experiments that search for the quantum chromodynamic (QCD) critical point. The first ever electron cooling based on the RF acceleration of electron beams was experimentally demonstrated on April 5, 2019 at LEReC at BNL. The first critical step in obtaining successful 3D non-magnetized cooling of the Au ion bunches in the RHIC cooling section was matching the electron beam energy with a relative error less than 5*10 -4 to the ion beam energy. Part of the LEReC beamline is a dipole magnet that bends the electron beam 180 degree. One of the most outstanding measurement challenges is that the dipole field is so low (≈200 G). Most of the existing NMR probes can only measure fields >400 G. Lower signal-to-noise ratio at low fields is requires the use of larger sample volumes. Working with CAYLAR, a NMR probe has been redesigned and optimized for these low field measurements with high resolution. We report the methods, challenges, and results for extensive magnetic field mappings of the 180 dipole magnet. A combination of NMR and Hall sensors has been successfully implemented to measure uniform field regimes inside the magnet center area and non-uniform field regimes at the magnet ends. Detailed measurement and mapping have been performed at five radii and five heights along the beam trajectory. Meanwhile a finite element magnetic modeling simulation of the magnet using Opera software has been performed. The calculated and measured data are compared, and the calculated data is good reference for measured data over long length mapping from magnet edge to center. The measured magnetic measurement data is directly useful for beam instrumentation, diagnostics and operation.

Magnetic Field Measurement↗
DOE OSTI2022

EEL: A MOOSE-based Application for 3D Electro-Chemo-Thermo-Mechanical Modeling of Solid-State Batteries

This report describes a tool for three dimensional, high fidelity, coupled electro-chemo-thermo-mechanical modeling of solid-state batteries. A complete version of the tool is available as open source software at https://github.com/hugary1995/eel.git. The theoretical framework of the tool revolves around an inf-sup statement on a total potential, comprising the Helmholtz free energy, the electrical kinetic potential, the chemical potential, the Fourier potential, the chemical reaction potential, and the external power expenditure. The tool uses the finite element framework of Multiphysics Object-Oriented Simulation Environment (MOOSE) and a variational formulation to solve the boundary value problem for solid-state battery incorporating the full set of multiphysics couplings. The variational formulation also enables a modular software architecture for the tool so additional physics can be easily included by specifying the new contribution to the total potential. The report discusses several method of manufactured solutions that were used to verify the implementation of different physics and the Butler-Volmer reaction kinetics at the electrolyte-electrode interface. Finally, the report discusses results from complete charge/discharge simulations of a quasi-1D and a 3D solid-state battery.

25 ENERGY STORAGE↗
DOE OSTI2021

Composites and additive manufacturing for high-field coil supports for stellarators

The utilization of composites layered on an additive manufactured substrate, for the production of coil supports for modular coils in small or middle size experimental stellarators is assessed. The focus of the study is a monolithic coil support comprising the coils of a half-period of a stellarator, somewhat similar to the ones in UST_2 and ARIES-CS stellarators. However, the concept may be applicable to quasi-monolithic coil supports (coil forms of NCSX type) or individual coil casings (W7-X type).Coil supports for stellarators require high precision, stiffness and strength for large contorted parts. Traditionally, monolithic coil supports are produced by steel casting/forging and final machining. This production method and material gives accurate, stiff and strong coil supports, but the method may be expensive for monolithic supports due to the geometrical complexity and required accuracy of the structure. In relation to those matters, this work investigates whether a monolithic coil support comprising an additive manufactured resin substrate, which is externally (outward from the coils) surrounded with a thick layer of fibre-reinforced resin, may achieve enough strength and stiffness under middle/high magnetic fields. Finite element calculations are produced to obtain the direction of the principal stresses and their values in compression and tension at different areas of the monolithic support, which is relevant for anisotropic materials. The feasibility of directional application of (carbon) fibres on the winding surface of the stellarator outward from the coils is experimentally tested on a scaled-down additively manufactured prototype of a monolithic support. The strength and stiffness of the composite structure appears sufficient for common magnetic fields in experimental stellarators, and the 3D-composite design and manufacturing was technically feasible.

Queral, Vicente↗