SEARCH · Engineering Papers

Results for “response surface methods”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Synthesis and Characterization of Polymeric Pore Former Microspheres using Modified Polymer Resin

Porosity in polymer parts is often desired to improve mechanical properties. Individual pores must be small enough that bulk mechanical properties are retained but significant enough to reduce the density of the material. Filler materials can be incorporated into polymer matrices and subsequently removed to create voids. The materials used to create the voids are known as pore-formers. The present research aims to create spherical pores in a silicone matrix by inserting and removing core materials or fillers post-cure. A method was developed to incorporate pore formers into an uncured polymer resin and remove the materials from the matrix. Polyethylene Glycol (PEG) and Polylactic Acid (PLA) microspheres are created and incorporated into silicone matrices as pore formers. Microspheres with controllable size distributions are formed via ultrasonic spray pyrolysis and emulsion formation. Microspheres are removed from the matrix by solvent extraction and calcination. PLA microspheres were modified by hydroxyl addition to the surface to study their solubility in water. PEG-PLA block copolymers are also synthesized and investigated, comparing their physical and chemical properties to the monomers to identify whether the block copolymer microspheres will be more easily removable than their monomeric microspheres. Moreover, other properties including thermal stability, rheology, and mechanical response are studied. The present dissertation provides a comprehensive understanding of the synthesis and characterization of PEG and PLA-based pore formers, toward the future development of feedstock for additive manufacturing porous structures.

Damptey, Ransford↗

DOE OSTI2020

Time-Domain Nuclear Magnetic Resonance Determination of Wettability Alteration: Analysis for Low-Salinity Water

Wettability has been shown to influence oil recovery. This property has become central to low-salinity (LSW) and smart (SWF) water flooding recovery mechanisms research. The challenge lies in the fact that oil recovery results from the combined effects of solid-liquid and liquid-liquid interactions. This demands methods that allow an independent interpretation of wettability alteration contributions. The primary objective of this work is to assess changes in wettability through the application of Time-Domain Nuclear Magnetic Resonance (TD-NMR) T 2 distribution and diffusion coefficient, starting with a well-controlled porous system, that is, glass beads, and then a model rock (Berea), in the presence of one phase, either oil or brine exclusively. Subsequently, two-phase fluid saturation was tested. For the glass beads, dimethyldichlorosilane was used to induce a hydrophobic response, as confirmed by contact angle experiments on slides of the same material. Sodium sulfate was used for its known positive influence on oil recovery during LSW and SWF. In cases where alteration of surface properties was expected, a leftward shift of the average T 2 distribution curve modes, accompanied by a reduction on the diffusion coefficient during the aging process was observed. The results of this work confirm that fluid-solid interactions during LSW and SWF, namely a shift in wettability, take place after the injection of low-salinity water.

42 ENGINEERING↗

DOE OSTI2021

Mechanical Loads on Spent Nuclear Fuel in the General 30 cm Package Drop Scenario

The U.S. Department of Energy Office of Nuclear Energy (DOE-NE) Spent Fuel and Waste Science & Technology (SFWST) research program is guided by the high-level goal of closing prioritized knowledge gaps related to spent nuclear fuel (SNF) storage and transportation, which are summarized by Saltzstein et al. (2020). One of the high-priority knowledge gaps is the identification and quantification of mechanical loads that are expected to affect SNF during normal conditions of transportation and storage to inform the range of physical SNF test programs. This report uses modeling and analysis methods to estimate the mechanical loads on spent nuclear fuel (SNF) in the general 30 cm package drop scenario. The drop scenario assumes impact limiters are in place in the transportation configuration and the impact surface is perfectly rigid. The goal of this analysis is to consider the universe of potential mechanical loading conditions that can happen to SNF and present the results in a manner that is useful for materials testing, decision making, and regulatory rule making purposes. This study uses validated finite element models and methods to perform a broad parametric study of key variables that can affect the mechanical loads on SNF during a hypothetical package free drop scenario. Physical drop test data from a cask and fuel assembly drop test campaign is the basis for model validation. Additionally, the results of the parametric study are used to inform a damage model, which uses multiple nonlinear regression to estimate the relationships between input variables and output response. The parametric finite element analyses consider thousands of input variable combinations, while the damage model estimates millions of combinations. The breadth of this study provides confidence that the potential range of mechanical loads that SNF might experience during the general 30 cm package drop scenario are characterized well enough to consider this knowledge gap closed. While this report documents the overall peak values calculated in this study, the 95 th percentile values, the histograms, and the observed trends are equally important. This study covered a large range of SNF temperatures, room temperature to 300°C, and burnups, 10 GWd/MTU to 62 GWd/MTU. Each temperature and burnup combination has a different cladding yield strain, so it is more meaningful to summarize the calculated cladding strain response as its factor of safety, which is defined relative to the yield strain. The factor of safety is calculated as the yield strain divided by the peak cladding strain. A factor of safety greater than unity indicates that the cladding remains below yield, whereas a value less than unity is indicative of plastic deformation. In all cases of this study a safety factor over 1.0 was calculated, although in the most limiting case at 300°C the safety factor was only 1.01, which suggests that yielding could occur when additional loads like rod internal pressure are included. When the temperature is restricted to 200°C the limiting safety factor increases to 1.28, which has significant margin to accommodate internal pressure and potential local cladding defects that could cause a local stress concentration. An important trend in the calculated fuel rod mechanical loads is that the 2 nd highest loaded fuel rod in an assembly tends to be significantly lower than the highest loaded rod. The implication is that even if one rod in an assembly experiences a failure the loads would have to be significantly higher to cause two or more rods to fail.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Enhancing Photosynthesis Simulation Performance in ESMs with Machine Learning-Assisted Solvers

When simulating vegetation dynamics, photosynthesis accounts for a large fraction of the computational cost in most Earth System Models (ESMs). This is largely since photosynthesis is represented as a system of nonlinear equations, and the solution requires the use of an initial guess followed by many iterations of the numerical solver to obtain a solution. We use machine learning (ML) to replicate the response surface of the model’s numerical solver to improve the choice of initial guess, therefore requiring fewer iterations to obtain a final solution. We implemented this test on the leaf-level calculations as well as at the canopy scale, and for both we observed fewer iterations of the photosynthesis solver when a ML-based initial guess was implemented. The model tested here is the Energy Exascale Earth System Model - Land Model (ELM). The ML-based algorithms used here are trained on simulations from the model itself and used only to improve the initial guess for the solver; therefore, the model maintains its own set of physics to obtain the final solution. This work shows novel ways to utilize ML-based methods to improve the performance of numerical solvers in ESMs.

Massoud, Elias [ORNL] (ORCID:0000000217725361)↗

DOE OSTI2020

A machine learning approach to emulation and biophysical parameter estimation with the Community Land Model, version 5

Abstract. Land models are essential tools for understanding and predicting terrestrial processes and climate–carbon feedbacks in the Earth system, but uncertainties in their future projections are poorly understood. Improvements in physical process realism and the representation of human influence arguably make models more comparable to reality but also increase the degrees of freedom in model configuration, leading to increased parametric uncertainty in projections. In this work we design and implement a machine learning approach to globally calibrate a subset of the parameters of the Community Land Model, version 5 (CLM5) to observations of carbon and water fluxes. We focus on parameters controlling biophysical features such as surface energy balance, hydrology, and carbon uptake. We first use parameter sensitivity simulations and a combination of objective metrics including ranked global mean sensitivity to multiple output variables and non-overlapping spatial pattern responses between parameters to narrow the parameter space and determine a subset of important CLM5 biophysical parameters for further analysis. Using a perturbed parameter ensemble, we then train a series of artificial feed-forward neural networks to emulate CLM5 output given parameter values as input. We use annual mean globally aggregated spatial variability in carbon and water fluxes as our emulation and calibration targets. Validation and out-of-sample tests are used to assess the predictive skill of the networks, and we utilize permutation feature importance and partial dependence methods to better interpret the results. The trained networks are then used to estimate global optimal parameter values with greater computational efficiency than achieved by hand tuning efforts and increased spatial scale relative to previous studies optimizing at a single site. By developing this methodology, our framework can help quantify the contribution of parameter uncertainty to overall uncertainty in land model projections.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Excitation of long electrostatic solitary waves in ion beam neutralization process

Unusually long electrostatic solitary waves (ESWs) are discovered in a particle-in-cell simulation study of the process of ion beam neutralization by electron emission from a filament. These ESWs are long because the density perturbation responses to the potential wells created by the ESWs are very small. Notably, the density perturbation is small because the trapped (positive) and untrapped (negative) electron density perturbations nearly compensate each other because of a non-Maxwellian electron velocity distribution in the beam.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Improved Methods for Modeling Functional Transition Metal Compounds in Complex Environments: Ground States, Excited States, and Spectroscopies

Rapid advances in energy applications require new theory and computational models to provide guidance for interpretation of experimental results and mechanistic understanding. New theory development is necessary to treat systems of increasing complexity, size, and relevance to real applications. Functional transition metal compounds, including molecules, clusters, nanoparticles, surfaces, and solids, provide particular promise for magnetic, optical, and catalytic applications. However, such systems can be exceptionally challenging to model. To make progress in understanding and designing transition metal compounds for energy applications, theory must be able to simulate such systems in complex environments, as well as simulate the spectra of such complex systems to provide direct connections with experiment. Leveraging the independent expertise of the team’s members, this project aimed to make inroads to the theoretical and computational challenges associated with studying transition metal compounds, their reaction chemistry, photophysics and photochemistry, and response to spectroscopic interrogation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Rapid Determination of Uranium Isotopic Abundance from Cotton Swipes: Direct Extraction via a Planar Surface Reader and Coupling to a Microplasma Ionization Source

The collection of solid particulates and liquids from surfaces by the use of cloth swipes is fairly ubiquitous. In such methods, there is a continuous concern regarding the ability to locate and quantitatively sample the analyte species from the material. In this effort, we demonstrate the initial coupling of an Advion Plate Express plate reader to a liquid sampling–atmospheric pressure glow discharge (LS-APGD) microplasma ionization source with an Orbitrap mass spectrometer to perform uranium isotopic analyses of solution residues on cotton swipes. The Plate Express employs a sampling probe head to engage and seal against the swipe surface. Subsequentially, the analyte residues are desorbed and transported within a 2% HNO 3 electrolyte flow to the ionization source. Quantitative recoveries were observed following a single 30 s extraction step, with the absolute mass sampled per extraction being ~100 ng. While the intrasample variability in the analytical responses for triplicate sampling of the same swipe yield ~30% RSD, this lack of precision is offset by the ability to determine isotope ratios for enriched uranium specimens with a precision of better than 10% RSD. Pooled, intersample precision ( n = 9) was found to be <5%RSD across the various sample compositions. Finally, 235 U/ 238 U determinations (ranging from 0.053 to 1.806) were accurate with errors of <10%, absolute. The 234 U- and 236 U-inclusive ratios were determined with similar accuracy in enriched samples. While the driving force for the effort is in the realm of nuclear nonproliferation efforts, the ubiquitous use of cloth swipes across many application areas could benefit from this convenient approach, including the use of versatile, reduced-format mass spectrometer systems.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2025

The importance of accounting for the Tolman correction to surface tension for nucleation and growth modeling of Fe clusters

Gibbs free energies of clusters are required for predictive modeling of cluster growth during condensation of a cooling vapor. Here, we present a straightforward method of calculating free energies of cluster formation using the data from molecular dynamics (MD) simulations. We apply this method to iron clusters having from 2 to 100 atoms. The energies obtained are verified by comparing to an MD-simulated equilibrium cluster size distribution in a sub-saturated vapor. We show that these free energies differ significantly from those obtained with a commonly used spherical cluster approximation, which relies on a surface tension coefficient of a flat surface, as it is used in the classical nucleation theory (CNT). We show that the spherical cluster approximation in CNT can be improved by using a cluster-size-dependent Tolman correction for the surface tension. The Tolman length and effective surface tension values were derived for iron clusters, and they significantly differ from the commonly used experimentally measured values. This improved approximation does not account for geometric magic number effects responsible for spikes and troughs in densities of neighbor cluster sizes. Nonetheless, it allows to more accurately model cluster formation from a cooling vapor. It better reproduces the condensation timeline, overall shape of the cluster size distribution, average cluster size, and the distribution width. In contrast, using a constant surface tension coefficient (as done in CNT) resulted in incorrect condensation dynamics and cluster size distributions. The analytical expression for cluster nucleation rate from CNT was updated to account for the size-dependence of cluster surface tension.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Specific conductivity and salinity of the Parker River, PIE LTER, Plum Island Sound MA, August-November 2022

This dataset contains specific conductivity and calculated salinity data of Parker River water at a tidal brackish wetland dominated by Typha angustifolia at the upper estuary of the Plum Island Sound in Newbury, Massachusetts (MA) within the Plum Island Ecosystems Long Term Ecological Research site (PIE LTER). Measurements were taken to evaluate temporal changes in surface water salinity in high frequency to characterize boundary conditions of soil and plant responses to changes in salinity. A PVC pipe was installed in a low elevation spot in the creek bank so that the bottom of the pipe sat on the sediment surface allowing flushing with water during flooding. Raw measurements were collected using an Onset HOBO U24-002 Saltwater Conductivity/Salinity data logger. The specific conductance and salinity measurements were corrected and calculated respectively using Onset’s HOBOware software and reference specific conductivity measurements taken in tandem with the first and last points recorded by the HOBO sensor. These reference measurements were taken using a HACH HQ14D Portable Conductivity Meter. Because of the installation design, only data one hour before and after high tide are used. Metadata files Typha_ctd_salinity_dd.csv and Typha_ctd_salinity_flmd.csv contain detailed information on data variables, sampling and QA/QC methods, and site location.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Specific conductivity and salinity of the Parker River, PIE LTER, Plum Island Sound MA, March-November 2023

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Exploring Whether Subsurface Fluid Production Can Minimize Triggered Seismicity in Geothermal Fields

Fluid injection and production related to energy recovery and other industrial operations alter the pressure and stress state of subsurface reservoirs which can lead to induced seismicity. Here, the primary goal is to investigate the hypothesis that the modulation of stress state in subsurface reservoirs through fluid injection/production operations can reduce the likelihood of inducing seismicity. Validation of this mitigation strategy will provide an active operational method to control induced seismicity in subsurface reservoir exploitations such as geothermal energy recovery and carbon storage. In this project, we analyze the relationship between extensive fluid production and the paucity of aftershock activity at the Coso Geothermal plant (CGP) following the 2019 Ridgecrest earthquake sequence, where high rates of aftershock triggering were expected. We developed a high-fidelity multiphase coupled thermo-hydro-mechanical (THM) model to simulate fluid injection/production activities at the CGP between 1986 and 2020. THM results of surface subsidence due to high rates of fluid production agree well with field observations from global positioning system (GPS) and interferometric synthetic aperture radar (InSAR). Subsequently, the variations of pore pressure, temperature and stress state were used to drive numerical earthquake simulations of the aftershock response to the 2019 Ridgecrest earthquakes. The earthquake simulations show that seismic quiescence may occur following the Ridgecrest event, depending on the initial stress state at the time of the mainshock and at the initiation of geothermal energy production. Seismic quiescence occurs in 20% of the cases we explored, where rates are decreased by 50% of the background rate in the two years following the mainshock. In circumstances where aftershock rates increase near the CGP following Ridgecrest, the average rate change is a factor of two larger than background rates from simulations with no operations. These findings indicate that unlike many other operations that bring faults closer to failure, operations at the CGP are acting to stabilize faults such that triggering in the current stress state is minimal.

58 GEOSCIENCES↗

DOE OSTI2023

Extent of interlocking and metallurgical bonding in friction riveting of aluminum alloy to steel

In this study, the joining of 6061-T6 aluminum alloy and DP590 steel using a M42 steel rivet via friction riveting technique is investigated. The surface morphology and microstructure characterization reveal the formation of an anchor zone that imparts mechanical interlock as well as the formation of metallurgical bonds at the interface of aluminum and steel. A combination of interlocking and bonding results in the achievement of a high load-carrying capacity of 5.7 kN during lap shear testing at room temperature. A finite element-based computational model was developed which accurately predicted the lap shear response of the joint. The model revealed that the metallurgical bond formed during fric-riveting adds 39% peak load strength to the joint. An extensive microstructural investigation, post-lap-shear fractography, and the modeling results, together provided insights on the joint failure mechanism. In conclusion, this study highlights that friction riveting is a promising method for aluminum-to-steel dissimilar joining, which is important for lighweighing automotive vehicles for energy efficiency.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Characterization of Temperature Heterogeneity in Utility-Scale Power Plant Boilers by Spatially Distributed Ultrasonic Measurements

In extreme environments, even hardened insertion sensors fail quickly. For such environments, we have developed an ultrasonic (US) method for measuring the spatial distribution of temperatures in solid materials and, specifically, across containments of extreme processes. By deploying US sensors in multiple locations, spatial heterogeneity in temperatures inside harsh environments may be characterized. In this note, we update our progress on testing the developed approach in an industrial setting and its performance in describing the temperature distribution inside a 500 MW (electrical) coalfired utility boiler of an electrical power generation plant. We use waveguides (WG) structured to contain echogenic features that produce a train of echoes in response to an excitation pulse. The time of flight between echoes encodes the information on the temperature distribution in the corresponding segment of the WG, which we reconstruct using its parametrization. Five waveguides were welded to the boiler’s heat exchange surface (waterwall) and different locations at the same boiler elevation and produced reliable US waveforms 18 months after the installation. For several weeks during the latest trial, we performed US measurements in five locations spanning half of the boiler’s 14-meter width. The measurements were acquired during regular operation, including load cycling to adjust for demand and intermittent contributions of renewable power sources, and during boiler shutdown for emergency maintenance and the subsequent restart. The interpretation of US waveforms by signal processing resulted in an accurate estimation of temperature distribution along the waveguides. We captured daily cyclical variations in demand-following load. Measurements performed in multiple locations revealed unexpected temperature variations across the boiler. We conclude that an array of developed sensors can provide responsive and spatial temperature measurements while maintaining functionality despite prolonged exposure to an extreme environment.

Walton, Kenneth↗

DOE OSTI2020

Coupling between plasmonic and photonic crystal modes in suspended three-dimensional meta-films

A complementary metal oxide semiconductor (CMOS) compatible fabrication method for creating three-dimensional (3D) meta-films is presented. In contrast to metasurfaces, meta-films possess structural variation throughout the thickness of the film and can possess a sub-wavelength scale structure in all three dimensions. Here we use this approach to create 2D arrays of cubic silicon nitride unit cells with plasmonic inclusions of elliptical metallic disks in horizontal and vertical orientations with lateral array-dimensions on the order of millimeters. Fourier transform infrared (FTIR) spectroscopy is used to measure the infrared transmission of meta-films with either horizontally or vertically oriented ellipses with varying eccentricity. Shape effects due to the ellipse eccentricity, as well as localized surface plasmon resonance (LSPR) effects due to the effective plasmonic wavelength are observed in the scattering response. The structures were modeled using rigorous coupled wave analysis (RCWA), finite difference time domain (Lumerical), and frequency domain finite element (COMSOL). The silicon nitride support structure possesses a complex in-plane photonic crystal slab band structure due to the periodicity of the unit cells. We show that adjustments to the physical dimensions of the ellipses can be used to control the coupling to this band structure. The horizontally oriented ellipses show narrow, distinct plasmonic resonances while the vertically oriented ellipses possess broader resonances, with lower overall transmission amplitude for a given ellipse geometry. We attribute this difference in resonance behavior to retardation effects. The ability to couple photonic slab modes with plasmonic inclusions enables a richer space of optical functionality for design of metamaterial-inspired optical components.

42 ENGINEERING↗

DOE OSTI2020

Accelerated Snowmelt Protocol to Simulate Climate Change Induced Impacts on Snowpack Dependent Ecosystems

Field studies that simulate the effects of climate change are important for a predictive understanding of ecosystem responses to a changing environment. Among many concerns, regional warming can result in advanced timing of spring snowmelt in snowpack dependent ecosystems, which could lead to longer snow-free periods and drier summer soils. Past studies investigating these impacts of climate change have manipulated snowmelt with a variety of techniques that include manual snowpack alteration with a shovel, infrared radiation, black sand and fabric covers. Within these studies however, sufficient documentation of methods is limited, which can make experimental reproduction difficult. Here, we outline a detailed plot-scale protocol that utilizes a permeable black geotextile fabric deployed on top of an isothermal spring snowpack to induce advanced snowmelt. The method offers a reliable and cost-effective approach to induce snowmelt by passively increasing solar radiation absorption at the snow surface. In addition, control configurations with no snowpack manipulation are paired adjacent to the induced snowmelt plot for experimental comparison. Past and ongoing deployments in Colorado subalpine ecosystems indicate that this approach can accelerate snowmelt by 14-23 days, effectively mimicking snowmelt timing at lower elevations. This protocol can be applied to a variety of studies to understand the hydrological, ecological, and geochemical impacts of regional warming in snowpack dependent ecosystems.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Mineralogical, magnetic and geochemical data constrain the pathways and extent of weathering of mineralized sedimentary rocks

The oxidative weathering of sulfidic rock can profoundly impact watersheds through the resulting export of acidity and metals. Weathering leaves a record of mineral transformation, particularly involving minor redox-sensitive phases, that can inform the development of conceptual and quantitative models. In sulfidic sedimentary rocks, however, variations in depositional history, diagenesis and mineralization can change or overprint the distributions of these trace minerals, complicating the interpretation of weathering signatures. Here we show that a combination of bulk mineralogical and geochemical techniques, micrometer-resolution X-ray fluorescence microprobe analysis and rock magnetic measurements, applied to drill core samples and single weathered fractures, can provide data that enable the development of a geochemically consistent weathering model. This work focused on one watershed in the Upper Colorado River Basin sitting within the Mesaverde Formation, a sedimentary sandstone bedrock with disseminated sulfide minerals, including pyrite and sphalerite, that were introduced during diagenesis and subsequent magmatic-hydrothermal mineralization. Combined analytical methods revealed the pathways of iron (Fe), carbonate and silicate mineral weathering and showed how pH controls element retention or release from the actively weathering fractured sandstone. Drill core logging, whole rock X-ray diffraction, and geochemical measurements document the progression from unweathered rock at depth to weathered rock at the surface. X-ray microprobe analyses of a 1-cm size weathering profile along a fracture surface are consistent with the mobilization of Fe(II) and Fe(III) into acidic pore water from the dissolution of primary pyrite, Fe-sphalerite, chlorite, and minor siderite and pyrrhotite. These reactions are followed by the precipitation of secondary minerals such as of goethite and jarosite, a Fe-(oxyhydr)oxide and hydrous Fe(III) sulfate, respectively. Microscale analyses also helped explain the weathering reactions responsible for the mineralogical transformations observed in the top and most weathered section of the drill core. For example, dissolution of feldspar and chlorite neutralizes the acidity generated by Fe and sulfide mineral oxidation, oversaturating the solution in both Fe-oxides. The combination of X-ray spectromicroscopy and magnetic measurements show that the Fe(III) product is goethite, mainly present either as a coating on fracture surfaces in the actively weathering region of the core or more homogeneously contained within the unconsolidated regolith at the top of the core. Low-temperature magnetic data reveal the presence of ferromagnetic Fe-sulfide pyrrhotite that, although it occurs at trace concentrations, could provide a qualitative proxy for unweathered sulfide minerals because the loss of pyrrhotite is associated with the onset of oxidative weathering. Pyrrhotite loss and goethite formation are detectable through room-temperature magnetic coercivity changes, suggesting that rock magnetic measurements can determine weathering intensity in rock samples at many scales. In conclusion, this work contributes evidence that the weathering of sulfidic sedimentary rocks follows a geochemical pattern in which the abundance of sulfide minerals controls the generation of acidity and dissolved elements, and the pH-dependent mobility of these elements controls their export to the ground- and surface-water.

58 GEOSCIENCES↗

DOE OSTI2021

Calibration and field deployment of low-cost sensor network to monitor underground pipeline leakage

Recent technological advances in methane detection have improved leak detection and repair. However, current methods to reliably measure methane concentrations rely on expensive instruments or demand significant labor input. There is interest in using affordable methane sensors that are responsive to ppmv level changes in methane concentrations in both urban and rural environments for monitoring underground natural gas pipeline leaks. This is especially relevant for situations where potentially significant leaks cannot be repaired immediately or smaller leaks that require long-term monitoring and further evaluation. In this work, a low-cost sensor unit, equipped with a metal oxide sensor, was designed, built, and calibrated over a wide range of methane concentrations and environmental conditions in preparation for field application. A network of these sensors was then installed at the test site and used to measure methane concentrations at ground level above known sub-surface natural gas emissions which emulated underground gas pipeline leaks. This low-cost sensor network measured over 4 days total for the two different known leakage rates. Results demonstrate that the sensors can continuously measure relative methane variability for extended periods but require calibration for a wide range of temperature and humidity conditions to properly determine absolute gas (i.e., methane) concentrations. Furthermore, when a regression analysis was conducted to evaluate the effects of meteorological parameters on methane concentration, air temperature and wind speed have strong impacts on the concentration. Overall, the network approach allows improved identification of leak location and monitoring of underground natural gas leaks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗