Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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.

At least 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Recent Thermal and Hydraulic Stimulation Results at Raft River, ID EGS Site

The Raft River geothermal system, located 160 km northwest of Salt Lake City, Utah is the site of a Department of Energy Enhanced Geothermal System (EGS) stimulation demonstration project. The target well RRG-9 ST-1, is roughly 1.6 km (1 mile) south of the power plant. This well has undergone a series of thermal and hydraulic stimulation treatments. Furthermore, the well passes through approximately 1,524 m MD (5,000 ft. MD) of laterally discontinuous Quaternary and Tertiary volcanoclastic and volcanic rocks that unconformably overlie the Precambrian basement. The Elba Quartzite formation located in the Precambrian basement is the target zone for the stimulation program. Fluid injected into the well exits through a fracture zone in the Elba Quartzite at 1,567 m TVD (5,140 ft. TVD). In June 2013 low-rate stimulation began at RRG-9 ST-1 using injection water from the geothermal plant. The average injection rate, temperature, and well head pressure of the low-rate stimulation was 163 Lpm (43 gpm), 34°C, and 1,931kPa (280 psig), respectively. This injection continued through late August of 2013 when wellhead pressures were intentionally increased to approximately 5,516 kPa (800 psig) at approximately 984 Lpm (260 gpm). From September 12 to September 24, 2013 cold well water was injected. The average surface injection temperature was approximately 12°C. On September 25, 2013 injection of the warmer plant water was resumed. This was continued through the rest of 2013 and into 2014. Between April 1, 2014 and April 4, 2014, a higher rate hydraulic stimulation was conducted. Colder water (~10°C) was injected at rates of 20 bpm to 30 bpm (3,178 Lpm to 4,770 Lpm) at wellhead pressures of 5,861 kPa (850 psig) over approximately 24 hours. Microseismic activity was recorded near the Tertiary-Precambrian contact. Subsequently, injection of water from the plant was resumed. As of mid-January, 2015, over 666 million liters (176 million gallons) have been injected into the well. The fate of this water is the subject of ongoing numerical analyses and continuous site monitoring. Various analyses have been implemented to monitor and assess the ongoing stimulation program. Modified Hall plot analysis shows increasing conductivity and/or a decreasing skin factor around the wellbore. Since April, 2014, the injectivity index has increased significantly from 0.28 Lpm/kPa (0.5 gpm/psig) to 0.93 Lpm/kPa (1.7 gpm/psig).

enhanced geothermal system↗
DOE OSTI2023

Annual cycle of aerosol properties over the central Arctic during MOSAiC 2019–2020 – light-extinction, CCN, and INP levels from the boundary layer to the tropopause

Abstract. The MOSAiC (Multidisciplinary drifting Observatory for the Study of Arctic Climate) expedition was the largest Arctic field campaign ever conducted. MOSAiC offered the unique opportunity to monitor and characterize aerosols and clouds with high vertical resolution up to 30 km height at latitudes from 80 to 90∘ N over an entire year (October 2019 to September 2020). Without a clear knowledge of the complex aerosol layering, vertical structures, and dominant aerosol types and their impact on cloud formation, a full understanding of the meteorological processes in the Arctic, and thus advanced climate change research, is impossible. Widespread ground-based in situ observations in the Arctic are insufficient to provide these required aerosol and cloud data. In this article, a summary of our MOSAiC observations of tropospheric aerosol profiles with a state-of-the-art multiwavelength polarization Raman lidar aboard the icebreaker Polarstern is presented. Particle optical properties, i.e., light-extinction profiles and aerosol optical thickness (AOT), and estimates of cloud-relevant aerosol properties such as the number concentration of cloud condensation nuclei (CCN) and ice-nucleating particles (INPs) are discussed, separately for the lowest part of the troposphere (atmospheric boundary layer, ABL), within the lower free troposphere (around 2000 m height), and at the cirrus level close to the tropopause. In situ observations of the particle number concentration and INPs aboard Polarstern are included in the study. A strong decrease in the aerosol amount with height in winter and moderate vertical variations in summer were observed in terms of the particle extinction coefficient. The 532 nm light-extinction values dropped from >50 Mm−1 close to the surface to <5 Mm−1 at 4–6 km height in the winter months. Lofted, aged wildfire smoke layers caused a re-increase in the aerosol concentration towards the tropopause. In summer (June to August 2020), much lower particle extinction coefficients, frequently as low as 1–5 Mm−1, were observed in the ABL. Aerosol removal, controlled by in-cloud and below-cloud scavenging processes (widely suppressed in winter and very efficient in summer) in the lowermost 1–2 km of the atmosphere, seems to be the main reason for the strong differences between winter and summer aerosol conditions. A complete annual cycle of the AOT in the central Arctic could be measured. This is a valuable addition to the summertime observations with the sun photometers of the Arctic Aerosol Robotic Network (AERONET). In line with the pronounced annual cycle in the aerosol optical properties, typical CCN number concentrations (0.2 % supersaturation level) ranged from 50–500 cm−3 in winter to 10–100 cm−3 in summer in the ABL. In the lower free troposphere (at 2000 m), however, the CCN level was roughly constant throughout the year, with values mostly from 30 to 100 cm−3. A strong contrast between winter and summer was also given in terms of ABL INPs which control ice production in low-level clouds. While soil dust (from surrounding continents) is probably the main INP type during the autumn, winter, and spring months, local sea spray aerosol (with a biogenic aerosol component) seems to dominate the ice nucleation in the ABL during the summer months (June–August). The strong winter vs. summer contrast in the INP number concentration by roughly 2–3 orders of magnitude in the lower troposphere is, however, mainly caused by the strong cloud temperature contrast. A unique event of the MOSAiC expedition was the occurrence of a long-lasting wildfire smoke layer in the upper troposphere and lower stratosphere. Our observations suggest that the smoke particles frequently triggered cirrus formation close to the tropopause from October 2019 to May 2020.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Geranium Photodiodes for Hard Xray Detection

This report summarizes the investigations engaged at UC Davis under subcontract B632083, as related to the LLNL-LDRD Project “High-Speed X-Ray Imager Arrays.” The first section entails work done within the Electrical and Computer Engineering (ECE) Department. The purpose of the work done so far is to understand photodiode design and which parameters are conducive to the design goals, which are: I: Maximize the detector quantum efficiency. The device must have a quantum efficiency of roughly 5% at 80 keV. II: Minimize temporal response. Collection time is limited to 1ns. All charge must be collected in that time. III: Minimize dark current. In this document, various photodiode parameters and designs are explored, and their impact on these goals are tabulated. In addition to these, various instruments that have been constructed for either testing these parameters or helping teammates with their experiments. The second section entails work done with the Materials Science Engineering (MSE) Department. In accordance with the proposed development of the germanium ROIC backside imagers, a critical step in the micromachining of these devices was the removal of the substrate on which the PIN junctions were epitaxially grown. A common technique for the micromachining of silicon wafers is the use of chemical etching alongside an “etch-stop” layer. In this technique, a chemical etchant will remove material until it reaches the etch-stop layer, at which the etch rate will drop significantly. This high selectivity (ratio of etch rates) allows any unevenness in the initial etching to be effectively smoothed away in the event that part of the wafer surface reaches the etch-stop layer sooner than another. An initial literature review found a paper by Divan et al. in which the authors used ion implantation to dope germanium wafers and achieve selectivity on the order of 100, sufficient for a reliable etch-stop.

42 ENGINEERING↗
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 OSTI2025

Formation and fate of freshwater on an ice floe in the Central Arctic

The melt of snow and sea ice during the Arctic summer is a significant source of relatively fresh meltwater. The fate of this freshwater, whether in surface melt ponds or thin layers underneath the ice and in leads, impacts atmosphere–ice–ocean interactions and their subsequent coupled evolution. Here, we combine analyses of datasets from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition (June–July 2020) for a process study on the formation and fate of sea ice freshwater on ice floes in the Central Arctic. Our freshwater budget analyses suggest that a relatively high fraction (58 %) is derived from surface melt. Additionally, the contribution from stored precipitation (snowmelt) outweighs by 5 times the input from in situ summer precipitation (rain). The magnitude and rate of local meltwater production are remarkably similar to those observed on the prior Surface Heat Budget of the Arctic Ocean (SHEBA) campaign, where the cumulative summer freshwater production totaled around 1 m during both. A relatively small fraction (10 %) of freshwater from melt remains in ponds, which is higher on more deformed second-year ice (SYI) compared to first-year ice (FYI) later in the summer. Most meltwater drains laterally and vertically, with vertical drainage enabling storage of freshwater internally in the ice by freshening brine channels. In the upper ocean, freshwater can accumulate in transient meltwater layers on the order of 0.1 to 1 m thick in leads and under the ice. The presence of such layers substantially impacts the coupled system by reducing bottom melt and allowing false bottom growth; reducing heat, nutrient, and gas exchange; and influencing ecosystem productivity. Regardless, the majority fraction of freshwater from melt is inferred to be ultimately incorporated into the upper ocean (75 %) or stored internally in the ice (14 %). Terms such as the annual sea ice freshwater production and meltwater storage in ponds could be used in future work as diagnostics for global climate and process models. For example, the range of values from the CESM2 climate model roughly encapsulate the observed total freshwater production, while storage in melt ponds is underestimated by about 50 %, suggesting pond drainage terms as a key process for investigation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A diagrammatic approach for automatically deriving analytical gradients of tensor hyper-contracted electronic structure methods

In this work, we introduce a diagrammatic approach to facilitate the automatic derivation of analytical nuclear gradients for tensor hyper-contraction (THC) based electronic structure methods. The automatically derived gradients are guaranteed to have the same scaling in terms of both operation count and memory footprint as the underlying energy calculations, and the computation of a gradient is roughly three times as costly as the underlying energy. The new diagrammatic approach enables the first cubic scaling implementation of nuclear derivatives for THC tensors fitted in molecular orbital basis (MO-THC). Furthermore, application of this new approach to THC-MP2 analytical gradients leads to an implementation, which is at least four times faster than the previously reported, manually derived implementation. Finally, we apply the new approach to the 14 tensor contraction patterns appearing in the supporting subspace formulation of multireference perturbation theory, laying the foundation for developments of analytical nuclear gradients and nonadiabatic coupling vectors for multi-state CASPT2.

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

Calculation of Friction Factors in Heat Pipes Using CFD in Support of the Sockeye Code

A heat pipe is two-phase heat transfer device which relies on surface tension and capillary pressure to provide a very efficient heat transfer mechanism. Currently, there is interest in designing micro nuclear power reactors using high temperature heat pipes to provide totally passive cooling. To support this, Sockeye is being developed as a heat pipe simulation capability. To support the development of Sockeye, individual phenomena must be tested. This work focuses on applying a CFD model to evaluate friction factors for use in Sockeye. A CFD model has been developed in Nek5000 which simulates each component fluid, i.e. sodium liquid and vapor, independently and captures the effects of the wick on friction factor. The model tests the effect of the wick permeability and the thickness of the wick occupied by each fluid component. The thickness of the wick occupied by each component can be directly related to the local vapor volume fraction from the Sockeye heat pipe model. Results from the CFD model indicate that friction factor in the liquid flow decreases with increasing wick permeability, while the friction factor in the vapor flow increases with increasing wick permeability. In the liquid, this is likely a direct result of the increased cross sectional area available for the flow. For the vapor, this is explained as the wick acting like a rough wall, while the increased area has only a negligible impact. Results for the tests of the thickness of the wick occupied by each fluid showed a thresholding behavior for both liquid and vapor. When the thickness of the wick occupied by each fluid was increased, it increased the effect of the wick, but only up to a certain point. Once a certain thickness was reached, the friction factors remained constant. When the fluid occupied only a very small fraction of the wick, friction factors returned to nominal values.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
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), recharge (Döll and Fiedler, 2008; Gleeson et al., 2016), lakes (Messager et al., 2016), 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. B., Graham, N. T., Yoon, J., Wild, T. B., Hejazi, M., Watson, D. J., & Vernon, C. R. (2025). Long-term hydro-economic analysis tool for evaluating global groundwater cost and supply: Superwell v1.1. Geoscientific Model Development, 18(5), 1737-1767. https://doi.org/10.5194/gmd-18-1737-2025 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_rec.csv contains all processed outputs, including aquifer properties like porosity, permeability, recharge, lake areas, aquifer thickness, and depth to groundwater. shapefiles.zip: contains all digitized GIS databases and shapefile for all aquifer properties prep_inputs.R and prep_inputs_recharge_lakes.R: R scripts that process the shapefiles to produce the aquifer_properties_rec.csv 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_rec.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_Recharge.png map_in_Grid_area_km.png map_in_Lake_area_km.png map_in_WHYClass.png Sample inputs sample_inputs.py: this script samples inputs from the aquifer_properties_rec 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_rec.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) Recharge: long-term annual averaged recharge rates (in m/yr; Döll and Fiedler, 2008; Gleeson et al., 2016) Grid_area: Area of the grid cell (in square meters) Lakes_area: Area of inland lakes (in square meters; Messager et al., 2016) 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 Recharge: Döll, P., & Fiedler, K. (2008). Global-scale modeling of groundwater recharge. Hydrol. Earth Syst. Sci., 12(3), 863-885. https://doi.org/10.5194/hess-12-863-2008; Gleeson, T., Befus, K. M., Jasechko, S., Luijendijk, E., & Cardenas, M. B. (2016). The global volume and distribution of modern groundwater. Nature Geoscience, 9(2), 161-167. https://doi.org/10.1038/ngeo2590 Inland Lakes: Messager, M. L., Lehner, B., Grill, G., Nedeva, I., & Schmitt, O. (2016). Estimating the volume and age of water stored in global lakes using a geo-statistical approach. Nature Communications, 7(1), 13603. https://doi.org/10.1038/ncomms13603 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. MultiSector Dynamics-Living, Intuitive, Value-adding, Environment. https://doi.org/10.57931/2484226 Contact Reach out to Hassan Niazi or open an issue in the superwell repository for questions or suggestions.

aquifer thickness↗
DOE OSTI2023

Light-trapping structures fabricated in situ for ultrathin III-V solar cells

Here, we describe a fully in situ method of fabricating light-scattering structures on III-V materials that generates a rough morphology via vapor phase etching and redeposition. Fully in situ methods support higher industrial throughput by utilizing the growth reactor to generate the light-trapping structures after device growth without removal from the reactor. We use HCl and PH 3 to etch and redeposit scattering morphologies on Ga 0.5 In 0.5 P in a dynamic hydride vapor phase epitaxy (D-HVPE) reactor. We show that the addition of PH 3 leads to redeposition during the vapor phase HCl etching of Ga 0.5 In 0.5 P and that HCl flow rate and time exposed to HCl-PH 3 each independently cause a linear increase in the redeposited feature size, indicating that redeposition proceeds by island growth in a III-Cl-limited, hydride-enhanced HVPE regime. Auger electron spectroscopy and scanning transmission electron microscopy with energy dispersive spectroscopy (STEM-EDS) reveal redeposition to be highly Ga-rich GaInP, i.e., Ga(In)P. The Ga-rich nature of the redeposition results from the higher thermodynamic driving force for Ga incorporation than for In during HVPE growth and the difference in the volatility of the III-Cl etch products. The resulting morphologies have high broadband scattering, as determined by normal specular reflectance and integrating sphere measurements, indicating effectiveness as light-scattering structures. In a 270-nm-thick GaAs photovoltaic device with a textured back surface, we achieve a 4.9% increase in short circuit current density (J SC ) without any loss in open-circuit voltage (V OC ) relative to a planar control using only a 60 s in situ texturing treatment.

14 SOLAR ENERGY↗
DOE OSTI2021

Topical Report – Findings on Subtask 2.7 – Wet ESP and Aerosol Testing at Coal Creek Station

Growing concerns over the impact of CO 2 emissions from combustion sources on global climate change have prompted numerous research and development projects aimed at developing cost-effective technologies for CO 2 capture. One family of technologies being demonstrated at pilot and full scale globally is postcombustion carbon capture (PCCC) systems that employ amine-based solvents. The captured CO 2 can be compressed and permanently stored underground or used for enhanced oil recovery. The proximity of North Dakota’s lignite-fired fleet of power plants to potential CO 2 storage options creates a unique atmosphere for PCCC within the state. However, the unique components present in lignite flue gas present a challenge for large-scale PCCC at North Dakota power plants by contributing to aerosol formation. Aerosols can negatively impact the long-term performance of amine-based solvents for CO 2 capture. Amine-based solvents are volatile, and flue gas particulate provides nucleation sites where amine vapors can condense as aerosols. Because aerosols cannot be easily captured at the column outlet using conventional technologies, the amine-laden aerosols escape the system and lead to accelerated solvent losses. Moreover, particulate components can chemically react with amines to form degradation products that can permanently deactivate the amine, cause fouling, and lead to hazardous emissions. Many of the elements that have been shown to catalyze solvent degradation are present in lignite coals and can exacerbate solvent replacement economics. Understanding this issue is critical to the implementation of solvent-based CO 2 capture systems as applied to lignite-fired generation systems. The Energy & Environmental Research Center (EERC) designed and carried out this project to fully characterize aerosol behavior with various control technologies installed to better optimize aerosol mitigation technology for CO 2 capture. To meet the goal of this project, the following objectives were identified: Determine the effectiveness of a wet electrostatic precipitator (WESP) on mitigating formation of problematic aerosols at Great River Energy’s Coal Creek Station, upstream of the PCCC system. Determine the effectiveness of the Mitsubishi Heavy Industries (MHI) proprietary amine emission reduction unit (AERU) as a postcapture solvent recovery system for reducing aerosol emissions and extending solvent life downstream of the PCCC system. Determine the impact of aerosols on the efficiency and degradation products of both commercial and advanced solvents within the PCCC system. Work was conducted at Coal Creek Station Unit 1 using a slipstream of flue gas from the outlet of the plant’s flue gas desulfurization (FGD) unit. Flue gas was routed through a pilot-scale FGD unit to remove SO 2 to very low levels (~1 ppm) and then through a direct contact cooler (DCC) to further cool the gas and to remove moisture. The gas exiting the DCC was then optionally routed through a WESP before passing to the CO 2 absorber columns. The MHI solvent was used to scrub CO 2 from the slipstream through a set of two absorber columns. The rich solvent was regenerated in a stripper column by heating to drive off captured CO 2 . Flue gas exiting the absorber column was routed to MHI’s proprietary AERU to recover entrained solvent. The system operated using a catch-and-release method where the CO 2 was separated to provide data on the process, but the captured CO 2 was released back into the host site stack. Particulate was measured, collected, and analyzed from multiple locations throughout the pilot-scale system. Unlike the performance observed in prior work, the inlet FGD and DCC did not remove significant particulate matter from the flue gas. This appears to be due to a difference in the nature of the particulate. The DCC seemed to increase particulate size and count, most likely owing to water condensing onto the surfaces of fly ash particles. When the WESP was operated, it achieved >95% particulate capture. Very little particulate matter or indications of solvent were detected at the AERU outlet. When operating with advanced KS 21 solvent, the particulate material at the AERU outlet was even further decreased. Solvent analysis showed that some species derived from flue gas and ash were slowly concentrating in the solvent over the duration of the test. The levels observed were reported to be within expected ranges and were not of concern to MHI. A high-level techno-economic assessment of installing CO 2 capture at Coal Creek Station suggested that, when using a standard monoethanolamine (MEA)-based solution with simple heat integration, the energy penalty to net generation would be 34%. The bulk of this was due to steam losses for regenerating solvent, followed by parasitic electrical demand for CO 2 compression and then by increased parasitic load for pushing flue gas through the absorber column. These demands could be decreased with a more advanced solvent that exhibits lower heat of regeneration and lower pressure drop than does a simple MEA solution. Further energy could be saved with more thorough heat integration to recover useful energy from the steam used for solvent regeneration. Installing a WESP was predicted to increase the cost of electricity by nearly $5/MWh. This would become cost-effective if solvent losses were roughly 10 times the baseline estimate when not using a WESP but could be returned to baseline by installing the WESP. Piping CO 2 for storage in more favorable geology could help with carbon capture and storage economics. Although storing off-site would necessitate construction of a CO 2 transport pipeline, the cost of this pipeline might be more than offset by reducing the number of wells required, the depths of the wells required, and the electrical demand for the CO 2 compressor. Additional factors that favor off-site storage costs include smaller expected CO 2 plume sizes, which translates to less monitoring and fewer landowner agreements. More detailed assessment of the specific geology in the region under and around Coal Creek Station would be needed to accurately assess the costs and benefits of different storage site options. This subtask was cofunded through the Energy & Environmental Research Center–U.S. Department of Energy Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE-FE0024233. Nonfederal funding was provided by the North Dakota Industrial Commission.

Strege, Joshua↗
DOE OSTI2023

Onshore U.S. Carbon Pipeline Deployment: Siting, Safety, and Regulation

Carbon capture, utilization, and storage (CCUS) technology has significant potential to reduce greenhouse gas (GHG) emissions and mitigate the impact of climate change, particularly in hard to decarbonize industrial and commercial sectors. CCUS involves capturing carbon dioxide (CO 2 ) from industrial processes or power generation and utilizing it for other purposes, such as enhanced oil recovery (EOR), or storing the captured CO 2 underground. CCUS technology can reduce the environmental impact of continued fossil fuel use while smoothing the transition to a low-carbon economy. CCUS can create new economic opportunities, such as the development of new industries and job creation, and can enhance energy security by diversifying energy sources. For these reasons, enabling CCUS has become a key objective of the Biden-Harris administration’s clean energy policy and has received bipartisan support. Despite its environmental and economic potential, CCUS faces multiple barriers to widespread deployment. One of the main challenges is the high cost and technical difficulty of implementing and operating large-scale CCUS infrastructure. CCUS remains a relatively expensive way to reduce carbon emissions (e.g., compared to solar photovoltaic technology’s displacement of coal generation). Additionally, financial incentives and supportive policies like those enacted to support solar photovoltaic development, especially at the state level, are inconsistent or nonexistent, which can discourage investment in CCUS projects. There are also technical challenges associated with safe and secure underground CO 2 storage and the development of new carbon utilization technologies. Public opposition to various aspects of CCUS technologies, ranging from concerns that CCUS will extend reliance on fossil fuels to CCUS infrastructure being sited in disadvantaged communities, is a growing challenge. This paper focuses on another significant barrier to broad CCUS deployment: the need for considerable expansion of the dedicated land-based CO 2 pipeline network in the United States to meet CCUS goals and the unique regulatory challenges to its development. To reach carbon emissions targets in the United States by 2050, CCUS technology will need to be supported by tens of thousands of miles of CO 2 pipelines. Estimates range from a minimum of roughly 29,000 pipeline miles (according to a 2020 Great Plains Institute study) to 66,000 pipeline miles (as per a 2021 Princeton University–led study). As of October 2022, however, the U.S. Department of Transportation (U.S. DOT) reports fewer than 5,400 miles of U.S. pipelines carrying CO 2 . This deficit—and what it means for the prospect of moving substantially larger quantities of CO 2 from source to use or storage—threatens to stifle the development of CCUS projects and technologies identified as an important tool to meet emissions targets. The current regulatory landscape facing CO 2 pipeline development can best be described as uncertain. At the federal level, the U.S. DOT Pipeline and Hazardous Materials Safety Administration (PHMSA) oversees safety regulation of pipelines transporting hazardous materials, including CO 2 upon commencement of operation. However, PHMSA’s definition of CO 2 as “a fluid consisting of more than 90 percent CO 2 molecules compressed to a supercritical state” has not been updated since its 1991 addition to the Federal Register. Because CO 2 can be transported in a gaseous, liquid, or supercritical state (indeed, the physical state of CO 2 can fluctuate within a single pipeline due to environmental changes), doubts persist about the extent of PHMSA’s purview—and raise questions about what, if anything, states should do to address this apparent gap. PHMSA has begun a major revision of its existing rules, but the agency does not expect a first draft before 2024. Economic oversight of CO 2 pipelines is even less clear. The Federal Energy Regulatory Commission (FERC) and Surface Transportation Board (STB)—which regulate the rates of interstate oil/natural gas and non-energy pipelines, respectively—have both declined jurisdiction over interstate CO 2 pipelines. This presumably leaves economic regulation to state and/or local governments, but few if any states have the laws or resources in place to oversee just and reasonable rates. Further, the interstate nature of CO 2 pipeline development creates questions around how different states should align their rate-making decisions. Onshore U.S. Carbon Pipeline Deployment: Siting, Safety, and Regulation Currently, regulatory responsibilities regarding CO 2 pipeline siting and permitting fall to state and local governments. The variety of laws and regulations across the country, however, creates a maze of requirements for pipeline developers to navigate. To secure necessary permits, most states require pipeline companies to be “common carriers” that provide transport service to the public at uniform rates. However, the specific definition of that term varies. Some states require clear evidence that a pipeline services the public, while others automatically deem any pipeline company transporting energy products or hazardous materials to be a “common carrier”—with little consideration for accessibility to third parties. Other states have eschewed common-carrier terminology entirely, placing private and publicly accessible pipelines on equal footing. Much like the variation in common-carrier requirements, laws governing eminent domain authority to secure rights-of-way (ROW) to commence construction on a planned pipeline route differ by state. Several states have no laws or rules governing CO 2 pipelines. In addition to creating questions about whether long-standing rules for other pipelines (e.g., natural gas or petroleum products) apply to CO 2 , this policy vacuum leaves local governments as the sole authority over sections of pipe within their boundaries. With dozens of counties along a given route, the probability of inconsistent regulation of the same pipeline is significant. Even in states with CO 2 pipeline laws in place, local regulatory attempts to address rising concerns over pipeline routing and safety have triggered lawsuits by pipeline companies seeking to delimit areas of federal, state, and local government responsibility. Meanwhile, legislators across the country have introduced bills to restrict the application of eminent domain to CO 2 pipeline projects, which could threaten a key means of securing ROW that companies cannot secure through negotiation with landowners. Taken separately, any of these regulatory issues—the narrow federal definition of CO 2 , FERC’s and STB’s decisions that CO 2 pipelines are not within their jurisdiction, and the considerable variation in state and local governments’ laws regulating CO 2 pipeline technologies—are extremely difficult to resolve. Adding the required scale of CO 2 pipeline expansion and the currently identified narrow window of time in which to reach climate target goals, the task becomes even more difficult—and raises a host of urgent questions for regulators. How should CO 2 be defined in federal regulations to ensure consistent safety standards across the country? What is the potential impact radius of a CO 2 pipeline rupture, and how should that inform local emergency response? In the absence of centralized federal oversight, what should state legislatures do to increase alignment for interstate CO 2 pipeline projects? This paper intends to serve as a primer for regulators and stakeholders who seek to better understand the regulatory challenges and opportunities facing this critical infrastructure.

42 ENGINEERING↗
DOE OSTI2023

TRACER-Coastal Urban Boundary-Layer Interactions with Convection (TRACER-CUBIC) Field Campaign Report

To better understand the complicated web of processes governing convective cloud life cycle and aerosol-convection interactions, the U.S. Department of Energy (DOE)’s Atmospheric Radiation Measurement (ARM) user facility supported deployment of a variety of advanced atmospheric measurement systems to the greater Houston, Texas, area from 1 October 2021 to 30 September 2022 as part of the Tracking Aerosol Convection Interactions Experiment (TRACER). Houston was selected as a study area because isolated convection and a variety of aerosol conditions are common in this region. This one-year ARM Mobile Facility (AMF) deployment featured a four-month intensive operational period (IOP) during summer 2022 (1 June–30 September). The ARM instrumentation was deployed at three sites along an east-west transect from La Porte, Texas to an ancillary site in a less-polluted rural region southwest of downtown Houston (Figure 1). At the La Porte Site, which is located near the Houston ship channel in an area that experiences significant pollution, the first ARM Mobile Facility (AMF1) was deployed. During the IOP, the ARM tethered balloon system (TBS) operated at the ancillary site. The second-generation C-Band Scanning ARM Precipitation Radar (CSAPR) operated near Pearland, Texas, roughly halfway between the Laporte and ancillary sites. As part of the TRACER- Coastal Urban Boundary-Layer Interactions with Convection (CUBIC) project, three boundary-layer profiling systems) were deployed along a north-south transect spanning from the University of Houston Coastal Center to the Aldine site north of downtown Houston (also blue dot in Figure 1) during the TRACER IOP. These systems included the National Oceanic and Atmospheric Administration (NOAA) National Severe Storms Laboratory CLAMPS2 (C2), which was deployed at the UHCC, the University of Wisconsin SPARC, which was deployed at the ARM CSAPR site near Pearland (orange diamond in middle of map in Figure 1), and the University of Oklahoma CLAMPS1 (C1), which was deployed at Aldine. These three systems have been successfully operated in various field campaigns, providing data sets that collectively offer new insights into atmospheric-boundary-layer (ABL) processes, sea-breeze (SB) circulations, and convection initiation (CI). For the TRACER IOP window, these systems ran continuously between 1 June and 26 September, 2022. Due to commitments to NOAA projects, the Doppler lidar at the UHCC site was not available until 24 June 2022. The CLAMPS and SPARC profiling systems are self-contained platforms that have benefited from several years of development and deployment. Instruments and data processing were maintained remotely, which made the 4-month deployment for the TRACER-CUBIC IOP period possible. The same basic instrument configuration comprises each system: a scanning Doppler wind lidar for flow characterization and passive profiler(s) for characterizing planetary-boundary-layer (PBL) thermodynamic properties. Each platform includes a Halo Streamline Doppler wind lidar, an Atmospheric Emitted Radiance Interferometer (AERI), and a surface meteorology station. CLAMPS1 and CLAMPS2 each also include a microwave radiometer (MWR, Figure 1d). The TRACER-CUBIC hypotheses included (i) Interactions of SB and urban circulations and how they affect the PBL structure in the Houston environment, causing spatially (horizontally and vertically) and temporally highly variable flow patterns, (ii) heat, moisture, and aerosol transport and mixing depend on these flow dynamics, and (iii) an improved understanding of the flow patterns and PBL structure are critical for investigating the processes leading to CI. To test these hypotheses, the project aimed at (i) characterizing SB circulations and their impacts on the diurnal evolution of the structure of the ABL, (ii) studying the evolution of Houston’s complex urban boundary layer, and (iii) identifying effects of urban-induced circulations on pre-convective environments. TRACER-CUBIC observations generally provide good coverage during the summer IOP. Initial screening of the data indicates a good number of cases with bay-breeze (BB) and/or SB signatures, local CI, and interesting boundary-layer features such as strong nocturnal low-level jets (LLJs, Table 1). The numbers listed in rows 3-5 in this table will be further updated as part of ongoing in-depth analyses and systematic identification of local circulations and CI events. More detailed information about the data availability and quality for each instrument is provided in the “readme” files that were submitted to the ARM Data Center along with each archived data sets. These “readme” files also provide instrument descriptions, information about the data collection and processing procedures, data formats, and any additional information relevant for further data analysis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The Odd Dark Matter Halos of Isolated Gas-rich Ultradiffuse Galaxies

Abstract We analyze circular velocity profiles of seven ultradiffuse galaxies (UDGs) that are isolated and gas-rich. Assuming that the dark matter halos of these UDGs have a Navarro–Frenk–White (NFW) density profile or a Read density profile (which allows for constant-density cores), the inferred halo concentrations are systematically lower than the cosmological median, even as low as −0.6 dex (about 5 σ away) in some cases. Alternatively, similar fits can be obtained with a density profile that scales roughly as 1/ r 2 for radii larger than a few kiloparsecs. Both solutions require the radius where the halo circular velocity peaks ( R max ) to be much larger than the median expectation. Surprisingly, we find an overabundance of such large- R max halos in the IllustrisTNG dark-matter-only simulations compared to the Gaussian expectation. These halos form late and have higher spins compared to median halos of similar masses. The inner densities of the most extreme among these late-forming halos are higher than their NFW counterparts, leading to a ∼1/ r 2 density profile. However, the two well-resolved UDGs in our sample strongly prefer lower dark matter densities in the center than the simulated ones. Comparing to IllustrisTNG hydrodynamical simulations, we also find a tension in getting both low enough circular velocities and high enough halo mass to accommodate the measurements. Our results indicate that the gas-rich UDGs present a significant challenge for galaxy formation models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Effect of electrode rinse solutions on the electrodialysis of concentrated salts

The data presented indicate that modest attention to the chemistry of the electrode rinse solution can yield as much as 56% improvement in process efficiency. When operating electrodialysis (ED) with salt concentrations greater than 10,000 mg/l (as NaCl) the amperage can be relatively high even at low voltages. To maximize the overall rate of ion transport from diluate to concentrate, there is a need to minimize resistances in the electrode cells since these can represent as much as 30% of the total resistances in the entire process and in this particular ED unit, represented resistance roughly equivalent to that within the membrane stack. This work describes the performance of a 10-cell pair, 200 cm 2 (0.02 m 2 ) per membrane, pilot ED unit operated in batch mode at 5 V potential. The feedstock to the stack was varied from 0.5% to 6% NaCl. Results from Volt-amp profiles (2.0–15 V) were used as the rationale for choosing the preferred electrode rinse solution, 90 g/l (kg/m 3 ) disodium sulfate at pH 12.5. The recommended solution of 30 g/l (kg/m 3 ) disodium sulfate at neutral pH was tested against stronger solutions (60, 90, and 120 g/l (kg/m 3 ) disodium sulfate) at neutral pH and with 1 g/l (kg/m 3 ) sodium hydroxide added to yield solutions at pH 12.5. There were stark differences in the performance of the ED unit between the solutions at pH 7 and those at pH 12.5. The slopes of the Volt-amp profiles improved indicating less resistance to ion flow at pH 12.5. Most of the difference can be attributed to the greater conductivity of sodium hydroxide compared to disodium sulfate, where even a modest addition of 1 g/l (kg/m 3 ) sodium hydroxide increased the conductivity of the solution between 8 and 20% depending on the concentration of disodium sulfate. Most notably, the initiation voltage was 2.39 ± 0.08 V with solutions at pH 7 and 1.99 V ± 0.04 V when the rinse solutions were around pH 12.5, a difference of around 0.42 V. Since the full batch runs were performed at 5.0 V, this 0.42 V differential represented a loss of 8% in process efficiency. A series of tests showed a shift in initiation voltage occurred between pH 11.5 and 12.5. Nernst equations for water electrolysis are coupled with a simple flux model that estimates the concentration of hydroxide and hydronium ions at the electrode surfaces. This model indicates that the potential demanded by the electrodes is directly linked to the current. Furthermore, there is a steep drop in the potential demanded by the anode when there is sufficient hydroxide available to generate oxygen from hydroxide as opposed to generation from water. The voltage demand is minimal at pH 12.5 and reflective of the reference standard voltage (just over 1.23 V predicted at 0.1 amps). However, the flux model predicts that the voltage demanded by the electrodes at pH 7 is 1.61 V at 0.1 amps (the lowest amperage measurable with the power pack). In conclusion, comparing this result to the reference standard (E 0 = -1.23 V) yields a satisfactory explanation for the observed difference below reference standard of 0.42 V.

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

Depth-resolved soil data from DOE-LM Riverton Wyoming site, April-September 2017

Averaged (from 3 to 5 replicates) and QA/QC:d data for soil cores taken roughly every month in the period April 5 to September 13 in 2017 at a location close to DOE Legacy Management well 855 at the Riverton, Wyoming floodplain site. The groundwater at this site exhibits persistent U, Mo, and sulfate plumes and is one of the field sites in focus for the SLAC Groundwater Quality SFA program. Cores were separated into 5-20 cm segments based on soil horizonation and each segment was subsampled for microbial analyses (data available at the NCBI Single Read Archive (SRA) Database, https://www.ncbi.nlm.nih.gov/sra/, BioProject ID# PRJNA626616) in addition to the soil analyses included in this dataset. The soil analyses dataset includes aerobic and anaerobic respiration rate estimates (MicroResp method), water content estimates (on collected samples, not in situ), total C, N, S, Ca, K, Mg, Na, Fe, As, Mo, U, and organic-C concentrations, water extractable concentrations of the same elements, and HCl-extractable Fe(II) concentrations.There are four files included in this data package: metadataFile.csv, avgDataFile.csv, methodFile.csv, and RVT2017soilData.xlsx. The excel file contains the same information as the individual csv-files (in separate tabs).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Edge‐Driven Fringe‐Field Effects, Reduced Screening, and Bandgap Widening in Graphene Nanoribbons Enable Single‑Molecule Sensitivity

Graphene nanoribbons (GNRs) offer promising platforms for single‐molecule sensing due to their quasi‐1D channels and discrete electronic states, providing superior sensitivity toward molecular perturbations. While prior studies emphasize smoother edges as essential for optimal performance, the potential benefits of controlled edge roughness remain largely unexplored. Additionally, most investigations focus on GNR arrays, leaving critical edge‐ and width‐dependent factors, including fringe fields, bandgap widening, interactions between adsorbing molecules and GNR atoms, density of states (DOS) suppression, and electrostatic screening lengths, and their collective impact on sensitivity, poorly understood. Here, in this work, we fabricated field‐effect transistors using individual GNRs (widths: 200–20 nm) and characterized their response to molecular adsorption with perfluorooctanoic acid as the model analyte. Narrower ribbons displayed significantly enhanced sensitivity, yielding a coverage‐normalized response of 116 ± 10 mV per molecule in 20 nm‐wide GNRs (from calibrated ensemble Dirac‐point shifts). Experimental and theoretical analyses reveal that this heightened sensitivity arises from stronger fringe fields, width‐dependent quantum confinement effects, reduced DOS, and increased edge roughness that facilitates molecular anchoring, enhanced orbital overlap, and improved charge transfer efficiency. Our findings challenge the conventional assumption that smoother edges inherently enhance sensor performance, demonstrating that controlled edge disorder substantially boosts molecular sensitivity in GNR sensors.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Metal-ion-controlled growth and nanoindentation response of 3D, bicontinuous Cu–Fe thin films

Bicontinuous, nanocomposite thin film morphologies depend largely on the deposition conditions applied during physical vapor deposition. With the introduction of high power impulse magnetron sputtering (HiPIMS), the range of potential morphologies achieved during deposition has been increased. In this work, we compare the deposition outcomes between traditional direct-current magnetron sputtering (DCMS) and HiPIMS for a thin film co deposit of Cu and Fe. Modular control of the columnarity, porosity and roughness was achieved by varying the Cu and Fe metal ion currents during deposition. The directionality of the nanostructured phase-separated morphology was also controlled as the ion current increased. At zero ion current for both Cu and Fe sputtered species during DCMS, the film exhibited lateral concentration modulations of Cu and Fe. The directionality of the Cu- and Fe-rich phases shifted to vertical concentration modulations at low ion currents of I Fe = 1A and I Cu = 0:1A and to lateral concentration modulations at relatively moderate ion currents of I Fe = 5A and I Cu = 2A: At high ion currents of I Fe = 18A and I Cu = 2A; a more randomized phase domain structure was observed on the nanoscale. This structural shift is rationalized using an interdiffusion model. Here, the role of different kinds of phase separated morphologies, achieved during DCMS deposition, on the mechanical properties has also been studied. Results indicated increase in hardness, indentation modulus and flow strength values with the increase in indentation strain rates. Bicontinuous Cu-Fe nanocomposites are found to be stronger than multilayer Cu-Fe samples.

36 MATERIALS SCIENCE↗
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↗