Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation output analysis”

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

Updated status of the ART neutronic fast reactor tools integration to the NEAMS Workbench

The Workbench initiative was launched in FY-2017 within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) Integration Product Line to facilitate the transition from conventional tools to high-fidelity tools. The Workbench provides a common user interface for model creation, real-time validation, execution, output processing, and visualization for integrated codes. The integration of the Argonne Reactor Computation (ARC) suite of codes into the NEAMS Workbench was initiated in FY-2017. The ARC codes contain both legacy codes like DIF3D and REBUS-3 that were developed with over 30 years of experience, and newer NEAMS additions like MC 2 -3 and PERSENT. The ARC integration into the NEAMS Workbench interface relies on the PyARC module which handles the pre- and post-processing of the native ARC codes input, and the runtime environment. The PyARC module together with the NEAMS Workbench interface are both released under Open Source Software licenses. This report describes the ARC capabilities available with the Workbench at the end of FY-2018, and a tutorial is provided together with the automatically-generated code documentation. Integrating the ARC codes into the Workbench benefits directly the Advanced Reactor Technology (ART) Campaign by providing a set of controlled, maintained and validated scripts to generate ARC inputs, which promotes best practices, reduces the learning curve, and facilitates project collaboration. The second benefit from this project is to improve the user experience with the ARC codes: the Workbench interface provides assistance for building an input through auto-completion, real-time validation, access to templates, and visualization. Finally, the ARC integration into the NEAMS Workbench enables new modeling capabilities, those will be extremely useful for advanced reactor design and analyses. The PyARC module facilitates and automatizes complex calculations and workflows for reactor analysis enabling geometrical perturbations, cross-section update through depletion, etc. The Dakota/PyARC coupling in the Workbench was also demonstrated to enable mathematical optimization and sensitivity analysis/uncertainty quantification (SA/UQ) techniques with ARC neutronic simulations. For demonstration purposes, the ARC codes were used through the Workbench for solving the Sodium-cooled Fast Reactor Uncertainty Analysis in Modelling (SFR-UAM) benchmark problems. The ARC codes are currently used at ANL through the Workbench by nuclear engineers for SFR core design analyses. Future efforts will focus on continuously adding modeling capabilities available with the ARC codes, training new users to continue building some experience, and on integrating high-fidelity codes such as PROTEUS into the Workbench.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Machine learning changes the rules for flux limiters

Learning to integrate non-linear equations from highly resolved direct numerical simulations has seen recent interest for reducing the computational load for fluid simulations. Here, we focus on determining a flux-limiter for shock capturing methods. Focusing on flux limiters provides a specific plug-and-play component for existing numerical methods. Since their introduction, an array of flux limiters has been designed. Using the coarse-grained Burgers' equation, we show that flux-limiters may be rank-ordered in terms of their log-error relative to high-resolution data. We then develop a theory to find an optimal flux-limiter and present flux-limiters that outperform others tested for integrating Burgers' equation on lattices with [Formula: see text], and 2x, 3x, 4x, and 8x coarse-grainings. We train a continuous piecewise linear limiter by minimizing the mean-squared misfit to six-grid point segments of high-resolution data, averaged over all segments. While flux limiters are generally designed to have an output of φ(r) = 1 at a flux ratio of r = 1, our limiters are not bound by this rule and yet produce a smaller error than standard limiters. Here we find that our machine learned limiters have distinctive features that may provide new rules-of-thumb for the development of improved limiters. Additionally, we use our theory to learn flux-limiters that outperform standard limiters across a range of values (as opposed to at a specific fixed value) of coarse-graining, number of discretized bins, and diffusion parameter. This demonstrates the ability to produce flux limiters that should be more broadly useful than standard limiters for general applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Temperature-Dependent Characterization, Modeling, and Switching Speed-Limitation Analysis of Third-Generation 10-kV SiC MOSFET

The temperature-dependent characteristics of the third-generation 10-kV/20-A SiC MOSFET including the static characteristics and switching performance are carried out in this paper. The steady-state characteristics, including saturation current, output characteristics, antiparallel diode, and parasitic capacitance, are tested. Here, a double pulse test platform is constructed including a circuit breaker and gate drive with >10-kV insulation and also a hotplate under the device under test for temperature-dependent characterization during switching transients. The switching performance is tested under various load currents and gate resistances at a 7-kV dc-link voltage from 25 to 125 C and compared with previous 10-kV MOSFETs. A simple behavioral model with its parameter extraction method is proposed to predict the temperature-dependent characteristics of the 10-kV SiC MOSFET. The switching speed limitations, including the reverse recovery of SiC MOSFET's body diode, overvoltage caused by stray inductance, crosstalk, heat sink, and electromagnetic interference to the control are discussed based on simulations and experimental results.

10-kV silicon carbide (SiC) metal–oxide–semiconduc↗
DOE OSTI2019

A Study into the Benefits of Using Discrete Event Simulation in TRU and Mixed Low-Level Waste (MLLW) Disposal at Advanced Mixed Waste Treatment Plant (AMWTP) - 19536

As part of the 1995 Settlement agreement reached between the Idaho State Government and the Idaho National laboratory, it was agreed that the approximately 65,000 cubic metres of waste stored at the AMWTP site was to be shipped out of the state to the waste isolation pilot plant (WIPP) in New Mexico or a similar disposal facility by the end of December 2018. In June 2016, the operation of AMWTP was taken over by Fluor Idaho, a subsidiary of Fluor Inc. To improve their capability of processing all of the varied waste contained on the extensive site, Fluor Idaho requested the assistance of DBD Ltd. We were tasked with producing an all-encompassing operational research (OR) model/discrete event simulation (DES) model of the site, which is able to produce live end-date predictions based on current inventory and waste tracking against milestones. The DES model, is capable of predicting years worth of plant operation in a matter of minutes including the variation caused by such complexities as resource/staff availability, equipment breakdown and maintenance, prioritisation of different waste types, and varying shift patterns of both staff and facilities. The DES model can be run from its own bespoke dashboard, that can be customised to suit the audience it is being used to inform, producing clear key performance indicator (KPI) graphics and high-level overviews for senior management, down to drum by drum breakdowns of the predicted feed plans for each individual facility on the AMWTP site for a work planner to analyse and compare to the current operations. Simulation can produce these types of outputs for very low cost and low risk, as it avoids experimentation to find potential alternatives to the status-quo that could result in loss of production or time. These methods can also be used to produce cost benefit analyses for suggested processing options feeding into optioneering studies and cost reduction strategies. As a tool for plant optimisation and tracking, discrete event simulation can provide detailed results quickly and reliably. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

Regional hydrologic response to climate change in the conterminous United States using high-resolution hydroclimate simulations

Despite the fact that Global Climate Model (GCM) outputs have been used to project hydrologic impacts of climate change using off-line hydrologic models for two decades, many of these efforts have been disjointed applications or at least calibrations have been focused on individual river basins and using a few of the available GCMs. This study improves upon earlier attempts by systematically projecting hydrologic impacts for the entire conterminous United States (US), using outputs from ten GCMs from the latest Coupled Model Intercomparison Project phase 5 (CMIP5) archive, with seamless hydrologic model calibration and validation techniques to produce a spatially and temporally consistent set of current hydrologic projections. The Variable Infiltration Capacity (VIC) model was forced with ten-member ensemble projections of precipitation and air temperature that were dynamically downscaled using a regional climate model (RegCM4) and bias-corrected to 1/24 (~4 km) grid resolution for the baseline (1966 2005) and future (2011 2050) periods under the Representative Concentration Pathway 8.5. Based on regional analysis, the VIC model projections indicate an increase in winter and spring total runoff due to increases in winter precipitation of up to 20% in most regions of the US. However, decreases in snow water equivalent (SWE) and snow-covered days will lead to significant decreases in summer runoff with more pronounced shifts in the time of occurrence of annual peak runoff projected over the eastern and western US. In contrast, the central US will experience year-round increases in total runoff, mostly associated with increases in both extreme high and low runoff. Furthermore, the projected hydrological changes described in this study have implications for various aspects of future water resource management, including water supply, flood and drought preparation, and reservoir operation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Local and Global Sensitivity Analysis of a Reactive Transport Model Simulating Floodplain Redox Cycling

Reactive transport models (RTMs) are essential tools that simulate the coupling of advective, diffusive, and reactive processes in the subsurface, but their complexity makes them difficult to understand, develop and improve without accompanying statistical analyses. Although global sensitivity analysis (SA) can address these issues, the computational cost associated with most global SA techniques limits their use with RTMs. In this study, we apply distance-based generalized sensitivity analysis (DGSA), a novel and computationally efficient method of global SA, to a floodplain-scale RTM and compare DGSA results to those from local SA. Our test case focuses on the impact of 17 uncertain environmental parameters on spatially and temporally variable redox conditions within a floodplain aquifer. The input parameters considered include flow and diffusion rates, geochemical reaction rates, and the spatial distribution of sediment facies. Sensitivity was evaluated for three distinct components of the model response, encompassing both multidimensional and categorical output. Parameter rankings differ between local SA and DGSA, due to nonlinear effects of individual parameters and interaction effects between parameters. DGSA results show that fluid residence time, which is controlled by aquifer permeability, generally exerts a stronger control on redox conditions than do geochemical reaction rates. Sensitivity indices also demonstrate that sulfate reduction is key for establishing and maintaining reducing conditions throughout the aquifer. Furthermore, these results provide insights into the key drivers of heterogeneous redox processes within floodplain aquifers, as well as the main sources of uncertainty when modeling complex subsurface systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Enabling Load Following Capability in the Transatomic Power MSR

This project is dedicated towards designing a fuel processing system that enables liquid-fueled molten salt reactors (MSR) to load follow by removing the dissolved xenon in the fuel salt. As one of the Gen-IV nuclear reactor concepts, the molten salt reactor receives increasing development interests in the recent years. One distinguishing feature of the liquid-fueled molten salt reactor is its improved ability to operate in a load-following mode by including the unique online fission product removal system. Load-following means that the reactor changes its power output based on the demand on the grid. Most of the current operating nuclear reactors have limited load-following ability and operate as the base load on the grid. Due to the rapid increase of solar energy, the requirement on load-following capacity is significantly increased because of the varying power output of the solar panels, yet the traditional load-following capacity is expected to decrease as the decarbonization of the grid continues. Therefore, the ability to perform load-following operation for the nuclear reactors will greatly enhance the resilience of the grid and make nuclear energy more economically competitive. This load-following feature is included in many commercial molten salt reactor designs, such as the designs by Transatomic Power, Terrestrial Energy, and Flibe Energy. Unfortunately, detailed analysis of the fuel processing system for commercial scale MSRs is still lacking, as well as how the fuel processing quantitively impacts the load-following operation. Moreover, experimental data for many of the underlying physics of fuel processing is limited. This project aims to pave the way for the fuel processing technology to advance to the commercial stage by performing combined experimental and simulation research. During the project period, four interconnected aspects of the development of the fuel processing system in liquid-fueled molten salt reactors are investigated. These aspects are the simulation and analysis of the fission product removal system, the fuel cycle simulation, the coupled reactor neutronics and thermal hydraulics transient simulation, and the gaseous fission product removal experiment. Multiphase CFD simulations are performed for components of the processing systems, and simplified air-water experiments are carried out to provide validation data. It is concluded that the CFD simulation can satisfactorily predict the system level performance of the components, and engineering models are constructed based on this success. Fuel cycle analysis is performed for two representative MSR design, the MSBR and the Transatomic Power MSR. Open-source code SaltProc is developed to incorporate the unique fuel processing system of the MSRs. It is concluded that the removal of xenon is essential for load-following operation in thermal spectrum MSR and Molten Salt Breeder Reactor. For the Transatomic Power MSR, the xenon poisoning effect is negligible due to its relatively fast neutron spectrum, though the overall fuel cycle economics still benefits from the removal of xenon. Coupled reactor neutronics and thermal hydraulics transient simulation is performed specifically for the Transatomic Power MSR. It is concluded that the reactor core design could perform power ramping fast enough to satisfy load-following operation. Combining the findings from each aspect, it is concluded that the load-following operation of a thermal neutron MSR is dependent upon the removal of xenon, which could be achieved for a commercial sized reactor using continuous inert gas sparging in a separate system with reasonable dimensions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1993

Computer program for analysis of high speed, single row, angular contact, spherical roller bearing, SASHBEAN. Volume 1: User's guide

The computer program SASHBEAN (Sikorsky Aircraft Spherical Roller High Speed Bearing Analysis) analyzes and predicts the operating characteristics of a Single Row, Angular Contact, Spherical Roller Bearing (SRACSRB). The program runs on an IBM or IBM compatible personal computer, and for a given set of input data analyzes the bearing design for it's ring deflections (axial and radial), roller deflections, contact areas and stresses, induced axial thrust, rolling element and cage rotation speeds, lubrication parameters, fatigue lives, and amount of heat generated in the bearing. The dynamic loading of rollers due to centrifugal forces and gyroscopic moments, which becomes quite significant at high speeds, is fully considered in this analysis. For a known application and it's parameters, the program is also capable of performing steady-state and time-transient thermal analyses of the bearing system. The steady-state analysis capability allows the user to estimate the expected steady-state temperature map in and around the bearing under normal operating conditions. On the other hand, the transient analysis feature provides the user a means to simulate the 'lost lubricant' condition and predict a time-temperature history of various critical points in the system. The bearing's 'time-to-failure' estimate may also be made from this (transient) analysis by considering the bearing as failed when a certain temperature limit is reached in the bearing components. The program is fully interactive and allows the user to get started and access most of its features with a minimal of training. For the most part, the program is menu driven, and adequate help messages were provided to guide a new user through various menu options and data input screens. All input data, both for mechanical and thermal analyses, are read through graphical input screens, thereby eliminating any need of a separate text editor/word processor to edit/create data files. Provision is also available to select and view the contents of output files on the monitor screen if no paper printouts are required. A separate volume (Volume-2) of this documentation describes, in detail, the underlying mathematical formulations, assumptions, and solution algorithms of this program.

Aggarwal, Arun K.↗
NASA NTRS2010

Trend Analysis and Detection of Changes in the Stratospheric Circulation

Increases in the circulation of the stratosphere appear to be a robust result of climate change in chemistry-climate models over decadal time scales. To date observations have yet to show a significant change in this circulation. It is important for the design of future observational missions to identify suitable atmospheric constituents and to determine the accuracy and length of record needed to identify a significant trend that can be attributed to circulation change. First, we determine what atmospheric variables can be used as proxies for stratospheric circulation changes. A few examples are changes in tropical lower stratospheric ozone, phase lag of the water vapor tape recorder, CO2, and SF6. Then, using both the Goddard Earth Observing System Chemistry-Climate Model (GEOS CCM) and observations from satellites and balloon soundings, we calculate the number of years needed to detect a significant trend, taking into account observational uncertainty. Model simulations will be evaluated to see how well they represent observed variability. In addition, the impacts of autocorrelation among the output or data and gaps in the observational record will be discussed.

Oman, Luke↗
DOE OSTI2015

SU-C-304-03: Experimental Investigation On the Accuracy of Plastic Scintillation Dosimeters in Small Fields

Purpose: To investigate the accuracy of the Exradin W1 (SI) and of an “in-house” plastic scintillation dosimeter (CHUQ PSD) in small radiation fields. Methods: Output factor (OF) measurements with the W1 and CHUQ PSD were performed for field sizes of 0.5 x 0.5, 1 x 1 and 2 x 2 cm{sup 2}. Both detectors were placed parallel to the central axis (CAX) in water. The spectrum discrimination calibration method was performed in each set-up to account for the Cerenkov (CRV) signal created in the fiber. The OFs were compared to the expected field factors in water derived using i) Monte Carlo (MC) simulations of an accurate accelerator model and ii) microLion (PTW) and D1V diode (SI) OFs. MC-derived correction factors were applied to both the microLion and D1V OFs. For the CHUQ PSD the calibration was repeated in water (// CAX), solid water (perpendicular to CAX) and under a shielded configuration. The signal was collected using a spectrometer (wavelength range = 185–1100 nm). Spectral analysis was performed to evaluate potential changes of the spectral distributions under the various calibration set-up configurations. Results: The W1 OFs presented an over-response for the 0.5 x 0.5 cm{sup 2} in the range of 3 – 4.1% relative to the expected field factor. The CHUQ PSD presented an under-response in the range of 1.5 – 2.7%, without accounting for volume averaging. The CRV spectra under the various calibration procedures appeared similar to each other and only minor changes were observed to the respective OFs. Conclusion: The W1 and CHUQ PSD can be used in small fields down to a 1 x 1 cm{sup 2} field size. Discrepancies were encountered between the two detectors for the smallest field size of 0.5 x 0.5 cm{sup 2} with the CHUQ PSD exhibiting a closer agreement to the expected field factor. Funding sources: 1) Alexander S. Onassis Public Benefit Foundation in Greece and 2) CREATE Medical Physics Research Training Network grant of the Natural Sciences and Engineering Research Council, Grant number: 432290, RGPIN-2014-06475.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

FTICR-MS Data from Multi-continent River Water and Sediment and from Coastal River Fresh and Saline Sediment Associated with: “Dissolved Organic Matter Functional Trait Relationships are Conserved Across Rivers”

This data package is associated with the publication “Dissolved Organic Matter Functional Trait Relationships are Conserved Across Rivers” submitted to PNAS (Stegen et al., 2023). The study aims to understand large-scale spatial structure of the dissolved organic matter (DOM) thermodynamic traits and inter-trait relationships by investigating (1) river water and sediments collected along 97 rivers spanning 3 continents and (2) coastal sediment collected from fresh and saline locations in Pacific and Gulf/Atlantic rivers. Sediment extracts and water samples were analyzed using ultrahigh resolution Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS). This dataset is comprised of three folders (1) Coastal, (2) WHONDR_S19S, and (3) Data_Dictionaries. Coastal contains (1) a subfolder with processed FTICR-MS data as csv files and sample collection metadata, (2) a subfolder with R scripts used to process the data and create associated figures, (3) a subfolder with the raw, unprocessed FTICR-MS data as .xml files, and (4) a readme file with more information about the dataset and instructions for using Formularity (https://omics.pnl.gov/software/formularity). WHONDRS_S19S contains (1) a csv file with processed FTICR data, (2) a csv with sample collection metadata, (3) a csv with sample geospatial data, (4) a csv with simulated lambda model outputs, (5) a subfolder with R scripts used to process the data and create associated figures, and (6) a readme file with more information regarding WHONDRS raw FTICR data and processing scripts. Data_Dictionaries contains data dictionaries for each csv file in the data package. The 97 global river corridors were part of a WHONDRS (https://whondrs.pnnl.gov) study. The raw, unprocessed FTICR-MS data with additional data can be found at doi:10.15485/1729719 for sediments and doi:10.15485/1603775 for water. This data package contains the processed data used in the associated manuscript. The coastal data has not been previously published, and this data package contains both the raw and processed data. Version 3 of this data package published February 2023 includes updates to the title of the manuscript, additional data and data dictionary and updated scripts linked to new analysis.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1988

Modal identification using single-mode projection filters and comparison with ERA and MLE results

The Single-Mode Projection Filter (SPF) is a newly developed algorithm for eigensystem parameter identification from both analytical results and test data. The SPF is formulated with a single mode only and practical for parallel processing implementation. Explicit formulations of SPF are derived for the multi-input multi-output (MIMO) system by using the orthogonal matrices of the controllability and observability matrices in the general sense. The modal parameters of SPF are initially obtained from an analytical model in modal space. The experimental data are then processed through SPF to update its modal parameters and to minimize a cost function defined by the norm of an error matrix. The updated modal parameters represent the characteristics of the test data. A two-dimensional global minimum optimization algorithm is developed and applied for the filter update by using the interval analysis method. The SPF is developed based on a single-mode subsystem and identifies only one modal frequency and one modal damping within a specified region. For an n-modes structure, n SPF can be implemented for parallel processing to reduce the computational burden. The SPF is applied to analyze the simulated data for the MAST beam structure. The estimated modal parameters are comparable to those from the Eigensystem Realization Algorithm (ERA) and repeated modal frequencies are identified. The modal analysis of the Spacecraft Control Laboratory Experiment (SCOLE) data is also performed by using the ERA and the Maximum Likelihood Estimate (MLE). The result shows that the first five modal frequencies are very close from ERA and MLE. However, there are slight disparities in the damping rates and the computational burdens are quite different among these two algorithms.

Huang, Jen-Kuang↗
NASA NTRS2003

Testing the Equivalence Principle in an Einstein Elevator: Detector Dynamics and Gravity Perturbations

We discuss specific, recent advances in the analysis of an experiment to test the Equivalence Principle (EP) in free fall. A differential accelerometer detector with two proof masses of different materials free falls inside an evacuated capsule previously released from a stratospheric balloon. The detector spins slowly about its horizontal axis during the fall. An EP violation signal (if present) will manifest itself at the rotational frequency of the detector. The detector operates in a quiet environment as it slowly moves with respect to the co-moving capsule. There are, however, gravitational and dynamical noise contributions that need to be evaluated in order to define key requirements for this experiment. Specifically, higher-order mass moments of the capsule contribute errors to the differential acceleration output with components at the spin frequency which need to be minimized. The dynamics of the free falling detector (in its present design) has been simulated in order to estimate the tolerable errors at release which, in turn, define the release mechanism requirements. Moreover, the study of the higher-order mass moments for a worst-case position of the detector package relative to the cryostat has led to the definition of requirements on the shape and size of the proof masses.

Hubbard, Dorthy↗
NASA NTRS2007

F-15B QuietSpike(TradeMark) Aeroservoelastic Flight Test Data Analysis

System identification or mathematical modelling is utilised in the aerospace community for the development of simulation models for robust control law design. These models are often described as linear, time-invariant processes and assumed to be uniform throughout the flight envelope. Nevertheless, it is well known that the underlying process is inherently nonlinear. The reason for utilising a linear approach has been due to the lack of a proper set of tools for the identification of nonlinear systems. Over the past several decades the controls and biomedical communities have made great advances in developing tools for the identification of nonlinear systems. These approaches are robust and readily applicable to aerospace systems. In this paper, we show the application of one such nonlinear system identification technique, structure detection, for the analysis of F-15B QuietSpike(TradeMark) aeroservoelastic flight test data. Structure detection is concerned with the selection of a subset of candidate terms that best describe the observed output. This is a necessary procedure to compute an efficient system description which may afford greater insight into the functionality of the system or a simpler controller design. Structure computation as a tool for black-box modelling may be of critical importance for the development of robust, parsimonious models for the flight-test community. Moreover, this approach may lead to efficient strategies for rapid envelope expansion which may save significant development time and costs. The objectives of this study are to demonstrate via analysis of F-15B QuietSpike(TradeMark) aeroservoelastic flight test data for several flight conditions (Mach number) that (i) linear models are inefficient for modelling aeroservoelastic data, (ii) nonlinear identification provides a parsimonious model description whilst providing a high percent fit for cross-validated data and (iii) the model structure and parameters vary as the flight condition is altered.

Kukreja, Sunil L.↗
NASA NTRS2015

Spatial Variability of Trace Gases During DISCOVER-AQ: Planning for Geostationary Observations of Atmospheric Composition

Results from an in-depth analysis of trace gas variability in MD indicated that the variability in this region was large enough to be observable by a TEMPO-like instrument. The variability observed in MD is relatively similar to the other three campaigns with a few exceptions: CO variability in CA was much higher than in the other regions; HCHO variability in CA and CO was much lower; MD showed the lowest variability in NO2All model simulations do a reasonable job simulating O3 variability. For CO, the CACO simulations largely under over estimate the variability in the observations. The variability in HCHO is underestimated for every campaign. NO2 variability is slightly overestimated in MD, more so in CO. The TX simulation underestimates the variability in each trace gas. This is most likely due to missing emissions sources (C. Loughner, manuscript in preparation).Future Work: Where reasonable, we will use these model outputs to further explore the resolvability from space of these key trace gases using analyses of tropospheric column amounts relative to satellite precision requirements, similar to Follette-Cook et al. (2015).

trace gas↗
NASA NTRSDate not supplied

An Interactive MATLAB Program for Fitting Transfer Functions to Frequency Responses

A computer program called FRFit (Frequency Response Fitting) is described for fitting single-input single-output transfer function models to empirical frequency response data. The program is interactive in that the user specifies ``elementary factors'' (gain, delay, pure differentiators and integrators, and first- and second-order zeros and poles) by entering numerical values or moving sliders in a graphical user interface. A nonlinear optimization can then be performed to obtain maximum likelihood estimates of transfer function parameters and uncertainties to provide feedback on the modeling and refine estimates. Several examples are discussed, including the identification of aircraft pitch dynamics from simulation data and data reported in the literature, approximating Theodorsen's function of unsteady aerodynamics, and obtaining a reduced-order model of a computational fluid dynamics code describing the unsteady aerodynamics around an aeroelastic wing. The program has some usefulness as a teaching aid, and can be applied to model structure determination, reduced-order modeling, preliminary analysis, and simple system identification problems. The program was written in MATLAB and is planned for public release through the NASA Software Catalog.

System identification↗
DOE OSTI2016

Environmental Characterization Scans Using Sodium-Iodide Gamma Spectroscopy to Determine Ra-226 Concentration in Surface Soil - 16326

The recent natural gas drilling boom has brought heightened attention to the issue of TENORM waste management. Radioactive residues from the extraction, treatment, and purification of materials containing uranium and thorium, or their progeny, pose a public health hazard primarily in the form of ground water contamination. Hydraulic fracturing can liberate formation waters from shale; the high salinity of which can cause radium to move with the water to the surface. The disposal and handling of this water has been a subject of public controversy. For soil remediation projects, radium-226 is generally one of the radionuclides of concern due to its alpha emissions. However performing a survey scan for Ra-226 is difficult because its photon emissions are low energy and have minimal probability of emission per decay (186.2 keV, 3.3% abundance). Additionally resolving its emission from U-235 (185.7 keV, 57% abundance) with field instrumentation is not feasible. Therefore detection of Ra-226 relies on the presence and photon emissions of its progeny. This paper outlines methodology for calculating the minimum detectable concentration of Ra-226 in characterization scans using a sodium-iodide detector array coupled with a multi-channel analyzer. Using progeny ingrowth as a detection proxy requires knowing decay time in order to determine progeny emission rates. Knowledge of site history, process, and prior survey results can be key in determining an accurate decay time. However in the absence of such information in-situ measurements can be used to estimate decay time by comparing observed progeny emission ratios against established values. These values are determined using Oak Ridge National Laboratory's ORIGEN-ARP program to generate photon energy spectra for Ra-226 and its progeny at specified time intervals. Once decay time is established, the detector's minimum detectable concentration can be calculated by synthesizing the expected emission spectra and Los Alamos National Laboratory's Monte Carlo N-Particle transport code to simulate detector response. Traditionally, one of the steps in calculating a system's minimum detectable concentration is to place an appropriate button source under the detector, record the response, and then move the source by a specified increment. This process is repeated until the detector's spatial efficiency has been mapped. Due to the impracticality of performing this with Ra-226 button sources of varying decay times, Monte Carlo N-Particle transport code was utilized to determine the detector's spatial efficiency. A dime sized source with the established emission spectra was placed into the center of a rectangle of soil. This rectangle, including the source, was used as a lattice element to fill the region of space below the detector. The simulation was designed to track the energy distribution of pulses created in the sodium iodide crystal using the F8 tally, and was setup to isolate each source; thereby correlating detector response to source position. The simulation output required post-processing to determine conversion factors from counts per second to disintegrations per second. Post-processing accounted for subtraction of the Compton continuum, and average the data over the scanning pattern to obtain the expected spatial sensitivity during a characterization study. This was performed for the following photo-peaks: 609 keV (Bi-214), and 1764 keV (Bi-214). For calculation of the scan minimum detectable concentration, the surveyor efficiency coefficient from NUREG 1507 was supplanted in favor of a scan speed deviation factor. This was used because the detector utilizes computerized data collection (spatial and spectral). This report finds that, for two shielded 10.16 cm x 10.16 cm x 40.64 cm sodium-iodide crystals with integrated multi-channel analyzers, the minimum detectable concentration in surface soils is less than 1 pCi per gram Ra-226 for a 1 second survey measurement. These results assume a volumetric source with a maximum depth of 15.24 cm and a minimum decay time of 21 days. Verification of the system's capabilities was performed by sending soil samples from scanned areas, and locations where in-situ measurements were performed, to an independent laboratory for analysis. The laboratory results demonstrated that scanned area was principally comprised of soil containing less than 1 pCi per gram Ra-226. This borders the system MDC and it underestimated the concentration of radium in the soil. The probable causes of this is escape of radon gas from the surface soil. Radon is the daughter of radium, and has a high mobility in soil. If the gas escapes, then the progeny used as the detection proxy will be present in lower than expected concentrations. Radionuclide identification using high volume sodium-iodide crystals in low level characterization scans is demonstrated to be possible. Further system refinement will focus on developing a methodology to estimate a radon escape factor for site specific use. Accounting for radon escape in the top 15.24 cm of soil will allow the system to better estimate Ra-226 concentrations in soil. This is a significant advancement over the gross counting statistics used in most studies to satisfy MARSSIM survey requirements, because it will increase stakeholder confidence. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2017

A Computational Simulation of Using Tungsten Gratings in Near-Field Thermophotovoltaic Devices

Near-field thermophotovoltaic (NFTPV) devices have received much attention lately as an alternative energy harvesting system, whereby a heated emitter exchanges super-Planckian thermal radiation with a photovoltaic (PV) cell to generate electricity. This work describes the use of a grating structure to enhance the power throughput of NFTPV devices, while increasing the energy conversion efficiency by ensuring that a large portion of the radiation entering the PV cell is above the band gap. The device contains a high-temperature tungsten grating that radiates photons to a room-temperature In 0.18 Ga 0.82 Sb PV cell through a vacuum gap of several tens of nanometers. Scattering theory is used along with the rigorous coupled-wave analysis (RCWA) to calculate the radiation energy exchange between the grating emitter and the TPV cell. A parametric study is performed by varying the grating depth, period, and ridge width in the range that can be fabricated using available fabrication technologies. It is found that the power output can be increased by 40% while improving the efficiency from 29.9% to 32.0% with a selected grating emitter as compared to the case of a flat tungsten emitter. Reasons for the enhancement are found to be due to the enhanced energy transmission coefficient close to the band gap. Furthermore, this work shows a possible way of improving NFTPV and sheds light on how grating structures interact with thermal radiation at the nanoscale.

42 ENGINEERING↗