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 757 records · Page 42

Rotational memory function of SPC/E water

Memory effects are essential for the dynamics of condensed materials and are responsible for non-exponential relaxation of correlation functions of dynamic variables through the memory function. Memory functions of dipole rotations for water have never been calculated directly from molecular dynamics simulations. We present here calculations of memory functions for single-dipole rotations and for the overall dipole moment of the sample for SPC/E water. The normalized memory functions for single-particle and collective dipole dynamics turn out to be nearly identical. This result validates theories of dielectric spectroscopy in terms of single-particle time correlation functions and the connection between the collective and single-particle relaxation times through the Kirkwood factor. The dielectric function in this formalism contains no new dynamic information that does not exist in the single-dipole correlation function. A short memory time, ≲1 fs, justifies the use of the mathematics of rotational diffusion to describe the dynamics of a single molecular dipole moment in bulk water. Here, an analytical equation for the rotational memory time is derived.

Asthagiri, Dilipkumar N. [Oak Ridge National Labor↗

Directed aging, memory, and nature’s greed

Disordered materials are often out of equilibrium and evolve very slowly in a rugged and tortuous energy landscape. This slow evolution, referred to as aging, is deemed undesirable as it often leads to material degradation. However, we show that aging also encodes a memory of the stresses imposed during preparation. Because of inhomogeneous local stresses, the material itself decides how to evolve by modifying stressed regions differently from those under less stress. Because material evolution occurs in response to stresses, aging can be “directed” to produce sought-after responses and unusual functionalities that do not inherently exist. Aging obeys a natural “greedy algorithm” as, at each instant, the material simply follows the path of most rapid and accessible relaxation. Our experiments and simulations illustrate directed aging in examples in which the material’s elasticity transforms as desired because of an imposed deformation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Exascale Dislocation Simulator

ExaDiS, or Exascale Dislocation Simulator, is a set of software modules written to enable numerical simulations of large groups of moving and interacting dislocations, line defects in crystals responsible for crystal plasticity. By tracking time evolution of sufficiently large dislocation ensembles, ExaDiS predicts plasticity response and plastic strength of crystalline materials. ExaDiS is built around a library of core functions for Discrete Dislocation Dynamics method specifically written to perform efficient computations on GPUs.

Bertin, Nicolas↗

From DNA to FBA: How to Build Your Own Genome-Scale Metabolic Model

Microbiological studies are increasingly relying on in silico methods to perform exploration and rapid analysis of genomic data, and functional genomics studies are supplemented by the new perspectives that genome-scale metabolic models offer. A mathematical model consisting of a microbe’s entire metabolic map can be rapidly determined from whole-genome sequencing and annotating the genomic material encoded in its DNA. Flux-balance analysis (FBA), a linear programming technique that uses metabolic models to predict the phenotypic responses imposed by environmental elements and factors, is the leading method to simulate and manipulate cellular growth in silico. However, the process of creating an accurate model to use in FBA consists of a series of steps involving a multitude of connections between bioinformatics databases, enzyme resources, and metabolic pathways. We present the methodology and procedure to obtain a metabolic model using PyFBA, an extensible Pythonbased open-source software package aimed to provide a platform where functional annotations are used to build metabolic models (http://linsalrob.github.io/PyFBA). Backed by the Model SEED biochemistry database, PyFBA contains methods to reconstruct a microbe’s metabolic map, run FBA upon different media conditions, and gap-fill its metabolism. The extensibility of PyFBA facilitates novel techniques in creating accurate genome-scale metabolic models.

genome annotation↗

Reliability assessment of thrust chamber cooling concepts using probabilistic analysis techniques

The reliability of OFHC (Oxygen Free High Conductivity) copper and NARloy-Z thrust chambers is assessed by applying probabilistic structural analysis techniques to incorporate design parameter variability and uncertainty. Thrust chambers specifically evaluated are the cylindrical test fixtures employed in a plug-nozzle configuration at the NASA Lewis Research Center. Direct sampling Monte Carlo simulations based on a simplified life prediction methodology established probability densities of firing cycles to structural failure. Simulated cyclic lives demonstrated modest agreement to experiment. Similarly, regions of high structural failure probability were determined using a limit state approach employing calculated cumulative distribution functions for effective stress response and an assumed material strength distribution. A probability of failure of 0.012 was calculated at the center of the coolant channel hot-gas-side wall for an OFHC milled channel. Structural response was found to be sensitive to the uncertainties in the thrust chamber thermal environment and the material's thermal expansion coefficient.

Rapp, Douglas C.↗

Reliability assessment of thrust chamber cooling concepts using probabilistic analysis techniques

The reliability of OFHC (Oxygen Free High Conductivity) copper and NARloy-Z thrust chambers is assessed by applying probabilistic structural analysis techniques to incorporate design parameter variability and uncertainty. Thrust chambers specifically evaluated are the cylindrical test fixtures employed in a plug-nozzle configuration at the NASA Lewis Research Center. Direct sampling Monte Carlo simulations based on a simplified life prediction methodology established probability densities of firing cycles to structural failure. Simulated cyclic lives demonstrated modest agreement to experiment. Similarly, regions of high structural failure probability were determined using a limit state approach employing calculated cumulative distribution functions for effective stress response and an assumed material strength distribution. A probability of failure of 0.012 was calculated at the center of the coolant channel hot-gas-side wall for an OFHC milled channel. Structural response was found to be sensitive to the uncertainties in the thrust chamber thermal environment and the material's thermal expansion coefficient.

Rapp, Douglas C.↗

Environmental effects on long term behavior of composite laminates

Model equations are presented for approximate methods simulating the long-term behavior of composite materials and structures in hot/humid service environments. These equations allow laminate property upgradings with time, and can account for the effects of service environments on creep response. These methodologies are illustrated for various individual and coupled temperature/moisture, longitudinal/transverse, and composite material type cases. Creep deformation is noted to rise dramatically for cases of matrix-borne, but not of fiber-borne, loading in hot, humid environments; the coupled influence of temperature and moisture is greater than a mere combination of their individual influences.

Singhal, S. N.↗

Active Turbine Tip Clearance Control Trade Space Analysis of an Advanced Geared Turbofan Engine

Tip clearance within the high pressure turbine of a gas turbine engine is a significant factor in engine performance and efficiency. In the pursuit of higher efficiency, aero-engine designs are migrating toward compact gas turbine (CGT) technology that seeks to increase the bypass ratio of the gas turbine engine without increasing the size of the fan, which is constrained by its underwing location. The reduced size of CGTs invoke concern over increased sensitivity of engine performance due to turbine tip clearance gap that makes an argument for advanced tip clearance mitigation and control techniques to be employed. This paper evaluates the tip clearance trade space for a conceptual geared turbofan engine with a CGT core. This is accomplished through a modeling and simulation approach that includes a sensitivity analysis of engine performance in response to high pressure turbine tip clearance as well as an evaluation of the sensitivity of tip clearance to various design parameters, including material properties and component cooling characteristics. Also included is a parametric study of actuators that provides preliminary requirements for implementation of active turbine tip clearance control actuation systems. The results produced from these studies are meant to be informative, with special emphasis on the demonstration of a systematic approach. The modeling approach appears to capture expected trends. The studies suggest that the tip clearance gap will have a greater impact on the new CGT engines and that a relatively slow, actively controlled actuation system may be sufficient as long as it has control authority to both open and close the tip clearance gap.

tip clearance control↗

Pulsed Sphere Time of Flight Neutron Spectrum Measurements: Blank Run and 1.8 Mean-Free-Path thick Polyethylene Moderating Sphere

From the late 1960s to 1985, Lawrence Livermore National Laboratory (LLNL) undertook a series of ‘pulsed sphere’ experiments for 32 materials involving 148 different experiments using 75 different spheres to measure the neutron leakage spectra for different materials from 14-MeV neutrons generated by 3 H(d, n) 4 He reactions induced by an incident D+ beam from the Insulated Core Transformer (ICT) accelerator at LLNL. This deuteron beam was focused to impinge on tritium loaded onto a titanium substrate within a low mass target assembly placed within the center of a spherical shell constructed from the material of interest. The neutron time of flight was measured at a detector outside a collimator. The program was intended to be comprehensive with numerous targets of various thicknesses spanning the periodic table from H to 239 Pu. Five different neutron detectors were utilized to avoid possible systematic errors due to differences in detector responses. Further, these measurements were designed with the aim that the simple geometry of the measurements could be easily simulated to validate Monte Carlo transport codes and nuclear data libraries. As a result, these benchmarks have been extensively used in the validation of ENDL, ENDF and other libraries. Two experiments selected as benchmark cases in this evaluation are: 1) target assembly without material sphere (blank run), and 2) target assembly with a hollow polyethylene sphere, 1.8 mean free path (mfp) or 16.56 cm thick.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Characterization of modeling and experimental data inconsistencies from burst testing for high-burnup commercial fuel rod applications

Fuel fragmentation, relocation, and dispersal are some of the largest issues remaining in the nuclear industry before rod-average burnup can be increased beyond 62 GWd/tU. The issue is primarily related to the potential for fuel to be dispersed into the reactor primary system, which may increase public risk. One way to prevent dispersal is to avoid cladding burst. The objective of this work is to support the high burnup safety case by evaluating cladding burst under high-burnup, full-length fuel rods and to identify uncertainties that could improve model predictions. The results of this analysis will evaluate realistic, prototypic loss-of-coolant accident (LOCA) conditions; support future cladding burst test designs; and inform the development of mechanistic material models. Realistic high-burnup operating conditions were implemented in the BISON fuel performance code to simulate steady-state and LOCA transient fuel rod evolution to the point at which cladding burst occurred. Parametric studies are performed to assess code response to changes in rod internal pressures and heating rates. Results were compared with simulated LOCA experiments to identify inconsistencies between commercial fuel rod analysis and experimental validation. The representative full-length fuel rod LOCA simulation results did not agree with cladding burst tests. Cladding burst tests indicated burst occurring well below (100–150 °C) those calculated in the full-length fuel rod LOCA analysis. Further investigation indicated that the cladding burst tests do not appear to be representative for full-length fuel rods. The inconsistency investigated in this work showed that the differences between the BISON simulation and the experiment's cladding burst conditions arise from an incomplete characterization of the cladding surface temperature, detailed rodlet characterization, lack of cladding strain measurements, and uncertainty in the cladding creep and failure models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Experimental verification of a crystal plasticity-based simulation framework for predicting microstructure and geometric shape changes: Application to bending and Taylor impact testing of Zr

This paper is concerned with experimental verification of a recently developed multi-scale simulation framework for plastic deformation of metallic materials from quasi-static to impact deformation conditions. The framework is a visco-plastic self-consistent (VPSC) polycrystalline model embedded in an implicit finite element method (FE-VPSC) to provide a microstructure-sensitive constitutive response at each material point. Each material point of the FEM model is a polycrystalline aggregate with crystallographic deformation mechanisms operating at the single crystal scale with their evolving activity based on a dislocation density-based hardening law and texture. Four beams and three cylinders machined in different orientations from a textured plate of high-purity zirconium are tested quasi-statically in 4-point bending and at speeds of 100 m/s, 170 m/s and 243 m/s during Taylor impact tests, respectively. The variation in dimensional changes resulting from different sample orientations in the plate with respect to loading directions is measured for each sample. Moreover, texture and twinning characterization is performed using electron backscattered diffraction (EBSD). The deformation processes and underlying evolution of microstructure are successfully simulated using the FE-VPSC framework. In doing so, the model parameters are optimized and validated across a broad range of strain rates and temperatures. Simulation results in terms of geometrical changes and microstructural evolution are compared with the experimental measurements. The model predicts anisotropic material flow resulting from the hard-to-deform crystallographic directions, the development of gradients in texture and twinning through the geometries, tension–compression asymmetry, as well as the extent of plasticity under impact.

42 ENGINEERING↗

Nanoindentation Stress–Strain for Fracture Analysis and Computational Modeling for Hardness and Modulus

Stress–Strain plots based on nanoindentation load–depth curves were obtained to study phenomena like internal fracture and ductile to brittle transitions. Fracture phenomena during the indentation process were analyzed based on the stress–strain plots. A transition from ductile to brittle fracture was observed on increasing the depth or load of indenter penetration. A new approach with shape factors in the fracture studies based on radial crack branching and micro-cracking was done. Hardness and modulus plots were fitted with polynomials. The fitting parameters were varied to obtain different hardness and modulus responses.

36 MATERIALS SCIENCE↗

Predicting forming limit diagrams for magnesium alloys using crystal plasticity finite elements

The crystal plasticity finite element (CPFE) method was used to predict forming limit strains of hexagonal close-packed (HCP) polycrystalline magnesium alloy sheets, namely AZ31 and ZE10, under an isothermal temperature condition. The strain rate-dependent uniaxial tensile test data along various loading directions and an initial texture were used to determine the constitutive parameters of the crystal plasticity model. A hybrid representative volume element approach, which combines the CPFE and the Marciniak–Kuczynski model, was developed to obtain the forming limit diagram of the magnesium alloys. The predicted forming limits were in excellent agreement with the Nakazima test results. Moreover, the microscopic responses such as slip/twin activation and deformation texture changes during various loading paths from uniaxial tension to the balanced biaxial tension were comprehensively analyzed and discussed from the CPFE results to understand the underlying micro-mechanism for the macro-mechanical responses.

36 MATERIALS SCIENCE↗

Defect identification in simulated Bragg coherent diffraction imaging by automated AI

X-ray Bragg coherent diffraction imaging is a powerful technique for operando and in situ materials characterization and provides a unique means of quantifying the influence of one-dimensional (1D) and two-dimensional (2D) material defects on material response. However, obtaining full images from raw x-ray diffraction data is nontrivial and computationally intensive, precluding real-time experimental feedback. Here, we present a machine learning approach to identify the presence of crystalline line defects (edge and screw) in samples from the raw, 2D, coherent diffraction data without the need for image reconstruction through iterative phase retrieval. Further, we compare different approaches to designing neural networks for this application and demonstrate the potential of automated ML (autoML) approaches.

36 MATERIALS SCIENCE↗

Mechanical-Electrochemical-Thermal Simulation of Lithium-Ion Cells

Models capture the force response for single-cell and cell-string levels to within 15%-20% accuracy and predict the location for the origin of failure based on the deformation data from the experiments. At the module level, there is some discrepancy due to poor mechanical characterization of the packaging material between the cells. The thermal response (location and value of maximum temperature) agrees qualitatively with experimental data. In general, the X-plane results agree with model predictions to within 20% (pending faulty thermocouples, etc.); the Z-plane results show a bigger variability both between the models and test-results, as well as among multiple repeats of the tests. The models are able to capture the timing and sequence in voltage drop observed in the multi-cell experiments; the shapes of the current and temperature profiles need more work to better characterize propagation. The cells within packaging experience about 60% less force under identical impact test conditions, so the packaging on the test articles is robust. However, under slow-crush simulations, the maximum deformation of the cell strings with packaging is about twice that of cell strings without packaging.

25 ENERGY STORAGE↗

Numerical investigation of laser-driven shock interaction with a deformable particle

A laser-driven shock propagating through an isolated particle embedded in a plastic (CH) target was studied using the radiation-hydrodynamic code FLASH. Preliminary simulations using IONMIX equations of state (EOS) showed significant differences in the shock Hugoniot of aluminum compared to experimental data in the low-pressure regime [O(10) GPa], resulting in higher streamwise compression and deformation of an aluminum particle. Hence, a simple modification to the ideal gas EOS was developed and employed to describe the target materials and examine the particle dynamics. The evolution of the pressure field demonstrated a complex wave interaction, resulting in a highly unsteady particle drag which featured two drag minima due to shock focusing at the rear end of the particle and rarefaction stretching due to laser shut-off. Although ~30% lateral expansion and ~25% streamwise compression were observed, the aluminum particle maintained considerable integrity without significant distortion. Additional simulations examined the particle response for a range of particle densities, sizes, and acoustic impedances. The data revealed that lighter particles such as aluminum gained significant momentum, reaching up to ~96% of the shocked CH's speed, compared to ~29% for the heavier tungsten particles. Despite the differences seen in the early stage of shock interaction, particles with varying acoustic impedances ultimately reached the same peak velocity. This identified particle-to-host density ratio is an important factor in determining the inviscid terminal velocity of the particle. In addition, the modified EOS model presented in this study could be used to approximate solid materials in hydrocodes that lack material strength models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Utilizing computer vision and artificial intelligence algorithms to predict and design the mechanical compression response of direct ink write 3D printed foam replacement structures

Additive Manufacturing (AM) of porous polymeric materials, such as foams, recently became a topic of intensive research due their unique combination of low density, impressive mechanical properties, and stress dissipation capabilities. Conventional methods for fabricating foams rely on complex and stochastic processes, making it challenging to achieve precise architectural control of structured porosity. In contrast, AM provides access to a wide range of printable materials, where precise spatial control over structured porosity can be modulated during the fabrication process enabling the production of foam replacement structures (FRS). Current approaches for designing FRS are based on intuitive understanding of their properties or an extensive number of finite element method (FEM) simulations. These approaches, however, are computationally expensive and time consuming. As such, in this work, we present a novel methodology for determining the mechanical compression response of direct ink write (DIW) 3D printed FRS using a simple cross-sectional image. By obtaining measurement data for a relatively small number of samples, an artificial neural network (ANN) was trained, and a computer vision algorithm was used to make inferences about foam compression characteristics from a single cross-sectional image. Finally, a genetic algorithm (GA) was used to solve the inverse design problem, generating the AM printing parameters that an engineer should use to achieve a desired compression response from a DIW printed FRS. The methods developed herein present an avenue for entirely autonomous design and analysis of additively manufactured structures using artificial intelligence.

36 MATERIALS SCIENCE↗

Defect Spirograph: Dynamical Behavior of Defects in Spatially Patterned Active Nematics

We report that topological defects in active liquid crystals can be confined by introducing gradients of activity. Here, we examine the dynamical behavior of two defects confined by a sharp gradient of activity that separates an active circular region and a surrounding passive nematic material. Continuum simulations are used to explain how the interplay among energy injection into the system, hydrodynamic interactions, and frictional forces governs the dynamics of topologically required self-propelling + 1/2 defects. Our findings are rationalized in terms of a phase diagram for the dynamical response of defects in terms of activity and frictional damping strength. Different regions of the underlying phase diagram correspond to distinct dynamical modes, namely immobile defects, steady rotation of defects, bouncing defects, bouncing-cruising defects, dancing defects, and multiple defects with irregular dynamics. These dynamic states raise the prospect of generating synchronized defect arrays for microfluidic applications.

36 MATERIALS SCIENCE↗