SEARCH · Search NASA

Results for “SPACE 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Systematic approach to define a prototype building model: Convenience Store

The US Department of Energy supports the development of commercial building energy codes and standards. To support commercial building energy research activities, continuous efforts have been made to convert 30 prototype building models, which cover 80% of commercial building floor space in the United States, into OpenStudio Standards prototype buildings. In this paper, we provide a systematic approach to define the building and system characteristics of a convenience store, which has the largest energy use intensity among the commercial building types not covered by current prototype models. This study will be used as a basis for developing the prototype convenience store building model.

Yoon, Yeobeom↗

DOE OSTI2021

RAVEN User Manual

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

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Epistatic Net allows the sparse spectral regularization of deep neural networks for inferring fitness functions

Abstract Despite recent advances in high-throughput combinatorial mutagenesis assays, the number of labeled sequences available to predict molecular functions has remained small for the vastness of the sequence space combined with the ruggedness of many fitness functions. While deep neural networks (DNNs) can capture high-order epistatic interactions among the mutational sites, they tend to overfit to the small number of labeled sequences available for training. Here, we developed Epistatic Net (EN), a method for spectral regularization of DNNs that exploits evidence that epistatic interactions in many fitness functions are sparse. We built a scalable extension of EN, usable for larger sequences, which enables spectral regularization using fast sparse recovery algorithms informed by coding theory. Results on several biological landscapes show that EN consistently improves the prediction accuracy of DNNs and enables them to outperform competing models which assume other priors. EN estimates the higher-order epistatic interactions of DNNs trained on massive sequence spaces-a computational problem that otherwise takes years to solve.

Aghazadeh, Amirali (ORCID:0000000302230873)↗

DOE OSTI2022

Carbon and energy cost impacts of electrification of space heating with heat pumps in the US

In order to meet climate goals, it will be necessary to significantly reduce the greenhouse gases emitted by homes. A key factor in the US is to reduce the on-site combustion of fossil fuels for heating end-uses and to replace this with use of electric heat pump technologies connected to a low-carbon grid. The replacement of natural gas furnaces with electric heat pumps is a key home decarbonization strategy. However, the potential for space heating electrification to reduce greenhouse gas emissions depends on the carbon dioxide equivalent (CO 2 e) content of the electricity used by the heat pump. This varies considerably depending on the source of electricity, with large state to state variability. Furthermore, household energy costs are likely to be impacted by the electrification of space heating, because retail energy prices for both natural gas and electricity in each state vary by factors of seven and four, respectively. Contractors, energy programs, government and building code officials, as well as consumers need clear indications of the likely CO 2 e and energy cost impacts of proposed electrification projects, because these will affect decarbonization choices and rationales around scaled heating electrification. Government and utility programs also need to be aware of the likely outcomes of any supported/incentivized measures. Here, in this paper, we investigate these effects by looking at new metrics to analyze the change in CO 2 e emitted and the cost to meet home heating loads when switching from a natural gas furnace to a heat pump for the contiguous 48 states of the mainland US.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2023

Energy transfer and scale dynamics in 2D and 3D laser-driven jets

We demonstrate a methodology for diagnosing the multiscale dynamics and energy transfer in complex HED flows with realistic driving and boundary conditions. The approach separates incompressible, compressible, and baropycnal contributions to energy scale-transfer and quantifies the direction of these transfers in (generalized) wavenumber space. We use this to compare the kinetic energy (KE) transfer across scales in simulations of 2D axisymmetric vs fully 3D laser-driven plasma jets. Using the FLASH code, we model a turbulent jet ablated from an aluminum cone target in the configuration outlined. We show that, in addition to its well known bias for underestimating hydrodynamic instability growth, 2D modeling suffers from significant spurious energization of the bulk flow by a turbulent upscale cascade. In 2D, this arises as vorticity and strain from instabilities near the jet's leading edge transfer KE upscale, sustaining a coherent circulation that helps propel the axisymmetric jet farther (⁠ ≈25% by 3.5 ns) and helps keep it collimated. In 3D, the coherent circulation and upscale KE transfer are absent. Here, the methodology presented here may also help with inter-model comparison and validation, including future modeling efforts to alleviate some of the 2D hydrodynamic artifacts highlighted in this study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Optimization of a High Bunch Charge ERL Injection Merger for PERLE

Delivery of high charge electron bunches into the main loop of an ERL (energy recovery linac) while preserving the emittance is challenging. This is because at the typical injection momentum, space charge forces still have a significant effect on the beam dynamics. In this work we consider the design of the merger for PERLE, an ERL test facility to be based at IJCLab in France. Previous simulations have shown that the baseline DC gun based injector can achieve the required emittance at the booster linac exit. The quality of the 500 pC bunches must then be preserved with space charge through the merger at total beam energy of 7 MeV keeping the emittance below 6 mm mrad. The beam dynamics in the merger were simulated using the code OPAL and optimised using a genetic algorithm. Three possible merger schemes were investigated. The goal of the optimisation was to minimise the emittance growth while also achieving the required Twiss parameters to match onto the spreader at the main linac exit. A three dipole solution is then examined in more detail.

Hounsell, B.↗

DOE OSTI2021

Design of a 10 MeV beamline at the Upgraded Injector Test Facility for e-beam irradiation

Electron beam irradiation near 10 MeV is suitable for wastewater treatment. The Upgraded Injector Test Facility (UITF) at Jefferson Lab is a CW superconducting linear accelerator capable of providing an electron beam of energy up to 10 MeV and up to 100 µA current. This contribution presents the beam transport simulations for a beamline to be used for the irradiation of wastewater samples at the UITF. The simulations were done using the code General Particle Tracer with the goal of obtaining an 8 MeV electron beam of radius (3-σ) of ~2.4 cm. The achieved energy spread is ~74.5 keV. The space charge effects were investigated when the bunch charge is varied to be up to 1000 times and the results showed that they do not affect the beam quality significantly.

Li, X.↗

DOE OSTI2021

Design of a 10 MeV beamline at the Upgraded Injector Test Facility for e-beam irradiation

Li, X.↗

DOE OSTI2020

Reduced Models for Chemical Kinetics derived from Parallel Ensemble Simulations of Stirred Reactors

Detailed chemical kinetic models are essential to improve the predictive capability of computational fluid dynamics (CFD) simulations and to impact the development of next generation combustion technologies. However, such detailed kinetic models are beyond the reach of CFD codes due to the computational cost of transporting additional chemical degrees of freedom. Reduced order models for combustion kinetics exploit the presence of lower dimensional manifolds in the state space and require far fewer number of scalars to be transported as part of CFD simulations. In this paper, we describe the software infrastructure developed to allow large ensembles of stirred reactor configurations with detailed chemical kinetics. The canonical reactor simulations are incorporated in a workflow with sparse grid and data reduction methodologies aimed at development of reduced order models. We present results from analyzing the results of the stirred reactor calculations for syngas mixture using principal component analysis (PCA) and neural network methodologies. The results show that nonlinear autoencoder model can surpass PCA in representing the original dataset with far fewer degrees of freedom.

Zhang, Pei↗

DOE OSTI2023

Improving the Estimation of the Atmospheric Water Vapor Pressure Using Interpretable Long Short-Term Memory Networks: Dataset, Python code, and trained models

Atmospheric water vapor pressure is an essential meteorological control on land surface and hydrologic processes. It is not as frequently observed as other meteorologic conditions, but often inferred through the August–Roche–Magnus formula by simply assuming dew point and daily minimum temperatures are equivalent or by empirically correlating the two temperatures using an aridity correction. The performance of both methods varies considerably across different regions and during different time periods; obtaining consistently accurate estimates across space and time remains a great challenge. We applied an interpretable Long Short-Term Memory (iLSTM) network conditioned on static, location specific attributes to estimate daily vapor pressure for 83 FLUXNET sites in the United States and Canada. This data package includes all raw data of the 83 FLUXNET sites, input data for model training/validation/test, trained models and results, and python codes for the manuscript "Improving the Estimation of the Atmospheric Water Vapor Pressure Using an Interpretable Long Short-term Memory Network". Specifically, it consists of five parts. - First, "1_Daymet_data_83sites.zip" includes raw data downloaded from Daymet for the 83 sites used in the paper according to their longitude and latitude, in which vapor pressure is used. It also includes a pre-processed CSV data file combining all data from the 83 sites which is specifically used for the paper. - Second, "2_Fluxnet2015_data_83sites.zip" includes raw half hourly data of the 83 sites downloaded from FLUXNET2015 data portal, pre-processed daily data of the 83 sites, a CSV file including combined pre-processed daily data of the 83 sites, and a CSV file including the information (site ID, site name, latitude, longitude, data available period) of the 83 sites. - Third, "3_MODIS_LAI_data_83sites_raw.zip" includes raw leaf area index (LAI) data downloaded from the AppEEARs data portal. - Fourth, "4_Scripts.zip" includes all scripts related to model training and post-processing of a trained model, and a jupyter notebook showing an example for model post-processing. Two typo errors in files titled "run2get_args.py" and "postprocess.py" were corrected on March 27, 2024 to avoid confusions. - Finally, "Trained_models_and_results.zip" includes three folders and three files with suffix ".npy", and each folder corresponds to one file with suffix ".npy" with the same title. Each of the three folders include all trained models associated with one iLSTM model configuration (35 models for each configuration, details are described in the paper). Each file with suffix ".npy" includes the post-processed results of the corresponding 35 models under one iLSTM model configuration.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Airfoil Computational Fluid Dynamics - 9k shapes, 2 AoA's

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 8,996 airfoil shapes, computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this data, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We fixed the linear deformations to be the mean over the database and sampled new shapes over a four-dimensional parameter space of higher-order perturbation. This sampling approaches allows for isolated analysis of non-linear airfoil shape deformations while holding other aspects (e.g., airfoil thickness) approximately constant. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, Reynolds number of 9M, and at two angles of attack, 4 deg. and 12 deg. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_9k_data.ipynb.

9k↗

DOE OSTI2025

Propagation of ion cyclotron emission in the DIII-D tokamak

Experiments in L- and H-mode plasmas performed on the DIII-D tokamak explored ion cyclotron emission (ICE) propagation via the recently upgraded ICE diagnostic. The distance between the plasma and the outer wall was scanned to alter the evanescent region, which was seen to impact the detection of runaway electron-driven whistler modes in a comparable frequency range to typical ICE harmonics (f≳10fci) [Heidbrink et al., Plasma Phys. Controlled Fusion 61, 014007 (2018)]. In L-mode plasmas, core ICE did not exhibit clear trends as the plasma-wall distance was altered. Instead, inadvertent changes to the fast ion density correlated with different ICE fine structures; the most spectral bands (spaced ∼150–200 kHz apart) were observed at the highest fast ion density, and a just single band when nfast decreased by ∼ 40%. The full-wave Petra-M code simulated core ICE propagation through the plasma and to ICE probes, agreeing with experimental findings that these are likely compressional fast waves, which propagate to probe locations with no dependence on plasma-wall gap. Edge ICE in H-mode plasmas was similarly unaffected by plasma-wall distance. Unlike in core ICE cases, there were negligible changes to the edge fast ion distribution and edge bulk plasma profiles. However, changes in ELM frequency seemed most closely tied to ICE behavior. Ultimately, core and edge ICE did not exhibit clear dependencies on the plasma-wall distance, and ICE harmonics were observed with even the largest plasma-wall gaps. This independence bodes well for similar diagnostics in future tokamaks, which might be placed further from high radiation areas without adversely impacting measurement capabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Symmetry-protected self-correcting quantum memory in three space dimensions

Whether self-correcting quantum memories can exist at nonzero temperature in a physically reasonable setting remains a great open problem. Furthermore, it has recently been argued that symmetry-protected topological (SPT) systems in three space dimensions subject to a strong constraint—that the quantum dynamics respect a 1-form symmetry—realize such a quantum memory. We illustrate how this works in Walker-Wang codes, which provide a specific realization of these desiderata. In this setting we show that it is sufficient for the 1-form symmetry to be enforced on a subvolume of the system. This strongly suggests that the SPT character of the state is not essential. We confirm this by constructing an explicit example with a trivial (paramagnetic) bulk that realizes a self-correcting quantum memory. We therefore show that the enforcement of a 1-form symmetry on a measure-zero subvolume of a three-dimensional system can be sufficient to stabilize a self-correcting quantum memory at nonzero temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Polarization–frequency hyperentangled photons: generation, characterization, and manipulation

Frequency-bin encoding is massively parallelizable and robust for optical fiber transmission. When coupled with an additional degree of freedom (DoF), the expansion of the Hilbert space allows for deterministic controlled operations between two DoFs within a single photon. Such capabilities, when combined with photonic hyperentanglement, are of great value for quantum communication protocols, including dense coding and single-copy entanglement distillation. In this talk, we present an all-fiber-coupled, ultrabroadband polarization–frequency hyperentangled source and conduct comprehensive quantum state tomography across multiple dense wavelength division multiplexing channels spanning the optical C+L-band (1530–1625 nm). In addition, we design and implement a high-fidelity controlled-NOT (cnot) operation between polarization and frequency DoFs by exploiting electro-optic phase modulation within a fiber Sagnac loop. Collectively, our hyperentangled source and two-qubit gate should unlock new opportunities for harnessing polarization–frequency resources in established telecommunication fiber networks for future quantum applications.

Lu, Hsuan-Hao↗

DOE OSTI2023

Pressurized Water Reactor Gadolinia Pin Location Optimization

This report presents the results of lattice optimization studies performed to find optimum locations for gadolinia burnable absorber (BA) rods in pressurized water reactor (PWR) lattice fuel designs. Initial excess reactivity suppression allows core designers to further improve operational economics by extending cycle length. Gadolinia BAs are commonly used in boiling water reactor assembly designs for this purpose. In recent years, gadolinia absorbers have been used in PWR designs owing to their longer effectiveness for reactivity suppression compared with common BAs used in PWR assemblies. This report examines the optimum gadolinia pin placement in 17 × 17 PWR lattices at different fuel and gadolinia concentrations for optimized lattice performance, using the SCALE/Polaris lattice physics code. The completed work is continuation of the Light Water Reactor LEU+ Lattice Optimization (ORNL/TM-2021/2366) project. An optimization driver called the metaheuristic optimization tool (MOT) is used to automate domain space exploration and optimization of the lattice designs. Heuristics from previous light-water reactor (LWR) lattice optimization studies were used to construct the objective function and define the domain space for optimization. This work successfully demonstrated that the optimization algorithms of MOT can generate feasible, nonproprietary PWR lattice designs with gadolinia.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

A Full-Stack Exploration of Language-Based Parallelism in Fortran 2023

This poster explores native parallel features in Fortran 2023 through the lens of supporting applications with libraries, compilers, and parallel runtimes. The language revision informally named Fortran 2008 introduced parallelism in the form of Single Program Multiple Data (SPMD) execution with two broad feature sets: (1) loop-level parallelism via do concurrent and (2) a Partitioned Global Address Space (PGAS) comprised of distributed “coarray” data structures. Fortran’s native parallelism has demonstrated high performance [1] and reduced the burden of inserting what sometimes amounts to more directives than code. Several compilers support both feature sets, typically by translating do concurrent into serial do loops annotated by parallel directives and by translating SPMD/PGAS features into direct calls to a communication library. Our research focuses primarily on two questions: (1) can the compiler’s parallel runtime library be developed in the language being compiled (Fortran) and (2) can we define an interface to the runtime that liberates compilers from being hardwired to one runtime and vice versa. We are answering these questions by developing the Parallel Runtime Interface for Fortran (PRIF) [2] and the Co-Array Fortran Framework of Efficient Interfaces to Network Environments (Caffeine) [3]. Caffeine is initially targeting adoption by LLVM Flang, a new open-source Fortran compiler developed by a broad community in industry, academia, and government labs. We are also exploring the use of these features in Inference-Engine, a deep learning library designed to facilitate neural network training and inference for high-performance computing applications written in modern Fortran.

Rasmussen, Katherine↗

DOE OSTI2021

MOA-2009-BLG-319Lb: A Sub-Saturn Planet inside the Predicted Mass Desert

We present an adaptive optics (AO) analysis of images from the Keck II telescope NIRC2 instrument of the planetary microlensing event MOA-2009-BLG-319. The ∼10 yr baseline between the event and the Keck observations allows the planetary host star to be detected at a separation of 66.5 ± 1.7 mas from the source star, consistent with the light-curve model prediction. The combination of the host star brightness and light-curve parameters yields host star and planet masses of M {sub host} = 0.524 ± 0.048M {sub ⊙} and m {sub p} = 67.3 ± 6.2M {sub ⊕} at a distance of D {sub L} = 7.1 ± 0.7 kpc. The star−planet projected separation is 2.03 ± 0.21 au. The planet-to-star mass ratio of this system, q = (3.857 ± 0.029) × 10{sup −4}, places it in the predicted “planet desert” at 10{sup −4} < q < 4 × 10{sup −4} according to the runaway gas accretion scenario of the core accretion theory. Seven of the 30 planets in the Suzuki et al. sample fall in this mass ratio range, and this is the third with a measured host mass. All three of these host stars have masses of 0.5 ≤ M {sub host}/M {sub ⊙} ≤ 0.7, which implies that this predicted mass ratio gap is filled with planets that have host stars within a factor of two of 1M {sub ⊙}. This suggests that runaway gas accretion does not play a major role in determining giant planet masses for stars somewhat less massive than the Sun. Our analysis has been accomplished with a modified DAOPHOT code that has been designed to measure the brightness and positions of closely blended stars. This will aid in the development of the primary method that the Nancy Grace Roman Space Telescope mission will use to determine the masses of microlens planets and their hosts.

47 OTHER INSTRUMENTATION↗

DOE OSTI2020

GW190814's Secondary Component with Mass 2.50–2.67 M ⊙ as a Superfast Pulsar

Using Stergioulas's RNS code for investigating fast pulsars with Equation of States (EOSs) on the causality surface (where the speed of sound equals that of light) of the high-density EOS parameter space satisfying all known constraints from both nuclear physics and astrophysics, we show that one possible explanation for the GW190814's secondary component of mass (2.50-2.67) M ⊙ is that it is a super-fast pulsar spinning faster than 971 Hz about 42% below its Kepler frequency. If confirmed, it would be the fastest pulsar with the highest mass observed presently. There is a large and physically allowed EOS parameter space below the causality surface where pulsars heavier than 2.50 M ⊙ are supported if they can rotate even faster with critical frequencies depending strongly on the high-density behavior of nuclear symmetry energy.

79 ASTRONOMY AND ASTROPHYSICS↗