Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite element analysis”

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

Real-space simulations of the hyperbolic plasmon polaritons in 1T′ tungsten ditelluride

Recent discoveries of hyperbolic surface plasmon polaritons (SPPs) in 1T′ WTe 2 have garnered significant attention in 2D materials and nanophotonics. In this study, we employ finite-element simulations to investigate the real-space characteristics of hyperbolic SPPs in thin WTe 2 flakes. Our results show that the SPPs exhibit a pronounced sensitivity to excitation energy and sample thickness. By analyzing the plasmonic field patterns, we extract key plasmonic parameters including plasmon wavelengths, hyperbolic angles, and plasmonic figures of merit. In addition, we examine SPP modes in stacked WTe 2 flakes with varying twist angles, demonstrating that the plasmonic field patterns can be effectively tuned by adjusting the twist angle. Notably, as the twist angle increases, the SPPs undergo a topological transition from open hyperbolic modes to closed elliptic modes. This twist-engineering capability offers promising potential for the development of tunable plasmonic devices based on WTe 2 .

2D materials↗
DOE OSTI2025

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

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Monolithic optoelectronic circuit design for on-chip terahertz applications

We demonstrate a monolithic coplanar stripline platform for on-chip terahertz (THz) generation, transmission, and detection, addressing key challenges of mode purity, bandwidth, and referencing. Capacitive coupling of the photoconductive generator switch enforces pure odd-mode propagation, increases THz field strength, and extends the operational frequency range, achieving 0.05–1.4 THz. Our architecture enables fully monolithic fabrication with amorphous silicon switches, in situ field referencing, and galvanic isolation between generation and detection. Finite-element simulations and experiments confirm that suppressing parasitic modes improves signal integrity, providing a robust platform for high-fidelity THz spectroscopy, ultrafast electronics, and nanoscale quantum materials research.

Capacitive coupling↗
DOE OSTI2025

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

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

Electromagnetic optics↗
DOE OSTI2026

AnisONet: A deep neural operator-based anisotropic permeability upscaler from pore to Darcy scale

Directional permeability variations, which govern directional fluid flow in porous media with anisotropy, are important to accurately predict flow behavior, reactive transport, and fluid–solid interactions for various processes such as enhanced geothermal systems, energy storage devices, and biological systems. However, the intricate architecture of porous media makes it difficult to predict directional permeabilities. In this work, we present a novel machine learning (ML) framework, AnisONet, built upon an integration of a convolutional neural network, Swin transformer, and the deep operator network architecture, designed to predict anisotropic permeability and upscale predictions to larger spatial domains. First, AnisONet was evaluated with three classes of two-dimensional (2D) porous media, including synthetic circular and elliptical grains and natural sandstone grains from micro-computed tomography images. A lattice Boltzmann model (LBM) was used to calculate directional permeabilities at every 10° angle, producing 19 data points per image of porous media. AnisONet is then trained to predict permeability as a function of rotation angle. AnisONet showed strong predictive capability of directional permeability. Second, we tested our model for five upscaling cases with a large image size in the finite-element method (FEM) for 2D Darcy flow with various permeability tensor construction methods. Overall, upscaled permeability tensors in FEM simulations produce a reasonably good match with LBM results, highlighting the importance of selecting appropriate tensor formation strategies for accurate permeability upscaling. AnisONet, as a directional permeability estimator, could be further developed for more complex geometries, with the potential to develop a foundational ML model for various applications in porous media.

42 ENGINEERING↗
DOE OSTI2025

Effect of toolpath in large-format additive manufacturing with bio-derived composites

There has been growing interest in integrating bio-derived composites into Large-format additive manufacturing (LFAM) feedstocks to reduce the use of petroleum-derived materials and reduce the overall carbon footprint of LFAM. However, these materials present unique challenges during manufacturing due to their variability, which can lead to unintended deformations and failures attributable to suboptimal process conditions. While numerical modelling has been extensively employed to simulate numerous manufacturing processes, its application in LFAM with bio-derived composites remains limited. This study addresses this gap by systematically developing a numerical model to simulate the LFAM process using bio-based materials, specifically wood fiber-reinforced polylactic acid (PLA/WF). Experimental investigations were conducted to characterise the thermal and mechanical properties of additively manufactured PLA/WF specimens. Numerical simulations were performed to predict temperature profiles and deformations during LFAM. The effect of varying infill patterns, internal structures, and tool paths on the temperature distribution and deformation of printed parts was explored using the developed model. This article aims to advance the utilisation of bio-derived composites in LFAM systems and provide a comprehensive understanding of the LFAM process. The findings offer valuable insights for optimising process parameters and enhancing the performance of LFAM with bio-based composites.

Large-format additive manufacturing↗
DOE OSTI2025

Flexible PCB Windings Size Optimization for Winding AC Resistance Minimization Under the Sinusoidal Voltage Regime

In this paper, optimization of the winding AC resistance of the flexible printed circuit board (FPCB) wire windings is performed. A one-dimensional model of the FPCB winding is introduced and an equation for the FPCB winding AC resistance and its low and medium frequency approximation are derived. The approximate FPCB winding AC resistance equation is used to derive optimum thickness of the FPCB winding conductor thickness at which minimum of winding AC resistance (global optimum) is achieved. The Finite Element Method analysis and experimental verification of derived equations (winding resistance, impedance, power loss, and temperature measurements) is performed in order validate derived model and winding resistance equation.

42 ENGINEERING↗
DOE OSTI2025

Electromagnetic Modeling of Lattice Structures in Additively Manufactured Electric Machines

Recent advancements in metal additive manufacturing (AM) show great potential to revolutionize the design and manufacturing of electromagnetic components used in the field of electrical engineering. Lattice structures directly printed by AM processes typically offer better structural performance with reduced weight, such as high stiffness, surface area, elongation, energy absorption, and porosity, than the solid counterpart. This paper aims to study the electromagnetic modeling of lattice structures used in additively manufactured magnetic cores or windings for electric machines. Three dimensional (3D) electromagnetic finite element (FE) analysis with high performance computing (HPC) shows the highest fidelity in predicting the electromagnetic performance of designs with lattice structures by preserving complex geometry details. FE-based homogenization methods have also been explored to potentially speed up concept design. A case study based on an additively manufactured axial-flux permanent magnet machine with a Hilbert pattern stator validates the discussed electromagnetic modeling approaches.

Chung, Shi-Uk [ANSYS, Inc.]↗
DOE OSTI2025

Design and Construction of a Prototype B1pF Large Aperture Rutherford Cable Superconducting Magnet for the EIC Interaction Region

Brookhaven National Laboratory (BNL) was chosen to host the international Electron-Ion Collider (EIC), which will collide high energy and highly polarized hadron and electron beams with a center of mass energy up to 140 GeV. The Interaction Region (IR) requires several large aperture, relatively high field superconducting dipole and quadrupole magnets, some of which are very closely spaced. B1pF is a large aperture (300 mm coil ID), medium field (4.2 T), 3 meter long superconducting dipole magnet. Its size is larger than the largest superconducting magnet in the RHIC accelerator (180 mm DX dipole, 4.2 T, 3.7 m long). As such, B1pF was chosen as the representative magnet to be prototyped to demonstrate the technological choice, design and construction details. Here, this paper describes the analyses and considerations which informed both the magnet and magnet tooling designs. This paper will summarize the construction results so far, in particular the plan for novel coil winding tooling and the test coil winding & curing program that was carried out prior to coil and magnet construction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

NSTXU Diagnostic Disruption Dynamic Loading Represented by Response Spectra

This article presents the results of transient dynamic simulations of loads due to disruption eddy currents on the NSTXU vacuum vessel. Dynamic loading at diagnostic mounting locations is expressed as response spectra derived from the time history results of the dynamic structural simulations of a variety of disruption scenarios. The disruption simulations draw on a history of the project assessments of worst case disruptions for specific components. Major efforts to assess disruption loading have included the vacuum vessel which is the major structural support for the machine, as well as the passive plates (PPs), high harmonic fast wave (HHFW) antenna, and centerstack casing. Each one of these efforts included transient electromagnetic simulations producing time-dependent eddy current Lorentz loads (and in some cases halo loads) which then were applied to time-dependent structural dynamic analyses intended to obtain the proper dynamic amplification factors. In some instances, the EM model and structural model were identical allowing direct transfer of EM forces to the structural model. In other cases, the EM and structural model were not identical and the vector potential (VP) transfer method was used. The results files from these analyses were available (or re-run) to post process in ANSYS Classic time history postprocessor. In conclusion, the ANSYS command is used to create response spectra from time history data at desired points on the vessel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A Modified DBC Substrate Improving Thermal Performance for Confined Space Applications

High-power modules use substrates to house the semiconductor device and for electrical insulation. These substrates are constructed with thermally conductive dielectric material sandwiched between two metals to extract heat from semiconductor chips. Thus, the required cooling performance of a power module is linked to the substrate’s thermal performance and can vary based on the substrate technologies. Here, in this study, five substrate technologies were evaluated for space-restricted applications: direct-bonded copper (DBC), an insulated metal substrate (IMS), a thermally annealed pyrolytic graphite (TPG)-based IMS, DBC-based double-sided cooling, and direct-bonded aluminum (DBA) where the heat sink is directly attached without thermal interface materials (TIMs). The finite element (FE) analysis results suggest that the popular DBC substrate has thermal performance better than that of the other substrates for space-constrained applications. To further improve the thermal performance, a modified DBC substrate was proposed where a copper block was added between the semiconductor and the DBC substrate to achieve heat spreading underneath the chip. The modified DBC performance was then compared with the aforementioned substrates and the results showed significant thermal performance improvement. The results were verified with experimental results where the proposed substrate showed 20% more loss handling capability compared to an identical DBC substrate.

42 ENGINEERING↗
DOE OSTI2026

Subsurface mechanical damage of fused silica glass during grinding by various sub-aperture tools with and without ultrasonics

The subsurface mechanical damage (SSD) depth after grinding fused silica glass with a comprehensive set of sub-aperture fixed abrasive grinding tools [cup, wheel, belt, pad, and rotary face mill (with and without ultrasonics)] and process parameters has been statistically measured using the taper wedge technique and evaluated. Consistent with a previously reported grinding model [J. Non-Cryst. Solids 352, 5601–5617 (2006) Crossref , Materials Science and Technology of Optical Fabrication (Wiley & Sons, 2018)], based on the sliding indentation fracture by sliding particles or asperities where the normal load per particle determines the depth of the fracture (and ultimately the overall SSD depth), the dominant factor controlling SSD depth was found to be the abrasive size regardless of the tool type and process conditions. Compared to full aperture grinding methods, the overall SSD depth was higher using the sub-aperture tools, likely due to the higher effective pressure and higher load per particle distribution. Here, in addition to abrasive size, a significant reduction in SSD depth was achieved by: (1) reducing the load distribution on the abrasive particles via increase in contact area and/or decrease in mechanical loading; (2) using a more compliant host tool medium; and (3) in what we believe is a more novel way, using ultrasonics. Combining low abrasive size, larger contact area, and a compliant host, the 6 µm diamond in a resin matrix (Trizact) on a foam pad led to very low SSD depth (~ 4.6 µm), relatively fast grinding rate (186 mm 3 /h), and little or no figure degradation. This grinding tool/process is an attractive choice for final grind, resulting in significantly reduced polish out (i.e., “grey out”) times. With the rotary face mill tool, the use of ultrasonics consistently led to a SSD depth reduction (ranging from 17%–34%). A new fracture mechanics-based model, to the best of our knowledge, where the relevant normal load is parallel to the feed direction, has been developed to explain how ultrasonics leads to lower SSD depth. The key factors, supported by finite element stress analysis and load measurements, are (1) the initiation of fractures at higher z heights during the tool’s ultrasonic vertical oscillations, thus propagating less deep into workpiece; (2) reduction in load (and therefore reduction in fracture propagation distance) due to smaller tool-workpiece feed direction contact area (again caused by higher heights relative to depth of cut); (3) upward movement of the tool during oscillation leads to fracturing toward the surface instead into the depth; and finally (4) at tool’s lowest point of oscillation cycle, there may not be enough time for the fracture to propagate to its full length.

Optics and optical instruments↗
DOE OSTI2025

Mesoscale Linear Elastic Modeling and Homogenization of Marine Energy Composites

The design of fiber-reinforced composite (FRC)-based components for marine energy applications necessitates a fundamental understanding of material properties and the resulting geometry to predict long-term performance. In this work, we present a modeling workflow to predict linear elastic and diffusive bulk properties at the mesoscale for an idealized geometry based on knowledge of fiber and resin properties. A parametric study was performed to identify the key model input parameters that influence bulk properties. Furthermore, we demonstrate how bulk properties can be leveraged in high-fidelity image-based simulations, where imperfections in tow geometry and voids captured during X-ray computed tomography imaging are explicitly represented within the simulation. Bulk properties of interest include moduli, Poisson’s ratios, hygroscopic swelling, diffusivity, and moisture uptake, which are key parameters for characterizing FRC performance within marine environments. Modeling predictions agreed well with experimental data, except for estimating swelling coefficients, likely due to crack accumulation as a function of moisture uptake. The mesoscale modeling workflow ultimately highlights a versatile framework for understanding the influence of material and geometric properties, which can be leveraged to rapidly assess new FRC-based components.

computational mechanics↗