Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interval analysis”

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

Freight Analysis Framework Version 5 (FAF5) Base Year 2017 Data Development Technical Report

The Freight Analysis Framework (FAF) integrates data from a variety of sources to create a comprehensive national picture of freight movements among states and major metropolitan areas by all modes of transportation. The latest of this data series is FAF5, which is the fifth generation FAF and is benchmarked on Commodity Flow Survey (CFS) 2017. Except for FAF1 that provided estimates for truck, rail, and water tonnage for calendar year 1998, later generations of FAF (FAF2 through FAF5) were built based on their benchmark year CFS data, for 2002, 2007, 2012, and 2017 respectively. The FAF is produced under a partnership between Bureau of Transportation Statistics (BTS) and Federal Highway Administration (FHWA). As a major data product of the FAF program, the FAF regional database provides a national picture of freight flows to, from, and within the United States (among regions and states), by commodity and mode for the base year, as well as for forecasts up to 30 years into the future in a 5-year interval. Additional FAF data products also include FAF network flows database, where truck movements are routed onto the national highway network, estimates of annual projections, and synchronized historical data series. This report is a technical document prepared to describe the data sources and methodologies applied in the process of building the FAF5 base-year 2017 regional database, released as FAF5.0 in February 2021. This report offers a description of the diverse data sources and modeling methods used in constructing the base year FAF5 regional database. The FAF5 base-year database is used as the base for development of forecasts and for assignment of truck flows on highway network. Similarly, the FAF5 base-year database will be used as the base to generate FAF5 annual estimates. In addition to this report, users are encouraged to refer to the FAF5 User’s Guide, which provides basic information of the data, including definitions of the data attributes, information on how to access the data and tool, as well as detailed data dictionary and code tables.

42 ENGINEERING↗
DOE OSTI2025

NOvA 2024 official data release (26.61E20 neutrino + 12.5E20 antineutrino)

This data release corresponds to the Bayesian 2024 analysis of NOvA $\nu_e$ appearance and $\nu_{\mu}$ disappearance data, corresponding to analysis described in https://arxiv.org/abs/2509.04361. File `NOvA_2024_data_histograms.root` contains data histograms for all the NOvA data samples. Exposure: * neutrino-enhanced beam: 26.61E20 protons on target * antineutrino-enhanced beam: 12.5E20 protons on target External constraints: * ss2th12=0.851, dm21=7.53e-5 are fixed at 2019 PDG values, with negligible effect on NOvA predictions. * ss2th13 & dm32: * RCDB1D: 1D constraint from Daya Bay for ss2th13, 0.0851+/-0.0024 * RCDB2D: Correlated 2D constraint from Daya Bay on ss2th13 & dm32, available in their official 2023 data release: https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.130.161802#supplemental An additional file containing the predictions (`NOvA_2024_prediction_with_systs_histograms.root`) for all channels’ signal and background components computed at the NOvA best-fit oscillation parameters and systematic pull terms. Numu samples include the no-oscillation case as well. Please note that any fits performed with these histograms are not expected to exactly reproduce the official NOvA results as parameterizations of the numerous systematic uncertainties considered in the official fits are not included in this release. The zip file also contains the 2D credible interval contours, with details in a README.md

Sztuc, Artur [University Coll. London] (ORCID:0000↗
DOE OSTI2025

Searching for Neutrino Tridents in the NOvA Near Detector

This dissertation presents a search for neutrino trident production in the NOvA near detector through the coherent ``dimuon" channel: $\nu_\mu +\hspace{1pt}\text{X} \rightarrow \nu_\mu + \mu^- + \mu^+ +\hspace{1pt}\text{X}$. Trident production is a rare, purely electroweak process with sensitivity to physics beyond the Standard Model. The theoretical background, motivation for studying the process, and previous experimental measurements are reviewed. The analysis uses data collected by the NOvA near detector (ND) from Fermilab's Neutrinos at the Main Injector (NuMI) beam between November 2014 and February 2024, corresponding to an exposure of $25.5\times 10^{20}$ protons on target. The ND is a segmented tracking calorimeter located 800~m from the beam target, receiving neutrinos with a mean energy of 2~GeV. A multi-pass background reduction strategy is implemented, including the development of a novel dimuon-specific tracking technique. Trident candidates are identified using a boost ed decision tree classifier trained on simulated signal and background events. Limited background Monte Carlo statistics necessitate the use of functional fits to sideband data, which are extrapolated to estimate backgrounds in the signal region. The unblinded data contain 9 trident-like events, with an estimated background of 5.66 $\pm$ 5.15 events. This yields a best fit estimate of 3.34 tridents compared to the Standard Model prediction of 4.66. A profiled Feldman-Cousins method is used to determine a 90\% confidence interval of [0,9.1] on the number of signal events, corresponding to an upper limit of 1.95$\times$ the Standard Model prediction. This result represents the lowest energy search for trident events to date, and the first experimental contribution to the process in 27 years.

Bowles, Reed Scott [Indiana U.]↗
DOE OSTI2020

The population of galaxies that contribute to the H i mass function

ABSTRACT We look at the contribution of different galaxy populations to the atomic hydrogen (H i) mass function (HIMF) and the H i density parameter, ΩH i, in the local Universe. Our analysis is based on a sample of 7857 H i selected galaxies selected from a volume common to the Sloan Digital Sky Survey (SDSS) and Arecibo Fast Legacy ALFA (40 per cent catalogue – α.40) surveys. We define different populations of galaxies in the colour (u − r)–magnitude (Mr) plane and compute the HIMF for each of them. Additionally we compute the HIMF for dark galaxies; these are undetected in SDSS and represent $\sim 2{{\ \rm per\ cent}}$ of the total sample. We find that the luminous red population dominates the total HIMF for $\log _{10}(M_{\rm{H\,{\small I}}}h^2_{70}/M_{\odot }) \ge 10.4$. The full red population – luminous and faint – represents about $\sim 17{{\ \rm per\ cent}}$ of the ΩH i budget, while that of the dark population is $\sim 3{{\ \rm per\ cent}}$. The HIMF about the knee, $\log _{10}(M_{\rm{H\,{\small I}}}h^2_{70}/M_{\odot }) \in [8,10.4]$, is dominated by the faint and luminous blue populations, the latter dominating at larger masses in this interval. Their total contribution to ΩH i is $\sim 55\!-\!70{{\ \rm per\ cent}}$, the variation depending on the definition of population. The dominant populations at the low-mass end, $\log _{10}(M_{\rm{H\,{\small I}}}h^2_{70}/M_{\odot }) \le 8.0$ are the faint blue and faint bluer populations, the latter’s dominance being sensitive to its definition. The full blue (blue–bluer luminous and faint) population represents $\sim 80{{\ \rm per\ cent}}$ of ΩH i. A bimodal HIMF suggested by our results is however not seen since the amplitude of the HIMF of the luminous red population is small compared to that of the luminous blue population.

Dutta, Saili↗
DOE OSTI2025

Lossy Compression: An Online Multi-Stage Technology for High-Fidelity Synchro- Waveform Measurements

Effective real-time monitoring and analysis of distributed grids necessitate the use of synchro-waveform measurements, which capture almost all high-frequency disturbances and transient phenomena. However, due to limitations in high-speed measurements and network bandwidth, it is challenging to transfer all high-fidelity synchro-waveforms losslessly and successfully. To cope with these challenges, a hybrid-based online multi-stage compression algorithm is proposed to significantly improve the compression efficiency for synchro-waveform measurements. Initially, the multiple discrete Wavelet transformation is deployed to deconstruct the waveform components. The delta encoding is further developed to decrease the magnitude. In conjunction with the Lempel-Ziv-Markov chain, the hybrid compression algorithm is implemented to achieve real-time compression for the synchro-waveform measurements. Moreover, an innovative error index that synergizes the time and frequency domain error and correlation is formulated to evaluate the waveform distortion. By integrating compression ratio, suitable parameters can be optimally selected. Finally, the simulation, laboratory experiments, as well as field tests across a spectrum of sampling frequencies and time intervals are conducted to substantiate the efficacy of the proposed method. Here, the outcomes demonstrated that a compression ratio of approximately 15.5 and 17.83 can be reached for 0.5 s and 1 s data under both offline and online scenarios, which equates to a substantial 93.5% to 94.39% reduction in data storage requirements.

High-fidelity synchro-waveform measurements↗
DOE OSTI2020

Best Practices for Shielding Analyses of Activated Metals and Spent Resins from Reactor Operation

This report investigated best practices for performing shielding evaluations of Type B waste packages, as defined in 10 CFR Part 71 regulations on packaging and transportation of radioactive material, or packages for which the contents are not defined before loading and may include a broad range of nuclides, geometries, and non-fuel materials. The following non-fuel waste streams were analyzed: Activated metals from decommissioned commercial power reactors, including Type 304 stainless steel, reactor vessel steel, and Inconel, Control blades from boiling water reactors, Neutron-activated corrosion products on surfaces of activated metals, and Spent resins from power plant operations. Measured elemental compositions, including major constituents and impurities, for steel and Inconel samples from commercial power reactors were used in activation calculations to determine radionuclide inventories in activated metals. For a simplified cask model, 60 Co contribution to the total external package dose rate at 30 days after shutdown varied from approximately 60% to 95%, depending on the activated metal, initial cobalt impurity concentration in the metal, and the thickness of the overpack gamma shield. Its maximum contribution to the total external package dose rate of approximately 100% was reached within the time interval of 2 to 5 years after shutdown and was maintained for up to 45 to 60 years after shutdown, depending on material, initial cobalt impurity concentration, and shield thickness. Thereafter, the 60 Co contribution to external package dose rate decreased with increasing decay time. Cobalt-60 is primarily produced by neutron reactions with the cobalt impurity in steel and Inconel. Other important radionuclides in activated metals contributing to package external dose rate are radionuclides with relatively short decay times, including 51 Cr, 59 Fe, 58 Co, and 54 Mn. These radionuclides may be represented as an equivalent 60 Co activity/source because 60 Co gamma ray emissions are bounding in terms of source strength and energy to other important radionuclides identified in the analyzed activated metals. Approaches for modeling the neutron-activated corrosion products that may be attached to activated reactor components were analyzed in this report. It was demonstrated that a surface source is more conservative than a uniform volumetric source for the treatment of neutron-activated corrosion products with respect to external package dose rates. An analysis of the maximum radionuclide loadings reported on spent resins identified the radionuclides 137 Cs, 60 Co, 134 Cs, 65 Zn, and 58 Co as the primary contributors to external package dose rate. For a resin cooled for 3.08 years, the external package dose rate was entirely produced by the reported 137 Cs and 60 Co inventory. The neutron sources from actinides found on spent resins or activated metals produced negligible dose rates and may be ignored in dose rate analyses. Effects of idealized waste material, source geometry, and spatial material/source distributions on external package dose rates were determined based on dose rate results for a simplified cask model under normal conditions of transport. Type 304 stainless steel, zirconium, and aluminum with adjusted mass densities based on a maximum content weight were analyzed for material modeling effects on external package dose rate. These materials produced identical external package dose rates within the statistical uncertainties of the dose rate estimates. Among four different source geometry configurations with homogeneous material of different mass densities, uniform volumetric source distribution, and the same source strength, the geometry configuration with lowest mass density (i.e., minimum self-shielding effect among the four cases) was most conservative. Spatial source distributions that better represent localized peak 60 Co activity values were more conservative than a uniform volumetric source distribution, assuming the same weight and total source strength per package. The increase in external dose rate caused by localized source peaks can be as much as the ratio between source peak activity density to average activity density, depending on the location of the activated metal with peak activity density. Therefore, the shielding analysis may be simplified if localized peak activities and the average activity per package can be measured/determined and documented at the time of cask loading. For simplicity, external package dose rates may be determined based on average source activity and a uniform volumetric source distribution. The dose rate results from that calculation model multiplied by the ratio of peak activity density to average activity density will produce maximum dose rate values for conservative estimates.

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

Sensitivity Analysis of Salt Storage Cavern Mechanical Integrity Test Parameters

This report examines the sensitivity of salt cavern Mechanical Integrity Tests (MIT) to uncertainties in key test parameters. MIT's are used by cavern operators to detect and quantify leak rates in access wells in underground salt storage caverns and involve a suite of measured and assumed parameters that have a direct impact on the sensitivity of the testing to detect actual leaks from the cavern storage system. Determining the sensitivity of the testing to these different parameters provides a basis for understanding the results from, and informing the design criteria for, this type of testing. Without fully understanding the sensitivity of the test to the testing parameters, it is possible that the test results may not accurately reflect the integrity of the cavern system; an actual leak may be missed, or an intact system may be interpreted as leaking. This report reviews the main parameters included in MITs and examines how selected changes in their values can impact test results. The deviations used in the sensitivity analyses were designed to be within the ranges believed to be similar to those which may be encountered during testing. The results show that small, plausible fluctuations in some of the parameters measured values can have a significant impact on the testing results. Of the parameters studied here, the sensitivity analyses showed the order of importance to be (from highest to lowest): nitrogen-oil interface depth measurement, well bore temperature, well head pressure, and finally the internal geometry of the testing interval.

42 ENGINEERING↗
DOE OSTI2024

Local-scale heterogeneity of soil thermal dynamics and controlling factors in a discontinuous permafrost region

In permafrost regions, the strong spatial and temporal variability in soil temperature cannot be explained by the weather forcing only. Understanding the local heterogeneity of soil thermal dynamics and their controls is essential to understand how permafrost systems respond to climate change and to develop process-based models or remote sensing products for predicting soil temperature. In this study, we analyzed soil temperature dynamics and their controls in a discontinuous permafrost region on the Seward Peninsula, Alaska. We acquired one-year temperature time series at multiple depths (at 5 or 10 cm intervals up to 85 cm depth) at 45 discrete locations across a 2.3 km 2 watershed. We observed a larger spatial variability in winter temperatures than that in summer temperatures at all depths, with the former controlling most of the spatial variability in mean annual temperatures. We also observed a strong correlation between mean annual ground temperature at a depth of 85 cm and mean annual or winter season ground surface temperature across the 45 locations. We demonstrate that soils classified as cold, intermediate, or warm using hierarchical clustering of full-year temperature data closely match their co-located vegetation (graminoid tundra, dwarf shrub tundra, and tall shrub tundra, respectively). We show that the spatial heterogeneity in soil temperature is primarily driven by spatial heterogeneity in snow cover, which induces variable winter insulation and soil thermal diffusivity. These effects further extend to the subsequent summer by causing variable latent heat exchanges. Finally, we discuss the challenges of predicting soil temperatures from snow depth and vegetation height alone by considering the complexity observed in the field data and reproduced in a model sensitivity analysis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Younger facial looks are associate with a lower likelihood of several age-related morbidities in the middle-aged to elderly

Abstract Background Looking older for one’s chronological age is associated with a higher mortality rate. Yet it remains unclear how perceived facial age relates to morbidity and the degree to which facial ageing reflects systemic ageing of the human body. Objectives To investigate the association between ΔPA and age-related morbidities of different organ systems, where ΔPA represents the difference between perceived age (PA) and chronological age. Methods We performed a cross-sectional analysis on data from the Rotterdam Study, a population-based cohort study in the Netherlands. High-resolution facial photographs of 2679 men and women aged 51.5–87.8 years of European descent were used to assess PA. PA was estimated and scored in 5-year categories using these photographs by a panel of men and women who were blinded for chronological age and medical history. A linear mixed model was used to generate the mean PAs. The difference between the mean PA and chronological age was calculated (ΔPA), where a higher (positive) ΔPA means that the person looks younger for their age and a lower (negative) ΔPA that the person looks older. ΔPA was tested as a continuous variable for association with ageing-related morbidities including cardiovascular, pulmonary, ophthalmological, neurocognitive, renal, skeletal and auditory morbidities in separate regression analyses, adjusted for age and sex (model 1) and additionally for body mass index, smoking and sun exposure (model 2). Results We observed 5-year higher ΔPA (i.e. looking younger by 5 years for one’s age) to be associated with less osteoporosis [odds ratio (OR) 0.76, 95% confidence interval (CI) 0.62–0.93], less chronic obstructive pulmonary disease (OR 0.85, 95% CI 0.77–0.95), less age-related hearing loss (model 2; B = −0.76, 95% CI −1.35 to −0.17) and fewer cataracts (OR 0.84, 95% CI 0.73–0.97), but with better global cognitive functioning (g-factor; model 2; B = 0.07, 95% CI 0.04–0.10). Conclusions PA is associated with multiple morbidities and better cognitive function, suggesting that systemic ageing and cognitive ageing are, to an extent, externally visible in the human face.

Mekić, Selma↗
DOE OSTI2021

Application of Point Precipitation Frequency Estimates to Watersheds

This report documents work sponsored by the U.S. Nuclear Regulatory Commission (NRC) at the Oak Ridge National Laboratory (ORNL) as part of the RES project, “Application of Point Precipitation Frequency Estimates to Watersheds.” This project was implemented as part of the Probabilistic Flood Hazard Assessment (PFHA) Research Program. The objective of the PFHA Research Program is to develop tools and guidance on the use of PFHA methods to risk-inform NRC’s licensing of new facilities as well as licensing and oversight of currently operating facilities as they relate to flooding hazards. Many nuclear power plants (NPPs) are located on or near rivers so riverine flooding hazards need to be considered in their design and operation. Probabilistic riverine flood models are important tools for realistic assessment of flooding risks. However, these models require areal estimates of the depth, duration, and frequency of rainfall distributed over the watershed, which are not often available. Point precipitation frequency estimates are more widely available. For example, the National Oceanic and Atmospheric Administration (NOAA) has published NOAA Atlas 14, which provides point precipitation frequency estimates for 5-minute through 60-day durations at average recurrence intervals of 1-year through 1,000-year. The research documented in this report addresses areal reduction factors (ARFs), which can be used to convert the widely available point precipitation frequency estimates, to estimates of areal precipitation frequency over a watershed. The most widely used ARF source is Technical Paper 29 (TP-29) published by the then U.S. Weather Bureau in 1958. However, both the methods and the underlying precipitation data used to produce TP-29 are seriously out of date. For example, due to the small gauge network available at the time of TP-29’s compilation, ARF estimates developed are only for watersheds smaller than about 400 square miles. Due to the relatively short record lengths of precipitation data available, frequency considerations could not be accurately determined. Other factors such as regional climate and seasonality were not addressed. Several newer methods have been published since TP-29 was developed and both the type and quantity of precipitation data have increased significantly, along with computational resources and analytical tools such as geographic information systems. This report reviewed and assessed the available precipitation products and methods for conducting ARF analysis. The work applied up-to-date precipitation data products and analysis methods with a novel watershed-based approach to investigate how ARF estimates vary across different methods, data sources, geographical locations, return periods, and seasons. The overall findings reported here regarding basic ARF trends are in line with other recent studies showing that ARFs decrease with increasing area, increase with increasing duration, and decrease with increasing return period. This study found significant differences among the available ARF methods. This work also found a strong geographical variability across different US hydrologic regions, suggesting that the ARF are specific to regional climate patterns and geographical characteristics and should not be applied arbitrarily to other locations. The results also reveal the importance of data record length, especially for high return level ARFs. The work reported in NUREG/CR-7271 will assist NRC staff in assessing different classes of ARF methods in conjunction with available rainfall data sets. It will also support the development of guidance for application of point precipitation data in PFHAs. It should be noted that the ARF values presented in this report for any location or region were developed for the purposes of comparing methods and investigating the factors that influence ARFs. They should not be considered official and should not be used in leu of a site-specific analysis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Uncertainty in Synthetic Tropical Cyclone Hazard and Risk Estimates: Insights from RAFT, CHAZ, MIT, STORM, and CLIMADA

We synthesize five complementary tropical cyclone (TC) hazard frameworks—RAFT (physics-based machine learning), CHAZ and MIT (statistical–dynamical), STORM (fully statistical), and CLIMADA (observation-driven resampling)—to characterize uncertainty in wind-related TC metrics relevant to energy applications. All datasets and the IBTrACS observational record are harmonized to a common 6-hourly, 2.5° grid. We compare basin-wide and coastal properties using consistent definitions for TC frequency, mean and maximum intensity, 24-hour intensification, and 6-hour translation speed, and quantify agreement with Pearson r, RMSE, and Kling–Gupta efficiency (KGE) alongside resampling-based confidence intervals. CLIMADA is included for basin context but excluded from coastal skill scoring because it resamples historical IBTrACS; if supplied with projected future tracks from an external hazard model, CLIMADA can be used to simulate future TC scenarios. Results show robust, cross-model signals: (i) a corridor of activity from the tropical Atlantic through the Caribbean into the Bahamas and western subtropical Atlantic; (ii) a meridional dipole in 24-hour intensification (low-latitude strengthening, subtropical weakening); and (iii) a transition from slower tropical motion to faster midlatitude translation. Coastal winds (mean and maximum) consistently cluster from the eastern Gulf into the Bahamas–western Atlantic transition. The largest structural spread occurs in the amplitude and footprint of lifetime maximum intensity and, secondarily, in translation speed; intensification exhibits similar central behavior across frameworks with variability in extremes. Translation speed shows the most uniform coastal agreement. These findings provide a decision envelope for wind-focused risk screening and clarify where uncertainty should be carried forward; wind-only results represent a lower bound on total hazard, motivating integration of surge and rainfall modules and a companion, asset-level damage analysis.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Understanding Solar Photovoltaic System Performance: An Assessment of 75 Federal Photovoltaic Systems

This report presents a performance analysis of 75 photovoltaic systems based on PV system production data collected as part of a FEMP Federal PV Performance Assessment project combined with co-incident insolation, and ambient temperature to analyze how actual performance compares with a performance model. FEMP collaborated with 17 Federal agencies and sub-agencies to collect the information required to analyze the performance of each system. The systems represent a total capacity of 30,714 kW and range in size from 1 kW to 4,043 kW, with an average size of 410 kW, and were installed between 2011 and 2020. The data is analyzed for Key Performance Indicators, Availability, Performance Ratio and Energy Ratio by comparing the measured production data to model production data. The System Advisor Model (SAM) combines a description of the system (such as inverter capacity, de-rating for temperature, balance-of-system efficiency) with environmental parameters (coincident solar and temperature data) to calculate predicted performance. The performance metrics are calculated by lining up the measured production data with the model estimate on an hour-by-hour, day-by-day, or month-by-month basis (depending on the interval resolution of the production data). A report with system description, photo of the system, special assumptions made for the site, graph of measured production and model production, table of key performance indicators, and links to O&M resources that might improve performance was produced and delivered to site and agency staff with a short on-line briefing.

14 SOLAR ENERGY↗
DOE OSTI2022

High-Resolution Wind Resource Data Set of the Greater Puerto Rico Region

In February 2022, the U.S. Department of Energy and six national laboratories launched the Puerto Rico Grid Resilience and Transitions to 100% Renewable Energy Study (PR100). PR100 aims to provide a comprehensive analysis of possible pathways for Puerto Rico's energy future, with a goal of 100% renewable energy by 2050. As a part of the renewable energy potential assessment in this project, we developed 20 years (2001-2020) of data using a numerical weather prediction (NWP) model for onshore and offshore wind resource assessment for the Puerto Rico. The research steps in developing the long-term wind resource data sets based on the NWP model were: 1. Model wind resource based on the Weather Research and Forecasting (WRF) model. 2. Develop WRF model configurations for Puerto Rico. 3. Test WRF with 11 different physics parameterizations for planetary boundary layer (PBL). 4. Assess WRF output from the different PBL schemes against observations. 5. Select a final model configuration which can produce the modeled wind speed with sufficient accuracy. 6. Produce 20 years wind resource data sets for Puerto Rico region. In the first stage of our framework for developing wind resource data, we developed a WRF model configuration using two nested domains (9 km and 3 km) to cover Puerto Rico and U.S. Virgin Islands and downscale the ERA5 reanalysis data (0.25 degrees x 0.25 degrees; hourly interval) to a 3-km domain. For the second stage, we implemented one-year simulations focused on using 11 different PBL physics parameterizations to find a combination of WRF physics parameterizations that could provide accurately modeled wind speed for Puerto Rico. We also analyzed the sensitivity of the modeled wind speed to PBL schemes for onshore and offshore locations. The WRF output resulting from the 11 WRF experiments using different PBL parameterizations were evaluated against observations obtained from the National Data Buoy Center (NDBC) as well as at hub height for a location for which measurements were available. A final model setup selected through the validation with observational data was used to produce 20 years of data with 3-km spatial and 5-minute temporal resolution. The WRF model output was post-processed to include wind profiles and basic atmospheric variables in a format that can be easily used for downstream modeling. The 20 years of wind resource data will be made available through NREL and support the estimation of wind energy development costs for the PR100 study.

17 WIND ENERGY↗
DOE OSTI2021

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

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

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Particle image velocimetry analysis with simultaneous uncertainty quantification using Bayesian neural networks

Particle image velocimetry (PIV) is an effective tool in experimental fluid mechanics for extracting flow fields from images. Recently, convolutional neural networks (CNNs) have been used to perform PIV analysis with accuracy on par with classical methods. Here we extend the use of CNNs to analyze PIV data while providing simultaneous uncertainty quantification on the inferred flow field. The method we apply in this paper is a Bayesian convolutional neural network (BCNN) which learns distributions of the CNN weights through variational Bayes. In order to demonstrate the utility of BCNNs for the PIV task, we compare the performance of three distinct BCNN models with simple architectures. The first network estimates flow velocity from image interrogation regions only. Our second model learns to infer velocity from both the image interrogation regions and interrogation region cross-correlation maps. Finally, our best performing network infers velocities from interrogation region cross-correlation maps only. We find that BCNNs using interrogation region cross-correlation maps as inputs perform better than those using interrogation windows only as inputs and discuss reasons why this may be the case. Additionally, we test the best performing BCNN on a full synthetic test image pair and a real image pair from the 1st International PIV Challenge. We show that ~98% of true particle displacements from the full synthetic image pair can be captured within the BCNN's 95% confidence intervals, and that the BCNN's performance on the real image pair is quantitatively similar to that of algorithms tested in the 1st International PIV Challenge. Finally, we show that BCNNs can be generalized to be used with multi-pass PIV algorithms with a moderate loss in accuracy, which may be overcome by future work on finetuning and training schemes. So to our knowledge, this is the first use of Bayesian neural networks to perform PIV.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Performance of Compact Pulsed Thermal Imaging System for In-Service Applications. Pulsed thermal tomography nondestructive examination of additively manufactured reactor materials and components

Additive manufacturing (AM) is an emerging method for cost-efficient fabrication of complex topology nuclear reactor parts from high-strength corrosion resistance alloys, such as stainless steel and Inconel. AM of metallic structures for nuclear energy applications is currently based on laser powder bed fusion (LPBF) process, which has the capability of melting metallic powder and net shaping the structures with relatively high precision. Some of the challenges with using LPBF method for nuclear manufacturing include the possibility of introducing pores into metallic structures. Integrity of AM structures needs to be evaluated nondestructively because material flaws could lead to premature failures due to creep in high temperature nuclear reactor environment. Currently, there exist limited capabilities to evaluate actual AM structures nondestructively. Pulsed Thermography (PT) imaging provides a capability for non-destructive evaluation (NDE) of sub-surface defects in arbitrary size structures. The PT method is based on recording material surface temperature transients with infrared (IR) camera following thermal pulse delivered on material surface with flash light. The PT method has advantages for NDE of actual AM structures because the method involves one-sided non-contact measurements and fast processing of large sample areas captured in one image. The data cube of PT measurements consists of surface temperature taken at sequential time intervals T(x,y,t). Material defects can be detected either by analyzing the thermograms T(x,y,t) data cube, or by using thermal tomography (TT) algorithm to obtain 3D spatial reconstruction of thermal effusivity e(x,y,z). To reduce the cost and enable in-service NDE in spatially constrained environment, it is highly desirable to develop PT with compact and inexpensive IR camera. Following initial qualification of an AM component for deployment in a nuclear reactor, a compact PT system can also be used for in-service nondestructive evaluation (NDE) applications. However, data cube obtained with PT based on compact IR camera suffers from strong thermal noises and loss of features due to relatively low sampling rate. In this report we describe two unsupervised machine learning (ML) algorithms for enhancement of PT images obtained with compact IR camera. In one approach, we introduce Sparse Coding Discrete Cosine Transform (SC/DCT) algorithm to remove additive white Gaussian noise (AWGN) from spatial thermal effusivity reconstructions. In another approach we introduce a Spatial Temporal Denoised Thermal Source Separation (STDTSS) ML algorithm to process thermograms. The STDTSS algorithm consists of spatial and temporal denoising using Gaussian and Savitzky–Golay filtering, followed by the matrix decomposition using Principal Component Analysis (PCA), and Independent Component Analysis (ICA) to automatically detect flaws. In the work described in this report, we constructed a compact PT system using a relatively small and low-cost FLIR A65 camera, consisting on uncooled microbolometer detector. Performance of SC/DCT algorithm was demonstrated on enhancing TT images of Inconel 718 AM plate. Performance of the STDTSS methods was investigated using thermography data obtained from imaging stainless steel 316L specimens produced with LPBF method with imprinted calibrated porosity defects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

SPRUCE Climate Warming and Elevated CO2 Rapidly Alter Peatland Soil Carbon Sources and Stability: Supporting Data

This data set reports a suite of complementary biogeochemical analyses of peat samples from the SPRUCE (Spruce and Peatland Responses Under Changing Environments) experiment. Results were collected using quantitative molecular analysis of bulk soil carbon to assess the stability of soil organic carbon following whole-ecosystem warming and exposure to elevated carbon dioxide concentrations (eCO2). Targeted soil organic carbon components include solvent-extractable compounds (alkanoic acids, alkanols, alkanes, steroids, and terpenoids), ester-bound hydrolysable biopolymers (cutin and suberin markers), lignin phenols, and pyrogenic carbon. Bulk peat samples were analysed by Soxhlet extraction and solid phase separation for solvent-extractable compounds, alkaline hydrolysis to extract hydrolysable biopolymers, copper (II) oxide oxidation to extract lignin phenols and benzene polycarboxylic acids (BPCAs) as an approximation of pyrogenic carbon. Samples were analysed by gas chromatography (GC) equipped with a flame ionization detector (GC-FID) and compound identification was performed on GC coupled to mass selective detector (MS) for solvent-extractable compounds, ester-bound hydrolysable biopolymers and lignin phenols, and high-performance liquid chromatograph (HPLC) for pyrogenic carbon. Results are presented in Ofiti et al. (accepted). The experimental work was conducted on samples collected in August 2018 at the SPRUCE climate manipulation experiment in northern Minnesota, 40 km north of Grand Rapids, in the USDA Forest Service Marcell Experimental Forest (MEF). Samples were collected and later analysed in a 10 cm increments over 0 to 50 cm depth and 25 cm intervals from 50 to 75 cm. Samples were analyzed for lignin phenols over 0 to 30 cm depth. This data set contains one file in comma separate (*.csv) format. This dataset contains data used to produce: Ofiti, N.O.E., Schmidt, M.W.I., Abiven, S., Hanson, P.J., Iversen, C.M., Wilson, R.M., Kostka, J.E., Wiesenberg, G.L.B., Malhotra, A. 2023. Climate warming and elevated CO2 rapidly alter peatland soil carbon sources and stability. Nat Commun 14, 7533. https://doi.org/10.1038/s41467-023-43410-z.

SPRUCE experiment, Marcell Experimental Forest, so↗
DOE OSTI2022

A multi-scale temperature-based strategy to map hydrologic exchange flows in highly dynamic systems

Mapping and quantifying hydrologic exchange flows (HEFs) is critical to environmental monitoring and remediation at contaminated sites; however, these objectives are challenging in highly dynamic systems, e.g., along dam-regulated rivers, where HEFs vary rapidly. Direct seepage measurements are labor-intensive and difficult to automate, whereas indirect (e.g., thermal) and remote sensing methods have potential to allow continuous monitoring with limited field effort. We present a preliminary assessment of a multi-scale temperature-based strategy for monitoring HEFs along the Hanford Reach of the Columbia River, in eastern WA, United States. Five thermal methods were assessed. First, a vertical temperature profile (VTP) was installed into the streambed. The VTP data were analyzed using a data assimilation algorithm designed for automated real-time estimation in dynamic systems. Second, a thermal infrared (TIR) camera was used in roving surveys to identify seeps. Third, a TIR camera was stationed at the VTP site to collect images at 1-h intervals. Together, the two TIR datasets provided a basis to assess the potential for drone-based TIR. Fourth, temperature was measured at the sediment/water interface to assess fiber-optic distributed temperature sensing. Fifth, imagery from the ECOSTRESS satellite mission was acquired to assess the potential of spaceborne thermal monitoring. Based on our preliminary assessment, VTP, TIR, and bed temperature measurements provide complementary spatial coverage, temporal sampling, and resolution; these methods have potential for long-term, automated monitoring of HEFs. The publicly available spaceborne imagery, however, proved inadequate because of insufficient spatial resolution and data gaps resulting from cloud cover and revisit frequency.

54 ENVIRONMENTAL SCIENCES↗