Search NASA⌕ Search

SEARCH · Search NASA

Results for “computer Codes Computer Programs”

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,117 records · Page 62

NEAMS Workbench Status and Capabilities

The mission of the US Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program is to develop, apply, deploy, and support state-of-the-art predictive modeling and simulation tools for the design and analysis of current and future nuclear energy systems. This is accomplished by using computing architectures that range from laptops to leadership-class facilities. The NEAMS Workbench is a new initiative that will facilitate the transition from conventional tools to high fidelity tools by providing a common user interface for model creation, review, execution, output review, and visualization for integrated codes. The Workbench can use common user input, including engineering-scale specifications that are expanded into application-specific input requirements through the use of customizable templates. The templating process can enable multifidelity analysis of a system from a common set of input data. Additionally, the common user input processor can provide an enhanced alternative application input that provides additional conveniences over the native input, especially for legacy codes. Expansion of the integrated codes and application templates available in the Workbench will broaden the NEAMS user community and will facilitate system analysis and design. Current and planned capabilities of the NEAMS Workbench are detailed here.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Uncertainty and Sensitivity Analysis of TREAT Transient Test #15

The U.S. Department of Energy (DOE) is preparing for the resumption of transient testing using the Transient Reactor Test Facility (TREAT) located at the Idaho National Laboratory (INL). In order to support the analysis of TREAT operation and experiment planning, the MAMMOTH reactor physics team has recently begun development of uncertainty and sensitivity analysis for the TREAT core simulated via MAMMOTH. MAMMOTH is a high-fidelity three-dimensional transient multi-physics code developed at INL under the DOE-NE Advanced Modeling and Simulation program (NEAMS) [1]. This research is part of the development presenting some initial uncertainty and sensitivity analysis of TREAT Transient Test #15 [2]. In this summary, we employ the Sobol sampling strategy implemented in RAVEN [3] to efficiently perform the sensitivity and uncertainty analysis. The idea is using Sobol functional decomposition with generalized polynomial chaos (gPC) expansions to estimate the responses of interest [4]. The gPC coefficients are then computed via the spectral projection method with sparse grid approach. In addition, this strategy includes evaluation of the second or higher order sensitivity indices, allowing the study of interactions between several input variables. The principal focus of this analysis is to quantify the uncertainties in the total power and maximum fuel temperatures with respect to uncertainties in the initial transient conditions, such as rated power, fuel temperature, graphite fuel density and control rod boron concentration. In addition, both the first order and second order time-dependent sensitivity indices are computed and analyzed. The results show that the interactions between initial control rod boron concentration and the rated power will contribute at most 1% to the total uncertainty of the power peak.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

VTK-m User's' Guide (V.1.7)

High-performance computing relies on ever finer threading. Advances in processor technology include ever greater numbers of cores, hyperthreading, accelerators with integrated blocks of cores, and special vectorized instructions, all of which require more software parallelism to achieve peak performance. Traditional visualization solutions cannot support this extreme level of concurrency. Extreme scale systems require a new programming model and a fundamental change in how we design algorithms. To address these issues we created VTK-m: the visualization toolkit for multi-/many-core architectures. VTK-m supports a number of algorithms and the ability to design further algorithms through a top-down design with an emphasis on extreme parallelism. VTK-m also provides support for finding and building links across topologies, making it possible to perform operations that determine manifold surfaces, interpolate generated values, and find adjacencies. Although VTK-m provides a simplified high-level interface for programming, its template-based code removes the overhead of abstraction. VTK-m simplifies the development of parallel scientific visualization algorithms by providing a framework of supporting functionality that allows developers to focus on visualization operations. Consider the listings in Figure 1.1 that compares the size of the implementation for the Marching Cubes algorithm in VTK-m with the equivalent reference implementation in the CUDA software development kit. Because VTK-m internally manages the parallel distribution of work and data, the VTK-m implementation is shorter and easier to maintain. Additionally, VTK-m provides data abstractions not provided by other libraries that make code written in VTK-m more versatile.This book includes contributions from the VTK-m community including the VTK-m development team and the user community.

97 MATHEMATICS AND COMPUTING↗

TH-AB-207A-07: Radiation Dose Simulation for a Newly Proposed Dynamic Bowtie Filters for CT Using Fast Monte Carlo Methods

Purpose: Dynamic bowtie filter is an innovative design capable of modulating the X-ray and balancing the flux in the detectors, and it introduces a new way of patient-specific CT scan optimizations. This study demonstrates the feasibility of performing fast Monte Carlo dose calculation for a type of dynamic bowtie filter for cone-beam CT (Liu et al. 2014 9(7) PloS one) using MIC coprocessors. Methods: The dynamic bowtie filter in question consists of a highly attenuating bowtie component (HB) and a weakly attenuating bowtie (WB). The HB is filled with CeCl3 solution and its surface is defined by a transcendental equation. The WB is an elliptical cylinder filled with air and immersed in the HB. As the scanner rotates, the orientation of WB remains the same with the static patient. In our Monte Carlo simulation, the HB was approximated by 576 boxes. The phantom was a voxelized elliptical cylinder composed of PMMA and surrounded by air (44cm×44cm×40cm, 1000×1000×1 voxels). The dose to the PMMA phantom was tallied with 0.15% statistical uncertainty under 100 kVp source. Two Monte Carlo codes ARCHER and MCNP-6.1 were compared. Both used double-precision. Compiler flags that may trade accuracy for speed were avoided. Results: The wall time of the simulation was 25.4 seconds by ARCHER on a 5110P MIC, 40 seconds on a X5650 CPU, and 523 seconds by the multithreaded MCNP on the same CPU. The high performance of ARCHER is attributed to the parameterized geometry and vectorization of the program hotspots. Conclusion: The dynamic bowtie filter modeled in this study is able to effectively reduce the dynamic range of the detected signals for the photon-counting detectors. With appropriate software optimization methods, the accelerator-based (MIC and GPU) Monte Carlo dose engines have shown good performance and can contribute to patient-specific CT scan optimizations.

60 APPLIED LIFE SCIENCES↗

RISMC Industry Application no. 1 (ECCS/LOCA) - Core Characterization Automation: Lattice Codes Interface for PHISICS/RELAP5-3D

The U.S. NRC is currently proposing rulemaking '10 CFR 50.46c' to revise the loss-of-coolant-accident (LOCA) and emergency core cooling system (ECCS) acceptance criteria to include the effects of higher burn-up on cladding performance. The key implications of this proposition are that the cladding performance cannot be evaluated in isolation anymore. Both cladding and ECCS performance need to be considered in a coupled manner. This also implies that models for cladding performance as well as LOCA methodologies will need to be updated for regulatory purposes. Furthermore, efforts need to be made to find the limiting case within the regulatory input space. The Risk-Informed Safety Margin Characterization (RISMC) Pathway, as part of the DOE Light Water Reactor Sustainability (LWRS) research and development program initialed several programs at the Idaho National Laboratory (INL) to support the industry in the transition to the new rule and to offer potential solutions to this LOCA/ECCS analysis. One of these programs is the Industry Application no. 1 (IA1). Its purpose is to provide to the plant owner/operator a vehicle to inform decisions to manage margins related to compliance with the proposed 10 CFR 50.46c rule. IA1's goal is to develop an Integrated Evaluation Model (IEM) to revisit how uncertainties are propagated across the physical disciplines and data involved, as well as well as how risks are evaluated in a LOCA under 10 CFR 50.46c. This IEM is called LOCA Tool-kit for the US (LOTUS) and it connects five major disciplines involved in LOCA analysis. The focus of LOTUS is to establish the automation interfaces among the five disciplines. LOTUS will utilize in a first step existing computer codes. The risk-informed margins management approach for IA1 will provide a means of quantifying the impact on the key LOCA analysis figures-of-merit like peak cladding temperature (PCT), equivalent cladding reacted (ECR), etc. of a change in LOCA analysis inputs. The information that the risk analysis and associated tools provide can then be used for decision-making and margin management. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

TNSL improvements in FUDGE [Slides]

FUDGE TNSL processing capabilities are improving. New development focuses on reading and writing GNDS, but FUDGE also supports writing to ACE and ENDL. Work provides a chance to compare in-depth with other codes. TNSL processing is mostly complete, but MCGIDI sampling has room for improvement. There are plans to revisit interpolation and to improve strategies for sampling coherent/incoherent elastic.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

SCALE Activities in FY23 [Slides]

This presentation shows that the SCALE 6.3 is available from RSICC. The lecture touches on new features including the new ENDF/B-VIII.0 data including covariances. It shows that the updated parallel infrastructure enables parallel capability on Windows. Production release with maintenance until 2026 at minimum addressing Code or data bugs, Performance, Ease of installation. There is a New Government Use Agreement (GUA) for SCALE 7.0 beta access. With site licenses available for non-commercial testing and feedback, handled through ORNL technology transfer.

97 MATHEMATICS AND COMPUTING↗

Investigating Fission Reaction Rate Ratio Sensitivities [Abstract]

Reaction rate ratios are a measurable parameter for reactor and criticality applications. A number of foil irradiations and fission chamber measurements have been performed for critical assemblies at Los Alamos National Laboratory starting in the 1950’s including (i) Godiva, a bare HEU spherical assembly; (ii) Flattop-25, a spherical assembly consisting of an HEU core and a natural uranium reflector; (iii) Jezebel, a bare 239 Pu assembly; and (iv) Flattop-Pu, a spherical assembly consisting of a 239Pu core and a natural uranium reflector. Fission ratio data for 238 U(n,f)/ 235 U(n,f), 237 Np(n,f)/ 235 U(n,f), 233 U(n,f)/ 235 U(n,f) and 239 Pu(n,f)/ 235 U(n,f) were obtained and reported. The EUCLID (Experiments Underpinned by Computational Learning for Improvements in nuclear Data) project at Los Alamos National Laboratory (LANL) aims to constrain nuclear data by using a suite of measurement types beyond k-effective. Recent investigations include the use of pulsed spheres for nuclear data validation and other measurement methods of interest. One focus of the work is to determine if other methods are complimentary to the critical experiments utilized for nuclear data validation. It is anticipated the investigations will help inform methods that may be utilized in machine learning algorithms for nuclear validation. In order to use a measurement type for nuclear validation, it is necessary to obtain cross-section sensitivities for parameters. This work looks at reaction rate ratio sensitivities with SENSMG and Monte Carlo N-Particle R Code Version 6.21.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Progress and status of the OpenMC Monte Carlo code

The present work describes the latest advances and progress in the development of the OpenMC Monte Carlo code, an open-source code originating from the Massachusetts Institute of Technology. First, an overview of the development workflow of OpenMC is given. Various enhancements to the code such as real-time XML input validation, state points, plotting, OpenMP threading, and coarse mesh finite difference acceleration are described. (authors)

97 MATHEMATICS AND COMPUTING↗

Development of the computer code for multi-calculations of source term of spent fuel - 15351

Many PWR spent fuels have been transferred between the neighboring nuclear power plants at Kori, Hanbit, Hanwul site, and from nuclear power plant to the KAERI(Korea Atomic Energy Research Institute) for the research purpose in Korea. In case of the PHWR type spent fuels, these have been transferred from the spent fuel pool to an on-site dry storage facility. Total radioactivity of spent fuels has to be calculated before transport in advance by the regulation requirement. In case of PWR, KN-12 transport cask can transfer 12 spent fuel assemblies at one time and in case of PHWR, Hi-STAR 63 transport cask can transfer 120 spent fuel assemblies at one time. Therefore, total radioactivity of all transported spent fuel assemblies that are transported generally are calculated by using the ORIGEN-ARP to meet this regulatory requirement. Basically, ORIGEN-ARP calculates the one fuel assembly at one time. SCALE code including this code provides the batch mode to perform the multi-calculations of the source term. In this study, in order to use the batch mode of ORIGEN-ARP, ARP-auto program was developed to generate input files and batch file of ORIGEN-ARP automatically and to make the summary file extracting the result from output files of ORIGEN-ARP. This code enhances the efficiency at the field. The working time to complete the calculation of source terms of 120 spent fuel assemblies take only 20 minutes. In order to confirm the credibility of this code, total radioactivity of 12 spent fuel assemblies of 17x17 Westinghouse type were calculated and compared to the previous result. The initial enrichments of these spent fuels are 3.46 ∼ 3.50 w/o and burnups are 29 ∼ 43 GWD/MTU. As the results, the relative errors are about 0.9 ∼ 3.0%. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Global Transformations for Legacy Parallel Applications via Structural Analysis and Rewriting

Performance and scalability optimization of large HPC applications is currently a labor-intensive, manual process with very low productivity. Major difficulties come from the disaggregated environment for HPC application development: the compiler is only involved in local decisions (core or multithreaded domain), while a library-based, communication-oriented programming model realizes whole-machine parallelism. Realizing any major global change in such a disaggregated environment is very difficult and involves changing large portions of the source code. We present semi-automated techniques, based on structural analysis and rewriting, for performing global transformations on an HPC application source code. We present two case studies using the Self-Consistent Field (SCF) standalone benchmark as well as the Coupled Cluster (CCSD) module (2.9 million lines of Fortran code), a key module of the NWChem computational chemistry application. We demonstrate how structural rewriting techniques can be used to automate transformations that affect multiple sections of the application's source code. We show that the transformations can be applied in a systematic fashion across the source code bases with minimal manual effort. These transformations improve the scalability of the SCF benchmark by more than two orders of magnitude and the performance of the full CCSD module by a factor of 4.0.

legacy parallel applications↗

STEPS: A Portable Numerical Simulation Toolkit for Electrical Power System Dynamic Studies

Numerical simulation is the key technique for large scale power system analysis. Redistribution of global renewable power via international interconnections requires new simulation tools to study the interconnected systems with different nominal frequencies as a whole. In this paper we introduce an open source simulation toolkit for electrical power systems (STEPS) which is hosted at Github. Its kernel is coded in C++ with major functions of power flow and electro-mechanical dynamic simulation. Flexible options are provided and configurable to improve power flow solution and dynamic simulation. Common devices and models are supported in STEPS for AC/DC hybrid system studies. Studies of interconnected systems with different nominal frequencies is supported in STEPS for research of international interconnection. Application program interfaces are provided and wrapped with Python to enable high-level interfaces for general applications. STEPS is thread safe and parallel computation is supported in both kernel and script levels to accelerate simulation. It is portable and works on Windows and GNU/Linux platforms. Cases from small to large scale systems are thoroughly tested to validate the toolkit with commercial packages as benchmarks.

42 ENGINEERING↗

An infeasible-start framework for convex quadratic optimization, with application to constraint-reduced interior-point and other methods

A framework is proposed for solving general convex quadratic programs (CQPs) from an infeasible starting point by invoking an existing feasible-start algorithm tailored for inequality-constrained CQPs. The central tool is an exact penalty function scheme equipped with a penalty-parameter updating rule. The feasible-start algorithm merely has to satisfy certain general requirements, and so is the updating rule. Under mild assumptions, the framework is proved to converge on CQPs with both inequality and equality constraints and, at a negligible additional cost per iteration, produces an infeasibility certificate, together with a feasible point for an (approximately) ℓ 1 -least relaxed feasible problem, when the given problem does not have a feasible solution. The framework is applied to a feasible-start constraint-reduced interior-point algorithm previously proved to be highly performant on problems with many more inequality constraints than variables (“imbalanced”). Numerical comparison with popular codes (OSQP, qpOASES, MOSEK) is reported on both randomly generated problems and support-vector machine classifier training problems. The results show that the former typically outperforms the latter on imbalanced problems. Finally, application of the proposed infeasible-start framework to other feasible-start algorithms is briefly considered, and is tested on a simplex iteration.

97 MATHEMATICS AND COMPUTING↗

A Review of Lightweight Thread Approaches for High Performance Computing

High-level, directive-based solutions are becoming the programming models (PMs) of the multi/many-core architectures. Several solutions relying on operating system (OS) threads perfectly work with a moderate number of cores. However, exascale systems will spawn hundreds of thousands of threads in order to exploit their massive parallel architectures and thus conventional OS threads are too heavy for that purpose. Several lightweight thread (LWT) libraries have recently appeared offering lighter mechanisms to tackle massive concurrency. In order to examine the suitability of LWTs in high-level runtimes, we develop a set of microbenchmarks consisting of commonlyfound patterns in current parallel codes. Moreover, we study the semantics offered by some LWT libraries in order to expose the similarities between different LWT application programming interfaces. This study reveals that a reduced set of LWT functions can be sufficient to cover the common parallel code patterns and that those LWT libraries perform better than OS threads-based solutions in cases where task and nested parallelism are becoming more popular with new architectures.

Castello, Adrian↗

Plasma Surface Interactions: Bridging from the Surface to the Micron Frontier through Leadership Class Computing (Final Technical Progress Report)

The performance demands on plasma-facing components (PFCs) in future fusion power plants are beyond the capability of current materials. In defining the plasma-surface interaction (PSI), it is evident that three coupled spatial regions influence PFC materials evolution and performance. These regions consist of (1) the edge and scrape-off layer region of the plasma, (2) the near-surface material response to extreme thermal and particle fluxes under the influence of, and feedback to, the plasma sheath, and (3) the structural material’s response to an intense, 14 MeV peaked neutron spectrum, which produces very high concentrations of transmuted elements through (n,p) and (n,α) reactions and structural material property degradation. The objective of this project is to develop and deploy validated, high-performance simulation tools capable of predicting the performance of tungsten-based plasma-facing components (PFCs) in a burning fusion plasma environment, which includes modeling surface morphology evolution in either erosion or re-deposition regimes, and the recycling of hydrogen species. This requires the development of a leadership-scale computational code to predict PFC behavior that is well integrated to a suite of multiscale modeling techniques to bridge the scales needed to address complex physical and computational issues at the plasma surface interface and the transition region below the surface where neutron damage processes in the bulk dominate material behavior. Successful completion of this project will provide PFC simulation tools to evaluate steady-state performance of tungsten-based PFC and divertor components in burning plasma environments. This will enable the identification of critical experiments to confirm whether practical PFC solutions exist for magnetic fusion energy beyond ITER. More specifically, the research activities within this project focus on two broadly defined research thrusts: 1) Developing a new simulation code, Xolotl-PSI to predict PFC operating lifetime and performance, that is specifically designed to take advantage of leadership class computing facilities; and 2) Integrating and applying discrete particle-based, as well as continuum-based, multiscale modeling techniques to provide scientific discovery of the mechanisms controlling PFC and bulk materials evolution under fusion plasma and 14-MeV neutron exposure; in this thrust, the challenge of scale-bridging between atomistic/microstructural modeling and the continuum-scale of PFC device performance is achieved through multiscale modeling techniques. Throughout the grant period, research at UMass Amherst focused on the atomic-scale analysis of the dynamics and kinetics (transport and reactions) of small mobile helium clusters near surfaces and grain boundaries in plasma-exposed tungsten, as well as the proper parameterization and incorporation of these processes into our SciDAC team’s continuum-level cluster dynamics code (Xolotl). Our work also included characterization of large-scale molecular-dynamics (MD) simulation results of helium implantation into tungsten based on the findings of our analysis and comparisons of Xolotl simulation predictions with the large-scale MD simulation results for Xolotl validation purposes; emphasis was placed on the impact of helium cluster dynamics/kinetics on tungsten surface morphology, near-surface tungsten structure, and helium retention. Furthermore, our work included initiation of a research program on evaluating the impact of structural defects in tungsten due to plasma exposure (such as helium nanobubbles) on tungsten’s thermophysical properties, focusing on lattice thermal conductivity. Our studies contributed to a fundamental understanding of the effects of plasma-surface interactions on the dynamical response of plasma-facing materials in nuclear fusion devices and the development of predictive hierarchical multiscale computational tools for the quantitative description of the material’s dynamical response to plasma exposure.

36 MATERIALS SCIENCE↗

On a stochastic approach to a code performance estimation

The main goal of an efficient profiling of software is to minimize the runtime overhead under certain constraints and requirements. The traces built by a profiler during the work, affect the performance of the system itself. One of important aspect of an overhead arises from the randomness of variability in the context in which the application is embedded, e.g., due to possible cache misses, etc. Such uncertainty needs to be taken into account in the design phase. In order to overcome these difficulties we propose to investigate this issue through the analysis of the probability distribution of the difference between profiler’s times for the same code. The approximating model is based on the finite normal mixtures within the framework of the method of moving separation of mixtures. We demonstrate some results for the MATLAB profiler using plotting of 3D surfaces by the function surf. The idea can be used for an estimating of a program efficiency.

97 MATHEMATICS AND COMPUTING↗

Experimental Characterization of OpenMP Offloading Memory Operations and Unified Shared Memory Support

The OpenMP specification recently introduced support for unified shared memory, allowing implementation to leverage underlying system software to provide a simpler GPU offloading model where explicit mapping of variables is optional. Support for this feature is becoming more available in different OpenMP implementations on several hardware platforms. A deeper understanding of the different implementation’s execution profile and performance is crucial for applications as they consider the performance portability implications of adopting a unified memory offloading programming style. This work introduces a benchmark tool to characterize unified memory support in several OepnMP compilers and runtimes, with emphasis on identifying discrepancies between different OpenMP implementations as to how they various memory allocation strategies interact with unified shared memory. The benchmark tool is used to characterize OpenMP compilers on three leading High Performance Computing platforms supporting different CPU and device architectures. The benchmark tool is used to assess the impact of enabling unified shared memory on the performance of memory-bound code, highlighting implementation differences that should be accounted for when applications consider performance portability across platforms and compilers.

Elwasif, Wael↗

OpenMP application experiences: Porting to accelerated nodes

As recent enhancements to the OpenMP specification become available in its implementations, there is a need to share the results of experimentation in order to better understand the OpenMP implementation’s behavior in practice, to identify pitfalls, and to learn how the implementations can be effectively deployed in scientific codes. We report on experiences gained and practices adopted when using OpenMP to port a variety of ECP applications, mini-apps and libraries based on different computational motifs to accelerator-based leadership-class high-performance supercomputer systems at the United States Department of Energy. Additionally, we identify important challenges and open problems related to the deployment of OpenMP. Through our report of experiences, we find that OpenMP implementations are successful on current supercomputing platforms and that OpenMP is a promising programming model to use for applications to be run on emerging and future platforms with accelerated nodes.

97 MATHEMATICS AND COMPUTING↗