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 505 records · Page 28

Simulated rarefied entry of the Galileo probe into the atmosphere of Jupiter

Flow properties and aerodynamics are computed with a direct simulation Monte Carlo (DSMC) method for rarefied entry of the Galileo Probe into the atmosphere of Jupiter. Accurate predictions of vehicle drag coefficients are needed in order to assess atmospheric properties from the onboard Atmosphere Structure Experiment where highly-sensitive accelerometers will measure the drag force to within 10-6 barr during the initial entry phase at high altitudes. The corresponding flow rarefraction extends from the free molecule limit to the near continuum transition regime (Re less than 1000). Simulation results indicate that C(sub D) varies from 2.1 at the free molecule limit down to 1.6 at Re(infinity) = 1,000. Temperatures, densities, and internal energies throughout the flow field were also computed at each altitude ranging from 735 km to 353 km above the 1 barr level in the Jovian atmosphere. Surface heating and temperatures of the probe were computed directly in the DSMC code by assuming radiative equilibrium. Material response was re-asssessed accurately during entry by accounting for conductivity, heat capacity, and pyrolysis which led to surface material mass efflux several times that of the freestream mass influx. The simulation also accounted for the quantum nature of the rotational energy mode of the dominant atmospheric species H2 through partial internal excitation in the freestream gas.

Hass, Brian L.↗

Comparison of Simulated and Experimentally Measured Detector Impulse Responses to Pulsed X-rays

The Gamma Array Simulation Toolkit (GAST) is a multi-physics software stack comprised of Geant4, MATLAB, and experimentally measured photomultiplier tube (PMT) impulse response waveforms that is used to model and inform the end-to-end performance of radiation detector systems. In order to validate the ability of GAST to predict the impulse response of a modeled detector configuration operated in current mode, experimental measurements using the Nevada National Security Site Transformational Diagnostics and Imaging 2 MeV endpoint energy electron linear accelerator (linac) with a tungsten target were performed. Seven different detector configurations consisting of unique combinations of scintillator material, size, reflectivity conditions, presence or absence of light guide, and PMT were evaluated. The resultant impulse response of each detector configuration was recorded and compared to a GAST simulation that emulated the respective linac experiment conditions. The comparison of the experimental and simulated detector impulse responses indicated that GAST can reliably predict detector impulse responses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Unstructured mesh-based multi-physics reactor transient simulation with iMC code

In this study, we present a time-dependent Monte Carlo simulation result about a complex fuel-loaded light water reactor system with thermal-hydraulic feedback using the KAIST iMC code. The intra-pin temperature distribution is calculated based on Monte Carlo-tallied detailed power distribution and the finite element heat transfer method. The temperature effect of fuel, absorber, and coolant are adequately reflected into the dynamic response by the on-the-fly Doppler broadened cross-section correction. For the time-dependent Monte Carlo transport simulation, the predictor-corrector quasi-static Monte Carlo method is used for this study. We tested the coupled analysis framework with a 17-by-17 fuel assembly loaded with the centrally-shielded burnable absorber (CSBA) fuel element. The temperature feedback effect on the material cross-section and the moderator density incurred negative reactivity in response to the external positive reactivity insertion, showing the inherent stability of the reactor system. The pin-wise power and temperature change during the system transient are presented and discussed. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

How relative defect migration energies drive contrasting temperature-dependent microstructural evolution in irradiated ceramics

Ceramic materials have become widely used in various fields of material science, and ceramic oxides such as cubic zirconia (c- ZrO 2 ) and magnesia (MgO) are candidate materials for nuclear energy applications. For the corresponding in-service conditions of these materials, there is a crucial need in studies at moderate or high temperatures of the physical phenomena underlying the damage buildup under irradiation. Here in the present work, we show, using x-ray diffraction, that these two materials exhibit a similar damage accumulation process under ion irradiation at fixed temperature. However, they display an unexpected opposite damaging rate to changes in the irradiation temperature. In fact, as the temperature is increased, the final damage state is reached earlier in c- ZrO 2 while it is delayed in MgO. Rate equation cluster dynamics simulations show that defect clustering is favored over defect recombination in c- ZrO 2 , but the situation is reversed for MgO, explaining the observed opposite response to temperature of the two materials. This contrasting behavior can be rationalized in terms of nonequivalent interstitial versus vacancy defect migration energies in MgO. We finally demonstrate that these results allow for a qualitative prediction of the evolution of the experimental irradiation-induced disorder with temperature, henceforth potentially reducing the cost in selecting and developing ad hoc materials.

36 MATERIALS SCIENCE↗

Predicting System Accidents with Model Analysis During Hybrid Simulation

Standard discrete event simulation is commonly used to identify system bottlenecks and starving and blocking conditions in resources and services. The CONFIG hybrid discrete/continuous simulation tool can simulate such conditions in combination with inputs external to the simulation. This provides a means for evaluating the vulnerability to system accidents of a system's design, operating procedures, and control software. System accidents are brought about by complex unexpected interactions among multiple system failures , faulty or misleading sensor data, and inappropriate responses of human operators or software. The flows of resource and product materials play a central role in the hazardous situations that may arise in fluid transport and processing systems. We describe the capabilities of CONFIG for simulation-time linear circuit analysis of fluid flows in the context of model-based hazard analysis. We focus on how CONFIG simulates the static stresses in systems of flow. Unlike other flow-related properties, static stresses (or static potentials) cannot be represented by a set of state equations. The distribution of static stresses is dependent on the specific history of operations performed on a system. We discuss the use of this type of information in hazard analysis of system designs.

Malin, Jane T.↗

Uncertainty quantification in multiscale simulation of woven fiber composites

We report woven fiber composites have been increasingly employed as light-weight materials in aerospace, construction, and transportation industries due to their superior properties. These materials possess a hierarchical structure that necessitates the use of multiscale simulations in their modeling. To account for the inherent uncertainty in materials, such simulations must be integrated with statistical uncertainty quantification (UQ) and propagation (UP) methods. However, limited advancement has been made in this regard due to the significant computational costs and complexities in modeling spatially correlated structural variations coupled at different scales. In this work, a non-intrusive approach is proposed for multiscale UQ and UP to address these limitations. We introduce the top-down sampling method that allows to model non-stationary and continuous (but not differentiable) spatial variations of uncertainty sources by creating nested random fields (RFs) where the hyperparameters of an ensemble of RFs is characterized by yet another RF. We employ multi-response Gaussian RFs in top-down sampling and leverage statistical techniques (such as metamodeling and dimensionality reduction) to address the considerable computational costs of multiscale simulations. We apply our approach to quantify the uncertainty in a cured woven composite due to spatial variations of yarn angle, fiber volume fraction, and fiber misalignment angle. Finally, our results indicate that, even in linear analysis, the effect of uncertainty sources on the material’s response could be significant.

42 ENGINEERING↗

The Influence of Mesh Density on the Impact Response of a Shuttle Leading-Edge Panel Finite Element Simulation

A study was performed to examine the influence of varying mesh density on an LS-DYNA simulation of a rectangular-shaped foam projectile impacting the space shuttle leading edge Panel 6. The shuttle leading-edge panels are fabricated of reinforced carbon-carbon (RCC) material. During the study, nine cases were executed with all possible combinations of coarse, baseline, and fine meshes of the foam and panel. For each simulation, the same material properties and impact conditions were specified and only the mesh density was varied. In the baseline model, the shell elements representing the RCC panel are approximately 0.2-in. on edge, whereas the foam elements are about 0.5-in. on edge. The element nominal edge-length for the baseline panel was halved to create a fine panel (0.1-in. edge length) mesh and doubled to create a coarse panel (0.4-in. edge length) mesh. In addition, the element nominal edge-length of the baseline foam projectile was halved (0.25-in. edge length) to create a fine foam mesh and doubled (1.0- in. edge length) to create a coarse foam mesh. The initial impact velocity of the foam was 775 ft/s. The simulations were executed in LS-DYNA version 960 for 6 ms of simulation time. Contour plots of resultant panel displacement and effective stress in the foam were compared at five discrete time intervals. Also, time-history responses of internal and kinetic energy of the panel, kinetic and hourglass energy of the foam, and resultant contact force were plotted to determine the influence of mesh density. As a final comparison, the model with a fine panel and fine foam mesh was executed with slightly different material properties for the RCC. For this model, the average degraded properties of the RCC were replaced with the maximum degraded properties. Similar comparisons of panel and foam responses were made for the average and maximum degraded models.

Jackson, Karen E.↗

Humidity Disrupts Structural and Chiroptical Properties of Chiral 2D Perovskites

Chiral two-dimensional (2D) hybrid organic–inorganic metal halide perovskite semiconductors have emerged as an exceptional material platform with many design opportunities for spintronic applications. However, a comprehensive understanding of changes to the crystal structure and chiroptical properties upon exposure to atmospheric humidity has not been established. We demonstrate phase degradation to the 1D (MBA)PbI 3 (MBA = methylbenzylammonium) and the hypothetical (MBA) 3 PbI 5 ·H 2 O hydrate phases, accompanied by a reduction and disappearance of the chiroptical response. Firstprinciple simulations show that water molecules preferentially locate at the interface between the organic cations and the inorganic framework, thereby disrupting the hydrogen bonding, impacting both the structural chirality and stability of the material. These findings provide critical insights into phase degradation mechanisms and their impact on chiroptical activity in chiral 2D perovskites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Shock compression response of an insensitive high explosive single crystal: 1,1-diamino-2,2-dinitroethene (FOX-7)

Insensitive high explosives (IHEs) – that do not compromise performance – are of immense interest as a safer alternative to conventional high explosives, such as pentaerythritol tetranitrate (PETN) and cyclotrimethylenetrinitramine (RDX). Despite the strong interest in using IHEs, shock compression experiments on IHE single crystals have not been reported. To address this need, plate impact experiments were conducted to measure wave profiles in 1,1-diamino-2,2-dinitroethene (FOX-7) single crystals – a representative IHE crystal – shocked to 21 GPa longitudinal stress. Particle velocity histories, measured using laser interferometry, show a clear two-wave structure (elastic-inelastic response) at modest stresses (<3.8 GPa). Wave profiles at higher stresses show a single (overdriven) wave. Measured shock velocities and wave profiles provide accurate Hugoniot data to 21 GPa. The measured wave profiles to 21 GPa show no sign of energy release due to chemical decomposition and constitute the first demonstration of an IHE single crystal insensitivity under plane shock compression. Numerical simulations using a phenomenological material model developed for FOX-7 showed good agreement with the measured wave profiles. The experimental results and continuum simulations presented here constitute a significant first step in gaining insight into the shock compression response of IHE single crystals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Simulation of Radiation-Induced DNA Damage With the Code RITRACKS

INTRODUCTION DNA damage is one of the most physiologically important effects of ionizing radiation. Clustered DNA damage events, like double-strand breaks (DSBs), have the most notable biological consequences. DNA damage types depend on both the track structure of the radiation and the spatial organization of the DNA. High linear energy transfer (LET) charged nuclei, found in galactic cosmic rays (GCR), are known to produce large numbers of complex DNA damage events. The human genome is packaged into chromatin, which can take on locus-dependent and cell type-dependent spatial conformations that correspond to epigenetic states, such as more open, extended structures in transcriptionally active chromatin. These epigenetic differences can affect DNA break patterns in response to ionizing radiation, potentially creating distinct DNA repair and signaling outcomes across the genome in different cells. MATERIAL AND METHODS The code RITRACKS (Relativistic Ion Tracks), which simulates stochastic radiation track structures and radiation chemistry, was used to model damage on isolated and histone-bound DNA by various types of ions and photons. The changes made to the code to perform radiation-induced DNA damage, and simulation results on single nucleosomes are given in our recent paper. In this work, the DNA building capabilities of RITRACKS have been extended to simulate more complex DNA structures build on the coarse-grain simulation framework meso-WLCsim. This code can sample generic chromatin fiber conformation ensembles based on the geometry of nucleosomes and mechanical properties of DNA. Using RITRACKS, we simulated the fragment length distributions (FLD) of irradiated DNA structures built using the chromatin conformations of WLCsim and obtained results representative of those obtained with Radiation-Induced Correlated Cleavage with sequencing (RICC-Seq) experiments [6]. We have also performed Fe ion and photon irradiations of K562, IMR90, BJ and RPE-1 cells at NSRL to experimentally validate results. Sample processing and data analysis are in progress and any available preliminary results will be discussed. DISCUSSION The recent updates in the code RITRACKS allow the calculation of several quantities such as the DNA damage yield and the FLD. This approach can be used to model epigenetic state-specific chromatin structure parameters to leverage the epigenetic state data available for many human cell types to infer relative DNA damage sensitivity among genomic loci.

I Plante↗

A modular methodology for time-domain stochastic seismic wave propagation

Presented here is a modular methodology for time-domain stochastic seismic wave propagation analysis. Presented methodology is designed to analyse uncertain seismic motions as an input, propagating through uncertain material. Traditional approach for uncertain wave propagation relies on models that include deep bedrock, local soil site, and their random process and random field information. Such models can become quite large and computationally intractable. The modular approach proposed herein features two step approach that allows separate consideration of the deep bedrock and local site along with corresponding random field information. In this work, the first step considers an auxiliary stochastic motions problem in the bedrock. Stochastic local site response can then be simulated in a reduced domain within certain depth from the surface. Application of uncertain seismic motions at depth, for local uncertain site response is done using stochastic effective forces developed through the Domain Reduction Method. By using Hermite polynomial chaos expansion to represent the non-Gaussian random field of material parameters and non-stationary random process of seismic motion, the proposed modular methodology is formulated using intrusive stochastic Galerkin approach, as seen in the Stochastic Elastic–Plastic Finite Element Method (SEPFEM). Developed modular methodology is illustrated using a 1-D stochastic seismic wave propagation analysis with three cases, and simulation results are also verified with results from conventional approach.

58 GEOSCIENCES↗

Investigation of a Macromechanical Approach to Analyzing Triaxially-Braided Polymer Composites

A macro level finite element-based model has been developed to simulate the mechanical and impact response of triaxially-braided polymer matrix composites. In the analytical model, the triaxial braid architecture is simulated by using four parallel shell elements, each of which is modeled as a laminated composite. The commercial transient dynamic finite element code LS-DYNA is used to conduct the simulations, and a continuum damage mechanics model internal to LS-DYNA is used as the material constitutive model. The material stiffness and strength values required for the constitutive model are determined based on coupon level tests on the braided composite. Simulations of quasi-static coupon tests of a representative braided composite are conducted. Varying the strength values that are input to the material model is found to have a significant influence on the effective material response predicted by the finite element analysis, sometimes in ways that at first glance appear non-intuitive. A parametric study involving the input strength parameters provides guidance on how the analysis model can be improved.

Goldberg, Robert K.↗

Modeling of Response of Thermoplastic Composites Under High Energy Impact

This work investigates the response of a thermoplastic panel under high energy ballistic impact. The composite material system is a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK3). Several high energy dynamic impact tests were performed where an aluminum projectile impacted various panels at different impact velocities. For this work, the test performed at an impact velocity of 73.2 m/s was investigated. The projectile rebounded upon impact, and the damage to the panel was characterized through ultrasonic and CT-scanning. Response of the panel under impact and the associated damage is discussed. Numerical simulation of the impact process was performed using the LS-DYNA® computer code [1] with a new material model MAT_COMPOSITE_DAMAGE_DGD (MAT 299), a modification over MAT 261 material model in terms of an improved representation of damage coupling and additional capability for strain rate effects. Delaminations were modeled using cohesive zone contact formulation in LS-DYNA®, wherein the mixed-mode fracture is captured via the Benzegaggh-Kenane (B-K) mode-dependent fracture energy interpolation law. Simulation of this impact test demonstrates a good match of results in terms of panel deflection and panel-projectile interactions. The overall extent of delaminations was under-predicted by the simulation. This work forms the basis to model tests performed at higher velocities, involving projectile penetration of panel.

Composite Materials↗

Generation of Synthetic Data for a Radiation Detection Algorithm Competition

This paper details the generation of synthetic radiation data using large-scale Monte Carlo transport models to evaluate radiation search detection algorithms. Modular 3-D Monte Carlo models spanning multiple city blocks were constructed, loosely based on downtown Knoxville, TN, containing buildings composed of multiple materials (brick, granite, and concrete), sidewalks, a four-lane road, side streets, parking lots, and grassy fields. Background and simulated source detector response calculations from these models were used to create synthetic list mode data sets for a 2"×4"×16" NaI(Tl) detector moving through a city street at a constant speed. For the background simulations, major isotopes were computed individually so that background composition and variability could be computed efficiently outside Monte Carlo. The source detector response included six simulated sources placed at 15 source locations. Detector response was developed to be periodic through the city street so that a detector path could begin at one end of the model and wrap around to the other. This framework allowed for the creation of diverse data sets, each with its own unique background and simulated source detector response. Synthetic data allows for high-quality labels, which are useful in developing data-driven radiation detection algorithms. This methodology was used to create synthetic data sets which were released as part of a public data competition to spur the development of new radiation detection algorithms for radiological search applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Development of ASME Division 5 Code proposal on temperature limits for simplified design methods

This report establishes proposed upper temperature limits for the ASME BP&V Code Section III, Division 5, Nonmandatory Appendix HBB-T design by elastic analysis provisions for bounding ratcheting strain and creep-fatigue damage in Class A high temperature nuclear reactor components. Limitations on the use of these design options are required because the design by elastic analysis methods rely on bounding theories that assume a non-unified, decoupled model of creep-plasticity. However, at high temperatures creep and plastic deformation become coupled and bounding theorems relying a decoupled material response may fail. The report describes a method for selecting appropriate upper temperature limits, demonstrates directly through a comparison to full inelastic simulations that the existing Code provisions can be nonconservative at high temperatures, and develops the requisite Code language required to implement the temperature limits in Section III, Division 5.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Experimental Study of Mechanical Response of Polyurethane Rubber at Varying Strain Rates for W80-4 System

This report describes an experimental study to determine the mechanical behavior of the polyurethane rubber material that was used in the W80-4 MCC Shock/Breach Phase 1 test series. Compression experiments were conducted on cylindrical specimens over a wide range of loading rates to characterize the material over the range of strain rates that were experienced in the shock/breach testing. Additionally, specimen diameter was varied to determine the effect of confinement on the material response and was found to be significant. This data is used to populate a material model to enable accurate analyses and finite element simulations of the shock/breach test series.

36 MATERIALS SCIENCE↗

Space radiation effects on the dimensional stability of a toughened epoxy graphite composite

The effects of a simulated space radiation environment on the dimensional stability of an elastomer-toughened epoxy-graphite composite were determined. The response of the material was characterized following exposure to radiation doses equivalent to geosynchronous orbit lifetimes ranging from 6 months to 30 years. The results show that radiation interacts with the epoxy matrix to embrittle the composite, beginning at relatively low total doses (10 to the 7th power rads). The embrittlement results in thermal expansion changes and significant laminate microcracking during thermal cycling. These property changes could limit the service life of this material in some spacecraft applications.

Sykes, G. F.↗

PRISMS-Fatigue computational framework for fatigue analysis in polycrystalline metals and alloys

Abstract The PRISMS-Fatigue open-source framework for simulation-based analysis of microstructural influences on fatigue resistance for polycrystalline metals and alloys is presented here. The framework uses the crystal plasticity finite element method as its microstructure analysis tool and provides a highly efficient, scalable, flexible, and easy-to-use ICME community platform. The PRISMS-Fatigue framework is linked to different open-source software to instantiate microstructures, compute the material response, and assess fatigue indicator parameters. The performance of PRISMS-Fatigue is benchmarked against a similar framework implemented using ABAQUS. Results indicate that the multilevel parallelism scheme of PRISMS-Fatigue is more efficient and scalable than ABAQUS for large-scale fatigue simulations. The performance and flexibility of this framework is demonstrated with various examples that assess the driving force for fatigue crack formation of microstructures with different crystallographic textures, grain morphologies, and grain numbers, and under different multiaxial strain states, strain magnitudes, and boundary conditions.

Chemistry↗