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,333 records · Page 74

Lattice Rotation Patterns and Strain Gradient Effects in Face-Centered-Cubic Single Crystals Under Spherical Indentation

Strain gradient effects are commonly modeled as the origin of the size dependence of material strength, such as the dependence of indentation hardness on contact depth and spherical indenter radius. However, studies on the microstructural comparisons of experiments and theories are limited. First, we have extended a strain gradient Mises-plasticity model to its crystal plasticity version and implemented a finite element method to simulate the load-displacement response and the lattice rotation field of Cu single crystals under spherical indentation. The strain gradient simulations demonstrate that the forming of distinct sectors of positive and negative angles in the lattice rotation field is governed primarily by the slip geometry and crystallographic orientations, depending only weakly on strain gradient effects, although hardness depends strongly on strain gradients. Second, the lattice rotation simulations are compared quantitatively with micron resolution, three-dimensional X-ray microscopy (3DXM) measurements of the lattice rotation fields under 100mN force, 100 mu m radius spherical indentations in < 111 >, < 110 >, and < 001 > oriented Cu single crystals. Third, noting the limitation of continuum strain gradient crystal plasticity models, two-dimensional discrete dislocation simulation results suggest that the hardness in the nanocontact regime is governed synergistically by a combination of strain gradients and source-limited plasticity. However, the lattice rotation field in the discrete dislocation simulations is found to be insensitive to these two factors but to depend critically on dislocation obstacle densities and strengths.

Gao, Y. F.↗

Dynamics of defects in Ce 3+ doped silica affecting its performance as protective filter in ultraviolet high-power lasers

We investigate defects forming in Ce 3+ -doped fused silica samples following exposure to nanosecond ultraviolet laser pulses and their relaxation as a function of time and exposure to low intensity light at different wavelengths. A subset of these defects are responsible for inducing absorption in the visible and near infrared spectral range, which is of critical importance for the use of this material as ultraviolet light absorbing filter in high power laser systems. The dependence of the induced absorption as a function of laser fluence and methods to most efficiently mitigate this effect are presented. Experiments simulating the operation of the material as a UV protection filter for high power laser systems were performed in order to determine limitations and practical operational conditions.

47 OTHER INSTRUMENTATION↗

Thermal performance of a mechanically attached ablator tile for on-orbit repair of shuttle TPS

The reusable surface insulation (RSI) material used in the shuttle thermal protection system is susceptible to damage. If any RSI tiles are damaged or lost during ascent, they must be repaired or replaced prior to entry. One approach to replacing a damaged or missing RSI tile consists of mechanically attaching a tile of ablation material in the place of the RSI tile. The thermal performance of this type of repair tile was evaluated in a simulated entry heating environment. The test specimen consisted of the ablator repair tile mechanically fastened to the strain isolation pad and surrounded by RSI tiles. The evaluation of the thermal performance was based on temperature response, surface integrity, and predicted flight performance. When the ablator tile protruded 1/8 inch above the surrounding RSI tiles, the forward facing steps caused significant inflow of hot gas down the ablator RSI joints and this inflow caused greatly increased back surface temperatures.

Thompkins, S. S.↗

Monte Carlo simulation of lunar megaregolith and implications

A realistic Monte Carlo model closely simulating the evolution of the lunar megaregolith over a large area of 67 million sq. km of the front surface of the moon is presented. Craters larger than 100 km in diameter observed over the entire surface of the moon and those less than 100 km lying in the referenced area are included in the simulation. A total of 21,664 craters are processed. The model predicts the average thickness of the megaregolith to be about 1.9-2.0 km. Curves for the variation of the regolith thickness across the simulated area are given and show that about 50% of the area is covered with regolith less than 1 km thick. The model produces crater structures similar to the ones observed in the lunar highlands, it partially supports the layering theory for crater structures that the variations in strength and density of target materials may be responsible for the observed differences in the morphologies of lunar craters, and rules out the possibility that all craters when formed are bowl-shaped with a fixed depth/diameter ratio characteristic of small craters.

Aggarwal, H. R.↗

Trajectory Based Heating and Ablation Calculations for MESUR Pathfinder Aeroshell

Based on the geometry of Mars Environment Survey (MESUR) Pathfinder aeroshell and an estimated Mars entry trajectory, two-dimensional axisymmetric time dependent calculations have been obtained using GIANTS (Gauss-Siedel Implicit Aerothermodynamic Navier-Stokes code with Thermochemical Surface Conditions) code and CMA (Charring Material Thermal Response and Ablation) Program for heating analysis and heat shield material sizing. These two codes are interfaced using a loosely coupled technique. The flowfield and convective heat transfer coefficients are computed by the GIANTS code with a species balance condition for an ablating surface, and the time dependent in-depth conduction with surface blowing is simulated by the CMA code with a complete surface energy balance condition. In this study, SLA-561V has been selected as heat shield material. The solutions, including the minimum heat shield thicknesses over aeroshell forebody, pyrolysis gas blowing rates, surface heat fluxes and temperature distributions, flowfield, and in-depth temperature history of SLA-561V, are presented and discussed in detail.

Chen, Y. K.↗

An Approximate Ablative Thermal Protection System Sizing Tool for Entry System Design

A computer tool to perform entry vehicle ablative thermal protection systems sizing has been developed. Two options for calculating the thermal response are incorporated into the tool. One, an industry-standard, high-fidelity ablation and thermal response program was integrated into the tool, making use of simulated trajectory data to calculate its boundary conditions at the ablating surface. Second, an approximate method that uses heat of ablation data to estimate heat shield recession during entry has been coupled to a one-dimensional finite-difference calculation that calculates the in-depth thermal response. The in-depth solution accounts for material decomposition, but does not account for pyrolysis gas energy absorption through the material. Engineering correlations are used to estimate stagnation point convective and radiative heating as a function of time. The sizing tool calculates recovery enthalpy, wall enthalpy, surface pressure, and heat transfer coefficient. Verification of this tool is performed by comparison to past thermal protection system sizings for the Mars Pathfinder and Stardust entry systems and calculations are performed for an Apollo capsule entering the atmosphere at lunar and Mars return speeds.

Dec, John A.↗

An Approximate Ablative Thermal Protection System Sizing Tool for Entry System Design

A computer tool to perform entry vehicle ablative thermal protection systems sizing has been developed. Two options for calculating the thermal response are incorporated into the tool. One, an industry-standard, high-fidelity ablation and thermal response program was integrated into the tool, making use of simulated trajectory data to calculate its boundary conditions at the ablating surface. Second, an approximate method that uses heat of ablation data to estimate heat shield recession during entry has been coupled to a one-dimensional finite-difference calculation that calculates the in-depth thermal response. The in-depth solution accounts for material decomposition, but does not account for pyrolysis gas energy absorption through the material. Engineering correlations are used to estimate stagnation point convective and radiative heating as a function of time. The sizing tool calculates recovery enthalpy, wall enthalpy, surface pressure, and heat transfer coefficient. Verification of this tool is performed by comparison to past thermal protection system sizings for the Mars Pathfinder and Stardust entry systems and calculations are performed for an Apollo capsule entering the atmosphere at lunar and Mars return speeds.

Dec, John A.↗

The Evaluation of Composite Energy Absorbers for Use in UAM EVTOL Vehicle Impact Attenuation

The demand for new personal air-taxi services is leading to the development of lightweight Vertical Take-off and Landing (eVTOL) vehicles with electric propulsion for the Urban Air Mobility (UAM) industry. Manufacturers (OEMs) are considering many different designs to develop a vehicle that is able to take-off, cruise, and land autonomously with seating arrangements ranging between 2 and 15 passengers. It is unclear at present how the eventual market will mature; however, one of the common design characteristics noted by many of the OEMs is the use of advanced materials such as composites. A test and analysis program was initiated at NASA Langley Research Center (LaRC) in 2018 to evaluate the impact attenuation capabilities of various composite material systems with the goal of eventual implementation into an eVTOL vehicle. A series of 3-inch diameter by 6-inch length tubular specimens were fabricated from different material systems which included both traditional carbon and hybrid woven layers of fibers. Additionally, a subset of specimens were filled with closed-cell polyisocyanurate foam to help both with stabilization and crush response. The ultimate goal of the test program was to design a specimen capable of limiting the sustained crush acceleration to 20 g through a stable crush progression. After a series of material tests, these specimens were evaluated under both static and dynamic conditions for impact energy attenuation characteristics and crush stability. Additionally, a series of simulation models were developed in parallel to the test efforts. It is anticipated that the models developed using the component level test efforts can be used to help guide the development of a design for use in full-scale eVTOL vehicle applications.

Littell, Justin↗

Overview of the U.S. DOE advanced fuels campaign

Advanced Fuel Campaign (AFC) has responsibility to develop advanced nuclear fuel technologies for the Fuel Cycle Research and Development (FCRD) program. AFC focus areas include 4 topics: advanced LWR fuels, actinide transmutation, research infrastructure development for the study of new nuclear fuels and advanced modeling and simulation. This paper reviews the latest technical advances in the 4 topics. Fuel and cladding materials under investigation include the following general classes of technologies: iron-based cladding, ceramic cladding, coated Zirconium alloy cladding, molybdenum-based cladding, high-density fuel, fully ceramic matrix fuel, and metallic alloy-based nuclear fuel. Concerning actinide transmutation fuel alloy additions to immobilize the lanthanide fission products, sodium-free annular fuel concepts, and cladding coatings and liners for mitigation of fuel-cladding chemical interaction have been investigated. Direct casting of metallic fuel into a zirconium liner was demonstrated. A Talos Transmission Electron Microscope was acquired for the study of advanced nuclear fuels and materials. The conceptual and preliminary designs for a static capsule, a flowing water loop, and a flowing sodium loop to support future transient testing needs have been completed. AFC interfaces with the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to develop multi-scale, multi-physics fuel performance codes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Mechanical Response of Porcine Liver Tissue under High Strain Rate Compression

In automobile accidents, abdominal injuries are often life-threatening yet not apparent at the time of initial injury. The liver is the most commonly injured abdominal organ from this type of trauma. In contrast to current safety tests involving crash dummies, a more detailed, efficient approach to predict the risk of human injuries is computational modelling and simulations. Further, the development of accurate computational human models requires knowledge of the mechanical properties of tissues in various stress states, especially in high-impact scenarios. In this study, a polymeric split-Hopkinson pressure bar (PSHPB) was utilized to apply various high strain rates to porcine liver tissue to investigate its material behavior during high strain rate compression. Liver tissues were subjected to high strain rate impacts at 350, 550, 1000, and 1550 s –1 . Tissue directional dependency was also explored by PSHPB testing along three orthogonal directions of liver at a strain rate of 350 s –1 . Histology of samples from each of the three directions was performed to examine the structural properties of porcine liver. Porcine liver tissue showed an inelastic and strain rate-sensitive response at high strain rates. The liver tissue was found lacking directional dependency, which could be explained by the isotropic microstructure observed after staining and imaging. Furthermore, finite element analysis (FEA) of the PSHPB tests revealed the stress profile inside liver tissue and served as a validation of PSHPB methodology. The present findings can assist in the development of more accurate computational models of liver tissue at high-rate impact conditions allowing for understanding of subfailure and failure mechanisms.

60 APPLIED LIFE SCIENCES↗

Probabilistic Modeling Of Ocular Biomechanics In VIIP: Risk Stratification

Visual Impairment and Intracranial Pressure (VIIP) syndrome is a major health concern for long-duration space missions. Currently, it is thought that a cephalad fluid shift in microgravity causes elevated intracranial pressure (ICP) that is transmitted along the optic nerve sheath (ONS). We hypothesize that this in turn leads to alteration and remodeling of connective tissue in the posterior eye which impacts vision. Finite element (FE) analysis is a powerful tool for examining the effects of mechanical loads in complex geometries. Our goal is to build a FE analysis framework to understand the response of the lamina cribrosa and optic nerve head to elevations in ICP in VIIP. To simulate the effects of different pressures on tissues in the posterior eye, we developed a geometric model of the posterior eye and optic nerve sheath and used a Latin hypercubepartial rank correlation coef-ficient (LHSPRCC) approach to assess the influence of uncertainty in our input parameters (i.e. pressures and material properties) on the peak strains within the retina, lamina cribrosa and optic nerve. The LHSPRCC approach was repeated for three relevant ICP ranges, corresponding to upright and supine posture on earth, and microgravity [1]. At each ICP condition we used intraocular pressure (IOP) and mean arterial pressure (MAP) measurements of in-flight astronauts provided by Lifetime Surveillance of Astronaut Health Program, NASA Johnson Space Center. The lamina cribrosa, optic nerve, retinal vessel and retina were modeled as linear-elastic materials, while other tissues were modeled as a Mooney-Rivlin solid (representing ground substance, stiffness parameter c1) with embedded collagen fibers (stiffness parameters c3, c4 and c5). Geometry creationmesh generation was done in Gmsh [2], while FEBio was used for all FE simulations [3]. The LHSPRCC approach resulted in correlation coefficients in the range of 1. To assess the relative influence of the uncertainty in an input parameter on the peak strains, we ranked and then normalized these coefficients, considering that normalized values 0.5 implied a substantial influence on the range of the peak strains in the optic nerve head (ONH). IOP and ICP were found to have a major influence on the peak strains in the ONH, as did optic nerve and LC stiffness. Interestingly, the stiffness of the sclera far from the scleral canal did not have a large influence on peak strains in ONH tissues; however, the collagen fiber stiffness in the peripapillary sclera and annular ring both influenced the peak strains within the ONH. We have created a physiologically relevant model that incorporated collagen fibers to study the effects of elevated ICP. Elevated ICP resulted in strains in the optic nerve that are not predicted to occur on earth: the upright or supine conditions. We found that IOP, ICP, lamina cribrosa stiffness and optic nerve stiffness had the highest association with these extreme strains in the ONH. These extreme strains may activate mechanosensitive cells that induce tissue remodeling and are a risk factor for the development of VIIP.

biomechanics↗

Dynamic Strength of Copper at High Pressures Using Pressure Shear Plate Experiments

Measurement of strength at high pressures is critical to understanding the behavior of materials subjected to extreme loading conditions. Recent developments have extended the pressure shear plate impact (PSPI) technique to the high-pressure regime. Such experiments provide a unique opportunity to extract the complete stress-strain behavior of materials at relatively high pressures. The modified technique includes a new fiber-optic heterodyne transverse velocity interferometer system and new analysis methods to account for the inelastic response of the anvil plates. In this study, PSPI experiments are conducted on oxygen-free high conductivity (OFHC) copper at pressures ranging from 10 to 43 GPa and at strain rates of ~10 5 s -1 . Complete stress-strain curves of copper are obtained using a hybrid methodology that relies on numerical simulations and particle velocity records gathered in these experiments. Experiments are conducted at similar pressures and strain rates using two types of anvils that undergo different levels of plastic strains: tool steel and tungsten carbide, to check the robustness of the hybrid method. The behavior of copper at these high pressures and strain rates is compared with previous literature data, and the effect of pressure and strain rate on the strength behavior of copper are discussed. Here, it was observed that the rate of change in yield strength change in pressure was 2.5 times the rate of change of shear modulus with pressure. The strength models that predict the response of copper at these extreme conditions are also examined and validated.

36 MATERIALS SCIENCE↗

Mechanisms of local stress sensing in multifunctional polymer films using fluorescent tetrapod nanocrystals

Nanoscale stress-sensing can be used across fields ranging from detection of incipient cracks in structural mechanics to monitoring forces in biological tissues. We demonstrate how tetrapod quantum dots (tQDs) embedded in block-copolymers act as sensors of tensile/compressive stress. Remarkably, tQDs can detect their own composite dispersion and mechanical properties, with a switch in optomechanical response when tQDs are in direct contact. Using experimental characterizations, atomistic simulations and finite-element analyses, we show that under tensile stress, densely-packed tQDs exhibit a photoluminescence peak shifted to higher energies (“blue-shift”) due to volumetric compressive stress in their core; loosely-packed tQDs exhibit a peak shifted to lower energies (“red-shift”) from tensile stress in the core. The stress-shifts result from the tQD's unique branched morphology in which the CdS arms act as antennas that amplify the stress in the CdSe core. Our nanocomposites exhibit excellent cyclability and scalability with no degraded properties of the host polymer. Colloidal tQDs allow sensing in many materials to potentially enable auto-responsive, smart structural nanocomposites that self-predict impending fracture.

36 MATERIALS SCIENCE↗

Damage initiation and evolution in Al-Si layered microstructures under shock loading conditions at atomic scales

Designing materials with microstructural features like multiphase interfaces have shown significant promise for usage in the next generation defense and nuclear applications. The dynamic response of these interfaces in extreme environments is observed to vary with the deformability of individual phases, which could subsequently alter the favored damage nucleation sites related to spall failure. Here, we investigate the role of spacing of interfaces in a nanocrystalline layered Aluminum-Silicon system on the shock wave propagation behavior, and microstructural and defect evolution using classical molecular dynamics simulations. The simulations are carried out with different hypothetical Al/Si microstructures including variations in the distribution of Si grains as layers in a nanocrystalline Al matrix. The molecular dynamics simulations suggest that spall failure is preferentially initiated at Al/Si interfaces. In addition, for the same system volume and same concentration of Si, a higher number of layers is marked by an increase in shock wave velocity and reduced resistance to spall failure.

36 MATERIALS SCIENCE↗

Material processing of convection-driven flow field and temperature distribution under oblique gravity

A set of mathematical formulation is adopted to study vapor deposition from source materials driven by heat transfer process under normal and oblique directions of gravitational acceleration with extremely low pressure environment of 10(exp -2) mm Hg. A series of time animation of the initiation and development of flow and temperature profiles during the course of vapor deposition has been obtained through the numerical computation. Computations show that the process of vapor deposition has been accomplished by the transfer of vapor through a fairly complicated flow pattern of recirculation under normal direction gravitational acceleration. It is obvious that there is no way to produce a homogeneous thin crystalline films with fine grains under such a complicated flow pattern of recirculation with a non-uniform temperature distribution under normal direction gravitational acceleration. There is no vapor deposition due to a stably stratified medium without convection for reverse normal direction gravitational acceleration. Vapor deposition under oblique direction gravitational acceleration introduces a reduced gravitational acceleration in vertical direction which is favorable to produce a homogeneous thin crystalline films. However, oblique direction gravitational acceleration also induces an unfavorable gravitational acceleration along horizontal direction which is responsible to initiate a complicated flow pattern of recirculation. In other words, it is necessary to carry out vapor deposition under a reduced gravity in the future space shuttle experiments with extremely low pressure environment to process vapor deposition with a homogeneous crystalline films with fine grains. Fluid mechanics simulation can be used as a tool to suggest most optimistic way of experiment with best setup to achieve the goal of processing best nonlinear optical materials.

Hung, R. J.↗

Sweep-tracing algorithm: in silico slip crystallography and tension-compression asymmetry in BCC metals

Abstract Direct Molecular Dynamics (MD) simulations are being increasingly employed to model dislocation-mediated crystal plasticity with atomic resolution. Thanks to the dislocation extraction algorithm (DXA), dislocation lines can be now accurately detected and positioned in space and their Burgers vector unambiguously identified in silico, while the simulation is being performed. However, DXA extracts static snapshots of dislocation configurations that by themselves present no information on dislocation motion. Referred to as a sweep-tracing algorithm (STA), here we introduce a practical computational method to observe dislocation motion and to accurately quantify its important characteristics such as preferential slip planes (slip crystallography). STA reconnects pairs of successive snapshots extracted by DXA and computes elementary slip facets thus precisely tracing the motion of dislocation segments from one snapshot to the next. As a testbed for our new method, we apply STA to the analysis of dislocation motion in large-scale MD simulations of single crystal plasticity in BCC metals. We observe that, when the crystal is subjected to uniaxial deformation along its [001] axis, dislocation slip predominantly occurs on the {112} maximum resolved shear stress plane under tension, while in compression slip is non-crystallographic (pencil) resulting in asymmetric mechanical response. The marked contrast in the observed slip crystallography is attributed to the twinning/anti-twinning asymmetry of shears in the {112} planes relatively favoring dislocation motion in the twinning sense while hindering dislocations from moving in the anti-twinning directions.

36 MATERIALS SCIENCE↗

Approach and rotation of reconnecting topological defect lines in liquid crystal

Topological defects are a universal concept across many disciplines, such as crystallography, liquid-crystalline physics, low-temperature physics, cosmology, and even biology. In nematic liquid crystals, topological defects called disclinations have been widely studied. For their three-dimensional (3D) dynamics, however, only recently have theoretical approaches dealing with fully 3D configurations been reported. Further, recent experiments have observed 3D disclination line reconnections, a phenomenon characteristic of defect line dynamics, but detailed discussions were limited to the case of approximately parallel defects. In this paper, we focus on the case of two disclination lines that approach at finite angles and lie in separate planes, a more fundamentally 3D reconnection configuration. Observing and analyzing such reconnection events, we find the square-root law of the distance between disclinations and the decrease of the interdisclination angle over time. We compare the experimental results with theory and find qualitative agreement on the scaling of distance and angle with time, but quantitative disagreement on distance and angle relative mobilities. To probe this disagreement, we derive the equations of motion for systems with reduced twist constant and also carry out simulations for this case. These, together with the experimental results, suggest that deformations of disclinations may be responsible for the disagreement. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗

Summary of Grizzly Development for Advanced Reactor Structural Materials

This summary report contains an overview of work performed during Fiscal Year 2019 under the Nuclear Engineering Advanced Modeling and Simulation (NEAMS) Fuels Product Line work package entitled "MS-19IN020106 - GRIZZLY - INL". Grizzly is a code based on the MOOSE framework that has been under development for several years funded by the US Department of Energy’s Light Water Reactor Sustainability (LWRS) program. Prior to 2019, Grizzly development was focused on development of capabilities needed to simulate the effects of aging mechanisms in light water reactor (LWR) components. The goal of this LWR-focused work has been to develop tools that can be used to assess the ability of existing LWRs to perform safety in long term operation (LTO) scenarios. This LWR-focused work has primarily focused on the effects of reactor pressure vessel (RPV) embrittlement and aging of reinforced concrete structures. Grizzly development to support LWR degradation needs is still ongoing, but starting in 2019, Grizzly’s capabilities are being expanded to also meet growing needs for a tool to address structural component integrity and degradation in advanced reactors. A simulation capability for advanced reactor components has many of the same basic requirements that one for LWRs does, such as the ability to solve coupled systems of partial differential equations that arise from the fact that many of the problems of interest involve multiple coupled physics. The materials issues also inherently involve multiple length scales, with the response at the engineering scale being strongly affected by mechanisms occurring at atomistic and mesoscopic scales.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗