Search NASA⌕ Search

SEARCH · Search NASA

Results for “Code development”

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 595 records · Page 33

Machine learning at the Spallation Neutron Source accelerator and target

We describe the ongoing efforts to apply Machine Learning techniques to improve the performance of our accelerator and target. Specially, we are looking to minimize halo beam losses in the absence of a proper physics model, automatically detect and log anomalies in the target support systems such as cooling, and detect and prevent errant beam pulses in the linac. We also describe the infrastructure we use to acquire and stream data to the GPU cluster for training, our code development cycle, and edge computing for model inference. To minimize halo beam losses, we use a Reinforcement Learning technique tested on a virtual accelerator. The target anomaly detection is trained on archived data using incomplete physics models and is made part of the existing target reporting system. The errant beam prevention analyzes beam current and beam phase waveforms as well as accelerator configuration data to predict errant pulses. We also develop continual learning to adapt to changes in the accelerator.

Accelerator Physics↗

Comparison of Generation of Higher-Order Neutron Scattering Cross Sections

The generation of high-order neutron scattering cross sections consistent with high-fidelity simulations remains an area of active research. Popular options include generating cross sections from continuous energy Monte Carlo calculations or generating cross sections from a deterministic neutron transport calculation with high-fidelity tabulated cross sections. Both options present challenges. Monte Carlo simulations can naturally process continuous energy cross section data and allow for the general description of anisotropic neutron scattering given a scattering law. However, Monte Carlo simulations are inherently related to particle weighting and it has been suggested that this may be unacceptable for generating high-order neutron scattering cross sections. Deterministic neutron transport calculations can easily calculate high-order moments of the flux to appropriately calculate higher-order neutron scattering cross sections but are generally limited by discretization of space, energy, and angle. In this work, the trade-offs between generation of high-order neutron scattering cross sections via Monte Carlo and deterministic neutron transport methods are investigated. The methods implemented in the Monte Carlo computer program Serpent and the deterministic fast reactor neutron cross section generator MC2-3 are compared. Cross sections resulting from these methods are used in Rattlesnake, a deterministic neutron transport code developed by Idaho National Laboratory (INL), and results are compared to a reference continuous energy Monte Carlo calculation. Whereas previous work investigating the effects of anisotropic neutron scattering has focused on light water reactor (LWR) simulations, this work focuses on high-order neutron scattering cross sections as they relate to fast reactor simulations. To investigate the consequences of the Serpent and MC2-3 methodologies, a test problem is developed. The test problem is a one-dimensional geometry with fast reactor materials designed to demonstrate deep penetration and exacerbate the effects of high-order neutron scattering. Based on the results of the deep penetration test problem, it is concluded that P3 neutron scattering cross sections are sufficient to describe anisotropic scattering in fast reactor materials. Any neutron scattering of order higher than P3 offers negligible change in the eigenvalue of the test problem. Additionally, it is determined that the methodology as implemented in Serpent is applicable for generating high-order neutron scattering cross sections through at least P3 in fast reactor materials.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Resiliency of Degraded Built Infrastructure

Infrastructure resiliency depends on the ability of infrastructure systems to withstand, adapt, and recover from chronic and extreme stresses. In this white paper, we address the resiliency of infrastructure assets and discuss improving infrastructure stability through development of our understanding of cement and concrete degradation. The resiliency of infrastructure during extreme events relies on the condition, adaptability, and recoverability of built infrastructure (roads, bridges, dams), which serves as the backbone of existing infrastructure systems. Much of the built infrastructure in the US has consistently been rated D+ by the American Society of Civil Engineers (ASCE). Aged infrastructure introduces risk to the system, since unreliable infrastructure increases the likelihood of failures under chronic and extreme stress and are particularly concerning when extreme events occur. To understand and account for this added risk from poor infrastructure quality, more research is needed on (i) how the changing environment alters the aging of new and existing built infrastructure and (ii) how degradation causes unique failure mechanisms. The aging of built infrastructure is based on degradation of the structural materials, such as concrete and steel supports, which causes failure. Current work in cement/concrete degradation is based on (i) the development of high strength and degradation resistance concrete mixtures, (ii) methods of assessing the age and reliability of existing structures, and (3) modeling of structural stability and the microstructural evolution of concrete/cement from degradation mechanisms (sulfide attack, carbonation, decalcification). Sandia National Laboratories (SNL) has made several investments in studying the durability and degradation of cement based materials, including using SNL-developed codes and methodologies (peridynamics, PFLOTRAN) to focus on chemo-mechanical fracture of cement for energy applications. Additionally, a recent collaboration with the University of Colorado Boulder has included fracture of concrete gravity dams, scaling the existing work to applications in full sized infrastructure problems. Ultimately, SNL has the experience in degradation of cementitious materials to extend the current research portfolio and answer concerns about the resilience of aging built infrastructure.

36 MATERIALS SCIENCE↗

Stockpile Assessment (Quad Chart)

Our goal was to characterize certain aspects of shaped charges. In order to determine the pressure field created by the jet and the jet velocity, I worked with and modified simulations using CTH, a code developed at Sandia. I was able to manipulate and modify certain variables in order to observe their effects and measure the pressure and velocities at different points in the simulation.

42 ENGINEERING↗

Institutional Computing: Annual Progress Report, Precision early universe simulations to constrain nuclear reactions and beyond standard model physics

The current one-year allocation (w20 qburst) is the first year in a two year allocation of the same name. The project has supported two publications in peer-review journals, another manuscript in preparation for publication, a white paper for the Astro2020 Decadal Survey, and a number of talks (see below). Code development continues on implementing a neutrino quantum kinetic equation solver into BURST, now called QBURST. While implementing coherent and forward-scattering processes into the dynamical neutrino evolution equations, multiple numerical issues have appeared. We have remedied a number of these issues through extensive debugging and testing, which include streamlining interpolation and implementing new ODE integration techniques. Further testing is required to solve the problem at a precision of better than one part in 10 6 . We expect that we will understand the numerical difficulties and implement solutions by end of year 2021.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Increasing the Usability of Sockeye: Implementing Customer-Driven Modeling Improvements, Features, and Documentation of Sockeye

A number of tasks for the heat pipe simulation tool Sockeye were identified based on customer needs and included both user documentation and code development: (a) Implement saturation-dependence in the capillary pressure model; (b) Provide an analytic utility to analyze the frozen startup limit; (c) Demonstrate custom friction closures; (d) Develop examples for different types of heat transfer sources and boundary conditions; (e) Demonstrate the capillary and sonic limits; and (f) Provide discretization details in the theory manual. These tasks were all completed, and results are given where applicable.

97 MATHEMATICS AND COMPUTING↗

Isotopic and Fuel Lattice Parameter Trends in Extended Enrichment and Higher Burnup LWR Fuel, Volume II: BWR

Commercial light water reactor (LWR) operators and fuel vendors in the United States are pursuing changes to nuclear fuel that include extended enrichment (EE) and accident-tolerant fuel (ATF) designs to further improve reactor safety and plant economics. Extended fuel enrichments above 5% 235 U pin enrichment and up to 10% 235 U are a subset of high assay low-enriched uranium (HALEU) that may be deployed in commercial US LWRs in the near term. ATF features, such as cladding coatings or alternative cladding materials, are designed to improve fuel system performance under accident conditions. One goal of EE is to improve fuel cycle economy by enabling fuel to be depleted to higher burnup than the typical current limits (62 gigawatt-days per metric ton of uranium [GWd/MTU]). Adoption of EE, ATF, and high burnup (HBU) fuels in the US commercial fleet requires a clear understanding of the effects on core physics parameters and used fuel isotopic content, as well as confidence in the accuracy of computer code predictions over an expanded range of materials, enrichment, and burnup. A thorough understanding of the applicability and adequacy of benchmark data (e.g., criticality, decay heat, isotopic content) for computer code validation is necessary to ensure that appropriate safety margins are maintained. As part of the US Nuclear Regulatory Commission (NRC) agreement number 31310019N0008, “SCALE Code Development, Assessment and Maintenance,” the effects of EE, ATF, and HBU are being assessed for selected representative LWR fuel designs. The project is divided into phases, and this report summarizes the findings of the Phase 1 work, which focuses on lattice physics parameter and used fuel isotopic changes for a conventional GE14 10 x 10 boiling water reactor (BWR) design with GNF-2 part length rod patterns to model a modern BWR assembly design.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Large Eddy Simulations of Turbulence below Antarctic Ice Shelves

Turbulence below ice shelves is of key importance for predicting ice-shelf melt rates and consequently the contribution of ice sheets to sea level rise. In this project we conducted Large-Eddy Simulations (LES) to improve our understanding of sub-ice-shelf ocean turbulence and the relationship between ice-shelf melt rates and ocean conditions. Over the course of the second and final year of this project (FY2020), we accomplished both major code developments for the PArallel Large-eddy simulation Model (PALM; Maronga et al., 2015) and conducted a suite of simulations which form the basis for a publication in preparation.

58 GEOSCIENCES↗

Statistical Distributions for Mesh Independent Solutions in ALEGRA

The representation of material heterogeneity (also referred to as "spatial variation") plays a key role in the material failure simulation method used in ALEGRA. ALEGRA is an arbitrary Lagrangian-Eulerian shock and multiphysics code developed at Sandia National Laboratories and contains several methods for incorporating spatial variation into simulations. A desirable property of a spatial variation method is that it should produce consistent stochastic behavior regardless of the mesh used (a property referred to as "mesh independence"). However, mesh dependence has been reported using the Weibull distribution with ALEGRA's spatial variation method. This report describes efforts towards providing additional insight into both the theory and numerical experiments investigating such mesh dependence. In particular, we have implemented a discrete minimum order statistic model with properties that are theoretically mesh independent.

36 MATERIALS SCIENCE↗

2D Magnetohydrodynamic Simulations of the Electrothermal Instability in Metallic Liners

The Virginia Tech (VT) Plasma Dynamics Laboratory Computational (PDCL) and Lawrence Livermore National Laboratory (LLNL) are performing two dimensional (2D) simulations of the electrothermal instability (ETI) using the LLNL multi-physics code Ares. Ares is a multi-physics arbitrary–Lagrangian-Eulerian (ALE) code developed by LLNL and is of particular use in studying magnetohydrodynamic (MHD) instabilities like the ETI due to its resistive MHD, magnetic diffusion, and radiative-hydrodynamics packages. Among its capabilities, it has the ability to model material strength, perform adaptive mesh refinement (AMR), and incorporate a wide variety of equations of state models and conductivity models. The 2D Ares simulation model created by VT-LLNL for studying the development and growth of the electrothermal instability has been configured with initial conditions based on the Mykonos Electrothermal Instability II (METI-II) experiments described by this grant and conducted by team members at the University of Nevada (UNR), the University of New Mexico (UNM), and Sandia National Laboratories. Previously, preliminary 2D Ares simulations of the ETI had been run to approximately 80ns. The rods in these simulations were initiated with sinusoidal perturbations at a similar order of magnitude to those measured on the aluminum rods used for the Mykonos experiment. This model has been improved by increasing the spatial resolution of the simulations and running the simulations further in time. In addition to the simulation run-times extending, the preliminary sinusoidal perturbation has been replaced with a perturbation derived from amplitude measurements by the experimental team, thereby correlating the simulation inputs better to the experimental runs. These new runs are capable of reaching 120ns of simulated time for the uncoated cases and to 200ns the 41 μm coated cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Profiling and Optimization [Slides]

Intended for Simulating Physics using Efficient and Effective code Development (SPEED) program. This is a LANL internal lecture series for computational science application developers.

97 MATHEMATICS AND COMPUTING↗

Jan-Dec 2021 IC Project: w20_earthquake (Yearly Report)

The simulations conducted under this IC project used an in-house code developed at Los Alamos National Laboratory called HOSS (Hybrid Optimization Software Suite) and its extended version called HOSS-FSIS (fluid-structure interaction solver). HOSS has been under development in LANL’s Geophysics Group since 2004 and was used to simulate the signal propagation in pieces of material with pre-existing defects and also to simulate the near source effects of above ground nuclear explosions. One of the key feature of HOSS is that it is superior to pure finite element method and discrete element method implementations because it allows for a natural transition from continuum to discontinuum when the material breaks or fractures.

58 GEOSCIENCES↗

Update of CEA DES Criticality-Safety Activities and Perspectives [Slides]

Neutronics staff at CEA DES are in charge of a whole set of modeling and simulation tools for neutronics: nuclear data evaluation and processing, transport codes development, calculation sequences edition, verification, validation and uncertainty analysis. Advanced user groups are in charge of specific nuclear analysis in reactor physics, fuel cycle, criticality-safety, radiation shielding and nuclear instrumentation for CEA and its partners, at all stages of nuclear facilities life cycle. CEA is willing to participate/collaborate on the following issues: detector design, experiments, analysis and evaluation, DH validation by combining TAGS measurements with Fission Yield variance-covariance matrices.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

LLNL Macroscopic Anisotropic Explosives Research At INL National Security Test Range: Shot Narrative

This series of experiments is exploring new frontiers in explosives science, technology, and engineering. The extensive list of experiment shots is necessary to sort out what works and what does not work in searching for an environment that supports anisotropic detonation properties in explosives. These experiments have been designed to serve as an economical and efficient first-cut determination as to whether there is any validity to an explosive analog to a series circuit of electrical diodes. That is, the detonation moves more easily in the “forward” direction than the “backward” direction. The basic data gleaned from these experiments will be used to design experiments that will be shot in the 10 kg tank at the High Explosives Application Facility (HEAF) at the Lawrence Livermore National Laboratory. In turn, those shots that will be heavy in diagostics and instrumentation, providing basic data that will be used by simulation code developers in creating high explosive models that account for air gap and oblique explosive surfaces undergoing sympathetic detonation.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Final Report DE-FE0031785 Mohsen Ahmadian, Ph.D. Demonstration of Proof of Concept of a Multiphysics Approach for Real-Time Remote Monitoring of Dynamic Changes in Pressure and Salinity in Hydraulically Fractured Networks

Hydraulic fracturing has evolved into a multistep process with varying flow rates, carrier fluids (e.g., gel or slickwater), proppant loadings, and proppant grain sizes. As a result, primary recovery from a hydraulically fractured tight-oil reservoir is often a tiny fraction of the original oil in place, ranging between 5 and 10%. As stated in the FOA1990, “part of this problem is due to the inability of current well completion processes to effectively stimulate the entire reservoir volume in contact with the wellbore. Innovative technologies are needed that can help improve the effectiveness of reservoir completion methods, maximize stimulated reservoir volumes, and optimize recovery over the entire producing life span of a well”. We first need to enhance the current fracture diagnostic techniques to improve a well-completion design. However, detecting and delineating a subsurface hydraulic fracture is extremely difficult because the induced fracture network is only fractionally propped, and these propped fractures are generally very thin. Microseismic and tiltmeter monitoring techniques can provide information on the fracture extent but provide little or no information on the movement and final distribution of proppant or production fluids. On the other hand, electromagnetic (EM) imaging has shown the capability to monitor proppant distribution throughout the fracture area, especially in the presence of Electrically Active Proppants (EAPs). A previous EM survey of hydraulic fracturing at the Devine Fracture Pilot Site (DFPS) and subsequent EM code developments demonstrated this survey as a robust technique to remotely interrogate the extent of the EAP-filled hydraulic fracture during its propagation. The objectives of the project were threefold: (1) to capitalize on the material properties of an EAP to demonstrate remote monitoring of relative changes in pressure, pressure, and flow that are commonly encountered during production from a hydraulically fractured reservoir; (2) to evaluate EM imaging tools, to achieve Objective 1 in near real-time; and (3) to develop a multi-physics joint inversion approach to precisely predict flow patterns and physiochemical changes within an EAP-filled fracture network. This research project was built upon our previous work at the Devine Test Site managed by the Bureau of Economic Geology (BEG) at The University of Texas at Austin (UT-Austin). It also leveraged a significant investment from the Advanced Energy Consortium (AEC) to address the DOE's interest in subsurface flow, containment, and characterization by multiple signals. This three-year and three-month project succeeded in demonstrating the feasibility of a real-time dynamic fluid flow mapping technique at Technology Readiness Level 5 (TRL5) by utilizing a commercially available surface-based Controlled-Source Electromagnetic (CSEM) method (Objectives 1, 2). We demonstrated that injections into an EAP-filled fracture could be successfully coupled with real-time electric field measurements on the surface, leading to remote monitoring of dynamic changes within the EAP-filled fracture. Furthermore, the observed electric field in our study is influenced by bottomhole pressure, flow rate, and salinity, which is demonstrated by comparing these parameters with the electrical field potentials. EM simulations solely based on assumptions of fracture conductivity changes during injection did not reproduce the whole measured electric field magnitudes. Preliminary estimates showed that including Streaming Potential (SP) in our geophysical model is likely needed to reduce the simulation misfit.

02 PETROLEUM↗

Efficient Multi-Scale, Thermo-Hydro-Mechanical Numerical Model for Simulation of Long-Term Stability of Rooms for Nuclear Waste Disposal

This project demonstrated the feasibility of developing a numerical modeling tool that provides an efficient and unique approach to analyze the geomechanical performance of a deep geological repository (DGR) for high-level nuclear waste. The numerical modeling tool is physics-based and models the thermo-hydro-mechanical (THM) processes during the required time (up to one million years) for a DGR. The unique capability of the developed code is explicit representation of fracturing and a discrete fracture network (DFN) in the model that includes the entire repository layout for simulated time of one million years. Therefore, it can potentially address one of the challenges of analysis of a DGR, which is efficient numerical simulation of a multi-scale, physics-based models.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Improvements in MCNP6.3 (Rev. 1) [Slides]

This presentation outlines improvements in the Monte Carlo N-Particle Code (MCNP) version 6.3. MCNP is an internationally renowned Monte Carlo particle transport code developed by Los Alamos National Laboratory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Updated Economic Model for Estimation of GDP Losses in the MACCS Offsite Consequence Analysis Code RDEIM Model Report for MACCS v4.2

This report updates the Regional Disruption Economic Impact Model (RDEIM) GDP-based model described in Bixler et al. (2020) used in the MACCS accident consequence analysis code. MACCS is the U.S. Nuclear Regulatory Commission (NRC) used to perform probabilistic health and economic consequence assessments for atmospheric releases of radionuclides. It is also used by international organizations, both reactor owners and regulators. It is intended and most commonly used for hypothetical accidents that could potentially occur in the future rather than to evaluate past accidents or to provide emergency response during an ongoing accident. It is designed to support probabilistic risk and consequence analyses and is used by the NRC, U.S. nuclear licensees, the Department of Energy, and international vendors, licensees, and regulators. The update of the RDEIM model in version 4.2 expresses the national recovery calculation explicitly, rather than implicitly as in the previous version. The calculation of the total national GDP losses remains unchanged. However, anticipated gains from recovery are now allocated across all the GDP loss types – direct, indirect, and induced – whereas in version 4.1, all recovery gains were accounted for in the indirect loss type. To achieve this, we’ve introduced new methodology to streamline and simplify the calculation of all types of losses and recovery. In addition, RDEIM includes other kinds of losses, including tangible wealth. This includes loss of tangible assets (e.g., depreciation) and accident expenditures (e.g., decontamination). This document describes the updated RDEIM economic model and provides examples of loss and recovery calculation, results analysis, and presentation. Changes to the tangible cost calculation and accident expenditures are described in section 2.2. The updates to the RDEIM input-output (I-O) model are not expected to affect the final benchmark results Bixler et al. (2020), as the RDEIM calculation for the total national GDP losses remains unchanged. The reader is referred to the MACCS revision history for other cost modelling changes since version 4.0 that may affect the benchmark. RDEIM has its roots in a code developed by Sandia National Laboratories for the Department of Homeland Security to estimate short-term losses from natural and manmade accidents, called the Regional Economic Accounting analysis tool (REAcct). This model was adapted and modified for MACCS. It is based on I-O theory, which is widely used in economic modeling. It accounts for direct losses to a disrupted region affected by an accident, indirect losses to the national economy due to disruption of the supply chain, and induced losses from reduced spending by displaced workers. RDEIM differs from REAcct in in its treatment and estimation of indirect loss multipliers, elimination of double-counting associated with inter-industry trade in the affected area, and that it is intended to be used for extended periods that can occur from a major nuclear reactor accident, such as the one that occurred at the Fukushima Daiichi site in Japan. Most input-output models do not account for economic adaptation and recovery, and in this regard RDEIM differs from its parent, REAcct, because it allows for a user-definable national recovery period. Implementation of a recovery period was one of several recommendations made by an independent peer review panel to ensure that RDEIM is state-of-practice. For this and several other reasons, RDEIM differs from REAcct.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗