Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite-Element”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1973

Finite Element Method for Thermal Analysis

A two- and three-dimensional, finite-element thermal-analysis program which handles conduction with internal heat generation, convection, radiation, specified flux, and specified temperature boundary conditions is presented. Elements used in the program are the triangle and tetrahedron for two- and three-dimensional analysis, respectively. The theory used in the program is developed, and several sample problems demonstrating the capability and reliability of the program are presented. A guide to using the program, description of the input cards, and program listing are included.

Heuser, J.↗
NASA NTRS1973

Cracked finite elements proposed for NASTRAN

The recent introduction of special crack-tip singularity elements, usually referred to as cracked elements, has brought the power and flexibility of the finite-element method to bear much more effectively on fracture mechanics problems. This paper recalls the development of two cracked elements and presents the results of some applications proving their accuracy and economy. Judging from the available literature on numerical methods in fracture mechanics, it seems clear that the elements described have been used more extensively than any others in practical fracture mechanics applications.

Aberson, J. A.↗
NASA NTRS1974

Dynamic analysis of a system of hinge-connected rigid bodies with nonrigid appendages

Equations of motion are derived for use in simulating a spacecraft or other complex electromechanical system amenable to idealization as a set of hinge-connected rigid bodies of tree topology, with rigid axisymmetric rotors and nonrigid appendages attached to each rigid body in the set. In conjunction with a previously published report on finite-element appendage vibration equations, this report provides a complete minimum-dimension formulation suitable for generic programming for digital computer numerical integration.

Likins, P. W.↗
NASA NTRS1973

An improved numerical process for solution of solid mechanics problems

This paper gives an overview of the development and status of an improved numerical process for the solution of solid mechanics problems. The proposed process uses a mixed formulation with the fundamental unknowns consisting of both stress and displacement parameters. The problem is formulated either by means of first-order partial differential equations or in a variational form by using a Hellinger-Reissner-type mixed variational principle. For presentation purposes, the components of a numerical process are characterized and the criteria for an ideal process are outlined. Commonly used finite-difference and finite-element procedures are examined in the light of these criteria and it is shown that they fall short in a number of ways. The proposed numerical process, on the other hand, satisfies most of the optimality criteria and appears to be particularly suited for use with the forthcoming generation computers.

Noor, A. K.↗
NASA NTRS1973

Dynamic analysis of a system of hinge-connected rigid bodies with nonrigid appendages

Equations of motion are derived for use in simulating a spacecraft or other complex electromechanical system amenable to idealization as a set of hinge-connected rigid bodies of tree topology, with rigid axisymmetric rotors and nonrigid appendages attached to each rigid body in the set. In conjunction with a previously published companion paper on finite-element appendage vibration equations, this paper provides a complete minimum-dimension formulation suitable for generic programming for digital computer numerical integration.

Likins, P. W.↗
NASA NTRS1973

Recent advances in numerical analysis of structural eigenvalue problems

A wide range of eigenvalue problems encountered in practical structural engineering analyses is defined, in which the structures are assumed to be discretized by any suitable technique such as the finite-element method. A review of the usual numerical procedures for the solution of such eigenvalue problems is presented and is followed by an extensive account of recently developed eigenproblem solution procedures. Particular emphasis is placed on the new numerical algorithms and associated computer programs based on the Sturm sequence method. Eigenvalue algorithms developed for efficient solution of natural frequency and buckling problems of structures are presented, as well as some eigenvalue procedures formulated in connection with the solution of quadratic matrix equations associated with free vibration analysis of structures. A new algorithm is described for natural frequency analysis of damped structural systems.

Gupta, K. K.↗
NASA NTRS1973

Optimization of structures to satisfy a flutter velocity constraint by use of quadratic equation fitting

Using the first and the second derivative of flutter velocity with respect to the parameters, the velocity hypersurface is made quadratic. This greatly simplifies the numerical procedure developed for determining the values of the design parameters such that a specified flutter velocity constraint is satisfied and the total structural mass is near a relative minimum. A search procedure is presented utilizing two gradient search methods and a gradient projection method. The procedure is applied to the design of a box beam, using finite-element representation. The results indicate that the procedure developed yields substantial design improvement satisfying the specified constraint and does converge to near a local optimum.

Motiwalla, S. K.↗
NASA NTRSDate not supplied

Battery Pack Shape Optimization using Transient Heat Conduction Coupled with Cell-Discharge Analysis

Battery electric systems exhibit significant time-dependence, especially when evaluated in the context of an aircraft mission profile with continually changing power demands. Additionally, when evaluating battery-powered aircraft concepts, it is important to accurately compute the temperature of the batteries and properly characterize the thermal response of the system. The temperature of the batteries has a significant impact on cell performance, in addition to safety considerations of maintaining battery temperatures below their operating limit. Because of these considerations, battery models for preliminary design and optimization of aircraft should include the capability to accurately compute the temperature distribution within the battery pack. Furthermore, battery pack designs should be as light-weight as possible to maximize the pack energy density, while also considering battery temperature limits. Here, we demonstrate a simultaneous trajectory and shape optimization of a battery pack concept, using a transient heat transfer finite element model coupled with a time-varying cell-discharge battery model to provide this capability. Including the transient finite element problem in the loop enables accurate temperatures that can be passed back to the cell discharge model, while the cell discharge model can supply the finite element model with time-varying heat boundary conditions to the finite element problem, further benefiting the fidelity of the thermal response of the batteries. We first demonstrate the coupling capability between the battery cell-discharge model and the transient finite-element heat transfer through an optimization which computes the optimal current profile for the battery pack while ensuring the battery temperatures remain below their operational limit. We then build on this optimization by adding shape optimization to the problem, which allows us to consider a composite objective function which also minimizes the mass of the battery pack, while also producing an optimal current discharge profile.

Optimization↗
DOE OSTI2025

Microwave Debinding of Ceramics Produced by Additive Manufacturing

Microwave-assisted debinding offers an innovative approach to processing 3D-printed engineering ceramics. Ceramic parts produced with photopolymer-based additive manufacturing techniques require thermal debinding to remove polymers that bind ceramic particles during the printing process. Microwave-assisted debinding significantly reduces processing times and offers energy savings over conventional thermal treatments due to higher heating rates and more uniform heating. Finite-element modeling of the debinding process allows for the optimization and prediction of temperature distribution, stress development, and potential defects, leading to improved process control and material quality.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Minimization of thermal deformation in crystal optics for high repetition-rate FEL

Minimizing thermal deformation in X-ray crystal optics is crucial for preserving coherence and wavefront in high-repetition-rate free-electron lasers (FELs). This study presents two approaches to reduce pulse-by-pulse transient thermal deformation in diamond crystals used in cavity-based X-ray FELs (CBXFELs): (1) cryogenic cooling with liquid nitrogen (LN₂), and (2) second-order correction via focusing optics. We revisit the temperature-dependent thermal-mechanical properties of diamond and silicon, implement a finite-element analysis (FEA) method to accelerate convergence to a quasisteady-state regime. Results show that LN₂-cooled diamond crystals meet the stringent deformation requirement of less than 15 pm RMS for the pulse at mJ scale at 1 MHz repetition frequency, and up to 1.5 mJ for 100 kHz. Second-order correction by using focusing elements within the cavity can reduce the impact of thermal deformation for both liquid nitrogen and water cooling.

free electron laser↗
DOE OSTI2025

Novel Design and Fabrication of a High Frequency Transient Heat Flux Sensor for Use in an RDE

Rotating detonation engine (RDE) combustion systems have been a topic of interest in the pressure gain combustion community for their benefits over traditional gas turbine engine combustors. However, cooling requirements for these engines are significantly higher and less predictable than non-detonating engines. To understand the high-speed heat transfer dynamics inside an RDE, a novel, high-frequency heat flux gage is presented. This study aims to design a robust, single-sided sensor that can withstand the high temperature and harsh environment of an RDE for extended durations. Sensor bench testing is performed using a hot plate as a heat source, and the sensor response is compared to a finite-element analysis (FEA) model. The sensor response is then tested inside a water-cooled RDE and the wall heat flux is compared to calorimetry data.

rotating detonation engines↗
DOE OSTI2026

Broad frequency tuning of a Nb$_{3}$Sn superconducting microwave cavity for dark matter searches

We demonstrate a novel broad-frequency tuning mechanism for superconducting microwave cavities designed for dark matter searches. Using a Nb$_3$Sn-coated cigar-shaped cavity operating at approximately 9 GHz, we achieve continuous frequency tuning exceeding 1 GHz by mechanically separating the two cavity halves: a "tuning-by-opening" technique. Finite-element method simulations predict that radiative losses do not degrade the quality factor even for large openings, as a closed cavity with an intrinsic quality factor of $10^7$ maintains this value for apertures up to 9 mm, corresponding to a tuning range from 9.0 to 7.5 GHz. Experimental validation using both copper ring spacers and a continuous sliding mechanism confirms $Q_0$ values exceeding the dark matter quality factor across the entire explored frequency range, despite mechanical imperfections and film non-uniformities. This tuning approach avoids inserting elements into the resonant volume, making it particularly suitable for high-Q superconducting cavities in axion haloscope experiments and readily applicable to REBCO-based implementations capable of operating in multi-tesla magnetic fields.

Maiello, D. [Padua U.; INFN, Padua] (ORCID:0009000↗
DOE OSTI2026

Optimal Hairpin Winding Configuration for EV Traction Motors to Enhance Active CMV Cancellation Capability of NPL.X Inverter

Active common-mode voltage (CMV) cancellation in neutral-point-less (NPL.X) inverter is an effective approach for reducing common-mode electromagnetic interference and alleviating voltage stress associated with high-speed switching in EV traction drives. However, the existing hairpin winding configuration often exhibits impedance imbalance due to slot leakage inductance and manufacturing constraints, limiting the effectiveness of active CMV cancellation. This paper investigates symmetric hairpin winding configurations for a dual-three-phase EV traction machine to improve impedance symmetry. Considering the manufacturing constraints of six-layer hairpin windings, four feasible symmetric winding configurations are analyzed. Two-dimensional finite-element analysis is performed to extract the impedance parameters of each configuration. The extracted impedance parameters are used to quantify the impedance mismatch between complementary phase pairs and evaluate its influence on active CMV cancellation. The results show that the fully interleaved winding configuration achieves perfectly balanced phase impedances with 0 % mismatch, whereas the mismatch in the other configurations ranges from approximately 2.4 % to 0.8 %. Consequently, the resulting total CMV overshoot is significantly suppressed to near zero voltage under high-voltage and high switching frequency operating conditions, demonstrating the importance of winding symmetry for maximizing the active CMV cancellation capability of the NPL.X inverter.

Lee, Kangbeen [Purdue Univ., West Lafayette, IN (U↗
DOE OSTI2026

In-Situ Magnetic Field Reconstruction in the MAGIS-100 Experiment

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.; Ferm↗
DOE OSTI2025

Physics-based stabilized finite element approximations of the Poisson–Nernst–Planck equations

We present and analyze two stabilized finite element methods for solving numerically the Poisson–Nernst–Planck equations. The stabilization we consider is carried out by using a shock detector and a discrete graph Laplacian operator for the ion equations, whereas the discrete equation for the electric potential need not be stabilized. Discrete solutions stemmed from the first algorithm preserve both maximum and minimum discrete principles. For the second algorithm, its discrete solutions are conceived so that they hold discrete principles and obey an entropy law provided that an acuteness condition is imposed for meshes. Remarkably the latter is found to be unconditionally stable. We validate our methodology through transient numerical experiments that show convergence toward steady-state solutions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Polynomial range estimation as a troubled-cell indicator for high-order methods

Two troubled-cell indicators based on polynomial range estimation methods are used to flag cells that may violate positivity constraints. One method uses interval extension, and the second uses the range enclosure property of the Bernstein polynomial basis. Furthermore, both methods reduce compute time for the positivity preserver by limiting its application to a subset of cells. The Bernstein polynomial method remains effective as the problem dimensionality increases. Interval extension applied to the internal energy equation permits the use of the troubled-cell indicators for rational functions, though performance suffers compared to directly applying the indicators to polynomial functions.

42 ENGINEERING↗
DOE OSTI2025

X-point effects on the ideal MHD modes in tokamaks in the description of dual-poloidal-region safety factor

The flux coordinates with dual-region safety factor (q) in the poloidal direction are developed in this work. The X-point effects on the ideal MHD modes in tokamaks are then analyzed using this coordinate system. Since the X-point effects mainly affect the edge region, the modes localized at the tokamak edge are particularly examined. Two types of modes are studied. The first is related to the conventional peeling or peeling-ballooning modes. The mode existence aligned with the local magnetic field in the poloidally core region as observed experimentally is confirmed. The X points are shown to contribute to a stabilizing effect for the conventionally treated modes with the surface-averaged q and with the tokamak edge portion truncated. The other is the axisymmetric modes localized in the vicinity of X points, which can affect the cross-field-line transport near the X points. The existence of axisymmetric modes points to the possibility of applying a toroidally axisymmetric resonant magnetic perturbation (RMP) in the X-point area for mitigating the edge localized modes, which can be an alternative to the current RMP design. The dual q description also has important implications for the existing non-axisymmetric RMP concept. It helps to understand why the RMP suppression of edge localized modes is difficult to achieve in the double-null tokamak configurations and points to the possibility of further improving the current RMP concept by considering the alignment to the local q.

Fourier analysis↗
DOE OSTI2024

Full-field quantitative visualization of shock-driven pore collapse and failure modes in PMMA

The dynamic collapse of pores under shock loading is thought to be directly related to hot spot generation and material failure, which is critical to the performance of porous energetic and structural materials. However, the shock compression response of porous materials at the local, individual pore scale is not well understood. This study examines, quantitatively, the collapse phenomenon of a single spherical void in PMMA at shock stresses ranging from 0.4 to 1.0 GPa. Using a newly developed internal digital image correlation technique in conjunction with plate impact experiments, full-field quantitative deformation measurements are conducted in the material surrounding the collapsing pore for the first time. The experimental results reveal two failure mode transitions as shock stress is increased: (i) the first in situ evidence of shear localization via adiabatic shear banding and (ii) dynamic fracture initiation at the pore surface. Numerical simulations using thermo-viscoplastic dynamic finite element analysis provide insights into the formation of adiabatic shear bands (ASBs) and stresses at which failure mode transitions occur. Further numerical and theoretical modeling indicates the dynamic fracture to occur along the weakened material inside an adiabatic shear band. Finally, analysis of the evolution of pore asymmetry and models for ASB spacing elucidate the mechanisms for the shear band initiation sites, and elastostatic theory explains the experimentally observed ASB and fracture paths based on the directions of maximum shear.

42 ENGINEERING↗