Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “source 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,261 records · Page 70

A Procedure for Evaluating Modified Pulsed-Neutron- Source Experiments in Subcritical Nuclear Reactors

A theoretical interpretation of the modified pulsed-neutron-source experiments of Sjöstrand, Gozani, and Garelis and Russell is given using exact steady-state Boltzmann equations. Here, interpretation is based on a phenomenological description of the experiments and is patterned after work done on the Garelis-Russell method by Corngold. The basic approximation made is that the fundamental prompt-mode decay constant is much larger than any delayed-neutron precursor decay constant. The theoretical interpretation allows the reactivities measured by the above three modified pulsed-source techniques to be easily calculated and compared to more conventional definitions of reactivity. The calculations can be performed by any standard source-iteration code that has been modified to solve the inhomogeneous problem. Experiments were performed on the Cornell University Critical Assembly and interpreted with the aid of the above theory. Calculations and experiments agree to within 20%. Sjöstrand's method is found to give the best result for this reactor.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Rotational Raman-Based Temperature Measurements in a High-Velocity Turbulent Jet

Spontaneous rotational Raman scattering spectroscopy is used to acquire the first ever high quality, spatially-resolved measurements of the mean and root mean square (rms) temperature fluctuations in turbulent, high-velocity heated jets. Raman spectra in air were obtained across a matrix of radial and axial locations downstream from a 50 mm diameter nozzle operating from subsonic to supersonic conditions over a wide range of temperatures and Mach numbers, in accordance with the Tanna matrix frequently used in jet noise studies. These data were acquired in the hostile, high noise (115 dB) environment of a large scale open air test facility at NASA Glenn Research Center (GRC). Temperature estimates were determined by performing nonlinear least squares fitting of the single shot spectra to the theoretical rotational Stokes spectra of N2 and O2, using a custom in-house code developed specifically for this investigation. The laser employed in this study was a high energy, long-pulsed, frequency doubled Nd:YAG laser. One thousand single-shot spectra were acquired at each spatial coordinate. Mean temperature and rms temperature variations were calculated at each measurement location. Excellent agreement between the averaged and single-shot temperatures was observed with an accuracy better than 2.5 percent for temperature, and rms variations in temperature between +/-2.2 percent at 296 K and +/-4.5 percent at 850 K. The results of this and planned follow-on studies will support NASA GRC's development of physics-based jet noise prediction, turbulence modeling and aeroacoustic source modeling codes.

Randy J Locke↗

Rotational Raman-Based Temperature Measurements in a High-Velocity Turbulent Jet

Spontaneous rotational Raman scattering spectroscopy is used to acquire the first ever high quality, spatially-resolved measurements of the mean and root mean square (rms) temperature fluctuations in turbulent, high-velocity heated jets. Raman spectra in air were obtained across a matrix of radial and axial locations downstream from a 50 mm diameter nozzle operating from subsonic to supersonic conditions over a wide range of temperatures and Mach numbers, in accordance with the Tanna matrix frequently used in jet noise studies. These data were acquired in the hostile, high noise (115 dB) environment of a large scale open air test facility at NASA Glenn Research Center (GRC). Temperature estimates were determined by performing nonlinear least squares fitting of the single shot spectra to the theoretical rotational Stokes spectra of N2 and O2, using a custom in-house code developed specifically for this investigation. The laser employed in this study was a high energy, long-pulsed, frequency doubled Nd:YAG laser. One thousand single-shot spectra were acquired at each spatial coordinate. Mean temperature and rms temperature variations were calculated at each measurement location. Excellent agreement between the averaged and single-shot temperatures was observed with an accuracy better than 2.5 percent for temperature, and rms variations in temperature between +/-2.2 percent at 296 K and +/-4.5 percent at 850 K. The results of this and planned follow-on studies will support NASA GRC's development of physics-based jet noise prediction, turbulence modeling and aeroacoustic source modeling codes.

Randy J Locke↗

KiteFAST

KiteFAST is a simulator for airborne wind energy systems based on NREL's OpenFAST whole turbine simulator. KiteFAST is a multi-physics engineering tool for modeling the coupled aero-servo-elastic dynamics of airborne wind energy systems i.e. a wind-energy kite. The engineering model is like NREL's OpenFAST tool for land-based or offshore horizontal-axis wind turbines atop towers, in that KiteFAST models the important physical phenomena and system couplings, including the wind/aerodynamic excitation and full-system dynamics response (fuselage, wings, tail, stabilizers, pylons, nacelles, rotors, tether, ground station, control surfaces, and controller) under both normal (for fatigue) and extreme (for ultimate) loading conditions. The configuration of the wind-energy kite is assumed to have one fuselage; two wings (starboard, port) attached to the fuselage; one vertical stabilizer attached the fuselage; two horizontal stabilizers (starboard, port) attached to the vertical stabilizer; one or more pylons per wing, each with two rotors (top, bottom); one ground-station-fixed tether with a bridal connection to the wings, and nominal circular motion of the energy kite in the port (left) direction. The general configuration of the control surfaces assumes one or more flaps per wing, two rudders on the vertical stabilizer, two elevators per horizontal stabilizer, and rotors with variable speed and rotor-collective blade pitch. KiteFAST was developed from (1) the existing MBDyn open-source multi-body dynamics code; (2) the existing/modified OpenFAST modules NWTC Subroutine Library, AirfoilInfo, InflowWind, MoorDyn, and TurbSim; and (3) new source code for a) generating MBDyn models of the energy kite (KiteMBDyn Preprocessor), b) computing quasi-steady energy-kite aerodynamics (KiteAeroDyn, KiteAeroDyn Driver, KiteVSM, ActuatorDisk), c) interfacing an energy-kite controller, and d) coupling the various KiteFAST modules together (KiteFASTMBD).

Jonkman, Jason↗

A Data Analysis Framework for Earth System Simulation within an In-Situ Infrastructure

This paper presents a generic procedure to implement a scalable and high performance data analysis framework for large-scale scientific simulation within an in-situ infrastructure. It demonstrates a unique capability for global Earth system simulations using advanced computing technologies (i.e ., automated code analysis and instrumentation), in-situ infrastructure (i.e ., ADIOS) and big data analysis engines (i.e ., SciKit-learn). This paper also includes a useful case that analyzes a globe Earth System simulations with the integration of scalable in-situ infrastructure and advanced data processing package. Finally, the in-situ data analysis framework can provides new insights on scientific discoveries in multiscale modeling paradigms.

97 MATHEMATICS AND COMPUTING↗

SALT3: An Improved Type Ia Supernova Model for Measuring Cosmic Distances

Abstract A spectral-energy distribution (SED) model for Type Ia supernovae (SNe Ia) is a critical tool for measuring precise and accurate distances across a large redshift range and constraining cosmological parameters. We present an improved model framework, SALT3, which has several advantages over current models—including the leading SALT2 model (SALT2.4). While SALT3 has a similar philosophy, it differs from SALT2 by having improved estimation of uncertainties, better separation of color and light-curve stretch, and a publicly available training code. We present the application of our training method on a cross-calibrated compilation of 1083 SNe with 1207 spectra. Our compilation is 2.5× larger than the SALT2 training sample and has greatly reduced calibration uncertainties. The resulting trained SALT3.K21 model has an extended wavelength range 2000–11,000 Å (1800 Å redder) and reduced uncertainties compared to SALT2, enabling accurate use of low- z I and iz photometric bands. Including these previously discarded bands, SALT3.K21 reduces the Hubble scatter of the low- z Foundation and CfA3 samples by 15% and 10%, respectively. To check for potential systematic uncertainties, we compare distances of low (0.01 < z < 0.2) and high (0.4 < z < 0.6) redshift SNe in the training compilation, finding an insignificant 3 ± 14 mmag shift between SALT2.4 and SALT3.K21. While the SALT3.K21 model was trained on optical data, our method can be used to build a model for rest-frame NIR samples from the Roman Space Telescope. Our open-source training code, public training data, model, and documentation are available at https://saltshaker.readthedocs.io/en/latest/ , and the model is integrated into the sncosmo and SNANA software packages.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluating wind speed and power forecasts for wind energy applications using an open-source and systematic validation framework

Building on the verification and validation work developed under the Second Wind Forecast Improvement Project, this work exhibits the value of a consistent procedure to evaluate wind power forecasts. We established an open-source Python code base tailored for wind speed and wind power forecast validation, WE-Validate. The code base can evaluate model forecasts with observations in a coherent manner. To demonstrate the systematic validation framework of WE-Validate, we designed and hosted a forecast evaluation benchmark exercise. We invited forecast providers in industry and academia to participate and submit forecasts for two case studies. We then evaluated the submissions with WE-Validate. Our findings suggest that ensemble means have reasonable skills in time series forecasting, whereas they are often inferior to single ensemble members in wind ramp forecasting. Adopting a voting scheme in ramp forecasting that allows ensemble members to detect ramps independently leads to satisfactory skill scores. Throughout this document, we also emphasize the importance of using statistically robust and resistant metrics as well as equitable skill scores in forecast evaluation.

17 WIND ENERGY↗

Development of the MCNP-ORIGEN activation automation tool

This paper introduces the MCNP-ORIGEN activation automation tool for streamlining the calculation of experiment source terms. This tool couples the Monte-Carlo radiation-transport solver, MCNP to the depletion tool, ORIGEN-S. To showcase its current capabilities, this paper presents an activation analysis exercise, which is conducted with Serpent, and the results are com- pared. The experiment consists in a 90%-enriched uranium sphere surrounded by light water. A non-fissile cylinder is placed in the water, representing an irradiation experiment. We conducted simulations for 2 different cylinder materials: iron and cobalt at two different temperatures: 300 K and 900 K. The exercise consists of an irradiation of 50 days at a constant power of 5 MW, followed by a decay of 50 days. These studies highlight the discrepancies between how the different tools handle nuclear data. Finally, the results underscore multiple areas of possible improvement.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Two XTE A01 Projects: A Multifrequency Study of Circinus X-1 and a Search for Microsecond Variability From Bright Galactic X-Ray Source

This final report describes the research of a single common portion of the above-named two projects, conducted by G. Jernigan, i.e., the theory for a new method, a variation of a Kolmogorov-Smirnov test, for determining the fastest variability present in an X-ray source. The current phase involves testing the newly developed code on real example sources (CYG X1). Unfortunately, there are no calibration sources for testing the code, which therefore required the development of an X-ray source simulation code. The goal is to evaluate the sensitivity of the code for the detection of a range of different types of variability (bursts, pulsations, etc.).

Jernigan, Garrett↗

The Particle Accelerator Simulation Code PyORBIT

The particle accelerator simulation code PyORBIT is presented. The structure, implementation, history, parallel and simulation capabilities, and future development of the code are discussed. The PyORBIT code is a new implementation and extension of algorithms of the original ORBIT code that was developed for the Spallation Neutron Source accelerator at the Oak Ridge National Laboratory. The PyORBIT code has a two level structure. The upper level uses the Python programming language to control the flow of intensive calculations performed by the lower level code implemented in the C++ language. The parallel capabilities are based on MPI communications. The PyORBIT is an open source code accessible to the public through the Google Open Source Projects Hosting service.

Accelerator Simulation↗

An independent analysis of bias sources and variability in wind plant pre-construction energy yield estimation methods

The wind resource assessment community has long had the goal of reducing the bias between wind plant pre-construction energy yield assessment (EYA) and the observed annual energy production (AEP). This comparison is typically made between the 50% probability of exceedance (P50) value of the EYA and the long-term corrected operational AEP (hereafter OA P50), and is known as the P50 bias. The industry has critically lacked an independent analysis of bias reduction investigated across multiple consultants to identify the greatest sources of uncertainty and variance in the EYA process and the best opportunities for uncertainty reduction. The present study addresses this gap by benchmarking consultant methodologies against each other and against operational data at a scale not seen before in industry collaborations. We consider data from 10 wind plants and evaluate discrepancies between eight consultancies in the steps taken from estimates of gross to net energy. Consultants tend to overestimate the gross energy produced at the turbines and then compensate by further overestimating downstream losses, leading to a mean P50 bias near zero, still with significant variability among the individual wind plants. Within our data sample, we find that consultant estimates of all loss categories, except environmental losses, tend to reduce the project-to-project variability of the P50 bias. The disagreement between consultants, however, remains flat throughout the addition of losses. Finally, we find that differences in consultants’ estimates of project performance can lead to differences up to $10/MWh in the levelized cost of energy for a wind plant.

Todd, Austin C.↗

XScan: An Integrated Tool for Understanding Open Source Community-Based Scientific Code

Many scientific communities have adopted community-based models that integrate multiple components to simulate whole system dynamics. The community software projects’ complexity, stems from the integration of multiple individual software components that were developed under different application requirements and various machine architectures, has become a challenge for effective software system understanding and continuous software development. The paper presents an integrated software toolkit called X-ray Software Scanner (in abbreviation, XScan) for a better understanding of large-scale community-based scientific codes. Our software tool provides support to quickly summarize the overall information of scientific codes, including the number of lines of code, programming languages, external library dependencies, as well as architecture-dependent parallel software features. The XScan toolkit also realizes a static software analysis component to collect detailed structural information and provides an interactive visualization and analysis of the functions. We use a large-scale community-based Earth System Model to demonstrate the workflow, functions and visualization of the toolkit. We also discuss the application of advanced graph analytics techniques to assist software modular design and component refactoring.

Zheng, Weijian↗

Economic assessment of seismic monitoring for underground hydrogen storage

Underground hydrogen storage (UHS) plays a key role in the energy landscape. However, like other subsurface engineering technologies, UHS may cause leakage into the groundwater or atmosphere and possibly induce local seismicity. To reduce these risks, seismic monitoring could be a viable technique to track the UHS plume, detect leakages, and locate induced seismicity events. Seismic monitoring has been proposed to safely monitor UHS, but research in this area is still new and requires field studies. Lab and theoretical studies have demonstrated the validity of seismic monitoring for UHS. Therefore, it is imperative to analyze the economic feasibility of seismic monitoring for UHS. Hence, we develop a cost model and open-source Python code for seismic monitoring that considers types of seismometers, comprehensive operational scenarios, detection thresholds, and long-term leakage monitoring. A case study is further provided to validate the cost model on reservoir simulations of UHS. We find that the levelized cost for a 10-year operating UHS site will range on the order of ∼0.003 $\$$/kg. The methods developed in this study could also be applied to the monitoring of groundwater, gas, and/or wastewater injection.

08 HYDROGEN↗

Emulation of radiation transport in 3D stochastic media using 1D planar Monte Carlo stochastic media radiation transport algorithms

A subset of stochastic media radiation transport problems involves those in which radiation is incident on a thin slab of stochastic material. Particle tracking in 3D for such problems is expensive, and 1D planar models lack accuracy because they only allow the material to change in one dimension. Therefore, we propose dimensional emulation, which through a slight modification allows existing 1D planar geometry stochastic media radiation transport models to reproduce results from the equivalent 3D models by allowing the material to change in all three dimensions, reproducing the fidelity of the 3D model for the low computational cost of the 1D planar model. In this work, we apply dimensional emulation to three Monte Carlo stochastic media radiation transport models: Chord Length Sampling (CLS), the Local Realization Preserving method (LRP), and a variant of Conditional Point Sampling (CoPS). For a common Markovian benchmark set, the 3D emulation variants of these algorithms are numerically verified to reproduce the results of the 3D variants within statistics while running 1.3 to 2 times faster in the implementation within Sandia National Laboratories open-source research code PlaybookMC. The 3D emulation variants are also shown to yield a 72%–92% reduction in error for the thin slab problems in comparison to the 1D benchmark. As a result, the 3D emulation variant of CLS and CoPS-1 are shown to reproduce 3D CLS results that were used to approximate results for a 3D spherical inclusion geometry benchmark set.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

MaterialsMap: A CALPHAD-based tool to design composition pathways through feasibility map for desired dissimilar materials, demonstrated with resistance spot welding joining of Ag-Al-Cu

Assembly of dissimilar metals can be achieved by different methods, for example, casting, welding, and additive manufacturing (AM). However, undesired phases formed in liquid-phase assembling processes due to solute segregation during solidification diminish mechanical and other properties of the processed parts. In the present work, an open-source software named MaterialsMap, has been developed based on the CALculation of Phase Diagrams (CALPHAD) approach. The primary objective of MaterialsMap is to facilitate the design of an optimal composition pathway for assembling dissimilar alloys with liquid-phases based on the formation of desired and undesired phases along the pathway. In MaterialsMap, equilibrium thermodynamic calculations are used to predict equilibrium phases formed at slow cooling rate, while Scheil-Gulliver simulations are employed to predict non-equilibrium phases formed during rapid cooling. Additionally, by combining these two simulations, MaterialsMap offers a thorough guide for understanding phase formation in various manufacturing processes, assisting users in making informed decisions during material selection and production. As a demonstration of this approach, a compositional pathway was designed from pure Al to pure Cu through Ag using MaterialsMap. The design was experimentally verified using resistance spot welding (RSW).

36 MATERIALS SCIENCE↗

Estimations of Mo X-pinch plasma parameters on QiangGuang-1 facility by L-shell spectral analyses

Plasma parameters of molybdenum (Mo) X-pinches on the 1-MA QiangGuang-1 facility were estimated by L-shell spectral analysis. X-ray radiation from X-pinches had a pulsed width of 1 ns, and its spectra in 2–3 keV were measured with a time-integrated X-ray spectrometer. Relative intensities of spectral features were derived by correcting for the spectral sensitivity of the spectrometer. With an open source, atomic code FAC (flexible atomic code), ion structures, and various atomic radiative-collisional rates for O-, F-, Ne-, Na-, Mg-, and Al-like ionization stages were calculated, and synthetic spectra were constructed at given plasma parameters. By fitting the measured spectra with the modeled, Mo X-pinch plasmas on the QiangGuang-1 facility had an electron density of about 10{sup 21} cm{sup −3} and the electron temperature of about 1.2 keV.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗