Search NASA⌕ Search

SEARCH · Search NASA

Results for “C 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 811 records · Page 45

Comparative study of decay heat calculations with FLUKA and MCNP/CINDER2008

In designing a high-power tungsten target, decay heat driven temperature rise in the spallation volume is a safety concern during maintenance and in loss of coolant accidents. As tungsten hydrates and becomes volatile in steam at above 800 ◦C, it is important to keep the target temperature below this threshold when active cooling is unavailable. Decay heat in a target is calculated with particle transport simulation codes combined with transmutation codes. The calculated decay heat usually differs depending on the nuclear cross sections and the decay particle transport models built in the code architecture. In this paper, we calculated decay heat of a water-cooled tungsten target using popular particle transport codes, FLUKA and MCNP6® paired with CINDER2008 and compared the results. The target-moderator-reflector (TMR) system is modeled with a water-cooled solid tungsten target, water premoderators, liquid hydrogen cold moderators and beryllium reflectors. Water-cooled stainless-steel shielding is modeled around the TMR system. The tungsten volume is clad with a thin layer of erosion/corrosion resistant material. This study provides information about the uncertainty range in decay heat prediction of high-power spallation targets for hazard analysis.

Lee, Yong Joong↗

Report on FY 2020 creep, fatigue and creep fatigue testing of Alloy 709 base metal at ORNL

The testing activities and research in support of ASME Code qualification of Alloy 709, an advanced austenitic steel, are being carried out at Oak Ridge National Laboratory (ORNL), Argonne National Laboratory, and Idaho National Laboratory. This report summarizes the status and results of FY 2020 planned testing at ORNL. Uniaxial tensile tests on the electroslag remelt (ESR) solution-annealed plate (heat number 58776-3RBC) with additional heat treatment were performed for the baseline mechanical properties evaluation for ASME code qualification of Alloy 709. The tensile properties were found to meet ASME SA-213 specifications and were comparable to the Nippon Steel NF709 data generated for an ASME Section I Code Case for seamless tubing. ORNL was tasked to carry out a subset of the Code Case testing for creep rupture. Creep rupture data from 18 tests were generated on a solution-annealed ESR plate with a solution annealing temperature of 1100°C. There are total of 19 intermediate and long-term creep rupture tests of Alloy 709 ongoing at ORNL. A preliminary fatigue design curve at 760°C was developed for Alloy 709, and the results show that the fatigue design curve of Alloy 709 is comparable to that of Alloy 800H at 760°C.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Case Study: Leveraging GenAI to Build AI-based Surrogates and Regressors for Modeling Radio Frequency Heating in Fusion Energy Science

This work presents a detailed case study on using Generative AI (GenAI) to develop AI surrogates for simulation models in fusion energy research. The scope includes the methodology, implementation, and results of using GenAI to assist in model development and optimization, comparing these results with previous manually developed models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

CSAPR2 CMAC 2.0 level c1 data.

Raw data from ARM precipitation radars must be corrected for atmospheric phenomena and instrument characteristics (e.g., attenuation, clutter) to retrieve precipitation properties. The Corrected Moments in Antenna Coordinates Version 2 (CMAC2) value-added product (VAP) is a set of algorithms and code that does such corrections, and it also retrieves precipitation quantities from the radar measurements. Similar to the X-SAPR radars at the ARM main site, C-SAPR2 radar data also needs corrections and improvements. CMAC 2.0 has been updated to work with the ARM C-SAPR 2 radar. Data and fields that have been processed to: correct for velocity aliasing, unfold and generate a cross-polarimetric phase difference that is monotonically increasing, removing impulses caused by non-uniform beam filling and phase shift on backscatter, recalculate specific differential phase using a 20-point Sobel filter on the aforementioned phase, correct for liquid path attenuation using the polarimetric signals, estimate rainfall rates at the gate using the specific attenuation. In addition, CMAC writes the data out into a community-standard format netCDF File using the CF/Radial conventions. The data are therefore compatible with new and existing National Center for Atmospheric Research (NCAR) tools such as RadXConvert for converting to a variety of popular file formats.

54 ENVIRONMENTAL SCIENCES↗

Code for multi-shape Gaussian process (GP) fitting with uncertainty quantification (UQ)

This is code associated with the publication “Nonparametric Multi-shape Modeling with Uncertainty Quantification,” authored by Hengrui Luo (Lawrence Berkeley National Laboratory) and Justin Strait (Los Alamos National Laboratory). The code is used to fit multiple-output Gaussian process (GP) models of planar closed curves to collections of ordered point sets, allowing for flexible nonlinear prediction of the underlying curve under dense or sparse point set samplings and with or without noise, as well as tractable uncertainty quantification. To do this, we employ use of a periodic kernel to account for the nonlinear input space of closed curves, and combine with coregionalization models to account for dependence both (i) between curve coordinates, and (ii) between pairs of curves. Functions in the code are capable of fitting these models, as well as performing additional tasks with the fitted curves such as (a) shape registration and alignment, (b) shape averaging, and (c) fitting for curve sub-populations / clusters.

Strait, Justin↗

RELAP5-3D Pre-test Prediction for High Temperature Test Facility Test PG-26

The RELAP5-3D computer code has been used to perform a pre-test prediction for Test PG-26 in the High Temperature Test Facility at Oregon State University. Test PG-26 will simulate a pressurized conduction cooldown transient following a trip of the primary coolant system gas circulator. The facility is initialized with a steady-state power of 700 kW and prototypic core inlet and outlet helium coolant temperatures of 259 and 687°C, respectively. Sensitivity calculations addressed the flow coastdown time, the reactor cavity cooling system flow rate and temperature, and decay power.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Review of the Evaluation of Building Energy Code Compliance in the United States

Building energy codes are essential tools for achieving energy efficiency in buildings. However, the full energy savings potential of these codes can only be realized if buildings are constructed in compliance with them. Therefore, evaluating building energy code compliance is crucial in bridging the gap between the energy efficiency requirements set by energy codes and the actualized energy savings achieved. An energy code compliance evaluation serves as a mechanism to assess construction practices, evaluate the effectiveness of code enforcement, identify gaps in compliance, and guide strategies for improvement through training and education. Conducting code compliance evaluation activities involves field studies that require careful design and significant resources. Historically, more emphasis has been placed on developing and adopting building energy codes, while efforts to evaluate compliance have been relatively limited and lacking consistent approaches. The passage of the 2009 American Recovery and Reinvestment Act (ARRA), which mandated that states create plans for achieving 90% compliance within eight years, stimulated the need for an energy code compliance evaluation. As a result, federal, state, and local governments, and utilities have invested in the development of methodologies and tools for code compliance evaluation studies. This paper reviews the code compliance evaluation studies conducted in the United States over the past three decades. It describes and compares the methodologies and metrics used to assess building energy code compliance, summarizes the general elements and steps involved in the evaluation process, and discusses common issues in these studies. Over time, code compliance evaluation methodologies have evolved from isolated development within individual states, regions, and utilities, to widely accepted protocols applicable across different states and local jurisdictions. There has been a transition in compliance metrics, shifting from historical compliance rates to energy-consumption-oriented approaches.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Advanced Ultra-Supercritical Component Test (ComTest) Project for 760*C Steam Conditions

ComTest is a $27M Department of Energy-funded project. Phase I, which began in November 2015, served to identify the technology gaps, as well as the scope and cost of required testing. Phase II, which was awarded in December 2018, includes an advanced manufacturing effort to complete U.S. based supply chain development for full commercial scale (800-850 MWe) AUSC components made of nickel-based alloys, components operating at up to 760°C. Completion of Phase II scheduled for September 30, 2021. Strategic Objectives of U.S. AUSC Project include: Power Plant Efficiency Improvements –Develop cost-effective, reliable technologies to improve the efficiency of new and existing high-temperature advanced generation power plants. Close gaps to achieve readiness for commercial scale demonstration of Advanced Ultra-Supercritical (AUSC) technology. Fabricate full-scale versions of key nickel-based alloy components. Validate capabilities of US supply chain for cost certainty. Support cross-cutting high-temperature generation technologies. Obtain ASME Code approval for new materials, components and processes. Increase power plant steam temperatures for higher cycle efficiency (Note: Average efficiency of US coal-fired fleet = 33% HHV; A-USC plant efficiency over 47% HHV at 1,400°F (760°C) steam temperature).

fireside corrosion↗

Data Analysis of Energy Code Compliance in Commercial Buildings (Rev. 1)

What percent of newly constructed commercial buildings comply with the energy code? How much energy and cost could be saved if the compliance increased? Which code requirements have both low compliance rates and high savings potential? These are the questions U.S. Department of Energy (DOE) is trying to answer through its Commercial Energy Code Field Study. Previous commercial studies have been very limited and did not result in a widely accepted and tested methodology. DOE’s goal is to create a standardized methodology that can produce actionable results at a reasonable study cost, that can be used by state and local governments and utilities and provides valuable information to policy makers. The field study team implemented the pilot methodology, compiling a data set of 230 office and retail buildings in two climate zones. The approach is based on identifying lost savings on a total energy cost basis rather than simply counting the quantity of measures that do not meet code. This report is focused on the analysis of the collected data. The primary goal was to analyze the data collected during the field study and determine the actual energy cost impact of each measure in a non-compliance situation. The energy impact results allowed for ranking the measures to identify which have the highest potential for lost savings. These results combined with the time required to verify each measure will allow future compliance verification to focus on measures that had a large impact on energy use over the life of the building and those that have the greatest savings recovery potential per verification hour.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

HiPACE++: A portable, 3D quasi-static particle-in-cell code

Modeling plasma accelerators is a computationally challenging task and the quasi-static particle-in-cell algorithm is a method of choice in a wide range of situations. In this work, we present the first performance-portable, quasi-static, three-dimensional particle-in-cell code HiPACE++. By decomposing all the computation of a 3D domain in successive 2D transverse operations and choosing appropriate memory management, HiPACE++ demonstrates orders-of-magnitude speedups on modern scientific GPUs over CPU-only implementations. The 2D transverse operations are performed on a single GPU, avoiding time-consuming communications. The longitudinal parallelization is done through temporal domain decomposition, enabling near-optimal strong scaling from 1 to 512 GPUs. HiPACE++ is a modular, open-source code enabling efficient modeling of plasma accelerators from laptops to state-of-the-art supercomputers.

97 MATHEMATICS AND COMPUTING↗

Effects of electron cyclotron current drive on magnetic islands in tokamak plasmas

The effects of the electron cyclotron current drive on magnetic islands in tokamak plasmas are studied using gyrokinetic simulations. By investigating the effects of different characteristics of the driven current, such as current density distribution and deposition location, the factors which can determine the suppression effect on the resistive tearing modes have been explored. It is found that an electron cyclotron wave (ECW) driven current with a larger peak value and more focused deposition region has a better stabilization effect. When the ECW-driven current is closer to the rational surface, it has a better stabilizing effect. These gyrokinetic toroidal code (GTC) linear simulations in the electron fluid limit of the tearing modes in the cylindrical geometry agree well with the magnetohydrodynamic codes. In conclusion, the optimal timing control of the current deposition on resistive tearing modes is demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Laser beam coupling with capillary discharge plasma for laser wakefield acceleration applications

One of the most robust approaches, demonstrated to date, of accelerating electron beams by laser-plasma sources is the utilization of plasma channels generated by the capillary discharges. Although the spatial structure of the installation is simple in principle, there may be some important effects caused by the open ends of the capillary, by the supplying channels etc., which require a detailed 3D modeling of the processes. In the present work, such simulations are performed using the code MARPLE. First, the process of capillary filling with cold hydrogen before the discharge is fired, through the side supply channels is simulated. Second, the simulation of the capillary discharge is performed with the goal to obtain a time-dependent spatial distribution of the electron density near the open ends of the capillary as well as inside the capillary. Lastly, to evaluate the effectiveness of the beam coupling with the channeling plasma wave guide and of the electron acceleration, modeling of the laser-plasma interaction was performed with the code INF&RNO.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Characterization of an imploding cylindrical plasma for electron transport studies using x-ray emission spectroscopy

We report here on the characterization of the conditions of an imploding cylindrical plasma by time-resolved x-ray emission spectroscopy. Knowledge about this implosion platform can be applied to studies of particle transport for inertial confinement fusion schemes or to astrophysical plasmas. A cylindrical Cl-doped CH foam within a tube of solid CH was irradiated by 36 beams (I total ~5×10 14 W/cm 2 , 1.5ns square pulse, and E total ~ 16.2kJ) of the OMEGA-60 laser to radially compress the CH toward the axis. The analysis of the time-resolved spectra showed that the compression can be described by four distinct phases, each presenting different plasma conditions. First the ablation of the cylinder is dominant; second, the foam is heated and induces a significant jump in emission intensities; third, the temperature and density of the foam reaches a maximum; and finally, the plasma expands. Ranges for the plasma temperature were inferred with the atomic physics code SCRAM (Spectroscopic Collisional-Radiative Atomic Model) and the experimental data have been compared to hydrodynamic simulations performed with the 2D code FLASH, which showed a similar implosion dynamic over time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Particle-in-cell modeling of plasma jet merging in the large-Hall-parameter regime

The merging process of magnetized plasma jets with parameters relevant to the plasma-jet-driven magneto-inertial fusion (PJMIF) design and the plasma liner experiment (PLX) is modeled by fully kinetic particle-in-cell (PIC) simulations in one and two spatial dimensions. Here, the modified two-stream instability is identified to be the main mechanism responsible for stopping the plasma jets and preventing species interpenetration. The electron and ion Hall parameters of the merged plasma are greater than unity, and the plasma β is close to unity, which is the desired characteristic of planned experiments at PLX. Our 2D PIC simulations validate the results of the radiation magneto-hydrodynamics code FLASH, which will be the primary tool for modeling various stages of future PJMIF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlinear light-output calibration of the oxygenated xylene scintillators used in OMEGA neutron time-of-flight spectrometers

Neutron time-of-flight (nTOF) spectrometers are essential instruments for measuring and evaluating the performance of inertial confinement fusion implosions. The neutron spectrometers utilized for the OMEGA laser include two liquid-based scintillators, each consisting of a large volume filled with xylene that is coupled to four photomultiplier tubes. Analysis of the signal from these detectors requires detailed knowledge of the scintillator’s light output, which is needed to fit the nTOF spectrum, from which the neutron energy spectrum is informed. The light output is nonlinearly proportional to the neutron energy, which, in turn, affects the interpretation of the neutron energy spectrum from a TOF signal. A recent campaign on OMEGA was performed to calibrate the xylene detectors and infer the shape of the light-output curve. The campaign utilized materials with increasing Z placed in the OMEGA target chamber to initiate scattering events with the 14 MeV fusion neutrons. This process leads to the production of backscatter neutrons of varying energies that appear as peaks in the nTOF data. Simulations using a neutron transport code were combined with the measured deuterium–tritium neutron yields to calculate the expected backscattered neutron yields from the well-known scattering cross sections of each material. Furthermore, the neutron-energy dependent light output of the scintillator inferred from the experiment is compared to the light-output curve simulated with a neutron transport code for the following neutron energies: 1.5, 2.5, 6, and 14 MeV.

47 OTHER INSTRUMENTATION↗

A Code-Agnostic Driver Application for Coupled Neutronics and Thermal-Hydraulic Simulations

While the literature has numerous examples of Monte Carlo and computational fluid dynamics (CFD) coupling, most are hard-wired codes intended primarily for research rather than as standalone, general-purpose applications. In this work, we describe an open source application, ENRICO, that enables coupled neutronic and thermal-hydraulic simulations between multiple codes that can be chosen at runtime (as opposed to a coupling between two specific codes). The application has been designed such that the control flow logic, domain mapping, nonlinear fixed-point iteration, solution transfers, and convergence checks are all agnostic to the underlying physics solvers used. Special emphasis has also been placed on enabling efficient execution on distributed-memory computing environments. The transfer of solution fields between solvers is performed in memory rather than through filesystem I/O. Additionally, solvers can be configured to run on overlapping or disjoint sets of processes. To date, coupling with the OpenMC and Shift Monte Carlo codes, the Nek5000 CFD code, and a simplified heat diffusion and subchannel solver has been implemented in ENRICO. We present results for coupled simulations of a single light-water reactor fuel assembly based on the NuScale reactor using various combinations of the physics solvers. For this problem, the coupled simulations are shown to converge in about four Picard iterations. A comparison of the heat source and temperature distributions computed by ENRICO using OpenMC coupled with Nek5000 and Shift coupled with Nek5000 illustrates remarkable agreement between the codes.

42 ENGINEERING↗

Computation of generalised magnetic coordinates asymptotically close to the separatrix

Integrals to calculate generalised magnetic coordinates from an input magnetic flux function asymptotically close to the separatrix are presented, and implemented in the GPEC/DCON code suite. These integrals allow characterisation of the magnetic equilibrium of a diverted tokamak, in magnetic coordinates, arbitrarily close to the last closed flux surface, avoiding the numerical issues associated with calculating diverging field-line integrals near a magnetic x-point. Finally, these methods may assist ongoing efforts to develop robust asymptotic equilibrium behaviour for spectral 3D MHD codes at the separatrix.

equilibrium edge truncation↗

Localizing resonant magnetic perturbations for edge localized mode control in KSTAR

An external 3D magnetic perturbation typically drives a resonant response at the rational surfaces from the core to the edge of tokamak plasmas, due to strong mode coupling and amplification. In this paper, we present a method to isolate the edge from core resonant fields using the ideal perturbed equilibrium code and to design an edge-localized resonant magnetic perturbation (RMP) for effective edge localized mode (ELM) control. A robust feature of the edge-localized RMP is the curtailed response to the field at the low-field-side (LFS) midplane, as opposed to typical RMPs which strongly resonate with the LFS fields. This emphasizes the importance of off-midplane coils to improve ELM control without provoking a large core response that could lead to devastating instabilities. The conceptual design of new ELM control coils based on the edge-localized RMP in KSTAR shows how this new insight can be utilized to enhance the efficiency of our ELM suppression capabilities. Simple window-pane coils matching the edge-localized resonant mode structure substantially expand in the ELM suppression window beyond the existing coil. Further optimization using the flexible optimized coils using space-curves code leads to additional enhancement in the edge-localized control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗