Search NASA⌕ Search

SEARCH · Search NASA

Results for “Source code instrumentation”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

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

Analysis of DESI × DES using the Lagrangian effective theory of LSS

In this work we use Lagrangian perturbation theory to analyze the harmonic space galaxy clustering signal of the Bright Galaxy Survey (BGS) and luminous red galaxies (LRGs) targeted by the dark energy spectroscopic instrument (DESI), combined with the galaxy-galaxy lensing signal measured around these galaxies using Dark Energy Survey Year 3 source galaxies. The BGS and LRG galaxies are extremely well characterized by DESI spectroscopy and, as a result, lens galaxy redshift uncertainty and photometric systematics contribute negligibly to the error budget of our “2 × 2-point” analysis. On the modeling side, this work represents the first application of the SPINOSAURUS code, implementing an effective field theory model for galaxy intrinsic alignments, and we additionally introduce a new scheme ( MAIAR ) for marginalizing over the large uncertainties in the redshift evolution of the intrinsic alignment signal. Furthermore, this is the first application of a hybrid effective field theory model for galaxy bias based on the A emulus 𝜈 simulations. Our main result is a measurement of the amplitude of the lensing signal, 𝑆 8 = 𝜎 8 ⁢(Ω 𝑚 /0.3) 0.5 = 0.85⁢0$^{+0.042}_{−0.050}$, consistent with values of this parameter derived from the primary cosmic microwave background. This constraint is artificially improved by a factor of 51% if we assume a more standard, but restrictive parametrization for the redshift evolution and sample dependence of the intrinsic alignment signal, and 63% if we additionally assume the nonlinear alignment model. Furthermore, we show that when fixing the cosmological model to the best-fit values from Planck PR4 there is > 5⁢𝜎 evidence for a deviation of the evolution of the intrinsic alignment signal from the functional form that is usually assumed in cosmic shear and galaxy-galaxy lensing studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Particle transport and heat loads in NIO1

NIO1 is a compact radio frequency ion source designed to generate a 60 kV-135 mA hydrogen negative ion beam and it aims at continuous operation, which implies a detailed thermo-mechanical analysis of the beam-facing components, in particular, the accelerator grids. A 3D analysis of the entire NIO1 beam has been performed for the first time with a fully 3D version of EAMCC, a relativistic particle tracking code for the calculation of the grid power deposition induced by particle impacts. According to the results presented in this paper, secondary and co-extracted electrons cause a non-negligible heat load on the grids, where different high-power density regions, within reasonable sustainable standard limits, are calculated.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Electron transport in gaseous detectors with a Python-based Monte Carlo simulation code

Understanding electron drift and diffusion in gases and gas mixtures is a topic of central importance for the development of modern particle detection instrumentation. The industry-standard MagBoltz code has become an invaluable tool during its 20 years of development, providing capability to solve for electron transport (‘swarm’) properties based on a growing encyclopedia of built-in collision cross sections. We have made a refactorization of this code from FORTRAN into Cython, and studied a range of gas mixtures of interest in high energy and nuclear physics. The results from the new open source PyBoltz package match the outputs from the original MagBoltz code, with comparable simulation speed. An extension to the capabilities of the original code is demonstrated, in implementation of a new Modified Effective Range Theory interface. We hope that the versatility afforded by the new Python code-base will encourage continued use and development of the MagBoltz tools by the particle physics community.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

SPIDER beam dump as diagnostic of the particle beam

The beam power produced by the negative ion source for the production of ion of deuterium extracted from RF plasma is mainly absorbed by the beam dump component which has been designed also for measuring the temperatures on the dumping panels for beam diagnostics. A finite element code has been developed to characterize, by thermo-hydraulic analysis, the sensitivity of the beam dump to the different beam parameters. The results prove the capability of diagnosing the beam divergence and the horizontal misalignment, while the entity of the halo fraction appears hardly detectable without considering the other foreseen diagnostics like tomography and beam emission spectroscopy.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Background gas density and beam losses in NIO1 beam source

NIO1 (Negative Ion Optimization 1) is a versatile ion source designed to study the physics of production and acceleration of H- beams up to 60 keV. In ion sources, the gas is steadily injected in the plasma source to sustain the discharge, while high vacuum is maintained by a dedicated pumping system located in the vessel. In this paper, the three dimensional gas flow in NIO1 is studied in the molecular flow regime by the Avocado code. The analysis of the gas density profile along the accelerator considers the influence of effective gas temperature in the source, of the gas temperature accommodation by collisions at walls, and of the gas particle mass. The calculated source and vessel pressures are compared with experimental measurements in NIO1 during steady gas injection.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Imaging suprathermal x-rays from a laboratory plasma jet using PIN-diode-based and scintillator-based 1D pinhole/coded aperture cameras

A PIN-diode-based 1D x-ray camera and a scintillator-based 1D x-ray camera, both with a microsecond to submicrosecond time resolution, have been developed to perform time-resolved imaging of transient, low-intensity, suprathermal x-rays associated with magnetohydrodynamic instabilities disrupting a plasma jet. These cameras have a high detection efficiency over a broad x-ray band, a wide field of view, and the capability to produce >50 time-resolved frames with a ≤1 μs time resolution. The x-ray images are formed by a pinhole or by a coded aperture placed outside a vacuum chamber in which the plasma jet is launched. The 1D imaging shows that the location of the x-ray source is either a few centimeters away from an inner disk electrode or near a spatially translatable metal frame that is 30–40 cm away from the electrode. Compared to a pinhole, a coded aperture increases the signal collection efficiency but also introduces unwanted artifacts.

Instruments & Instrumentation↗
DOE OSTI2018

High-resolution Gamma Spectroscopy Measurements of Pressurized Water Reactor Spent Nuclear Fuel Rods

International nuclear safeguards use nondestructive assay techniques to verify declared characteristics of spent nuclear fuel (SNF) assemblies. The declared characteristics usually include the initial enrichment, assembly average burnup, and irradiation history. The primary nondestructive assay tool used by the International Atomic Energy Agency has been the fork detector irradiated fuel measurement system, which uses gross gamma and neutron correlations; however, there have been several research publications within the last decade on the use of high-resolution gamma spectroscopy for the determination of SNF characteristics. References [1]–[6] are published journal articles on this subject from 2006–2016. This report describes measurements of 25 individual SNF rods using high-resolution gamma spectroscopy. Details of the measurement geometry, measurement settings, and rod characteristics are documented. These spectra are intended for the validation of simulation codes that calculate SNF inventory and source term characteristics. The spectrum archive is available to evaluate future spectrum analysis methods. The archive contains full axial scan spectra of each SNF rod. This allows more detailed comparison between simulation and nondestructive analysis than is possible on assemblies.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Balmer-α spectrum measurements of the LHD one-third ion source

Wavelength spectra of Balmer-α light from plasmas in the extraction region of the Large Helical Device-R&D negative ion source, or the LHD one-third ion source have exhibited a blue shift as a negative bias voltage was applied to the plasma grid. The blue shift increased as the negative bias voltage with respect to the local plasma potential was increased. The measured spectra were compared with the velocity distributions of surface reflected hydrogen atoms calculated by atomic collisions in amorphous target code. The arc power and the source H{sub 2} pressure also affected the shift and broadening in the observed Balmer-α spectra. The possibility of identifying the negative hydrogen ions produced at the low work function plasma grid surface by high resolution spectroscopy is discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

System Advisor Model (SAM) General Description (Version 2017.9.5)

This document describes the capabilities of the System Advisor Model (SAM) developed and distributed by the U.S. Department of Energy's National Renewable Energy Laboratory. The document is for potential users and others wanting to learn about the model's capabilities. SAM is a techno-economic computer model that calculates performance and financial metrics of renewable energy projects. Project developers, policy makers, equipment manufacturers, and researchers use graphs and tables of SAM results in the process of evaluating financial, technology, and incentive options for renewable energy projects. SAM simulates the performance of photovoltaic, concentrating solar power, solar water heating, wind, geothermal, biomass, and conventional power systems. The financial models are for projects that either buy and sell electricity at retail rates (residential and commercial) or sell electricity at a price determined in a power purchase agreement (PPA). SAM's simulation tools facilitate parametric and sensitivity analyses, Monte Carlo simulation and weather variability (P50/P90) studies. SAM can also read input variables from Microsoft Excel worksheets. For software developers, the SAM software development kit (SDK) makes it possible to use SAM simulation modules in their applications written in C/C plus plus, C sharp, Java, Python, MATLAB, and other languages. NREL provides both SAM and the SDK as free downloads at https://sam.nrel.gov. SAM is an open source project, so its source code is available to the public. Researchers can study the code to understand the model algorithms, and software programmers can contribute their own models and enhancements to the project. Technical support and more information about the software are available on the website.

47 OTHER INSTRUMENTATION↗
DOE OSTI2019

Simulations of particle interactions in a high-current RFQ

High-intensity beams of composite ions are susceptible to particle interactions that do not affect single ions, such as stripping of electrons and charge exchange. Beam dynamics simulation techniques used in accelerator design do not often include these particle interactions, and so the transmission of the real beam can be significantly different from the design simulations. We have modified the IMPACT-T code to include particle interactions as part of the beam dynamics simulations, including space-charge effects and beam losses. We have carried out validation simulations of interactions in a molecular hydrogen ion beam through an existing high-intensity deuteron RFQ design, as beam commissioning is often carried out using molecular hydrogen beams, to avoid deuteron-deuteron interactions in the accelerator leading to neutron activation. The new code framework successfully models the interactions within the beam, tracks the different source and product ions through the accelerator, and calculates the collective effects such as space charge from all particles together. Finally, we also discuss other simulation codes that could be modified to include similar interaction effects, and possible applications to other types of accelerators.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Sub-Sampled Imaging for STEM: Maximising Image Speed, Resolution and Precision Through Reconstruction Parameter Refinement

Sub-sampling during image acquisition in scanning transmission electron microscopy (STEM) has been shown to provide a means to increase the overall speed of acquisition while at the same time providing an efficient means to control the dose, dose rate and dose overlap delivered to the sample. In this paper, we discuss specifically the parameters used to reconstruct sub-sampled images and highlight their effect on inpainting using the beta-process factor analysis (BPFA) methodology. The selection of the main control parameters can have a significant effect on the resolution, precision and sensitivity of the final inpainted images, and here we demonstrate a method by which these parameters can be optimised for any image in STEM. As part of this paper, we also provide a link to open source code and a tutorial on its use, whereby these parameters can be tested for any datasets. When coupled with the hardware necessary to rapidly sub-sample images in STEM, this approach can have significant implications for imaging beam sensitive materials and dynamic processes.

47 OTHER INSTRUMENTATION↗
DOE OSTI2015

ABSORPTION MEASURE DISTRIBUTION IN Mrk 509

In this paper we model the observed absorption measure distribution (AMD) in Mrk 509, which spans three orders of magnitude in ionization level with a single-zone absorber in pressure equilibrium. AMD is usually constructed from observations of narrow absorption lines in radio-quiet active galaxies with warm absorbers. We study the properties of the warm absorber in Mrk 509 using recently published broadband spectral energy distribution observed with different instruments. This spectrum is an input in radiative transfer computations with full photoionization treatment using the titan code. We show that the simplest way to fully reproduce the shape of AMD is to assume that the warm absorber is a single zone under constant total pressure. With this assumption, we found theoretical AMD that matches the observed AMD determined on the basis of the 600 ks reflection grating spectrometer XMM-Newton spectrum of Mrk 509. The softness of the source spectrum and the important role of the free–free emission breaks the usual degeneracy in the ionization state calculations, and the explicit dependence of the depths of AMD dips on density open a new path to the density diagnostic for the warm absorber. In Mrk 509, the implied density is of the order of 10{sup 8} cm{sup −3}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Implications of 3D Thermal Structure on 1D Atmospheric Retrieval

Using the atmospheric structure from a 3D global radiation-hydrodynamic simulation of HD 189733b and the open-source Bayesian Atmospheric Radiative Transfer (BART) code, we investigate the difference between the secondary-eclipse temperature structure produced with a 3D simulation and the best-fit 1D retrieved model. Synthetic data are generated by integrating the 3D models over the Spitzer , the Hubble Space Telescope ( HST ), and the James Web Space Telescope ( JWST ) bandpasses, covering the wavelength range between 1 and 11 μ m where most spectroscopically active species have pronounced features. Using the data from different observing instruments, we present detailed comparisons between the temperature–pressure profiles recovered by BART and those from the 3D simulations. We calculate several averages of the 3D thermal structure and explore which particular thermal profile matches the retrieved temperature structure. We implement two temperature parameterizations that are commonly used in retrieval to investigate different thermal profile shapes. To assess which part of the thermal structure is best constrained by the data, we generate contribution functions for our theoretical model and each of our retrieved models. Our conclusions are strongly affected by the spectral resolution of the instruments included, their wavelength coverage, and the number of data points combined. We also see some limitations in each of the temperature parametrizations, as they are not able to fully match the complex curvatures that are usually produced in hydrodynamic simulations. The results show that our 1D retrieval is recovering a temperature and pressure profile that most closely matches the arithmetic average of the 3D thermal structure. When we use a higher resolution, more data points, and a parametrized temperature profile that allows more flexibility in the middle part of the atmosphere, we find a better match between the retrieved temperature and pressure profile and the arithmetic average. The Spitzer and HST simulated observations sample deep parts of the planetary atmosphere and provide fewer constraints on the temperature and pressure profile, while the JWST observations sample the middle part of the atmosphere, providing a good match with the middle and most complex part of the arithmetic average of the 3D temperature structure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Forward Modeling of Large-scale Structure: An Open-source Approach with Halotools

Here, we present the first stable release of Halotools (v0.2), a community-driven Python package designed to build and test models of the galaxy-halo connection. Halotools provides a modular platform for creating mock universes of galaxies starting from a catalog of dark matter halos obtained from a cosmological simulation. The package supports many of the common forms used to describe galaxy-halo models: the halo occupation distribution (HOD), the conditional luminosity function (CLF), abundance matching, and alternatives to these models that include effects such as environmental quenching or variable galaxy assembly bias. Satellite galaxies can be modeled to live in subhalos, or to follow custom number density profiles within their halos, including spatial and/or velocity bias with respect to the dark matter profile. Here, the package has an optimized toolkit to make mock observations on a synthetic galaxy population, including galaxy clustering, galaxy-galaxy lensing, galaxy group identification, RSD multipoles, void statistics, pairwise velocities and others, allowing direct comparison to observations. Halotools is object-oriented, enabling complex models to be built from a set of simple, interchangeable components, including those of your own creation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

Cosmic Microwave Background Anisotropy Numerical Solution (CMBAns) I: An Introduction to $C_l$ Calculation

Cosmological Boltzmann codes are often used by researchers for calculating the CMB angular power spectra from different theoretical models, for cosmological parameter estimation, etc. Therefore, the accuracy of a Boltzmann code is of utmost importance. Different Markov Chain Monte Carlo based parameter estimation algorithms typically require 10^3 - 10^4 iterations of Boltzmann code. This makes the time complexity of such codes another critical factor. In the last two decades, several Boltzmann packages, such as CMBFAST, CAMB, CMBEasy, CLASS etc., have been developed. In this paper, we present a new cosmological Boltzmann code, CMBAns, that can be used for accurate calculation of the CMB power spectrum. At present, CMBAns is developed for a flat background matrix. It is mostly written in the C language. However, we borrowed the concept of class from C++. This gives researchers the flexibility to develop their own independent package based on CMBAns, without an in-depth understanding of the source code. We also develop multiple stand-alone facilities which can be directly compiled and run on a given parameter set. In this paper, we discuss all the mathematical formulation, approximation schemes, integration methods etc., that are used in CMBAns. The package will be made available through github for public use in the near future.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Study of Data-Driven Mesh-Model Optimization in System Thermal-Hydraulic Simulation

Quantification of nuclear power plant safety risk requires a systematic and yet practical approach to identification of accident scenarios, assessment of their likelihood and consequences. Instrumental to this goal is risk-informed safety margin characterization (RISMC) framework, whose realization requires computationally robust and affordable methods for sufficiently accurate simulation of complex multi-dimensional physical phenomena, such as turbulent and multi-phase flow. The CFD-like correlation-based codes with 3D simulation capability (e.g., GOTHIC) ensure computational efficiency using coarse mesh size and the sub-grid phenomena in the boundary layer that can be captured by adequate constitutive correlations. However, the error sources and user effects on the selection of mesh size and models lead to unpredictable simulation error, while rich High-Fidelity (HF) data from experiments and numerical simulation using validated code or DNS are not fully explored. It would be useful to have a “smart” data-driven multi-scale framework in which the low-resolution models can be “taught” to emulate high-resolution models. The objective of this work is to develop and evaluate a physical-based data-driven mesh-model optimization approach (Optimized Mesh/Model Information System, OMIS) to estimate the simulation error and give advice on the optimized selection of coarse mesh size and Low-Fidelity (LF) models for System Thermal Hydraulic (STH) simulation to achieve accuracy comparable to that of HF models. This approach takes advantages of computational efficiency of coarse-mesh simulation and application of Machine Learning (ML) algorithms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Alternative modeling methods for plasma-based Rf ion sources

Rf-driven ion sources for accelerators and many industrial applications benefit from detailed numerical modeling and simulation of plasma characteristics. For instance, modeling of the Spallation Neutron Source (SNS) internal antenna H{sup −} source has indicated that a large plasma velocity is induced near bends in the antenna where structural failures are often observed. This could lead to improved designs and ion source performance based on simulation and modeling. However, there are significant separations of time and spatial scales inherent to Rf-driven plasma ion sources, which makes it difficult to model ion sources with explicit, kinetic Particle-In-Cell (PIC) simulation codes. In particular, if both electron and ion motions are to be explicitly modeled, then the simulation time step must be very small, and total simulation times must be large enough to capture the evolution of the plasma ions, as well as extending over many Rf periods. Additional physics processes such as plasma chemistry and surface effects such as secondary electron emission increase the computational requirements in such a way that even fully parallel explicit PIC models cannot be used. One alternative method is to develop fluid-based codes coupled with electromagnetics in order to model ion sources. Time-domain fluid models can simulate plasma evolution, plasma chemistry, and surface physics models with reasonable computational resources by not explicitly resolving electron motions, which thereby leads to an increase in the time step. This is achieved by solving fluid motions coupled with electromagnetics using reduced-physics models, such as single-temperature magnetohydrodynamics (MHD), extended, gas dynamic, and Hall MHD, and two-fluid MHD models. We show recent results on modeling the internal antenna H{sup −} ion source for the SNS at Oak Ridge National Laboratory using the fluid plasma modeling code USim. We compare demonstrate plasma temperature equilibration in two-temperature MHD models for the SNS source and present simulation results demonstrating plasma evolution over many Rf periods for different plasma temperatures. We perform the calculations in parallel, on unstructured meshes, using finite-volume solvers in order to obtain results in reasonable time.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗