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 1,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Finite Element Method for Capturing Ultra-relativistic Shocks

While finite element methods are used extensively by researchers solving computational fluid dynamics in fields other than astrophysics, their use in astrophysical fluid simulations has been predominantly overlooked. Current simulations using other methods such as finite difference and finite volume (based on finite difference) have shown remarkable results, but these methods are limited by their fundamental properties in aspects that are important for simulations with complex geometries and widely varying spatial and temporal scale differences. We have explored the use of finite element methods for astrophysical fluids in order to establish the validity of using such methods in astrophysical environments. We present our numerical technique applied to solving ultra-relativistic (Lorentz Factor Gamma >> 1) shocks which are prevalent in astrophysical studies including relativistic jets and gamma-ray burst studies. We show our finite element formulation applied to simulations where the Lorentz factor ranges up to 2236 and demonstrate its stability in solving ultra-relativistic flows. Our numerical method is based on the Flowfield Dependent Variation (FDV) Method, unique in that numerical diffusion is derived from physical parameters rather than traditional artificial viscosity methods. Numerical instabilities account for most of the difficulties when capturing shocks in this regime. Our method results in stable solutions and accurate results as compared with other methods.

Richardson, G. A.↗
NASA NTRS2003

Application of Probabilistic Analysis to Aircraft Impact Dynamics

Full-scale aircraft crash simulations performed with nonlinear, transient dynamic, finite element codes can incorporate structural complexities such as: geometrically accurate models; human occupant models; and advanced material models to include nonlinear stressstrain behaviors, laminated composites, and material failure. Validation of these crash simulations is difficult due to a lack of sufficient information to adequately determine the uncertainty in the experimental data and the appropriateness of modeling assumptions. This paper evaluates probabilistic approaches to quantify the uncertainty in the simulated responses. Several criteria are used to determine that a response surface method is the most appropriate probabilistic approach. The work is extended to compare optimization results with and without probabilistic constraints.

Lyle, Karen H.↗
NASA NTRS2008

Effect of Microscopic Damage Events on Static and Ballistic Impact Strength of Triaxial Braid Composites

In previous work, the ballistic impact resistance of triaxial braided carbon/epoxy composites made with large flat tows (12k and 24k) was examined by impacting 2 X2 X0.125" composite panels with gelatin projectiles. Several high strength, intermediate modulus carbon fibers were used in combination with both untoughened and toughened matrix materials. A wide range of penetration thresholds were measured for the various fiber/matrix combinations. However, there was no clear relationship between the penetration threshold and the properties of the constituents. During some of these experiments high speed cameras were used to view the failure process, and full-field strain measurements were made to determine the strain at the onset of failure. However, these experiments provided only limited insight into the microscopic failure processes responsible for the wide range of impact resistance observed. In order to investigate potential microscopic failure processes in more detail, quasi-static tests were performed in tension, compression, and shear. Full-field strain measurement techniques were used to identify local regions of high strain resulting from microscopic failures. Microscopic failure events near the specimen surface, such as splitting of fiber bundles in surface plies, were easily identified. Subsurface damage, such as fiber fracture or fiber bundle splitting, could be identified by its effect on in-plane surface strains. Subsurface delamination could be detected as an out-of-plane deflection at the surface. Using this data, failure criteria could be established at the fiber tow level for use in analysis. An analytical formulation was developed to allow the microscopic failure criteria to be used in place of macroscopic properties as input to simulations performed using the commercial explicit finite element code, LS-DYNA. The test methods developed to investigate microscopic failure will be presented along with methods for determining local failure criteria that can be used in analysis. Results of simulations performed using LS-DYNA will be presented to illustrate the capabilities and limitations for simulating failure during quasi-static deformation and during ballistic impact of large unit cell size triaxial braid composites.

Littell, Justin D.↗
DOE OSTI2016

Thermo-mechanical simulations of CO{sub 2} laser–fused silica interactions

CO{sub 2} laser heating of silica glass is used in many scientific and industrial applications. Particularly, localized CO{sub 2} laser heating of silica glass has demonstrated its ability to mitigate surface damage on optics used for high power laser applications. To develop such applications, the control of temperature, heat affected area, and resulting mechanical stresses are critical. Therefore, it is necessary to understand the silica transformation, the material ejection, and the thermo-mechanical stresses induced by the laser heating and subsequent cooling. In this paper, we detail the development of comprehensive thermo-mechanical numerical simulations of these physical processes, based on finite-element method. The approach is developed for 2D or 3D cases to tackle the case of a moving beam at the surface of the sample, and we particularly discuss the choice of the different parameters based on bibliographic inputs. The thermal and mechanical numerical results have been compared to different dedicated experimental studies: infrared thermography measurements at the surface of the irradiated area, optical profilometry measurements of the laser-processed sites, and photo-elastic measurements. Very consistent results are obtained between numerical and experimental results for the description of the temperature gradients, the material ejection, and the residual stresses.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1989

Numerical simulation of controlled large space structures

Large Space Structures do not have much damping, which necessitates the installation of a controller onto the structure. If the controller is improperly designed, the structure may become unstable and be destroyed. Since Large Space Structures are extremely expensive pieces of hardware, new controllers must not be tested first on the structure. They must first be tested in computer simulations. Until now, the usual procedure for simulating controlled Large Space Structures is to compute a reduced order modal representation of the structure and then apply the controller. However, this procedure entails modal truncation error. A new software package which is free from this error is currently under development within the Center for Space Construction. The more accurate finite element representation of the structure is used in the simulation, instead of the less accurate reduced order modal representation. This software also features an efficient matrix storage scheme, which effectively deals with the asymmetric system matrices which occur when control is added to the structure. Also, an integration algorithm was chosen so that the simulation is a reliable indicator of system stability or instability. The software package is fairly general in nature. Linearity of the finite element model and of the controller is the only assumption made. Actuator dynamics, sensor dynamics, noise, and disturbances can be handled by the package. In addition, output feedback of displacement, velocity, and/or acceleration signals can be simulated. Kalman state estimation was also implemented. This software was tested on a finite element model of a real Large Space Structure: The Mini-Mast Truss. Mini-Mast is a testbed at NASA-Langley which is currently under development. A 714 degree of freedom finite element model was computed, and a 19 state controller was designed for it. Torque wheel dynamics were added to the model, and the entire closed loop system was simulated with the software package.

Quan, Ralph↗
DOE OSTI2019

Thermal and Thermomechanical Modeling to Design a Gallium Oxide Power Electronics Package: Preprint

The power electronics industry is currently undergoing a major transition by replacing silicon devices with wide-bandgap devices. Gallium oxide devices have the potential to accelerate this transition by offering comparable or even superior performance to other wide-bandgap devices, but at a much lower cost. Recent breakthroughs include demonstration of a laboratory-scale gallium oxide transistors and diodes; however, a functional power electronics package for these devices is yet to be developed. In this paper, the research methodology in designing an electronics package for gallium oxide devices is outlined. Finite element-based thermal and thermomechanical modeling simulations were conducted to realize a package design that meets the combined target of minimal thermal resistance and improved reliability. Different package designs that include various material combinations and cooling configurations were explored, and their thermal and thermomechanical performance are reported. Furthermore, the short-circuit withstanding capabilities of gallium oxide devices were studied and compared with silicon carbide.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

Damage development in neutron-irradiated concrete in a test reactor: Hygro-thermal and mechanical simulations

This paper reports the development of a 3D mesoscale hygro-thermal-mechanical simulation approach to predict damage in concrete irradiated in a test reactor. This framework, developed in MOOSE, considers the effects of elevated temperature, moisture content, and high neutron fluence (energy threshold, E > 0.1 MeV) on the mortar and aggregates separately. The first-stage simulation implements hygro-thermal analysis to determine the temperature and RH inside the specimen as a function of imposed radiation energy. These are used as inputs to the second stage, which considers radiation-induced volumetric expansion (RIVE) of aggregates, and creep, shrinkage, and stress-strain response of mortar to predict the expansion, stresses, and damage in specimens made using different coarse aggregates and subjected to different irradiation times. The irradiation time-dependent damage in the mortar is expressed using an isotropic damage parameter. This multi-physics model serves as a predictive tool for damage quantification in concrete due to neutron irradiation.

36 MATERIALS SCIENCE↗
NASA NTRS2009

Aeroelastic-Acoustics Simulation of Flight Systems

This paper describes the details of a numerical finite element (FE) based analysis procedure and a resulting code for the simulation of the acoustics phenomenon arising from aeroelastic interactions. Both CFD and structural simulations are based on FE discretization employing unstructured grids. The sound pressure level (SPL) on structural surfaces is calculated from the root mean square (RMS) of the unsteady pressure and the acoustic wave frequencies are computed from a fast Fourier transform (FFT) of the unsteady pressure distribution as a function of time. The resulting tool proves to be unique as it is designed to analyze complex practical problems, involving large scale computations, in a routine fashion.

Gupta, kajal K.↗
DOE OSTI2022

An E and B gyrokinetic simulation model for kinetic Alfvén waves in tokamak plasmas

The gyrokinetic particle simulation is a powerful tool for studies of transport, nonlinear phenomenon, and energetic particle physics in tokamak plasmas. While most gyrokinetic simulations make use of the scalar and vector potentials, a new model (GK-E&B) has been developed by using the E and B field in a general form and has been implemented in simulating kinetic Alfvén waves in uniform plasma. In our work, the Chen et al. GK-E&B model has been expressed, in general, tokamak geometry using the local orthogonal coordinates and general tokamak coordinates. Its reduction for uniform plasma is verified, and the numerical results show good agreement with the original work. The theoretical dispersion relation and numerical results in the local model in screw pinch geometry are also in excellent agreement. Numerical results show excellent performance in a realistic parameter regime of burning plasmas with high values of β/(M e k$^{2}_{⊥}$ρ$^{2}_{i}$), which is a challenge for traditional methods due to the “cancellation” problem. As one application, the GK-E&B model is implemented with kinetic electrons in the local single flux surface limit. With the matched International Tokamak Physics Activity-Toroidicity-induced Alfvén Eigenmodes parameters adopted, numerical results show the capability of the GK-E&B in treating the parallel electron Landau damping for realistic tokamak plasma parameters. As another application, the global GK-E&B model has been implemented with the dominant electron contribution in the cold electron limit. Its capability in simulating the finite E || due to the finite electron mass is demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Development of Predictive Model for Accurate Rupture Time from Multi-Axial Creep in Alloy 709 with Physics-Based Simulations

A physics-based model is developed to predict multiaxial creep behavior in Alloy 709 (A709), an advanced austenitic stainless steel intended for high-temperature applications such as Sodium Fast Reactors (SFRs). Compared to conventional stainless steels like 316H, A709 offers superior high-temperature performance; however, comprehensive data on its multiaxial creep response remain limited. To address this gap, a crystal plasticity finite element (CPFE) framework is used to simulate the deformation and failure mechanisms of A709 under multiaxial loading conditions. The model incorporates an extended Hu-Cocks dislocation creep formulation that accounts for precipitation effects, along with the Sham–Needleman model to capture grain boundary cavitation-driven failure. These advanced constitutive models enable a detailed understanding of the interplay between microstructural evolution and macroscopic creep response. Furthermore, the study evaluates the predictive accuracy of various effective stress measures in estimating creep rupture life, leveraging simulated multiaxial creep data. The findings provide critical insights into the applicability of different stress measures for engineering design and life prediction of A709 components operating under complex loading conditions. This work contributes to improving the reliability of high-temperature structural components by advancing predictive modeling capabilities for advanced austenitic steels.

Alloy 709↗
NASA NTRSDate not supplied

Simulating High Energy Dynamic Impact of IM7/PEKK Continuous Fiber Laminated Thermoplastic Composite using Open Hole Coupon Experiments

The use of advanced thermoplastic composites (TPC) for structural applications in the aerospace and automotive industries has grown increasingly popular due to their high performance, high manufactured part output, and sustainability. Enhancing verification simulations of the progressive failure in TPCs at the coupon scale is needed to better validate advanced composite material models at larger length scales, such as element level panels. MAT_213, a tabulated composite material model, has been applied to simulate high velocity dynamic impact (HEDI) and is integrated in the LS-DYNA explicit finite element (FE) software. Previous efforts using MAT_213 have calibrated HEDI simulated failures with experimental data. Traditionally simulations have relied on a structured mesh to represent quasi-isotropic panels. For this study, a continuous fiber unidirectional tape, IM7/PEKK, was analyzed through a series of verification studies of notched laminate coupons, then validation analyses of quasi-isotropic panels under HEDI were compared with experimental data. Experiments using unidirectional TPC were performed to obtain tabulated stress vs strain curves using Digital Image Correlation (DIC) for input into MAT_213. Next, notched coupon experiments were used as an intermediate step to calibrate the damage and failure behavior in MAT_213. Lastly, investigations into differences between HEDI simulations using a fiber aligned mesh were compared with experimental failure patterns and damage sizes. The use of experimentally generated material behavior along with selective mesh layout improved HEDI predictions with minimal material parameter calibration.

Polymer matrix composites↗
DOE OSTI2012

A combined three-dimensional finite element and scattering matrix method for the analysis of plane wave diffraction by bi-periodic, multilayered structures

A three-dimensional finite element method (FEM) for the analysis of plane wave diffraction by a bi-periodic slab is described and implemented. A scattering matrix formalism based on the FEM allows the efficient treatment of light reflection and transmission by multilayer bi-periodic structures, and the computation of Bloch modes of three-dimensional arrays. Numerical simulations, which show the accuracy and flexibility of the FEM, are presented.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Monitoring of mechanical changes in a pipe assembly with complex geometry using multi-mode acoustic signals

Guided wave ultrasonic testing has been widely used for nondestructive and long-range evaluation of diverse mechanical structures, such as pipes and plates. In general, a large number of guided wave modes exist even in the simplest geometries, and they exhibit highly dispersive characteristics. Thus, sophisticated excitation methods and elaborate modal analyses need to be implemented to excite, acquire, and process well-defined pure modes as part of such a measurement method. If the interrogated structure has complex geometry, for example, with cross-sectional variations or multiple joints, the utilization of pure modes becomes significantly more difficult. To overcome this difficulty, this paper proposes the use of multi-mode, multi-frequency guided waves for long-range inspection of structures with arbitrarily complex geometries. Furthermore, using broadband chirp signals and multi-mode time–frequency analysis in a prototypical complex pipe assembly that includes multiple pipe diameters, elbows, and flange connections, we experimentally show sensitive detection of mechanical changes (both material addition and loss). In addition, we demonstrate that the complexity of the structure can be modeled in the digital domain with sufficient detail to allow accurate numerical simulation of the acoustic response using this method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Temperature compensated transfer function for probability of detection (POD) in high temperature ultrasonic NDE using low temperature signals - 249

Ultrasonic piezoelectric based transducers offer the potential for in-coolant NDT measurements and imaging required at high temperature for the safe operation of advanced small modular reactors (SMRs). Current work estimates a significant 97% reduction in the probability of detection (POD) for an artificial flaw size of 2{lambda} due to degradation of the transducer performance at elevated temperature. The source of variability for POD in the current work is the temperature dependence of material matrix of piezoelectric material PZT-5A used in the sensor. The physics based model developed using finite element method quantifies reduction in the scattering amplitude for the side drilled hole (SDH) from 15 C to 195 C. Using model assisted probability of detection (MAPOD) approach, sensor temperature effect on the POD is demonstrated for SDHs in the low carbon steel specimen. Scattering amplitudes for the room temperature simulated data is validated with the experimental data. A temperature correction factor and transfer factors are reported to transform the simulated temperature effect from physics based model to the experimental data for estimation of high temperature POD using low temperature ultrasonic measurements. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2015

Impact Testing and Simulation of a Sinusoid Foam Sandwich Energy Absorber

A sinusoidal-shaped foam sandwich energy absorber was developed and evaluated at NASA Langley Research Center through multi-level testing and simulation performed under the Transport Rotorcraft Airframe Crash Testbed (TRACT) research project. The energy absorber, designated the "sinusoid," consisted of hybrid carbon- Kevlar® plain weave fabric face sheets, two layers for each face sheet oriented at +/-45deg with respect to the vertical or crush direction, and a closed-cell ELFOAM(TradeMark) P200 polyisocyanurate (2.0-lb/ft3) foam core. The design goal for the energy absorber was to achieve an average floor-level acceleration of between 25- and 40-g during the full-scale crash test of a retrofitted CH-46E helicopter airframe, designated TRACT 2. Variations in the design were assessed through quasi-static and dynamic crush testing of component specimens. Once the design was finalized, a 5-ft-long subfloor beam was fabricated and retrofitted into a barrel section of a CH-46E helicopter. A vertical drop test of the barrel section was conducted onto concrete to evaluate the performance of the energy absorber prior to retrofit into TRACT 2. Finite element models were developed of all test articles and simulations were performed using LSDYNA ®, a commercial nonlinear explicit transient dynamic finite element code. Test analysis results are presented for the sinusoid foam sandwich energy absorber as comparisons of load-displacement and acceleration-time-history responses, as well as predicted and experimental structural deformations and progressive damage for each evaluation level (component testing through barrel section drop testing).

Jackson, Karen E.↗
DOE OSTI2013

Preserving spherical symmetry in axisymmetric coordinates for diffusion problems

Persevering symmetric solutions, even in the under-converged limit, is important to the robustness of production simulation codes. We explore the symmetry preservation in both a continuous nodal and a mixed finite element method. In their standard formulation, neither method preserves spherical solution symmetry in axisymmetric (RZ) coordinates. We propose two methods, one for each family of finite elements, that recover spherical symmetry for low-order finite elements on linear or curvilinear meshes. This is a first step toward understanding achieving symmetry for higher-order elements. (authors)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2017

Constitutive Modeling of the Facesheet to Core Interface in Honeycomb Sandwich Panels Subject to Mode I Delamination

A new cohesive zone traction-separation law, which includes the effects of fiber bridging, has been developed, implemented with a finite element (FE) model, and applied to simulate the delamination between the facesheet and core of a composite honeycomb sandwich panel. The proposed traction-separation law includes a standard initial cohesive component, which accounts for the initial interfacial stiffness and energy release rate, along with a new component to account for the fiber bridging contribution to the delamination process. Single cantilever beam tests on aluminum honeycomb sandwich panels with carbon fiber reinforced polymer facesheets were used to characterize and evaluate the new formulation and its finite element implementation. These tests, designed to evaluate the mode I toughness of the facesheet to core interface, exhibited significant fiber bridging and large crack process zones, giving rise to a concave downward concave upward pre-peak shape in the load-displacement curve. Unlike standard cohesive formulations, the proposed formulation captures this observed shape, and its results have been shown to be in excellent quantitative agreement with experimental load-displacement and apparent critical energy release rate results, representative of a payload fairing structure, as well as local strain fields measured with digital image correlation.

composites↗