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,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

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↗
DOE OSTI2014

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↗
DOE OSTI2025

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↗
DOE OSTI2018

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↗
DOE OSTI2022

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↗
DOE OSTI2015

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↗
DOE OSTI2019

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↗
DOE OSTI2018

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↗
DOE OSTI2015

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↗
DOE OSTI2015

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↗
DOE OSTI2016

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↗
DOE OSTI2018

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↗
DOE OSTI2022

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↗
DOE OSTI2024

Mixed Stochastic-Deterministic Approach for Many-Body Perturbation Theory Calculations

Here, we present an approach for GW calculations of quasiparticle energies with quasi-quadratic scaling by approximating high-energy contributions to the Green's function in its Lehmann representation with effective stochastic vectors. The method is easy to implement without altering the GW code, converges rapidly with stochastic parameters, and treats systems of various dimensionality and screening response. Our calculations on a 5.75° twisted MoS 2 bilayer show how large-scale GW methods include geometry relaxations and electronic correlations on an equal basis in structurally nontrivial materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

PyOMP: Multithreaded Parallel Programming in Python

We know that Python is a widely used language in scientific computing. When the goal is high performance, however, Python lags far behind low-level languages such as C and Fortran. To support applications that stress performance, Python needs to access the full capabilities of modern CPUs. That means support for parallel multithreading. In this paper, we describe PyOMP, a system that enables OpenMP in Python. Programmers write code in Python with OpenMP, Numba generates code that compiles to LLVM, and the resulting programs run with performance that approaches that from code written with C and OpenMP. In this paper we provide an update on the PyOMP project and explain how to install it and use it to write parallel multithreaded code in Python.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

GPU Lifetimes on Titan Supercomputer: Survival Analysis and Reliability

The Cray XK7 Titan was the top supercomputer system in the world for a long time and remained critically important throughout its nearly seven year life. It was an interesting machine from a reliability viewpoint as most of its power came from 18,688 GPUs whose operation was forced to execute three rework cycles, two on the GPU mechanical assembly and one on the GPU circuitboards. We write about the last rework cycle and a reliability analysis of over 100,000 years of GPU lifetimes during Titan’s 6-year-long productive period. Using time between failures analysis and statistical survival analysis techniques, we find that GPU reliability is dependent on heat dissipation to an extent that strongly correlates with detailed nuances of the cooling architecture and job scheduling. We describe the history, data collection, cleaning, and analysis and give recommendations for future supercomputing systems. We make the data and our analysis codes publicly available.

Ostrouchov, George↗
DOE OSTI2020

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↗
DOE OSTI2016

Endoluminal ultrasound applicators for MR-guided thermal ablation of pancreatic tumors: Preliminary design and evaluation in a porcine pancreas model

Purpose: Endoluminal ultrasound may serve as a minimally invasive option for delivering thermal ablation to pancreatic tumors adjacent to the stomach or duodenum. The objective of this study was to explore the basic feasibility of this treatment strategy through the design, characterization, and evaluation of proof-of-concept endoluminal ultrasound applicators capable of placement in the gastrointestinal (GI) lumen for volumetric pancreas ablation under MR guidance. Methods: Two variants of the endoluminal applicator, each containing a distinct array of two independently powered transducers (10 × 10 mm 3.2 MHz planar; or 8 × 10 × 20 mm radius of curvature 3.3 MHz curvilinear geometries) at the distal end of a meter long flexible catheter assembly, were designed and fabricated. Transducers and circulatory water flow for acoustic coupling and luminal cooling were contained by a low-profile polyester balloon covering the transducer assembly fixture. Each applicator incorporated miniature spiral MR coils and mechanical features (guiding tips and hinges) to facilitate tracking and insertion through the GI tract under MRI guidance. Acoustic characterization of each device was performed using radiation force balance and hydrophone measurements. Device delivery into the upper GI tract, adjacent to the pancreas, and heating characteristics for treatment of pancreatic tissue were evaluated in MR-guided ex vivo and in vivo porcine experiments. MR guidance was utilized for anatomical target identification, tracking/positioning of the applicator, and MR temperature imaging (MRTI) for PRF-based multislice thermometry, implemented in the real-time RTHawk software environment. Results: Force balance and hydrophone measurements indicated efficiencies of 48.8% and 47.8% and −3 dB intensity beam-widths of 3.2 and 1.2 mm for the planar and curvilinear transducers, respectively. Ex vivo studies on whole-porcine carcasses revealed capabilities of producing ablative temperature rise (ΔT > 15 °C) contours in pancreatic tissue 4–40 mm long and 4–28 mm wide for the planar transducer applicator (1–13 min sonication duration, ∼4 W/cm{sup 2} applied acoustic intensity). Curvilinear transducers produced more selective heating, with a narrower ΔT > 15 °C contour length and width of up to 1–24 mm and 2–7 mm, respectively (1–7 min sonication duration, ∼4 W/cm{sup 2} applied acoustic intensity). Active tracking of the miniature spiral coils was achieved using a Hadamard encoding tracking sequence, enabling real-time determination of each coil’s coordinates and automated prescription of imaging planes for thermometry. In vivo MRTI-guided heating trials in three pigs demonstrated capability of ∼20 °C temperature elevation in pancreatic tissue at 2 cm depths from the applicator, with 5–7 W/cm{sup 2} applied intensity and 6–16 min sonication duration. Dimensions of thermal lesions in the pancreas ranged from 12 to 28 mm, 3 to 10 mm, and 5 to 10 mm in length, width, and depth, respectively, as verified through histological analysis of tissue sections. Multiple-baseline reconstruction and respiratory-gated acquisition were demonstrated to be effective strategies in suppressing motion artifacts for clear evolution of temperature profiles during MRTI in the in vivo studies. Conclusions: This study demonstrates the technical feasibility of generating volumetric ablation in pancreatic tissue using endoluminal ultrasound applicators positioned in the stomach lumen. MR guidance facilitates target identification, device tracking/positioning, and treatment monitoring through real-time multislice PRF-based thermometry.

60 APPLIED LIFE SCIENCES↗