Search NASA⌕ Search

SEARCH · Search NASA

Results for “Open Source Open Science”

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,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Terrain-Influenced Winds and Fire-Fire Interactions in Wildland Fire Simulations [Dissertation]

Ensemble-based approaches to prescribed fire planning cannot be supported by computational fluid dynamics based models like FIRETEC and the Wildland-Urban Interface Fire Dynamics Simulator (WFDS) because they are too computationally expensive and cannot leverage large eddy simulation approaches like CAWFE and WRF-SFIRE because they have too coarse of resolution. QUIC-Fire was developed to fill this gap but it cannot currently address complex terrain, that is typical for instance in the Western United States. This dissertation describes a variety of improvements made to QUIC-Fire and its various incorporated algorithms in an effort to make it a viable tool in simulating wildland and prescribed fires on terrain. The modifications made to QUIC-Fire are described in three chapters. The first chapter describes the extension of the diagnostic wind model QUIC-URB, the wind engine of QUIC-Fire, to a terrain-following coordinate system. The terraininfluenced winds it generates are analyzed and compared. In particular, this chapter presents the mathematical derivation of the wind solver leading to a linear system of equations that are solved through the successive over-relaxation method. The model is validated against a standard test used in previous works (the Askervein Hill) and against a new dataset from measurements in the Socorro Mountains, New Mexico. The terrain-following implementation captures the correct phenomenology for the isolated Askervein Hill, with a wind speed up at the top of the hill. The model agrees well with measurements on the upwind side of the peak, but overestimates speed-up on the downwind side of the hill. This is due to the inability of the model to generate flow separation and wake-eddy dynamics. On a common laptop, the divergence-free wind field is obtained in 6 s, making the solver appealing for coupled fire-atmosphere simulations. The Socorro Mountain is highly complex, with many cliff faces, peaks, and valleys. Although the model captures the magnitude and direction of inlet and outlet areas of the domain, it performs rather poorly in the valley region and in the regions near the steep cliffs. Hence, the model shows good agreement with data in areas of open sloped terrain but lacks in areas where flow separation and thermally driven effects may be present (neither effect is addressed in this work). In the second chapter the implementation of the terrain-following version of QUIC-URB into QUIC-Fire, and the necessary changes needed to include terrain are described. No changes to the underlying fire spread algorithm are made other than what is required to correctly account for the inclusion of terrain. Previously published FIRETEC results that use five different topographies that share the same centerline profile are compared to simulation results from the modified QUIC-Fire that use the same topographies and fuels. QUIC-Fire results show overall similar behaviors in terms of how the topographies affect fire shapes and trends in spread rates. Due to the terrain-following version of QUIC-URB being unable to generate flow separations at the crest of hills, fire spread rates in these regions across all topographies are over-predicted when compared to FIRETEC. Lateral fire growth shows similar trends with FIRETEC between topographies but does not capture the increase in spread due to a diagonal interface between grassland and forested fuel region of the domain. These results suggest that there are three algorithms within QUIC-Fire that could use improvement: how flame tilt angle is accounted for, the incorporation of non-local drag effects, and the inclusion of the wake-eddy parameterizations that are used in QUIC-URB. Lastly, the third chapter describes a modification to the initial guess used for the calculation of the QUIC-URB mass-conserved wind solution during fire simulations. The modification is aimed at improving fire-fire interactions in QUIC-Fire simulations. The modification consists of using the solution from the previous timestep as the starting point for the calculation of the solution for the next timestep. Fire-fire interactions is greatly improved by the change but a new source of error is introduced. Due to how plumes are modelled in QUIC-Fire the new solution contains errors where gaps in the plume structure are present. However, these errors are mostly limited to the upper atmosphere, where they do not affect fire behavior at the surface, and their magnitude isn’t significant enough to discount the amount of new fire phenomenology now captured in QUIC-Fire with the change.

58 GEOSCIENCES↗
DOE OSTI2022

Automated identification of local contamination in remote atmospheric composition time series

Abstract. Atmospheric observations in remote locations offer a possibility of exploring trace gas and particle concentrations in pristine environments. However, data from remote areas are often contaminated by pollution from local sources. Detecting this contamination is thus a central and frequently encountered issue. Consequently, many different methods exist today to identify local contamination in atmospheric composition measurement time series, but no single method has been widely accepted. In this study, we present a new method to identify primary pollution in remote atmospheric datasets, e.g., from ship campaigns or stations with a low background signal compared to the contaminated signal. The pollution detection algorithm (PDA) identifies and flags periods of polluted data in five steps. The first and most important step identifies polluted periods based on the derivative (time derivative) of a concentration over time. If this derivative exceeds a given threshold, data are flagged as polluted. Further pollution identification steps are a simple concentration threshold filter, a neighboring points filter (optional), a median, and a sparse data filter (optional). The PDA only relies on the target dataset itself and is independent of ancillary datasets such as meteorological variables. All parameters of each step are adjustable so that the PDA can be “tuned” to be more or less stringent (e.g., flag more or fewer data points as contaminated). The PDA was developed and tested with a particle number concentration dataset collected during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition in the central Arctic. Using strict settings, we identified 62 % of the data as influenced by local contamination. Using a second independent particle number concentration dataset also collected during MOSAiC, we evaluated the performance of the PDA against the same dataset cleaned by visual inspection. The two methods agreed in 94 % of the cases. Additionally, the PDA was successfully applied to a trace gas dataset (CO2), also collected during MOSAiC, and to another particle number concentration dataset, collected at the high-altitude background station Jungfraujoch, Switzerland. Thus, the PDA proves to be a useful and flexible tool to identify periods affected by local contamination in atmospheric composition datasets without the need for ancillary measurements. It is best applied to data representing primary pollution. The user-friendly and open-access code enables reproducible application to a wide suite of different datasets. It is available at https://doi.org/10.5281/zenodo.5761101 (Beck et al., 2021).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Low-cost Manufacturing of Semitransparent CdTe PV for Building Integration

Solar has been demonstrated to be a robust renewable energy source, constituting a significant portion of the United States’ renewable energy portfolio. Despite its growth across residential, commercial, and utility sectors over the past decades, it remains a small fraction of the overall energy infrastructure. Challenges persist in fully harnessing solar power to meet the nation's escalating energy demands. Among these challenges lies the hurdle of efficiently distributing large quantities of solar-generated electricity to densely populated regions with the highest energy needs and costs. Traditional utility-scale arrays demand extensive land, a luxury often unavailable in metropolitan areas. Consequently, installations must be situated at a distance, necessitating additional infrastructure for electricity transmission to service areas. While metropolitan landscapes lack sprawling open spaces suitable for conventional utility-scale solar deployment, they offer a different resource: windows. Semitransparent photovoltaic window technology has the potential to not only bolster the grid's energy capacity, but to also provide HVAC and economic advantages to building owners. However, commercial availability of building-integrated photovoltaic windows remains limited. Silicon based photovoltaics currently dominate the solar market but adapting them for use in windows poses a variety of engineering and economic challenges such as relatively low power density, high costs associated with custom manufacturing, and aesthetic considerations. Addressing these challenges, this project explored the use of laser ablation patterning to manufacture cost-effective, high-efficiency, semitransparent Cadmium Telluride photovoltaic modules. Results showcased the potential of this methodology in developing photovoltaic windows and other innovative semitransparent PV applications. The ablation manufacturing technique demonstrated great versatility in achieving different patterns and levels of visible light transmission, and the power loss due to ablation was nearly directly proportional to the amount of material removed. Furthermore, the manufacturing process for Cadmium Telluride modules already has established advantages in material and energy efficiency, and the conversion of a standard submodule to semitransparent essentially requires a single additional process step, ensuring scalability. It is important to note that during this work, Toledo Solar experienced substantial organizational upheaval stemming from a lawsuit with First Solar. An external investigation led the board of directors to remove and replace the previous management team, and several other members of the staff elected to depart as well, including the then acting Principle Investigator on this project. The remaining Toledo Solar team attempted to recover from the disruption and deliver on the remaining tasks, but upon its own review, the Department of Energy ruled the project in default and terminated the contract in December 2023.

14 SOLAR ENERGY↗
DOE OSTI2018

Rapid growth of organic aerosol nanoparticles over a wide tropospheric temperature range

Nucleation and growth of aerosol particles from atmospheric vapors constitutes a major source of global cloud condensation nuclei (CCN). The fraction of newly formed particles that reaches CCN sizes is highly sensitive to particle growth rates, especially for particle sizes <10 nm, where coagulation losses to larger aerosol particles are greatest. Currnet results show that some oxidation products from biogenic volatile organic compounds are major contributors to particle formation and initial growth. However, whether oxidized organics contribute to particle growth over the broad span of tropospheric temperatures remains an open question, and quantitative mass balance for organic growth has yet to be demonstrated at any temperature. Here, in experiments performed under atmospheric conditions in the Cosmics Leaving Outdoor Droplets (CLOUD) chamber at the European Organization for Nuclear Research (CERN), we show that rapid growth of organic particles occurs over the range from – 25 ° C to 25 ° C. The lower extent of autoxidation at reduced temperatures is compensated by the decreased volatility of all oxidized molecules. This is confirmed by particle-phase composition measurements, showing enhanced uptake of relatively less oxygenated products at cold temperatures. We can reproduce the measured growth rates using an aerosol growth model based entirely on the experimentally measured gas-phase spectra of oxidized organic molecules obtained from two complementary mass spectrometers. We report that the growth rates are sensitive to particle curvature, explaining widespread atmospheric observations that particle growth rates increase in the single-digit-nanometer size range. Our results demonstrate that organic vapors can contribute to particle growth over a wide range of tropospheric temperatures from molecular cluster sizes onward.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Distinct Fracture Mineralogy That is Out of Equilibrium With Modern Groundwaters Provides Important Context for Subsurface Life

Rock fracture surfaces in the crust are essential habitat for microorganisms. Fracture‐groundwater interfaces provide physical substrates for biofilm growth and are sources of carbon, nutrients, and electron donors and acceptors. To better understand geochemical processes impacting fracture surfaces and the subsurface microbiome, we identified fractures in archived rock cores from the Soudan formation, which is known to host saline groundwaters and isolated microbial communities dependent on rock‐water interactions. Cores with open fractures were thin sectioned and studied via electron microprobe and synchrotron X‐ray fluorescence microprobe. Most fracture surfaces had mineralogy distinct from that of the bulk rock. Chlorite minerals were abundant on fracture surfaces and had elemental compositions suggesting deposition during late‐stage hydrothermal alteration. Fracture‐lining chlorites likely limit access to iron oxide and sulfide minerals that are active in subsurface biogeochemical cycles. Calcium‐rich rinds were also observed along fracture edges. These rinds were too thin and poorly ordered to be identified via light microscopy or X‐ray diffraction; however, Ca K‐edge micro‐X‐ray absorption near‐edge structure spectroscopy identified them as carbonates, minerals not observed in the bulk rock. Thermodynamic modeling shows that carbonate precipitation is largely unfavorable in Soudan groundwaters, indicating that fracture edge conditions differed from those in modern water samples. Because of the low carbon concentrations in Soudan groundwaters, carbonate rinds likely play an important role in subsurface carbon cycling and may mark fracture surfaces that once hosted biofilms. Overall, this study suggests that fracture alteration can both play an active role in and suppress rock‐water interactions essential to subsurface life.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Evaluation of Data Catalog Software for Hanford Site Environmental Datasets

Environmental information and data underpin achievement of the U.S. Department of Energy (DOE) Office of Environmental Management (EM) mission at the Hanford Site. The Hanford Environmental Data Management (HEDM) Program is the DOE Richland Operations Office (RL) approach to develop and implement a formal program for managing environmental data and the associated records, materials, and systems at the Hanford Site. The current project, contract, organization, and contractor-specific efforts at managing environmental data sets are insufficient to provide orderly, long-term, site-wide access. A vital element to be created within the HEDM program plan is a catalog of data sources, called the Hanford Environmental Information and Data Index (HEIDI), that will enable long-term access and retrievability for the multiple independent sources of data that might otherwise be difficult to discover. This report compares leading open source and commercial data catalog platforms using criteria to assess the functionality needed to develop the HEIDI catalog of Hanford data sources that connects and exchanges data with established Hanford Local Area Network (HLAN) enterprise information technology systems. Proprietary platforms evaluated included ArcGIS Enterprise Sites, Junar, OpenDataSoft, and Socrata, and non-proprietary platforms included Energy Data eXchange (EDX), Comprehensive Knowledge Archive Network (CKAN), and DKAN (a Drupal-based open data portal based on CKAN). Capabilities supporting data discoverability, retrieval, and archival, as well as metadata standard requirements and integration into the HLAN were rated as either failing to meet requirements (F), meeting requirements (M), or exceeding requirements by delivering additional desired features (E). The lowest rating for any capability area was assigned as the overall rating for the platform. These findings enable DOE-RL and the contractors implementing the HEDM plan to focus on candidate tools likely to meet the requirements for implementing HEIDI. All of the platforms receiving an overall rating of ‘F’ were unable to be deployed on Hanford infrastructure or within dedicated cloud resources. A propriety software-as-a-service (SaaS) model of delivering a data catalog (e.g., found in software such as Junar and OpenDataSoft) favors consistency across customers at the expense of customization and configurable roles that are needed for Hanford work. Hosting data on a shared commercial platform places limits on dataset size (maximum of 240 Mb for OpenDataSoft), a significant limitation for HEIDI implementation. EDX, a government data catalog based on CKAN, received the ‘F’ rating due to an inability to incorporate authentication from HLAN into the system. Among platforms rated ‘M’ or ‘E’, only the Socrata platform had a SaaS delivery model. In contrast to other SaaS platforms, Socrata provided custom roles and gateways that allow local datasets to be incorporated into an online catalog. Socrata also complies with the Federal Risk and Authorization Management Program, a significant benefit for cloud-based management of Hanford data. The other platforms rated ‘M’ or ‘E’, ArcGIS Enterprise Sites, CKAN, and DKAN, provide fully self-hosted options, allowing for greater control and flexibility with the HEIDI catalog. These widely used tools have supportive communities of practice, extensive customization options, and demonstrated deployments that provide evidence that they can meet requirements, often deliver additional desired features, and work well with federal government systems. Completely customized alternatives built on a collection of applications were not evaluated because achieving similar performance to CKAN or DKAN requires substantial resources, especially in the absence of the active communities that have grown to support these tools. ArcGIS Enterprise Sites, Socrata, CKAN, and DKAN were evaluated as strong candidates for successful implementation with HEIDI.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Real-Time Partitioning of Diurnal Stem CO 2 Efflux into Local Stem Respiration and Xylem Transport Processes

The apparent respiratory quotient (ARQ) of tree stems, defined as the ratio of net stem CO 2 efflux (E S_CO2 ) to net stem O 2 influx (E S_O2 ), offers insights into the balance between local respiratory CO 2 production and CO 2 transported via the xylem. Traditional static chamber methods for measuring ARQ can introduce artifacts and obscure natural diurnal variations. Here, we employed an open flow-through stem chamber with ambient air coupled with cavity ring-down spectrometry, which uses the molecular properties of CO 2 and O 2 molecules to continuously measure E S_CO2 , E S_O2 , and ARQ, at the base of a California cherry tree (Prunus ilicifolia) during the 2024 growing season. Measurements across three stem chambers over 3–11-day periods revealed strong correlations between E S_CO2 and E S_O2 and mean ARQ values ranging from 1.3 to 2.9, far exceeding previous reports. Two distinct diurnal ARQ patterns were observed: daytime suppression with nighttime recovery, and a morning peak followed by gradual decline. Partitioning E S_CO2 into local respiration and xylem-transported CO 2 indicated that the latter can dominate when ARQ exceeds 2.0. Furthermore, transported CO 2 exhibited a higher temperature sensitivity than local respiration, with both processes showing declining temperature sensitivity above 20 °C. These findings underscore the need to differentiate stem CO 2 flux components to improve our understanding of whole-tree carbon cycling.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Fluid-Kinetic Coupling: Advanced Discretizations for Simulations on Emerging Heterogeneous Architectures (LDRD FY20-0643)

Plasma physics simulations are vital for a host of Sandia mission concerns, for fundamental science, and for clean energy in the form of fusion power. Sandia's most mature plasma physics simulation capabilities come in the form of particle-in-cell (PIC) models and magnetohydrodynamics (MHD) models. MHD models for a plasma work well in denser plasma regimes when there is enough material that the plasma approximates a fluid. PIC models, on the other hand, work well in lower-density regimes, in which there is not too much to simulate; error in PIC scales as the square root of the number of particles, making high-accuracy simulations expensive. Real-world applications, however, almost always involve a transition region between the high-density regimes where MHD is appropriate, and the low-density regimes for PIC. In such a transition region, a direct discretization of Vlasov is appropriate. Such discretizations come with their own computational costs, however; the phase-space mesh for Vlasov can involve up to six dimensions (seven if time is included), and to apply appropriate homogeneous boundary conditions in velocity space requires meshing a substantial padding region to ensure that the distribution remains sufficiently close to zero at the velocity boundaries. Moreover, for collisional plasmas, the right-hand side of the Vlasov equation is a collision operator, which is non-local in velocity space, and which may dominate the cost of the Vlasov solver. The present LDRD project endeavors to develop modern, foundational tools for the development of continuum-kinetic Vlasov solvers, using the discontinuous Petrov-Galerkin (DPG) methodology, for discretization of Vlasov, and machine-learning (ML) models to enable efficient evaluation of collision operators. DPG affords several key advantages. First, it has a built-in, robust error indicator, allowing us to adapt the mesh in a very natural way, enabling a coarse velocity-space mesh near the homogeneous boundaries, and a fine mesh where the solution has fine features. Second, it is an inherently high-order, high-intensity method, requiring extra local computations to determine so-called optimal test functions, which makes it particularly suited to modern hardware in which floating-point throughput is increasing at a faster rate than memory bandwidth. Finally, DPG is a residual-minimizing method, which enables high-accuracy computation: in typical cases, the method delivers something very close to the $L^2$ projection of the exact solution. Meanwhile, the ML-based collision model we adopt affords a cost structure that scales as the square root of a standard direct evaluation. Moreover, we design our model to conserve mass, momentum, and energy by construction, and our approach to training is highly flexible, in that it can incorporate not only synthetic data from direct-simulation Monte Carlo (DSMC) codes, but also experimental data. We have developed two DPG formulations for Vlasov-Poisson: a time-marching, backward-Euler discretization and a space-time discretization. We have conducted a number of numerical experiments to verify the approach in a 1D1V setting. In this report, we detail these formulations and experiments. We also summarize some new theoretical results developed as part of this project (published as papers previously): some new analysis of DPG for the convection-reaction problem (of which the Vlasov equation is an instance), a new exponential integrator for DPG, and some numerical exploration of various DPG-based time-marching approaches to the heat equation. As part of this work, we have contributed extensively to the Camellia open-source library; we also describe the new capabilities and their usage. We have also developed a well-documented methodology for single-species collision operators, which we applied to argon and demonstrated with numerical experiments. We summarize those results here, as well as describing at a high level a design extending the methodology to multi-species operators. We have released a new open-source library, MLC, under a BSD license; we include a summary of its capabilities as well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Soil properties and root characteristics across four lowland Panamanian forests from 0 - 1 m soil depths

Objectives:Fine roots significantly influence ecosystem-scale cycling of nutrients, carbon (C), and water, yet there is limited understanding of how fine root traits vary across and within tropical forests, some of Earth's most C-rich ecosystems. The biomass of fine roots can impact soil carbon storage, as root mortality is a primary source of new carbon to soils. A positive relationship has been observed between fine root biomass and soil carbon stocks in Panama (Cusack et al 2018). Beyond biomass, root characteristics like specific root length (SRL) could also influence soil carbon, as roots with higher SRL are less dense and thinner, potentially decomposing more easily or promoting soil aggregation. Understanding the effects of root morphology and tissue quality on soil carbon storage and with soil properties in general can improve predictions of landscape-scale carbon patterns. We aggregated new data of root biomass, morphology and nutrient content at 0-10 cm, 10-20 cm, 20-50 cm and 50-100 cm depth increments across four distinct lowland Panamanian forests and paired with already published datasets (Cusack et al 2018; Cusack and Turner 2020) of soil chemistry from the same sites and soil depths to explore relationship between soil carbon stocks and root characteristics.Datasets included:The datasets provided include .csv and .xlsx files for fine root characteristics and soil chemistry from four different forests across 0-10 cm, 10-20 cm, 20-50 cm, and 50-100 cm depth increments. Root characteristics include live fine root biomass, dead fine root biomass, coarse root biomass, specific root length, root diameter, root tissue density, specific root area, root %N, root %C, and root C/N ratio. Soil chemistry data includes total carbon (TC), dissolved organic carbon (DOC), bulk density, total phosphorus (TP), available phosphorus (AEM Pi), and various Mehlich-extractable elements such as aluminum, calcium, iron, potassium, manganese, phosphorus, and zinc. Nitrogen content measures include ammonium, nitrate, total dissolved nitrogen (TDN), dissolved inorganic nitrogen (DIN), and dissolved organic nitrogen (DON). The dataset also includes total exchangeable bases (TEB) and effective cation exchange capacity (ECEC) in both centimoles of charge per kilogram and micromoles of charge per gram. The soil chemistry data was obtained from Cusack et al (2018) and Cusack and Turner (2020) and paired with root characteristics data for the same depth increments and sites. Additionally, a .kml file is provided with coordinates for all 32 plots included in the study across four forests (n = 8 plots per site). Root data was averaged across these 8 plots per site and soil data was collected in one pit in each site. This dataset serves as baseline data before a throughfall exclusion experiment, Panama Rainforest Changes with Experimental Drying (PARCHED), was implemented. No special software is needed to open these files.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Bi-Modal Model for Neutron Emissions from PuO{sub 2} and MOX Holdup

The measurement of uranium and plutonium holdup in plants during process activity and for decommissioning is important for nuclear safeguards and material control. The amount of plutonium and uranium holdup in glove-boxes, pipes, ducts, and other containers has been measured for several decades using both neutron and gamma-ray techniques. For the larger containers such as hot cells and glove-boxes that contain processing equipment, the gamma-ray techniques are limited by self-shielding in the sample as well as gamma absorption in the equipment and associated shielding. The neutron emission is more penetrating and has been used extensively to measure the holdup for the large facilities such as the MOX processing and fabrication facilities in Japan and Europe. In some case the totals neutron emission rates are used to determine the holdup mass and in other cases the coincidence rates are used such as at the PFPF MOX fabrication plant in Japan. The neutron emission from plutonium and MOX has 3 primary source terms: 1) Spontaneous fission (SF) from the plutonium isotopes, 2) The (α,n) reactions from the plutonium alpha particle emission reacting with the oxygen and other impurities, and 3) Neutron multiplication (M) in the plutonium and uranium as a result of neutrons created by the first two sources. The spontaneous fission yield per gram is independent of thickness, whereas, the above sources 2) and 3) are very dependent on the thickness of the deposit. As the effective thickness of the deposit becomes thin relative to the alpha particle range, the (α,n) reactions and neutrons from multiplication (M) approach zero. In any glove-box, there will always be two primary modes of holdup accumulation, namely direct powder contact and non-contact by air dispersal. These regimes correspond to surfaces in the glove-box that have come into direct contact with the process MOX powder versus surface areas that have not had direct contact with the powder. The air dispersal of PuO{sub 2} particles has been studied for several decades by health physicists, because the primary health hazard of plutonium is breathing the airborne particles. The air dispersal mechanism results from the smaller particles in the top layer of powder that are lifted into the air by the electrostatic charge buildup from the alpha decay process, and the air convection carries the particles to new more distant locations. If there is open plutonium powder in a glove-box, the surfaces at more distant locations will become contaminated over time. The range of an alpha particle in a solid or powder is a function of the particle energy, the material density, and the atomic number A of the material. The average energy of a plutonium alpha particle is ∼5.2 MeV and the range in air is ∼37 mm. The range in other materials can be estimated via the Bragg-Kleenman equation. For plutonium, A is 94, and the typical density for a single particle is ∼11.5 g/cm{sup 3}, but for a powder, the density would be less because of the air packing fraction. The significance of the small diameter is that the range of the alpha particle is ∼50 μm for powder density 2.5 and significantly less for a single particle with density 11.5, so the thin deposit of separate small particles will have a greatly reduced (α,n) yield. The average alpha transit length to the surface in the isolated MOX particle would be < 2.5 μm; whereas, the range of the alpha particle is much longer. Thus, most of the alpha particles would escape from the MOX particle and be absorbed by the walls and air. The air dispersal particles will have access to a large surface area that includes the walls, whereas, the powder contact surface area will be orders of magnitude smaller. Thus, the vast majority of the glove-box surface area does not produce the full (α,n) reaction neutron yield, even from the O{sub 2} in the PuO{sub 2} as well as any impurity contamination such as H{sub 2}O. To obtain a more quantitative estimate of the neutron (α,n) yields as a function of holdup deposit thickness, we have used MCNPX calculations to estimate the absorption of alpha particles in PuO{sub 2} holdup deposits. The powder thickness was varied from 0.1 μm to 5000 μm and the alpha particle escape probability was calculated. As would be expected, as the thickness approaches zero, the escape probability approaches 1.0, and as the thickness gets much greater than the alpha particle range (∼50 μm), the escape probability becomes small. Typically, the neutron holdup calibration measurement are performed using sealed containers of thick MOX that has all 3 sources of neutrons [SF, (α,n), and M], and no significant impurities. Thus, the calibration counting rates need to include corrections for M and (α,n) yields that are different for the holdup compared with the calibration samples. If totals neutron counting is used for the holdup measurements, the variability of the (α,n) term needs to be considered.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Data from: "Moisture rivals temperature in limiting photosynthesis by trees establishing beyond their cold-edge range limit under ambient and warmed conditions"

This archive contains data files that were used to draw conclusions in “Moisture rivals temperature in limiting photosynthesis by trees establishing beyond their cold-edge range limit under ambient and warmed conditions”, by Moyes et al., 2015. All field research was completed in common garden plots set up as part of the Alpine Treeline Warming Experiment (ATWE) on Niwot Ridge, Colorado, USA.There are two main data file formats in this archive: comma-separated values (.csv), and Microsoft Excel (.xls and .xlsx). .xlsx files can be read using Microsoft Excel and Google Sheets, and .csv files can be read using any simple text editor program, such as TextEdit (Mac) and Notepad (Windows). This .pdf data user’s guide can be read using Adobe Acrobat Reader, or any other compatible software. Seedling photographs and their corresponding leaf area-processed images are available in .jpg/.JPG image format, and can be opened using Preview (Mac) and Photos (Windows). To provide additional spatial context, two types of geospatial files are also published in this data package: ESRI shapefiles (.shp) and .kml files. Shapefiles are compatible with any GIS software able to read the file type (such as QGIS or ESRI’s ArcGIS suite), and .kml files can be opened with Google Earth or Google Maps. Figures 3 and 4 in the publication contain data from Moyes et al. 2013. This publication is cited in the References section in this archive, and data files can be accessed via the Alpine Treeline Warming Experiment project portal on ESS-DIVE. ---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- Climate change is altering plant species distributions globally, and warming is expected to promote uphill shifts in mountain trees. However, at many cold-edge range limits, such as alpine treelines in the western United States, tree establishment may be colimited by low temperature and low moisture, making recruitment patterns with warming difficult to predict.- We measured response functions linking carbon (C) assimilation and temperature- and moisture-related microclimatic factors for limber pine (Pinus flexilis) seedlings growing in a heating × watering experiment within and above the alpine treeline. We then extrapolated these response functions using observed microclimate conditions to estimate the net effects of warming and associated soil drying on C assimilation across an entire growing season.- Moisture and temperature limitations were each estimated to reduce potential growing season C gain from a theoretical upper limit by 15–30% (c. 50% combined). Warming above current treeline conditions provided relatively little benefit to modeled net assimilation, whereas assimilation was sensitive to either wetter or drier conditions.- Summer precipitation may be at least as important as temperature in constraining C gain by establishing subalpine trees at and above current alpine treelines as seasonally dry subalpine and alpine ecosystems continue to warm.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

A Microalgae-Based Platform for the Beneficial Re-use of Carbon Dioxide Emissions from Power Plants

This project sought to address the technical and economic barriers to carbon dioxide capture and utilization using microalgae. Operating data were collected in 2016 and 2017 during cultivation of Scenedesmus acutus at Duke Energy’s East Bend Station – a coal-fired power plant located in northern Kentucky – using flue gas as the CO 2 source. Algae were grown in a 1200 L “cyclic flow” photobioreactor (PBR) designed by the University of Kentucky. A key finding was that the harvested algae contained only very low concentrations of heavy metals (As, Cd, Hg, Se), consistent with heavy metals incorporation from the supplied nutrients. This indicates that algal biomass produced from coal-derived flue gas would be suitable for a variety of applications, including the production of bioplastics, use as fertilizer, etc. A lifecycle assessment showed that the UK-designed PBR employed in this work qualifies as a net CO 2 capture technology. Indeed, over a 30-year period, net CO 2 capture would equate to 43% of the targeted amount, i.e., the amount captured from the supplied flue gas. A techno-economic analysis indicated that the minimum production cost of Scenedesmus acutus biomass in the US is in the order of 875 USD/ton, excluding the cost of capital. While this figure is not too dissimilar to values reported for open raceway ponds in similar scenarios, it emphasizes that for current cultivation technology any pathway to economic viability will require applications for which algal boimass can be sold at prices in excess of 1,000 USD/ton. Currently, such applications represent relatively small markets, such as pigments (e.g., astaxanthin) and nutraceuticals (ω-3 unsaturated fatty acids), as well as nutritional supplements (whole algae) for human consumption and for use in pet food. Consequently, the commercialization of large-scale algae-based CO 2 capture and utilization will require the development of new technologies to reduce the cost of algae production and/or the development of new, high-value applications for algal biomass. One of the more promising applications for algal biomass is in the production of bioplastics. In this work, the potential for the algae grown in the East Bend PBR for the production of bioplastic with adequate mechanical properties was clearly shown. Positive features of the produced biomass included a high protein content and a composition that was generally more homogeneous than biomass grown in open ponds (in which many species may be present). The best candidate for further review, after incorporation into ethylene-vinyl acetate, was a lipid and sugar extracted material, which showed the highest extension values with comparable load values to other UK-derived samples. It also demonstrated extension benefits against the Algix Bloom product (currently offered commercially), even though it contained agglomerates which generally exert a negative effect on mechanical properties. This leads to speculation that with enhanced milling that exists on the commercial scale, the sugar and fat extracted product may be even more competitive. This, in turn, points to the need for additional work in order to assess the properties of such optimized algae-based plastics and the price point they can command.

09 BIOMASS FUELS↗
DOE OSTI2018

Self-Navigating Robotic Assistant for Long-term Wide Area Floor Contamination Monitoring - 18213

For Routine Contamination Testing (RCT) of an operational nuclear facility or one undergoing Decontamination and Decommissioning (D and D), workers utilize specified travel corridors (enclosed hallways or areas designated with tape or rope in open spaces) to minimize the risk of spreading contamination. Radiation health technicians must continually monitor these areas for contamination to minimize worker dose and prevent its spread to uncontrolled areas. For many sites, it is impractical for the technician to adequately survey such areas. The Nuclear Robotics Group (NRG) at UT Austin has developed a self-navigating robotic platform equipped with radiation detectors to augment the technician to improve both survey coverage and frequency. While maintaining a controlled velocity and floor-to-sensor distance, the system takes radiation measurements to identify and/or localize contamination sources. Readings above background levels are reported to a remote workstation to allow technicians to then respond accordingly. The completion of these mundane tasks by a robot allows the trained technicians to focus instead on operational tasks, D and D, or responding to alarms. The robot generates a map using a Lidar range sensor with established Simultaneous Localization and Mapping (SLAM) techniques. Range data is then compared to previously generated maps for robot localization and autonomous navigation. To maximize coverage, the robot executes a Complete Coverage Path Planning (CCPP) algorithm using a static map that is discretized into cells to assure comprehensive yet efficient coverage. As the system executes the generated plan, it builds a radiation map of the area for the technician to validate and verify the results. If contamination is detected, the robot will minimize the spread of contamination by stopping, sounding an alarm, and alerting the remote technician. The proposed system allows for long-term deployment with minimal human intervention. It is robust to both highly dynamic obstacles (i.e. personnel) and low dynamic changes (i.e. temporally stored cart) in the environment. This is accomplished through continuous spatiotemporal updating of the static map, and the utilization of a supervisory controller that allows the system to reliably respond to unexpected events or states (i.e. low battery, closed door, etc.). While the supervisory controller can respond to temporary complications, the system also recognizes the environment can change slowly over time necessitating that the CCPP automatically be updated to account for observed 'permanent' changes to the environment. To do this, the robot uses intermediate probabilistic data to make intelligent planning decisions offline while still maintaining local reactive obstacle avoidance and safety features. The self-updating of the map allows the robot to run longer without the need to remake the map, leading to more consistent monitoring data and less maintenance from the human technician. This paper will review the current system hardware, the CCPP algorithm, the supervisory controller, the probabilistic algorithm for identifying permanent changes, and recent experiments at UT Austin to experimentally validate the reported sensor coverage. (authors)

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

COVID19 Disease Map, a computational knowledge repository of virus–host interaction mechanisms

We need to effectively combine the knowledge from surging literature with complex datasets to propose mechanistic models of SARS-CoV-2 infection, improving data interpretation and predicting key targets of intervention. Here, we describe a large-scale community effort to build an open access, interoperable and computable repository of COVID-19 molecular mechanisms. The COVID-19 Disease Map (C19DMap) is a graphical, interactive representation of disease-relevant molecular mechanisms linking many knowledge sources. Notably, it is a computational resource for graph-based analyses and disease modelling. To this end, we established a framework of tools, platforms and guidelines necessary for a multifaceted community of biocurators, domain experts, bioinformaticians and computational biologists. The diagrams of the C19DMap, curated from the literature, are integrated with relevant interaction and text mining databases. We demonstrate the application of network analysis and modelling approaches by concrete examples to highlight new testable hypotheses. This framework helps to find signatures of SARS-CoV-2 predisposition, treatment response or prioritisation of drug candidates. Such an approach may help deal with new waves of COVID-19 or similar pandemics in the long-term perspective.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Mondo: integrating disease terminology across communities

Precision medicine aims to enhance diagnosis, treatment, and prognosis by integrating multimodal data at the point of care. However, challenges arise due to the vast number of diseases, differing methods of classification, and conflicting terminological coding systems and practices used to represent molecular definitions of disease. This lack of interoperability artificially constrains the potential for diagnosis, clinical decision support, care outcome analysis, as well as data linkage across research domains to support the development or repurposing of therapeutics. There is a clear and pressing need for a unified system for managing disease entities⁠—including identifiers, synonyms, and definitions. To address these issues, we created the Mondo disease ontology—a community-driven, open-source, unified disease classification system that harmonizes diverse terminologies into a consistent, computable framework. Mondo integrates key medical and biomedical terminologies, including Online Mendelian Inheritance in Man (OMIM), Orphanet, Medical Subject Headings (MeSH), National Cancer Institute Thesaurus (NCIt), and more, to provide a comprehensive and accurate representation of disease concepts with fully provenanced and attributed links back to the sources. Mondo can be used as the handle for curation of gene–disease associations utilized in diagnostic applications, research applications such as computational phenotyping, and in clinical coding systems in clinical decision support by pointing the clinician to the numerous knowledge resources linked to the Mondo identifier. Mondo's community-centric approach, stewarded by the Monarch Initiative's expertise in ontologies, ensures that the ontology remains adaptable to the evolving needs of biomedical research and clinical communities, as well as the knowledge providers.

biomedical informatics↗
DOE OSTI2022

Cabana: A Performance Portable Library for Particle-Based Simulations

Particle-based simulations are ubiquitous throughout many fields of computational science and engineering, spanning the atomistic level with molecular dynamics (MD), to mesoscale particle-in-cell (PIC) simulations for solid mechanics, device-scale modeling with PIC methods for plasma physics, and massive N-body cosmology simulations of galaxy structures, with many other methods in between (Hockney & Eastwood, 1989). While these methods use particles to represent significantly different entities with completely different physical models, many low-level details are shared including performant algorithms for short- and/or long-range particle interactions, multi-node particle communication patterns, and other data management tasks such as particle sorting and neighbor list construction. Cabana is a performance portable library for particle-based simulations, developed as part of the Co-Design Center for Particle Applications (CoPA) within the Exascale Computing Project (ECP) (Alexander et al., 2020). The CoPA project and its full development scope, including ECP partner applications, algorithm development, and similar software libraries for quantum MD, is described in (Mniszewski et al., 2021). Cabana uses the Kokkos library for on-node parallelism (Edwards et al., 2014; Trott et al., 2022), enabling simulation on multi-core CPU and GPU architectures, and MPI for GPU-aware, multi-node communication. Cabana provides particle simulation capabilities on almost all current Kokkos backends, including serial execution, OpenMP (including OpenMP-Target for GPUs), CUDA (NVIDIA GPUs), HIP (AMD GPUs), and SYCL (Intel GPUs), providing a clear path for the coming generation of accelerator-based exascale hardware. Cabana builds on Kokkos by providing new particle data structures and particle algorithms resulting in a similar execution policy-based, node-level programming model that is intended to be used in addition to the core Kokkos library within an application. Cabana is designed as an application and physics agnostic, but particle-specific toolkit which can either be used to generate a new application, or to be used as needed in existing applications at various levels of invasiveness including through interfaces that wrap user memory in existing data structures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Coastal Ohio Wind Project

The Coastal Ohio Wind Project intends to address problems that impede deployment of wind turbines in the coastal and offshore regions of Northern Ohio. The project evaluates different wind turbine designs and the potential impact of offshore turbines on migratory and resident birds by developing multidisciplinary research, which involves wildlife biology, electrical and mechanical engineering, and geospatial science. Firstly, the project conducts cost and performance studies of two- and three-blade wind turbines using a turbine design suited for the Great Lakes. The numerical studies comprised an analysis and evaluation of the annual energy production of two- and three-blade wind turbines to determine the levelized cost of energy. This task also involved wind tunnel studies of model wind turbines to quantify the wake flow field of upwind and downwind wind turbine-tower arrangements. The experimental work included a study of a scaled model of an offshore wind turbine platform in a water tunnel. The levelized cost of energy work consisted of the development and application of a cost model to predict the cost of energy produced by a wind turbine system placed offshore. The analysis found that a floating two-blade wind turbine presents the most cost effective alternative for the Great Lakes. The load effects studies showed that the two-blade wind turbine model experiences less torque under all IEC Standard design load cases considered. Other load effects did not show this trend and depending on the design load cases, the two-bladed wind turbine showed higher or lower load effects. The experimental studies of the wake were conducted using smoke flow visualization and hot wire anemometry. Flow visualization studies showed that in the downwind turbine configuration the wake flow was insensitive to the presence of the blade and was very similar to that of the tower alone. On the other hand, in the upwind turbine configuration, increasing the rotor blade angle of attack reduced the wake size and enhanced the vortices in the flow downstream of the turbine-tower compared with the tower alone case. Mean and rms velocity distributions from hot wire anemometer data confirmed that in a downwind configuration, the wake of the tower dominates the flow, thus the flow fields of a tower alone and tower-turbine combinations are nearly the same. For the upwind configuration, the mean velocity shows a narrowing of the wake compared with the tower alone case. The downwind configuration wake persisted longer than that of an upwind configuration; however, it was not possible to quantify this difference because of the size limitation of the wind tunnel downstream of the test section. The water tunnel studies demonstrated that the scale model studies could be used to adequately produce accurate motions to model the motions of a wind turbine platform subject to large waves. It was found that the important factors that affect the platform is whether the platform is submerged or surface piercing. In the former, the loads on the platform will be relatively reduced whereas in the latter case, the structure pierces the wave free surface and gains stiffness and stability. The other important element that affects the movement of the platform is depth of the sea in which the wind turbine will be installed. Furthermore, the wildlife biology component evaluated migratory patterns by different monitoring systems consisting of marine radar, thermal IR camera and acoustic recorders. The types of radar used in the project are weather surveillance radar and marine radar. The weather surveillance radar (1988 Doppler), also known as Next Generation Radar (NEXRAD), provides a network of weather stations in the US. Data generated from this network were used to understand general migratory patterns, migratory stopover habitats, and other patterns caused by the effects of weather conditions. At a local scale our marine radar was used to complement the datasets from NEXRAD and to collect additional monitoring parameters such as passage rates, flight paths, flight directions, and flight altitudes of nocturnal migrating species. Our work focused on the design and development of custom built marine radar that used t-bar and parabolic dish antennas. The marine radar used in the project was Furuno (XANK250) which was coupled with a XIR3000B digitizing card from Russell Technologies for collection of the radar data. The radar data was processed by open source radR processing software using different computational techniques and methods. Additional data from thermal IR imaging cameras were collected to detect heat emitted from objects and provide information on movements of birds and bats, data which we used for different animal flight behavior analysis. Lastly, the data from the acoustic recorders were used to provide the number of bird calls for assessing patterns and peak passage rates during migration. The development of the geospatial database included collection of different data sources that are used to support offshore wind turbine development. Many different data sets were collected and organized using initial version of web-based repository software tools that can accommodate distribution of rectified pertinent data sets such as the lake depth, lake bottom engineering parameters, extent of ice, navigation pathways, wind speed, important bird habitats, fish efforts and other layers that are relevant for supporting robust offshore wind turbine developments. Additional geospatial products developed during the project included few different prototypes for offshore wind farm suitability which can involve different stakeholders and participants for solving complex planning problems and building consensus. Some of the prototypes include spatial decision support system (SDSS) for collaborative decision making, a web-based Participatory Geographic Information System (PGIS) framework for evaluating importance of different decision alternatives using different evaluation criteria, and an Android application for collection of field data using mobile and tablet devices . In summary, the simulations of two- and three-blade wind turbines suggested that two-bladed machines could produce comparable annual energy as the three-blade wind turbines but have a lighter tower top weight, which leads to lower cost of energy. In addition, the two-blade rotor configuration potentially costs 20% less than a three blade configuration that produces the same power at the same site. The cost model analysis predicted a potential cost savings of approximately 15% for offshore two-blade wind turbines. The foundation design for a wind turbine in Lake Erie is likely to be driven by ice loads based on the currently available ice data and ice mechanics models. Hence, for Lake Eire, the cost savings will be somewhat smaller than the other lakes in the Great Lakes. Considering the size of cranes and vessels currently available in the Great Lakes, the cost optimal wind turbine size should be 3 MW, not larger. The surveillance data from different monitoring systems suggested that bird and bat passage rates per hour were comparable during heavy migrations in both spring and fall seasons while passage rates were significantly correlated to wind directions and wind speeds. The altitude of migration was higher during heavy migrations and higher over water relative to over land. Notable portions of migration on some spring nights occurred parallel the shoreline, often moving perpendicular to southern winds. The birds approaching the Western basin have a higher propensity to cross than birds approaching the Central basin of Lake Erie and as such offshore turbine development might be a better option further east towards Cleveland than in the Western basin. The high stopover density was more strongly associated with migration volume the following night rather than the preceding night. The processed mean scalar wind speeds with temporal resolutions as fine as 10-minute intervals near turbine height showed that August is the month with the weakest winds while December is the month, which typically has the strongest winds. The ice data suggests that shallow western basin of Lake Erie has higher ice cover duration many times exceeding 90 days during some winters.

17 WIND ENERGY↗