Search NASA⌕ Search

SEARCH · Search NASA

Results for “sample average approximation”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Vitrification of Hanford Tank Wastes for Condensate Recycle and Feed Composition Changeover Testing

During the vitrification of Hanford Site nuclear waste at the Waste Treatment and Immobilization Plant (WTP), which is a part of the safe and efficient retrieval, treatment, and disposal mission of the U.S. Department of Energy Office of River Protection, the offgas condensate generated from the waste-to-glass conversion is currently planned to be concentrated by evaporation in the Effluent Management Facility (EMF). This concentrated condensate can then be recycled back to the incoming waste and vitrified. To test the recycle process, an apparatus was designed and built to mimic the EMF evaporator and was then used to concentrate a volume of condensate produced during the vitrification of a sample of Hanford tank 241-AP-107 (referred to herein as AP-107) waste in a continuous laboratory-scale melter (CLSM). The concentrated condensate was added to an additional sample of AP-107 waste, to mimic one round of the recycle process, and the combined solution was vitrified, producing a second round of recycle condensate. In the current study, the EMF test apparatus was used to concentrate the second-round recycle condensate under evaporation conditions (at 45 °C and 1.4 psia) designed to emulate EMF operation. The condensate was successfully concentrated by a factor of ~10 while retaining over 95 % of the technetium-99 ( 99 Tc), Cs, and I inventories in the concentrate. Another portion of AP-107 waste was retrieved by Washington River Protection Solutions, LLC (WRPS) and transferred to Pacific Northwest National Laboratory (PNNL), where it was pretreated and then combined with the second-round recycle AP-107 condensate concentrate and glass-forming chemicals (GFCs) to form the two-time recycle AP-107 melter feed, approximating a second round to the recycling action to be performed at the WTP. A portion of AP-105 waste was also retrieved by WRPS and provided to PNNL for pretreatment and combining with GFCs to form AP-105 melter feed. The two-time recycle AP-107 and AP-105 melter feeds were processed consecutively in the CLSM. The CLSM run proceeded for 13.63 hours, producing 9.70 kg of glass for an average glass production rate of 1464 kg m 2 d -1 during the two-time recycle AP-107 feed charging and 1568 kg m 2 d -1 during the AP-105 feed charging. The rate during AP-107 charging was essentially equivalent to the rate when processing no-recycle AP-107 feed and lower than that achieved when processing one-time recycle AP-107 feed. However, all rates were within the potential range of variability when processing melter feeds with similar composition in the CLSM. Likewise, the rate during AP-105 charging was higher than the previous rate processing AP-105, but within the potential CLSM range. The cold-cap characteristics changed from the typically thin AP-107 cold cap to a foamy-edged cold cap as previously seen with AP-105 shortly after transitioning to the AP-105 melter feed. The glass produced during the CLSM run was within 10 % of its target composition for the primary glass components. The CaO and Li 2 O targets varied by more than 1 wt% between the two-time recycle AP-107 and AP-105 glass targets and it took about 2 turnovers of the CLSM glass inventory to reach a relative chemical steady state in the glass for CaO and Li 2 O after the melter feed inputs were switched. The 99 Tc and total cesium content in the melter feeds were maintained at concentrations expected to be experienced at the WTP. During the CLSM run, while processing the two-time recycle AP-107 melter feed at a relative chemical steady state, the 99 Tc/Cs ratio was 10, and 34% of 99 Tc and 74% of Cs were retained in the glass. These values were higher than those measured in the CLSM run with one-time recycle AP-107 melter feed. After the transition to processing the AP-105 melter feed, when the production reached a relative chemical steady state, the 99 Tc/Cs ratio was 77 while 44% of 99 Tc was retained in the glass. The C's retention during this time frame reached 200% due to the excess C's in the glass after the target content decreased to 15% of its initial level in the two-time recycle AP-107 melter feed to the lower target in the Ap-105 melter feed. While iodine was below inductively coupled plasma mass spectrometry analytical reporting limits in the melter feed and glass samples, it was detected in quantities above the analytical reporting limits in the liquid and filter samples collected from the CLSM offgas treatment system. The behavior of iodine in the CLSM offgas treatment system followed a similar pattern to those of 99 Tc and Cs.

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

Fiscal Year 2022. Filtration of Hanford Tank 241-AP-101 Supernatant at 16 °C

Bench-scale filtration testing of ~9 liters of supernatant from Hanford waste tank 241-AP-101, chilled to 16 °C, was conducted using a backpulse dead-end filter (BDEF) filtration system equipped with a feed vessel and a Mott inline filter Model 6610 (Media Grade 5) in the hot cells of the Radiochemical Processing Laboratory at Pacific Northwest National Laboratory. This was done to assess the performance of the anticipated third feed to the Tank Side Cesium Removal (TSCR) system. The as-received samples were diluted to the target sodium concentration and transferred to 1.5-liter polyethylene bottles and held at 16 °C for approximately 1 week prior to filtration. The feed was filtered through the BDEF system at a targeted flux of 0.065 gpm/ft 2 to match the prototypic operation of the TSCR system. During filtration, the differential pressure required to effect filtration at 0.065 gpm/ft 2 increased little over the filtration campaign and never reached 2 psid (the TSCR action limit). This indicates that the TSCR filter should perform well when processing AP-101 supernatant. After completing filtration of the AP-101 feed, the filter was cleaned. Solids concentrated from the backpulse solutions displayed sodium nitrate-type phases, aluminum and silicon phases reported as cancrinite or nitrate-cancrinite, a mixed chromium-aluminum oxide, iron oxides, and Ca-bearing phases (calcite). Scanning electron microscopy analyses showed that the average particle size was 0.5 micron.

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

A hybrid 3D/2D field response calculation for liquid argon detectors with PCB based anode plane

Liquid Argon Time Projection Chamber (LArTPC) technology is commonly utilized in neutrino detector designs. It enables detailed reconstruction of neutrino events with high spatial precision and low energy threshold. Its field response (FR) model describes the time-dependent electric currents induced in the anode-plane electrodes when ionization electrons drift nearby. An accurate and precise FR is a crucial input to LArTPC detector simulations and charge reconstruction. Established LArTPC designs have been based on parallel wire planes. It allows accurate and computationally economic two-dimensional (2D) FR models utilizing the translational symmetry along the direction of the wires. Recently, novel LArTPC designs utilize electrodes formed on printed circuit board (PCB) in the shape of strips with through holes. The translational symmetry is no longer a good approximation near the electrodes and a new FR calculation that employs regions with three dimensions (3D) has been developed. Extending the 2D models to 3D would be computationally expensive. Fortuitously, the nature of strips with through holes allows for a computationally economic approach based on the finite-difference method (FDM). In this paper, we present a new software package pochoir that calculates LArTPC field response for these new strip-based anode designs. This package combines 3D calculations in the volume near the electrodes with 2D far-field solutions to achieve fast and precise field response computation. We apply the resulting FR to simulate and reconstruct samples of cosmic-ray muons and 39 Ar decays from a Vertical Drift (VD) detector prototype operated at CERN. We find the difference between real and simulated data within 5%. Current state-of-the-art LArTPC software requires a 2D FR which we provide by averaging over one dimension and estimate that variations lost in this average are smaller than 7%.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Correcting for filter-based aerosol light absorption biases at the Atmospheric Radiation Measurement program's Southern Great Plains site using photoacoustic measurements and machine learning

Abstract. Measurement of light absorption of solar radiation by aerosols is vital for assessing direct aerosol radiative forcing, which affects local and global climate. Low-cost and easy-to-operate filter-based instruments, such as the Particle Soot Absorption Photometer (PSAP), that collect aerosols on a filter and measure light attenuation through the filter are widely used to infer aerosol light absorption. However, filter-based absorption measurements are subject to artifacts that are difficult to quantify. These artifacts are associated with the presence of the filter medium and the complex interactions between the filter fibers and accumulated aerosols. Various correction algorithms have been introduced to correct for the filter-based absorption coefficient measurements toward predicting the particle-phase absorption coefficient (Babs). However, the inability of these algorithms to incorporate into their formulations the complex matrix of influencing parameters such as particle asymmetry parameter, particle size, and particle penetration depth results in prediction of particle-phase absorption coefficients with relatively low accuracy. The analytical forms of corrections also suffer from a lack of universal applicability: different corrections are required for rural and urban sites across the world. In this study, we analyzed and compared 3 months of high-time-resolution ambient aerosol absorption data collected synchronously using a three-wavelength photoacoustic absorption spectrometer (PASS) and PSAP. Both instruments were operated on the same sampling inlet at the Department of Energy's Atmospheric Radiation Measurement program's Southern Great Plains (SGP) user facility in Oklahoma. We implemented the two most commonly used analytical correction algorithms, namely, Virkkula (2010) and the average of Virkkula (2010) and Ogren (2010)–Bond et al. (1999) as well as a random forest regression (RFR) machine learning algorithm to predict Babs values from the PSAP's filter-based measurements. The predicted Babs was compared against the reference Babs measured by the PASS. The RFR algorithm performed the best by yielding the lowest root mean square error of prediction. The algorithm was trained using input datasets from the PSAP (transmission and uncorrected absorption coefficient), a co-located nephelometer (scattering coefficients), and the Aerosol Chemical Speciation Monitor (mass concentration of non-refractory aerosol particles). A revised form of the Virkkula (2010) algorithm suitable for the SGP site has been proposed; however, its performance yields approximately 2-fold errors when compared to the RFR algorithm. To generalize the accuracy and applicability of our proposed RFR algorithm, we trained and tested it on a dataset of laboratory measurements of combustion aerosols. Input variables to the algorithm included the aerosol number size distribution from the Scanning Mobility Particle Sizer, absorption coefficients from the filter-based Tricolor Absorption Photometer, and scattering coefficients from a multiwavelength nephelometer. The RFR algorithm predicted Babs values within 5 % of the reference Babs measured by the multiwavelength PASS during the laboratory experiments. Thus, we show that machine learning approaches offer a promising path to correct for biases in long-term filter-based absorption datasets and accurately quantify their variability and trends needed for robust radiative forcing determination.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

High resolution variability in wet deposition in the southeastern United States

Rainwater chemistry is determined by atmospheric pollutants and particles which vary spatially and temporally. Industrial and agricultural activities and meteorological events (e.g. sea breezes, severe weather, blowing dust) alter atmospheric particle and trace gas compositions. These gases and particles are scavenged by cloud and rain droplets that drive wet deposition. During an Intensive Operation Period (IOP) from April to October 2021, rainwater was collected at higher frequency intervals, usually daily, after precipitation events at three locations on the Savannah River Site (SRS). The farthest locations were separated by approximately 20 km. The mean concentration (μeq/L) of seven ions followed the Cl⁻ > SO 4 2− > Na⁺ > NO 3 ⁻ > K⁺ > Mg 2+ > Ca 2+ downward trend. Ion concentrations were compared to background ion concentrations from the National Atmospheric Deposition Program (NADP). The high frequency monthly averaged SRS data compared well with the monthly averaged NADP background but demonstrated extensive variability. In some months in 2021, the high frequency data compared better to the NADP site near the coast while in other months inland sites compared better. Strong spatial variability for ion concentrations was observed across SRS which was attributed to localized impacts in rainfall spatial variability. High frequency measurements allowed for comparison to regional weather patterns indicating influences from the Atlantic Ocean, Gulf of Mexico, and cities. This can account for spatial variability in the wet deposition flux. Sea breezes, Saharan dust, and anthropogenic sources were shown to impact wet deposition flux variability. Higher frequency precipitation chemistry sampling at numerous locations better captures ion concentration variability and improves measurement representativeness.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

High temperature neutron diffraction investigation of U 3 Si 2

The U-Si system is actively undergoing studies due to its promise as a component of an accident tolerant nuclear fuel. In this work, the crystal structure of the U 3 Si 2 compound in the U-Si system was investigated as a function of temperature from room temperature to 1298 K for a sample of approximately 42.02 at.% Si composition and 1373 K for a sample with approximately 42.10 at.% Si using high temperature time-of-flight neutron diffraction on the High-Pressure Preferred Orientation (HIPPO) diffractometer at Los Alamos Neutron Science Center (LANSCE). The simultaneous Rietveld refinement of five histograms from the five HIPPO detector backs (40°, 60°, 90°, 120°, and 145°) provided datasets for the lattice parameters, anisotropic atomic displacement parameters, and atomic positions as a function of temperature. To explore the possibility of a homogeneity range, two sample compositions were analyzed, a 40.02 at.% Si (near stoichiometric) and a potentially hyperstoichiometric (40.10 at.% Si) sample. While minor differences in the anisotropic atomic displacement parameters between the two samples were observed, over the entire investigated temperature range no additional phases were found indicative of a single phase sample and the ability of the U 3 Si 2 compound to accommodate nonstoichiometry, suggesting that the U 3 Si 2 compound is not a line compound confirming previous determinations. However, the differences observed in the average coefficient of linear thermal expansion identified between the two compositions warrant future investigation.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Classical combinatorial optimization scaling for random Ising models on 2D heavy-hex graphs

Motivated by near term quantum computing hardware limitations, combinatorial optimization problems that can be addressed by current quantum algorithms and noisy hardware with little or no overhead are used to probe capabilities of quantum algorithms such as the quantum approximate optimization algorithm. In this study, a specific class of near term quantum computing hardware defined combinatorial optimization problems, Ising models on heavy-hex graphs both with and without geometrically local cubic terms, are examined for their classical computational hardness via empirical computation time scaling quantification. Specifically the time-to-solution (TTS) metric using the classical heuristic simulated annealing is measured for finding optimal variable assignments (ground states), as well as the time required for the optimization software Gurobi to find an optimal variable assignment. Because of the sparsity of these Ising models, the classical algorithms are able to find optimal solutions efficiently even for large instances (i.e. 100 000 spin variables). The Ising models both with and without geometrically local cubic terms exhibit average-case linear-time or weakly quadratic scaling when solved exactly using Gurobi, and the Ising models with no cubic terms show evidence of exponential-time TTS scaling when sampled using simulated annealing. These findings point to the necessity of developing and testing more complex, namely more densely connected, optimization problems in order for quantum computing to ever have a practical advantage over classical computing. Our results are another illustration that different classical algorithms can indeed have exponentially different running times, thus making the identification of the best practical classical technique important in any quantum computing vs. classical computing comparison.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Model Parameter Development for Complex Materials: Species-Specific Diffusion Barriers in 316 Stainless Steel from Systematic DFT Calculations

Vacancy-mediated diffusion barriers in 316 stainless steel have been systematically calculated using density functional theory to provide essential parameters for mesoscale microstructure evolution models. A statistical sampling approach employing 210 nudged elastic band calculations across multiple special quasi-random structures captures the effects of local chemical environments in this concentrated alloy. The computational methodology addresses challenges specific to chemically disordered systems, including proper magnetic treatment throughout multi-step calculations and validation against experimental structural properties. The calculated activation barriers reveal clear species-dependent diffusion behavior with the hierarchy Ni >> Fe ˜ Cr >> Mo. Nickel exhibits the highest barriers (0.74–1.31 eV, mean 1.045 eV), confirming its role as the slowest-diffusing major component. Iron and chromium show similar moderate barriers averaging 0.587 eV and 0.522 eV, respectively. Remarkably, molybdenum demonstrates exceptionally low barriers (0.12–0.28 eV, mean 0.194 eV), suggesting much higher mobility than previously recognized and potentially significant implications for precipitation kinetics and microstructure evolution. The barrier ranges remain consistent across different 316 SS compositions, supporting parameter transferability for modeling applications. The overall mean barrier of 0.64 eV provides a practical approximation for phase field simulations, while species-specific values enable detailed treatments of diffusion-controlled processes. This systematic approach establishes a validated framework for generating diffusion parameters in other concentrated alloys where experimental data are limited, while providing the first systematic set of species-specific barriers for predictive modeling of 316 stainless steel microstructure evolution.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Monitoring Potential Transport of Radioactive Contaminants in Shallow Ephemeral Channels: FY2019

Desert Research Institute (DRI) conducted a field assessment of the potential for contaminated soil to be transported from the Smoky Site Contamination Area (CA) because of storm runoff. This activity supported U.S. Department of Energy (DOE) Environmental Management Nevada Program (EM NV) efforts to establish post-closure monitoring plans for the Smoky Site Soils Corrective Action Unit (CAU) 550. The work was intended to confirm the likely mechanism of transport and determine the meteorological conditions that might cause the movement of contaminated soils, as well as determine the particle size fraction most closely associated with transported radionuclide-contaminated soils. These data will facilitate the design of the appropriate post-closure monitoring program. In 2011, DRI installed a meteorological monitoring station on the west side of the Smoky Site CA and a hydrologic (runoff) monitoring station within the CA, near the east side. The meteorological station collected air temperature, wind speed, wind direction, relative humidity, precipitation, solar radiation, barometric pressure, soil temperature, and soil water content data. The maximum, minimum, and average or total values (as appropriate) for each of these parameters were recorded for each 10-minute interval. The maximum, minimum, and average water depth in the flume installed at the hydrologic station were also recorded for each 10-minute interval. This report presents the data collected from these stations during fiscal year (FY) 2019. During the FY2019 reporting period, the warmest months were June, July, and August and the coldest were December and February. The highest monthly solar radiation (i.e., sunshine) values were recorded in June, July, and August, whereas December through February had the lowest solar radiation values. Monthly mean wind speeds were highest in the spring (April) and early summer (June). Winds were predominantly from the north, except in February, July, August, and September when winds were mostly from the west. The monthly average relative humidity ranged from 17 percent to approximately 60 percent. Humidity was lowest in the summer and highest in the winter. Monthly total precipitation ranged from 0.0 during August and September to 2.22 inches (in) (56.39 millimeters [mm]) in February. Total precipitation for FY2019 was 7.83 in (198.88 mm). During the reporting period, a single, nonzero, flow event was recorded at the flume. This occurred over two days between March 6 and 7, 2019, in response to 1.16 in (29.46 mm) of precipitation falling within 9.3 hours on March 6, 2019. This event produced a peak water depth of 0.39 in (0.99 centimeters [cm]) in the flume, and runoff was measured for approximately 23.5 hours. However, the maximum recorded flow depth was less than the limit for flow discharge and velocity that can be calculated for the six-inch (15 cm) throat of a Parshall flume. The runoff event produced insufficient flow velocities through the natural channel to cause local erosion and transport of bedload materials. No bedload samples were collected for radiological analysis. Observing meteorological and environmental conditions that lead to storm runoff will help identify the parameters and threshold conditions that should be incorporated into post-closure monitoring plans for this and similar Soils CAU sites. Routine monitoring of meteorological and hydrologic parameters at the Smoky Site CA has been discontinued at the request of EM NV, effective at the conclusion of FY2019. Environmental data acquisition was completed and the meteorological and runoff monitoring stations were removed from the Smoky Site on October 15, 2019.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Calibration and Localization of Optically Pumped Magnetometers Using Electromagnetic Coils

In this paper, we propose a method to estimate the position, orientation, and gain of a magnetic field sensor using a set of (large) electromagnetic coils. We apply the method for calibrating an array of optically pumped magnetometers (OPMs) for magnetoencephalography (MEG). We first measure the magnetic fields of the coils at multiple known positions using a well-calibrated triaxial magnetometer, and model these discreetly sampled fields using vector spherical harmonics (VSH) functions. We then localize and calibrate an OPM by minimizing the sum of squared errors between the model signals and the OPM responses to the coil fields. We show that by using homogeneous and first-order gradient fields, the OPM sensor parameters (gain, position, and orientation) can be obtained from a set of linear equations with pseudo-inverses of two matrices. The currents that should be applied to the coils for approximating these low-order field components can be determined based on the VSH models. Computationally simple initial estimates of the OPM sensor parameters follow. As a first test of the method, we placed a fluxgate magnetometer at multiple positions and estimated the RMS position, orientation, and gain errors of the method to be 1.0 mm, 0.2°, and 0.8%, respectively. Lastly, we calibrated a 48-channel OPM array. The accuracy of the OPM calibration was tested by using the OPM array to localize magnetic dipoles in a phantom, which resulted in an average dipole position error of 3.3 mm. The results demonstrate the feasibility of using electromagnetic coils to calibrate and localize OPMs for MEG.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

NNSS Soils Monitoring: Plutonium Valley (CAU 366) FY2019

Desert Research Institute (DRI) conducted a field assessment of the potential for radionuclide-contaminated soil to be transported from the Plutonium Valley Dispersion Sites Contamination Area (CA) because of both wind and storm water runoff. This activity supported U.S. Department of Energy (DOE) Environmental Management Nevada Program (EM NV) efforts to establish post-closure monitoring plans for the Plutonium Valley Dispersion Sites Corrective Action Unit (CAU) 366. The DRI task was intended to identify the likely mechanism(s) of transport and determine the meteorological conditions that might cause the movement of radionuclide-contaminated soils. The emphasis of the work was on collecting sediment transported by channelized storm runoff and measuring airborne dust concentrations and associated wind conditions. These data will facilitate an appropriate post-closure monitoring program. In 2011, DRI installed meteorological monitoring stations north and south of the Plutonium Valley CA—known as Pu Valley North and Pu Valley South, respectively—as well as a fluvial sediment sampling station within the CA. Since installation, temperature, wind speed, wind direction, relative humidity, precipitation, solar radiation, barometric pressure, soil temperature, volumetric soil moisture content, and airborne particulate concentrations have been collected at both meteorological stations. The maximum, minimum, and average or total (as appropriate) for each of these parameters were recorded for each 10-minute interval. The sediment sampling station included an automatically activated sampling pump with collection bottles for suspended sediment, which was activated when sufficient flow was present in the channel, and passive traps for bedload material that was transported down the channel during runoff events. This report presents data collected from these stations during fiscal year (FY) 2019. During the FY2019 (October 1, 2018, through September 30, 2019) reporting period, the warmest month was July and the coldest month was February. Monthly total solar radiation was highest in June or July (depending on the station) and lowest in December. Monthly mean wind speeds were highest during February. At the Pu Valley North station, winds were commonly northerly, northwesterly, and southerly. At the Pu Valley South station, the wind direction was more variable and winds frequently came from most directions, excluding westerly and northwesterly winds. Winds above 15.0 miles per hour (mph) (24.1 kilometers per hour [km/hr]) were frequently southerly at both stations, and less-frequent strong winds were northerly (at the Pu Valley South station) or northwesterly (at the Pu Valley North station). Monthly average relative humidity ranged from 18 percent in August to 62 percent in January and February. Monthly total precipitation ranged from zero at both stations in September to 2.23 inches (in) (56.6 millimeters [mm]) during February at the Pu Valley South station and 2.09 in (53.1 mm) during March at the Pu Valley North station. From October 1, 2018, through September 30, 2019, the total precipitation was 8.23 in (209 mm) and 8.15 in (207 mm) at the Pu Valley South and Pu Valley North stations, respectively. The largest daily precipitation totals of 1.46 in (37.1 mm) and 1.14 in (29.0 mm) occurred at the Pu Valley South and Pu Valley North stations, respectively, on March 6, 2019. However, the ISCO autosampler was not activated during any storm event in FY2019 because the channel water depth measurements were insufficient for suspended sediment sample collection. Additionally, because no significant flow through the channel was indicated, no bedload samples were collected. Light breezes of 0.0 mph to 5.0 mph (0.0 km/hr to 8.0 km/hr) occurred most frequently (approximately 55 percent to 57 percent of the time). The frequency of occurrence diminished exponentially as the wind speed increased such that winds in excess of 25.0 mph (40.2 km/hr) occurred less than 0.01 percent of the time. Winds in excess of 15.0 mph (24.1 km/hr) were most commonly southerly, which is likely controlled by the topography of the valley because topographic highs define the east and west sides of the valley and converge toward the north.The concentrations of PM2.5 (particulate matter with an aerodynamic diameter ≤2.5 micrometers [μm]) and PM10 (particulate matter with an aerodynamic diameter ≤10 μm) in the air generally increase as wind speed increases. However, the particle profiler at the Pu Valley South station regularly malfunctioned, which resulted in poor data quality and a limited period of record within FY2019. Therefore, discussion of the annual dust observations is mostly limited to the Pu Valley North station. At the Pu Valley North station, significant increases in windblown dust concentrations were observed when wind speeds exceeded 15.0 mph (24.1 km/hr), especially during southerly winds. When all wind speeds were considered, high dust concentrations at the Pu Valley North station were frequently associated with northerly, northwesterly, or southerly winds. The majority of PM10 transport at the Pu Valley North station occurred when winds were blowing northerly and northwesterly because of the high frequency of winds up to 15.0 mph (24.1 km/hr). The ratio of PM10 to PM2.5 is a qualitative indicator of the proximity of dust sources to the observation point and the values of the PM10 to PM2.5 ratio tend to be higher nearer to the source. The PM10 to PM2.5 ratio increased from 1.2 to 2.8 with higher wind speeds at the Pu Valley North station. This increase suggests that some locally sourced PM10 was present in the air when local winds exceeded 15.0 mph (24.1 km/hr). Local dust suspension is interpreted during major dust events when peaks in the PM10 to PM2.5 ratio occur simultaneously with the maximum wind speed and the peak PM10 concentration. Analysis of the five major dust events in FY2019 revealed that northerly winds were dominant in three events and southerly winds were dominant in two events. Local suspension of dust within Plutonium Valley is suspected to have occurred during two of the five major events (April 9, 2019, and September 28 to 29, 2019), and a nonlocal source of dust is presumed for the three other major dust events (October 14, 2018; May 16, 2019; and September 2, 2019). The three suspected regional dust transport events exhibited winds that typically did not exceed the 15.0 mph (24.1 km/hr) wind speed threshold for local dust transport and featured PM10 to PM2.5 ratios that indicate large-scale regional dust transport. Similar event characteristics (e.g., wind speed and PM10 concentration) were observed at air monitoring stations approximately 70 mi (113 km) to the northwest within the Tonopah Test Range (TTR). Therefore, these three regional events may not reflect transport from the Plutonium Valley CA. Routine monitoring of meteorological and hydrologic parameters at Plutonium Valley has been discontinued at the request of EM NV, effective at the conclusion of FY2019. Environmental data acquisition has been completed, and both Pu Valley North and Pu Valley South meteorological stations within Plutonium Valley as well as the ISCO autosampler station within the CA were removed on October 14, 2019.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

California Trees Seasonally Use Augmented Water Sources: Water Isotope Tracking in a Groundwater‐Dependent Ecosystem

Sustainable groundwater management must account for the needs of groundwater dependent ecosystems. To understand the relationship of ecosystems and seasonal water use, we studied the stable isotope composition (δ 18 O and δ 2 H) of water in streamside trees in a semi-arid streamside environment (Livermore, California, USA). We sampled seven trees at two sites every other month from April 2024 through April 2025 for tree xylem stable water isotope signatures. These data were compared to potential source waters: precipitation, imported surface water, soil water and regional groundwaters. Large daily precipitation events were found to be isotopically similar to regional groundwater and were thus treated as one water source. A Bayesian mixing model using stable water isotopes was used to determine the ratios of these three potential source waters (small daily precipitation events, groundwater/large daily precipitation events and imported water) present in tree xylem water. On average, tree water sources include 32% imported water (SD = 9%), 35% small daily precipitation events (SD = 10%) and 33% groundwater (SD = 3%), with significant seasonal variation (t summer-winter = 30.8, p < 0.01), particularly drawing more imported water (more than 55%) in the summer. While small daily precipitation events contribute only 10% of the total precipitation in our dataset, it represents a third of water used by trees. In addition, while these ecosystems are designated as groundwater dependent ecosystems, the trees use approximately one third imported water and even more during dry summer months. This approach provides water managers with a practical tool for quantifying ecosystem water needs, supporting data-driven decisions and regulatory compliance. While California's Sustainable Groundwater Management Act emphasizes supporting groundwater dependent ecosystems, it allows flexibility in demonstrating benefits. In conclusion, our methodology provides a way to document that management actions (such as managed aquifer recharge with imported water) deliver measurable co-benefits to GDEs.

Geosciences↗
DOE OSTI2021

Developing Upscaling Approach for Swarming Hydraulic Fractures Observed at Hydraulic Fracturing Test Site through Multiscale Simulations

We show that this work aims to address a challenge posed by recent observations of tightly spaced hydraulic fractures in core samples from the hydraulic fracturing test site (HFTS) in the Middle Wolfcamp Formation. Many fractures in retrieved cores have subfoot spacing, which is at odds with conventional models in which usually one hydraulic fracture is initiated per cluster. Models assuming a single fracture at each cluster, although a common practice, often predict excessive fracture propagation that is inconsistent with microseismic observation. Here, we aim to develop a numerical approach to effectively account for densely spaced hydraulic fractures in field-scale simulations. Because it is impractical to explicitly model all aforementioned fractures, we develop a new upscaling law that enables existing simulation tools to predict reservoir response to fracture swarms. The upscaling law is derived based on an energy equivalence argument and validated through multiscale simulations using a high-fidelity code, GEOS. The swarming fractures are first modeled with a spacing that is much smaller than the cluster spacing; these fractures are then approximated by an upscaled, single fracture based on the proposed upscaling law. The upscaled fracture is shown to successfully match the energy input rate and produce the total fracture aperture and average propagation length of the explicitly simulated swarm. Afterward, the upscaling approach is further implemented in 3D field-scale simulations and validated against the HFTS microseismic data of a horizontal well. Our results show that hydraulic fracture swarming can significantly affect fracture propagation behaviors compared with the propagation of single fractures as assumed by conventional modeling approaches. Under the considered situations, the conventional treatment yields fast propagation speed that far exceeds that indicated by the microseismic data. We also illustrate that this discrepancy can be reduced readily through the implementation of the upscaling law. Our results demonstrate the importance of accounting for the fracture swarming effect in field-scale simulations and the efficacy of this approach to enable realistic predictions of reservoir responses to fracture swarms, without the need to model tightly spaced fractures individually.

02 PETROLEUM↗
DOE OSTI2021

Computationally Efficient Multiscale Neural Networks Applied to Fluid Flow in Complex 3D Porous Media

Abstract The permeability of complex porous materials is of interest to many engineering disciplines. This quantity can be obtained via direct flow simulation, which provides the most accurate results, but is very computationally expensive. In particular, the simulation convergence time scales poorly as the simulation domains become less porous or more heterogeneous. Semi-analytical models that rely on averaged structural properties (i.e., porosity and tortuosity) have been proposed, but these features only partly summarize the domain, resulting in limited applicability. On the other hand, data-driven machine learning approaches have shown great promise for building more general models by virtue of accounting for the spatial arrangement of the domains’ solid boundaries. However, prior approaches building on the convolutional neural network (ConvNet) literature concerning 2D image recognition problems do not scale well to the large 3D domains required to obtain a representative elementary volume (REV). As such, most prior work focused on homogeneous samples, where a small REV entails that the global nature of fluid flow could be mostly neglected, and accordingly, the memory bottleneck of addressing 3D domains with ConvNets was side-stepped. Therefore, important geometries such as fractures and vuggy domains could not be modeled properly. In this work, we address this limitation with a general multiscale deep learning model that is able to learn from porous media simulation data. By using a coupled set of neural networks that view the domain on different scales, we enable the evaluation of large ( $$>512^3$$ > 512 3 ) images in approximately one second on a single graphics processing unit. This model architecture opens up the possibility of modeling domain sizes that would not be feasible using traditional direct simulation tools on a desktop computer. We validate our method with a laminar fluid flow case using vuggy samples and fractures. As a result of viewing the entire domain at once, our model is able to perform accurate prediction on domains exhibiting a large degree of heterogeneity. We expect the methodology to be applicable to many other transport problems where complex geometries play a central role.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Simultaneous Optical Transmission Spectroscopy of a Terrestrial, Habitable-zone Exoplanet with Two Ground-based Multiobject Spectrographs

Investigating the atmospheres of rocky exoplanets is key to performing comparative planetology between these worlds and the terrestrial planets that reside in the inner solar system. Terrestrial exoplanet atmospheres exhibit weak signals, and attempting to detect them pushes at the boundaries of what is possible for current instrumentation. We focus on the habitable-zone terrestrial exoplanet LHS 1140b. Given its 25-day orbital period and 2 hr transit duration, capturing transits of LHS 1140b is challenging. We observed two transits of this object, approximately 1 yr apart, which yielded four data sets thanks to our simultaneous use of the IMACS and LDSS3C multiobject spectrographs mounted on the twin Magellan telescopes at Las Campanas Observatory. We present a jointly fit white light curve, as well as jointly fit 20 nm wavelength-binned light curves from which we construct a transmission spectrum. Binning the joint white light-curve residuals to 3-minute time bins gives an rms of 145 ppm; binning down to 10-minute time bins gives an rms of 77 ppm. Our median uncertainty in R{sub p}{sup 2}/R{sub s}{sup 2} in the 20 nm wavelength bins is 260 ppm, and we achieve an average precision of 1.3× the photon noise when fitting the wavelength-binned light curves with a Gaussian process regression. Our precision on R{sub p}{sup 2}/R{sub s}{sup 2} is a factor of four larger than the feature amplitudes of a clear, hydrogen-dominated atmosphere, meaning that we are not able to test realistic models of LHS 1140b’s atmosphere. The techniques and caveats presented here are applicable to the growing sample of terrestrial worlds in the Transiting Exoplanet Survey Satellite era, as well as to the upcoming generation of ground-based giant segmented mirror telescopes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Automated co-adding and energy calibration of large array microcalorimeter data with zero sample knowledge

State-of-the-art microcalorimeter spectrometers now contain large detector arrays with hundreds of individual pixels. Each individual pixel outputs a unique and non-linear response with respect to deposited energy. This work describes a pattern-recognition algorithm to combine these responses into a single energy-calibrated histogram, referred to as co-adding pixels. Photo-peaks from different pixels are matched together based upon how well the match aligns the centroids and heights of neighboring peaks. This usually results in around 100 co-adding calibration points from 30 to 300 keV for a several day acquisition of plutonium items with masses between 0.5 and 10 grams. An additional algorithm energy-calibrates this co-added spectrum using the fluoresced K x-ray emissions from a tantalum absorber and inherent x-ray escape peaks from the tin absorbers. Both algorithms operate without knowledge of the source and are fully automated. This work presents results from the acquisitions of high and low burnup plutonium, 10% enriched uranium, a 153 Gd calibration source, and a 57 Co+ 166m Ho calibration source. In all measurements, resolution defined as the full-width at half-maximum (FWHM) of photo-peaks is preserved between the individual pixel and co-added spectra at around 65 eV for incident photon energies between 60 and 208 keV. The energy calibration algorithm is approximate and yields a calibration curve off by an average of around 200 eV for incident photon energies between 60 and 208 keV.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Deep learning to estimate permeability using geophysical data

Time-lapse electrical resistivity tomography (ERT) is a popular geophysical method to estimate three-dimensional (3D) permeability fields from electrical potential difference measurements. Traditional inversion and data assimilation methods are used to ingest this ERT data into hydrogeophysical models to estimate permeability. Due to ill-posedness and the curse of dimensionality, existing inversion strategies provide poor estimates and low resolution of the 3D permeability field. Recent advances in deep learning provide us with powerful algorithms to overcome this challenge. This paper presents a deep learning (DL) framework to estimate the 3D subsurface permeability from time-lapse ERT data. To test the feasibility of the proposed framework, we train DL-enabled inverse models on simulation data. Each measurement in both synthetic and field data is standardized by removing the mean and scaling the time-series to unit variance. This pre-processing step is necessary to bring simulation data closer to field observations. Subsurface process models based on hydrogeophysics are used to generate this synthetic data. Training performed on limited simulation data resulted in the DL model over-fitting. An advanced data augmentation based on mixup is implemented to generate additional training samples to overcome this issue. This mixup technique creates weakly labeled (low-fidelity) samples from strongly labeled (high-fidelity) data. The weakly labeled training data is then used to develop DL-enabled inverse models and reduce over-fitting. As both time-lapse ERT (1133048 features/realization) and 3D permeability (585453 features/realization) data samples are from a high-dimensional space, principal component analysis (PCA) is employed to reduce dimensionality. Encoded ERT and encoded permeability are generated using the trained PCA estimators. A deep neural network is then trained to map the encoded ERT to encoded permeability. This mixup training and unsupervised learning allowed us to build a fast and reasonably accurate DL-based inverse model under limited simulation data. Results show that proposed weak supervised learning can capture salient spatial features in the 3D permeability field. Quantitatively, the average mean squared error (in terms of the natural log) on the strongly labeled training, validation, and test datasets is less than 0.5. The R 2 -score (global metric) is greater than 0.75, and the percent error in each cell (local metric) is less than 10%. Finally, an added benefit in terms of computational cost is that the proposed DL-based inverse model is at least O(10 4 ) times faster than running a forward model once it is trained. Data generation, DL model training, and hyperparameter tuning to identify optimal neural network architectures utilized high-performance computing resources while the DL inference is performed on a standard laptop. Approximately, O(10 5 ) processor hours are used for generating data and DL tuning and training. We acknowledge that the data generation and DL model development are expensive. But once a DL model is trained, it can be re-used for inversion rapidly for the given system, with set physics and domain. Note that traditional inversion may require multiple forward model simulations (e.g., in the order of 10 to 1000), which are very expensive. This computational savings ≈ O(10 5 ) – O(10 7 )) makes the proposed DL-based inverse model attractive for subsurface imaging and real-time ERT monitoring applications due to fast and yet reasonably accurate estimations of permeability field.

58 GEOSCIENCES↗
DOE OSTI2021

Large-scale Gravitational Lens Modeling with Bayesian Neural Networks for Accurate and Precise Inference of the Hubble Constant

We investigate the use of approximate Bayesian neural networks (BNNs) in modeling hundreds of time delay gravitational lenses for Hubble constant (H 0 ) determination. Here, our BNN was trained on synthetic Hubble Space Telescope quality images of strongly lensed active galactic nuclei with lens galaxy light included. The BNN can accurately characterize the posterior probability density functions (PDFs) of model parameters governing the elliptical power-law mass profile in an external shear field. We then propagate the BNN-inferred posterior PDFs into an ensemble H 0 inference, using simulated time delay measurements from a plausible dedicated monitoring campaign. Assuming well-measured time delays and a reasonable set of priors on the environment of the lens, we achieve a median precision of 9.3% per lens in the inferred H 0 . A simple combination of a set of 200 test lenses results in a precision of 0.5 km s –1 Mpc –1 (0.7%), with no detectable bias in this H 0 recovery test. The computation time for the entire pipeline—including the generation of the training set, BNN training and H0 inference—translates to 9 minutes per lens on average for 200 lenses and converges to 6 minutes per lens as the sample size is increased. Being fully automated and efficient, our pipeline is a promising tool for exploring ensemble-level systematics in lens modeling for H 0 inference.

79 ASTRONOMY AND ASTROPHYSICS↗