Search NASA⌕ Search

SEARCH · Search NASA

Results for “System Code”

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,765 records · Page 98

Parallelized real-time physics codes for plasma control on DIII-D

A real-time safe multi-threading library was developed on the DIII-D plasma control system to optimize the real-time TORBEAM and real-time STRIDE physics codes. These physics codes are crucial for future fusion power plant operation as they provide information about electron cyclotron wave propagation and heating as well as inform about ideal plasma stability limits. The real-time TORBEAM code executed consistently in under 20 ms while the real-time STRIDE code computes in 100 ms. The multi-threading library developed in this work can be applied to other real-time physics-based codes that will be crucial for the next generation of fusion devices.

DIII-D↗

Integrated State Estimation and Contingency Analysis Software Implementation using High Performance Computing Techniques

Power system simulation tools are traditionally developed in sequential mode and codes are optimized for single core computing only. However, the increasing complexity in the power grid models requires more intensive computation. The traditional simulation tools will soon not be able to meet the grid operation requirements. Therefore, power system simulation tools need to evolve accordingly to provide faster and better results for grid operations. This paper presents an integrated state estimation and contingency analysis software implementation using high performance computing techniques. The software is able to solve large size state estimation problems within one second and achieve a near-linear speedup of 9,800 with 10,000 cores for contingency analysis application. The performance evaluation is presented to show its effectiveness.

Chen, Yousu↗

The statistical analysis techniques to support the NGNP fuel performance experiments

This paper describes the development and application of statistical analysis techniques to support the Advanced Gas Reactor (AGR) experimental program on Next Generation Nuclear Plant (NGNP) fuel performance. The experiments conducted in the Idaho National Laboratory’s Advanced Test Reactor employ fuel compacts placed in a graphite cylinder shrouded by a steel capsule. The tests are instrumented with thermocouples embedded in graphite blocks and the target quantity (fuel temperature) is regulated by the He–Ne gas mixture that fills the gap volume. Three techniques for statistical analysis, namely control charting, correlation analysis, and regression analysis, are implemented in the NGNP Data Management and Analysis System for automated processing and qualification of the AGR measured data. The neutronic and thermal code simulation results are used for comparative scrutiny. The ultimate objective of this work includes (a) a multi-faceted system for data monitoring and data accuracy testing, (b) identification of possible modes of diagnostics deterioration and changes in experimental conditions, (c) qualification of data for use in code validation, and (d) identification and use of data trends to support effective control of test conditions with respect to the test target. Analysis results and examples given in the paper show the three statistical analysis techniques providing a complementary capability to warn of thermocouple failures. It also suggests that the regression analysis models relating calculated fuel temperatures and thermocouple readings can enable online regulation of experimental parameters (i.e. gas mixture content), to effectively maintain the fuel temperature within a given range.

99 GENERAL AND MISCELLANEOUS↗

Capturing Correlation Effects in Positron Binding to Atoms and Molecules

A major challenge in contemporary electronic structure theory involves the development of methods to describe in a balanced manner the contribution of correlation effects to energy differences. This challenge can be even greater for multicomponent systems containing more than one type of quantum particle. In the present work, we describe a flexible code for carrying out self-consistent field and configuration interaction (CI) calculations on multicomponent systems and use it to generate trial wave functions for use in diffusion Monte Carlo (DMC) calculations of the positron affinity of Be, Be 2 , Be 4 , Mg, CS 2 , and benzene. The resulting positron affinities (PAs) are in good agreement with the best values from the literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tritium Mitigation/Control for Advanced Reactor System

A tritium removal facility, which is similar to the design used for tritium recovery in fusion reactors, is proposed in this study for fluoride-salt-cooled high-temperature reactors (FHRs) to result in a two-loop FHR design with the elimination of an intermediate loop. Using this approach, an economic benefit can potentially be obtained by removing the intermediate loop, while the safety concern of tritium release can be mitigated. In addition, an intermediate heat exchanger (IHX) that can yield a similar tritium permeation rate to the production rate of 1.9 Ci/day in a 1,000 MWe PWR needs to be designed to prevent the residual tritium that is not captured in the tritium removal system from escaping into the power cycle and ultimately the environment. The main focus of this study is to aid the mitigation of tritium permeation issue from the FHR primary side to significantly reduce the concentration of tritium in the secondary side and the process heat application side (if applicable). The goal of the research is to propose a baseline FHR system without the intermediate loop. The specific objectives to accomplish the goals are: To estimate tritium permeation behavior in FHRs; To design a tritium removal system for FHRs; To meet the same tritium permeation level in FHRs as the tritium production rate of 1.9 Ci/day in 1,000 MWe PWRs; To demonstrate economic benefits of the proposed FHR system via comparing with the three-loop FHR system. The objectives were accomplished by designing tritium removal facilities, developing a tritium analysis code, and conducting an economic analysis. In the fusion reactor community, tritium extraction has been widely investigated and researched. Borrowing the experiences from the fusion reactor community, a tritium control and mitigation system was proposed. Based on mass transport theories, a tritium analysis code was developed, and the tritium behaviors were analyzed using the developed code. Tritium removal facilities were designed and laboratory-scale experiments were proposed for the validation of the proposed tritium removal facilities.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

CPMC-Lab: A Matlab package for Constrained Path Monte Carlo calculations

We describe CPMC-Lab, a Matlab program for the constrained-path and phaseless auxiliary-field Monte Carlo methods. These methods have allowed applications ranging from the study of strongly correlated models, such as the Hubbard model, to ab initio calculations in molecules and solids. The present package implements the full ground-state constrained-path Monte Carlo (CPMC) method in Matlab with a graphical interface, using the Hubbard model as an example. The package can perform calculations in finite supercells in any dimensions, under periodic or twist boundary conditions. Importance sampling and all other algorithmic details of a total energy calculation are included and illustrated. This open-source tool allows users to experiment with various model and run parameters and visualize the results. It provides a direct and interactive environment to learn the method and study the code with minimal overhead for setup. Furthermore, the package can be easily generalized for auxiliary-field quantum Monte Carlo (AFQMC) calculations in many other models for correlated electron systems, and can serve as a template for developing a production code for AFQMC total energy calculations in real materials. Several illustrative studies are carried out in one- and two-dimensional lattices on total energy, kinetic energy, potential energy, and charge- and spin-gaps.

Computer Science↗

Design-to-Deployment Continuum Platform for Microscopes and Computing Ecosystems

Science ecosystems with networked computing systems and physical instruments are increasingly being deployed with a goal to achieve the productivity promised by AI-supported remote automation. In support of these efforts, the virtual infrastructure twins (VITs) have been successfully utilized to develop the orchestration codes for these ecosystems without requiring physical access to expensive instruments, such as electron microscopes. Currently, the utility of such a VIT is severely limited by the computing capacity and capability of the computing system used as its host. Furthermore, codes developed on the VIT typically need to be transferred and refactored for production use, particularly, on high-performance systems with accelerators. In response, we develop a design-to-deployment continuum platform wherein a VIT runs natively on the ecosystem's own computing system, and thereby facilitates the continual in-situ testing and transition of codes for production use. Here, we describe the development and testing of software for remote microscope steering and GPU-based image reconstruction using this platform on a multi-GPU computing system networked to Nion microscopes. We demonstrate a continual transition of steering and reconstruction codes developed under VIT platform to production ecosystem deployment.

Al-Najjar, Anees [Oak Ridge National Laboratory (O↗

Programmable In Situ System for Iterative Workflows

We describe an in situ system for solving iterative problems. We specifically target inverse problems, where expensive simulations are approximated using a surrogate model. The model explores the parameter space of the simulation through iterative trials, each of which becomes a job managed by a parallel scheduler. Our work extends Henson [1], a cooperative multi-tasking system for in situ execution of loosely coupled codes.

Lohrmann, Erich↗

Critical review and analysis of hydrogen safety data collection tools

The wider adoption of hydrogen in multiple sectors of the economy requires that safety and risk issues be rigorously investigated. Quantitative Risk Assessment (QRA) is an important tool for enabling safe deployment of hydrogen fueling stations and is increasingly embedded in the permitting process. QRA requires reliability data, and currently hydrogen QRA is limited by the lack of hydrogen specific reliability data, thereby hindering the development of necessary safety codes and standards [1]. Four tools have been identified that collect hydrogen system safety data: H2Tools Lessons Learned, Hydrogen Incidents and Accidents Database (HIAD), National Renewable Energy Lab's (NREL) Composite Data Products (CDPs), and the Center for Hydrogen Safety (CHS) Equipment and Component Failure Rate Data Submission Form. This work critically reviews and analyzes these tools for their quality and usability in QRA. It is determined that these tools lay a good foundation, however, the data collected by these tools needs improvement for use in QRA. Areas in which these tools can be improved are highlighted, and can be used to develop a path towards adequate reliability data collection for hydrogen systems.

08 HYDROGEN↗

Petascale solvers for anisotropic PDEs in atmospheric modelling on GPU clusters

Memory bound applications such as solvers for large sparse systems of equations remain a challenge for GPUs. Fast solvers should be based on numerically efficient algorithms and implemented such that global memory access is minimised. To solve systems with trillions ( O ( 10 12 ) ) unknowns the code has to make efficient use of several million individual processor cores on large GPU clusters. We describe the multi-GPU implementation of two algorithmically optimal iterative solvers for anisotropic PDEs which are encountered in (semi-) implicit time stepping procedures in atmospheric modelling. In this application the condition number is large but independent of the grid resolution and both methods are asymptotically optimal, albeit with different absolute performance. In particular, an important constant in the discretisation is the CFL number; only the multigrid solver is robust to changes in this constant. We parallelise the solvers and adapt them to the specific features of GPU architectures, paying particular attention to efficient global memory access. We achieve a performance of up to 0.78 PFLOPs when solving an equation with 0.55 · 10 12 unknowns on 16384 GPUs; this corresponds to about 3% of the theoretical peak performance of the machine and we use more than 40% of the peak memory bandwidth with a Conjugate Gradient (CG) solver. Although the other solver, a geometric multigrid algorithm, has a slightly worse performance in terms of FLOPs per second, overall it is faster as it needs less iterations to converge; the multigrid algorithm can solve a linear PDE with half a trillion unknowns in about one second.

Computer Science↗

The HANDE-QMC Project: Open-Source Stochastic Quantum Chemistry from the Ground State Up

Building on the success of Quantum Monte Carlo techniques such as diffusion Monte Carlo, alternative stochastic approaches to solve electronic structure problems have emerged over the past decade. The full configuration interaction quantum Monte Carlo (FCIQMC) method allows one to systematically approach the exact solution of such problems, for cases where very high accuracy is desired. The introduction of FCIQMC has subsequently led to the development of coupled cluster Monte Carlo (CCMC) and density matrix quantum Monte Carlo (DMQMC), allowing stochastic sampling of the coupled cluster wave function and the exact thermal density matrix, respectively. In this Article, we describe the HANDE-QMC code, an open-source implementation of FCIQMC, CCMC and DMQMC, including initiator and semistochastic adaptations. We describe our code and demonstrate its use on three example systems; a molecule (nitric oxide), a model solid (the uniform electron gas), and a real solid (diamond). An illustrative tutorial is also included.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multilevel Morphology of Complex Nanoporous Materials

This work exploits gas adsorption and small-angle X-ray scattering (SAXS) to determine the morphology of complex nanoporous materials. We resolve multiple classes of porosity including previously undetected large-scale texture that significantly compromises the canonical interpretation of gas adsorption. Specifically, a UVM-7 class mesoporous silica was synthesized that has morphological features on three length scales: macropores due to packing of 150 nm globules, 1.9 nm radius spherical mesopores inside the globules, and >7 nm pockets on and between the globules. The total and external surface areas, as well as the mesopore volume, were determined using gas adsorption (α s -plot) and SAXS. A new approach was applied to the SAXS data using multilevel fitting to determine the surface areas on multiple length scales. The SAXS analysis code is applicable to any two-phase system and is freely available to the public. The total surface area determined by SAXS was 12% greater than that obtained by gas adsorption. The macropore interfacial area, however, is only 30% of the external surface area determined by the α s -plot. The overestimation of the external surface area by the αs-plot method is attributed to capillary condensation in nanoscale surface irregularities. Here, the discrepancy is resolved assuming that the macropore–globule interfaces harbor fractally distributed nooks and crannies, which lead to gas adsorption at pressures above the mesopore filling pressure.

36 MATERIALS SCIENCE↗

Numerical prediction of the thermodynamic properties of ternary Al-Ni-Hf alloys

Thermodynamic properties of ternary Al-Hf-Ni system, such as {sup ex}G, μ{sub Al}, μ{sub Ni} and μ{sub Zr} at 1373K were predicted on the basis of thermodynamic properties of binary systems included in the investigated ternary system. The idea of predicting {sup ex}G values was regarded as the calculation of excess Gibbs energy values inside a certain area (a Gibbs triangle) unless all boundary conditions, that is values of {sup ex}G on all legs of the triangle are known. {sup ex}G and L{sub ijk} ternary interaction parameters in the Muggianu extension of the Redlich-Kister formalism are calculated numerically using Wolfram Mathematica 9 software.

36 MATERIALS SCIENCE↗

Chemical insight from density functional modeling of molecular adsorption: Tracking the bonding and diffusion of anthracene derivatives on Cu(111) with molecular orbitals

We present a method of analyzing the results of density functional modeling of molecular adsorption in terms of an analogue of molecular orbitals. This approach permits intuitive chemical insight into the adsorption process. Applied to a set of anthracene derivates (anthracene, 9,10-anthraquinone, 9,10-dithioanthracene, and 9,10-diselenonanthracene), we follow the electronic states of the molecules that are involved in the bonding process and correlate them to both the molecular adsorption geometry and the species’ diffusive behavior. We additionally provide computational code to easily repeat this analysis on any system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Commissioning of a Varian Clinac iX 6 MV photon beam using Monte Carlo simulation

Monte Carlo modelling of a linear accelerator is the first and most important step in Monte Carlo dose calculations in radiotherapy. Monte Carlo is considered today to be the most accurate and detailed calculation method in different fields of medical physics. In this research, we developed a photon beam model for Varian Clinac iX 6 MV equipped with MilleniumMLC120 for dose calculation purposes using BEAMnrc/DOSXYZnrc Monte Carlo system based on the underlying EGSnrc particle transport code. Monte Carlo simulation for this commissioning head LINAC divided in two stages are design head Linac model using BEAMnrc, characterize this model using BEAMDP and analyze the difference between simulation and measurement data using DOSXYZnrc. In the first step, to reduce simulation time, a virtual treatment head LINAC was built in two parts (patient-dependent component and patient-independent component). The incident electron energy varied 6.1 MeV, 6.2 MeV and 6.3 MeV, 6.4 MeV, and 6.6 MeV and the FWHM (full width at half maximum) of source is 1 mm. Phase-space file from the virtual model characterized using BEAMDP. The results of MC calculations using DOSXYZnrc in water phantom are percent depth doses (PDDs) and beam profiles at depths 10 cm were compared with measurements. This process has been completed if the dose difference of measured and calculated relative depth-dose data along the central-axis and dose profile at depths 10 cm is ≤ 5%. The effect of beam width on percentage depth doses and beam profiles was studied. Results of the virtual model were in close agreement with measurements in incident energy electron 6.4 MeV. Our results showed that photon beam width could be tuned using large field beam profile at the depth of maximum dose. The Monte Carlo model developed in this study accurately represents the Varian Clinac iX with millennium MLC 120 leaf and can be used for reliable patient dose calculations. In this commissioning process, the good criteria of dose difference in PDD and dose profiles were achieve using incident electron energy 6.4 MeV.

61 RADIATION PROTECTION AND DOSIMETRY↗

An approach to model reactor core nodalization for deterministic safety analysis

Adopting good nodalization strategy is essential to produce an accurate and high quality input model for Deterministic Safety Analysis (DSA) using System Thermal-Hydraulic (SYS-TH) computer code. The purpose of such analysis is to demonstrate the compliance against regulatory requirements and to verify the behavior of the reactor during normal and accident conditions as it was originally designed. Numerous studies in the past have been devoted to the development of the nodalization strategy for small research reactor (e.g. 250kW) up to the bigger research reactor (e.g. 30MW). As such, this paper aims to discuss the state-of-arts thermal hydraulics channel to be employed in the nodalization for RTP-TRIGA Research Reactor specifically for the reactor core. At present, the required thermal-hydraulic parameters for reactor core, such as core geometrical data (length, coolant flow area, hydraulic diameters, and axial power profile) and material properties (including the UZrH{sub 1.6}, stainless steel clad, graphite reflector) have been collected, analyzed and consolidated in the Reference Database of RTP using standardized methodology, mainly derived from the available technical documentations. Based on the available information in the database, assumptions made on the nodalization approach and calculations performed will be discussed and presented. The development and identification of the thermal hydraulics channel for the reactor core will be implemented during the SYS-TH calculation using RELAP5-3D{sup ®} computer code. This activity presented in this paper is part of the development of overall nodalization description for RTP-TRIGA Research Reactor under the IAEA Norwegian Extra-Budgetary Programme (NOKEBP) mentoring project on Expertise Development through the Analysis of Reactor Thermal-Hydraulics for Malaysia, denoted as EARTH-M.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The Influence of the Hall Term on the Development of Magnetized Laser-Produced Plasma Jets

We present 2D axisymmetric simulation results describing the influence of the Hall term on laser-produced plasma jets and their interaction with an applied magnetic field parallel to the laser axis. Bending of the poloidal B-field lines produces an MHD shock structure surrounding a conical cavity, and a jet is produced from the convergence of the shock envelope. Both the jet and the conical cavity underneath it are bound by fast MHD shocks. We compare the MHD results generated using the extended-MHD code Physics as an Extended-MHD Relaxation System with an Efficient Upwind Scheme (PERSEUS) with MHD results generated using GORGON and find reasonable agreement. We then present extended-MHD results generated using PERSEUS, which show that the Hall term has several effects on the plasma jet evolution. A hot low-density current-carrying layer of plasma develops just outside the plume, which results in a helical rather than a purely poloidal B-field, and reduces magnetic stresses, resulting in delayed flow convergence and jet formation. The flow is partially frozen into the helical field, resulting in azimuthal rotation of the jet. The Hall term also produces field-aligned current in strongly magnetized regions. In particular, we find the influence of Hall physics on this problem to be scale-dependent. In conclusion, this points to the importance of mitigating the Hall effect in a laboratory setup, by increasing the jet density and system dimensions, in order to avoid inaccurate extrapolation to astrophysical scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Correlation consistent effective core potentials for late 3 d transition metals adapted for plane wave calculations

We construct a new modification of correlation consistent effective core potentials (ccECPs) for late 3d elements Cr–Zn with Ne-core that are adapted for efficiency and low energy cut-offs in plane wave calculations. The decrease in accuracy is rather minor, so that the constructions are in the same overall accuracy class as the original ccECPs. The resulting new constructions work with energy cut-offs at or below ≈400 Ry and, thus, make calculations of large systems with transition metals feasible for plane wave codes. We also provide the basic benchmarks for atomic spectra and molecular tests of this modified option that we denote as ccECP-soft.

Chemistry↗