Search NASA⌕ Search

SEARCH · Search NASA

Results for “material response 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,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Transient Thermal Analysis of Calorimeters Used in Characterization of the ACRR Radiation Environments

Silicon calorimeters have been used for active radiation dosimetry in the central cavity of the Annular Core Research Reactor (ACRR) for over a decade. Recently, there has been interest in using other materials for calorimetry to accurately measure the prompt gamma-ray energy deposition in the mixed neutron and gamma-ray environment. The calorimeters used in the ACRR use a thermocouple (TC) to measure the change in temperature of specific materials in the radiation environment. The temperature change is related to the instantaneous dose received by the material in a pulse-transient operation. SOLIDWORKS Simulation and ANSYS Mechanical were used to model the calorimeter and analyze the thermal behavior under pulse-transient conditions. This report compares the results from modeling to experimental results for selected calorimeter materials and radiation environments. These materials include bismuth, tin, zirconium, and silicon. Calorimeters assembled with each material were irradiated in the ACRR central cavity in the free- field, LB44, CdPoly, and PLG radiation environments. The neutronics code Monte-Carlo N- Particle (MCNP) was used to calculate the neutron and gamma-ray response of the calorimeter materials at the experimental locations in the central cavity. Different response tallies were used and found to give different results for the gamma-ray energy deposition. It was determined that performing the neutron/gamma-ray/electron transport in MCNP using the *F8 electron tally gave the overall best agreement with the experimental results. The *F8 tally, however, is much more computationally intensive than the neutron/gamma-ray transport calculations. Also, this report contains parametric analyses that examine the ways to improve the current design of the calorimeters. One finding from the parametric analysis was that the TC should be placed closer to the outer radius of the disks to obtain a measurement closer to the maximum temperature of the disk. Also, the parametric analysis showed that the most dominant mechanism of heat loss in the calorimeters is conduction through the alumina posts. In future designs, the conduction should be minimized to reduce the effect of heat loss on the measurements.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1993

Current capabilities for simulating the extreme distortion of thin structures subjected to severe impacts

The explicit transient dynamics technology in use today for simulating the impact and subsequent transient dynamic response of a structure has its origins in the 'hydrocodes' dating back to the late 1940's. The growth in capability in explicit transient dynamics technology parallels the growth in speed and size of digital computers. Computer software for simulating the explicit transient dynamic response of a structure is characterized by algorithms that use a large number of small steps. In explicit transient dynamics software there is a significant emphasis on speed and simplicity. The finite element technology used to generate the spatial discretization of a structure is based on a compromise between completeness of the representation for the physical processes modelled and speed in execution. That is, since it is expected in every calculation that the deformation will be finite and the material will be strained beyond the elastic range, the geometry and the associated gradient operators must be reconstructed, as well as complex stress-strain models evaluated at every time step. As a result, finite elements derived for explicit transient dynamics software use the simplest and barest constructions possible for computational efficiency while retaining an essential representation of the physical behavior. The best example of this technology is the four-node bending quadrilateral derived by Belytschko, Lin and Tsay. Today, the speed, memory capacity and availability of computer hardware allows a number of the previously used algorithms to be 'improved.' That is, it is possible with today's computing hardware to modify many of the standard algorithms to improve their representation of the physical process at the expense of added complexity and computational effort. The purpose is to review a number of these algorithms and identify the improvements possible. In many instances, both the older, faster version of the algorithm and the improved and somewhat slower version of the algorithm are found implemented together in software. Specifically, the following seven algorithmic items are examined: the invariant time derivatives of stress used in material models expressed in rate form; incremental objectivity and strain used in the numerical integration of the material models; the use of one-point element integration versus mean quadrature; shell elements used to represent the behavior of thin structural components; beam elements based on stress-resultant plasticity versus cross-section integration; the fidelity of elastic-plastic material models in their representation of ductile metals; and the use of Courant subcycling to reduce computational effort.

Key, Samuel W.↗
DOE OSTI2017

Coupling Experiments and Simulations to Understand Metallic Fuel Behavior

The development of metal alloy fuels for minor actinide transmutation in fast reactor systems is of high importance. However, it faces numerous technical challenges such as: i) fuel behavior has to be benign during core off-normal events; ii) fuel has to maintain integrity to high burnup; iii) low minor actinide fabrication loss under remote handling conditions; and iv) fuel has to lend itself to low-loss recycling process. Limited experimental fuel performance data is currently available to support the selection of a fuel type for a transmutation system. Metallic fuels are being considered for application in fast reactor environment due to the high burn-up, high fissile and fertile density capability, inherent safety characteristics, favorable thermal response, and relatively easy recycling process. The uranium-plutonium-zirconium (U-Pu- Zr) based alloys are considered as one of the most promising metallic fuels for such application. In this work, we present results of microstructural characterization of Pu-Zr and U-Pu-Zr alloys and provide much needed linkage between thermal properties and microscopic scale observation. A firm understanding of microstructural evolution in reactor materials is imperative and will lead to deeper understanding of in-reactor degradation behavior, creating a new generation of radiation-resistant materials and fuels with the potential for significant increases in burn-up.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Application of GOTHIC to Groundwater Transport Analysis - 20152

Migration of pollutants and hazardous wastes, potentially containing radioactive isotopes, via groundwater transport is a concern at most waste cleanup sites. Predictive analysis can be used to evaluate mitigating actions intended to minimize impact on the environment and public exposure. GOTHIC is a multipurpose thermal hydraulics code that is used extensively in the nuclear industry for design, licensing and operation evaluations. It combines the capabilities of typical one-dimensional system codes and the essential features of Computational Fluid Dynamics (CFD) codes for three-dimensional analysis. There are other codes that are specifically developed for groundwater transport analysis and the results presented here are consistent with prior analyses. However, GOTHIC has some unique features that offer advantages for applications related to nuclear waste. Most importantly, it has been developed and maintained under a Quality Assurance program in compliance with the requirements of 10CFR50 Appendix B [2] and applicable portions of ASME NQA-1 [3] since 1995. Available GOTHIC capabilities that make the code especially useful for groundwater transport of nuclear materials include: - Tracking of any number of tracer elements for contaminants and other species of interest; - Radioactive decay and progeny of tracer elements; - Adsorption/desorption of tracer elements; - Tracking of any number of dissolved gases; - Release and absorption for dissolved gases; - Vapor phase tracking; - Non-Newtonian fluid modeling. The general porous body modeling approach makes GOTHIC well suited to groundwater transport analysis. The multi-region modeling approach used by GOTHIC simplifies model construction for regions of varying hydrologic characteristics and focuses the computational effort on regions of particular interest while simultaneously capturing the macroscopic response and any feedback effects across the larger domain. The applicability of GOTHIC to groundwater transport applications is demonstrated by comparing code results with available analytic or semi-analytic solutions for groundwater behavior. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2013

Uncertainty Quantification of Calculated Temperatures for the AGR-1 Experiment

This report documents an effort to quantify the uncertainty of the calculated temperature data for the first Advanced Gas Reactor (AGR-1) fuel irradiation experiment conducted in the INL’s Advanced Test Reactor (ATR) in support of the Next Generation Nuclear Plant (NGNP) R&D program. Recognizing uncertainties inherent in physics and thermal simulations of the AGR-1 test, the results of the numerical simulations can be used in combination with the statistical analysis methods to improve qualification of measured data. Additionally, the temperature simulation data for AGR tests can be used for validation of the fuel transport and fuel performance simulation models. The crucial roles of the calculated fuel temperatures in ensuring achievement of the AGR experimental program objectives require accurate determination of the model temperature uncertainties. The report is organized into three chapters. Chapter 1 introduces the AGR Fuel Development and Qualification program and provides overviews of AGR-1 measured data, AGR-1 test configuration and test procedure, and thermal simulation. Chapters 2 describes the uncertainty quantification procedure for temperature simulation data of the AGR-1 experiment, namely, (i) identify and quantify uncertainty sources; (ii) perform sensitivity analysis for several thermal test conditions; (iii) use uncertainty propagation to quantify overall response temperature uncertainty. A set of issues associated with modeling uncertainties resulting from the expert assessments are identified. This also includes the experimental design to estimate the main effects and interactions of the important thermal model parameters. Chapter 3 presents the overall uncertainty results for the six AGR-1 capsules. This includes uncertainties for the daily volume-average and peak fuel temperatures, daily average temperatures at TC locations, and time-average volume-average and time-average peak fuel temperatures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Rapid activation of non-oriented mechanophores via shock loading and spallation

Mechanophores (MPs), stimuli-responsive molecules that respond chromatically to mechanochemical reactions, are important for understanding the coupling between mechanics and chemistry as well as in engineering applications. However, the atomic-level understanding of their activation originates from gas phase studies or under simple linear elongation forces directly on molecules or polymer chains containing MPs. The effect of many-body distortions, pervasive in condensed-phase applications, is not understood. Therefore, we performed large-scale molecular dynamics (MD) simulations of a poly(methyl methacrylate)-spiropyran copolymer under dynamic mechanical loading and studied the activation of the MP under various conditions from dynamical compression to tension during unloading. Detailed analysis of the all-atom MD trajectories shows that the MP blocks experience significant many-body intramolecular distortion that can significantly decrease the activation barrier as compared with when deformation rates are slow relative to molecular relaxation time scales. Here we find that the reactivity of MPs under material compression states is governed by many-body effects of intramolecular torsions, whereas under tension, the reactions are governed by tensile stresses.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Injection data analysis using material balance time for CO 2 storage capacity estimation in deep closed saline aquifers

Estimating the ultimate storage capacity of deep saline aquifers is important to address the formation potential to store the envisioned large volumes of CO 2 . Injection data (i.e. injection rate, bottomhole pressure, and cumulative injected volume of CO 2 ) are routinely recorded during storage operations. These data contain valuable information on the subsurface (e.g. the reservoir pore volume and the formation storage capacity) that can be extracted. In this paper, we present a two-step graphical technique to infer the pore volume and the ultimate storage capacity of closed saline aquifers by analyzing the available injection data. First, the pore volume is inferred through adapting the concept of the material balance time. Material balance time is an approximate superposition time function developed to interpret production data from oil and gas wells operating at variable pressure/rate conditions during the boundary-dominated flow period. Using material balance techniques, the ultimate storage capacity is then estimated through linear extrapolation of the average pressure trend to the maximum allowable pressure the formation can withstand. The average pressure is not available in practice, but is can be obtained from the injection data. Two approaches are presented in this study to calculate the average pressure; namely the rigorous and the approximate approaches. Unlike the rigorous approach, the approximate approach does not require a prior knowledge of some reservoir properties (e.g. relative permeability, absolute permeability, formation porosity and thickness) to calculate the average pressure. To investigate its potential and reliability in analyzing CO 2 injection data, the proposed technique is applied to four synthetic cases representing different well operating conditions. Results indicate that the approximate approach consistently overestimates the actual (simulated) storage capacity as compared to the rigorous approach. The agreement - between the inferred and the simulated reservoir pore volume, and between the analytical and numerical estimates of storage capacity - validates the potential application of the technique to CO 2 storage in closed saline aquifers. The technique is further substantiated through application to a field data set utilized from a commercial-scale geological storage (CGS) project. Finally, field data interpretation shows that the proposed technique can be utilized to identify the degree of hydraulic continuity and reservoir compartmentalization within a target formation by interpreting the corresponding pressure and rate responses.

02 PETROLEUM↗
DOE OSTI2017

Modeling slip system strength evolution in Ti-7Al informed by in-situ grain stress measurements

In this work, far-field high-energy X-ray diffraction microscopy is used to asses the evolution of slip system strengths in hexagonal close-packed (HCP) Ti-7Al during tensile deformation in-situ. The following HCP slip system families are considered: basal $\langle$a$\rangle$, prismatic $\langle$a$\rangle$, pyramidal$\langle$a$\rangle$, and first-order pyramidal $\langle$c + a$\rangle$. A 1 mm length of the specimen's gauge section, marked with fiducials and comprised of an aggregate of over 500 grains, is tracked during continuous deformation. The response of each slip system family is quantified using ‘slip system strength curves’ that are calculated from the average stress tensors of each grain over the applied deformation history. These curves, which plot the average resolved shear stress for each slip system family versus macroscopic strain, represent a mesoscopic characterization of the aggregate response. A short time-scale transient softening is observed in the basal $\langle$a$\rangle$, prismatic $\langle$a$\rangle$, and pyramidal $\langle$a$\rangle$ slip systems, while a long time-scale transient hardening is observed in the pyramidal $\langle$c + a$\rangle$ slip systems. These results are used to develop a slip system strength model as part of an elasto-viscoplastic constitutive model for the single crystal behavior. A suite of finite element simulations is performed on a virtual polycrystal to demonstrate the relative effects of the different parameters in the slip system strength model. The model is shown to accurately capture the macroscopic stress-strain response using parameters that are chosen to capture the mesoscopic slip system responses.

36 MATERIALS SCIENCE↗
DOE OSTI2018

A continuum dislocation-based model of wedge microindentation of single crystals

Recent Electron Backscatter Diffraction (EBSD) experiments have revealed the emergence of heterogeneous dislocation microstructures forming under a wedge indenter in fcc crystals, where micro-meter dislocation patterns challenge the predictions of traditional models of plasticity. In order to explain the formation of these features and develop a relationship between the force-displacement curve and the dislocation substructure, here we present a model of wedge indentation based on the continuum theory of dislocations. The model accounts for large deformation kinematics through the multiplicative split of the deformation gradient tensor, where the incompatible plastic component of deformation results from the flux of dislocations on different and interacting slips systems. Constitutive equations for dislocation fluxes are determined from a dissipative variational principle. As a result, each dislocation density satisfies an initial-boundary value problem with convective-diffusive character, which is coupled to the macroscopic stress and displacement fields governing the deformation process. Solution to the self-consistent continuum formulation is found using the finite element method. Computer simulations mimic the experimental conditions of wedge micro-indentation experiments into Ni single-crystals used by Kysar et al. (2010a). A comparison of overall dislocation density distribution and macroscopic mechanical response shows good overall agreement with the experimental results in terms of the detailed features of dislocation patterns and lattice rotations as well as the macroscopic force-displacement response.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effect of lattice orientation on compressive properties of selective laser sintered nylon lattice coupons

As AM lattices become more popular in medical devices, it is important to consider how lattice design parameters affect the mechanical integrity and performance of a device. This research investigated the effect of lattice orientation on the compressive mechanical response of five common lattice geometries: Hexagonal Honeycomb (Hex), Diamond, Voronoi Tessellation Method (VTM), sheet-based Gyroid, and Face-Centered-Cubic (FCC). Samples were tested at two relative densities and eight orientations, printed in Nylon 12 (PA2200) on an EOS P396. The mechanical response was compared to simulated finite element analysis (FEA) response for Hex and FCC lattices. Hex lattices displayed significant orientation dependance. Diamond lattices demonstrated some orientation dependance when at a lower relative density. Gyroid, FCC, and VTM all displayed minimal orientation dependance. Increasing lattice relative density appeared to reduce anisotropy among tested orientations of all lattice geometries except Hex. FEA simulations were able to produce a trend response that matched the mechanical response trends based on the investigated orientations. Furthermore, these results support the need for multi axial lattice mechanical test considerations. It also bolsters the need for lattice FEA models to have a sufficient comparator when identifying “worst-case” lattice test scenarios.

36 MATERIALS SCIENCE↗
DOE OSTI2024

In-situ HRTEM observations of intermediate phase transformation in lattice reorientation in HCP rhenium

In-situ high-resolution transmission electron microscopy (HRTEM) is performed to investigate the deformation behavior of hexagonal close-packed rhenium (Re) which is compressed along the $\langle1\overline{1}00\rangle$ direction. Atomistic simulations are also conducted to better understand the deformation mechanisms. Two types of lattice reorientation are observed during compression. The first type involves the reorientation of one lattice by ~90° around $\langle1\overline{1}00\rangle$, which is accomplished by the formation of an intermediate face-center-cubic (FCC) phase at the interface. This transformation sequence can be described as {$1\overline{1}00$} matrix → {$111$} FCC → ($0001$) twin . In the second type, a new grain is formed but does not satisfy any known twin relationship with the matrix, and an intermediate FCC phase is also formed. The transformation sequence can be described as {$1\overline{1}00$} matrix → {$111$} FCC → ($0001$) grain . Mechanisms responsible for the observed lattice reorientation and sequential phase transitions are analyzed by conducting lattice correspondence analyses on the simulation results. Strain accommodation is also analyzed to explain the mechanisms for lattice reorientation and the intermediate phase transformations. In conclusion, the results provide new insight into the deformation behavior of HCP metals

36 MATERIALS SCIENCE↗
DOE OSTI2019

Time-Dependent Complete Active Space Embedded in Polarizable Force Field

Electron correlation and environmental effects play important roles in electron dynamics and spectroscopic observables of chemical systems in condensed phase. In this paper, we present a time-dependent complete active space configuration interaction (TD-CASCI) approach embedded in a polaraizable force filed, MMPol. The present implementation of TD-CASCI/MMPol utilizes a direct matrix-vector contraction, allowing studies of large systems. This scheme is used to study the solvatochromic shift of coumarin 153 in methanol. The TD-CASCI/MMPol approach captures the double excitation character in the excited state wave function and accurately predicts the solvatochromic red-shift of coumarin 153 dye within the experimental range, outperforming linear response time-dependent density functional theory. The effect of using different reference orbitals for the TD-CASCI/MMPol simulation is also investigated, highlighting the need for an unbiased treatment of all electronic states in the energy range of interest.

Liu, Hongbin↗
DOE OSTI2022

Experimental study of electromagnetic wave scattering from a gyrotropic gaseous plasma column

We experimentally demonstrate the controlled scattering of incident transverse-electric electromagnetic waves from a gyrotropic magnetized plasma cylindrical discharge. Scattered electromagnetic waves can bend left and right by changing the external magnetic field of a plasma rod. Measured scattered wavefronts are in good agreement with electromagnetic simulations. A gyrotropic response is observed for incident wave frequencies ranging from 3.5 to 5.6 GHz for conditions corresponding to a ratio of cyclotron frequency to plasma frequency, ω ce /ω p ≈ 0.16. In conclusion, the observation of a gyrotropic response from cylindrical plasma discharges paves the way for their use as building blocks for future devices such as magnetized plasma photonic crystals, topological insulators, plasma metamaterials, non-reciprocal waveguide structures, and other devices, which require a tunable gyrotropic response from centimeter to meter-scale materials with application-specific geometry.

42 ENGINEERING↗
DOE OSTI2014

Study of the Durability of Doped Lanthanum Manganite and Cobaltite Cathode Materials under ''Real World'' Air Exposure Atmospheres

The overall objective of the program is to develop and validate mechanisms responsible for the overall structural and chemical degradation of lanthanum manganite as well as lanthanum ferrite cobaltite based cathode when exposed to “real world” air atmosphere exposure conditions during SOFC systems operation. Of particular interest are the evaluation and analysis of degradation phenomena related to and responsible for (a) products formation and interactions with air contaminants, (b) dopant segregation and oxide exolution at free surfaces, (c) cation interdiffusion and reaction products formation at the buried interfaces, (d) interface morphology changes, lattice transformation and the development of interfacial porosity and (e) micro-cracking and delamination from the stack repeat units. Reaction processes have been studied using electrochemical and high temperature materials compatibility tests followed by structural and chemical characterization. Degradation hypothesis has been proposed and validated through further experimentation and computational simulation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Updating the Building Science Advisor (BSA): A Tool to Assist in the Design of Durable Building Envelopes

Predicting the moisture durability of building envelope components remains challenging due to multiple influencing factors, including material selection, assembly positioning, local climate conditions, air tightness, interior environment, and construction quality. Building codes increasingly emphasize energy efficiency through enhanced insulation and tighter envelopes but offer limited guidance on moisture durability considerations. Consequently, builders face uncertainty, particularly as new materials and assemblies enter the market.The Building Science Advisor (BSA) is a free, web-based expert system developed to address these challenges by providing actionable insights into the moisture durability and energy efficiency of both new and retrofit wall designs. Recently updated, we are now providing version 3.0 of the tool. BSA features significant user interface improvements, enhancing navigation and user interaction through a refreshed, intuitive design. Additionally, the tool incorporates a newly developed database containing pre-simulated wall assembly cases, significantly reducing response times and improving the accuracy of moisture durability assessments. Furthermore, the updated BSA includes moisture content as a performance criterion, providing users with a more comprehensive understanding of moisture-related durability risks. These enhancements enable rapid, reliable assessments tailored to specific climate zones and local building practices. BSA continues to offer targeted guidance on wall retrofit scenarios and delivers access to an expanded library of location-specific building science resources.This paper describes these key updates, highlighting the enhanced features, expanded capabilities, and overall improvements to user experience and educational content. The paper includes a demonstration that illustrates how the revised BSA effectively supports practitioners in designing durable, energy-efficient building envelope assemblies.

Salonvaara, Mikael [ORNL] (ORCID:0000000318991554)↗
NASA NTRS1987

Structural modeling for multicell composite rotor blades

Composite material systems are currently good candidates for aerospace structures, primarily for the design flexibility they offer, i.e., it is possible to tailor the material and manufacturing approach to the application. A working definition of elastic or structural tailoring is the use of structural concept, fiber orientation, ply stacking sequence, and a blend of materials to achieve specific performance goals. In the design process, choices of materials and dimensions are made which produce specific response characteristics, and which permit the selected goals to be achieved. Common choices for tailoring goals are preventing instabilities or vibration resonances or enhancing damage tolerance. An essential, enabling factor in the design of tailored composite structures is structural modeling that accurately, but simply, characterizes response. The objective of this paper is to present a new multicell beam model for composite rotor blades and to validate predictions based on the new model by comparison with a finite element simulation in three benchmark static load cases.

Rehfield, Lawrence W.↗
NASA NTRS2016

Overview of the CHarring Ablator Response (CHAR) Code

An overview of the capabilities of the CHarring Ablator Response (CHAR) code is presented. CHAR is a one-, two-, and three-dimensional unstructured continuous Galerkin finite-element heat conduction and ablation solver with both direct and inverse modes. Additionally, CHAR includes a coupled linear thermoelastic solver for determination of internal stresses induced from the temperature field and surface loading. Background on the development process, governing equations, material models, discretization techniques, and numerical methods is provided. Special focus is put on the available boundary conditions including thermochemical ablation and contact interfaces, and example simulations are included. Finally, a discussion of ongoing development efforts is presented.

Amar, Adam J.↗
NASA NTRSDate not supplied

Icarus Application to Dragonfly Heatshield

Icarus, an in-house material response solver developed at NASA Ames, is applied to perform 2-D bondline-temperature analysis on the heatshield surface of the latest Dragonfly entry system design for a mission to Titan. Two axisymmetric domains are used: the near-shoulder region and the entire heatshield. The heatshield consists of a PICA TPS bonded to a layered system of Aluminum honeycomb wrapped by an M55J carbon-fiber facesheet. The backshell, which is also included for more accurate analysis, consists of similar materials layered. Icarus simulations using both domains indicate that the in-plane thermal conductivity property of the M55J facesheet is observed to play a dominant role in bondline temperature. Similarly, the maximum bondline temperature is found around the PICA-tile-interface region instead of the near-shoulder region or the stagnation point based on the current trajectory using the orthotropic properties of M55J. Therefore, these findings indicate the significance of 2-D or higher-dimensional material modeling analysis in fully understanding how bondline temperature behaves across the heatshield and obtaining a basis for TPS design locations. Verification at different radial locations ahead of the shoulder, where thermal conduction is close to 1-D, between Icarus (2-D), Icarus (1-D), and FIAT is performed and good agreement is observed. With the entire heatshield domain, the total mass loss of the aeroshell materials is approximately 3% of their initial mass potentially due to thermal decomposition of phenolic resins inside PICA under the current high-heating environment.

TPS material modeling↗