Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer Codes”

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,315 records · Page 73

Effective field theory for dark matter absorption on single phonons

Single phonon excitations, with energies in the 1–100 meV range, are a powerful probe of light dark matter (DM). Utilizing effective field theory, we derive a framework to compute DM absorption rates into single phonons starting from general DM-electron, proton, and neutron interactions. We apply the framework to a variety of DM models: Yukawa coupled scalars, axionlike particles with derivative interactions, and vector DM coupling via gauge interactions or Standard Model electric and magnetic dipole moments. We find that GaAs or Al 2 O 3 targets can set powerful constraints on a U ( 1 ) B − L model, and targets with electronic spin ordering are similarly sensitive to DM coupling to the electron magnetic dipole moment. Lastly, we make the code, (an extension of the existing code which computes general DM–single phonon scattering rates), publicly available. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

SPACE code simulation of cold leg small break LOCA in the ATLAS integral test

SPACE code is a system analysis code for pressurized water reactors. This code uses a two-fluid and three-field model. For a few years, intensive validations have been performed to secure the prediction accuracy of models and correlations for two-phase flow and heat transfer. Recently, the code version 1.0 was released. This study is to see how well SPACE code predicts thermal hydraulic phenomena of an integral effect test. The target experiment is a cold leg small break LOCA in the ATLAS facility, which has the same two-loop features as APR1400. Predicted parameters were compared with experimental observations. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Barrier elision for production parallel programs

Large scientific code bases are often composed of several layers of runtime libraries, implemented in multiple programming languages. In such situation, programmers often choose conservative synchronization patterns leading to suboptimal performance. In this paper, we present context-sensitive dynamic optimizations that elide barriers redundant during the program execution. In our technique, we perform data race detection alongside the program to identify redundant barriers in their calling contexts; after an initial learning, we start eliding all future instances of barriers occurring in the same calling context. We present an automatic on-the-fly optimization and a multi-pass guided optimization. We apply our techniques to NWChem - a 6 million line computational chemistry code written in C/C++/Fortran that uses several runtime libraries such as Global Arrays, ComEx, DMAPP, and MPI. Our technique elides a surprisingly high fraction of barriers (as many as 63%) in production runs. This redundancy elimination translates to application speedups as high as 14% on 2048 cores. Our techniques also provided valuable insight about the application behavior, later used by NWChem developers. Overall, we demonstrate the value of holistic context-sensitive analyses that consider the domain science in conjunction with the associated runtime software stack.

97 MATHEMATICS AND COMPUTING↗

RAVEN User Manual

RAVEN is a generic software framework to perform parametric and probabilistic analysis based on the response of complex system codes. The initial development was aimed to provide dynamic risk analysis capabilities to the Thermo-Hydraulic code RELAP-7, currently under development at the Idaho National Laboratory (INL). Although the initial goal has been fully accomplished, RAVEN is now a multi-purpose probabilistic and uncertainty quantification platform, capable to agnostically communicate with any system code. This agnosticism includes providing Application Programming Interfaces (APIs). These APIs are used to allow RAVEN to interact with any code as long as all the parameters that need to be perturbed are accessible by inputs files or via python interfaces. RAVEN is capable of investigating the system response, and investigating the input space using Monte Carlo, Grid, or Latin Hyper Cube sampling schemes, but its strength is focused toward system feature discovery, such as limit surfaces, separating regions of the input space leading to system failure, using dynamic supervised learning techniques. The development of RAVEN has started in 2012, when, within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the need to provide a modern risk evaluation framework became stronger. RAVEN principal assignment is to provide the necessary software and algorithms in order to employ the concept developed by the Risk Informed Safety Margin Characterization (RISMC) program. RISMC is one of the pathways defined within the Light Water Reactor Sustainability (LWRS) program. In the RISMC approach, the goal is not just the individuation of the frequency of an event potentially leading to a system failure, but the closeness (or not) to key safety-related events. Hence, the approach is interested in identifying and increasing the safety margins related to those events. A safety margin is a numerical value quantifying the probability that a safety metric (e.g. for an important process such as peak pressure in a pipe) is exceeded under certain conditions. The initial development of RAVEN has been focused on providing dynamic risk assessment capability to RELAP-7, currently under development at the INL and, likely, future replacement of the RELAP5-3D code. Most the capabilities that have been implemented having RELAP-7 as principal focus are easily deployable for other system codes. For this reason, several side activaties are currently ongoing for coupling RAVEN with software such as RELAP5-3D, etc. The aim of this document is the explanation of the input requirements, focalizing on the input structure.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A Python Package for Fuel Cycle and Waste Management Analysis - 19545

The future of nuclear power utilization in the United States may involve a transition from the once-through light water reactor fuel cycle to a so-called advanced fuel cycle. An advanced fuel cycle may offer improvements in natural resource utilization, reactor safety, and environmental impact. In addition, adoption of an advanced fuel cycle will result in changes in the amount, composition, and disposal of the nuclear wastes they produce. Characteristics intrinsic to the wastes themselves, such as mass and radioactivity at specific times, were used as waste management metrics in the Department of Energy 2014 Fuel Cycle Evaluation and Screening (FCES) study. Binning of these metrics informed fuel cycle according to relative potential benefit. However, waste properties alone cannot qualify management strategy. Therefore, comparing fuel cycles from the perspective of waste management can provide additional insights. To that end, an object-oriented Python package has been developed to perform fuel cycle waste management calculations. The package, named nwpy (short for Nuclear Waste analysis in Python), mimics the mass flow operations reported in the FCES study and interfaces with the ORIGEN-S code to assess metric data as reported in the FCES. In addition, this Python package determines the loading of waste streams into forms suitable for final disposal in a geological repository and uses such data for repository performance assessment. This work describes the Python package methodology, assumptions, and capabilities. As input data, the package requires information about the fuel cycle such as the mass of fuel discharged, the cooling time between reactor discharge and subsequent operations, and the separations process efficiencies and waste forms it produces. In addition, the isotopic composition of the reactor discharge stream, such as that available from the FCES database, is required. With that, it can perform prescribed mass flow and decay operations throughout the back end of the fuel cycle, the latter carried out by the ORIGEN-S code distributed with the SCALE6.2 software suite. The result is time-dependent data characterizing the waste on a per-package basis. Package loading for waste streams depends on the stream form (i.e. UNF vs. HLW, or glass HLW from aqueous reprocessing vs. ceramic HLW from electrochemical processing) and is based on literature review and engineering assumption. As an example of the application of the code, calculations for the areal footprint of direct-contact repositories for different fuel cycles is shown. Development is ongoing, and the package is expected to grow incorporating capabilities for waste management analyses such as waste transportation and radionuclide transport through the biosphere. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Extension of mutation testing for the requirements and design faults - 322

Mutation Testing is a fault-based software testing technique that is considered to be the most efficient in detecting faults. It focuses on the software code level and follows a set of rules which are called mutation operators to systematically seed faults into the source code. Each seeded fault results is a new version of the software, which is called a mutant. Mutation testing promises to be effective in identifying adequate test data which can be used to find faults in the source code. However, traditional mutation testing is not geared towards the detection of requirements and design faults. To cover the full spectrum of possible faults, the mutation framework has been extended from the code level to the requirements and design level. Twenty-nine defects that appear in requirements and design documents were identified to define appropriate mutants. Then corresponding mutation operators that model the defects were developed. For each mutation operator, the strategy to generate a mutant and the evaluation of the number of possible mutants were defined. Existing code level cost reduction techniques were tailored for requirements and design mutation operators that result in a large number of mutants. With the extended mutation testing framework, the full spectrum of possible faults will be tested and identified to ensure the safety and functionality of the software. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

SU-F-J-72: A Clinical Usable Integrated Contouring Quality Evaluation Software for Radiotherapy

Purpose: To introduce the Auto Contour Evaluation (ACE) software, which is the clinical usable, user friendly, efficient and all-in-one toolbox for automatically identify common contouring errors in radiotherapy treatment planning using supervised machine learning techniques. Methods: ACE is developed with C# using Microsoft .Net framework and Windows Presentation Foundation (WPF) for elegant GUI design and smooth GUI transition animations through the integration of graphics engines and high dots per inch (DPI) settings on modern high resolution monitors. The industrial standard software design pattern, Model-View-ViewModel (MVVM) pattern, is chosen to be the major architecture of ACE for neat coding structure, deep modularization, easy maintainability and seamless communication with other clinical software. ACE consists of 1) a patient data importing module integrated with clinical patient database server, 2) a 2D DICOM image and RT structure simultaneously displaying module, 3) a 3D RT structure visualization module using Visualization Toolkit or VTK library and 4) a contour evaluation module using supervised pattern recognition algorithms to detect contouring errors and display detection results. ACE relies on supervised learning algorithms to handle all image processing and data processing jobs. Implementations of related algorithms are powered by Accord.Net scientific computing library for better efficiency and effectiveness. Results: ACE can take patient’s CT images and RT structures from commercial treatment planning software via direct user input or from patients’ database. All functionalities including 2D and 3D image visualization and RT contours error detection have been demonstrated with real clinical patient cases. Conclusion: ACE implements supervised learning algorithms and combines image processing and graphical visualization modules for RT contours verification. ACE has great potential for automated radiotherapy contouring quality verification. Structured with MVVM pattern, it is highly maintainable and extensible, and support smooth connections with other clinical software tools.

60 APPLIED LIFE SCIENCES↗

Yield Of Radioactive Products Of Transfer Reactions Induced By {sup 6}He Ions On {sup 197}Au Target

The excitation functions and isomeric cross section ratios of Tl isotopes obtained at irradiation of {sup 197}Au by a {sup 6}He-ions beams with energies from 15 to 60 MeV are analyzed. The EMPIRE-2.18 code which takes into account the branching ratios of gamma-transitions of discrete gamma-cascade in a residual nucleus in the calculations of the cross sections of the isomeric and ground nuclear states production is used.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Electron transfer in extended systems: characterization by periodic density functional theory including the electronic coupling

In this paper, we describe a new computer implementation of electron transfer (ET) theory in extended systems treated by periodic density functional theory (DFT), including the calculation of the electronic coupling transition element V AB . In particular, the development opens up the full characterization of electron transfer in the solid state. The approach is valid for any single-determinant wavefunction with localized character representing the electronic structure of the system, from Hartree–Fock (HF) theory, to density functional theory (DFT), hybrid DFT theory, DFT+U theory, and constrained DFT (cDFT) theory. The implementation in CP2K reuses the high-performance functions of the code. The computational cost is equivalent to only one iteration of an HF calculation. We present test calculations for electron transfer in a number of systems, including a 1D-model of ferric oxide, hematite Fe2O 3 , rutile TiO 2 , and finally bismuth vanadate BiVO 4 .

14 SOLAR ENERGY↗

Validation of a plant dynamics code for 4S - Test analysis of natural circulation behavior

A plant transient dynamics code for a sodium-cooled fast reactor was developed by Toshiba. The code is used to evaluate the safety performance of Super-Safe, Small, and Simple reactor (4S) for Anticipated Operational Occurrences (AOOs), Design Basis Accident (DBA) and Beyond DBA (BDBA). The code is currently undergoing verification and validation (V and V). As one of the validation, test analysis of the Shutdown Heat Removal Test (SHRT)-17 performed in the Experimental Breeder Reactor (EBR)-II was conducted. The SHRT-17 is protected loss of flow test. The purpose of this validation is to confirm capability of the code to simulate natural circulation behavior of the plant. As a result, good agreements are shown between the analytical results and the measured data which were available from instrumented subassembly. The detailed validation result of the natural circulation behavior is described in this paper. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

RADIATION TRANSPORT FOR EXPLOSIVE OUTFLOWS: OPACITY REGROUPING

Implicit Monte Carlo (IMC) and Discrete Diffusion Monte Carlo (DDMC) are methods used to stochastically solve the radiative transport and diffusion equations, respectively. These methods combine into a hybrid transport-diffusion method we refer to as IMC-DDMC. We explore a multigroup IMC-DDMC scheme that in DDMC, combines frequency groups with sufficient optical thickness. We term this procedure ''opacity regrouping''. Opacity regrouping has previously been applied to IMC-DDMC calculations for problems in which the dependence of the opacity on frequency is monotonic. We generalize opacity regrouping to non-contiguous groups and implement this in SuperNu, a code designed to do radiation transport in high-velocity outflows with non-monotonic opacities. We find that regrouping of non-contiguous opacity groups generally improves the speed of IMC-DDMC radiation transport. We present an asymptotic analysis that informs the nature of the Doppler shift in DDMC groups and summarize the derivation of the Gentile-Fleck factor for modified IMC-DDMC. We test SuperNu using numerical experiments including a quasi-manufactured analytic solution, a simple 10 group problem, and the W7 problem for Type Ia supernovae. We find that opacity regrouping is necessary to make our IMC-DDMC implementation feasible for the W7 problem and possibly Type Ia supernova simulations in general. We compare the bolometric light curves and spectra produced by the SuperNu and PHOENIX radiation transport codes for the W7 problem. The overall shape of the bolometric light curves are in good agreement, as are the spectra and their evolution with time. However, for the numerical specifications we considered, we find that the peak luminosity of the light curve calculated using SuperNu is ∼10% less than that calculated using PHOENIX.

79 ASTRONOMY AND ASTROPHYSICS↗

Code Verification Capabilities and Assessments in Support of ASC V&V Level 2 Milestone #6035

This document provides a summary of the code verification activities supporting the FY17 Level 2 V&V milestone entitled “Deliver a Capability for V&V Assessments of Code Implementations of Physics Models and Numerical Algorithms in Support of Future Predictive Capability Framework Pegposts.” The physics validation activities supporting this milestone are documented separately. The objectives of this portion of the milestone are: 1) Develop software tools to support code verification analysis; 2) Document standard definitions of code verification test problems; and 3) Perform code verification assessments (focusing on error behavior of algorithms). This report and a set of additional standalone documents serve as the compilation of results demonstrating accomplishment of these objectives.

97 MATHEMATICS AND COMPUTING↗

Verification of Monte Carlo transport codes against measured small angle p-, d-, and t-emission in carbon fragmentation at 600 MeV/nucleon

Momentum spectra of hydrogen isotopes have been measured at 3.5° from 12 C fragmentation on a Be target. Momentum spectra cover both the region of fragmentation maximum and the cumulative region. Differential cross sections span five orders of magnitude. The data are compared to predictions of four Monte Carlo codes: QMD, LAQGSM, BC, and INCL++. There are large differences between the data and predictions of some models in the high momentum region. The INCL++ code gives the best and almost perfect description of the data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Influence of Magnetically Conjugate Fragments of Auroral Emission Images on the Accuracy of Determining E{sub av} of Precipitating Electrons

Results of magnetic conjugation of image fragments of auroral emissions at different altitudes along the magnetic field lines and preliminary results of evaluation of their influence on the accuracy of remote mapping of energy characteristics of precipitating electrons are presented. The results are obtained using the code of tracing being an integral part of the software Vector M intended for calculation of accompanying, geophysical, and astronomical information for the center of mass of a space vehicle (SV) and remote observation of aurora by means of Aurovisor-VIS/MP imager onboard the SV Meteor-MP to be launched.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Calculation of electron-impact rotationally elastic total cross sections for NH{sub 3}, H{sub 2}S, and PH{sub 3} over the energy range from 0.01 eV to 2 keV

This paper report results of calculation of the total cross section Q{sub T} for electron impact on NH{sub 3}, H{sub 2}S, and PH{sub 3} over a wide range of incident energies from 0.01 eV to 2 keV. Total cross sections Q{sub T} (elastic plus electronic excitation) for incident energies below the ionization threshold of the target were calculated using the UK molecular R-matrix code through the Quantemol-N software package and cross sections at higher energies were derived using the spherical complex optical potential formalism. The two methods are found to give self-consistent values where they overlap. The present results are, in general, found to be in good agreement with previous experimental and theoretical results.

74 ATOMIC AND MOLECULAR PHYSICS↗

Spent fuel pool analysis using TRACE code

The storage requirements of Spent Fuel Pools have been analyzed with the purpose to increase their rack capacities. In the past, the thermal limits have been mainly evaluated with conservative codes developed for this purpose, although some works can be found in which a best estimate code is used. The use of best estimate codes is interesting as they provide more realistic calculations and they have the capability of analyzing a wide range of transients that could affect the Spent Fuel Pool. Two of the most representative thermal-hydraulic codes are RELAP-5 and TRAC. Nowadays, TRACE code is being developed to make use of the more favorable characteristics of RELAP-5 and TRAC codes. Among the components coded in TRACE that can be used to construct the model, it is interesting to use the VESSEL component, which has the capacity of reproducing three dimensional phenomena. In this work, a thermal-hydraulic model of the Maine Yankee spent fuel pool using the TRACE code is developed. Such model has been used to perform a licensing calculation and the results obtained have been compared with experimental measurements made at the pool, showing a good agreement between the calculations predicted by TRACE and the experimental data. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Software Quality Metrics for Geant4: An Initial Assessment - Paper 122

In the context of critical applications, such as shielding and radiation protection, ensuring the quality of simulation software they depend on is of utmost importance. The assessment of simulation software quality is crucial both to determine its adoption in experimental applications and to guarantee reproducibility of outcome over time. Geant4 is a simulation system used in a wide range of experimental applications, some of which, such as radiation protection and biomedical studies, may involve critical use cases and would especially profit from an objective assessment of the software quality. In this study, we present the first results from an ongoing analysis of Geant4 code based on established software metrics from the perspective of software engineering. The analysis evaluates the current status of the code to quantify its characteristics with respect to documented quality standards; further appraisals concern evolutions over a series of release distributions. The quantitative assessment of the software is initially focused on a subset of Geant4 packages, which play a key role in scientific applications and are representative of different software development processes. We describe the chosen software metrics that map the entities of a software system to numeric values through mathematical definitions, quantifying software attributes ranging from functionality to maintainability. To measure such metrics, we exploited a set of software metrics tools, i.e. programs that implement a set of software metrics definitions. They allow software engineers to assess a software system by extracting the required attributes from the software and providing the corresponding metric values. To correlate software metrics values in a well-defined way, we used software quality standard models, such as the ISO/IEC 9126 standard, which aggregate numerical values to aid quality analysis and assessment. With this work we aim to create a baseline for evaluating correlations between software quality embedded in the Geant4 development process and simulation observables produced by Geant4-based applications. Furthermore, we intend to highlight what metrics are most effective at evaluating radiation transport software quality. We are confident that the result will provide constructive guidance both to improve key software tools, such as Geant4, and to estimate their contribution in risk analysis concerning sensitive applications. The remainder of this paper is structured as follows. We first discuss the approach adopted to fulfil our assessment. Then, we present the setup of our experiments and the way we collected data. We provide an initial interpretation of a subset of data resulting from measurements made on the selected Geant4 packages. As the analysis is still ongoing at the time of submitting this abstract, some results will be documented in the final paper and in the conference presentation. With this work we aim to create a proper Geant4 data set for evaluating its software characteristics and identifying both strengths and weaknesses. In this initial study, we considered the geometry and processes packages, which play a major role in most simulation applications. The data analysis is in progress at the time of writing; its complexity is related to the differences among software metrics tools, to the large amount of data to be evaluated and to the pioneering character of this project. We are confident that the results will provide valuable indications for the Geant4 user community and constructive guidance to improve key Geant4 packages. Other Monte Carlo codes could profit from the methodology developed in this pilot project for similar assessments of their software quality. At the end of this study, we intend to create a baseline for evaluating correlations between software quality embedded in the Geant4 development process and simulation observables produced by Geant4-based applications. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

Modelling an induced thermal plume with data from electrical resistivity tomography and distributed temperature sensing: a case study in northeast Italy

Groundwater tracer tests are often used to improve aquifer characterization, but they present several disadvantages, such as the need to pour solutions or dyes into the aquifer system and alteration of the water’s chemical properties. Thus, tracers can affect the groundwater flow mechanics and data interpretation becomes more complex, hindering effective study of ground heat pumps for low enthalpy geothermal systems. This paper presents a preliminary methodology based on a multidisciplinary application of heat as a tracer for defining the main parameters of shallow aquifers. The field monitoring techniques electrical resistivity tomography (ERT) and distributed temperature sensing (DTS) are noninvasive and were applied to a shallow-aquifer test site in northeast Italy. The combination of these measurement techniques supports the definition of the main aquifer parameters and therefore the construction of a reliable conceptual model, which is then described through the numerical code FEFLOW. This model is calibrated with DTS and validated by ERT outcomes. The reliability of the numerical model in terms of fate and transport is thereby enhanced, leading to the potential for better environmental management and protection of groundwater resources through more cost-effective solutions.

54 ENVIRONMENTAL SCIENCES↗