Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Control variates”

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 1,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Deep learning-based spatio-temporal estimate of greenhouse gas emissions using satellite data

Accurate estimation of greenhouse gases (GHGs) emissions is very important for developing mitigation strategies to climate change by controlling and reducing GHG emissions. This project aims to develop multiple deep learning approaches to estimate anthropogenic greenhouse gas emissions using multiple types of satellite data. NO2 concentration is chosen as an example of GHGs to evaluate the proposed approach. Two sentinel satellites (sentinel-2 and sentinel-5P) provide multiscale observations of GHGs from 10-60m resolution (sentinel-2) to ~kilometer scale resolution (sentinel-5P). Among multiple deep learning (DL) architectures evaluated, two best DL models demonstrate that key features of spatio-temporal satellite data and additional information (e.g., observation times and/or coordinates of ground stations) can be extracted using convolutional neural networks and feed forward neural networks, respectively. In particular, irregular time series data from different NO 2 observation stations limit the flexibility of long short-term memory architecture, requiring zero-padding to fill in missing data. However, deep neural operator (DNO) architecture can stack time-series data as input, providing the flexibility of input structure without zero-padding. As a result, the DNO outperformed other deep learning architectures to account for time-varying features. Overall, temporal patterns with smooth seasonal variations were predicted very well, while frequent fluctuation patterns were not predicted well. In addition, uncertainty quantification using conformal inference method is performed to account for prediction ranges. Overall, this research will lead to a new groundwork for estimating greenhouse gas concentrations using multiple satellite data to enhance our capability of tracking the cause of climate change and developing mitigation strategies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A comparative investigation of the effect of gas type on foam strength and flow behavior in tight carbonates

Foam generation technique has been practically applied to overcome gas mobility issue during enhanced oil recovery (EOR) processes. Incorporated with the utilized gas, the capability of foaming agent to generate persistent foam determines the flow performance under reservoir conditions. Here, in this study, we first identified the surfactants with excellent phase behavior and foaming properties and, second, evaluated their bulk foam performance with different gases, including N 2 , CO 2 , and hydrocarbon gas (CH 4 ) under high-pressure and high-temperature conditions. In bulk, the foam generated with the high-density gas, i.e., CO 2 , exhibited the greatest foamability but the lowest foam stability. With amidoamine oxide-based surfactant, CH 4 -foam was found to be more stable than N 2 -foam, while the opposite was observed for sulfobetaine-based surfactant. The most effective surfactant was later used to perform flow experiments on two tight carbonate rocks, including Minnesota Norther Cream Buff and Edwards limestone, with permeabilities of ~0.8 and 16 mD, respectively. The core-flooding experiments showed that the optimum quality of CO 2 -foam was lower than that of N 2 - and CH 4 -foams. For a given flow rate and fixed foam quality, the steady-state strength followed the descending order of CH 4 -, N 2 -, and CO 2 -foam. However, in general, the difference in foam strength between N 2 - and CH 4 -foams was insignificant when the total flow rate was varied. We also found that the variation in rock permeability had a greater effect on the strength of the CO 2 -foam compared to those generated with N 2 and CH 4 . In the presence of oil, the shear-thinning behavior of all foams was found more prominent. The increase in oil saturation had a more adverse effect on the strength of N 2 - and CH 4 -foams, particularly in the higher permeability rock. The results indicated that the foam strength, which is key in the effectiveness of foam EOR applications, is greatly influenced by several factors such as gas density, gas solubility in water, the presence of oil, and rock permeability.

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

Bio‐Imaging quorum sensing signal molecules in a soil‐mimic gel

Rhizosphere, the narrow region surrounding plant roots directly influenced by root exudates and the root associated microbiome, plays an important role in plant productivity and the rhizosphere is well known for stimulating microbial metabolic activities. How microbial communities interact to form stable, metabolically interconnected functional communities is an area of intense interest. The question remains on how microbially produced secreted molecules that function as intercellular communication signals shape the structure and function of microbial communities. Our goal is to detect and quantify signal molecules produced by microbes or plants in the root‐soil environment, spatially and temporarily. As a proof‐of‐concept, we are imaging diffusible extracellular microbial metabolites involved in a bacteria‐bacteria communication process called quorum‐sensing. Quorum‐sensing relies on the accumulation of high concentrations of signal molecules in the environment to control bacterial gene expression, influencing rhizosphere colonization and plant health. Local concentration of quorum sensing molecules are determined using aptamer based sensors. Aptamers that can specifically bind the desired signal molecules are selected through SELEX and immobilized on the surface of nano‐porous membranes. Binding of the signal molecules and aptamer covered surface results in changes in surface charge distribution and steric hinderance and thus modifying the transmembrane ionic transport. Changes in transmembrane impedance can be measured through electrochemical impedance spectroscopy methods to monitor local concentrations of signal molecules. Sensor responses were determined for different concentrations of signal molecules and results showed detection of C4‐HSL in a soil‐mimic solution with K D of 10 nM. We demonstrate that the quorum‐sensing signal molecule C4‐homoserine lactone and (C4‐HSL) can rapidly diffuse in a soil‐mimic gel, providing evidence for our use of a soil‐mimic for further aptasensor development and validation. The sensors were then inserted into soil mimic gel for monitoring C4‐HSL diffusion and the resulted impedance changes were used to determine the C4‐HSL concentration at different positions in the gel. Measurements of local C4‐HSL concentrations variation and numerical solution of diffusion equation were used to create 4D images of C4‐HSL molecule diffusion in the soil‐mimic gel.

Jiang, Nianyu↗
DOE OSTI2025

Multivariate environmental and trait-based controls of transpiration in the Central Amazon Rainforest

Tropical forest tree mortality is increasing due to more severe droughts, yet our understanding of how tree traits and life strategies are linked to drought stress has been limited by measurement scarcity. The BIONTE (BIOmass and NuTrient Experiment) near Manaus, Brazil hosts one of the world’s largest sap flow installations, with sensors in 90 canopy trees across a wood density gradient monitored since June 2022. The 2023 El Niño drought provided a unique opportunity to evaluate how water availability impacts tree transpiration. An interpretable machine learning framework was used to study the complex interactions between transpiration and multiple environmental variables such as soil water availability and vapor pressure deficit (VPD), and how these interactions vary with wood density and individual trees. We found varying responses of transpiration from different trees during the El Niño drought. Transpiration generally increased with temperature, with stronger effects in wetter areas and in trees with low to medium wood density. However, this response was modulated by stomatal sensitivity to VPD, which constrained transpiration under high atmospheric demand, particularly in intermediate-moisture area. The inflection in transpiration rate at high temperatures (>32°C) underscores the role of stomatal and hydraulic regulation in limiting water loss and protecting trees from excessive evaporative demand. Analysis of soil water contribution to transpiration revealed unimodal patterns in wetter area, with peak contributions near 0.45 cm 3 cm -3 of surface soil water and declining or flat responses beyond that threshold, suggesting a shift from water- to energy-limited transpiration. In contrast, drier areas exhibited limited transpiration sensitivity to soil water conditions and minimal trait-based variation in VPD responses, indicating supply-limited conditions. Despite higher wood density trees being generally more resilient, this study shows diverse tree drought resilience, prompting further investigation into the specific traits and dynamics between environmental variables in regulating transpiration and other physiological processes in trees.

Drought↗
DOE OSTI2020

Hydraulic characterization of a fault zone from fracture distribution

A quantitative assessment of how faults control the migration of geofluids is critical in many areas of geosciences. Here we integrated geological fieldwork, quantitative analysis of the fractures distribution and numerical modeling to build a geometrical representation of a fault zone and to characterize its hydraulic properties. Our target is a fault located in the Majella Mountain (Italy). We collected 21 scan lines across the fault profile in order to characterize its architecture. The numerical modeling of the fracture network of the damage zones and their hydraulic parameters was performed using both commercial (Move ® ) and open source software (dfnWorks and PFLOTRAN). Move ® was used to build a representative model of the fault zone using fracture spacing as a proxy, and to model the hydraulic parameters of the different fault domains. dfnWorks and PFLOTRAN were employed to infer the hydraulic parameters of the damage zones of the fault and then upscale these properties to an equivalent continuum domain, suitable for fluid flow simulations through the whole fault zone. Our findings show how even in a relatively small area it is possible to describe changes in terms of hydraulic properties of a fault zone and to build models capable to represent these variations.

58 GEOSCIENCES↗
DOE OSTI2023

Characterizing Current THD’s Dependency on Solar Irradiance and Supraharmonics Profiling for a Grid-Tied Photovoltaic Power Plant

The rapidly increasing distributed energy resources (DERs) in power systems are now getting interconnected to set community grid structures, where power quality will be a major concern. The grid-to-grid (G2G) bidirectional power transfer among the distribution microgrid will not be considered commercially feasible unless the upstream harmonics are under the limits. The aggregation of such harmonics, measured as total harmonic distortion (THD), is feared to be beyond tolerable limits with the progression of rooftop grid-tied PV-like installations. Hence, this THD needs to be characterized with DER generation end variables. In this work, the photovoltaic (PV) DERs’ dependency on environment variables such as irradiance was profiled in the context of generating and injecting harmonics into the grid. A mathematical model of a grid-tied three-phase PV DER was developed as part of this correlation characterization, matching the fundamental unit structure of a 1.4 MW solar canopy located on the Florida International University (FIU) Miami campus. To determine the qualitative association with produced THD patterns, the model was evaluated with various irradiance settings. A real-time digital simulation (RTDS) platform was used to verify it. Following this confirmation, sets of data from power quality meters at the point of common coupling and FIU field sensors were utilized to validate further the correlation model. The results showed that the grid current’s THD exhibited a high correlation with the irradiance profile and its variation over time. The early morning and late afternoon periods of the day, associated with a low irradiance, constantly had higher harmonics generated from the PV DER. The midday THD was rather rational with partial shadings, hence a geolocation-dependent factor. These findings were verified by an RTDS and validated by real field data. In quantifying the THD injected by a single DER at a high-frequency (2–150 kHz) supraharmonics (SH) level, a 3% peak increment in magnitude was observed from the high-fixed to the low-fixed irradiance profile. The correlation characteristics depicted that the hybrid microgrid suffered from a daytime-dependent harmonic insertion from the grid-tied DER. This is a global problem unless specific measures are taken to mitigate the harmonics. The electrically notorious higher-frequency SH was found to increase proportionally. The G2G power transfer can be limited because of the higher THD in the early morning and late afternoon, which will also worsen because the numbers of grid-tied PV DERs (i.e., rooftop solar and industrial solar) are likely to increase rapidly soon. The community grid structure can thus have a controlled harmonics filtration setup purposefully designed to address the findings of this work, which also fall within the scope of our future research.

14 SOLAR ENERGY↗
DOE OSTI2023

Modeling Short-Range Order in Disordered Rocksalt Cathodes by Pair Distribution Function Analysis

Pair distribution function (PDF) analysis is a powerful technique for the characterization of short-range order (SRO) in disordered materials. Accurate interpretation of experimental PDF data is critically reliant on the development of structural models that can account for local variations in site occupancies and bond lengths. To this end, we outline an approach to model SRO using first-principles calculations based on the cluster-expansion formalism. These methods are validated on neutron scattering data from two disordered rocksalt oxyfluorides, Li 1.3 Mn 0.4 Ti 0.3 O 1.7 F 0.3 and Li 1.3 Mn 0.4 Nb 0.2 Ti 0.1 O 1.7 F 0.3 . For each composition, we demonstrate that an average structure without any SRO fails to reproduce several key features in the experimental PDF. To pinpoint the origin of the suspected SRO in these materials, configurational and displacive effects were separately investigated using two disparate models. Special quasi-random structures were relaxed using density functional theory to account for local changes in bond lengths while maintaining a near-random ionic configuration. This leads to slightly improved accuracy but still misrepresents asymmetry in the first few peaks of the PDF. Monte Carlo simulations were performed to model configurational SRO on a fixed lattice, which by itself is shown to have a minimal influence on the PDF. Instead, we find that it is the bond length relaxations within environments created by SRO which controls the details of the PDF, thereby highlighting the subtle but important coupling between configurational and displacive SRO in disordered materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Nuclear Criticality Safety Assessment of Criticality Control Containers without Moderation Control at the Waste Isolation Pilot Plant

The Waste Isolation Pilot Plant (WIPP) provides for safe, permanent disposal of government-owned transuranic (TRU) and TRU mixed wastes. Receipt and disposal of waste at the WIPP site began in March 1999. The Sandia report, Consideration of Nuclear Criticality When Disposing of Transuranic Waste at the Waste Isolation Pilot Plant, addressed potential nuclear criticality safety issues based on the projected inventory characteristics known at the time [1]. As designs for inventory, waste forms, and disposal packages have changed, new analyses have been performed, and updates have been made to address any potential effects to the WIPP safety basis. New analyses performed include Saylor 2017 [2] and Brickner 2019 [3], which address certain waste containers with specified loadings under post-closure conditions. Both examined several hypothetical scenarios and included analyses to bound (from a criticality potential standpoint) credible configurations that could occur at WIPP during the repository regulatory post-closure disposal time period for feature, event, and process (FEP) considerations—10,000 years. During this post-closure period at WIPP, the screening of FEPs is governed by the risk-based standards and implementing regulations of the US Environmental Protection Agency (EPA) (i.e., 40 CFR 191 and 40 CFR 194, respectively) [4,5]. An FEP screening can be based on either a low-consequence or low-probability rationale. A low-probability rationale includes either (a) a qualitative rationale that the FEP is not credible or (b) a quantitative demonstration that the probability is less than 10-4 in 104 years. In this evaluation, a qualitative lowprobability rationale of not credible is used by demonstrating that bounding configurations of the waste are not critical. The demonstration of subcriticality is through quantitative calculations, but a probability of criticality is not evaluated. Rather, the rationale for this evaluation is that bounding configurations with an effective neutron multiplication factor (keff) well below the upper subcriticality limit (USL) make criticality incredible. Reference [2] documented a nuclear criticality assessment of the WIPP repository for disposal of dilute surplus plutonium materials using the Dilute and Dispose Approach and packaging in criticality control overpacks (CCOs). The CCO is the waste disposal container recently designed to allow for up to 380 fissile gram equivalent (FGE) 239 Pu per drum, which is a higher fissile loading than typical waste containers. The CCO consists of a criticality control container (CCC) positioned by upper and lower plywood spacers within a standard 55 gal drum. The CCC is used to establish a geometry control for fissile materials during transportation and WIPP emplacement operations. The current WIPP waste acceptance criteria for CCO payloads limit beryllium to less than or equal to 1% by weight of the waste contents and require the waste form to be non-machine compacted. Reference [2] considered two scenario progressions—room closure from salt creep, hereafter referred to as the reconfigured dry scenario, and flooding with brine, hereafter referred to as the reconfigured wet scenario. The subsequent drying out of the reconfigured wet scenarios was also considered. For all scenarios, subcriticality was maintained when 50 g of B 4 C (acting as a neutron absorber) per CCC was intermixed within the plutonium disposition waste form. The analysis used a waste form description that limits the amount of moderation that could be present within the waste form (i.e., it limits the amount of water and polyethylene that could be present based on planned processing conditions). This analysis to evaluate increased limits on the amount of moderation that could be present was performed as a companion to Reference [2] to address concerns associated with verifying moisture and/or plastic contents of waste materials following packaging of dilute surplus plutonium in the CCO. To that end, this analysis used the models and methods from Reference [2] to evaluate a more generic base waste form consisting of water and polyethylene that is more similar (and nearly identical) to the generic waste forms utilized in other models/analyses supporting the TRU Package Transporter Model II (TRUPACTII) safety analysis [6] (all are without moderation controls). The waste form in this analysis uses a base mixture of 75% water and 25% polyethylene, the total amount of which is varied to determine the optimum moderation to fissile material (H/Pu) ratio. The fissile loading is maintained at up to 380 FGE 239 Pu (modeled as PuO 2 ) per CCO with an additional 545 g of beryllium (to bound the 1% by weight contents restriction) and 50 g of B 4 C intermixed per CCO. The beryllium content (1% by weight) is based on the total allowed waste weight (this does not include packaging and container weights). Figures ES-1 and ES-2 display summary results, showing that with this model including 50 g of B 4 C per CCO, the system keff remains under 0.85 for all moderator amounts and provides a significant margin against post-closure criticality under postulated bounding conditions for compaction. Figure ES-1 compares an infinite model with a room model at the initial emplacement spacing and under full radial compaction. Full radial compaction places each CCC in direct contact and does not credit any anticipated spacing associated with current post-closure geomechanical modeling of the repository [7]. The effects of variations in the H/Pu ratio were evaluated by varying the amount of the water/polyethylene component of the waste model, with fissile loading maintained at 380 239 Pu FGE. Similarly, Figure ES-2 illustrates how various amounts of B 4 C per CCO influence k eff at different radial compactions, all at the H/Pu ratio of 200 (in the room array model). Therefore, while the results from Saylor 2017 [2] modeled more realistic process limits associated with packaging of dilute surplus plutonium, this analysis demonstrates that limits on moderation (plastic and water content) are not necessary to ensure subcriticality in the WIPP repository, provided the requisite B 4 C absorber is present.

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

Surface and Buried Thermal, and RGB Unexploded Ordnance Data Collection

This document provides a description of a data collection campaign of unexploded ordnance (UXOI) set. The dataset captures a controlled UAV imaging campaign designed to support detection of UXO across varied environmental conditions. Data were collected during three campaigns in Norris and Northeast Knoxville, Tennessee, using RGB, and thermal sensors mounted on Parrot UKR. In total, the dataset contains 9925 images, 26 full-motion video, and approximately 81.99 GB of data, collected across late spring/summer conditions, every hour during sunlight, and multiple surface contexts, including tall grass, short grass, gravel, as well as buried in sand, and other gravel mixtures. The collection was designed to capture thermal and visual variability relevant to UXO detection in agricultural land, bare earth, and subsurface. Review of the imagery showed that ordnance was most detectable during periods of changing solar input, especially approximately 10-60 minutes after sunrise, approximately 20-60 minutes after sunset, and 2-3 min after cloud cover interrupted prolonged solar heating. These conditions increased thermal contrast because many ordnance items retained or released heat differently than the surrounding vegetation and ground surface. This dataset provides a useful resource for developing and evaluating airborne UXO detection methods under realistic field conditions. All ordnance used in the study was inert, and thermal behavior may differ from that of live ordnance. In addition, variation in ordnance type, composition, and placement introduced differences in thermal response that should be considered when interpreting results.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2026

Validation of an integrated modeling framework for investigating 3D plasma responses in tokamak plasmas

As contemporary experimental tokamaks are pushed toward reactor-relevant operation, they provide essential testbeds for demonstrating and exploring ELM control strategies for deployment in future fusion pilot plants. Accurate predictions of full plasma responses are essential to guide and optimize these demonstrations. This paper introduces and validates an integrated modeling framework over a historical range of DIII-D operational space. The integrated modeling framework uses only scalar plasma parameters and optional reference boundary to self-consistently and flexibly scan through tokamak operational space and estimate the corresponding ELM-suppression relevant plasma response over 3D coil phase space. The framework generates tightly converged equilibria that satisfy a target set of plasma parameters (I p , β N , l i ), with kinetic profiles constrained by an EPED(NN)-computed pedestal and empirical core model. The plasma response of these modeled equilibria is calculated with GPEC. Validated against 55 distinct DIII-D equilibria, the framework consistently reproduces experimentally constrained equilibrium pressure, q, and other representative profiles using the time-varying information from only evolving scalar plasma parameters; plasma response validation is performed for the n=3 perturbation. In particular, the plasma response of a DIII-D discharge scanning q 95 in search of ELM suppression windows is quantitatively reproduced. By accurately reproducing experimental equilibria and plasma responses across wide parameter variation, the framework supports potential for synthetic parameter scans in key operational, stability, and plasma response dimensions to investigate RMP ELM-suppression experiments and inform predictive RMP scenario optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Differences in substrate use linked to divergent carbon flow during litter decomposition

ABSTRACT Discovering widespread microbial processes that create variation in soil carbon (C) cycling within ecosystems may improve soil C modeling. Toward this end, we screened 206 soil communities decomposing plant litter in a common garden microcosm environment and examined features linked to divergent patterns of C flow. C flow was measured as carbon dioxide (CO2) and dissolved organic carbon (DOC) from 44-days of litter decomposition. Two large groups of microbial communities representing ‘high’ and ‘low’ DOC phenotypes from original soil and 44-day microcosm samples were down-selected for fungal and bacterial profiling. Metatranscriptomes were also sequenced from a smaller subset of communities in each group. The two groups exhibited differences in average rate of CO2 production, demonstrating that the divergent patterns of C flow arose from innate functional constraints on C metabolism, not a time-dependent artefact. To infer functional constraints, we identified features – traits at the organism, pathway or gene level – linked to the high and low DOC phenotypes using RNA-Seq approaches and machine learning approaches. Substrate use differed across the high and low DOC phenotypes. Additional features suggested that divergent patterns of C flow may be driven in part by differences in organism interactions that affect DOC abundance directly or indirectly by controlling community structure.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Arctic extreme precipitation changes from 1980 to 2022 in response to sea ice decline and enhanced atmospheric rivers

Arctic extreme precipitation (EP) broadly impacts permafrost degradation, glacier and snow cover changes, and ice sheet mass balance as well as ecosystems. However, investigation of EP spatiotemporal variations over the Arctic remains challenging, and their primary drivers are still poorly understood. Performance estimation of three state-of-the-art reanalysis products (Climate Forecast System Reanalysis (CFSR), European Centre for Medium-Range Weather Forecasts Reanalysis version 5 (ERA-5), and Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2)) against gauge-based precipitation observations reveals that MERRA-2 outperforms other reanalysis for Arctic EP changes. Based on MERRA-2 data, both annual EP amount and occurrence days averaged over the Arctic show statistically significant positive trends during 1980–2022 (3.37 ± 1.03 mm dec −1 and 0.42 ± 0.17 d dec −1 , respectively), with the most pronounced increase in the autumn. Spatial heterogeneity in annual and seasonal EP trends is found across the Arctic, with the largest positive annual trends of 30 mm dec −1 over the Bering Sea and the Denmark Strait. The significant EP increase is closely associated with intensified atmospheric rivers (ARs) and widespread decline in sea ice concentration (SIC). Specially, SIC and ARs are responsible for 12% and 50% of Arctic EP inter-annual variance, respectively, while ARs directly contribute 28.3% of the total annual EP amounts. These findings explain the mechanistic controls on Arctic EP, providing critical insights for projecting polar weather and climate extremes and their impacts on the Arctic environment.

atmospheric rivers↗
DOE OSTI2024

Seasonal Cost-Benefit Analysis of Automated Distribution Feeder Upgrades with Advanced Mitigation Technologies

The increasing deployment of distributed solar photovoltaics (DPV) to meet clean energy goals can trigger adverse grid operation issues, such as voltage excursions and the violation of thermal loading constraints of the power delivery elements (e.g., lines and transformers) on the evolving electricity infrastructure. Such integration issues would require distribution upgrades with associated costs to mitigate them and to maintain reliable and resilient grid operating conditions. Traditional distribution network upgrade approaches use a specific single snapshot analysis that is overly conservative. This study considers a multi-time point analysis to capture both moderate (probable bounds) and extreme grid operating conditions using time points such as minimum load with minimum photovoltaics (PV), maximum load with maximum PV, maximum load with minimum PV, and minimum load with maximum PV. Further, this study investigates seasonal variation impacts and associated distribution upgrade costs for a spring season case (March, representing a low load and high PV scenario) and a summer case (July, representing a high load and high PV scenario). Such seasonal analysis will allow system operators to characterize upgrade requirements and associated costs across various periods. Because the spatial distribution of DPV can impact upgrade and associated costs, this study investigates three common DPV deployment scenarios - randomly deployed, close to the substation, and far from the substation - at different penetration levels. Apart from spatial distribution impacts, this project evaluates the techno-economic impacts of the nodal photovoltaic penetration factor (NPPF) for generating the various DPV deployment scenarios at increasing penetration levels. This project investigates the impact of varying nodal PV-to-load ratios using conservative and extreme NPPF values of 3 and 10, respectively. This study investigates the deployment of traditional infrastructure upgrade strategies, such as installing new voltage regulating equipment, transformers and lines replacements, and the activation of advanced inverter functionality (e.g., autonomous volt/VAR) in expanding PV hosting capacity. Existing DPV systems are assumed to operate with the legacy unity power factor, and we considered the possibility of retrofitting such systems with the activation of volt/VAR control as integration standards and regulations continue to evolve to allow such functions. The cost-benefit analysis metrics used in study include distribution upgrade costs, average cost per watt of the upgrade cost, average marginal cost per watt of the upgrade cost, and power losses.

14 SOLAR ENERGY↗
DOE OSTI2024

Large reductions in Permian Basin methane intensity shown in multi-year comparison of aerially-visible methane emissions

Spanning the US states of Texas and New Mexico, the Permian Basin has been a hotspot of methane emissions from oil and natural gas activity 1–4, although studies disagree over the magnitude of these emissions. The most comprehensive measurement campaigns published were conducted in 2019 1–5 .There have been large changes in the energy industry since then, including in the prices of oil and gas, both state and federal regulatory environments, investor and activist pressure over methane emissions, and the adoption of new technologies and policies by energy operators. Understanding how any or all of these might influence methane emissions is important for policy makers, oil and gas operators, and other stakeholders. We characterize the time evolution of Permian Basin methane emissions using a series of comprehensive aerial surveys conducted every year from 2020-2023 and compare them to the 2019 results cited above. To maintain comparability, all the data sets are from surveys using Insight M point source methane sensing technology. The scope of these surveys expanded over time: from 33-46% of wells, oil production, and gas production in 2020 to 60-65% in 2021, to 84% of wells and over 90% of both oil and gas production in 2023. These surveys by Insight M also include hundreds of gas processing plants and compressor stations as well as 1000s of km of gathering and transmission pipelines. Considering only the aerially detected portion of emissions (typically the majority of the total in such surveys 4), we find reductions of more than 70% in methane emissions intensity compared to the 2019 New Mexico-only Insight M survey, with variation depending on the year 3,4. Notably, although sources below 100 kg/hr contributed less than 10% of aerially measured emissions the 2019 New Mexico survey 3, these smaller sources constitute a larger proportion of total aerially measured emissions (although not the majority) in 2020-2023. Production facilities and gathering pipelines are responsible for the larges shares of total emissions, followed by compressor stations and gas processing plants. Permian methane emissions were also measured in a comprehensive 2019 Permian-wide survey by the Carbon Mapper team 2. That analysis led to a lower total emissions estimate at the time 4. These new Insight M-based emission rates are still roughly 30-70% lower than the aerially measured portion of the 2019 Carbon Mapper-based estimates 4. Further work is needed to harmonize these surveys in space and time to create the most intercomparable numbers possible 5. Additional analysis is needed to compare our findings to the more spatially constrained 2020, 2021, and 2023 Carbon Mapper surveys in the Permian 4,6. The evidence is strong from these two survey teams that emissions intensity has declined significantly since 2019. Reasons for this trend are currently unclear but point to possible success of emissions control programs. Future work investigating frequency, source, and operator-specific intensities could provide insights into the causes of this promising trend.

methane, oil and gas, data science, remote sensing↗
DOE OSTI2026

Monitoring of Liquid Metal Reactor Heater Zones with Recurrent Neural Network Learning of Temperature Time Series

Advanced high-temperature fluid reactors (ARs), such as sodium fast reactors (SFRs) and molten salt cooled reactors (MSCRs) utilize high-temperature fluids at ambient pressure. To melt the fluid during reactor startup and prevent fluid freezing during cooldown, the thermal–hydraulic systems of such ARs include heater zones consisting of specific heaters with controllers, temperature sensors, and thermal insulation. The failure of heater zones due to insulation material degradation or improper installation, resulting in parasitic heat losses, can lead to fluid freezing. The detection of faults using a heat-transfer model is difficult because of a lack of knowledge of the experimental details. Data-driven machine learning of heater zone temperature time series offers a viable alternative. In this study, we benchmarked the performance of recurrent neural networks (RNNs) in an analysis of heat-up transient temperature time series of heater zones installed on a liquid sodium vessel. The RNN models include long short-term memory (LSTM) and gated recurrent unit (GRU) networks, as well as their bi-directional variants, BiLSTM and BiGRU. Anomalous temperature points were designated using a percentile-based threshold applied to residual fluctuations in the detrended temperature time series. Additionally, the impact of the exponentially weighted moving average (EWMA) method on detection accuracy was examined. The RNN models’ performance was assessed using precision, recall, and F 1 score metrics. Results demonstrated that RNN models effectively detect anomalies in temperature time series with the best models for each heater zone achieving F 1 scores of over 93%. To explain the variations in RNN model performance across different heater zones, we used Kullback–Leibler (KL) divergence to quantify the relative entropy between training and testing data, and the Detrended Fluctuation Analysis (DFA) to assess long-range temporal correlations. For datasets with strong long-range correlations and minimal relative entropy between training and testing data, GRU is the best-performing model. When the data exhibits weaker long-term correlations and a significant relative entropy between training and testing distributions, BiGRU shows the best performance. For the data sets with intermediate values of both KL divergence and DFA, the best performance is obtained with LSTM and BiLSTM, respectively.

gated recurrent unit↗
DOE OSTI2026

Effects of Gasoline Composition and Operating Parameters on the Response of End-Gas Autoignition to Nitric Oxide in a Lean-Burn DI-SI Engine

Lean operation of spark-ignition engines can lead to engine thermal efficiency gains and lower NOx emissions due to reduced combustion temperatures. Yet, lean operation could still face challenges in end-gas autoignition and knock generation due to higher intake pressures and trapped NO in the residual gas. Here, this study evaluates the impact of NO on end-gas autoignition for two gasoline fuels with similar octane rating but different composition: high cycloalkane fuel (HCA) and high olefin fuel (HO). Experiments were performed at stoichiometric and lean (λ = 2) conditions and at two engine speeds of 1400 rpm and 2000 rpm. Accompanying chemical kinetics simulations in CHEMKIN revealed that the mechanisms controlling the effect of NO on autoignition are similar λ = 2 and λ = 1, with NO + HO 2 = NO 2 + OH being the main pathway for enhancing reactivity by promoting low-temperature heat release (LTHR). The compositionally different fuels reacted differently to NO seeding and engine speed, and differences were augmented at λ = 2 compared to λ = 1 as the end-gas autoignition shifted to the low temperature regime. HO, which has inherent low temperature chemistry, was strongly impacted by engine speed at low NO seeding levels, with no noticeable peak of LTHR detected at 2000 rpm. On the other hand, LTHR of HCA was marginally affected by shortened residence time at higher engine speed as NO + HO 2 reaction was not greatly affected by shorter time scales, since HO 2 production was sustained even at 2000 rpm to support OH generation from NO + HO 2 . Contrary to HO, HCA exhibited greater sensitivity to NO seeding, as the increased OH production at higher NO concentrations offset the OH-quenching effect of cyclopentane, which accounts for 28.6% of HCA’s composition. Consequently, a sensitivity analysis revealed that fuels with weak inherent low-temperature chemistry, like HCA, are likely to be more sensitive to variations in NO concentration and charge temperature, whereas fuels with strong low temperature chemistry are more sensitive to variations in end-gas λ and intake pressure.

Energy↗
DOE OSTI2022

Intra-canopy leaf trait variation facilitates high leaf area index and compensatory growth in a clonal woody encroaching shrub

Leaf trait variation enables plants to utilize large gradients of light availability that exist across canopies of high leaf area index (LAI), allowing for greater net carbon gain while reducing light availability for understory competitors. While these canopy dynamics are well understood in forest ecosystems, studies of canopy structure of woody shrubs in grasslands are lacking. To evaluate the investment strategy used by these shrubs, we investigated the vertical distribution of leaf traits and physiology across canopies of Cornus drummondii, the predominant woody encroaching shrub in the Kansas tallgrass prairie. We also examined the impact of disturbance by browsing and grazing on these factors. Our results reveal that leaf mass per area (LMA) and leaf nitrogen per area (Na) varied approximately threefold across canopies of C. drummondii, resulting in major differences in the physiological functioning of leaves. High LMA leaves had high photosynthetic capacity, while low LMA leaves had a novel strategy for maintaining light compensation points below ambient light levels. The vertical allocation of leaf traits in C. drummondii canopies was also modified in response to browsing, which increased light availability at deeper canopy depths. As a result, LMA and Na increased at lower canopy depths, leading to a greater photosynthetic capacity deeper in browsed canopies compared to control canopies. This response, along with increased light availability, facilitated greater photosynthesis and resource-use efficiency deeper in browsed canopies compared to control canopies. Furthermore, our results illustrate how C. drummondii facilitates high LAI canopies and a compensatory growth response to browsing—both of which are key factors contributing to the success of C. drummondii and other species responsible for grassland woody encroachment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Accuracy of HVAC Load Predictions: Validation of EnergyPlus and DOE-2 using FLEXLAB Measurements

The aim of the project reported here was to better understand the level of accuracy of three building energy simulation (BES) engines (‘engines’) — EnergyPlus™, DOE-2.1e, and DOE-2.2 — by identifying and investigating significant deviations between the performance predicted by these engines and actual performance as measured in the FLEXLAB® test facility at Lawrence Berkeley National Laboratory (LBNL). The specific test conditions included some of those prescribed in ANSI/ASHRAE Standard 140 - Standard Method of Test for the Evaluation of Building Energy Analysis Computer Programs. Detailed measurements of FLEXLAB performance, including indoor temperatures and heat fluxes and air-flow and water flow rates and temperatures in the Heating, Ventilating and Air Conditioning (HVAC) system, together with hourly weather data, were recorded and used in analyzing the simulation results from EnergyPlus v8.8, DOE-2.2 v3.65 and DOE-2.1e v127. These engines are commonly used in the United States for building energy code compliance, federal, state, and utility incentives programs, as well as energy efficient design of new buildings and energy retrofit of existing buildings. Seven conventional overhead mixing ventilation scenarios were tested and each engine was found to have a similar level of agreement with the measurements of space-level heating and sensible cooling loads. These results provide useful information regarding the accuracy of these engines in predicting the cooling and heating load elements of whole building energy performance. This information is intended for practitioners who are concerned about transitioning between simulation tools with different engines and for managers of utility programs leveraging these tools for evaluating and/or projecting measure savings to be incentivized under their programs. The results of the comparisons of simulated and measured performance indicate that the predictions from all three engines are not significantly different. The 24-hour average value of the absolute mean bias indicates the likely magnitude of the error in any particular case. The average mean bias is reduced by cancelation of overprediction in one case by underprediction in another. The daytime absolute mean biases, which may be more important for both energy performance and occupant comfort, are ~6%, presumably because of the greater complexity involved in simulating in the presence of solar radiation. EnergyPlus typically overpredicts the cooling load and/or underpredicts the heating load by ~1.5% and the DOE-2 engines typically underpredict the cooling load by approximately the same amount. The Root Mean Square Error is relatively more sensitive to shorter term variations in the difference between predicted and measured loads; the three engines have similar values, ~10%, suggesting that the uncertainties in their predictions of peak loads may also be similar in magnitude. The implication of these results is that users, both designers and program analysts, can use EnergyPlus, DOE-2.1e, or DOE-2.2 to model conventional commercial buildings equipped with overhead mixing ventilation with a similar level of confidence. Further work is required to better understand the variability in the level of agreement between the engine predictions and FLEXLAB measurements, where a particular engine will agree well with FLEXLAB in some cases and not so well in others and another engine will agree or disagree in different cases. As the sources of this variability are identified and eliminated or reduced significantly, it is recommended that the experimental capabilities and methods developed in the study reported here should be applied to validating heating and cooling load calculations for spaces with different types of furniture and miscellaneous loads. These methods should then be applied to low energy space conditioning systems in EnergyPlus including, in particular, radiant slab and radiant ceiling panel cooling and heating systems and ‘mixed mode’ systems that combine mechanical cooling and natural ventilation systems, focusing on controls, including control of thermal mass. The work reported here addresses the conventional method of heating and cooling occupied spaces; other methods, such as the use of radiant heating and cooling systems have the potential to provide equivalent occupant comfort, or better, with lower energy consumption. These systems are addressed more explicitly in EnergyPlus but there is a need for empirical validation to give users the same level of confidence in modeling these systems that they have, or should have, in modeling conventional systems, based on the results presented here.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗