Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reference Values”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Linking resource availability to pantropical forest canopy resistance and resilience to cyclone disturbance

Statement of purpose: Tropical cyclones are intensifying and occurring at higher latitudes in recent decades, but the mechanisms underpinning the resistance (ability to withstand disturbance-induced change) and resilience (pace of return to pre-disturbance reference values) of tropical forests to cyclones remains largely unexplored at the pantropical scale. We conducted a meta-analysis to investigate the role of soil resource availability (i.e., total soil phosphorus concentration) in mediating site-level forest canopy resistance and resilience to cyclones pan-tropically. We evaluated cyclone-induced and post-cyclone litterfall mass (g/m2/day), phosphorus (P) and nitrogen (N) fluxes (mg/m2/day), as well as concentrations (mg/g) across 73 case studies in Australia, Guadeloupe, Hawaii, Mexico, Puerto Rico, and Taiwan. The dataset zip file includes three data and two metadata files: - The compiled Litterfall Mass Flux data from tropical forests across the globe prior to and after varying tropical cyclone disturbances are provided in Litterfall_Mass.csv. This data file also includes site location, geographical characteristics, elevation, soil phosphorus concentration, geology, and several variables related to each tropical cyclone disturbance. - The compiled Litterfall Nitrogen and Phosphorus Flux data from tropical forests across the globe prior to and after varying tropical cyclone disturbances are provided in Litterfall_Nutrients.csv. This data file also includes site location, geographical characteristics, elevation, soil phosphorus concentration, geology, and several variables related to each tropical cyclone disturbance. - Tropical cyclone track data compiled from HURDAT2 and IBTrACS databases and used as input in the HURRECON model (https://github.com/hurrecon-model/HurreconR) to generate wind data is provided in hurdat2-1851-2019-052520.txt. - The metadata file (Metadata_Meta-analysis_Litterfall-Mass.pdf) has the complete information on each variable included in the Litterfall_Mass.csv dataset, the data sources, and data processing information. - The metadata file (Metadata_Meta-analysis_Litterfall-Nutrients.pdf) has the complete information on each variable included in the Litterfall_Nutrients.csv dataset, the data sources, and data processing information.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Temperature Effects on Droplet Oscillation Decay with Application to Fuel Property Measurement

Observation of oscillation decay in droplets has been shown to be an effective approach in determining physical properties such as viscosity and surface tension for emerging biofuels using μl quantities. Herein this work extends the approach to higher temperatures relevant to fuel injection conditions for internal combustion engines. Experiments are conducted that use high-resolution strobed imaging of moving, heated, μm-sized, fuel droplets to capture shape oscillation decay through image processing and analysis. Two fuels are investigated, iso-butanol and a primary reference fuel (PRF 84), which is a mixture of iso-octane and n-heptane. A piezoelectric droplet generator is used to generate a continuous train of single droplets, which are given an initial perturbation and subsequently undergo damped oscillations, captured using strobed imaging. Surface tension and viscosity calculated at four different temperatures ranging from 30°C-55°C using the frequency and decay time associated with the fundamental mode are found to be within ~ 10% of reference values obtained from the literature. Complementary numerical simulations are performed that utilize a volume-of-fluid approach to track transient droplet oscillation phenomena along with heat and mass transfer in a 2D axisymmetric domain. Simulations, where droplet size and temperature can be independently varied, capture an expected decrease in oscillation frequency and increase in decay time, with increase in fuel temperature. The simulations are further used to study the relative contribution to deviations in surface tension and viscosity predictions due to heat and mass transfer from the droplet, as well as viscous effects violating the inviscid flow assumption in the droplet oscillation theory. Mass loss effects are found to be negligible. Temperature change due to heat transfer is found to have the next highest sensitivity, particularly for the more volatile PRF 84, which has a lower viscosity and associated Ohnesorge number. For isobutanol, which has a higher viscosity and Ohnesorge number, viscous effects contribute the most to deviation in fuel property predictions.

09 BIOMASS FUELS↗
DOE OSTI2018

Subsoil Potassium in Central Iowa Soils: Status and Future Challenges

A reevaluation of plant-available K within central Iowa soils is needed, especially if crop residues are harvested for animal feed or cellulosic feedstock. This study evaluated plant-available K concentrations to depths of 120 cm (n = 2433) within soil profiles at 10 experimental sites. Data from 0- to 5-, 5- to 15-, 0- to 15-, 15- to 30-, 30- to 60-, 60- to 90-, and 90- to 120-cm depths were compiled and evaluated. Based on over 400 samples collected from each depth increment below 15 cm, a reference value of 100 mg kg –1 appears common for central Iowa subsoils. Sixty-four percent of all samples had soil-test K levels below the 120 mg kg –1 threshold for very low nutrient status for crop production in Iowa. Only 28% of surface soil samples (<15-cm depth) were considered to have at least optimum nutrient status with concentrations greater than 160 mg kg –1 . As K is not an environmentally sensitive nutrient, its management in central Iowa soils has generally been neglected in recent years, particularly when compared with N and P management initiatives. Building awareness of soil K profile concentrations, rather than focusing only on surface concentrations, will be necessary to meet nutrient requirements for 21st century agricultural production.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Risks to Ecological Receptors Posed by Contaminants of Potential Concern in the Lower Three Runs Cooling Ponds and Canals

The upper portion of Lower Three Runs includes several ponds, reservoirs, and canals that were formerly used as a cooling system for nuclear production reactors. This area was divided into nine exposure areas (EAs) for the assessment of environmental contamination resulting from past reactor operations and other industrial processes. A tiered screening process identified several contaminants of potential concern including aluminum, cyanide, lead, manganese, mercury, DDD, DDE, and DDT. Risks posed by these contaminants to ecological receptors (river otter, belted kingfisher, raccoon, and blue heron) were assessed using contaminant exposure models that estimated contaminant intake resulting from ingestion of food, water, and sediment/ soil and compared these intakes with toxicity reference values (TRVs). The contaminant exposure models showed that the TRVs were not exceeded in the otter model, exceeded by aluminum in EA 7 (Pond 2 and associated canals) in the raccoon model, and exceeded by mercury in EAs 2, 3 (Pond B), 6 (Par Pond), and 8 (Ponds 4 and 5 and Canal to Pond C) in both the kingfisher and blue heron models. Hazard quotients (total exposure dose divided by the TRV) were 2.8 for aluminum and 1.7- 3.6 for mercury. The primary route of exposure for aluminum was the ingestion of soil, and the primary route of exposure for mercury was the ingestion of mercury contaminated fish. Elevated levels of mercury in fish were at least partly the result of the aerial deposition of mercury onto Lower Three Runs and its watershed. The atmospheric deposition of mercury creates pervasive contamination in fish throughout the Savannah River basin. Another possible source of mercury was the discharge of mercury contaminated Savannah River water into the Lower Three Runs cooling ponds and canals during previous years of reactor operation. This contamination originated from industries located upstream of the SRS. The aluminum exceedance for the raccoon was likely the result of naturally high aluminum levels in SRS soils rather than SRS operations. Aluminum exceedances have previously been observed in relatively undisturbed background locations as well as areas affected by SRS operations. Aluminum exceedances are more likely with the raccoon than the other receptors because it consumes more soil as a result of its feeding habits. Sensitivity analysis showed that model uncertainty can be reduced by adequate sampling of key variables (e.g., fish and sediments). Although sediment samples were collected from all EAs, fish samples were not collected from three EAs and some analytes (pesticides and cyanide) were not measured in fish. Water-to-fish concentration ratios were used to estimate contaminant levels in fish when direct measurements from fish were unavailable; however, such estimates are potentially less accurate than direct measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Validation of Spent Fuel Dose Rate Calculations Using NDA Measurement Campaign Data

The predicted radiation gamma dose rate is a key parameter to determine if a spent nuclear fuel assembly is self-protected from misuse as well as an indicator for worker/inspector safety. Recent studies showed that using state-of-the-art techniques for dose rate calculations, including the one adopted for this study, the numerical calculations of the gamma dose rates are three times lower than those from the most commonly-cited reference values for spent LWR fuel used by domestic and international agencies. With only a few publically available experimental data points worldwide for this type of benchmark study, the gamma flux measurements from the Spent Fuel Non Destructive Assay (NDA) measurement campaign can be used to help validate the physics approach used in calculating gamma dose rates. Specifically, this data would validate the gamma flux calculation from spent fuel, which is a key component of these multi-physics techniques for dose rate assessments; dose rate measurements were not part of the NDA campaign. These data are primarily the passive gamma measurements of the spent nuclear fuel assemblies from the Swedish Nuclear Fuel and Waste Management Company (SKB).

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

DECOVALEX-2019 (Task E Final Report)

The DECOVALEX Project is an on-going international research collaboration, established in 1992, to advance the understanding and modeling of coupled Thermal (T), Hydrological (H), Mechanical (M) and Chemical (C) processes in geological in geological systems. DECOVALEX was initially motivated by the recognition that prediction of these coupled effects is an essential part of the performance and safety assessment of geologic disposal systems for radioactive waste and spent nuclear fuel. Later it was realized that these processes also play a critical role in other subsurface engineering activities, such as subsurface CO 2 storage, enhanced geothermal systems, and unconventional oil and gas production through hydraulic fracturing. Research teams from many countries (e.g., Canada, China, Czech Republic, Finland, France, Germany, Japan, Republic of Korea, Spain, Sweden, Switzerland, Taiwan, United Kingdom, and the United States) various institutions have participated in the DECOVALEX Project over the years, providing a wide range of perspectives and solutions to these complex problems. These institutions represent radioactive waste management organizations, national research institutes, regulatory agencies, universities, as well as industry and consulting groups. This document is the final report of Task E which was proposed and coordinated by Andra, the National Radioactive Waste Management Agency in France, presenting the technical definitions of the problems studied, approaches applied, achievements made and outstanding issues for future research. The purpose of Task E of the DECOVALEX-2019 project is to investigate upscaling THM modelling from small-scale experiments (some cubic meters) to full-scale experiments (some ten cubic meters) and finally to the scale of the waste repository (cubic kilometers). To achieve this aim, the data of two in-situ heating experiments performed by Andra (the French National Radioactive Waste Management Agency) in the Meuse/Haute-Marne Underground Research Laboratory (MHM URL) have formed the basis for the understanding of the THM behavior of the COx at different scales. The first experiment provided the reference values of the THM parameters by means of a calibration exercise and they were used for a blind prediction and an interpretative analysis of the second one.

58 GEOSCIENCES↗
DOE OSTI2020

Transfer Factors for the FRMAC Assessment Manual and Turbo FRMAC to Improve Radiological Dose Assessment

The Turbo FRMAC analysis tool is used to perform complex calculations to quickly evaluate radiological consequences and aid in decision making during an emergency response by assessing impacts to the public, workers, and the food supply. Turbo FRMAC calculations are based on methods established by the Federal Radiological Monitoring and Assessment Center (FRMAC). To be able to assess impacts, input data called transfer factors that describe radionuclide uptake by local plants and animals are required. During the code application exercises in late 2016, identifying, finding, and validating nonstandard transfer factors proved to be time consuming and diverted the teams’ activity away from other critical tasks. As a result, a task was undertaken to dramatically expand the list of available transfer factors (food and non-food) and incorporate these factors into Turbo FRMAC. This will ultimately result in improved efficiency of the assessment team to perform calculations during times when the FRMAC is activated and provide more defensible, vetted data from which to calculate results. As a result of these and subsequent exercises, transfer factors were needed for the following items: bell peppers, Christmas tree, deer, flowers, fresh cucumbers, tomatoes, grapefruit, lichens, mushrooms, oranges, snap beans, squash, strawberries, sugarcane for sugar and seed, sweet corn, tea, tobacco, tree bark, and watermelon. To expand the applicability of the information tables, generic transfer data were also provided for common categories and recommendations were made for expanding the list of chemicals based on chemical similarity. The data presented in this report were compiled from recent literature with most of the data encompassing the period from 2000 to 2018. The following radionuclides were targeted during the literature search: elements associated with reactor accidents or nuclear detonations (Sr 89/90, Cs 134/137, Ce 141/144, Ru 103/106, I 129/131/133, Pu 238/239, Am 241, Zr 95, and Nb 95) and elements associated with industrial accidents or dirty bombs (Ir 192 and Co 60). Data from other elements were evaluated if they were identified during the literature search. In addition, reports were evaluated that were recommended by the research consultants. For each plant or animal transfer factor, the goal was to determine the geometric mean, the geometric standard deviation, the minimum, the maximum, and the number of measurements used. If only one measurement was available that was presented as the mean. Concentration ratios with large geometric standard deviations (GSDs) were generally the result of a paucity of measurements or a few disparate measurements. Measurement disparity was observed for data from different soil types. The geometric mean (GM) is a good reference value for planning and responses purposes, but the location-specific concentrations are unlikely to be similar to the model results. Recommendations for further work include: segregating the data to reflect the influence of soil type, developing approaches to incorporate animal data based on aggregated transfer measurements, and including data for foliar deposition on plants.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2024

Initial Results of the Nuclear Accident Dosimetry Intercomparison at the Armed Forces Radiobiology Research Institute's TRIGA Reactor

This report presents the initial results for IER-602, the international nuclear accident dosimetry intercomparison performed at the Armed Forces Radiobiology Research Institute (AFRRI). Ten groups participated (six DOE laboratories, three international groups, and one DOD group) in a week-long exercise in late June 2024. Two additional groups observed the exercise (one DOE laboratory, one university). Three days of irradiations were performed, and participants reported results within 24 hours of the irradiations. Results were compared to the reference values and performance criteria established by IER-484 AFRRI ER1 Dosimetry Characterization and ANSI/HPS N13.3-2013, respectively. Passing performance varied from 17% to 100%, with an average of 67%.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2026

Machine Learning Vacancy Formation Energy in Nickel-Based Superalloys

Thermal vacancies play a critical role in high-temperature Ni-based superalloys and influence various properties such as creep resistance, oxidation, etc. This study systematically investigates the impact of commonly used transition metals (Cr, Co, Fe), refractory metals (Nb, Ta, Mo, W) and other elements (Al, Cu, Ti, Mn) on the thermodynamic stability of 36 binary, 20 ternary, 11 quaternary, 9 quinary, and 3 senary FCC Ni-based alloys covering various elemental combinations. Density functional theory-based studies on Ni-X binary alloys show that higher concentrations of Cr, Nb, Ta, Al, and Ti introduce significant lattice distortions and broaden the distribution of vacancy formation energies (standard deviation up to 0.15 eV). These elements partially donate electrons, reducing their self-consistent chemical potentials relative to single-element reference values and lowering vacancy formation energies, while Co, Fe, Mo, and W show lower charge localization. These trends extend from 3-6 element alloys, where Cr, Nb, and Ta-rich compositions have low-energy states (~0.5 eV) that increase vacancy concentrations. Finally, graph neural network models are developed to screen over 5000 virtual alloys. Eleven leading compositions are identified with mean vacancy formation energy higher than 1.75 eV and vacancy concentration ~2 orders of magnitude lower than pure Ni at 1000 K. These results provide valuable guidelines to achieve controlled defect engineering in structural alloys.

DFT↗
DOE OSTI2025

Anion exchange membrane test protocol validation

This study presents the validation of protocols for measuring ion exchange capacity (IEC) and alkaline stability of anion exchange membranes (AEMs) for low-temperature water electrolysis. While protocols are often tested within individual laboratories, their results across multiple laboratories with varying equipment, environmental conditions, and personnel qualification remain unverified. The validation involved Los Alamos National Laboratory (LANL), National Renewable Energy Laboratory (NREL), and University of Oregon (UO) using the same commercially available AEM to assess reproducibility and reliability of the protocols under diverse conditions. For the IEC protocol, results across laboratories were consistent within ±10% of the NMR-determined reference value. The alkaline stability protocol could pose greater challenges due to factors such as variations in sample collection timing, preservation methods, and analytical techniques, but consistent test results for percentage IEC loss were demonstrated across institutions. These results highlight the reliability and applicability of the protocols, emphasizing the importance of validation to ensure consistency in diverse research environments.

08 HYDROGEN↗
DOE OSTI2024

Techno-Economic Analysis of Cement Decarbonization Techniques: Oxygen Enrichment vs. Hydrogen Fuel

The Paris Agreement aims to limit global warming, and one of the most polluting sectors is heavy industry, where cement production is a significant contributor. This work briefly explores some alternatives, recycling, reducing clinker content, waste heat recovery, and carbon capture, discussing their advantages and drawbacks. Then, it examines the economic viability and benefits of increasing oxygen concentration in the primary burning air from 21 to 27 vol.%, which could improve clinker production by 7%, and the production of hydrogen through PEM electrolysis to make up 5% of the fuel thermal fraction, considering both in a cement plant producing 3000 tons of clinker per day. This analysis used reference values from Secil, an international company for cement and building materials, to determine the required scale of the oxygen and hydrogen production, respectively, and calculate the CAPEX of each approach. It is concluded that oxygen enrichment can provide substantial fuel savings for a relatively low cost despite a possible significant increase in NOx emissions. However, hydrogen production at this scale is not currently economically viable.

Domingues, Bruno C.↗
DOE OSTI2023

Electron Scattering from 1-Methyl-5-Nitroimidazole: Cross-Sections for Modeling Electron Transport through Potential Radiosensitizers

In this study, we present a complete set of electron scattering cross-sections from 1-Methyl-5-Nitroimidazole (1M5NI) molecules for impact energies ranging from 0.1 to 1000 eV. This information is relevant to evaluate the potential role of 1M5NI as a molecular radiosensitizers. The total electron scattering cross-sections (TCS) that we previously measured with a magnetically confined electron transmission apparatus were considered as the reference values for the present analysis. Elastic scattering cross-sections were calculated by means of two different schemes: The Schwinger multichannel (SMC) method for the lower energies (below 15 eV) and the independent atom model-based screening-corrected additivity rule with interferences (IAM-SCARI) for higher energies (above 15 eV). The latter was also applied to calculate the total ionization cross-sections, which were complemented with experimental values of the induced cationic fragmentation by electron impact. Double differential ionization cross-sections were measured with a reaction microscope multi-particle coincidence spectrometer. Using a momentum imaging spectrometer, direct measurements of the anion fragment yields and kinetic energies by the dissociative electron attachment are also presented. Cross-sections for the other inelastic channels were derived with a self-consistent procedure by sampling their values at a given energy to ensure that the sum of the cross-sections of all the scattering processes available at that energy coincides with the corresponding TCS. This cross-section data set is ready to be used for modelling electron-induced radiation damage at the molecular level to biologically relevant media containing 1M5NI as a potential radiosensitizer. Nonetheless, a proper evaluation of its radiosensitizing effects would require further radiobiological experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Representation of the QM Subsystem for Long-Range Electrostatic Interaction in Non-Periodic Ab Initio QM/MM Calculations

In QM/MM calculations, it is essential to handle electrostatic interactions between the QM and MM subsystems accurately and efficiently. To achieve maximal efficiency, it is convenient to adopt a hybrid scheme, where the QM electron density is used explicitly in the evaluation of short-range QM/MM electrostatic interactions, while a multipolar representation for the QM electron density is employed to account for the long-range QM/MM electrostatic interactions. In order to avoid energy discontinuity at the cutoffs, which separate the short- and long-range QM/MM electrostatic interactions, a switching function should be utilized to ensure a smooth potential energy surface. In this study, we benchmarked the accuracy of such hybrid embedding schemes for QM/MM electrostatic interactions using different multipolar representations, switching functions and cutoff distances. For test systems (neutral and anionic oxyluciferin in MM (aqueous and enzyme) environments), the best accuracy was acquired with a combination of QM electrostatic potential (ESP) charges and dipoles and two switching functions (long-range electrostatic corrections (LREC) and Switch) in the treatment of long-range QM/MM electrostatics. It allowed us to apply a 10Å distance cutoff and still obtain QM/MM electrostatics/polarization energies within 0.1 kcal/mol and time-dependent density functional theory (TDDFT)/MM vertical excitation energies within 10−3 eV from theoretical reference values.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

A Novel Machine Learning Algorithm for Cloud Detection Using AERI Measurement Data

Infrared hyperspectral remote sensing has been widely used in the field of meteorology. Many scientists have carried out research on inversion methods of meteorological elements such as thermodynamic profile, boundary layer height, cloud base height, etc. In this study, a method based on machine learning for cloud detection using ground-based infrared hyperspectral radiation data is proposed. The features of outliers, the cloudy and cloud-free data of Atmospheric Emitted Radiance Interferometer (AERI) radiation are extracted. The “reference values” of cloudy and cloud-free are determined based on the observation data of Vaisala CL31 ceilometer within the time range of 8 min before the corresponding time of AERI. A support vector machine (SVM) algorithm is used for training. The dataset comes from the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site and North Slope Alaska (NSA) site from 2015 to 2017, and the ARM West Antarctic Radiation Experiment (AWARE) site in 2016 is also analyzed. The instruments used in this paper include AERI, ceilometer, etc. The experimental results reveal that the agreement of cloud detection results between the proposed algorithm and ceilometer is about 93% at each site. However, for high clouds or optically thin clouds, the agreement will decrease.

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

Bias corrections of GOSAT SWIR XCO 2 and XCH 4 with TCCON data and their evaluation using aircraft measurement data

We describe a method for removing systematic biases of column-averaged dry air mole fractions of CO 2 (XCO 2 ) and CH 4 (XCH 4 ) derived from short-wavelength infrared (SWIR) spectra of the Greenhouse gases Observing SATellite (GOSAT). We conduct correlation analyses between the GOSAT biases and simultaneously retrieved auxiliary parameters. We use these correlations to bias correct the GOSAT data, removing these spurious correlations. Data from the Total Carbon Column Observing Network (TCCON) were used as reference values for this regression analysis. To evaluate the effectiveness of this correction method, the uncorrected/corrected GOSAT data were compared to independent XCO 2 and XCH 4 data derived from aircraft measurements taken for the Comprehensive Observation Network for TRace gases by AIrLiner (CONTRAIL) project, the National Oceanic and Atmospheric Administration (NOAA), the US Department of Energy (DOE), the National Institute for Environmental Studies (NIES), the Japan Meteorological Agency (JMA), the HIAPER Pole-to-Pole observations (HIPPO) program, and the GOSAT validation aircraft observation campaign over Japan. These comparisons demonstrate that the empirically derived bias correction improves the agreement between GOSAT XCO 2 /XCH 4 and the aircraft data. Finally, we present spatial distributions and temporal variations of the derived GOSAT biases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Kinetic Monte Carlo Simulation of Hydrogen Diffusion in Tungsten

Nuclear fusion is considered as the most promising sustainable energy source. If achieved, it could provide electricity to the mankind for millions of years. During the reaction, the fusion ash is continuously removed from the magnetic gap, the divertor region, and the divertor facing component is subjected to intense plasma bombardment. It is proposed to use tungsten as the divertor facing material, because it has low physical sputtering yield and superior thermomechanical properties. However, the high mobility of hydrogen in tungsten and the high bombardment flux bring the hydrogen isotopes deep into the tungsten tile. The presence of hydrogen in tungsten tends to weaken the gram boundaries and produce surface blistering. The ejection of tungsten particulates into the plasma core will cool the reaction and cease the fusion reactions. Therefore, it is important to have a detailed understanding of the hydrogen diffusion in tungsten bulk. The KMC method is a stochastic computer simulation to predict the time evolution of a certain system. Typically, at a given time, the system will evolve through a few possible paths with known probabilities. In each step, the KMC method analyzes the entire system and randomly determines the system evolving events based on all transition rates of all possible events and advance the physical time. This process is repeated till the simulation of the entire physical time is completed. In previous work, the authors developed a kinetic Monte Carlo (KMC) algorithm and program to simulate hydrogen diffusion process on tungsten reconstructed surface. In this article, this KMC method is extended to model the hydrogen diffusion in tungsten bulk. The diffusion coefficients are calculated and fitted to the Arrhenius equation. The result matches the experiment measurement very well. KMC diffusion simulation needs the locations of the hydrogen trapping sites, diffusion paths and the corresponding diffusion energy barriers. The hydrogen interstitial sites in tungsten crystal have been thoroughly studied by density function theory (DFT) and molecular dynamics (MD). Johnson and Carter used DFT to determine the preferred location of the hydrogen interstitial sites m tungsten crystal, which are the tetrahedral sites (T-sites) in body-centered-cubic (bcc) lattice surfaces. All T-sites are located on the surfaces of a lattice cell. There are four T-sites per surface and twelve T-sites per lattice cell. The hydrogen atoms diffuse by traveling from one T-site to one of its nearest four T-site neighbors. Thus, each T-site is surrounded by two neighbors on the same surface, and one neighbor on each of the two adjacent surfaces perpendicular to the surface holding the center T-site. The diffusion barrier between two neighboring T-sites was calculated to be 0.38 eV. Knowing the diffusion path, energy barrier, and the crystal structure, the KMC simulation is readily to be performed. Based on other researchers' DFT calculations, the hydrogen stable interstitial site in tungsten bulk is determined as the tetrahedral site on the surface of the bee lattice. Each T-site has 4 neighbors and the energy barrier between two adjacent T-sites are estimated to be 0.38 eV. Based on those facts and the KMC theory, a MATLAB program is developed to model the hydrogen diffusion in tungsten bulk. In the MATLAB program, substrate temperature, hydrogen concentration, physical time to be simulated, and the number of cycles are all customizable. To compare the KMC determined diffusion coefficients with experiment results in the literature, the simulations are performed for temperature 300 K - 2500 K and low hydrogen concentration (0.1%), mimicking the experiment condition. The calculated diffusion coefficients match the reference values measured by the experiment under low hydrogen concentration condition. The results are also fitted to the Arrhenius equation. The pre-exponential factor and the average activation energy in Arrhenius equation are obtained. The good agreement between simulation and experiment proves that the developed KMC algorithm truly reflect the well-established theory of interstitial hydrogens in tungsten. The developed program could easily be used to model the diffusion process and predict the diffusion coefficients under other concentrations and temperatures. This work is part of the ongoing development of hydrogen diffusion in tungsten integrating surface diffusion (previous work), bulk diffusion (this work), and adsorption/desorption on surface (future work)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Adaptive Sparse-Grid Stochastic Collocation Uncertainty Quantification Convergence for Multigroup Diffusion

Advanced methods in uncertainty quantification for numerical models in computational physics continue to gain acceptance in nuclear modeling. Previously, the efficiency of sparse-grid stochastic collocation in comparison with Monte Carlo for uncertainty quantification through convergence studies was demonstrated. An adaptive method for anisotropic sparse grid collocation and its convergence compared to previous efforts is considered. The physical system in consideration is a two-dimensional quarter-core reactor, consisting of 5 materials distributed in 121 regions. The two-group neutron diffusion criticality approximation is solved. Vacuum boundaries are applied on the top and right, and reflecting boundaries on the bottom and left. The material cross sections, neutron multiplication factors, and diffusion coefficients are potential sources of uncertainty. Uniformly-distributed uncertainty within 5% of the reference values is introduced for material properties. The distribution of the uncertain parameters make up the uncertainty space Γ included in R{sup N}, where N is the number of uncertain parameters. The input parameter uncertainties are independently distributed. As expected, the convergence of Monte Carlo on the benchmark as a function of computational solves is nearly linear and lethargic compared to the other methods. The static isotropic index set, in which each input dimension is treated with equal importance, converges much more efficiently than Monte Carlo, but lacks the benefit of emphasizing some dimensions over others. The anisotropic case is most efficient, where first-order sensitivity information was used to determine anisotropic weights for polynomials to use in the polynomial chaos expansion. Interestingly, the adaptive algorithm converges slightly less efficiently than the anisotropic case, while still performing better than the isotropic case. Because the adaptive case has to determine the importance of each input dimension, it occasionally predicts inaccurately and wastes computation solves compared to the more ideal anisotropic case. However, in circumstances when an ideal anisotropy is not understood a priori, the adaptive algorithm may be optimally efficient in representing the original model. Given the potential loss of efficiency when using a poorly-chosen anisotropic weighting, the adaptive algorithm may be an effective choice for initial uncertainty quantification in problems with a low-dimensionality input space. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

A Study on Reconstruction of Intra-Pin Information in Pin-Wise Reactor Analysis

Along with nodal calculation, accurate pin-wise calculations with parallel computational platforms have been recently developed as alternative tools to nodal methods to maximize core performance and minimize uncertainties used in reactor safety analysis. Therefore, it is necessary to have a method for accurately reconstructing intra-pin information comparable to whole-core transport calculation while being significantly more computationally efficient. In this study, a method to reconstruct intra-pin information (intra-pin reconstruction) is proposed. This method is based on combining the local intra-pin form function with the global homogeneous distribution. The form functions are generated from a single pin calculation with reflective boundary conditions using Serpent 2, a Monte Carlo neutron transport code. The global homogeneous flux was reconstructed using the analytical solution of the diffusion equation due to its superior accuracy and speed compared to other methods. The study here is limited to a two-dimensional two-group problem of rectangular geometric fuel pin element. The pin-wise boundary constraints (BCs) and reference intra-pin power used for reconstruction and comparison were obtained from the Serpent 2 reference calculation. From our study, it is indicated that maximum relative error (MRE) between reconstructed and reference values is within 5%. In addition, a sensitivity study was also performed to assess the feasibility of the method when the proposed method is coupled with a pin-wise calculation. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗