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,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

High-speed nanoscale characterization of dewetting via dynamic transmission electron microscopy

The dewetting of thin films can occur in either the solid or the liquid state for which different mass transport mechanisms are expected to control morphological changes. Traditionally, dewetting dynamics have been examined on time scales between several seconds to hours, and length scales ranging between nanometers and millimeters. The determination of mass transport mechanisms on the nanoscale, however, requires nanoscale spatial resolution and much shorter time scales. This study reports the high-speed observation of dewetting phenomena for kinetically constrained Ni thin films on crystalline SrTiO{sub 3} substrates. Movie-mode Dynamic Transmission Electron Microscopy (DTEM) was used for high-speed image acquisition during thin film dewetting at different temperatures. DTEM imaging confirmed that the initial stages of film agglomeration include edge retraction, hole formation, and growth. Finite element modeling was used to simulate temperature distributions within the DTEM samples after laser irradiation with different energies. For pulsed laser irradiation at 18 μJ, experimentally observed hole growth suggests that Marangoni flow dominates hole formation in the liquid nickel film. After irradiation with 13.8 μJ, however, the observations suggest that dewetting was initiated by nucleation of voids followed by hole growth through solid-state surface diffusion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Development of Graphite Thermal and Mechanical Modeling Capabilities in Grizzly

Nuclear-grade graphites are used extensively in the core designs of multiple types of advanced nuclear reactors. In the reactor environment, graphite is exposed over long durations to extreme conditions, including high temperatures, radiation and potentially molten salt and oxygen. Exposure to these conditions can cause several degradation mechanisms in graphite, including nonuniform volumetric strains induced by irradiation and thermal expansion, which lead to stresses that can compromise the performance of graphite components. Evaluating component integrity, predicting component performance over the reactor lifetime, and developing design standards all require robust tools for predicting fracture initiation and propagation in graphite structural components in nuclear reactors. This report documents progress in an ongoing effort to develop modeling and simulation tools in the Grizzly code for predicting the performance of graphite exposed to reactor conditions. Recent developments include a set of thermal and mechanical models that now include the IG-110, NBG-18, and H-451 graphite grades. Improvements have also been made to a nonlinear damaged plasticity model applicable to predicting damage under tension and compression to quasibrittle materials, including graphite. In addition, enhancements have been made to the extended finite element method implementation targeted at simulating graphite fracture. These include new capabilities for crack nucleation in the interior of a solid body, improved treatment of crack nucleation on free surfaces, and more robust modeling approaches for crack growth approaching free surfaces or other cracks.

36 MATERIALS SCIENCE↗
DOE OSTI2023

TReactMech v4.217

TReactMech couples geomechanical processes (poroelasticity, failure, and inelastic strain) with multiphase nonisothermal flow (derived from TOUGH2) and reactive geochemical transport. At its core is the reactive-transport code TOUGHREACT v4.13. TReactMech is an efficient hybrid parallel simulator, solving the geomechanics using finite elements and MPI/PetSc, the multiphase flow using integrated finite difference and MPI/PETSc, and the reactive chemistry using OpenMP. The advantages of TReactMech are in its multiphase flow capabilities (e.g., supercritical CO2, supercritical water, air) and parallel geomechanics including full 3-D stress tensor, shear and tensile failure, coupled to porosity and permeability changes. It is backwardly compatible with TOUGH2 and TOUGHREACT v4.13, allowing for easier transitions between the codes. TReactMech can be used to simulate many natural and engineered subsurface systems, including geothermal reservoirs, borehole heat exchangers, geologic carbon sequestration, geologic storage of nuclear waste, groundwater resources, weathering, sediment diagenesis, seafloor hydrothermal circulation, hydrofracturing in unconventional reservoirs, and injection/production-induced surface deformation.

Sonnenthal, Eric↗
DOE OSTI2015

Seismic design of a nuclear waste repository against an earthquake induced rock shear load - 15346

Concepts for isolation of spent nuclear fuel have been studied by organizations in countries that employ nuclear power. These concepts typically involve a deep geological repository in bedrock, with processed radioactive waste encased in deposition canisters emplaced in underground excavations and surrounded by clay-based buffer. This multi-barrier approach is intended to provide long-term isolation of the waste package from the biosphere. According to many statistics, there is a tendency toward increasing numbers of earthquakes in the Korean peninsula each year. If an earthquake induced rock shear load occurs at a repository site, the engineered barrier system (EBS) can be sheared, resulting in a pressure increased in the buffer material. This is expected to result in damage to the canisters. Therefore, because the damage, which can be fatal to a disposal canister, can be given in case an earthquake occurs, the method protecting the canisters safely is needed. In this paper, 'an earthquake-resistance type buffer' was proposed using a method for protection from earthquake damage. To aid in achieving this objective, we have developed a technique for estimating the behavior of the EBS affected by earthquake induced rock shear using an analysis method, such as finite element method. As a first step, experimental tests are performed for understanding the behavior of the buffer material. Comparing the results of the numerical simulations to the experimental results can be used to verity the reliability of the analysis model. Then we developed the real scale analysis model of a nuclear waste repository. The main parameter having an effect on the earthquake-resistant performance was analyzed and an earthquake-resistance proof type buffer was designed. To improve the earthquake-resistance performance, we designed the buffer by putting the low- density buffer between the canister and buffer. The simulation results show that the earthquake-resistant performance was improved by about 65%. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1991

A combined finite element-boundary element formulation for solution of axially symmetric bodies

A new method is presented for the computation of electromagnetic scattering from axially symmetric bodies. To allow the simulation of inhomogeneous cross sections, the method combines the finite element and boundary element techniques. Interior to a fictitious surface enclosing the scattering body, the finite element method is used which results in a sparce submatrix, whereas along the enclosure the Stratton-Chu integral equation is enforced. By choosing the fictitious enclosure to be a right circular cylinder, most of the resulting boundary integrals are convolutional and may therefore be evaluated via the FFT with which the system is iteratively solved. In view of the sparce matrix associated with the interior fields, this reduces the storage requirement of the entire system to O(N) making the method attractive for large scale computations. The details of the corresponding formulation and its numerical implementation are described.

Collins, Jeffrey D.↗
DOE OSTI2025

Finite Element Analysis of the TRUST Nonlinear Dynamics Testbed

This paper builds on prior work conducted within the Los Alamos National Laboratory (LANL) Testbeds to Reduce Uncertainty in Simulations and Tests (TRUST) program. Specifically, it builds on finite element (FE) modeling efforts for the TRUST program’s Nonlinear Dynamics (ND) testbed. Historically, the FE model for the ND testbed has exclusively utilized a Lanczos eigensolver that linearly extracts the system’s natural frequencies. This paper investigates the Abaqus 2024’s explicit dynamic solver, which implements a central difference explicit solver. The central difference method used in Abaqus can capture nonlinear material responses in structural dynamic simulations making it suitable for the ND FE model.

42 ENGINEERING↗
NASA NTRS2002

P1 Nonconforming Finite Element Method for the Solution of Radiation Transport Problems

The simulation of radiation transport in the optically thick flux-limited diffusion regime has been identified as one of the most time-consuming tasks within large simulation codes. Due to multimaterial complex geometry, the radiation transport system must often be solved on unstructured grids. In this paper, we investigate the behavior and the benefits of the unstructured P(sub 1) nonconforming finite element method, which has proven to be flexible and effective on related transport problems, in solving unsteady implicit nonlinear radiation diffusion problems using Newton and Picard linearization methods. Key words. nonconforrning finite elements, radiation transport, inexact Newton linearization, multigrid preconditioning

Kang, Kab S.↗
NASA NTRS1988

An internal variable constitutive model for the large deformation of metals at high temperatures

The advent of large deformation finite element methodologies is beginning to permit the numerical simulation of hot working processes whose design until recently has been based on prior industrial experience. Proper application of such finite element techniques requires realistic constitutive equations which more accurately model material behavior during hot working. A simple constitutive model for hot working is the single scalar internal variable model for isotropic thermal elastoplasticity proposed by Anand. The model is recalled and the specific scalar functions, for the equivalent plastic strain rate and the evolution equation for the internal variable, presented are slight modifications of those proposed by Anand. The modified functions are better able to represent high temperature material behavior. The monotonic constant true strain rate and strain rate jump compression experiments on a 2 percent silicon iron is briefly described. The model is implemented in the general purpose finite element program ABAQUS.

Brown, Stuart↗
DOE OSTI2015

A stochastic approach to uncertainty in the equations of MHD kinematics

The magnetohydrodynamic (MHD) kinematics model describes the electromagnetic behavior of an electrically conducting fluid when its hydrodynamic properties are assumed to be known. In particular, the MHD kinematics equations can be used to simulate the magnetic field induced by a given velocity field. While prescribing the velocity field leads to a simpler model than the fully coupled MHD system, this may introduce some epistemic uncertainty into the model. If the velocity of a physical system is not known with certainty, the magnetic field obtained from the model may not be reflective of the magnetic field seen in experiments. Additionally, uncertainty in physical parameters such as the magnetic resistivity may affect the reliability of predictions obtained from this model. By modeling the velocity and the resistivity as random variables in the MHD kinematics model, we seek to quantify the effects of uncertainty in these fields on the induced magnetic field. We develop stochastic expressions for these quantities and investigate their impact within a finite element discretization of the kinematics equations. We obtain mean and variance data through Monte Carlo simulation for several test problems. Toward this end, we develop and test an efficient block preconditioner for the linear systems arising from the discretized equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Finite Element Analysis and Machine Learning Guided Design of Carbon Fiber Organosheet-Based Battery Enclosures for Crashworthiness

Carbon fiber composite can be a potential candidate for replacing metal-based battery enclosures of current electric vehicles (E.V.s) owing to its better strength-to-weight ratio and corrosion resistance. However, the strength of carbon fiber-based structures depends on several parameters that should be carefully chosen. Here, in this work, we implemented high throughput finite element analysis (FEA) based thermoforming simulation to virtually manufacture the battery enclosure using different design and processing parameters. Subsequently, we performed virtual crash simulations to mimic a side pole crash to evaluate the crashworthiness of the battery enclosures. This high throughput crash simulation dataset was utilized to build predictive models to understand the crashworthiness of an unknown set. Our machine learning (ML) models showed excellent performance (R 2 > 0.97) in predicting the crashworthiness metrics, i.e., crush load efficiency, absorbed energy, intrusion, and maximum deceleration during a crash. We believe that this FEA-ML work framework will be helpful in down select process parameters for carbon fiber-based component design and can be transferrable to other manufacturing technologies.

36 MATERIALS SCIENCE↗
NASA NTRS2004

Multi-Terrain Impact Testing and Simulation of a Composite Energy Absorbing Fuselage Section

Comparisons of the impact performance of a 5-ft diameter crashworthy composite fuselage section were investigated for hard surface, soft soil, and water impacts. The fuselage concept, which was originally designed for impacts onto a hard surface only, consisted of a stiff upper cabin, load bearing floor, and an energy absorbing subfloor. Vertical drop tests were performed at 25-ft/s onto concrete, soft-soil, and water at NASA Langley Research Center. Comparisons of the peak acceleration values, pulse durations, and onset rates were evaluated for each test at specific locations on the fuselage. In addition to comparisons of the experimental results, dynamic finite element models were developed to simulate each impact condition. Once validated, these models can be used to evaluate the dynamic behavior of subfloor components for improved crash protection for hard surface, soft soil, and water impacts.

Fasanella, Edwin L.↗
NASA NTRS2007

Multi-Terrain Impact Testing and Simulation of a Composite Energy Absorbing Fuselage Section

Comparisons of the impact performance of a 5-ft diameter crashworthy composite fuselage section were investigated for hard surface, soft soil, and water impacts. The fuselage concept, which was originally designed for impacts onto a hard surface only, consisted of a stiff upper cabin, load bearing floor, and an energy absorbing subfloor. Vertical drop tests were performed at 25-ft/s onto concrete, soft-soil, and water at NASA Langley Research Center. Comparisons of the peak acceleration values, pulse durations, and onset rates were evaluated for each test at specific locations on the fuselage. In addition to comparisons of the experimental results, dynamic finite element models were developed to simulate each impact condition. Once validated, these models can be used to evaluate the dynamic behavior of subfloor components for improved crash protection for hard surface, soft soil, and water impacts.

Fasanella, Edwin L.↗
NASA NTRSDate not supplied

Applicability of Loads Estimation Techniques Using Sparse Acceleration Sensor Data to Spacecraft Structural Health Monitoring

The use of structural health monitoring systems on spacecraft structures can play a crucial role in ensuring the safety, reliability, and longevity of the structure by gathering and analyzing onboard sensor data. Of specific importance is monitoring for excessive loading at critical interfaces as any off-nominal structural excitations experienced by spacecraft structures can cause early unpredicted high structural life consumption or damage. The availability and cost of flight-certified sensors along with the size of spacecraft structures and allowable payload mass drives the need for a method to estimate loads using sparsely-located sensors. Numerous approaches such as physics-based, statistical learning, and physics-enhanced statistical learning algorithms have gained popularity among structural prognostics applications. However, developing noise-robust prediction models to assess loads and structural life predictions from a sparse multi-sensor data acquisition system can be a challenging task. This paper discusses the evaluation of physics-based versus machine-learning algorithms for predicting loads and structural life at mission critical locations on the spacecraft structure using a finite element loads analysis with the application of simulated noise and noise reduction techniques. To estimate the loads from accelerations, the physics-based algorithm leverages a loads transformation matrix from a Craig-Bampton reduced finite element model. A System Equivalent Reduction Expansion Process (SEREP) and a pseudo-inverse approach are considered to expand from the onboard sensor degrees of freedom to the Craig-Bampton model degrees of freedom. The machine learning algorithm provides a data driven solution/mapping of the sensor accelerations to the loads at the mission critical locations using a high dimensionality analysis. Although these strategies produce comparable loads prediction without noise, the limitations of these strategies with incorporating simulated noise and noise reduction techniques with low signal to noise ratio signals are evaluated. The study demonstrates the immense potential of statistical learning algorithms for sparse structural prognostic models and enhancing signal denoising techniques. These findings also highlight the need for noise-resilient prognostic models and low-noise data acquisition systems onboard spacecraft structures.

Spacecraft Structural Health Monitoring↗
DOE OSTI2018

Influences of Quench Cooling Rate on Microstructure and Corrosion Resistance of Al-Cu-Mg Alloy Based on the End-Quenching Test

To investigate the effects of the quench cooling rate on corrosion resistance of Al-Cu-Mg alloy, an end-quenching test was conducted and the microstructures at different cooling rates were observed by SEM and TEM. Additionally, the corrosion resistance was characterized by an intergranular corrosion test and electrochemical test. Moreover, the finite element method was applied to simulate the end quenching process. The results indicate that the actual end quenching process can be approximated as one-dimensional heat transfer, and the cooling rate varies at different cooling distances. By affecting the microstructures, decreasing the cooling rate leads to a decline in the corrosion properties. Low cooling rates coarsen the constituent particles and grain boundary particles, resulting in a wide precipitation-free zone and an increase in the intensity of corrosion reactions. A high cooling rate concentrates on the intragranular precipitant, which can reduce the pitting depth and represents a conversion from localized corrosion to general corrosion.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Compact spherical neutron polarimeter using high-T{sub c} YBCO films

We describe a simple, compact device for spherical neutron polarimetry measurements at small neutron scattering angles. The device consists of a sample chamber with very low (<0.01 G) magnetic field flanked by regions within which the neutron polarization can be manipulated in a controlled manner. This allows any selected initial and final polarization direction of the neutrons to be obtained. We have constructed a prototype device using high-T{sub c} superconducting films and mu-metal to isolate regions with different magnetic fields and tested device performance in transmission geometry. Finite-element methods were used to simulate the device’s field profile and these have been verified by experiment using a small solenoid as a test sample. Measurements are reported using both monochromatic and polychromatic neutron sources. The results show that the device is capable of extracting sample information and distinguishing small angular variations of the sample magnetic field. As a more realistic test, we present results on the characterization of a 10 μm thick Permalloy film in zero magnetic field, as well as its response to an external magnetic field.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1994

Analysis of new composite architectures

Global/local method in conjunction with a special macro finite element is used to gain computational efficiency in the simulation of textile composites behavior. Results are included to demonstrate the effectiveness of the method. Also, 2D finite elements are used to study boundary effects in plain weave composite specimens subjected to extension, shear, and flexure loads. Effective extension, shear, and flexural moduli were found to be quite sensitive to specimen size. For extension and flexure loads stress distributions were affected by a free surface, but the free surface boundary effect did not appear to propagate very far into the interior. For shear load, the boundary effect appeared to propagate much further into the interior. The report is in the form of two technical articles, the first describes the global/local method while the second describes the use of the 2D finite elements.

Whitcomb, John D.↗
DOE OSTI2016

Incorporating physically-based microstructures in materials modeling: Bridging phase field and crystal plasticity frameworks

Here, the mechanical properties of materials systems are highly influenced by various features at the microstructural level. The ability to capture these heterogeneities and incorporate them into continuum-scale frameworks of the deformation behavior is considered a key step in the development of complex non-local models of failure. In this study, we present a modeling framework that incorporates physically-based realizations of polycrystalline aggregates from a phase field (PF) model into a crystal plasticity finite element (CP-FE) framework. Simulated annealing via the PF model yields ensembles of materials microstructures with various grain sizes and shapes. With the aid of a novel FE meshing technique, FE discretizations of these microstructures are generated, where several key features, such as conformity to interfaces, and triple junction angles, are preserved. The discretizations are then used in the CP-FE framework to simulate the mechanical response of polycrystalline α-iron. It is shown that the conformal discretization across interfaces reduces artificial stress localization commonly observed in non-conformal FE discretizations. The work presented herein is a first step towards incorporating physically-based microstructures in lieu of the overly simplified representations that are commonly used. In broader terms, the proposed framework provides future avenues to explore bridging models of materials processes, e.g. additive manufacturing and microstructure evolution of multi-phase multi-component systems, into continuum-scale frameworks of the mechanical properties.

36 MATERIALS SCIENCE↗