Search NASA⌕ Search

SEARCH · Search NASA

Results for “distributed surface roughness”

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.

251 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

A multi criteria analog model for assessing the vulnerability of rural catchments to road spills of hazardous substances

Road spills of hazardous substances are common in developing countries due to increasing industrialization and traffic accidents, and represent a serious threat to soils and water in catchments. There is abundant literature on equations describing the wash-off of pollutants from roads during a storm event and there are a number of watershed models incorporating those equations in storm water quality algorithms that route runoff and pollution yields through a drainage system towards the catchment outlet. However, methods describing catchment vulnerability to contamination by road spills based solely on biophysical parameters are scarce. These methods could be particularly attractive to managers because they can operate with a limited amount of easily collectable data, while still being able to provide important insights on the areas more prone to contamination within the studied watershed. The purpose of this paper was then to contribute with a new vulnerability model. To accomplish the goal, a selection of medium properties appearing in wash-off equations and routing algorithms were assembled and processed in a parametric framework based on multi criteria analysis to define the watershed vulnerability. However, parameters had to be adapted because wash-off equations and water quality models have been developed to operate primarily in the urban environment while the vulnerability model is meant to run in rural watersheds. The selected parameters were hillside slope, ground roughness (depending on land use), soil permeability (depending on soil type), distance to water courses and stream density. The vulnerability model is a spatially distributed algorithm that was prepared to run under the IDRISI Selva software, a GIS platform capable of handling spatial and alphanumeric data and execute the necessary terrain model, hydrographic and thematic analyses. For illustrative purposes, the vulnerability model was applied to the legally protected Environmental Protection Area (APA), located in the Uberaba region, state of Minas Gerais, Brazil. In this region, the risk of accidents causing chemical spills is preoccupying because large quantities of dangerous materials are transported in two important distribution highways while the APA is fundamental for the protection of water resources, the riverine ecosystems and remnants of native vegetation. In some tested scenarios, model results show 60% of vulnerable areas within the studied area. The most sensitive parameter to vulnerability is soil type. To prevent soils from contamination, specific measures were proposed involving minimization of land use conflicts that would presumably raise the soil's organic matter and in the sequel restore the soil's structural functions. Additionally, the present study proposed the preservation and reinforcement of riparian forests as one measure to protect the quality of surface water. - Highlights: • A multi criteria analog model was developed to assess rural catchment vulnerability along roads. • Model parameters were defined by analogy with urban wash-off equations and routing algorithms. • The model mixes up various biophysical and socio-economic parameters. • Model application was based on a scenario analysis. • The study is focused on the Environmental Protection Area of Uberaba River, Brazil.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Quantification of Off-Channel Inundated Habitat for Pacific Chinook Salmon ( Oncorhynchus tshawytscha ) along the Sacramento River, California, Using Remote Sensing Imagery

Off-channel areas are one of the most impacted aquatic habitats by humans globally, as extensive agricultural and urban development has limited them to roughly 10% of historical extent. This is also true for California’s Sacramento River Valley, where historically frequent widespread inundation has been reduced to a few off-channel water bodies along the mid-Sacramento River. This remaining shallow-water habitat provides crucial ecological benefits to multiple avian and fish species, but especially to floodplain-adapted species such as Chinook salmon (Oncorhynchus tshawytscha). Characterizing spatiotemporal off-channel dynamics, including inundation extent and residence time, is fundamental to better understanding the intrinsic value of such habitats and their potential to support recovery actions. Remote sensing techniques have been increasingly used to map surface water at regional and local scales, with improved resolutions. As such, this study maps off-channel inundation areas and describes their temporal dynamics by analyzing pixel-based time- series of multiple water indices, modified Normalized Difference Water Index (mNDWI) and the Automated Water Extraction Index (AWEI), generated from LandSat-8 and Sentinel-2 data between 2013–2021. Quantified off-channel area was similar with each water index and method used, but improved performance was associated with Sentinel-2 products and AWEI index to identify wetted areas under lower mainstem discharges. Results indicate an uneven distribution of off-channel habitat in the study area, with limited inundated areas in upstream reaches (<16% of total off-channel area for greater flows). In addition, much less habitat exists for flows under 400 m 3 /s, an important migration cue for endangered winter-run Chinook salmon, limiting juvenile access to areas with enhanced rearing conditions. Off-channel habitat residence times averaged between 7 and 16 days, primarily defined by the rate of receding flows, with rapid flow recession providing marginal off-channel habitat. This study shows reasonable performance of moderate resolution LandSat-8 and Sentinel-2 remote sensing imagery to characterize shallow-water inundated habitat in higher-order rivers, and as a method to inform restoration and native fish recovery efforts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Impact of induced seismic events on seal integrity, Texas Gulf Coast

Recent publications have suggested that large-scale CO 2 injection could trigger earthquakes and that even small- to moderate-sized earthquakes may threaten the seal integrity of the injection zone, and potentially damage buildings and other surface structures. In this study, we compared seal thickness to estimated fault displacement due to a single hypothetical seismic event in a selected area of the Texas Gulf Coast comprising an offshore strip of state waters along two Texas counties. To evaluate the slip generated by a single seismic event, we compiled well log information on shale/sand sequences and seismic information on fault geometric characteristics of a section of Lower Miocene age. The section is thousands of feet thick and is overlain and underlain by marine shales (Amph. B and Anahuac, respectively) that are relatively easy to correlate between wells. The Amph. B. shale is the secondary and ultimate seal for all injection intervals in the Lower Miocene. Given its thickness, no realistic seismic event or small series of seismic events will offset it significantly. However, this may not be true of smaller local primary seals. An analysis of geophysical logs of a total of 71 wells yielded a total of 2,871 sand / shale binary intervals. An analysis of the dedicated 3D seismic survey counted 723 fault traces at five roughly horizontal horizons within the Lower Miocene Fault displacement estimated using the product of the fault length times an uncertain multiplier coefficient assumed to follow a triangular distribution with a 10 -3 to 10 -5 range and a mode of 8 × 10 -5 . We then compared estimated single-event fault displacements to seal thicknesses by means of a Monte-Carlo analysis. Only 1.8% of thickness/displacement pairs display a displacement greater than 20% of the seal thickness. Only 0.26% of the pairs result in a displacement of half the seal thickness and only 0.05% of thickness/displacement pairs result in a clear seal rupture. The next step was to compare the magnitude of the event generated by such a displacement to documented magnitudes of “large” earthquakes generated by waterflooding and fluid disposal. Based on this analysis, we conclude that seismicity that may arise from CO 2 injection appears not to be a serious complication for CO 2 storage integrity, at least in the Gulf Coast area.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Early Failure of Lithium–Sulfur Batteries at Practical Conditions: Crosstalk between Sulfur Cathode and Lithium Anode

Lithium–sulfur (Li–S) batteries are one of the most promising next-generation energy storage technologies due to their high theoretical energy and low cost. However, Li–S cells with practically high energy still suffer from a very limited cycle life with reasons which remain unclear. Here, through cell study under practical conditions, it is proved that an internal short circuit (ISC) is a root cause of early cell failure and is ascribed to the crosstalk between the S cathode and Li anode. The cathode topography affects S reactions through influencing the local resistance and electrolyte distribution, particularly under lean electrolyte conditions. The inhomogeneous reactions of S cathodes are easily mirrored by the Li anodes, resulting in exaggerated localized Li plating/stripping, Li filament formation, and eventually cell ISC. Manipulating cathode topography is proven effective to extend the cell cycle life under practical conditions. The findings of this work shed new light on the electrode design for extending cycle life of high-energy Li–S cells, which are also applicable for other rechargeable Li or metal batteries.

25 ENERGY STORAGE↗
DOE OSTI2024

Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins

This repository of global hydrogeologic datasets contains aquifer properties on 0.5° scale, including depth to groundwater (Fan et al., 2013), aquifer thickness (de Graaf et al., 2015), WHYMap aquifer classes (Richts et al., 2011), porosity and permeability (Gleeson et al., 2014), digitized and geo-processed from their respective sources. Globally gridded aquifer properties could be used independently to estimate global groundwater availability or used as critical inputs to the superwell model to simulate groundwater extraction and provide estimates of pumped volumes and unit costs under user-specific scenarios. Key resources related to this data are: Niazi, H., Ferencz, S., Graham, N., Yoon, J., Wild, T., Hejazi, M., Watson, D., & Vernon, C. (2024; In-prep). Long-term Hydro-economic Assessment Tool for Evaluating Global Groundwater Cost and Supply: Superwell v1. Geoscientific Model Development. superwell model repository which uses this data to simulate groundwater extraction and provides estimates of the global extractable volumes and unit-costs ($/km3) of accessible groundwater production under user-specified extraction scenarios. Repository Overview Main output: aquifer_properties.csv contains all processed outputs, including aquifer properties like porosity, permeability, aquifer thickness, and depth to groundwater. shapefiles.zip: contains all digitized GIS databases and shapefile for all aquifer properties prep_inputs.R: R script that processes the shapefiles to produce the aquifer_properties file plot_inputs.R: R script for plotting the maps and conducting preliminary analysis on the available groundwater volume basin_to_country_mapping.csv, basin_country_region_mapping.csv and continent_county_mapping.csv provide the mapping between continents, 32 energy-economic macro regions, countries, and water basins for post-processing aquifer_properties.csv Maps: Each map visualizes the spatial distribution of one of the aquifer properties across the globe map_in_Porosity.png map_in_Permeability.png map_in_Aquifer_thickness.png map_in_Depth_to_water.png map_in_Grid_area_km.png map_in_WHYClass.png Sample inputs sample_inputs.py: this script samples inputs from the aquifer_properties dataset, ensuring the sampled and original inputs maintain the same distributions sampled_data_100.csv contains 100 sampled data points and sampled_data_100.png compares their distributions Dataset Overview The main outputs are consolidated in a comprehensive aquifer_properties.csv file and include the following fields: GridCellID: Unique identifier for each (roughly 0.5°) grid cell Continent: Continent name Country: Country name GCAM_basin_ID: Identifier for GCAM hydrologic basin Basin_long_name: Full name of the basin WHYClass: Hydrogeologic classification based on WHYMap aquifer classes (Richts et al., 2011) Porosity: Soil porosity (%) (Gleeson et al., 2014) Permeability: Soil permeability (in square meters; Gleeson et al., 2014) Aquifer_thickness: Thickness of the aquifer (in meters; de Graaf et al., 2015) Depth_to_water: Depth to groundwater (in meters; Fan et al., 2013) Grid_area: Area of the grid cell (in square meters) Key References The datasets are digitized versions of global hydrogeologic properties from the following key literature sources: Depth to Groundwater: Fan, Y., Li, H., & Miguez-Macho, G. (2013). Global Patterns of Groundwater Table Depth. Science, 339(6122), 940-943. https://doi.org/10.1126/science.1229881 Aquifer Thickness: de Graaf, I. E. M., Sutanudjaja, E. H., van Beek, L. P. H., & Bierkens, M. F. P. (2015). A high-resolution global-scale groundwater model. Hydrol. Earth Syst. Sci., 19(2), 823-837. https://doi.org/10.5194/hess-19-823-2015 Porosity and Permeability: Gleeson, T., Moosdorf, N., Hartmann, J., & van Beek, L. P. H. (2014). A glimpse beneath earth's surface: GLobal HYdrogeology MaPS (GLHYMPS) of permeability and porosity. Geophysical Research Letters, 41(11), 3891-3898. https://doi.org/10.1002/2014GL059856 Aquifer classes: Richts, A., Struckmeier, W. F., & Zaepke, M. (2011). WHYMAP and the Groundwater Resources Map of the World 1:25,000,000. In J. A. A. Jones (Ed.), Sustaining Groundwater Resources: A Critical Element in the Global Water Crisis (pp. 159-173). Springer Netherlands. https://doi.org/10.1007/978-90-481-3426-7_10 Cite as Niazi, H., Watson, D., Hejazi, M., Yonkofski, C., Ferencz, S., Vernon, C., Graham, N., Wild, T., & Yoon, J. (2024). Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins. MSD-LIVE Data Repository. https://doi.org/10.57931/2307831 Contact Reach out to Hassan Niazi or open an issue in superwell repository for questions or suggestions.

aquifer thickness↗
DOE OSTI2016

Additive Manufacturing of Advanced High Temperature Masking Fixtures for EBPVD TBC Coating

The purpose of this Manufacturing Demonstration Facility (MDF) technical collaboration project between Praxair Surface Technologies, Inc. (PST) and Oak Ridge National Laboratory (ORNL) was to develop an additive manufacturing process to fabricate next generation high temperature masking fixtures for coating of turbine airfoils with ceramic Thermal Barrier Coatings (TBC) by the Electron Beam Physical Vapor Deposition (EBPVD) process. Typical masking fixtures are sophisticated designs and require complex part manipulation in order to achieve the desired coating distribution. Fixtures are typically fabricated from high temperature nickel (Ni) based superalloys. The fixtures are fabricated from conventional processes by welding of thin sheet material into a complex geometry, to decrease the weight load for the manipulator and to reduce the thermal mass of the fixture. Recent attempts have been made in order to fabricate the fixtures through casting, but thin walled sections are difficult to cast and have high scrap rates. This project focused on understanding the potential for fabricating high temperature Ni based superalloy fixtures through additive manufacturing. Two different deposition processes; electron beam melting (EBM) and laser powder bed fusion were evaluated to determine the ideal processing route of these materials. Two different high temperature materials were evaluated. The high temperature materials evaluated were Inconel 718 and another Ni base alloy, designated throughout the remainder of this document as Alloy X, as the alloy composition is sensitive. Inconel 718 is a more widely utilized material for additive manufacturing although it is not currently the material utilized for current fixtures. Alloy X is the alloy currently used for the fixtures, but is not a commercially available alloy for additive manufacturing. Praxair determined it was possible to build the fixture using laser powder bed technology from Inconel 718. ORNL fabricated the fixture geometry using the EBM technology in order to compare deposition features such as surface roughness, geometric accuracy, deposition rate, surface and subsurface porosity, and material quality. It was determined that the laser powder bed technology was ideal for the geometry and requirements of the fixture set by Praxair, and Praxair moved forward with the purchase of a laser powder bed system. The subsequent portion of the project focused on determining the ideal processing parameters for alloy X for the laser powder bed system using ORNL’s Renishaw laser powder bed system. Praxair supplied gas atomized powders of alloy X material with properties specified by ORNL. ORNL printed text cube arrays in order to determine the ideal combination of laser powder and laser travel speed in order to maximize material density, improve surface quality, and maintain geometric accuracy. Additional powder supplied by Praxair was used to fabricate a full-scale fixture component.

42 ENGINEERING↗
DOE OSTI2022

In-situ Characterizations of Microstructural Degradation of Perovskite Solar Cells

Rapid progress has been achieved in perovskite solar cells (PSCs), and their efficiencies have improved from 3.8 % to 24.2 % in less than a decade. With low-cost processing, PSCs have shown exciting photovoltaic (PV) properties, such as effective optical absorption, a long carrier lifetime, and unique defect tolerance. While recent studies demonstrated improved stability up to 100 days, PSC technology is still challenged to meet the stringent industry requirements for commercialization. Despite considerable efforts, the underlying physical mechanisms for the inferior stability of PSCs are not well understood. One reason for this divergence is that many established measurement techniques (e.g., quantum efficiency, photoluminescence) probe the properties on length scales far greater than that of electronic and/or structural inhomogeneity (i.e., < 1 μm near grain boundaries) and therefore characterize convoluted and/or averaged properties. Ion/electron beam-based techniques have been extensively used to access the microstructures of PSCs, enabling atomic/nanoscale structural, chemical, optical, and electrical characterizations. For example, focused ion beam (FIB) milling produces an atomically smooth surface that minimizes the artifacts attributed to the surface roughness. FIB techniques can also create a well-defined cross-section of PSCs without mechanical damage in a physical cleaving sample preparation. While powerful, there are some concerns about possible beam damage of inorganic-organic perovskites via chemical-bond breakage and local heating. This project aims to comprehensively understand how the microstructural/interfacial properties of PSCs (e.g., Methylammonium Lead Iodide [MAPbI 3 ]) are modified under the irradiating ion beams. Specifically, we investigate the sub-surface properties of PSCs before and after Ar-ion beam injections. Kelvin probe force microscopy (KPFM) measures the contact potential differences (CPDs). Photoluminescence (PL) microscopy in conjunction with Finite-Difference Time-Domain (FDTD) simulations infers the formation of a “dead layer” (< 15 nm) on the subsurface of MAPbI 3 during Ar+ milling processes while preserving the initial bulk properties. The x-ray photoemission spectroscopy (XPS) confirms this modified surface is a lead-rich and iodine-deficient surface. We initiate customizing in-situ measurement setup while measuring the local optical and electrical properties of PSC under thermal (cooling, heating) and light stressors. Our results provide in-depth knowledge of the ion-beam impact on metal-halide perovskites and how this modified sub-surface impacts their properties under accelerated stressors of light and heat. Intensive Monte Carlo simulations of an electron beam interacting with PSCs provide the beam energy distribution in PSCs, proposing possible measurement conditions of using e-beam with minimizing beam damage. Our in-situ measurement platform can accommodate the diverse architecture of PSC devices for studying deterioration mechanisms under mixed environmental stressors.

14 SOLAR ENERGY↗
DOE OSTI2019

Accounting for hydrological and microbial processes on greenhouse gas budgets from river systems (Final Report)

River systems have a disproportionately large impact on carbon and nitrogen cycling on this planet. For instance, despite their much smaller size in area, it is estimated that global carbon dioxide emission from river systems exceed the emission from lakes and reservoirs, and are on par with the estimated oceanic and land carbon sinks. Closing carbon and nitrogen budgets in river systems requires knowledge of end member greenhouse gas (GHG) products (e.g. carbon dioxide, nitrous oxide, methane), which typically have not been measured in river systems. Moreover, while the production of these gasses is thought to be microbially derived, little is known about microbial activities in sediments of large rivers, despite their potentially high influence on biogeochemical budgets. Typically, when we think of rivers we visualize the channels that carve through the earth, but this only represents a small part of the hydrologic exchange that a river participates in. Rivers are more than conduits; they are dynamically exchanging water and materials with surrounding terrestrial and subsurface environments. In fact, 96% of the carbon dioxide produced within river systems occurs in hyporheic zones, or sediment regions beneath and alongside a river where groundwater and river water mix. In fact, it is widely recognized that microbial processes in hyporheic zones, also known as the river’s liver, strongly influence the fate of carbon, nutrients, and contaminants transported through watersheds. Despite this general knowledge, observational studies linking nitrogen and carbon gas production to microbial activity, and studies of overall gas emissions, are limited across river systems, but especially absent from larger rivers. To identify microbial processes responsible for GHG emissions from large rivers, and to quantify their environmental drivers and dependencies, we selected the Columbia River as a representative field site. This is the largest river west of the North American Continental Divide, and produces a drainage and discharge roughly double of the Nile in Egypt. This river is also an important source for agricultural irrigation and generates more hydroelectric power than any other North American river. Dams manage the release of water in the river, with these changes in river discharge patterns strongly impacting microbial activities in hyporheic zones. Specifically, we target a region down-river from the Priest Rapids Dam (WA) associated with the PNNL SFA, where water discharge from the dam is highly controlled, resulting in 1.5 m fluctuations in river stage daily. This dynamic river stage creates an ideal research opportunity for investigating hydrological perturbations in GHG production, uptake, and emission from large river systems. In this exploratory 1-year project, we conducted data-driven field research to provide observations needed to support predictions about GHG flux from river ecosystems. Using porewater dialysis and surface flux measurements, we observed where and when GHGs were produced and is the emission impacted by river stage. Using meta-omics tools, we determined what microbial processes control production and consumption of GHG in these systems. Finally, we provided key observations to parameterize models that could predict the effect of environmental conditions on GHG emission in the Columbia River near shore environment, at the scale that the modeled processes operate. Closing knowledge gaps in the linkages between river stage, microbial processes, and GHG production and emission, as well as the spatial distribution and ecosystem dependencies of these processes, is necessary to better parameterize land surface models. We provide some of the first, detailed GHG measurements paired microbial activity from large river systems. Our findings have direct relevance to many other large rivers with flows altered by hydroelectric energy, reservoirs, snowmelt, tides, or agriculture.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

THE APPLICATION OF MAMMOTH FOR A DETAILED TIGHTLY COUPLED FUEL PIN SIMULATION WITH A STATION BLACKOUT

Accurate calculation of desired quantities to predict fuel behavior requires the solution of interlinked equations representing different physics. Traditional fuels performance codes often rely on internal empirical models for the pin power density and a simplified boundary condition on the cladding edge. These simplifications are performed because of the difficulty of coupling applications or codes on differing domains and mapping the required data. To demonstrate an approach closer to first principles, the neutronics application Rattlesnake and the thermal hydraulics application RELAP-7 were coupled to the fuels performance application BISON under the master application MAMMOTH. A single fuel pin was modeled based on the dimensions of a Westinghouse 17x17 fuel rod. The simulation consisted of a depletion period of 1343 days, roughly equal to three full operating cycles, followed by a station blackout (SBO) event. The fuel rod was depleted for 1343 days for a near constant total power loading of 65.81 kW. After 1343 days the fission power was reduced to zero (simulating a reactor shut-down). Decay heat calculations provided the time-varying energy source after this time. For this problem, Rattlesnake, BISON, and RELAP-7 are coupled under MAMMOTH in a split operator approach. Each system solves its physics on a separate mesh and, for RELAP-7 and BISON, on only a subset of the full problem domain. Rattlesnake solves the neutronics over the whole domain that includes the fuel, cladding, gaps, water, and top and bottom rod holders. Here BISON is applied to the fuel and cladding with a 2D axi-symmetric domain, and RELAP-7 is applied to the flow of the circular outer water channel with a set of 1D flow equations. The mesh on the Rattlesnake side can either be 3D (for low order transport) or 2D (for diffusion). BISON has a matching ring structure mesh for the fuel so both the power density and local burn up are copied accurately from Rattlesnake. At each depletion time step, Rattlesnake calculates a power density, fission density rate, burn-up distribution and fast flux based on the current water density and fuel temperature. These are then mapped to the BISON mesh for a fuels performance solve. BISON calculates the fuel temperature and cladding surface temperature based upon the current power density and bulk fluid temperature. RELAP-7 then calculates the fluid temperature, water density fraction and water phase velocity based upon the cladding surface temperature. The fuel temperature and the fluid density are then passed back to Rattlesnake for another neutronics calculation. Six Picard or fixed-point style iterations are preformed in this manner to obtain consistent tightly coupled and stable results. For this paper a set of results from the detailed calculation are provided for both during depletion and the SBO event. We demonstrate that a detailed calculation closer to first principles can be done under MAMMOTH between different applications on differing domains.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Urban wash-off of tire wear particles

Tire wear particles (TWPs) are an important class of microplastics due to their toxicity and abundance. Because most TWPs are generated on impervious road surfaces, urban wash-off is the critical first phase of waterborne transport from their zone of production to stormwater drainage. However, little is known about the driving factors behind their mobilization. In this study, we use a rainfall simulator to investigate how surface roughness, rainfall intensity, and surface slope affect wash-off behaviors of TWPs. We also analyze how the size and shape of mobilized TWPs change over the course of simulated storm events. We found that low surface roughness, high rainfall intensity (most significant factor), and low slope result in the most rapid conveyance of TWP load. On average, large particles (>1000 µm) travelled faster than small particles (<125 µm). Particle shape explained a very small amount of variance in TWP wash-off velocity but was found to be more important under higher surface roughness conditions. In addition to wash-off velocity, we found similar conditions controlled the percent mobilization of TWPs. Low surface roughness and high rainfall intensity resulting in higher TWP wash-off rates is consistent with mineral sediment wash-off behavior. Conversely, low surface slope and large particle size leading to faster conveyance is directly opposed to mineral sediment wash-off. Our findings suggest drag-dominated flow and that sufficient runoff depth is the most important parameter governing TWP wash-off. These findings are important first steps to understanding wash-off behaviors of TWPs and informing future modeling efforts and mitigation strategies.

13 HYDRO ENERGY↗
DOE OSTI2015

The effect of dry and wet deposition of condensable vapors on secondary organic aerosols concentrations over the continental US

The effect of dry and wet deposition of semi-volatile organic compounds (SVOCs) in the gas phase on the concentrations of secondary organic aerosol (SOA) is reassessed using recently derived water solubility information. The water solubility of SVOCs was implemented as a function of their volatility distribution within the WRF-Chem regional chemistry transport model, and simulations were carried out over the continental United States for the year 2010. Results show that including dry and wet removal of gas-phase SVOCs reduces annual average surface concentrations of anthropogenic and biogenic SOA by 48 and 63% respectively over the continental US. Dry deposition of gas-phase SVOCs is found to be more effective than wet deposition in reducing SOA concentrations (−40 vs. −8% for anthropogenics, and −52 vs. −11% for biogenics). Reductions for biogenic SOA are found to be higher due to the higher water solubility of biogenic SVOCs. The majority of the total mass of SVOC + SOA is actually deposited via the gas phase (61% for anthropogenics and 76% for biogenics). Results are sensitive to assumptions made in the dry deposition scheme, but gas-phase deposition of SVOCs remains crucial even under conservative estimates. Considering reactivity of gas-phase SVOCs in the dry deposition scheme was found to be negligible. Further sensitivity studies where we reduce the volatility of organic matter show that consideration of gas-phase SVOC removal still reduces average SOA concentrations by 31% on average. We consider this a lower bound for the effect of gas-phase SVOC removal on SOA concentrations. A saturation effect is observed for Henry's law constants above 10 8 M atm −1 , suggesting an upper bound of reductions in surface level SOA concentrations by 60% through removal of gas-phase SVOCs. Other models that do not consider dry and wet removal of gas-phase SVOCs would hence overestimate SOA concentrations by roughly 50%. Assumptions about the water solubility of SVOCs made in some current modeling systems ( H * = H * (CH 3 COOH); H * = 10 5 M atm −1 ; H * = H * (HNO 3 )) still lead to an overestimation of 35%/25%/10% compared to our best estimate.

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

Synthesis and Characterization of Pd-based Nanomaterials

Bimetallic nanoparticles (BNPs) consist of two different types of metals or alloys that are bonded together. Unique properties such as optical, electronic, thermal, and catalytic effects differ for each type of BNP. Important BNPs range from Au-Pd, Ag- Pt, Au-Pt, and Ag-Ni. Pd bimetallic nanoparticles are of interest due to their many applications such as catalysis and sensing. Bimetallic catalysts have increase reaction rates and have improved catalyst stability through the geometry and ligand distribution. Pd nanoparticles are considered to be a strong catalyst due to their high activity at low temperatures and high tolerance to moisture. The catalytic properties of bimetallic nanoparticles depend on the structural properties such as size and shape. Core-shell, hollow structure, and multi-shell alloy are three possible structures nanoparticles can form as bimetallic catalysts. BNPs can be synthesized through different methods to control the size, shape, and structure. To obtain different morphologies, a variety of methods can be performed. Different methods can range from the usage of the glancing angle deposition (GLAD) to the galvanic replacement reaction, but the methods all depend on the properties of the metals. The galvanic displacement reaction was the method used to obtain Pd-based nanoparticles. This reaction is best know for obtaining hollow shaped NPs. To determine what redox process was preformed, the activity series of metals was used. From the activity series of metals, silver (Ag) was selected to preform Pd-based nanoparticles. Objectives: Synthesize Ag nanoparticles and Ag-Pd nanoparticles to understand the morphology. Characterize the synthesized nanoparticles using scanning electron microscopy (SEM), phase analysis light scattering (PALS), dynamic light scattering (DLS), energy dispersive X-ray spectroscopy (EDS), and UV-Vis spectroscopy. Results: In the UV-Vis spectrum, the Ag-Pd bimetallic NP's plasmon band decreased as the volume of palladium increased. The surface charge increases as the concentration of palladium increases. The Pd{sup 2+} ions interact with the sodium citrate surface, and decrease the negative charge. Conclusion: Ag-Pd nanoparticles were successfully created and stabilized with sodium citrate. The addition of Pd decreased the prominent plasmon band of the Ag nanoparticles. The SEM analysis showed that Ag nanoparticles had a well-defined structure, while the Ag-Pd nanoparticles showed hollow and rough structure. The EDX analysis confirmed the presence of silver and palladium. This material can be used in many industrial and research fields such as organic synthesis, fuel cells, and environmental sensing and remediation.

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

Evaluating the Feasibility of Measuring and Analyzing Strain Caused by Changes in Barometric Pressure in the Vadose Zone

Changes in barometric pressure propagate into the subsurface where they can affect water level measurements in wells and cause vertical strain. Previous strain sensors were developed to measure vertical strain, but they were limited to measurements at a single depth, making it difficult to evaluate strains from a migrating pressure wave. Technology known as Coherence-length-gated Microwave Photonics Interferometry (CMPI) uses optical fiber sensors to measure strain at multiple locations. It has the potential to detect strain caused by variations in barometric pressure at multiple depths in the subsurface. The goal of this study is to evaluate the feasibility of using this technology to record the changes in vertical strain in the vadose zone from small fluctuations in air pressure. Propagation of these air pressure fluctuations is affected by hydrogeological properties such as water content and permeability. A second goal of this study is to evaluate how these properties affect the air pressure distribution, strain, and gas-phase diffusivity with a long-term goal of using strain to monitor the vadose zone. The focus of this study is on a suite of laboratory experiments that used a sand-filled column with an open head space. The CMPI fiber was packaged and embedded along the axis of the column to measure the strain at multiple locations while air pressure transducers were installed through the wall of the column. Small periodic fluctuations in air pressure, similar to barometric pressure fluctuations, were created in the head space of the column using an audio driven speaker driven at 4 Hz using a sinusoidal signal created by a function generator. Then, the fluctuations propagated along the column where they were measured by the transducers and recorded as functions of time. Initial tests were conducted using dry sand, but then additional tests were conducted after injecting water into the column, which changed the hydrologic and mechanical properties. Another experiment was conducted after creating a thin barrier to air flow, which was designed to simulate a thin layer of saturated ground during rainfall. A final experiment was conducted by changing the grain size of the material to create a heterogeneity in the upper half of the column. The suite of six experiments was conducted to highlight effects of variations in water content and heterogeneities on air pressure and vertical strain. Numerical simulations were also created to evaluate the laboratory data and to provide a baseline analysis for the results. The numerical simulations used the governing equations to linear poroelasticity and two-phase flow with boundary conditions representing the experiments. The conceptual model for strain caused by barometric pressure recognizes that a fluctuating barometric pressure causes fluctuations in the air pore pressure that decrease in amplitude and lag in time with increase depth. The pressure distribution results in two different loads that cause strain: 1.) a mechanical load as the barometric pressure acts on the ground surface; and 2.) pressure loading in the pore space. The results show that the total vertical strain is a contribution of both of these loads, and this is the strain that was observed in the laboratory. The lab experiments show that both air pressure and strain vary as sinusoids that propagate with depth. The amplitude of the air pressure and strain decrease as a function of depth in all the conditions that were evaluate. For example, in the dry sand, the amplitude of the air pressure decreases from 25 Pa at the head space to 5 Pa at a depth of 70 cm, whereas the amplitude of the strain decreases from 0.15 to 0.06 µε over the same depth interval. In many cases the amplitudes decrease is approximately negative exponential functions of depth. When water was injected into the sand and the saturation increased, the amplitude of the air pressure increased in the partially saturated sand but decreased sharply to zero where the soil was saturated. The strain decreased with depth, but the transition from partially to fully saturated conditions had little effect on the strain profile, even though it had a major effect on the air pressure. This apparently is because strain is caused by variations in both water and air pressure. The fluctuating air pressure caused the water pressure to fluctuate, which caused strain throughout the column. Phase delay of the air pressure and the strain both increase as linear to bi-linear functions of depth. The phase delays indicate that pressure propagates at a velocity of 7 to 8 m/s in dry or partially saturated sand, but it drops by more than an order of magnitude in the vicinity of the capillary fringe. The phase delays indicate the strain propagates faster than the pressure. The data indicate that the strain propagates from 1.2 to more than 4x faster than the pressure. Switching from sand to silt in the column caused the velocity of the pressure to drop by roughly half for both pressures, but curiously, it had little effect on the velocity of the strain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Carbon Dioxide Shuttling Thermochemical Storage Using Strontium Carbonate

Phase I concludes with significant progress made towards the SunShot ELEMENTS goals of high energy density, high power density, and high temperature by virtue of a SrO/SrCO 3 based material. A detailed exploration of sintering inhibitors has been conducted and relatively stable materials supported by YSZ or SrZO 3 have been identified as the leading candidates. In 15 cycle runs using a 3 hour carbonation duration, several materials demonstrated energy densities of roughly 1500 MJ/m 3 or greater. The peak power density for the most productive materials consistently exceeded 40 MW/m 3 —an order of magnitude greater than the SOPO milestone. The team currently has a material demonstrating nearly 1000 MJ/m 3 after 100 abbreviated (1 hour carbonation) cycles. A subsequent 8 hour carbonation after the 100 cycle test exhibited over 1500 MJ/m 3 , which is evidence that the material still has capacity for high storage albeit with slower kinetics. Kinetic carbonation experiments have shown three distinct periods: induction, kinetically-controlled, and finally a diffusion-controlled period. In contrast to thermodynamic equilibrium prediction, higher carbonation temperatures lead to greater conversions over a 1 hour periods, as diffusion of CO 2 is more rapid at higher temperatures. A polynomial expression was fit to describe the temperature dependence of the linear kinetically-controlled regime, which does not obey a traditional Arrhenius relationship. Temperature and CO 2 partial pressure effects on the induction period were also investigated. The CO 2 partial pressure has a strong effect on the reaction progress at high temperatures but is insignificant at temperatures under 900°C. Tomography data for porous SrO/SrCO 3 structures at initial stage and after multiple carbonation/decomposition cycles have been obtained. Both 2D slices and 3D reconstructed representations have been obtained from the raw data. Porosity and surface area have been computed based on tomography. The two parameters, permeability and Forchheimer constants, in the Darcy’s law for linear and non-linear range have been obtained from the fluid dynamics simulation. The structure is nearly isotropic from the hydrodynamics point of view. Thermal conductivity at elevated temperature has been obtained from the 3-D pore-scale simulations. The computed thermal conductivities of the structure compare well with an existing correlation. Computational models have been developed for solving the energy and mass transport equations with chemical reaction on the surface and inside the particles. The chemical kinetics or rate law is needed on the particle scale; however it is generally measured using a lumped model based on the average of reacted as well as unreacted mass. Many simulations have been conducted on a 2-D model to determine the effects of those pore scale parameters on the heat release based on the measured data. The simulation results using 2-D model show that average volumetric heat flow rate of 4 MW/m 3 in the first 5 minutes can be obtained. The insight gained from the 2-D simulation will be useful for the next stage simulation based on realistic 3-D structures. A probabilistic analysis has been conducted on the baseline system of recovering heat from the hot CO 2 evolved from the endothermic reaction and compressing the cooled CO 2 into a storage tank at 20-21 bar. Using a Monte Carlo simulation using 2,000,000 samples and both uniform and triangular probability distributions of four variables (heat exchanger pinch points, compressor efficiency, and power block efficiency), exergetic efficiency of 90% appears to be very likely.

14 SOLAR ENERGY↗
DOE OSTI2025

Investigating Spatial Variability of Aerosol, Cloud Condensation Nuclei, and Ice Nucleating Particles in Mountainous Terrain Field Campaign Report

The U.S. Department of Energy Atmospheric System Research (ASR)-supported Surface Atmosphere Integrated Field Laboratory (SAIL) campaign in the East River Watershed (ERW) of the Upper Colorado River Basin in southwestern Colorado ran from fall 2021 to spring 2023. Two monitoring sites were deployed in the ERW as part of SAIL. The two sites were the Aerosol Observation System (AOS) located on Crested Butte Ski Mountain, and the second ARM Mobile Facility (AMF2), located at the Rocky Mountain Biological Laboratory in Gothic, Colorado. To gain a more comprehensive understanding of aerosols in complex, mountainous terrain, Handix Scientific deployed SAIL-Net, a distributed network of six measurement nodes spanning the domain of the SAIL research area from October 2021 to July 2023. Each node measured aerosol particles between 140 nm and 3.4 μm in diameter using a small portable optical particle spectrometer (POPS; Gao et al. 2016), cloud condensation nuclei (CNN) using a miniature CCN counter (CloudPuck), and ice nucleating particles (INP) using the time-resolved aerosol filter sampler (TRAPS; Creamean et al. 2018). Our approach was similar to other studies that aimed to better characterize and understand aerosols and gas-phase pollutants using networks of lower-cost sensors (Caubel et al. 2019, Kelly et al. 2021, Asher et al. 2022). Such studies have identified neighborhood-level variations in pollutant concentrations (Schneider et al. 2017, Popoola et al. 2018, Caubel et al. 2019). Small-scale variations such as this are poorly represented in models and poorly measured by a single monitoring system (Caubel et al. 2019). Previous work has shown the representation error (the ability of measurements to represent a larger area) increases with complex orography, leading to decreases in model accuracy (Schutgens et al. 2017). The overall goal of SAIL-Net was to improve our understanding of the variability of aerosol in the ERW, thus increasing our knowledge of aerosol-cloud interactions in this region and informing the usefulness of distributed networks of measurements for future studies. We met this goal by answering the following science questions: 1. What is the aerosol temporal variability, and how does aerosol inhomogeneity vary seasonally? Is there significant seasonal variability in sources, or are short-term meteorological conditions the most important determining factor in sources for cloud nuclei? 2. What is the aerosol spatial variability? What are the aerosol characteristics at cloud base, presumably the particles most representative of those acting as cloud nuclei? 3. How should measurement networks be designed to capture aerosol-cloud interactions, and what do they need to measure? Can a single measurement site accurately represent aerosol properties in regions of complex terrain? SAIL-Net consisted of six measurement nodes spread across the ERW near Crested Butte, Colorado. The primary objective in site placement was to select locations that captured the vertical variation in aerosol properties while also spanning the domain of the SAIL campaign. The elevation of the sites ranged from roughly 2750 m along the valley floor of the ERW to approximately 3500 m near the top of Crested Butte Mountain, which is one of the taller peaks in the ERW. The farthest distance between sites was 14 km, while the closest two sites were approximately 1 km apart. Two of the sites were collocated with the ARM SAIL sites; our instruments sat on top of one of the trailers at AOS and another one of our sites was located in a meadow just above AMF2.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Proposal and evaluation of subordinate standard solar irradiance spectra for applications in solar energy systems

Reference solar irradiance spectra are needed to specify key parameters of solar technologies such as photovoltaic cell efficiency, in a comparable way. The IEC 60904-3 and ASTM G173 standards present such spectra for Direct Normal Irradiance (DNI) and Global Tilted Irradiance (GTI) on a 37 degrees tilted sun-facing surface for one set of clear-sky conditions with an air mass of 1.5 and low aerosol content. The IEC/G173 standard spectra are the widely accepted references for these purposes. Hence, the authors support the future replacement of the outdated ISO 9845 spectra with the IEC spectra within the ongoing update of this ISO standard. The use of a single reference spectrum per component of irradiance is important for clarity when comparing and rating solar devices such as PV cells. However, at some locations the average spectra can differ strongly from those defined in the IEC/G173 standards due to widely different atmospheric conditions and collector tilt angles. Therefore, additional subordinate standard spectra for other atmospheric conditions and tilt angles are of interest for a rough comparison of product performance under representative field conditions, in addition to using the main standard spectrum for product certification under standard test conditions. This simplifies the product selection for solar power systems when a fully-detailed performance analysis is not feasible (e.g. small installations). Also, the effort for a detailed yield analyses can be reduced by decreasing the number of initial product options. After appropriate testing, this contribution suggests a number of additional spectra related to eight sets of atmospheric conditions and tilt angles that are currently considered within ASTM and ISO working groups. The additional spectra, called subordinate standard spectra, are motivated by significant spectral mismatches compared to the IEC/G173 spectra (up to 6.5%, for PV at 37 degrees tilt and 10-15% for CPV). These mismatches correspond to potential accuracy improvements for a quick estimation of the average efficiency by applying the appropriate subordinate standard spectrum instead of the IEC/G173 spectra. The applicability of these spectra for PV performance analyses is confirmed at five test sites, for which subordinate spectra could be intuitively selected based on the average atmospheric aerosol optical depth (AOD) and precipitable water vapor at those locations. The development of subordinate standard spectra for DNI and concentrating solar power (CSP) and concentrating PV (CPV) is also considered. However, it is found that many more sets of atmospheric conditions would be required to allow the intuitive selection of DNI spectra for the five test sites, due in particular to the stronger effect of AOD on DNI compared to GTI. In conclusion, the matrix of subordinate GTI spectra described in this paper are recommended to appear as an option in the annex of future standards, in addition to the obligatory use of the main spectrum from the ASTM G173 and IEC 60904 standards.

14 SOLAR ENERGY↗