Search NASA⌕ Search

SEARCH · Search NASA

Results for “extended finite element method”

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 289 records · Page 16

Monitoring of contact stresses in advanced propulsion systems

This paper discusses a new and innovative numerical procedure for evaluating and monitoring contact stresses in the mating surfaces of propulsion systems. The procedure is intended to extend the range of applicability and the accuracy of contact stress analyses. Originally proposed by Burton Paul, et. al., (1981) the procedure is based upon fundamental solutions from the theory of elasticity. It is an iterative numerical procedure. It is believed that the method has distinct advantages over the classical Hertz method, the finite-element method and over existing approaches with the boundary element method. Several examples are discussed. Favorable comparisons are obtained with classical results. Applications with mating element of propulsion systems are discussed.

Somprakit, Paisan↗

A weighted shifted boundary method for immersed moving boundary simulations of Stokes' flow

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. The surrogate domain is constructed so as to avoid cut cells and the associated problematic implementation and numerical integration issues. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions: hence the name of the method, that shifts the location and values of the boundary conditions. Here, in this article, we extend the SBM to the simulation of incompressible Stokes flow, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach allows to drastically reduce spurious pressure oscillations in time, which are produced if the total volume of active fluid were to change abruptly over a time step. The proposed Weighted SBM (W-SBM) exactly preserves states of hydrostatic equilibrium, and induces small mass and momentum conservation errors, which converge as the grid is refined. This is in analogy to cutFEMs and related unfitted approaches, which rely on an affine representation of cut boundaries. We demonstrate the robustness and accuracy of the proposed method with an extensive suite of two-dimensional tests.

97 MATHEMATICS AND COMPUTING↗

Scale–Resolving Simulations of a Fundamental Trailing–Edge Cooling Slot Using a Discontinuous–Galerkin Spectral–Element Method

The accurate prediction of turbulent mixing in high-pressure turbines that incorporate various airfoil surface-cooling strategies is becoming increasing critical to the design of modern gas turbine engines where the quest for improved efficiency is driving compressor overall pressure ratios and turbine inlet temperatures to much higher levels than ever before. In the present paper, a recently developed computational capability for accurate and efficient scale–resolving simulations of turbomachinery is extended to study the turbulent mixing mechanism of a simplified abstraction of an airfoil trailing-edge cooling slot - a plane wall jet with finite lip thickness discharging into an ambient flow. The computational capability is based on an entropy– stable, discontinuous–Galerkin approach that extends to arbitrarily high orders of spatial and temporal accuracy. The numerical results show that the present simulations capture the trends observed in the experiments. Discrepancies between the simulations and experiments are believed to be due to differences in the inflow profiles and tunnel side–wall effects. The thick lip configuration leads to a thicker wake and higher unsteadiness in the wall jet compared to the thin lip. A detailed comparison of the turbulent flowfields is presented to highlight differences arising due to lip thickness variations.

Scale Resolving↗

Numerical studies of nonlinear ultrasonic guided waves in uniform waveguides with arbitrary cross sections

Nonlinear guided waves have been investigated widely in simple geometries, such as plates, pipe and shells, where analytical solutions have been developed. This paper extends the application of nonlinear guided waves to waveguides with arbitrary cross sections. The criteria for the existence of nonlinear guided waves were summarized based on the finite deformation theory and nonlinear material properties. Numerical models were developed for the analysis of nonlinear guided waves in complex geometries, including nonlinear Semi-Analytical Finite Element (SAFE) method to identify internal resonant modes in complex waveguides, and Finite Element (FE) models to simulate the nonlinear wave propagation at resonant frequencies. Two examples, an aluminum plate and a steel rectangular bar, were studied using the proposed numerical model, demonstrating the existence of nonlinear guided waves in such structures and the energy transfer from primary to secondary modes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Determining failure properties of as-received and hydrided unirradiated Zircaloy-4 from ring compression tests

An integrated experimental and computational approach was developed to determine the plastic strain at fracture and to investigate the fracture energy of as-received and hydrided unirradiated Zircaloy-4 nuclear cladding tubes in the hoop direction at room temperature from ring compression testing (RCT). This work builds on previous methods in the literature that were developed to obtain the mechanical properties of nuclear cladding before fracture (Young’s modulus, yield stress, and strain-hardening parameters) from experimental RCT results and extends the characterization of this material class to include material failure properties. Furthermore, a ductile damage approach was developed and implemented in a validated finite element model to predict material failure, evaluate fracture energy, and quantify the plastic strain at fracture initiation. The hydrogen content of the Zircaloy-4 specimens was varied up to 610 wppm when the failure became brittle. As expected, increasing hydrogen content caused embrittlement of the Zircaloy-4 samples tested and the plastic strain at fracture and fracture energy to decrease.

36 MATERIALS SCIENCE↗

RF wave simulation for cold edge plasmas using the MFEM library

A newly developed generic electro-magnetic (EM) simulation tool for modeling RF wave propagation in SOL plasmas is presented. The primary motivation of this development is to extend the domain partitioning approach for incorporating arbitrarily shaped SOL plasmas and antenna to the TORIC core ICRF solver, which was previously demonstrated in the 2D geometry [S. Shiraiwa, et. al., “HISTORIC: extending core ICRF wave simulation to include realistic SOL plasmas”, Nucl. Fusion in press], to larger and more complicated simulations by including a 3D realistic antenna and integrating RF rectified sheath potential model. Such an extension requires a scalable high fidelity 3D edge plasma wave simulation. We used the MFEM [http://mfem.org], open source scalable C++ finite element method library, and developed a Python wrapper for MFEM (PyMFEM), and then a radio frequency (RF) wave physics module in Python. This approach allows for building a physics layer rapidly, while separating the physics implementation being apart from the numerical FEM implementation. An interactive modeling interface was built on pScope [S Shiraiwa, et. al. Fusion Eng. Des. 112, 835] to work with an RF simulation model in a complicated geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Advanced electric motor technology flux mapping

Design of electric motors which fulfill the needs of Thrust Vector Control (TVC) actuators used in large rocket propelled launch vehicles is covered. To accomplish this end the methodology of design is laid out in some detail. In addition a point design of a motor to fulfill the requirements of a certain actuator specified by MSFC is accomplished and reported upon. In the course of this design great stress has been placed on ridding the actuator of internally generated heat. To conduct the heat out of the motor use is made of the unique properties of the in house MSFC designed driving electronics. This property is that as along as they are operated in a quasi-linear mode the electronics nullify the effects of armature inductance as far as the phase of the armature current versus the rotor position is concerned. Actually the additional inductance due to the extended end turns in this design is of benefit because in the shorted armature failure mode the armature current in the fault (caused by the rotor flux sweeping past the armature) is diminished at a given rotor speed when compared to a more conventional motor with lower inductance. The magnetic circuit is analyzed using electromagnetic finite element methods.

Doane, George B., III↗

Probabilistic finite elements for transient analysis in nonlinear continua

The probabilistic finite element method (PFEM), which is a combination of finite element methods and second-moment analysis, is formulated for linear and nonlinear continua with inhomogeneous random fields. Analogous to the discretization of the displacement field in finite element methods, the random field is also discretized. The formulation is simplified by transforming the correlated variables to a set of uncorrelated variables through an eigenvalue orthogonalization. Furthermore, it is shown that a reduced set of the uncorrelated variables is sufficient for the second-moment analysis. Based on the linear formulation of the PFEM, the method is then extended to transient analysis in nonlinear continua. The accuracy and efficiency of the method is demonstrated by application to a one-dimensional, elastic/plastic wave propagation problem. The moments calculated compare favorably with those obtained by Monte Carlo simulation. Also, the procedure is amenable to implementation in deterministic FEM based computer programs.

Liu, W. K.↗

Enhanced accuracy through ensembling of randomly initialized auto-regressive models for dynamical systems

Computational mechanics simulations using traditional finite element methods (FEM) require prohibitively expensive computational resources for real-time engineering applications, design optimization, and digital twin implementations. While machine learning (ML) surrogate models offer significant computational speedups, autoregressive ML models for time-dependent mechanical systems suffer from error accumulation that compromises long-term prediction reliability - a critical concern for engineering applications where accuracy over extended time horizons is essential for safety and performance assessments. Here, we propose a deep ensemble framework specifically designed to address this challenge in computational mechanics applications, where multiple ML surrogate models with random weight initializations are trained in parallel and their predictions aggregated during inference. This approach leverages statistical diversity to maximize information gain from a fixed set of training data and to mitigate error propagation, while maintaining the computational efficiency that makes ML surrogates attractive for engineering practice. We validate the framework on three representative problems spanning critical areas of computational mechanics: stress field evolution in heterogeneous microstructures under complex loading (relevant to advanced materials design and composite analysis), planetary-scale shallow water dynamics (applicable to environmental and geotechnical engineering), and Gray-Scott reaction-diffusion systems (relevant to mass transport and chemical process engineering). Across all test cases, the ensemble approach demonstrates consistent error reduction of 15-33% compared to individual models. The codes for this work are available on GitHub (https://github.com/Graham-Brady-Research-Group/AutoregressiveEnsemble_SpatioTemporal_Evolution).

autoregressive prediction↗

Compositional modeling in porous media using constant volume flash and flux computation without the need for phase identification

The paper deals with the numerical solution of a compositional model describing compressible two-phase flow of a mixture composed of several components in porous media with species transfer between the phases. The mathematical model is formulated by means of the extended Darcy's laws for all phases, components continuity equations, constitutive relations, and appropriate initial and boundary conditions. The splitting of components among the phases is described using a new formulation of the local thermodynamic equilibrium which uses volume, temperature, and moles as specification variables. The problem is solved numerically using a combination of the mixed-hybrid finite element method for the total flux discretization and the finite volume method for the discretization of transport equations. A new approach to numerical flux approximation is proposed, which does not require the phase identification and determination of correspondence between the phases on adjacent elements. The time discretization is carried out by the backward Euler method. The resulting large system of nonlinear algebraic equations is solved by the Newton–Raphson iterative method. We provide eight examples of different complexity to show reliability and robustness of our approach.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Semiclassical states on Lie algebras

The effective technique for analyzing representation-independent features of quantum systems based on the semiclassical approximation (developed elsewhere) has been successfully used in the context of the canonical (Weyl) algebra of the basic quantum observables. Here, we perform the important step of extending this effective technique to the quantization of a more general class of finite-dimensional Lie algebras. The case of a Lie algebra with a single central element (the Casimir element) is treated in detail by considering semiclassical states on the corresponding universal enveloping algebra. Restriction to an irreducible representation is performed by “effectively” fixing the Casimir condition, following the methods previously used for constrained quantum systems. We explicitly determine the conditions under which this restriction can be consistently performed alongside the semiclassical truncation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Effect of type of load on stress analysis of thin-walled ducts

The standard procedure for qualifying the design of duct (pipe) systems in the Space Shuttle Main Engine (SSME) has been fairly well defined. However, since pipe elbows are quite common and important in the SSME duct systems, a clear understanding of the detailed stress profile of the components is necessary for accurate structural and life assessments. This study was initiated to predict the stress profile at/near the tangent point along the cross section of the duct under various types of loads. Also, this study was further extended to understand the stiffening effect on stresses due to pressure at the tangent point. The intention of this study was to identify the importance of selecting proper locations for mounting strain gauges and to utilize the obtained results to anchor dynamic models for accurate structural and life assessments of the SSME ducts under a dynamic environment. The finite element method was utilized in this study.

Min, J. B.↗

Recent developments of the NESSUS probabilistic structural analysis computer program

The NESSUS probabilistic structural analysis computer program combines state-of-the-art probabilistic algorithms with general purpose structural analysis methods to compute the probabilistic response and the reliability of engineering structures. Uncertainty in loading, material properties, geometry, boundary conditions and initial conditions can be simulated. The structural analysis methods include nonlinear finite element and boundary element methods. Several probabilistic algorithms are available such as the advanced mean value method and the adaptive importance sampling method. The scope of the code has recently been expanded to include probabilistic life and fatigue prediction of structures in terms of component and system reliability and risk analysis of structures considering cost of failure. The code is currently being extended to structural reliability considering progressive crack propagation. Several examples are presented to demonstrate the new capabilities.

Millwater, H.↗

High-Fidelity Micromechanics Model Enhanced for Multiphase Particulate Materials

This 3-year effort involves the development of a comprehensive micromechanics model and a related computer code, capable of accurately estimating both the average response and the local stress and strain fields in the individual phases, assuming both elastic and inelastic behavior. During the first year (fiscal year 2001) of the investigation, a version of the model called the High-Fidelity Generalized Method of Cells (HFGMC) was successfully completed for the thermo-inelastic response of continuously reinforced multiphased materials with arbitrary periodic microstructures (refs. 1 and 2). The model s excellent predictive capability for both the macroscopic response and the microlevel stress and strain fields was demonstrated through comparison with exact analytical and finite element solutions. This year, HFGMC was further extended in two technologically significant ways. The first enhancement entailed the incorporation of fiber/matrix debonding capability into the two-dimensional version of HFGMC for modeling the response of unidirectionally reinforced composites such as titanium matrix composites, which exhibit poor fiber/matrix bond. Comparison with experimental data validated the model s predictive capability. The second enhancement entailed further generalization of HFGMC to three dimensions to enable modeling the response of particulate-reinforced (discontinuous) composites in the elastic material behavior domain. Next year, the three-dimensional version will be generalized to encompass inelastic effects due to plasticity, viscoplasticity, and damage, as well as coupled electromagnetothermomechanical (including piezoelectric) effects.

Pindera, Marek-Jerzy↗

Methods Developed by the Tools for Engine Diagnostics Task to Monitor and Predict Rotor Damage in Real Time

Tools for Engine Diagnostics is a major task in the Propulsion System Health Management area of the Single Aircraft Accident Prevention project under NASA s Aviation Safety Program. The major goal of the Aviation Safety Program is to reduce fatal aircraft accidents by 80 percent within 10 years and by 90 percent within 25 years. The goal of the Propulsion System Health Management area is to eliminate propulsion system malfunctions as a primary or contributing factor to the cause of aircraft accidents. The purpose of Tools for Engine Diagnostics, a 2-yr-old task, is to establish and improve tools for engine diagnostics and prognostics that measure the deformation and damage of rotating engine components at the ground level and that perform intermittent or continuous monitoring on the engine wing. In this work, nondestructive-evaluation- (NDE-) based technology is combined with model-dependent disk spin experimental simulation systems, like finite element modeling (FEM) and modal norms, to monitor and predict rotor damage in real time. Fracture mechanics time-dependent fatigue crack growth and damage-mechanics-based life estimation are being developed, and their potential use investigated. In addition, wireless eddy current and advanced acoustics are being developed for on-wing and just-in-time NDE engine inspection to provide deeper access and higher sensitivity to extend on-wing capabilities and improve inspection readiness. In the long run, these methods could establish a base for prognostic sensing while an engine is running, without any overt actions, like inspections. This damage-detection strategy includes experimentally acquired vibration-, eddy-current- and capacitance-based displacement measurements and analytically computed FEM-, modal norms-, and conventional rotordynamics-based models of well-defined damages and critical mass imbalances in rotating disks and rotors.

Baaklini, George Y.↗

A general mapping procedure for variable area duct acoustics

A general mapping procedure is described and applied to the study of noise propagation in variable area ducts. The mapping provides a boundary fitted co-ordinate system which is ideal for the finite difference solution of acoustic fields with irregular boundaries, without the burden of large matrices required by finite element methods. The procedure is first described in general and then applied to a particular two-dimensional geometry under current experimental investigation. This method should be ideally suited to the study of high frequency noise propagation in variable area ducts and in cases where the far field is included in the calculation procedure. Moreover, the current approach can be directly extended to three-dimensions, resulting in numerical calculation over a rectangular parallelepiped in the transformed plane.

White, J. W.↗

Numerical Differentiation for Adaptively Refined Finite Element Meshes

Postprocessing of point-wise data is a fundamental process in many fields of research. Numerical differentiation is a key operation in computational electromagnetics. In the case of data obtained from a finite element method with automatic mesh refinement much work needs still to be done. This paper addresses some issues in differentiating data obtained from a finite element electromagnetic code with adaptive mesh refinement, and it proposes a methodology for deriving the electric field given the magnetic field on a mesh of linear triangular elements. The procedure itself is nevertheless more general and might be extended for numerically differentiating any point-wise solution based on triangular meshes.

Borgioli, Andrea↗

Finite element analysis of an extended end-plate connection using the T-stub approach

Beam-to-column end-plate bolted connections are usually used as moment-resistant connections in steel framed structures. For this joint type, the deformability is governed by the deformation capacity of the column flange and end-plate under tension and elongation of the bolts. All these elements around the beam tension flange form the tension region of the joint, which can be modeled by means of equivalent T-stubs. In this paper a beam-to-column end-plate bolted connection is substituted with a T-stub of appropriate effective length and it is analyzed using the commercially available finite element software ABAQUS. The performance of the model is validated by comparing the behavior of the T-stub from the numerical simulation with the behavior of the connection as a whole. The moment-rotation curve of the T-stub obtained from the numerical simulation is compared with the behavior of the whole extended end-plate connection, obtained by numerical simulation, experimental tests and analytical approach.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗