Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite Element Simulation”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Low spontaneous Brillouin scattering in anti-resonant hollow-core fibers in GHz frequency range

Brillouin light scattering (BLS) is a powerful experimental tool that can be used to gain insights into the fundamental and applied properties of matter, like dispersions of quasiparticles in a solid, as well as their spatiotemporal dynamics. Many applications of light scattering favor the use of optical fibers in place of free-space optics. In this study, we compare the performance of anti-resonant hollow core fibers to that of conventional solid core fused silica fibers for BLS experiments in the GHz frequency range. Conventional fibers are barely suitable for low-noise measurements because of the spontaneous scattering of photons on various phononic modes present in the core and cladding. In the case of the hollow-core fiber, we identify a range of discrete phononic modes and associate them with the various acoustic modes of the structure surrounding the hollow core using finite-element numerical simulations. The measured relative intensity of the spontaneous BLS signal from these modes is orders of magnitude smaller than that of a solid-core fiber, making anti-resonant hollow-core fibers one of the best solutions for the single-mode light guidance for BLS and potentially other low-noise photonic experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Design and Analysis of an Integrated Additively Manufactured Test Article for Plasma-Facing Components

Current plasma-facing components (PFCs) used in helium-cooled divertor modules are complex structures with tungsten tile, steel sleeve components, and cartridges, all assembled in a helium-cooled multiple jet (HEMJ) structure. The goal of this project is to simplify the complex PFC design using additive manufacturing techniques to create a single integrated tungsten test article. Apart from the flexibility this opens up in exploring a wide array of geometries for the article, having a single integrated article significantly reduces the number of joints and parts in the article, thus reducing chances of leaks. A process called electron beam melting has shown to produce very high-density samples and unique geometries, enabling HEMJ or similar designs. To validate and optimize this novel design, the model underwent a series of computational fluid dynamics and finite element analysis simulations to replicate steady-state heat flux in the divertors. The simulations presented in this study consider a steady-state base heat flux of 5 MW/m 2 , with water serving as the coolant. Future research will explore the use of helium as a coolant, simulate edge-localized-mode conditions, and include experimental validation. Since 3D-printed tungsten is anisotropic, the build direction versus build plane of the article are taken into consideration for the test article strength. Because of the high operating temperatures and low ductility of tungsten, thermal creep and brittle fracture are important failure mechanisms to consider. In conclusion, the cap is evaluated with various flow velocities and nozzle diameters, and an optimal design choice is made for which this cap will survive the divertor conditions with a conservative safety margin.

Divertor↗
DOE OSTI2020

Shield Design for 50 kW Three-Phase Wireless Charging System

In this paper, magnetic field emissions (MFE) from a 50 kW bipolar coil-based three-phase wireless charging system (WCS) is investigated. Then, a shield design is proposed to suppress the MFE below the safety limit identified by the International Commission of Non-Ionized Radiation Protection (ICNIRP) guidelines. Bipolar coil-based three-phase (3φ) WCSs provide extremely high power-density; therefore, the 3φ WCSs are highly promising for high-power fast-wireless-charging of the electric vehicles. For the high-power EV application, it is a critical challenge to design a shield to suppress the MFE below the safety limit. The traditional aluminum shield is found ineffective for the bipolar 3φ-WCS. The proposed shield is designed with high-permeability magnetic material, i.e. ferrite, nanocrystalline material, etc. The proposed shield is simulated in finite element analysis (FEA) and the FEA results show that the proposed magnetic shield effectively suppresses the magnetic field emission below the ICNIRP limit. The generalized design of the proposed shield can be used for even higher power bipolar coil based WCSs.

Mohammad, Mostak↗
DOE OSTI2023

Additively Manufactured Copper Windings with Hilbert Structure

In this paper, finite element method simulations and frequency analysis of additively manufactured windings with a Hilbert–cross section pattern were performed over a wide frequency range. The current density distributions for the proposed winding conductor and for a square-wire winding conductor at selected frequencies were compared. An approximately 25% winding loss reduction at higher frequencies and an approximately 5% increase at low frequencies were observed for the analyzed Hilbert structure.

Wojda, Rafal↗
DOE OSTI2023

Design and Optimization of High-frequency Transformer for Isolated Single-Stage Three-Phase AC/DC Converter using Bidirectional Switches

This paper presents the design and optimization of the high-frequency transformer for an isolated single-stage threephase AC/DC converter enabled by the 4H-SiC BiDirectional Field Effect Transistor (BiDFET). Non-dominated Sorting Genetic Algorithm II (NSGA-II) is used to perform a multiobjective optimization to minimize loss, weight, and volume. A 20 kW, 50 kHz hardware prototype based on the selected optimal point is characterized by a power amplifier. The finite element analysis simulation results are also presented to validate the core magnetic flux density and winding electric field strength. The experimental prototype was tested at 800 Vdc/480 Vrms at 10 kW rated power.

BiDFET, High-frequency transformer, Multiobjective↗
DOE OSTI2020

Comparison of Magnetic Field Emission from Unipolar and Bipolar Coil-Based Wireless Charging Systems

In this study, the leakage magnetic field (LMF) from the unipolar (circular, rectangular, etc.) and bipolar (double-D) coil-based wireless charging systems (WCSs) are investigated and compared. The unipolar and bipolar pads, which are the two most widely adopted charging-pads, have significantly different magnetic field emission characteristics. In this paper, the spatial-distribution and vector-pattern of the leakage magnetic field are studied and compared for similar 11 kW unipolar and bipolar coil-based WCSs. The field emissions are simulated using finite element analysis (FEA) and tested experimentally with an 11 kW bipolar coil-based WCS prototype. The results show that the leakage magnetic field of the unipolar and bipolar pads have widely different vector-orientation and spatial distributions of the leakage magnetic field. The results show that the unipolar coil based WCS has a dominant vertical and the bipolar pads coil-based WCS has dominant horizontally aligned leakage magnetic field. Considering the shielding, the aluminum plate is found highly effective for the unipolar pad and almost ineffective for bipolar pads.

Mohammad, Mostak↗
DOE OSTI2020

Computationally Efficient AC Resistance Model for Stator Winding With Rectangular Conductors

The resistance of the electric machine’s windings is an important parameter when analyzing the machine’s performance, as conduction losses in the winding are often the largest loss component. However, this resistance changes with frequency due to the skin and proximity effects, making it challenging to accurately estimate these losses. In this article, a computationally efficient, general method of calculating the ac resistance of an electric machine winding with rectangular conductors is presented. The method uses a 1-D magnetoquasi-static model of the stator slot region to develop boundary transfer relations between individual turns. These relations are then used to create an impedance matrix from which the ac resistance can be calculated. This technique is more flexible than previous solutions since it can accommodate different winding configurations such as fractional-pitch windings. It is shown that a fractional-pitch winding has a lower ac resistance than a full-pitch winding. Finally, finite-element analysis simulations and experimental results are used to validate the proposed method.

42 ENGINEERING↗
DOE OSTI2022

Laminated Permanent Magnets Enable Compact Magnetic Components in Current-Source Converters

Magnetic components have become one of the primary barriers to high power density converters, especially to current-source converters (CSCs). In CSCs, the inductors/transformers have a predominantly dc flux to store energy, which can be offset by standard permanent magnets (PMs). However, eddy currents in standard PMs induce significant losses and thermal stress for medium-frequency applications. This article proposes using laminated PMs to offset the dc flux while reducing eddy currents, leading to significant reductions in the size, cost, and losses of inductors/transformers. Furthermore, the laminated PMs’ optimal location, orientation, and distribution are investigated to generalize this approach for maximum benefits. Here, three-dimensional finite-element analysis simulation and hardware experiments are presented to validate the effectiveness of the proposed approach in a medium-frequency transformer (MFT) for CSCs. Compared with standard PMs, the proposed use of laminated PMs reduces aggregate core-plus-PM losses by 85% in experiments, and thus, relaxes the MFT thermal design. Finally, the proposed approach is experimentally validated in a 40 kVA flyback-type MFT for a soft-switching solid-state transformer. Compared to the traditional design without any PMs, the proposed design increases the saturation current by 46% while inducing only 5% more losses, leading to significant savings in magnetics cost and size.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2020

Mechanism-Based Energy Regularization in Computational Modeling of Quasibrittle Fracture

Abstract Quasibrittle materials are featured by a strain-softening constitutive behavior under many loading scenarios, which could eventually lead to localization instability. It has long been known that strain localization would result in spurious mesh sensitivity in finite element (FE) simulations. Previous studies have shown that, for the case of fully localized damage, the mesh sensitivity can be mitigated through energy regularization of the material constitutive law. However, depending on the loading configuration and structural geometry, quasibrittle structures could exhibit a complex damage process, which involves both localized and diffused damage patterns at different stages of loading. This study presents a generalized energy regularization method that considers the spatial and temporal evolution of damage pattern. The method introduces a localization parameter, which describes the local damage pattern. The localization parameter governs the energy regularization of the constitutive model, which captures the transition from diffused to localized damage during the failure process. The method is cast into an isotropic damage model, and is further extended to rate-dependent behavior. The energy regularization scheme is directly incorporated into the kinetics of damage growth. The model is applied to simulate static and dynamic failures of ceramic specimens. It is shown that the present model is able to effectively mitigate the spurious mesh sensitivity in FE simulations of both types of failure. The present analysis demonstrates the essential role of mechanism-based energy regularization of constitutive relation in FE simulations of quasibrittle fracture.

Mechanics↗
DOE OSTI2021

Design of a Single Layer Metamaterial for Pressure Offloading of Transtibial Amputees

While using a prosthesis, transtibial amputees may experience pain and discomfort brought on by large pressure gradients at the interface between the residual limb and prosthetic socket. Current prosthetic interface solutions attempt to alleviate these pressure gradients by using soft homogenous liners to reduce and distribute pressures. This research investigates an additively manufactured metamaterial inlay with a tailored mechanical response to reduce peak pressure gradients around the limb. The inlay uses a hyperelastic behaving metamaterial (US10244818) comprised of triangular pattern unit cells which can be 3D printed with walls of various thicknesses controlled by draft angles. The hyperelastic material properties are modeled using a Yeoh 3rd order model. The 3rd order coefficients can be adjusted and optimized, which corresponds to a change in the unit cell wall thickness to create an inlay that can meet the unique offloading needs of an amputee. Finite element analysis simulations evaluated the pressure gradient reduction from: 1) a common homogenous silicone liner, 2) a prosthetist's inlay prescription that utilizes three variations of the metamaterial, and 3) a metamaterial solution with optimized Yeoh 3rd order coefficients. When compared to a traditional homogenous silicone liner for two unique limb loading scenarios, the prosthetist prescribed inlay and the optimized material inlay can achieve equal or greater pressure gradient reduction capabilities. These results show the potential feasibility of implementing this metamaterial as a method of personalized medicine for transtibial amputees by creating a customizable interface solution to meet the unique performance needs of an individual patient.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Electrode Mesoscale as a Collection of Particles: Coupled Electrochemical and Mechanical Analysis of NMC Cathodes

Battery electrodes are composed of polydisperse particles and a porous, composite binder domain. These materials are arranged into a complex mesostructure whose morphology impacts both electrochemical performance and mechanical response. We present image-based, particle-resolved, mesoscale finite element model simulations of coupled electrochemical-mechanical performance on a representative NMC electrode domain. Beyond predicting macroscale quantities such as half-cell voltage and evolving electrical conductivity, studying behaviors on a per-particle and per-surface basis enables performance and material design insights previously unachievable. Voltage losses are primarily attributable to a complex interplay between interfacial charge transfer kinetics, lithium diffusion, and, locally, electrical conductivity. Mesoscale heterogeneities arise from particle polydispersity and lead to material underutilization at high current densities. Particle-particle contacts, however, reduce heterogeneities by enabling lithium diffusion between connected particle groups. While the porous composite binder domain (CBD) may have slower ionic transport and less available area for electrochemical reactions, its high electrical conductivity makes it the preferred reaction site late in electrode discharge. Mesoscale results are favorably compared to both experimental data and macrohomogeneous models. This work enables improvements in materials design by providing a tool for optimization of particle sizes, CBD morphology, and manufacturing conditions.

42 ENGINEERING↗
DOE OSTI2021

Effect of workpiece curvature on the tool influence function during hemispherical sub-aperture tool glass polishing

The influence of workpiece curvature on the tool influence function spot during polishing of fused silica glass with cerium oxide slurry, while using a rotating hemispherical pad-foam tool for a wide variety of process conditions (tool displacement, inclination angle, and rotation rate), has been investigated. (Workpiece curvature ranged from 500 mm radius concave to 43 mm radius convex.) The TIF spot decreases in diameter and increases in the peak removal rate on more convex workpieces. In contrast, the TIF spot increases both in diameter and peak removal rate on more concave workpieces. For the range of workpiece curvatures investigated, both the spot size and the peak removal rate changed significantly, as much as 2 times. An elastic sphere–sphere contact mechanics model, which utilizes both a modified displacement (that leads to a change in the applied load) as well as a mismatch factor (that influences the pressure distribution shape), has been developed. Furthermore, the model was validated using both offline load-displacement measurements and finite-element analysis simulations. The model quantitatively describes the measured change in the relative contact diameter and relative pressure distribution, as well as semiquantitively describes the change in the relative volumetric removal rate on a large variety of TIF spots. The change in the volumetric removal rate for convex workpieces is a result of the balance between a decreasing spot size (reducing removal) and an increasing peak pressure (increasing removal), which usually results in relatively small changes in volumetric removal. In the case of concave workpieces, the volumetric removal rate change is also governed by a similar balance, but the spot size increase contribution dominates, resulting in a significant increase in volumetric removal rate. Understanding these trends can enable methods to add greater determinism during the fabrication of freeform optics by adjusting polishing parameters (such as dwell time) while the tool translates along a workpiece surface with different local curvatures.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Self-assembled multifunctional nanostructures for surface passivation and photon management in silicon photovoltaics

This work reports the fabrication and characterization of multifunctional, nanostructured passivation layers formed using a self-assembly process that provide both surface passivation and improved light trapping in crystalline silicon photovoltaic (PV) cells. Scalable block copolymer self-assembly and vapor phase infiltration processes are used to form arrays of aluminum oxide nanostructures (Al 2 O 3 ) on crystalline silicon without substrate etching. The Al 2 O 3 nanostructures are characterized using scanning electron microscopy (SEM), transmission electron microscopy (TEM), and spectroscopic ellipsometry. Injection-level dependent photoconductance measurements are used to determine the effective carrier lifetime of the samples to confirm the nanostructures successfully passivate the Si surface. Finite element method simulations and reflectance measurement show that the nanostructures increase the internal rear reflectance of the PV cell by suppressing the parasitic optical losses in the metal contact. An optimized morphology of the structures is identified for their potential use in PV cells as multifunctional materials providing surface passivation, photon management, and carrier transport pathways.

14 SOLAR ENERGY↗
DOE OSTI2020

Semi-Empirical Material Model for Hydrogen Uptake Kinetics by 1,4-bis(phenylethynyl)benzene (DEB)-based Getters

We report on development of a semi-empirical model describing the kinetics of H 2 /D 2 reaction with 1,4-bis(phenylethynyl)benzene (DEB)-based getters in pellet and powder form. The physical basis of the model is solid-state diffusion of DEB plus reaction of DEB with H atoms on the surface of nanometer-scale Pd catalyst particles. The purpose of this material model is to facilitate more accurate predictions of performance for DEB-based getter assemblies. These predictions will be accomplished by integrating the model presented here into finite element, reaction/diffusion simulations of getter assembly hydrogen uptake performance. As of Fall 2020, this integration is underway.

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

Initial framework for engineering-scale statistical creep-fatigue modeling

This report describes the integration of new solid and interface-cohesive mechanics systems into MOOSE. The purpose of these new systems is to support the ability of MOOSE to run full-field crystal plasticity finite element method simulations of key material processes in high temperature metallic materials. These simulations could be used to help accurately predict the performance of key high temperature structural materials in future advanced nuclear reactor components. Previous work implemented preliminary versions of many of these systems in MOOSE Apps. The current work reports on their integration into the main MOOSE tensor mechanics module along with associated improvements to the basic formulations and numerical implementations. Finally, the report provides an example of the full-field crystal plasticity simulations now possible in MOOSE, including examples of realistic geometries requiring millions of degrees of freedom to resolve the microstructural features and macroscale geometry.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Design and Scoping Tests on Alloy 617 Using Notched Specimen Geometry to Validate Methods for Multiaxial Stress Relaxation

In FY 2022, a development effort was initiated at the US Department of Energy’s Oak Ridge National Laboratory (ORNL) to examine the multiaxial stress-relaxation behavior and multiaxial stress-state effect on the creep-fatigue (CF) performance for Alloy 617 at elevated temperatures. This effort supported the development of the design rules in the American Society of Mechanical Engineers (ASME) Boiler and Pressure Vessel Code (BPVC), Section III, Division 5. In this work, two types of the notch specimen geometries were designed. An inelastic constitutive model was used to investigate the stress triaxiality and elastic follow-up effects on the notch specimens under CF deformation in the finite element (FE) simulations. The numerical study demonstrated that the stress triaxiality and elastic follow-up caused by the notches both play significant roles in the stress-relaxation behavior. In addition, Alloy 617 CF experiments were designed, and the testing on the notch specimens is ongoing. CF test failure data will be generated on the specimens with notches under various conditions in FY 2023, and the results will be used to validate methods for accounting for the multiaxial stress relaxation effect in the design code.

36 MATERIALS SCIENCE↗
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↗