Search NASA⌕ Search

SEARCH · Search NASA

Results for “Code Assessment”

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,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Detailed modeling of electron emission for transpiration cooling of hypersonic vehicles

Electron transpiration cooling (ETC) is a recently proposed approach to manage the high heating loads experienced at the sharp leading edges of hypersonic vehicles. Computational fluid dynamics (CFD) can be used to investigate the feasibility of ETC in a hypersonic environment. In this paper, a modeling approach is presented for ETC, which includes developing the boundary conditions for electron emission from the surface, accounting for the space-charge limit effects of the near-wall plasma sheath. The space-charge limit models are assessed using 1D direct-kinetic plasma sheath simulations, taking into account the thermionically emitted electrons from the surface. The simulations agree well with the space-charge limit theory proposed by Takamura et al. for emitted electrons with a finite temperature, especially at low values of wall bias, which validates the use of the theoretical model for the hypersonic CFD code. The CFD code with the analytical sheath models is then used for a test case typical of a leading edge radius in a hypersonic flight environment. The CFD results show that ETC can lower the surface temperature of sharp leading edges of hypersonic vehicles, especially at higher velocities, due to the increase in ionized species enabling higher electron heat extraction from the surface. Finally, the CFD results also show that space-charge limit effects can limit the ETC reduction of surface temperatures, in comparison to thermionic emission assuming no effects of the electric field within the sheath.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Impact of ejecta morphology and composition on the electromagnetic signatures of neutron star mergers

The electromagnetic transients accompanying compact binary mergers (γ-ray bursts, afterglows and ‘macronovae’) are crucial to pinpoint the sky location of gravitational wave sources. Macronovae are caused by the radioactivity from freshly synthesized heavy elements, e.g. from dynamic ejecta and various types of winds. We study macronova signatures by using multidimensional radiative transfer calculations. We then employ the radiative transfer code SUPERNU and state-of-the-art LTE opacities for a few representative elements from the wind and dynamical ejecta (Cr, Pd, Se, Te, Br, Zr, Sm, Ce, Nd, U) to calculate synthetic light curves and spectra for a range of ejecta morphologies. The radioactive power of the resulting macronova is calculated with the detailed input of decay products. We assess the detection prospects for our most complex models, based on the portion of viewing angles that are sufficiently bright, at different cosmological redshifts (z). The brighter emission from the wind is unobscured by the lanthanides (or actinides) in some of the models, permitting non-zero detection probabilities for redshifts up to z = 0.07. We also find that the nuclear mass model and the resulting radioactive heating rate are crucial for the detectability. While for the most pessimistic heating rate (from the finite range droplet model) no reasonable increase in the ejecta mass or velocity, or wind mass or velocity, can possibly make the light curves agree with the observed near-infrared excess after GRB130603B, a more optimistic heating rate (from the Duflo–Zuker model) leads to good agreement. We conclude that future reliable macronova observations would constrain nuclear heating rates, and consequently help constrain nuclear mass models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

A continuous symmetry breaking measure for finite clusters using Jensen-Shannon divergence

A quantitative measure of symmetry breaking is introduced that allows the quantification of which symmetries are most strongly broken due to the introduction of some kind of defect in a perfect structure. The method uses a statistical approach based on the Jensen-Shannon divergence. The measure is calculated by comparing the transformed atomic density function with its original. Software code is presented that carries the calculations out numerically using Monte Carlo methods. The behavior of this symmetry breaking measure is tested for various cases including finite size crystallites (where the surfaces break the crystallographic symmetry), atomic displacements from high symmetry positions, and collective motions of atoms due to rotations of rigid octahedra. Finally, the approach provides a powerful tool for assessing local symmetry breaking and offers new insights that can help researchers understand how different structural distortions affect different symmetry operations.

atomic & molecular structure↗
DOE OSTI2025

Dataset: "Widespread Drought-driven Declines in Streamflows and Water quality in the Upper Colorado River Basin (1998-2022)"

This data package contains the associated data and scripts for Nagamoto, E., Ombadi, M., Ciulla, F. et al. Widespread drought-driven declines in streamflows and water quality in the Upper Colorado River Basin during 1998-2022. Commun Earth Environ 7, 734 (2026). https://doi.org/10.1038/s43247-026-03890-5. This purpose of this study was to investigate the impact of the 21st century drought on water quantity and quality at catchments throughout the Upper Colorado River Basin (UCRB). We used stream flow, water temperature, specific conductance, air temperature, precipitation, and catchment attribute data for over 200 sites in the UCRB, collected from the National Water Information System using Basin3D (Varadharajan, 2023), GAGESII (Falcone, 2010), and the Google Earth Engine. We identified years of severe drought between 1998 and 2022 using the Standardized Precipitation Evaporation Index (SPEI), then calculated the relative change percentage of the stream flow, water temperature, and specific conductance from drought versus non-drought years. We used the attribute information from GAGESII to investigate what physical traits of catchments are associated streamflow vulnerability (greater relative change) or resilience to drought. We used land cover data from the National Land Cover Database (USGS, 2024) to assess any changes to physical attributes that may not be represented in the static attributes information in GAGESII. To increase data availability, we modeled stream temperature using methods from Willard, 2023. While the study period is water years 1998 to 2022, the raw water quantity and quality data extends to 1950 and the meteorological data extends to 1980. The data and code can be downloaded via the UCRB_drought.zip. Within the zip, the files are organized as follows: - INPUTS: Contains all input data used in UCRB_Drought_Workflow.ipynb - OUTPUTS: Contains all intermediate data created from UCRB_Drought_Workflow.ipynb as well as final products including the calculated Standardized Evapotranspiration Index (SPEI) - climatic_variables: The code used to collect meteorologic data from Google Earth Engine - feature_importance: The code used for the catchment attributes analysis - preprocessing: Code used in UCRB_Drought_Workflow_Preprocessing.ipynb - pyeto: Code used in UCRB_Drought_Workflow_Preprocessing.ipynb - calculations: Code used in UCRB_Drought_Workflow_Impacts.ipynb - plotting: Code used in UCRB_Drought_Workflow_Impacts.ipynb - README.md - UCRB_Drought_Workflow_Preprocessing.ipynb: The code used to prep raw data for the analysis - UCRB_Drought_Workflow_Impact.ipynb: The code which uses the prepped raw data for analysis, and plots all figures - requirements_ucrb-drought_v2.yml: The requirements file to create a virtual environment and Jupyter Lab kernel to run the code The INPUTS folder is organized into the following major directories and sub-directories. The "RDC_WT_SC_RAW" folder contains raw data for streamflow, water temperature, and specific conductance in a ".h5" file. The "NLCD_RAW" folder contains ".csv" files with annual land cover percentages for counties within the UCRB. The "MET_RAW" folder contains a ".csv" file with monthly meteorological data (air temperature and precipitation) for the sites in the UCRB which was obtained from code in the climatic_variables folder. The "GAGESII" folder contains ".csv" files with physical catchment attribute variables for catchments across the country. The "WT_LSTM_data" folder contains ".csv" files with calculated WT (Willard, 2023) and the associated RMSEs. The "Upper_Colorado_River_Basin_Boundary" folder contains geographic data including a shapefile for plotting in the UCRB_Drought_Workflow.ipynb. The "RESERVOIRS_RAW" folder contains ".csv" files for each reservoir in the UCRB with daily reservoir storage. There are also two files in the INPUTS folder that have combined reservoir storage data and reservoir metadata. The OUTPUTS folder is organized into the following major directories and sub-directories. The "RDC_WT_SC_data" folder contains a folder "Water_year" with the associated cleaned data, metadata, and data availability information in ".csv" files, a folder "Median_Relchange" with the relative change comparing drought to non-drought years in ".csv" files, and a folder "Peak95_Min5_Relchange" that has ".csv" files for the relative change in peak (95th %) and minimum (5th %) variables. The "NLCD_data" folder contains the difference in land cover from the beginning to end of the study period and the percentage of the county that is within UCRB bounds can be found in Nagamoto et al (2025)). The "MET_data" folder contains separated monthly air temperature and precipitation data and the calculated PET in ".csv" files. The "SPEI_data" folder contains ".csv" files with calculated SPEI values (one restricted to the study period and the other with information from the entire MET data period). The "Paper_Tables" folder contains two ".csv" files containing site information and data availability and information about the GAGESII trait aggregated categories. The base directory includes the file “flmd.csv” for a list and description of all files and the file “dd.csv” for data dictionaries. Scripts for preprocessing, analysis, and figure generation are located in the associated GitHub repository found at [https://github.com/iNAIADS/drought-impacts/tree/develop/UCRB-drought]. UPDATE 1: Title and code file updated to match submitted manuscript 10-15-2025. UPDATE 2: Code and data files updated to match revised manuscript 3-4-2026. UPDATE 3: Code and data files updated to match revised manuscript 6-7-2026. ** NOTE: DD and FLMD have not been updated yet. UPDATE 4: Added associated Manuscript information and DD and FLMD have been updated. To cite this code, please use the following BibTeX: @misc{nagamoto2025drought, author = {Emily Nagamoto and Fabio Ciulla and Mohammad Ombadi and Jared Willard and Rosemary Carroll and Charuleka Varadharajan}, title = {Dataset: "Widespread Drought-driven Declines in Streamflows and Water quality in the Upper Colorado River Basin (1998-2022)"}, year = {2025}, doi = {10.15485/2551894}, publisher = {ESS-DIVE Repository}, url = {https://data.ess-dive.lbl.gov/datasets/doi:10.15485/2551894} }

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Response Technical Tools User Guide [Manual]

Developed by the U.S. Nuclear Regulatory Commission for assessment of reactor accidents, Response Technical Manual (RTM) is a paper-based report which contains simple methods for estimating possible accident scenarios and relevant consequences for different kinds of radiological events. Based on RTM, a software called Response Technical Tools (RTT) was developed by Sandia National Laboratories to convert the paper manual into an automated and easy-to-use code. In particular, the RTT focuses on the nuclear power plant severe accidents and is informed by state-of-the-art analyses and software programs, such as MELCOR and MAAP. RTT evaluations can be used to track and predict, at a very coarse level, the progression of a severe accident in nuclear power plant. The RTT allows a user to track the progress of an accident in a nuclear power plant from the point of initiation through the point of containment breach and release to the environment.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2020

Involute Working Group – Validation of CFD Turbulence Models for Steady-State Safety Analysis

There are three research reactors in the world having fuel plates curved as circle-involute (a spiral generated around a circle): The Oak Ridge National Laboratory (ORNL) High Flux Isotope Reactor (HFIR) located in Tennessee, U.S.A. [ORNL, 2018]; The Laue-Langevin Institute (ILL) High Flux Reactor (RHF) located in Grenoble, France [ILL, 2018]; and, The Technical University of Munich (TUM) Research Neutron Source Heinz Maier-Leibnitz (FRM II) located in Garching, Germany [TUM, 2018]. All three reactors are currently using Highly Enriched Uranium (HEU, 235U/U ≥ 20 wt. %) as fuel and all three are actively engaged in activities to convert to Low Enriched Uranium (LEU, 235U/U < 20 wt. %) fuel. More information on these reactors can be found in Section two. For various reasons, these reactors have expressed interest in using Computational Fluid Dynamics (CFD) tools to perform their Steady-State Thermal-Hydraulic (SSTH) safety calculations. Using CFD tools for SSTH analysis represents generally a significant departure from traditional methods and requires thorough verification and validation (V&V) to be accepted by regulators. This is why Argonne National Laboratory (ANL) and the involute reactors formed an informal group to help each other in this endeavor. This so-called Involute Working Group (IWG) aims at supporting the qualification of CFD tools for SSTH safety calculations. Activities include benchmarking, code-to-code comparison, Verification and Validation (V&V) as well as technical support. More information on the IWG and the reasons why involute-plate reactors are interested in CFD can be found in Section three. The present study aims at testing three turbulence models commonly implemented in commercial software (namely K-Epsilon, K-Omega and Reynolds-Stress) and assessing their relevance for involute-plate reactor application. An experiment performed by Gambill and Bundy in 1961 was performed to measure heat transfer coefficients in conditions prototypic of the HFIR reactor [Gambill, 1961]. More details on this experiment are provided in Section four. Based on the information available, CFD models of this experiment have been created and the experimental results have been compared to the ones obtained with these CFD models. CFD models and the results of this comparison are provided in Section five.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2018

A Detailed Assessment of the Finite Element Analysis of Reinforced Concrete Structures Impacted by Large Missiles

This document consolidates the work performed by Sandia National Laboratories and the US Nuclear Regulatory Commission in participation of Program IRIS: “Improving the Robustness of the Assessment Methodologies for Structures Impacted by Missiles”. Three round-robin benchmark exercises on improving the robustness of the assessment of structures impacted by large missiles at medium to high velocities were organized by either the IAGE Subgroup on Ageing of Concrete Structures of the Organization for Economic Co-operation and Development Nuclear Energy Agency (NEA) or Électricité de France (EDF). The objectives of the exercises were to develop guidance for conducting impact analyses including issues related to computer codes, modeling approaches, and analysis techniques. The full project was comprised of three phases: Phase I, impact of walls; Phase II, impact of larger structures; and Phase III, transmission of shock and vibration to internal components.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Research to Develop Flood Barrier Testing Strategies for Nuclear Power Plants

This report presents results of research performed by Idaho National Laboratory (INL) for the Nuclear Regulatory Commission (NRC) to identify and develop strategies for testing nuclear power plant (NPP) flood barriers. This research project began with a literature review that documented the results of previous NRC-funded research, work performed by other government agencies, reports generated by industry organizations, NRC licensee submittals, information from testing facilities, and information from commercial firms engaged in NPP decommissioning. This literature review was then used to develop a simplified flood protection categorization scheme and related glossaries. Using this categorization scheme and glossaries, a survey of currently used flood barriers was conducted, focusing on both permanent and temporary flood barriers. External flood barriers such as levees and berms were outside the scope of this survey. This research project also assessed the current state of NPP decommissioning for potential flood barrier harvesting opportunities, surveyed capabilities of domestic and international flood testing facilities, and explored technical challenges to harvesting and testing of flood barriers. Other factors that should be considered in developing flood barrier testing strategies were also explored and are documented in this report (e.g., selection of flood barriers for testing, codes and standards for flood barrier testing, potential alternatives to harvesting such as in-situ testing, testing performance criteria, and testing parameters).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Raman Microscopy Analysis of Wyoming CarbonSAFE Pilot Well Thin Sections for Mineralogy and Organic Matter Characterization

Scanning confocal Raman microspectroscopy (RMS) was used to analyze thin sections for inorganic mineral content and for evaluation of organic material (OM). Thin sections were selected from a stratigraphic test well core in the Powder River Basin, Wyoming under the Wyoming CarbonSAFE project. The well was analyzed with a variety of rock and fluid characterization techniques to determine the feasibility of a commercial-scale CO2 storage site. RMS complements other analyses, including traditional petrography, SEM, porosity and permeability, and XRD For this aspect of the study, RMS is especially important in the evaluation of OM in sealing lithologies. At prospective geologic CO2 storage sites, it is imperative to assess the unconventional oil and gas potential of seals to ensure that any future development would not compromise the integrity of the seals. Whole-slide mineralogical surveys were performed on thin sections from various shale and sandstone formations. Surveys were analyzed with Direct Classical Least Squares to identify and quantify minerals. The location and concentration of minerals was color-coded and overlaid on optical images for visualization of the distribution of minerals. Dense hyperspectral Raman mapping of OM was performed on five thin sections. Eleven spectral parameters diagnostic of organic type and thermal maturity were used to train a Partial Least Squares (PLS) calibration against a set of artificially matured samples spanning the pre- to mid-oil window. The PLS was applied to the study set and a post-mature set. Additionally, the PLS was applied to each point in hyperspectral maps for visualization of trends in maturity across sample sets and discrimination of organic matter types within a given map. In inorganic surveys on thin sections, a total of 14 unique inorganic minerals were identified in Raman spectra including quartz, dolomite, calcite, hematite and anhydrite. Shale thin sections tended to be dominated by organic material. OM was often observed mixed with inorganic minerals. Sand- and mudstones were dominated by inorganic minerals. The PLS extrapolated the post-mature set to reflectances >1.2%. The study set ranged from very immature to postmature in the median of map fit-peak parameters. However, point maturity maps indicate that matrix OM in all study samples is immature and that discreet organic particles selected for mapping, which may be inertinites, bias medians towards more-mature. The work demonstrates the capabilities of RMS to perform both whole-slide mineralogy and OM analysis with applications to formation evaluation in oil & gas, carbon sequestration and mining. Here, analysis of sealing formations in the well indicates high levels of immature organic matter that would not be a viable target for future oil production that could compromise the CO2 storage site. The work brings together diverse disciplines from geology and petrography to analytical chemistry, big data and microscopy.

Myers, Grant↗
DOE OSTI2016

PFLOTRAN Integrated Performance Assessment at Waste Isolation Pilot Plant

The Waste Isolation Pilot Plant (WIPP) is developed by the U.S. Department of Energy (DOE) for the deep geologic disposal of transuranic (TRU) waste in southeast New Mexico. The purpose of performance assessment (PA) is to estimate probability of potential radionuclide releases from the repository to the environment for a regulatory period of 10,000 years after its closure. The use of PFLOTRAN has been proposed to enhance PA with more physically realistic flow and transport models in 3D. PFLOTRAN is a state-of-art massively parallel subsurface flow and reactive transport code. This paper introduces the progress of implementation of PFLOTRAN in the PA and the preliminary results for undisturbed repository conditions. (authors)

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

Validation of single-fluid and two-fluid magnetohydrodynamic models of the helicity injected torus spheromak experiment with the NIMROD code

We present a comparison study of 3-D pressureless resistive MHD (rMHD) and 3-D presureless two-fluid MHD models of the Helicity Injected Torus with Steady Inductive helicity injection (HIT-SI). HIT-SI is a current drive experiment that uses two geometrically asymmetric helicity injectors to generate and sustain toroidal plasmas. The comparable size of the collisionless ion skin depth d{sub i} to the resistive skin depth predicates the importance of the Hall term for HIT-SI. The simulations are run with NIMROD, an initial-value, 3-D extended MHD code. The modeled plasma density and temperature are assumed uniform and constant. The helicity injectors are modeled as oscillating normal magnetic and parallel electric field boundary conditions. The simulations use parameters that closely match those of the experiment. The simulation output is compared to the formation time, plasma current, and internal and surface magnetic fields. Results of the study indicate 2fl-MHD shows quantitative agreement with the experiment while rMHD only captures the qualitative features. The validity of each model is assessed based on how accurately it reproduces the global quantities as well as the temporal and spatial dependence of the measured magnetic fields. 2fl-MHD produces the current amplification (I{sub tor}/I{sub inj}) and formation time τ{sub f} demonstrated by HIT-SI with similar internal magnetic fields. rMHD underestimates (I{sub tor}/I{sub inj}) and exhibits much a longer τ{sub f}. Biorthogonal decomposition (BD), a powerful mathematical tool for reducing large data sets, is employed to quantify how well the simulations reproduce the measured surface magnetic fields without resorting to a probe-by-probe comparison. BD shows that 2fl-MHD captures the dominant surface magnetic structures and the temporal behavior of these features better than rMHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Assessment of advanced solvent-based post-combustion CO 2 capture processes using a bi-objective optimization technique

We report the optimized performance of two advanced CO 2 capture processes is compared to that of a monoethanolamine (MEA) baseline for a gas-powered CO 2 capture retrofit of an existing coal-fired facility. The advanced temperature-swing processes utilize piperazine and mixed-salt solvents. The mixed-salt treatment involves the use of ammonia for CO 2 absorption and potassium carbonate primarily to control ammonia slip. The processes are represented in terms of energy duty requirements within a modular heat integration code developed for CO 2 capture modeling and optimization. The model includes a baseload coal plant, a gas-fired subsystem containing gas turbines and a heat recovery steam generator (HRSG), and a CO 2 capture facility. A formal bi-objective optimization procedure is applied to determine the design (e.g., detailed HRSG components and pressure levels, gas turbine capacity, CO 2 capture capacity) and time-varying operations of the facility to simultaneously maximize net present value (NPV) and minimize total capital requirement (TCR), while meeting a maximum CO 2 emission intensity constraint. For a realistic scenario constructed using historical data, optimization results indicate that both advanced processes outperform MEA in both objectives, and the mixed-salt process in turn outperforms the piperazine process. Specifically, for the scenario considered, the base case mixed-salt process achieves 16% greater NPV and 14% lower TCR than the MEA process, and 10% greater NPV and 5% lower TCR than the piperazine process. A five-case sensitivity study of the mixed-salt process indicates that it is competitive with the piperazine process and consistently outperforms the MEA process.

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

User Friendly Web-Based Tool to Assess the Energy Efficiency and Durability of Residential Wall Retrofits

According to the U.S. Department of Energy Windows and Building Envelope Research and Development Roadmap for Emerging Technologies, building envelope wall energy loss in the United States accounts for about 5.9 x 1010 kWh or 2 quads of energy annually, costing homeowners and occupants billions of dollars. Enclosure retrofits targeting these losses can save significant energy, reduce greenhouse gas emissions, and save occupants millions of dollars over time. Older homes, built before 1992 when DOE’s Building Energy Codes Program was established, represent approximately 68 percent of the residential building stock in the country, often having significant air leakage and inadequate insulation. Homes with little to no air sealing or insulation have heating and cooling losses that can represent a substantial portion of utility bills.High-performance building envelope retrofit systems are rarely selected for retrofit applications. Current solutions are expensive and/or unfamiliar to many designers, builders, contractors, and code officials and therefore are perceived as risky. The dominant perceived risk is durability specifically related to condensation and moisture accumulation in the building envelope component.The Building Science Advisor (BSA) is a rule-based expert system web-based tool that was originally developed to assist building professionals in designing energy efficient and durable wall systems for new construction. With the present focus being placed on upgrading the existing building stock, a retrofit module has been developed that, based on the location, existing construction, and planned retrofit strategy, provides recommendations on how to address the retrofit in a manner that will perform in accordance with IECC 2018 building code and not create a durability problem. This paper will describe the development of this tool and demonstrate its features and capabilities.

Desjarlais, Andre Omer↗
DOE OSTI2015

Water Governance Process Assessment: Evaluating the Link between Decision Making Processes and Outcomes in the Columbia River Basin

Academics and practitioners agree that in water governance, the quality of a decision making process should influence the quality of the outcome and the degree to which it is accepted by interested parties. However, finding a feasible way to evaluate and then improve the quality of a decision making process has proven elusive. Systematically collecting evidence of a link between process and outcome is also challenging. In my dissertation, I developed a synthesis framework for evaluating and improving water governance decision making to address these two challenges. The synthesis framework, which I call the Water Governance Process Assessment (Water GPA), draws upon 22 existing frameworks rooted in resilience, adaptive governance, and good governance. From these frameworks, I identified and provided a way to evaluate four characteristics critical to good water governance decision making processes: 1) accountability, 2) inclusivity, and 3) information, and 4) context. I applied the Water GPA framework to the recent reviews of the Columbia River Treaty by the United States and Canada. I collected data for the case studies through semi-structured interviews and surveys of process participants from the federal agencies, Tribal and First Nations, state/provincial governments, local governments, stakeholder interests, and citizens. I coded and analyzed the interviews using the qualitative analysis software QSR NVivo and the characteristics identified in the Water GPA framework. I used the two case study applications to demonstrate how to use the framework. I identified what aspects of the four process categories served as barriers and building blocks for good water governance in each water governance process. I also gleaned lessons learned and recommendations including some for determining process leadership, ensuring meaningful engagement and inclusivity, addressing sovereignty issues, setting decision criteria, sharing decision authority, allocating resources in future processes in the US and BC portions of the Columbia River Basin and similar basins. In my case studies, I also investigated which characteristics of the water governance decision making process influenced the direct outcomes of those processes (the decisions) as well as other non-target outcomes (such as trust, co-produced science, new coalitions, etc.). All four characteristics outlined in the Water Governance Process Assessment (Water GPA), played some role in the development of the content and/or support of the two case study decision documents (US Regional Recommendation and BC Provincial Decision). Generally, when the characteristic of the process was done well, it improved the legitimacy and acceptance of the decision. At the same time performing poorly in one area did not necessarily torpedo the process or decision. In both case studies, it appears that the interplay between the process characteristics (that is where two or more characteristics converge) had a greater influence over the decision. Further work is needed to clearly identify what characteristics of a process are most influential in different situations. The Water GPA is one useful tool for this effort.

13 HYDRO ENERGY↗
DOE OSTI2017

Complexity Reduction in Large Quantum Systems: Fragment Identification and Population Analysis via a Local Optimized Minimal Basis

We present, within Kohn-Sham Density Functional Theory calculations, a quantitative method to identify and assess the partitioning of a large quantum mechanical system into fragments. We then introduce a simple and efficient formalism (which can be written as generalization of other well-known population analyses) to extract, from first principles, electrostatic multipoles for these fragments. The corresponding fragment multipoles can in this way be seen as reliable (pseudo-) observables. By applying our formalism within the code BigDFT, we show that the usage of a minimal set of in-situ optimized basis functions is of utmost importance for having at the same time a proper fragment definition and an accurate description of the electronic structure. With this approach it becomes possible to simplify the modeling of environmental fragments by a set of multipoles, without notable loss of precision in the description of the active quantum mechanical region. Furthermore, this leads to a considerable reduction of the degrees of freedom by an effective coarsegraining approach, eventually also paving the way towards efficient QM/QM and QM/MM methods coupling together different levels of accuracy.

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

Non-local thermodynamic equilibrium 1.5D modeling of red giant stars

Spectra for two-dimensional (2D) stars in the 1.5D approximation are created from synthetic spectra of one-dimensional (1D) non-local thermodynamic equilibrium (NLTE) spherical model atmospheres produced by the PHOENIX code. The 1.5D stars have the spatially averaged Rayleigh-Jeans flux of a K3-4 III star while varying the temperature difference between the two 1D component models (ΔT {sub 1.5D}) and the relative surface area covered. Synthetic observable quantities from the 1.5D stars are fitted with quantities from NLTE and local thermodynamic equilibrium (LTE) 1D models to assess the errors in inferred T {sub eff} values from assuming horizontal homogeneity and LTE. Five different quantities are fit to determine the T {sub eff} of the 1.5D stars: UBVRI photometric colors, absolute surface flux spectral energy distributions (SEDs), relative SEDs, continuum normalized spectra, and TiO band profiles. In all cases except the TiO band profiles, the inferred T {sub eff} value increases with increasing ΔT {sub 1.5D}. In all cases, the inferred T {sub eff} value from fitting 1D LTE quantities is higher than from fitting 1D NLTE quantities and is approximately constant as a function of ΔT {sub 1.5D} within each case. The difference between LTE and NLTE for the TiO bands is caused indirectly by the NLTE temperature structure of the upper atmosphere, as the bands are computed in LTE. We conclude that the difference between T {sub eff} values derived from NLTE and LTE modeling is relatively insensitive to the degree of the horizontal inhomogeneity of the star being modeled and largely depends on the observable quantity being fit.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

MCNP Simulation of Discrete Gamma-ray Spectra for PGNAA Applications

Prompt Gamma Neutron Activation (PGNAA) systems use the neutron-induced gamma-ray spectra from test objects in order to identify their compositions. The U.S. Army in particular uses these systems to identify the contents of field-recovered munitions in order to determine the presence of chemical warfare materiel (CWM) or explosives. These field-recovered munitions can have a variety of different fills, including experimental fills, have a large range of sizes, and be assessed in different field conditions. All these variations can have a significant effect on the neutron-induced gamma-ray spectrum and therefore a large impact on automatic fill identification algorithms. In order to develop robust identification algorithms, parametric studies of the effects of munition size, fill composition, and assessment conditions are very useful but time-consuming and expensive in the laboratory. Monte Carlo modeling of the system response to these variations is a very attractive alternative, provided the modeling accurately reflects real gamma spectra and can be performed quickly and inexpensively when compared to experiment. At Idaho National Laboratory we have used the Monte Carlo N-Particle (MCNP6) code [2] to calculate the response of the Portable Isotopic Neutron Spectroscopy (PINS) system to a number of simulated chemical warfare munitions as well as performed laboratory measurements with the system. Two versions of PINS system were used; one excited the test object with neutrons from a 252C f source and one with neutrons from a portable deuterium-deuterium (DD) neutron generator. The sets of measurements and simulations with the two different neutron sources provided experimental data for development of automatic fill identification algorithms as well as providing a good test of the agreement between the MCNP calculations and experiment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Metagenomic analysis of the microbiota in the highly compartmented hindguts of six wood- or soil-feeding higher termites

Termites are important contributors to carbon and nitrogen cycling in tropical ecosystems. Higher termites digest lignocellulose in various stages of humification with the help of an entirely prokaryotic microbiota housed in their compartmented intestinal tract. Previous studies revealed fundamental differences in community structure between compartments, but the functional roles of individual lineages in symbiotic digestion are mostly unknown. Furthermore, we conducted a highly resolved analysis of the gut microbiota in six species of higher termites that feed on plant material at different levels of humification. Combining amplicon sequencing and metagenomics, we assessed similarities in community structure and functional potential between the major hindgut compartments (P1, P3, and P4). Cluster analysis of the relative abundances of orthologous gene clusters (COGs) revealed high similarities among woodand litter-feeding termites and strong differences to humivorous species. However, abundance estimates of bacterial phyla based on 16S rRNA genes greatly differed from those based on protein-coding genes. In conclusion, the community structure and functional potential of the microbiota in individual gut compartments are clearly driven by the digestive strategy of the host. The metagenomics libraries obtained in this study provide the basis for future studies that elucidate the fundamental differences in the symbiont-mediated breakdown of lignocellulose and humus by termites of different feeding groups. Furthermore, the high proportion of uncultured bacterial lineages in all samples calls for a reference-independent approach for the correct taxonomic assignment of protein-coding genes.

59 BASIC BIOLOGICAL SCIENCES↗