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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

The evolution history of superalloy powders during hot consolidation and plastic deformation

Highlights: • Element segregations in a powder metallurgy superalloy at different states are identified. • The surface energy drives precipitation of oxides, carbides, and γ' particles on powder surface during hot consolidation. • Prior particle boundary structure is distorted and more uniformly refined grains are obtained after hot extrusion. - Abstract: Powder metallurgy superalloys are widely used to produce critical rotating parts in turbine engines. However, the issues like what is the driving force of precipitation on prior particle boundary (PPB) during consolidation and how the PPB structure evolves during plastic deformation still remain to be resolved. In this study, with the characterization of the microstructure evolution from powder, consolidation, to plastic deformation, the segregation of alloying elements on the particle boundaries at elevated temperature was analyzed by equilibrium thermodynamic calculations, and the distortion of PPB structure under large plastic deformation was quantitively simulated by finite element (FE) methods. The results indicated that surface energy serves as the main driving force of PPB precipitation and hot extrusion could distort the PPB structure and yield more uniformly refined grains.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Simulations of Two Heavily Confined PBX 9502 Cook-off Experiments.

Finite element models of two cook-off experiments were constructed using the Aria code.These experiments provided thermocouple measurements that were used as boundary conditions to carry out finite element (FE) simulations and comparisons of the FE thermal response with internal temperature measurements in the PBX 9502. These FE models contain a pressure dependent PBX 9502 thermal decomposition model developed by Hobbs’ et al.. Comparisons of the FE results and the experiments are presented. These comparisons assess the predictive capabilities of Hobbs’ model under the experimental conditions. Limiting cases of the PBX 9502 decomposition model are considered. Two parameterizations associated with the pressure dependency of the PBX 9502 decomposition model were carried out to fit the experiments using internal temperature measurements. Two inverse heat conduction problems were implemented to estimate temperatures as functions of time at locations on the outer surface of the PBX 9502. Thermocouple conduction error studies were carried out to assess the comparisons of the FE model thermal response with the experiments.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS2004

Stress Distribution During Deformation of Polycrystalline Aluminum by Molecular-Dynamics and Finite-Element Modeling

In this paper, a multiscale modelling strategy is used to study the effect of grain-boundary sliding on stress localization in a polycrystalline microstructure with an uneven distribution of grain size. The development of the molecular dynamics (MD) analysis used to interrogate idealized grain microstructures with various types of grain boundaries and the multiscale modelling strategies for modelling large systems of grains is discussed. Both molecular-dynamics and finite-element (FE) simulations for idealized polycrystalline models of identical geometry are presented with the purpose of demonstrating the effectiveness of the adapted finite-element method using cohesive zone models to reproduce grain-boundary sliding and its effect on the stress distribution in a polycrystalline metal. The yield properties of the grain-boundary interface, used in the FE simulations, are extracted from a MD simulation on a bicrystal. The models allow for the study of the load transfer between adjacent grains of very different size through grain-boundary sliding during deformation. A large-scale FE simulation of 100 grains of a typical microstructure is then presented to reveal that the stress distribution due to grain-boundary sliding during uniform tensile strain can lead to stress localization of two to three times the background stress, thus suggesting a significant effect on the failure properties of the metal.

Yamakov, V.↗
NASA NTRS2013

FE Simulation of SMA Seal for Mars Sample Return

Several NASA rovers and landers have been on Mars and performed successful in-situ exploration. Returning Martian samples to Earth for extensive analysis is of great interest to the planetary science community. Current Mars sample return architecture would require leaving the acquired samples on Mars for years before being retrieved by subsequent mission. Each sample would be sealed securely to keep its integrity. A reliable seal technique that does not affect the integrity of the samples and uses a simple low-mass tool is required. The shape memory alloy (SMA) seal technique is a promising candidate. A study of the thermal performances of several primary designs of a SMA seal for sample tubes by finite element (FE) simulation are presented in this paper. The results show sealing the sample tube by SMA plugs and controlling the sample temperature below the allowed temperature level are feasible.

seal↗
DOE OSTI2019

Unsupervised machine learning based on non-negative tensor factorization for analyzing reactive-mixing

Analysis of reactive-diffusion simulations representing complex mixing processes requires a large number of independent model runs. For each high-fidelity model simulation, the model inputs are varied and the predicted mixing behavior is represented by temporal and spatial changes in species concentration. It is then required to discern how the model inputs (such as diffusivity, dispersion, anisotropy, and velocity field properties) impact the mixing process. This task is challenging and typically involves interpretation of large model outputs representing temporal and spatial changes of species concentration within the model domain. Yet, the task can be automated and substantially simplified by applying Machine Learning (ML) methods. Here, we present an application of an unsupervised ML method (called NTF k ) using Non-negative Tensor Factorization (NTF) coupled with a custom clustering procedure based on k-means to reveal the temporal and spatial features in product concentrations. An attractive and unique aspect of the proposed ML method is that it ensures the extracted features are non-negative, which is important to obtain a meaningful deconstruction of the mixing processes. The ML methodology is applied to a large set of high-resolution finite-element model simulations representing anisotropic reaction-diffusion processes in perturbed vortex-based velocity fields. The applied finite-element method ensures that spatial and temporal species concentration are always non-negative, even in the case of high anisotropic contrasts. The simulated reaction is a fast irreversible bimolecular reaction A+B→C , where species A and B react to form species C. The reactive-diffusion model input parameters that control mixing include properties of the velocity field (such as vortex structures), anisotropic dispersion, and molecular diffusion. We demonstrate the applicability of the ML feature extraction method to produce a meaningful deconstruction of model outputs to discriminate between different physical processes impacting the reactants, their mixing, and the spatial distribution of the product C. The introduced ML analysis allowed us to identify additive temporal and spatial features that characterize mixing behavior. The application of the proposed NTF k approach is not limited to reactive-mixing. NTFk can be readily applied to any observed or simulated datasets that can be represented as tensors (multi-dimensional arrays) and have separable latent signatures or features.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2000

Comparison of Nonlinear Random Response Using Equivalent Linearization and Numerical Simulation

A recently developed finite-element-based equivalent linearization approach for the analysis of random vibrations of geometrically nonlinear multiple degree-of-freedom structures is validated. The validation is based on comparisons with results from a finite element based numerical simulation analysis using a numerical integration technique in physical coordinates. In particular, results for the case of a clamped-clamped beam are considered for an extensive load range to establish the limits of validity of the equivalent linearization approach.

Rizzi, Stephen A.↗
DOE OSTI2018

Room-Temperature Charpy Impact Property of 3D-Printed 15-5 Stainless Steel

In this study, the room-temperature Charpy impact property of 3D-printed 15-5 stainless steel was investigated by a combined experimental and finite element modeling approach. The experimentally measured impact energy is 10.85 ± 1.20 J/cm{sup 2}, which is comparable to the conventionally wrought and non-heat treated 15-5 stainless steel. In parallel to the impact test experiment, a finite element model using the Johnson–Cook material model with damage parameters was developed to simulate the impact test. The simulated impact energy is 10.46 J/cm{sup 2}, which is in good agreement with the experimental data. The fracture surface from the experimentally tested specimen suggests that the 3D-printed specimens undergo predominately brittle fracture.

36 MATERIALS SCIENCE↗
NASA NTRS1989

Optical laboratory solution and error model simulation of a linear time-varying finite element equation

The use of simplified error models to accurately simulate and evaluate the performance of an optical linear-algebra processor is described. The optical architecture used to perform banded matrix-vector products is reviewed, along with a linear dynamic finite-element case study. The laboratory hardware and ac-modulation technique used are presented. The individual processor error-source models and their simulator implementation are detailed. Several significant simplifications are introduced to ease the computational requirements and complexity of the simulations. The error models are verified with a laboratory implementation of the processor, and are used to evaluate its potential performance.

Taylor, B. K.↗
DOE OSTI2022

Thermal Design and Optimization of High- Power Wireless Charging System

In this paper, the thermal design and optimization of a high-power wireless charging system (WCS) is proposed. An integrated electromagnetic and thermal co-optimization is essential to design a high power-density WCS pad. This study presents the thermal analysis and the potential optimization scopes for a polyphase WCS pad. The coil and core causes most of the power losses of a WCS pad and causes thermal hotspot in the pad. In this paper, the thermally conductive epoxy is introduced to balance the effective volumetric loss density in the pad and mitigate the thermal hotspots in the coil and core. The proposed design is simulated through finite element analysis (FEA) and tested experimentally for a 50 kW three-phase WCS. Simulation and experimental results show that the conventional design has localized hotspots in the coil and core. The proposed thermal design mitigates the thermal hotspots without increasing the coil or core volume.

Mohammad, Mostak↗
DOE OSTI2018

Fourier amplitude sensitivity test applied to dynamic combined finite-discrete element methods-based simulations

Fracture propagation plays a key role for a number of applications of interest to the scientific community, from dynamic fracture processes like spallation and fragmentation in metals to failure of ceramics, airplane wings, etc. Simulations of material deformation and fracture propagation rely on accurate knowledge of material characteristics such as material strength and the amount of energy being dissipated during the fracture process. Within the combined finite-discrete element method (FDEM) framework material fracture behavior is typically described through a parametrized softening curve, which defines a stress-strain relationship unique to each material. We apply the Fourier amplitude sensitivity test to explore how each of these parameters influences the simulated damage processes and to determine the key input parameters that have the most impact on the model response. We present several sensitivity numerical experiments for the simulation of a split Hopkinson pressure bar (SHPB) test for weathered granite samples using different combinations of model parameters. We validate the obtained results against SHPB experimental data. The experiments show that the model is mostly sensitive to parameters related to tensile and shear strengths, even in the presence of other parameter perturbations. Lastly, the results suggest that the specification of tensile and shear strengths at the interfaces dominate the stress-time history of the FDEM simulation of SHPB test.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Data-driven Whitney forms for structure-preserving control volume analysis

Control volume analysis models physics via the exchange of generalized fluxes between subdomains. Here, we introduce a scientific machine learning framework adopting a partition of unity architecture to identify physically-relevant control volumes, with generalized fluxes between subdomains encoded via Whitney forms. The approach provides a differentiable parameterization of geometry which may be trained in an end-to-end fashion to extract reduced models from full field data while exactly preserving physics. The architecture admits a data-driven finite element exterior calculus allowing discovery of mixed finite element spaces with closed form quadrature rules. An equivalence between Whitney forms and graph networks reveals that the geometric problem of control volume learning is equivalent to an unsupervised graph discovery problem. The framework is developed for manifolds in arbitrary dimension, with examples provided for H(div) problems in $\mathbb{R}$ establishing convergence and structure preservation properties. Finally, we consider a lithium-ion battery problem where we discover a reduced finite element space encoding transport pathways from high-fidelity microstructure resolved simulations. The approach reduces the 5.89M finite element simulation to 136 elements while reproducing pressure to under 0.1% error and preserving conservation.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1981

Engine dynamic analysis with general nonlinear finite-element codes. I Overall approach and development of bearing damper element

NASA-sponsored research on engine dynamic simulation using general finite element nonlinear time transient computer codes available on the open market is reviewed. The approach taken was to develop software packages to model engine components which are not typically found on dynamical structures and are therefore not already computer codes. The software package developed for squeeze-film bearing dampers is outlined, and the results of a parametric study of damper pressure for a variety of specified circular orbits are presented for both long-bearing and short-bearing solutions. The data from a four-degree-of-freedom rotor-damper-stator model under conditions of small rotor unbalance through large rotor unbalance are also given.

Adams, M. L.↗
DOE OSTI2019

Modelling and simulation of a porous core photonic crystal fibre for terahertz wave propagation

A porous core photonic crystal fibre based on conventional hexagonal lattice cladding is proposed for propagating terahertz radiation. The structure is designed and theoretically investigated using full vectorial finite element method. Simulation results show that at 300 µm core diameter, with a high porosity of 85%, and an operating frequency of 1.3 THz, the proposed fibre reduces the bulk absorption loss of cyclic olefin copolymer (TOPAS) by about 81%, which corresponds to an ultra-low effective material loss value of 0.039 cm{sup −1}. Furthermore, the proposed fibre shows near zero dispersion coefficient of 0.47 ps/THz/cm with an extremely small variation of 0.05 over a broad 1.3 THz bandwidth; with confinement and bending losses investigated and found to be negligibly low. It is anticipated that the proposed waveguide can potentially be used for short range transmission of terahertz radiation in the communication window.

42 ENGINEERING↗
DOE OSTI2016

Simulation of CNT-AFM tip based on finite element analysis for targeted probe of the biological cell

Carbon nanotubes (CNTs) are potentially ideal tips for atomic force microscopy (AFM) due to the robust mechanical properties, nano scale diameter and also their ability to be functionalized by chemical and biological components at the tip ends. This contribution develops the idea of using CNTs as an AFM tip in computational analysis of the biological cell’s. Finite element analysis employed for each section and displacement of the nodes located in the contact area was monitored by using an output database (ODB). This reliable integration of CNT-AFM tip process provides a new class of high performance nanoprobes for single biological cell analysis.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
NASA NTRS2003

Residual Strength Analyses of Monolithic Structures

Finite-element fracture simulation methodology predicts the residual strength of damaged aircraft structures. The methodology uses the critical crack-tip-opening-angle (CTOA) fracture criterion to characterize the fracture behavior of the material. The CTOA fracture criterion assumes that stable crack growth occurs when the crack-tip angle reaches a constant critical value. The use of the CTOA criterion requires an elastic- plastic, finite-element analysis. The critical CTOA value is determined by simulating fracture behavior in laboratory specimens, such as a compact specimen, to obtain the angle that best fits the observed test behavior. The critical CTOA value appears to be independent of loading, crack length, and in-plane dimensions. However, it is a function of material thickness and local crack-front constraint. Modeling the local constraint requires either a three-dimensional analysis or a two-dimensional analysis with an approximation to account for the constraint effects. In recent times as the aircraft industry is leaning towards monolithic structures with the intention of reducing part count and manufacturing cost, there has been a consistent effort at NASA Langley to extend critical CTOA based numerical methodology in the analysis of integrally-stiffened panels.In this regard, a series of fracture tests were conducted on both flat and curved aluminum alloy integrally-stiffened panels. These flat panels were subjected to uniaxial tension and during the test, applied load-crack extension, out-of-plane displacements and local deformations around the crack tip region were measured. Compact and middle-crack tension specimens were tested to determine the critical angle (wc) using three-dimensional code (ZIP3D) and the plane-strain core height (hJ using two-dimensional code (STAGS). These values were then used in the STAGS analysis to predict the fracture behavior of the integrally-stiffened panels. The analyses modeled stable tearing, buckling, and crack branching at the integral stiffener using different values of critical CTOA for different material thicknesses and orientation. Comparisons were made between measured and predicted load-crack extension, out-of-plane displacements and local deformations around the crack tip region. Simultaneously, three-dimensional capabilities to model crack branching and to monitor stable crack growth of multiple cracks in a large thick integrally-stiffened flat panels were implemented in three-dimensional finite element code (ZIP3D) and tested by analyzing the integrally-stiffened panels tested at Alcoa. The residual strength of the panels predicted from STAGS and ZP3D code compared very well with experimental data. In recent times, STAGS software has been updated with new features and now one can have combinations of solid and shell elements in the residual strength analysis of integrally-stiffened panels.

Forth, Scott↗
DOE OSTI2023

On Practical Aspects of Variational Consistency in Contact Dynamics

Usage of contact mechanics methodologies is a pervasive modeling requirement in dynamic simulations. While for some trivial problems, solutions taken from analytical geometry are available, use of a finite element framework is common to achieve formulation generality. This work explores two dynamic contact formulations: one based on the traditional node-to-segment (NTS) approach, and a variationally consistent segment-to-segment (STS) mortar formulation. The NTS formulation employed here enforces the constraints kinematically (i.e., the interpenetration is enforced to the solver tolerance), whereas the mortar approach uses Lagrange multipliers to enforce the contact constraints. Both approaches are implemented in the open-source finite element framework Multiphysics Object-Oriented Simulation Environment (MOOSE). The results highlight two relevant contact-interface-related dynamic phenomena in finite element simulations. First, stabilization of contact constraints is discussed, taking into account the evolution of the total energy in a benchmark problem. Second, the influence of finite element discretization on both of the aforementioned contact formulations is analyzed by exercising a large-deformation example with continuous relative sliding. Variationally consistent contact approaches such as the mortar formulation lead to improved energy preservation and avoid spurious excitation of the system's frequencies. This is especially relevant in settings where inertia and vibrations are of importance.

42 ENGINEERING↗
NASA NTRS1995

Verification of a three-dimensional resin transfer molding process simulation model

Experimental evidence was obtained to complete the verification of the parameters needed for input to a three-dimensional finite element model simulating the resin flow and cure through an orthotropic fabric preform. The material characterizations completed include resin kinetics and viscosity models, as well as preform permeability and compaction models. The steady-state and advancing front permeability measurement methods are compared. The results indicate that both methods yield similar permeabilities for a plain weave, bi-axial fiberglass fabric. Also, a method to determine principal directions and permeabilities is discussed and results are shown for a multi-axial warp knit preform. The flow of resin through a blade-stiffened preform was modeled and experiments were completed to verify the results. The predicted inlet pressure was approximately 65% of the measured value. A parametric study was performed to explain differences in measured and predicted flow front advancement and inlet pressures. Furthermore, PR-500 epoxy resin/IM7 8HS carbon fabric flat panels were fabricated by the Resin Transfer Molding process. Tests were completed utilizing both perimeter injection and center-port injection as resin inlet boundary conditions. The mold was instrumented with FDEMS sensors, pressure transducers, and thermocouples to monitor the process conditions. Results include a comparison of predicted and measured inlet pressures and flow front position. For the perimeter injection case, the measured inlet pressure and flow front results compared well to the predicted results. The results of the center-port injection case showed that the predicted inlet pressure was approximately 50% of the measured inlet pressure. Also, measured flow front position data did not agree well with the predicted results. Possible reasons for error include fiber deformation at the resin inlet and a lag in FDEMS sensor wet-out due to low mold pressures.

Fingerson, John C.↗