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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

FY23 Simulation of Elastic-Plastic Failure Propagation

This memo summarizes the simulation of ductile failure propagation work conducted under the ASC project “V&V of Ductile Failure” conducted during FY 23. Physically, the failure propagation consists of crack propagation in the material. In the numerical setting—specifically in a finite element model—propagation can be accomplished through element death when critical conditions occur locally at an element that is then deleted from the simulation. The validation of the finite element models is evaluated by direct comparison between the experimental and simulation results regarding the rate of crack growth and its influence on the load-deflection response of the specimens tested. This work considers two geometries that display stable crack propagation under displacement-controlled conditions. The first geometry consists of hat specimens loaded in compression with nominally identical geometries but made with three different materials: Steel A286, Al 7075-T651 and 304L stainless steel. The three materials represent a range of ductility values that affect the response and crack propagation within the specimen. The crack induced propagates under an essentially mode-II type of deformation. The second geometry consists of a pre-cracked 304L stainless steel compact tension test specimen loaded so as to induce a mode-I deformation at the crack.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Sensitivity analysis of high resolution gamma-ray detection for safeguards monitoring at natural uranium conversion facilities

Under the policies proposed by recent International Atomic Energy Agency (IAEA) circulars and policy papers, implementation of safeguards exists when any purified aqueous uranium solution or uranium oxides suitable for isotopic enrichment or fuel fabrication exists. Under IAEA Policy Paper 18, the starting point for nuclear material under safeguards was reinterpreted, suggesting that purified uranium compounds should be subject to safeguards procedures no later than the first point in the conversion process. In response to this technical need, a combination of simulation models and experimental measurements were employed in previous work to develop and validate gamma-ray nondestructive assay monitoring systems in a natural uranium conversion plant (NUCP). In particular, uranyl nitrate (UO 2 (NO 3 ) 2 ) solution exiting solvent extraction was identified as a key measurement point (KMP). Passive nondestructive assay techniques using high resolution gamma-ray spectroscopy were evaluated to determine their viability as a technical means for drawing safeguards conclusions at NUCPs, and if the IAEA detection requirements of 1 significant quantity (SQ) can be met in a timely manner. Building upon the aforementioned previous validation work on detector sensitivity to varying concentrations of uranyl nitrate via a series of dilution measurements, this work investigates detector response parameter sensitivities to gamma-ray signatures of uranyl nitrate. The full energy peak efficiency of a detection system is dependent upon the sample, geometry, absorption, and intrinsic efficiency parameters. Perturbation of these parameters translates into corresponding variations of the 185.7 keV peak area of the 235 U in uranyl nitrate. Such perturbations in the assayed signature impact the quality or versatility of the safeguards conclusions drawn. Given the potentially high throughput of uranyl nitrate in NUCPs, the ability to assay 1 SQ of material requires uncertainty «1%. Accounting for material self-shielding properties, pipe thickness, and source-detector orientation is instrumental in determining the robustness of gamma-ray detection in the process monitoring of uranyl nitrate in NUCPs. Monte Carlo models and ray-tracing models were employed to determine the sensitivity of the detected 185.7 keV photon to self-shielding properties, pipe thickness, and source-detector geometry. Considering the implementation of the detection of 1 SQ, diversion of 1 SQ becomes essentially undetectable given the systematic uncertainty, in addition to considerations such as propagating uncertainties due to pipe offset/position, as well as minor variations in pipe thickness. Consequently, pipe thickness was the most sensitive variable in affecting full energy efficiency of the 185.7 keV signature peak with up to 8% variation in efficiency for ±0.5 mm changes in Schedule 40 304L stainless steel piping. Furthermore, computation of the attenuation correction factor of the uranyl nitrate solution [CF(AT) (i.e. εsample)] using Parker's method using with the approximation for the geometrical factor κ≈π/4 was validated through experimental, Monte Carlo and ray-tracing calculations for a uranyl nitrate filled transfer pipe segment. Furthermore, quantifying sensitivity in detector position, as well as voiding effects due to bubbly flow or laminar flow with an air gap in the uranyl nitrate becomes increasingly important as considerations from (static) design-scale measurements translate into (dynamic) field operations tests.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2017

Efficient and Extensible Quasi-Explicit Modular Nonlinear Multiscale Battery Model: GH-MSMD

Complex physics and long computation time hinder the adoption of computer aided engineering models in the design of large-format battery cells and systems. A modular, efficient battery simulation model -- the multiscale multidomain (MSMD) model -- was previously introduced to aid the scale-up of Li-ion material and electrode designs to complete cell and pack designs, capturing electrochemical interplay with 3-D electronic current pathways and thermal response. Here, this paper enhances the computational efficiency of the MSMD model using a separation of time-scales principle to decompose model field variables. The decomposition provides a quasi-explicit linkage between the multiple length-scale domains and thus reduces time-consuming nested iteration when solving model equations across multiple domains. In addition to particle-, electrode- and cell-length scales treated in the previous work, the present formulation extends to bus bar- and multi-cell module-length scales. We provide example simulations for several variants of GH electrode-domain models.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2025

Fluoride-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model Updates

This work presents the latest improvements to, and investigations performed with, the Fluoride-Cooled High-Temperature Pebble-Bed Reactor reference plant models for the United States Nuclear Regulatory Commission. These models, developed with the Comprehensive Reactor Analysis Bundle, or BlueCRAB, serve as the foundation for the future development of detailed design evaluation models based on license applications. BlueCRAB is the code suite proposed for non-light-water reactor systems safety analyses, and it incorporates various simulation tools developed by the Nuclear Energy Advance Modeling and Simulation program, including the Griffin code for reactor physics, the Pronghorn and SAM codes for core thermal fluids, the BISON code for solid conduction and fuel performance, and the SAM code for system analysis. The primary objective of this work is to assess the level of readiness of BlueCRAB for modeling fluoride-cooled high-temperature pebble-bed reactors. We first developed numerical models in BlueCRAB that include the key physics for this technology, ensuring an adequate level of fidelity for modeling the core performance during accident scenarios. This was followed by simulation of transient scenarios, two loss-of-forced-cooling events (one protected and one unprotected), and two control rod withdrawal events (one delayed and one prompt supercritical reactivity insertion). The analysis includes comparisons between the 2-D thermal fluid porous media models in Pronghorn and SAM, comparisons between coupled Pronghorn-Griffin and coupled SAM-Griffin models for two loss-of-forced cooling events and one control rod withdrawal event, and comparisons between SAM single-solve and domain-overlapping approaches for multi-scale thermal fluid coupling. In addition, we performed comparisons between 3-D, 2-D, and 0-D neutronic models for the two control rod withdrawal scenarios with Pronghorn-Griffin. The results show that the BlueCRAB models led to physically intuitive solutions for the scenarios examined. The changes in the various scalar and vector fields such as the neutron flux, power, temperatures, densities, pressures, and velocities are all within the expected ranges, and their distributions can be explained from the system response of the transients and the geometric and material variations. Several comparisons suggest that the porous media models in Pronghorn and SAM can lead to similar solutions, even though they are based on different methodologies. The simulations demonstrate that there are differences between the various levels of fidelity, and it is advisable to have flexible tools that can cover the breadth and depth of needs that may arise in future technical evaluations. We believe that BlueCRAB’s capabilities represent a significant asset for confirmatory analyses aimed at resolving important safety questions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

An efficient computational approach to model statistical correlations in photon counting x-ray detectors

Purpose: To introduce and evaluate an increment matrix approach (IMA) describing the signal statistics of energy-selective photon counting detectors including spatial–spectral correlations between energy bins of neighboring detector pixels. The importance of the occurring correlations for image-based material decomposition is studied. Methods: An IMA describing the counter increase patterns in a photon counting detector is proposed. This IMA has the potential to decrease the number of required random numbers compared to Monte Carlo simulations by pursuing an approach based on convolutions. To validate and demonstrate the IMA, an approximate semirealistic detector model is provided, simulating a photon counting detector in a simplified manner, e.g., by neglecting count rate-dependent effects. In this way, the spatial–spectral correlations on the detector level are obtained and fed into the IMA. The importance of these correlations in reconstructed energy bin images and the corresponding detector performance in image-based material decomposition is evaluated using a statistically optimal decomposition algorithm. Results: The results of IMA together with the semirealistic detector model were compared to other models and measurements using the spectral response and the energy bin sensitivity, finding a good agreement. Correlations between the different reconstructed energy bin images could be observed, and turned out to be of weak nature. These correlations were found to be not relevant in image-based material decomposition. An even simpler simulation procedure based on the energy bin sensitivity was tested instead and yielded similar results for the image-based material decomposition task, as long as the fact that one incident photon can increase multiple counters across neighboring detector pixels is taken into account. Conclusions: The IMA is computationally efficient as it required about 10{sup 2} random numbers per ray incident on a detector pixel instead of an estimated 10{sup 8} random numbers per ray as Monte Carlo approaches would need. The spatial–spectral correlations as described by IMA are not important for the studied image-based material decomposition task. Respecting the absolute photon counts and thus the multiple counter increases by a single x-ray photon, the same material decomposition performance could be obtained with a simpler detector description using the energy bin sensitivity.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Microscale modeling of effective mechanical and electrical properties of textiles

Here, a computational framework for assisting in the development of novel textiles is presented. Electronic textiles are key in the rapidly growing field of wearable electronics for both consumer and military uses. There are two main challenges to the modeling of electronic textiles: the discretization of the textile microstructure and the interaction between electromagnetic and mechanical fields. A director-based beam formulation with an assumed electrical current is used to discretize the fabric at the level of individual fibrils. The open-source package FEniCS was used to implement the finite element model. Contact integrals were added into the FEniCS framework so that multiphysics contact laws can be incorporated in the same framework, leveraging the code generation and automated differentiation capabilities of FEniCS to produce the tangents needed by the implicit solution method. The computational model is used to construct and determine the mechanical, thermal, and electrical properties of a representative volume elements of a plain woven textile. Dynamic relaxation to solve the mechanical fields and the electrical and thermal fields is solved statically for a given mechanical state. The simulated electrical responses are fit to a simplified Kirchhoff network model to determine effective resistances of the textile.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Micromechanical simulation of damage progression in carbon phenolic composites

Carbon/phenolic composites are used extensively as ablative insulating materials in the nozzle region of solid rocket motors. The current solid rocket motor (RSRM) on the space shuttle is fabricated from woven rayon cloth which is carbonized and then impregnated with the phenolic resin. These plies are layed up in the desired configuration and cured to form the finished part. During firing, the surface of the carbon/phenolic insulation is exposed to 5000 F gases from the rocket exhaust. The resin pyrolizes and the material chars to a depth which progresses with time. The rate of charring and erosion are generally predictable, and the insulation depth is designed to allow adequate safety margins over the firing time of the motor. However, anomalies in the properties and response of the carbon/phenolic materials can lead to severe material damage which may decrease safety margins to unacceptable levels. Three macro damage modes which were observed in fired nozzles are: ply lift, 'wedge out', and pocketing erosion. Ply lift occurs in materials with plies oriented nearly parallel to the surface. The damage occurs in a region below the charred material where material temperatures are relatively low - about 500 F. Wedge out occurs at the intersection of nozzle components whose plies are oriented at about 45 deg. The corner of the block of material breaks off along a ply interface. Pocketing erosion occurs in material with plies oriented normal to the surface. Thermal expansion is restrained in two directions resulting in large tensile strains and material failure normal to the surface. When a large section of material is removed as a result of damage, the insulation thickness is reduced which may lead to failure of the nozzle due to excessive heating of critical components. If these damage events cannot be prevented with certainty, the designer must increase the thickness of the insulator thus adding to both weight and cost. One of the difficulties in developing a full understanding of these macro damage mechanisms is that the loading environment and the material response to that environment are extremely complex. These types of damage are usually only observed in actual motor firings. Therefore, it is difficult and expensive to evaluate the reliability of new materials. Standard material tests which measure mechanical and thermal properties of test specimens can only provide a partial picture of how the material will respond in the service environment. The development of the ANALOG test procedure which can combine high heating rates and mechanical loads on a specimen will improve the understanding of the interactive effects of the various loads on the system. But a mechanistic model of material response which can account for the heterogeneity of the material, the progression of various micromechanical damage mechanisms, and the interaction of mechanical and thermal stresses on the material is required to accurately correlate material tests with response to service environments. A model based on fundamental damage mechanisms which is calibrated and verified under a variety of loading conditions will provide a general tool for predicting the response of rocket nozzles. The development of a micromechanical simulation technique was initiated and demonstrated to be effective for studying across-ply tensile failure of carbon/phenolic composites.

Slattery, Kerry T.↗
DOE OSTI2016

Atomistic Conversion Reaction Mechanism of WO 3 in Secondary Ion Batteries of Li, Na, and Ca

Abstract Intercalation and conversion are two fundamental chemical processes for battery materials in response to ion insertion. The interplay between these two chemical processes has never been directly seen and understood at atomic scale. Here, using in situ HRTEM, we captured the atomistic conversion reaction processes during Li, Na, Ca insertion into a WO 3 single crystal model electrode. An intercalation step prior to conversion is explicitly revealed at atomic scale for the first time for Li, Na, Ca. Nanoscale diffraction and ab initio molecular dynamic simulations revealed that after intercalation, the inserted ion–oxygen bond formation destabilizes the transition‐metal framework which gradually shrinks, distorts and finally collapses to an amorphous W and M x O (M=Li, Na, Ca) composite structure. This study provides a full atomistic picture of the transition from intercalation to conversion, which is of essential importance for both secondary ion batteries and electrochromic devices.

He, Yang↗
DOE OSTI2025

Dispatch Manager for NEML2 Constitutive Model Calculations Embedded in MOOSE

This report describes the extended capabilities of the NEML2 constitutive modeling library, including a flexible and efficient work dispatching system designed to leverage both CPU and GPU resources. This enhancement addresses one of the primary computational challenges in large-scale simulations: the ability to distribute and execute batches of material model evaluations across heterogeneous computing devices. The new dispatch system introduces a modular set of dispatcher and scheduler classes that coordinate the flow of data and execution between devices. The dispatcher is responsible for efficiently packaging work, managing device-specific memory operations, and synchronizing results. This modularity allows for extensibility, making it straightforward to integrate additional computing backends in the future. From an implementation standpoint, the dispatcher system interfaces seamlessly with NEML2's existing models. They handle device-aware tensor operations, optimize memory transfers, and support asynchronous execution when applicable. This design ensures that batches of material points can be evaluated concurrently, substantially improving throughput compared to previous single-device or serial implementations. These improvements not only enhance the raw performance of NEML2 but also improve its usability in multiscale and high-fidelity simulations, where the simultaneous evaluation of large material point batches is critical. Benchmarks included in the report demonstrate the system’s scalability, highlighting its effectiveness when leveraging modern GPU architectures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Nanoprecipitates to Enhance Radiation Tolerance in High-Entropy Alloys

The growth of advanced energy technologies for power generation is enabled by the design, development, and integration of structural materials that can withstand extreme environments, such as high temperatures, radiation damage, and corrosion. High-entropy alloys (HEAs) are a class of structural materials in which suitable chemical elements in four or more numbers are mixed to typically produce single-phase concentrated solid solution alloys (CSAs). Many of these alloys exhibit good radiation tolerance like limited void swelling and hardening up to relatively medium radiation doses (tens of displacements per atom (dpa)); however, at higher radiation damage levels (>50 dpa), some HEAs suffer from considerable void swelling limiting their near-term acceptance for advanced nuclear reactor concepts. In this study, we developed a HEA containing a high density of Cu-rich nanoprecipitates distributed in the HEA matrix. The Cu-added HEA, NiCoFeCrCu 0.12 , shows excellent void swelling resistance and negligible radiation-induced hardening upon irradiation up to high radiation doses (i.e., higher than 100 dpa). The void swelling resistance of the alloy is measured to be significantly better than NiCoFeCr CSA and austenitic stainless steels. Density functional theory simulations predict lower vacancy and interstitial formation energies at the coherent interfaces between Cu-rich nanoprecipitates and the HEA matrix. The alloy maintained a high sink strength achieved via nanoprecipitates and the coherent interface with the matrix at a high radiation dose (~50 dpa). From our experiments and simulations, the effective recombination of radiation-produced vacancies and interstitials at the coherent interfaces of the nanoprecipitates is suggested to be the critical mechanism responsible for the radiation tolerance of the alloy. Furthermore, the materials design strategy based on incorporating a high density of interfaces can be applied to high-entropy alloy systems to improve their radiation tolerance.

36 MATERIALS SCIENCE↗
DOE OSTI2024

High-purity germanium semiconductor modeling in the detector response function toolkit

In this study, we have extended the detector response function toolkit (DRiFT) to provide modeling capabilities of semiconductor sensors. DRiFT provides realistic nuclear instrumentation response by post-processing Monte-Carlo N-particle (MCNP®) radiation transport outputs. MCNP® is capable of modeling radiation transport in complex environments, but has limited detector physics and readout electronics modeling capabilities. Semiconductor detector response can be calculated with a high-fidelity for a flexible range of environments by utilizing MCNP® to simulate radiation interactions inside of detector volumes, and then using DRiFT to model charge transport and signal formation in the semiconductor, as well as the readout electronics. DRiFT models charge transport in the semiconductor, the preamplifier, shaping amplifier, pulse pile-up, and electronic noise to generate detector response. The semiconductor application in DRiFT can model a range of semiconductor materials, shapes, and sizes; and is demonstrated here for a large volume coaxial high-purity germanium (HPGe) detector. Here, we compare detector response functions of a coaxial HPGe detector with measurement of 60 Co, 133 Ba, and 137 Cs at varying count rates, and we conduct a parameter study to demonstrate the effect of changing parameters in the DRiFT simulation. The HPGe detector response function shows excellent agreement with measurements of difference sources with varying dead times and count rates.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

RHEGAL: Resistive heating gas enclosure loadframe for in situ neutron scattering

In situ neutron scattering is a powerful tool to reveal materials atomic structural response such as phase transformation, lattice straining, and texture under external stimuli. The advent of a high flux neutron source such as the Spallation Neutron Source (SNS) allows fast measurement in even non-equilibrium conditions, i.e., phase transformation in steels. However, the commercial fast heating apparatus such as commercial physical simulation equipment is not designed for in situ neutron scattering, which limits its application to in situ materials research by using neutrons. Here we present a resistive heating gas enclosure loadframe (RHEGAL) for non-equilibrium phase transformation studies by using in situ neutron scattering, which takes the advantage of high flux neutron sources like SNS. RHEGAL enables fast resistive heating of metal samples to 1200 °C at a rate up to 60 °C/s in an inert atmosphere. It provides both horizontal and vertical positions for scattering optimization. The mechanical loading capability also allows in situ high temperature tension above the oxidation temperature limit. The optimized translucent neutron scattering window by silicon allows both reflection and transmission measurements, making this equipment applicable for neutron diffraction, small angle scattering, and imaging. Furthermore, to demonstrate the fast heating capability, the phase transformations of an example of advanced high strength steel heated at 3 °C/s and 30 °C/s were measured with the VULCAN engineering diffractometer, and the different phase transformation kinetics by neutron diffraction were presented.

47 OTHER INSTRUMENTATION↗
DOE OSTI2017

Process Integration and Optimization of ICME Carbon Fiber Composites for Vehicle Lightweighting: A Preliminary Development

Process integration and optimization is the key enabler of the Integrated Computational Materials Engineering (ICME) of carbon fiber composites. In this paper, automated workflows are developed for two types of composites: Sheet Molding Compounds (SMC) short fiber composites, and multi-layer unidirectional (UD) composites. For SMC, the proposed workflow integrates material processing simulation, microstructure representation volume element (RVE) models, material property prediction and structure preformation simulation to enable multiscale, multidisciplinary analysis and design. Processing parameters, microstructure parameters and vehicle subframe geometry parameters are defined as the design variables; the stiffness and weight of the structure are defined as the responses. For multi-layer UD structure, this work focuses on the discussion of different design representation methods and their impacts on the optimization performance. Challenges in ICME process integration and optimization are also summarized and highlighted. Two case studies are conducted to demonstrate the integrated process and its application in optimization.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Enabling Nuclear and Solar Thermal Propulsion with the Computational Laboratory: Faster and Cheaper Materials Characterization

Nuclear and solar thermal propulsion offer enhanced efficiencies for in-space travel which may open up the outer Solar System to both crewed and automated craft. However, these advanced propulsion systems require cutting-edge materials which are often not easily tested in the laboratory. Ab initio thermodynamics and other computational methods offer a way to test materials at high pressures and temperatures more quickly, and with reduced cost. Here, several case studies are presented which highlight the utility of modern simulation techniques. First, an extensive thermodynamic study was performed at the ab initio level to elucidate the role of hot hydrogen flow on the erosion of four refractory carbides to examine their appropriateness for use as channel coatings. A detailed analysis of reaction products informs estimates of erosion which are validated in comparison with heritage NERVA data. Specific material suggestions are made for both nuclear and solar designs. Secondly, we examined the mechanical behavior and chemical stability of various components of an idealized nuclear cermet core design using techniques spanning from the atomic to microscale. Next, the mechanical and chemical behavior of tungsten (a potential matrix material) and the stability of uranium mononitride (a potential ceramic fuel element), are characterized in the presence of hydrogen and at elevated temperatures with ab initio thermodynamics. Finally, the mechanical response and failure of tungsten under tensile strain is simulated at the micron scale with a combined finite element and dislocation dynamics approach. The resultant stress-strain data agrees well with experiments and leads to a workflow which can incorporate ab initio thermodynamic data into micron-scale simulations and thus provide real-world relevant estimates of erosion and stress-strain behavior. The power of this framework and our plans to use it in the future will be discussed.

William C Tucker↗
DOE OSTI2024

Giant electric field-induced second harmonic generation in polar skyrmions

Electric field-induced second harmonic generation allows electrically controlling nonlinear light-matter interactions crucial for emerging integrated photonics applications. Despite its wide presence in materials, the figures-of-merit of electric field-induced second harmonic generation are yet to be elevated to enable novel device functionalities. Here, we show that the polar skyrmions, a topological phase spontaneously formed in PbTiO3/SrTiO3 ferroelectric superlattices, exhibit a high comprehensive electric field-induced second harmonic generation performance. The second-order nonlinear susceptibility and modulation depth, measured under non-resonant 800 nm excitation, reach ~54.2 pm V −1 and ~664% V −1 , respectively, and high response bandwidth (higher than 10 MHz), wide operating temperature range (up to ~400 K) and good fatigue resistance (>10 10 cycles) are also demonstrated. Through combined in-situ experiments and phase-field simulations, we establish the microscopic links between the exotic polarization configuration and field-induced transition paths of the skyrmions and their electric field-induced second harmonic generation response. Our study not only presents a highly competitive thin-film material ready for constructing on-chip devices, but opens up new avenues of utilizing topological polar structures in the fields of photonics and optoelectronics.

Nonlinear optics↗
DOE OSTI2023

A Novel Spatio-Temporal Regime Tracking Method for Impact Simulations

In this proposal, we present a novel method of tracking the rheological regimes activated during impact cratering events that will allow researchers to gain new insights into cratering mechanics. Rheology describes the stress-strain response of rocks to different conditions. Planetary impact cratering events often occur on too large of a scale to be feasibly captured in controlled experiments. Instead, these dynamic events are primarily studied using multi-physics codes equipped with complex material models that enable calculations of the impact event at scale. However, determining which physical processes are required for the problem of interest is challenging. Because the dominant rheological regimes change with space and time during crater formation, it is difficult to link numerical simulations with observable features of craters at the end of the event. The basis of this work is the implementation of numerical flags that track the activation of each rheological regime throughout impact simulations. We demonstrate this with the ‘Rock Model’ implemented in the CTH shock-physics code. We will use this rheology tracking method to ’zoom in’ on a specific region within an event and track the conditions the rock experiences over time. This work will develop community benchmarks to validate and distribute our implemented rheological models. We will focus on improving the melt models used by the planetary impact modeling community by developing an EOS-aware rheological transition from solid to melt. Through analysis of cell and tracer-particle based tracking data, we will study the effects of different rheologies on modeled outcomes, particularly on the volume and distributions of impacts melts. We will also use this method to link observable features with the rheological mechanisms responsible for them. The deliverables (peer-reviewed papers) from the proposed work are (i) tests of the implemented rheologic processes and demonstrations of the tracking flags; (ii) the first calculations of the spatio-temporal evolution of the dominant rheologies during impact cratering events; and (iii) application to delivery of impactor iron during basin-scale impacts.

58 GEOSCIENCES↗
NASA NTRS2012

Characterization and Analysis of Triaxially Braided Polymer Composites under Static and Impact Loads

In order to design impact resistant aerospace components made of triaxially-braided polymer matrix composite materials, a need exists to have reliable impact simulation methods and a detailed understanding of the material behavior. Traditional test methods and specimen designs have yielded unrealistic material property data due to material mechanisms such as edge damage. To overcome these deficiencies, various alternative testing geometries such as notched flat coupons have been examined to alleviate difficulties observed with standard test methods. The results from the coupon level tests have been used to characterize and validate a macro level finite element-based model which can be used to simulate the mechanical and impact response of the braided composites. In the analytical model, the triaxial braid unit cell is approximated by using four parallel laminated composites, each with a different fiber layup, which roughly simulates the braid architecture. In the analysis, each of these laminated composites is modeled as a shell element. Currently, each shell element is considered to be a smeared homogeneous material. Simplified micromechanics techniques and lamination theory are used to determine the equivalent stiffness properties of each shell element, and results from the coupon level tests on the braided composite are used to back out the strength properties of each shell element. Recent improvements to the model include the incorporation of strain rate effects into the model. Simulations of ballistic impact tests have been carried out to investigate and verify the analysis approach.

Goldberg, Robert K.↗
DOE OSTI2016

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

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

36 MATERIALS SCIENCE↗