Search NASA⌕ Search

SEARCH · Search NASA

Results for “Code assessment”

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,405 records · Page 78

Using Likwid and Byfl to Benchmark Hardware Performance

This paper outlines a benchmarking study conducted during my internship at LANL, focusing on CPU (Computer Processing Unit) and program performance assessment. The primary goal was to gather memory access data using three methods across five polybench kernels The data gathered would then be used to compare and contrast to one another and calculate operational intensity for performance comparisons. Benchmarking tools like Byfl and Likwid were employed, with Byfl offering hardware-independent data through LLVM compiler communication and Likwid directly interacting with computer hardware. The study considered various benchmarking factors, including optimization levels, Big O notation ((n)), CPU diversity and specific kernel equations. Big O notation was utilized to simplify code complexity, with detailed breakdwons of operations and memory components for each polybench application. Specific O(n) equations enabled nuanced kernel compariosns, facilitating the identification of performance variations. CPU efficiency assessments were conducted using Likwid tests on two CPUs. The central focus on code optimization aimed at achieving higher speeds and reduced memory usage through streamlined code. Future work propsoes creating a roofline model, synthesizing benchmarking data into a comprehensive data graph to assist in optimizing code and improving hardware performance. The potential impact on the laboratory or national mission was underscored, emphasizing the importance of optimizing applications and hardware to conserve resources and accelerate program execution. The specific relevance to LANL’s operations in math-intensive fields such as Nuclear Fission, Space Exploration, and Nanotechnology highlights the necessity of efficient benchmarking for resource conservation and proram speed. Overall, this study contributes to the understanding of CPU and program performance, providing insights for future optimization efforts in a laboratory setting

97 MATHEMATICS AND COMPUTING↗

Initial Probabilistic Evaluation of Reactor Pressure Vessel Fracture with Grizzly and Raven

The Grizzly code is being developed with the goal of creating a general tool that can be applied to study a variety of degradation mechanisms in nuclear power plant components. The first application of Grizzly has been to study fracture in embrittled reactor pressure vessels (RPVs). Grizzly can be used to model the thermal/mechanical response of an RPV under transient conditions that would be observed in a pressurized thermal shock (PTS) scenario. The global response of the vessel provides boundary conditions for local models of the material in the vicinity of a flaw. Fracture domain integrals are computed to obtain stress intensity factors, which can in turn be used to assess whether a fracture would initiate at a pre-existing flaw. These capabilities have been demonstrated previously. A typical RPV is likely to contain a large population of pre-existing flaws introduced during the manufacturing process. This flaw population is characterized stastistically through probability density functions of the flaw distributions. The use of probabilistic techniques is necessary to assess the likelihood of crack initiation during a transient event. This report documents initial work to perform probabilistic analysis of RPV fracture during a PTS event using a combination of the RAVEN risk analysis code and Grizzly. This work is limited in scope, considering only a single flaw with deterministic geometry, but with uncertainty introduced in the parameters that influence fracture toughness. These results are benchmarked against equivalent models run in the FAVOR code. When fully developed, the RAVEN/Grizzly methodology for modeling probabilistic fracture in RPVs will provide a general capability that can be used to consider a wider variety of vessel and flaw conditions that are difficult to consider with current tools. In addition, this will provide access to advanced probabilistic techniques provided by RAVEN, including adaptive sampling and parallelism, which can dramatically decrease run times.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Pilot's Automated Weather Support System (PAWSS) concepts

The FAA and NASA are investigating concepts for providing automated support to pilots in obtaining and using weather information. The pilot's Automated Weather Support System (PAWSS) activity is using computer technology to assist the pilot in acquiring and displaying weather information in the cockpit, and to assist in assessing and applying weather data in flight management and en route support. The PAWSS concepts include data link communications with a short update cycle of new weather data, color-coded map depictions with alpha-numeric data available for backup and detailed data presentation, automated monitoring of weather trends with alerting for hazardous conditions, and expert system technology applications in weather data assessments and flight management applications.

Crabill, N. L↗

Verification, Validation, and Uncertainty Quantification in Thermal Hydraulics, Freeman Scholar Lecture (2019)

Engineering problems are generally solved by analytical models or computer codes. These models, in addition to conservation equations, also include many empirical relationships and approximate numerical methods. Each of these components contributes to the uncertainty in the prediction. A systematic approach to judge the applicability of the code to the intended application is needed. It starts from verification of implementation of formulation in the code, identification of important phenomena, finding relevant tests with quantified uncertainty for these phenomena, and validation of the code by comparing predictions with the relevant test data. The relevant tests must address phenomena as expected in the intended application. In case of small size or limited condition tests, the scaling analyses are needed to assess the relevancy of the tests. Finally, a statement of uncertainty in the prediction is needed. Systematic approaches are described to aggregate uncertainties from different components of the code for intended application. Here, verification, validation, and uncertainty quantifications (VVUQs) are briefly described.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Experiments suitable for wind tunnel wall interference assessment/correction

Three experiments suitable for wall interference assessment and evaluation of proposed correction methods are presented. The experiments are: (1) a series of airfoil tests using a newly designed transonic flow facility that employs side-wall boundary layer suction and upper- and lower-wall shaping; (2) tests on a swept airfoil section spanning a solid-wall wind tunnel with fixed contouring on all four walls; and (3) tests on a swept wing of aspect ratio 3 mounted in a solid-wall wind tunnel with fixed flat walls. Each of the experiments provides data on the airfoil sections as well as on the wind tunnel walls. All the experiments were performed in solid wall wind tunnels corrected for boundary layer displacement effects. Although the experiments were performed primarily to evaluate computer code performance, it is believed that they also provide information that can be used to evaluate methods for assessing and correcting wall interference effects.

Marvin, J. G.↗

LSPRAY-IV: A Lagrangian Spray Module

LSPRAY-IV is a Lagrangian spray solver developed for application with parallel computing and unstructured grids. It is designed to be massively parallel and could easily be coupled with any existing gas-phase flow and/or Monte Carlo Probability Density Function (PDF) solvers. The solver accommodates the use of an unstructured mesh with mixed elements of either triangular, quadrilateral, and/or tetrahedral type for the gas flow grid representation. It is mainly designed to predict the flow, thermal and transport properties of a rapidly vaporizing spray. Some important research areas covered as a part of the code development are: (1) the extension of combined CFD/scalar-Monte- Carlo-PDF method to spray modeling, (2) the multi-component liquid spray modeling, and (3) the assessment of various atomization models used in spray calculations. The current version contains the extension to the modeling of superheated sprays. The manual provides the user with an understanding of various models involved in the spray formulation, its code structure and solution algorithm, and various other issues related to parallelization and its coupling with other solvers.

Raju, M. S.↗

Introducing PRIMRE's MRE Software Knowledge Hub (February 2021)

This paper focuses on the role of the Marine Renewable Energy (MRE) Software Knowledge Hub on the Portal and Repository for Information on Marine Renewable Energy (PRIMRE). The MRE Software Knowledge Hub provides online services for MRE software users and developers, and seeks to develop assessments and recommendations for improving MRE software in the future. Online software discovery platforms, known as the Code Hub and the Code Catalog, are provided. The Code Hub is a collection of open-source MRE software that includes a landing page with search functionality, linked to files hosted on the MRE Code Hub GitHub organization. The Code Catalog is a searchable online platform for discovery of useful (open-source or commercial) software packages, tools, codes, and other software products. To gather information about the existing MRE software landscape, a software survey is being performed, the preliminary results of which are presented herein. Initially, the data collected in the MRE software survey will be used to populate the MRE Software knowledge hub on PRIMRE, and future work will use data from the survey to perform a gap analysis and develop a vision for future software development. Additionally, as one of PRIMRE's roles is to support development of MRE software within project partners, a silo of knowledge relating to best practices has been gathered. An early draft of new guidance developed from this knowledge is presented.

gap analysis↗

Quantifying mean, variability, and uncertainty in indoor radon exposure in Pennsylvania using random forest and quantile regression forest models

Radon is a naturally occurring radioactive gas that poses a serious health risk as the primary cause of lung cancer in non-smokers. Despite the well-known adverse association with health outcomes, current radon exposure assessments are limited to county-level or average-level estimates, which fail to capture regional variability. This study uses Machine Learning models, including Random Forest (RF) and Quantile Regression Forest (QRF), to estimate the indoor radon concentrations at the ZCTA (Zip code tabulation area)-level and characterize uncertainties in model estimates. Incorporating geological, meteorological, and building-specific data, the models aim to improve radon risk assessment by capturing mean exposure, variability, and extreme concentration levels. Processed radon test data (n = 718,111) were analyzed using average, variability, and quantile prediction methods. Models that estimate the average radon exposure at the ZCTA-level can yield promising model-fit results, but they do not capture the underlying variability of indoor radon exposure within a ZCTA. We utilize volatility analyses to identify characteristics indicative of high variability of indoor radon exposure. We also show that a QRF model can be used to estimate upper quantiles of residential radon exposure, thereby uncovering localized areas of elevated exposure that were not apparent in mean estimates. The results highlighted the need for a deep characterization of exposure risk and show that regions with moderate average exposure levels could still harbor extreme outliers with implications for evaluating health risks. Utilizing multiple radon exposure models allows for a deeper characterization of radon risk within a geographic area and can better identify high-risk areas. The results from this study provide a foundation for developing mitigation strategies and examining associations between radon exposure and health outcomes at fine scales. Future research should extend the geographic scope and incorporate additional environmental risk factors to establish a comprehensive framework for risk assessment.

Lee, Heechan [ORNL]↗

Proposed Semi-Analytic Benchmark for Coupled Neutronics/Thermal-Hydraulics Feedback in Multiphysics Simulations

The past decade has seen considerable advancements in multiscale and multiphysics simulations that account for nonlinear effects of feedback between different physical processes within a single calculation. One area of particular interest has been multiphysics simulations that include both neutronics and thermal-hydraulics effects to accurately model the thermal feedback effects present in nuclear reactors. During this time, many groups have successfully modeled thermal feedback effects by coupling a variety of different types of neutron transport and thermal-hydraulic solvers together. However, formal verification and validation of these coupled code systems has proven challenging due to a lack of benchmark problems with well-characterized reference solutions. One approach is to confirm gross energy balance using bulk coolant properties (outlet temperature/enthalpy) for code verification purposes. Alternatively, code-to-code comparisons of results for a predefined reference/challenge problem have been used to establish reasonability of local results (neutron flux, heating, enthalpy rise) for various implementations, for example the VERA Core Physics Progression Problems. While community solutions are valuable for assessing the relative performance and fidelity of different implementations, they are not guaranteed to provide a measure of the absolute accuracy of a given implementation. In fact, these code-to-code comparisons can give a false sense of confidence by masking the effects of errors or approximations that are common to all of the solution methods tested. This is a particular concern when many methods derive from a common approach, such as Picard iteration. To this end, an analytical benchmark would be very valuable for verifying individual coupled physics implementations, as well as for providing a formal assessment of algorithm accuracy and sensitivity to model approximations such as mesh size and number of convergence iterations. However, creation of an analytical benchmark is difficult due to the complex nature of the eigenvalue form of the neutron transport equation (which does not permit a simple analytical solution for non-trivial cases) as well as the variety of different empirical closure relationships used in thermal-hydraulic modeling. In this paper we propose a benchmark for coupled neutronics/thermal-hydraulics simulations. The proposed benchmark overcomes the challenges associated with a multiphysics benchmark problem by using a fixed-source (rather than eigenvalue) formulation of the neutron transport equation and limiting fluid properties to the single-phase regime where the density and temperature of the coolant are related to enthalpy by a known relationship, while still preserving significant nuclear/thermal feedback effects within the simplified problem. The new benchmark is easy to simulate using existing coupled physics implementations, and allows for a direct comparison with a reference solution, as provided in this summary.

42 ENGINEERING↗

Performance Assessment of a Generic Nuclear Waste Repository in Shale - 18392

A performance assessment (PA) is an important component of a comprehensive safety analysis for a nuclear waste repository. In a PA, probabilistic simulations of the total repository system are performed, and results are evaluated against performance metrics. Uncertainty and sensitivity analyses help prioritize additional research and model development. The United States Department of Energy (DOE) has been developing a state of the art PA simulation software tool-kit, Geologic Disposal Safety Assessment (GDSA) Framework, that couples increasingly higher fidelity models of subsystem processes into total system PA simulations. Over the past several years, PAs of generic repositories in three geologic media (salt, shale, and crystalline rock) have demonstrated ongoing developments in capability. The current PA of a nuclear waste repository in a generic shale formation showcases GDSA Framework, including capabilities in domain discretization (Cubit), multi-physics simulations (PFLOTRAN), uncertainty and sensitivity analysis (Dakota), and visualization (Paraview). The generic shale reference case considers the disposal of 22,000 metric tons heavy metal of commercial spent nuclear fuel (SNF) in a generic shale formation. PA simulations account for the thermal load and radionuclide inventory of the waste form, components of the engineered barrier system including waste package and buffer, and components of the natural barrier system including the host rock shale and underlying and overlying stratigraphic units. Two repository layouts are considered, one for emplacement of waste packages containing 12 pressurized water reactor (PWR) assemblies, and one for 4-PWR waste package emplacement. Model domains are half-symmetry and contain between 7 and 22 million grid cells. Grid refinement captures the detail of individual waste packages, emplacement drifts, access drifts, and shafts. Simulations are run in a high-performance computing (HPC) environment on 512 or 2048 processes. The governing equations describing coupled heat and fluid flow and reactive transport are solved with PFLOTRAN, an open-source, massively parallel multiphase flow and reactive transport code. Additional simulated processes include waste package degradation; waste form dissolution; radioactive decay and ingrowth in aqueous, solid, and adsorbed phases; and a simple biosphere based on a pumping well. Simulations are run to 106 y; dose and radionuclide concentrations are observed within aquifers at a point approximately 5 km downgradient of the repository. Dakota is used to sample likely ranges of input parameters including waste form and waste package degradation rates and properties of engineered and natural materials to quantify uncertainty in predicted concentrations and sensitivity to input parameters. Given the assumptions of the reference case, results of 12-PWR and 4-PWR simulations are very similar. Dose at the pumping well 5-km downgradient of the repository does not exceed 10-9 Sv/y. I-129 is the largest contributor to dose. The dose due to Cl-36, the next largest contributor, is orders of magnitude smaller. At the pumping well, the concentration of I-129 5-km downgradient does not exceed 10-14 mol/L. Rank correlation coefficients indicate that the maximum concentration of I-129 within the aquifer is sensitive to shale (repository host rock) porosity and aquifer permeability. These results affirm that HPC-capable codes can be used to simulate important multi-physics couplings directly in a total system performance assessment of a geologic repository, and can be used in prioritization of future research and model development. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Evalulation of the Specification for Safety-Related Steel Structures for Nuclear Facilities, ANSI/AISC N690-18, for Application to Nuclear Power Plants

This report describes the assessment of the recently published ANSI/AISC N690-18, Specification for Safety-Related Steel Structures for Nuclear Facilities. This work was performed for the U.S. Nuclear Regulatory Commission for potential endorsement of N690-18 for use in the design of nuclear power plants. Currently the NRC endorses the use of an earlier version: N690-1994 (R2004), and thus, a technical review of the updated N690 Specification, which contains significant changes, is warranted. The N690-18 Specification reflects advances and improvements in steel design gained from experience, analytical studies, and experimental data over the years. The N690-18 specification also contains a new Appendix N9 that addresses the design of safety-related steel-plate composite (SC) walls which are a relatively new type of structural member for which previously there has not been any U.S. code or standard that governs its design. The technical review of N690-18 was performed in order to assess the adequacy of the provisions for use by the NRC in updating its regulatory guidance for nuclear power plants. The research developed the technical basis for acceptance of the new specification and also identified areas where additional staff guidance is needed for the design of safety-related steel structures at nuclear power plants. The results of the research have determined that the N690-18 specification is a significant enhancement from the N690-1994 (R2004) edition. The new specification contains substantial upgrades in the design of structural steel and includes new provisions for design of SC wall sections, which can benefit the nuclear power industry as well as the NRC review process. All of the improvements, if endorsed by the NRC with certain regulatory positions, will allow clearer and more consistent regulatory guidance, and will avoid the need for justification by licensees and related reviews each time certain methods or approaches are used.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Evaluation of the Specification for Safety-Related Steel Structures for Nuclear Facilities, ANSI/AISC N690-18, for Application to Nuclear Power Plants

This report describes the assessment of the recently published ANSI/AISC N690-18, Specification for Safety-Related Steel Structures for Nuclear Facilities. This work was performed for the U.S. Nuclear Regulatory Commission for potential endorsement of N690-18 for use in the design of nuclear power plants. Currently the NRC endorses the use of an earlier version: N690-1994 (R2004), and thus, a technical review of the updated N690 Specification, which contains significant changes, is warranted. The N690-18 Specification reflects advances and improvements in steel design gained from experience, analytical studies, and experimental data over the years. The N690-18 specification also contains a new Appendix N9 that addresses the design of safety-related steel-plate composite (SC) walls which are a relatively new type of structural member for which previously there has not been any U.S. code or standard that governs its design. The technical review of N690-18 was performed in order to assess the adequacy of the provisions for use by the NRC in updating its regulatory guidance for nuclear power plants. The research developed the technical basis for acceptance of the new specification and also identified areas where additional staff guidance is needed for the design of safety-related steel structures at nuclear power plants. The results of the research have determined that the N690-18 specification is a significant enhancement from the N690-1994 (R2004) edition. The new specification contains substantial upgrades in the design of structural steel and includes new provisions for design of SC wall sections, which can benefit the nuclear power industry as well as the NRC review process. All of the improvements, if endorsed by the NRC with certain regulatory positions, will allow clearer and more consistent regulatory guidance, and will avoid the need for justification by licensees and related reviews each time certain methods or approaches are used.

42 ENGINEERING↗

Large-Eddy Simulation Sensitivities to Variations of Configuration and Forcing Parameters in Canonical Boundary-Layer Flows for Wind Energy Applications

Here, the sensitivities of idealized Large-Eddy Simulations (LES) to variations of model configuration and forcing parameters on quantities of interest to wind power applications are examined. Simulated wind speed, turbulent fluxes, spectra and cospectra are assessed in relation to variations of two physical factors, geostrophic wind speed and surface roughness length, and several model configuration choices, including mesh size and grid aspect ratio, turbulence model, and numerical discretization schemes, in three different code bases. Two case studies representing nearly steady neutral and convective atmospheric boundary layer (ABL) flow conditions over nearly flat and homogeneous terrain were used to force and assess idealized LES, using periodic lateral boundary conditions. Comparison with fast-response velocity measurements at five heights within the lowest 50 m indicates that most model configurations performed similarly overall, with differences between observed and predicted wind speed generally smaller than measurement variability. Simulations of convective conditions produced turbulence quantities and spectra that matched the observations well, while those of neutral simulations produced good predictions of stress, but smaller than observed magnitudes of turbulence kinetic energy, likely due to tower wakes influencing the measurements. While sensitivities to model configuration choices and variability in forcing can be considerable, idealized LES are shown to reliably reproduce quantities of interest to wind energy applications within the lower ABL during quasi-ideal, nearly steady neutral and convective conditions over nearly flat and homogeneous terrain.

17 WIND ENERGY↗

Large-eddy simulation sensitivities to variations of configuration and forcing parameters in canonical boundary-layer flows for wind energy applications

Abstract. The sensitivities of idealized large-eddy simulations (LESs) to variations of model configuration and forcing parameters on quantities of interest to wind power applications are examined. Simulated wind speed, turbulent fluxes, spectra and cospectra are assessed in relation to variations in two physical factors, geostrophic wind speed and surface roughness length, and several model configuration choices, including mesh size and grid aspect ratio, turbulence model, and numerical discretization schemes, in three different code bases. Two case studies representing nearly steady neutral and convective atmospheric boundary layer (ABL) flow conditions over nearly flat and homogeneous terrain were used to force and assess idealized LESs, using periodic lateral boundary conditions. Comparison with fast-response velocity measurements at 10 heights within the lowest 100 m indicates that most model configurations performed similarly overall, with differences between observed and predicted wind speed generally smaller than measurement variability. Simulations of convective conditions produced turbulence quantities and spectra that matched the observations well, while those of neutral simulations produced good predictions of stress, but smaller than observed magnitudes of turbulence kinetic energy, likely due to tower wakes influencing the measurements. While sensitivities to model configuration choices and variability in forcing can be considerable, idealized LESs are shown to reliably reproduce quantities of interest to wind energy applications within the lower ABL during quasi-ideal, nearly steady neutral and convective conditions over nearly flat and homogeneous terrain.

17 WIND ENERGY↗

TH-A-19A-10: Fast Four Dimensional Monte Carlo Dose Computations for Proton Therapy of Lung Cancer

Purpose: To develop and validate a fast and accurate four dimensional (4D) Monte Carlo (MC) dose computation system for proton therapy of lung cancer and other thoracic and abdominal malignancies in which the delivered dose distributions can be affected by respiratory motion of the patient. Methods: A 4D computer tomography (CT) scan for a lung cancer patient treated with protons in our clinic was used to create a time dependent patient model using our in-house, MCNPX-based Monte Carlo system (“MC{sup 2}”). The beam line configurations for two passively scattered proton beams used in the actual treatment were extracted from the clinical treatment plan and a set of input files was created automatically using MC{sup 2}. A full MC simulation of the beam line was computed using MCNPX and a set of phase space files for each beam was collected at the distal surface of the range compensator. The particles from these phase space files were transported through the 10 voxelized patient models corresponding to the 10 phases of the breathing cycle in the 4DCT, using MCNPX and an accelerated (fast) MC code called “FDC”, developed by us and which is based on the track repeating algorithm. The accuracy of the fast algorithm was assessed by comparing the two time dependent dose distributions. Results: The error of less than 1% in 100% of the voxels in all phases of the breathing cycle was achieved using this method with a speedup of more than 1000 times. Conclusion: The proposed method, which uses full MC to simulate the beam line and the accelerated MC code FDC for the time consuming particle transport inside the complex, time dependent, geometry of the patient shows excellent accuracy together with an extraordinary speed.

60 APPLIED LIFE SCIENCES↗

Bumper Implementation of LMEEM

The Hypervelocity Impact Technology (HVIT) group at the NASA Johnson Space Center (JSC) maintains the Bumper 3 computer program (referred to as “Bumper”) to perform Micrometeoroid and Orbital Debris (MMOD) risk assessments for spacecraft. To perform its calculations, this program requires detailed mathematical models of the space environments and can include those models representing orbital debris, micrometeoroids, and lunar ejecta. The lunar ejecta environment is a projectile environment specific to the lunar surface. Due to the Moon’s lack of an atmosphere, meteoroids frequently strike the lunar surface. When this happens, the impact ejects material radially away from the impact point, which can pose a potential impact hazard to surface-based hardware. This lunar ejecta environment was previously described in 1969, in NASA Special Publication SP-8013, which was used as the basis for a model that will be referred to in this document as the “SP-8013 environment model,” or simply, “SP-8013.” Since around 2020, this model has been used in MMOD risk assessments for the Human Landing System (HLS) program and other lunar surface programs. A new lunar meteoroid ejecta environment was developed via a computer code named Lunar Meteoroid Ejecta Engineering Model (LMEEM), produced by the Natural Environments Branch (NEB) at Marshall Space Flight Center (MSFC). A global lunar surface environment definition using LMEEM is to be provided in the Cross-Program Design Specification for Natural Environments (DSNE) Revision J in 2024. Contact the NASA MSFC Natural Environments Branch for an in-depth discussion of the LMEEM environment model and its derivation. This environment model will supersede SP-8013 for MMOD risk assessments. To perform such assessments, Bumper must read data from text files output by LMEEM and use that data in risk calculations. The purpose of this document is to: •Demonstrate verification that the Bumper code has properly implemented the LMEEM environment model (Section 3) •Discuss the LMEEM environment as understood through an initial environment file (Section 4) •Compare the calculated risk due to the SP-8013 and LMEEM models (Section 5)

Bumper↗

Application of a computational glass model to the shock response of soda-lime glass

This article details the implementation and application of the glass-specific computational constitutive model by Holmquist and Johnson [1] to simulate the dynamic response of soda-lime glass under high rate and high pressure shock conditions. The predictive capabilities of this model are assessed through comparison of experimental data with numerical results from computations using the CTH shock physics code. The formulation of this glass model is reviewed in the context of its implementation within CTH. Using a variety of experimental data compiled from the open literature, a complete parameterization of the model describing the observed behavior of soda-lime glass is developed. Simulation results using the calibrated soda-lime glass model are compared to flyer plate and Taylor rod impact experimental data covering a range of impact and failure conditions spanning an order of magnitude in velocity and pressure. In conclusion, the complex behavior observed in the experimental testing is captured well in the computations, demonstrating the capability of the glass model within CTH.

36 MATERIALS SCIENCE↗

Data driven modeling of plastic deformation

In this paper the application of machine learning techniques for the development of constitutive material models is being investigated. A flow stress model, for strain rates ranging from 10 –4 to 10 12 (quasi-static to highly dynamic), and temperatures ranging from room temperature to over 1000 K, is obtained by beginning directly with experimental stress-strain data for Copper. An incrementally objective and fully implicit time integration scheme is employed to integrate the hypo-elastic constitutive model, which is then implemented into a finite element code for evaluation. Accuracy and performance of the flow stress models derived from symbolic regression are assessed by comparison to Taylor anvil impact data. The results obtained with the free-form constitutive material model are compared to well-established strength models such as the Preston-Tonks-Wallace (PTW) model and the Mechanical Threshold Stress (MTS) model. Here, preliminary results show candidate free-form models comparing well with data in regions of stress-strain space with sufficient experimental data, pointing to a potential means for both rapid prototyping in future model development, as well as the use of machine learning in capturing more data as a guide for more advanced model development.

36 MATERIALS SCIENCE↗