Search NASA⌕ Search

SEARCH · Search NASA

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A Unified Design Theory for Multi-Port Polyphase Transformers Enabling Scalable Power-Multiplexed EV Fleet Charging Systems

This paper presents a unified analytical design theory for multi-port polyphase transformers, targeting scalable and isolated high-power Electric Vehicle (EV) fleet charging systems with power multiplexing capability. As fleet electrification accelerates, conventional one-to-one charger architectures face significant challenges in infrastructure cost, peak power demand, and low utilization of installed power electronics. Power-multiplexed charging architectures, which dynamically distribute power from a shared pool of converter modules across multiple vehicles, have emerged as a promising solution. However, such architectures require scalable, isolated multi-port power interfaces capable of routing energy among multiple inputs and outputs, whose design remains complex and dependent on iterative modeling. To address this gap, the proposed theory provides closed-form expressions for self-inductance, leakage inductance, and mutual coupling terms for arbitrary multi-phase, multi-port transformer structures. The formulation enables direct synthesis of isolated multi-input and multi-output resonant converter systems without reliance on geometry-specific finite-element analysis or extensive parameter extraction. This capability is particularly critical for power-multiplexed systems, where modular converter structures must interface with multiple vehicles while maintaining galvanic isolation and flexible power allocation. The effectiveness of the proposed framework is demonstrated through the design of a 360 kW multi-phase system operating over a 700–900 VDC input and 400–1250 VDC output range. PLECS simulation results confirm accurate prediction of system behavior and validate the applicability of the approach to multi-port, power-multiplexed charging scenarios. The proposed method significantly reduces design complexity while enabling scalable, cost-effective, and fully utilized EV fleet charging infrastructure.

Asa, Erdem [ORNL] (ORCID:0000000190884812)↗
DOE OSTI2025

Automated Redox Titrations via Interdigitated Electrode Arrays: Application to the Mediated Electron Transfer Interrogation of Charge and Rate on Electrodeposited Polymers

Mediated electron transfer (MET) plays a crucial role in energy storage and conversion technologies such as redox targeting flow batteries (RTFBs), yet its experimental investigation often requires labor-intensive and low-throughput setups. To address this, we developed a microfabricated interdigitated electrode array (IDA) platform that enables automated, high-throughput electrochemical redox titration measurement to be performed to study the MET process. Our redox titration method enables simultaneous measurement of the charge capacity and rate of MET processes on a material or surface. Automated redox titration (ART) facilitates systematic investigation of the MET process across a broad parameter space, exemplified through the study of polypyrrole (PPy) and a pyrene-4,5,9,10-tetrone azo group-based polymer (PTAP), both redox-active polymers relevant to various energy storage applications. Using PPy as a model material, 500 redox titration measurements were conducted within 50 h, varying the electrode gap widths, polymer charging potentials, voltammetric scan rates, and electrolyte concentrations. Finite-element simulations confirmed the electrochemical responses and elucidated the kinetics of the MET reactions. Our automated methodology was further tested with PTAP, revealing a surprising charging potential dependence on the rate of MET. The automation, flexibility, and scalability of our redox titration platform pave the way not only for advanced studies of MET processes relevant to RTFBs, but also with implications in the understanding of next-generation energy storage materials, molecular electrocatalysis, and biosensing.

electrochemical analysis↗
NASA NTRS1976

Optimization of continuous one-dimensional structures under steady harmonic excitation

The paper examines some questions relating to the optimal design of continuous one-dimensional structures driven by harmonically oscillating loads. Optimal-control methods are applied to a cantilever bar driven sinusoidally by an axial force at its tip to illustrate the minimum-weight design of one-dimensional structures under dynamic excitation. Realistic constraints are imposed during the optimizations, including a maximum allowable stress amplitude at any point along the bar and a minimum cross-sectional area. It is shown that in the absence of damping, the design space may contain many disjoint feasible regions, and multiple optima can exist. Detailed solutions are obtained for continuous bars with an excitation frequency less than, and then greater than, the fundamental free-vibration frequency. It is found that above a certain excitation frequency, two or more arcs with different constraints characterize the optimal designs. It is concluded that when more than two different constrained arcs characterize the optimal solution, the continuum approach may be impractical, and finite-element approximations may offer the only alternative.

Johnson, E. H.↗
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

Benchmarking of massively parallel phase-field codes for directional solidification

We present a detailed benchmark comparing two state-of-the-art phase-field implementations for simulating alloy solidification under experimentally relevant conditions. The study investigates the directional solidification of Al-3wt%Cu under high-velocity solidification conditions and SCN-0.46wt% camphor under microgravity conditions from National Aeronautics and Space Administration (NASA) DECLIC-DSI-R experiments. Both codes, one employing finite-difference discretization with uniform mesh and GPU-acceleration (GPU-PF) and the other one employing finite-element discretization with adaptive-mesh and CPU-parallelization (PRISMS-PF), solve the same quantitative phase-field formulation that incorporates an anti-trapping current for the solidification of dilute alloys. We evaluate the predictions of each code for dendritic morphology, primary spacing, and tip dynamics in both 2D and 3D, as well as their numerical convergence and computational performance. While existing benchmark problems have primarily focused on simplified or small-scale simulations, they do not reflect the computational and modeling challenges posed by employing experimentally relevant time and length scales. Our results provide a practical framework for assessing phase-field code performance as well as validating and facilitating their application in integrated computational materials engineering (ICME) workflows that require integration with realistic experimental data.

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

Two-dimensional finite-element temperature variance analysis

The finite element method is extended to thermal analysis by forming a variance analysis of temperature results so that the sensitivity of predicted temperatures to uncertainties in input variables is determined. The temperature fields within a finite number of elements are described in terms of the temperatures of vertices and the variational principle is used to minimize the integral equation describing thermal potential energy. A computer calculation yields the desired solution matrix of predicted temperatures and provides information about initial thermal parameters and their associated errors. Sample calculations show that all predicted temperatures are most effected by temperature values along fixed boundaries; more accurate specifications of these temperatures reduce errors in thermal calculations.

Heuser, J. S.↗
NASA NTRS1972

Computer program: Jet 3 to calculate the large elastic plastic dynamically induced deformations of free and restrained, partial and/or complete structural rings

A user-oriented FORTRAN 4 computer program, called JET 3, is presented. The JET 3 program, which employs the spatial finite-element and timewise finite-difference method, can be used to predict the large two-dimensional elastic-plastic transient Kirchhoff-type deformations of a complete or partial structural ring, with various support conditions and restraints, subjected to a variety of initial velocity distributions and externally-applied transient forcing functions. The geometric shapes of the structural ring can be circular or arbitrarily curved and with variable thickness. Strain-hardening and strain-rate effects of the material are taken into account.

Wu, R. W.↗
NASA NTRS1974

Development of solution techniques for nonlinear structural analysis

Nonlinear structural solution methods in the current research literature are classified according to order of the solution scheme, and it is shown that the analytical tools for these methods are uniformly derivable by perturbation techniques. A new perturbation formulation is developed for treating an arbitrary nonlinear material, in terms of a finite-difference generated stress-strain expansion. Nonlinear geometric effects are included in an explicit manner by appropriate definition of an applicable strain tensor. A new finite-element pilot computer program PANES (Program for Analysis of Nonlinear Equilibrium and Stability) is presented for treatment of problems involving material and geometric nonlinearities, as well as certain forms on nonconservative loading.

Vos, R. G.↗
NASA NTRS1975

Application of NASTRAN for stress analysis of left ventricle of the heart

Knowing the stress and strain distributions in the left ventricular wall of the heart is a prerequisite for the determination of the muscle elasticity and contractility in the process of assessing the functional status of the heart. NASTRAN was applied for the calculation of these stresses and strains and to help in verifying the results obtained by the computer program FEAMPS which was specifically designed for the plane-strain finite-element analysis of the left ventricular cross sections. Adopted for the analysis are the true shape and dimensions of the cross sections reconstructed from multiplanar X-ray views of a left ventricle which was surgically isolated from a dog's heart but metabolically supported to sustain its beating. A preprocessor was prepared to accommodate both FEAMPS and NASTRAN, and it has also facilitated the application of both the triangular element and isoparameteric quadrilateral element versions of NASTRAN. The stresses in several crucial regions of the left ventricular wall calculated by these two independently developed computer programs are found to be in good agreement. Such confirmation of the results is essential in the development of a method which assesses the heart performance.

Pao, Y. C.↗
NASA NTRS1974

Finite-element analysis of vibration and flutter of cantilever anisotropic plates

The hybrid stress finite element method is used to study the effects of filament angle and orthotropicity parameter on the vibration and flutter characteristics of cantilevered anisotropic plates. The results indicate a generally strong lack of monotonic dependence on filament angle. Also, the natural frequency in certain cases involving the first few modes, can become relatively insensitive to both filament angle and orthotropicity parameter for a range of filament angle beyond 70 deg. Values of critical dynamic pressure are obtained by a modal approach, in which the mode shapes are obtained by the hybrid stress method. Convergence of the modal method is rather rapid for the configurations analyzed, and a comparison of the method with an exact solution for the case of an isotropic simply supported plate shows that use of six modes gives excellent agreement.

Rossettos, J. N.↗
NASA NTRS1975

Finite-element formulations for problems of large elastic-plastic deformation

An Eulerian finite element formulation is presented for problems of large elastic-plastic flow. The method is based on Hill's variational principle for incremental deformations, and is ideally suited to isotropically hardening Prandtl-Reuss materials. Further, the formulation is given in a manner which allows any conventional finite element program, for 'small strain' elastic-plastic analysis, to be simply and rigorously adapted to problems involving arbitrary amounts of deformation and arbitrary levels of stress in comparison to plastic deformation moduli. The method is applied to a necking bifurcation analysis of a bar in plane-strain tension. The paper closes with a unified general formulation of finite element equations, both Lagrangian and Eulerian, for large deformations, with arbitrary choice of the conjugate stress and strain measures. Further, a discussion is given of other proposed formulations for elastic-plastic finite element analysis at large strain, and the inadequacies of some of these are commented upon.

Mcmeeking, R. M.↗
NASA NTRS1976

Finite-element analysis of dynamic fracture

Applications of the finite element method to the two dimensional elastodynamics of cracked structures are presented. Stress intensity factors are computed for two problems involving stationary cracks. The first serves as a vehicle for discussing lumped-mass and consistent-mass characterizations of inertia. In the second problem, the behavior of a photoelastic dynamic tear test specimen is determined for the time prior to crack propagation. Some results of a finite element simulation of rapid crack propagation in an infinite body are discussed.

Aberson, J. A.↗
NASA NTRS1977

Three dimensional finite-element analysis of finite-thickness fracture specimens

The stress-intensity factors for most of the commonly used fracture specimens (center-crack tension, single and double edge-crack tension, and compact), those that have a through-the-thickness crack, were calculated using a three dimensional finite-element elastic stress analysis. Three-dimensional singularity elements were used around the crack front. The stress intensity factors along the crack front were evaluated by using a force method, developed herein, that requires no prior assumption of either plane stress or plane strain. The calculated stress-intensity factors from the present analysis were compared with those from the literature whenever possible and were generally found to be in good agreement. The stress-intensity factors at the midplane for all specimens analyzed were within 3 percent of the two dimensional plane strain values. The stress intensity factors at the specimen surfaces were considerably lower than at the midplanes. For the center-crack tension specimens with large thickness to crack-length ratios, the stress-intensity factor reached a maximum near the surface of the specimen. In all other specimens considered the maximum stress intensity occurred at the midplane.

Raju, I. S.↗
NASA NTRS1977

Simplified and refined structural modeling for economical flutter analysis and design

A coordinated use of two finite-element models of different levels of refinement is presented to reduce the computer cost of the repetitive flutter analysis commonly encountered in structural resizing to meet flutter requirements. One model, termed a refined model (RM), represents a high degree of detail needed for strength-sizing and flutter analysis of an airframe. The other model, called a simplified model (SM), has a relatively much smaller number of elements and degrees-of-freedom. A systematic method of deriving an SM from a given RM is described. The method consists of judgmental and numerical operations to make the stiffness and mass of the SM elements equivalent to the corresponding substructures of RM. The structural data are automatically transferred between the two models. The bulk of analysis is performed on the SM with periodical verifications carried out by analysis of the RM. In a numerical example of a supersonic cruise aircraft with an arrow wing, this approach permitted substantial savings in computer costs and acceleration of the job turn-around.

Ricketts, R. H.↗