Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “load uncertainty”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Regulatory Testing of RPP-WTP HLW Glasses to Support Delisting Compliance, VSL-04R4780-1, Rev. 0 (Sep 2004)

The primary goal of the testing described in this report was to collect data to demonstrate compliance of the immobilized high-level waste (IHLW) glasses with delisting requirements. The collected data will be used to support a petition to delist the IHLW glasses destined for the national disposal facility. The Delisting Data Quality Objectives (DQO) (Cook and Blumenkranz 2003) identified a list of (16) inorganic constituents of potential concern (COPCs) and their associated limits for delisting. These COPCs can be divided into three groups, Cases 1, 2, and 3, based on their Toxicity Characteristic Leaching Procedure (TCLP) responses versus their respective delisting limits. To briefly summarize, Case 1 COPCs are those that, when loaded at their highest expected concentration in Waste Treatment Plant (WTP) glasses, are not expected to leach at their respective delisting limits when the glasses are exposed to the TCLP. Case 2 COPCs may reach the delisting limits in TCLP leachates of WTP glasses if loaded to concentrations near their maximum expected concentrations in glass. Finally, Case 3 COPCs are components that are likely to be present in concentrations sufficient to exceed their respective delisting limits in TCLP leachates of some possible glasses. The test objective was to show that, for (i) the expected range of inorganic contents in the waste feed to the Hanford WTP high level waste (HLW) vitrification facility, (ii) the expected range of glass product compositions, and (iii) the Case 1 and Case 2 COPCs identified by the DQO, the IHLW glasses meet all the relevant requirements for delisting. For Case 3 COPCs (i.e., Cd), the testing was to demonstrate the relationship between glass composition and TCLP cadmium (Cd) release, and then employ the results to develop TCLP-composition response models. During WTP operations, TCLP-composition models can be used to predict, within the required statistical uncertainties, TCLP responses of IHLW production glasses that are within the compositional region used to develop the model.

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

Statistical Correlation of Heliostat Pointing Deviation With Wind

This work was carried out as part of the Heliostat Consortium (HelioCon) Field Deployment subtask with the aim to develop a reduced order model framework for correlating wind speed and pointing deviation of a heliostat facet. There are only sparse field measurements of heliostat pointing deviations and accompanying wind conditions published in the literature. Heliostat test standards, such as IEC 62862-4-3, propose a suite of tests including laser pointing repeatability at wind speeds below 4 m/s, and provide technical requirements for heliostat slope and tracking deviations in coarse average wind speed bins of 4 m/s, 6 m/s, and 8 m/s. In addressing the gap of the variation of heliostat pointing deviation with wind speed, field measurements of laser pointing on a grid target and wind conditions were analyzed in this study at the Third-Party Metrology Platform at the National Laboratory of the Rockies (NLR) Flatirons Campus. Horizontal pointing deviations were found to follow a logarithmic relationship with peak wind speed, whereas vertical pointing deviations follow an exponential relationship with peak wind speed. Both horizontal and vertical pointing deviations also follow a second order polynomial relationship, as expected from the proportionality of elastic loads and deformations with the square of wind speed. The results indicate that heliostat facet pointing deviations in the vertical direction increase at a faster rate than in the horizontal direction with increasing wind speed over the tested range, however these are dependent on the heliostat structural design. Next steps are recommended for additional field measurements to confirm a linear relationship of pointing deviation with applied moment on a heliostat facet, and to distinguish between gravity-induced and wind-induced pointing deviations at different elevation angles. The derived correlations in the preliminary analysis in this report serve as a case study for heliostat developers and plant operators to estimate the wind-induced pointing deviations and their variation with peak gust wind speed. Next steps in future work would recommend higher resolution and longer duration datasets for different elevation angles and wind directions to reduce uncertainties and variance of collected laser beam spot data and their correlations with bin-averaged wind speed.

17 WIND ENERGY↗
DOE OSTI2023

Off-policy deep reinforcement learning with automatic entropy adjustment for adaptive online grid emergency control

Electric overloading conditions and contingencies put modern power systems at risk of voltage collapse and blackouts. Load shedding is crucial to maintain voltage stability for grid emergency control. However, the rule- or model-based schemes rely on accurate dynamic system models and face considerable challenges in adapting to various operating conditions and uncertain event occurrences. Here, to address these issues, this paper proposes a novel deep reinforcement learning (DRL)-based voltage stability control algorithm with automatic entropy adjustment (AEA) for grid emergency control. Various dynamic network components for complex system operations are modeled to construct the DRL environment. An off-policy soft actor-critic architecture is developed to maximize the expected reward and policy entropy simultaneously. The AEA mechanism is proposed to facilitate the policy maximum entropy procedure, and the proposed method can automatically provide effective discrete and continuous actions against various fault scenarios. Our approach accomplishes high sampling efficiency, scalability, and auto-adaptivity of the control policies under high uncertainties. Comparative studies with the existing DRL-based control methods in IEEE benchmarks indicate salient performance improvement of the proposed method for dynamic system emergency control.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Utilization of ACE nuclear data file toolkit ACEtk to calculate relative sensitivity coefficients of point-kinetics parameters

Sensitivity and uncertainty methods are quintessential for nuclear criticality safety and experiment design. This type of analysis relies on calculations of sensitivity coefficients; sensitivity coefficients of the effective neutron multiplication factor with respect to some nuclear data are predominantly calculated and used. As a part of the Laboratory Directed Research & Development project EUCLID (Experiments Underpinned by Computational Learning for Improvements in nuclear Data) at Los Alamos National Laboratory, sensitivity coefficients of many radiation detector measurement responses with respect to nuclear data were investigated. Specifically, this paper outlines a method to calculate point-kinetics parameters relative sensitivity coefficients with respect to nuclear data. Point-kinetics parameters such as the prompt neutron decay constant, effective delayed neutron fraction, and neutron generation time are especially important to experimenters and reactor operators designing systems with dynamic neutron populations. This method couples capabilities of the ACE (A Compact ENDF) nuclear data file toolkit, ACEtk, with the ability to load cross sections into the radiation transport code Monte Carlo N-Particle (MCNP). In conclusion, key aspects of optimizing this method for a particular application and sensitivity profiles of the Jezebel criticality experiment are examined and discussed.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2022

High-Fidelity Measurements for Flattop-HEU Benchmark Reevaluation

Flattop was first built in the 1950’s at Los Alamos National Laboratory. Flattop-HEU is composed of a sphere of highly enriched uranium (HEU) surrounded by a thick spherical natural uranium (NU) reflector. The reflector is composed of three parts: a stationary hemisphere and two movable quarter spheres. For fine control of the reactivity of the system, there are three control rods of natural uranium located in voids in the stationary hemisphere. The final components that make Flattop a useful critical assembly are the glory hole and mass adjustment pieces. These pieces can be loaded in various configurations into the glory hole and the core pedestal to control the known worth of the system. The glory hole and mass adjustment pieces are mostly small pieces of HEU with some mass adjustment pieces fabricated from NU. This allows for the irradiation of samples to a specified level. To better document the system, Flattop was evaluated and included in the International Criticality Safety Benchmark Evaluation Project (ICSBEP) handbook. The original benchmark evaluation of Flattop-HEU was written in 1999 based on an experiment completed in the 1960’s. This original evaluation was written to provide a single diameter that defined critical mass; however, as computational capabilities have increased, the focus for benchmark evaluations has shifted to include detailed modelswith all physical dimensions. Thus, as Flattop is a lynchpin in critical experiment work, the benchmark is being reevaluated at current standards. This summary discusses some of the largest known uncertainties from the evaluation and the high-fidelity measurements taken to reduce these uncertainties.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Voltage positioning using co-optimization of controllable grid assets in radial networks

With increasing penetration of solar PV, some distribution feeders are experiencing highly variable net-load flows and even reverse flows. To optimize distribution systems under such conditions, the scheduling of mechanical devices, such as OLTCs and capacitor banks, needs to take into account forecasted solar PV and actual grid conditions. However, these legacy switching assets are operated on a daily or hourly timescale, due to the wear and tear associated with mechanical switching, which makes them unsuitable for real-time control. Therefore, there is a natural timescale- separation between these slower mechanical assets and the responsive nature of inverter-based resources. In this paper, we present a network admissible convex formulation for holistically scheduling controllable grid assets to position voltage optimally against solar PV. An optimal hourly schedule is presented that utilizes mechanical resources to position the predicted voltages close to nominal values, while minimizing the use of inverter-based resources (i.e., DERs), making them available for control at a faster time-scale (after the uncertainty reveals itself). A convex, inner approximation of the OPF problem is adapted to a mixed-integer linear program that minimizes voltage deviations from nominal (i.e., maximizes voltage margins). Here, the resulting OPF solution respects all the network constraints and is, hence, robust against modeling simplifications. Simulation based analysis on IEEE distribution feeders validates the approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Probing the atmospheric boundary layer with integrated remote-sensing platforms during the American WAKE ExperimeNt (AWAKEN) campaign

The American WAKE ExperimeNt (AWAKEN) collaboration is an observational-based field campaign in northern Oklahoma intended to analyze the potential influence of onshore wind farms and their collective wakes on wind power production, turbine structural loads, and on the atmospheric boundary layer (ABL). Focusing on the ABL effects, the University of Oklahoma and the Lawrence Livermore National Laboratory collected continuous high-resolution kinematic and thermodynamic profile measurements during 2022 and Summer 2023. The deployment strategy for these campaigns is detailed first, followed by an initial comparison of data from two sites in the AWAKEN domain: a near-farm site to examine collective wake impacts on the ABL, and a far-field site remaining outside the wind farm-waked region. Here, we summarize the datasets available and demonstrate the benefits of these observations and multiple value-added products (VAPs) for investigation of ABL features observed during AWAKEN. We also highlight examples of preliminary analyses, including ABL height detection and nocturnal low-level jet examination, which are produced using novel VAPs based on optimal estimation to retrieve deeper Doppler lidar wind profiles than previously resolved, along with their uncertainty. By including the near-farm and far-field site in these analyses, we identified a pattern of stronger lower-atmospheric mixing at the near-farm site than the far-field site, motivating deeper investigation into the relationship between wind farms and general ABL characteristics. Future analysis will delve deeper into this relationship by examining other ABL characteristics, such as atmospheric stability and convection.

17 WIND ENERGY↗
DOE OSTI2025

Sensitivity-based voltage constraints for optimal power flow in low-voltage distribution feeders

The optimal power flow (OPF) problem for distribution systems can include network details down to the low-voltage (LV) points of interconnection of individual customers. This paper addresses the implementation of voltage magnitude constraints, and sets forth a practicable approach for capturing the effects on voltage from the switching behavior of loads (e.g., heat pumps, air conditioners, water heaters, or pool pumps) and from the variability of renewable generation (e.g., rooftop solar). The proposed method adjusts the OPF voltage constraints based on forecasts of load and generation upper and lower bounds, in conjunction with sensitivity factors derived from the power flow equations. An illustrative OPF formulation is also provided, which incorporates transformer models that include core loss. We demonstrate that accurate modeling of these LV network components is critical to avoid voltage violations at customer points of interconnection. Furthermore, the ideas are validated through numerical case studies on a realistic distribution feeder.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

MECHANISTIC MULTIPHYSICS MODELING OF CLADDING RUPTURE IN NUCLEAR FUEL RODS DURING LOSS-OF-COOLANT ACCIDENT CONDITIONS

The Loss of Coolant Accident (LOCA) is a design basis accident that is included as part of the safety analysis of nuclear power plants. As the nuclear industry desires to increase the cycle length and discharge burnup of existing nuclear power plants they must demonstrate safe operation during a LOCA on high burnup fuel. During a LOCA transient on high burnup fuel rods, the rods may undergo a process known as fuel fragmentation, relocation, and dispersal (FFRD). To permit dispersal, the cladding encapsulating the fuel must rupture with an opening size large enough to allow the fragmented fuel particles to release. Current licensing tools used by industry and the United States Nuclear Regulatory Commission are limited in geometric fidelity and materials that can be analyzed. These simulation tools generally employ a quasi-two-dimensional (1.5D or Layered1D) or 2D-RZ axisymmetric geometric representations exclusively. While a valid approach under some instances, there are times when important physics have an asymmetric behavior in the fuel rod. Examples include fuel fragmentation, thermal-hydraulic boundary conditions, and cladding rupture, all of which are important for LOCA analysis. As industry pursues burnup extensions it must be demonstrated that fuel dispersal can be mitigated or eliminated. To do this, an understanding of the rupture opening after cladding failure is required. This work presents the development of a model for predicting the size and location of the rupture opening in failed fuel rods during LOCA conditions using advanced modeling and simulation tools. In order to supply the rupture model with appropriate boundary conditions, improvements to fuel fragmentation, axial relocation and oxidation modeling were required. First, the eXtended Finite Element Method (XFEM) is used to mechanistically predict the number of fuel fragments that form due to material strength randomization, criteria for strength randomization, mesh density, power ramping rates and irradiation effects. These predictions with associated uncertainty were compared to empirical correlations developed for UO2 verifying that they can be used with increased confidence in subsequent axial relocation analyses. Secondly, a new first-of-its-kind Layered2D computational framework was developed that provides the ability to apply azimuthally varying boundary conditions while providing discrete layers to track fuel movement during the LOCA. An existing fuel axial relocation model developed for Layered1D was extended to work within the Layered2D framework. A large sensitivity study was performed on the initial version of the model to identify modeling parameters of particular importance, with the emissivity used for radiation after blowdown being the primary source of uncertainty. Then, a simplistic approach to incorporate mechanical degradation of the cladding due to oxidation was also developed to investigate the impact of reduced cladding thickness on predictions of the time to failure of cladding tubes. It was found that the cladding will typically rupture prior to a reduction in thickness that is sufficient to impact the rupture behavior. A model was then developed for predicting cladding rupture that transfers the cladding surface temperatures, rod internal and external pressures, fast neutron fluence, and fast neutron flux from a more detailed Layered1D, Layered2D, or 2D-RZ analysis to a 3D cladding only analysis. Comparisons of the rupture model to a few experiments indicate reasonable predictions. The rupture model was then applied to two accident tolerant fuel concepts (FeCrAl and Cr-coated Zircaloy) where it predicted that both ATF concepts would have smaller rupture openings and delayed rupture times than the standard Zircaloy-4 cladding material under identical loading conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Three-dimensional microstructure-explicit and void-explicit mesoscale simulations of detonation of HMX at millimeter sample size scale

Fully three-dimensional (3D) microstructure-explicit and void-explicit mesoscale simulations of the shock-to-detonation (SDT) process of pressed granular HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine) are performed. The overall size scale of the models is up to 3 × 3 × 15 mm3, with ∼30 000 grains and 206 265 voids. The models account for the heterogeneous material microstructure, constituent distribution, constituent morphology, and voids. Loading conditions considered involve piston velocities in the range of 600–1200 m/s or pressures in the range of 4–8 GPa. The focus is on analyzing the SDT process and the effects of microstructure and voids on the run-to-detonation distance (RDD). Companion two-dimensional (2D) simulations are also carried out to assess the differences between 2D and 3D. Statistically equivalent microstructure sample sets (SEMSSs) are generated and used for both 2D and 3D, allowing the prediction of the statistical and probabilistic Pop plots (PPs). The predictions are in general agreement with trends in available experimental data in the literature. It is found that both the microstructure (heterogeneous grain size, morphology, and size distribution) and voids significantly affect the RDD and the PPs. These effects are systematically delineated and quantified via the use of SEMSSs with different combinations of attributes. A recently developed probabilistic formulation for the PPs is used to characterize the results, allowing uncertainties in the relations between the shock pressure and RDD arising from material heterogeneities to be quantified. The probabilistic formulation is further used to quantify the confidence levels in the ranked order of influences of different combinations of microstructure and voids on the PPs.

Miller, Christopher (ORCID:0000000206844381)↗
DOE OSTI2023

Evaluation of thermostat location for multizone commercial building performance

In multi-zone buildings, it is often found that a single shared thermostat controls more than one conditioned zones. Although these shared zones are supposed to have similar thermal needs (e.g., cooling and heating load), in reality, they are not mainly due to different orientations, sizes of windows, occupancy, space types, etc. This can cause unnecessary energy waste or thermal discomfort for the occupants. How to quantify this impact in multizone buildings remains a research gap. Therefore, this study aims to evaluate the impact of different sensor (i.e., thermostat) locations for multizone commercial buildings through a comprehensive modeling study. Here, two different scenarios for the sensor locations were selected to evaluate the impact in terms of energy and thermal comfort. The scenario (1) is that one to five sensors distributed among the five zones, but the sensor readings from selected zones will be used for no-sensor zones, which is no-mean sensor scenario. The scenario (2) is that one to five sensors distributed among the five zones, but the average temperature from the shared zones will be used for each of the shared zones, which is a mean sensor scenario. The uncertainty analysis was performed for different sensor location scenarios.(a)The major findings from an energy perspective, for scenario (1), the differences of cooling energy go as high as 17% more or 12% less, compared with the baseline. For heating energy consumption, the discrepancies go as high as 51% more or 52% less, compared with baseline. For site energy consumption, the discrepancies go as high as 3.2% more or 3.2% less, compared with baseline. For fan energy consumption, the discrepancies go as high as 3.2% more or as low as 1.0% less, compared with baseline. For scenario (2), the discrepancies of cooling energy go as high as 3% more, or 0.5% less, compared with the baseline. For heating energy consumption, the discrepancies, are go as high as 10.1% more or as low as 3.0% less, compared with baseline. For site energy consumption, the discrepancies go as high as 1.3% more or 0.3% less, compared with baseline. For fan energy consumption, the discrepancies go as high as 3.1% more or as low as 1.0% less, compared with baseline.(b) In terms of the indoor thermal comfort, for the no-mean-sensor scenarios, the discrepancies of unmet hours for cooling mode can be as high as 1,200 h, compared with the baseline. The discrepancies of unmet hours for heating mode can be as high as 740 h, compared with the baseline. For the mean-sensor scenarios, the discrepancies of unmet hours for cooling mode can be as high as 750 h, compared with the baseline. The discrepancies of unmet hours for heating mode can be as high as 50 h, compared with the baseline.

42 ENGINEERING↗
DOE OSTI2025

A scoping study of far-SOL main-wall protection limiters for steady-state operation of compact pilot plant tokamaks

We present a novel method for handling steady-state heat fluxes incident on the main wall of pilot plant-scale magnetic fusion devices, based on the utilization of protection limiters in the far scrape-off layer (SOL). This method helps avoid large plasma-wall gaps, without excessively compromising blanket performance. We present an optimization algorithm for determining the appropriate size and scale of these protection limiters given (1) probability distributions of SOL plasma parameters and (2) assumed risk tolerance. As part of this optimization, we have developed an analytic description of parallel heat fluxes across limiter shadows, and an objective cost function (the ‘Far-SOL Marginal Cost’) to quantify the impact that different main-wall thermal management design choices have on reactor capital cost. Applying the model to a midscale fusion pilot plant concept shows that making use of far-SOL protection limiters can reduce capital costs on the order of $500 M, relative to naively increasing the plasma-wall gap. Our analysis demonstrates that the far-SOL power decay length is the highest-leverage plasma assumption for thermal loading of the first wall, and the primary cost driver for main wall thermal management. The relative cost efficiency of protection limiters increases as assumptions on the far-SOL heat flux become more pessimistic. The concepts described in this paper motivate the further development of far-SOL protection limiters as part of larger efforts to design economical core-edge-wall compatible solutions for a fusion pilot plant.

Design under uncertainty↗
DOE OSTI2022

Calibration of cloud and aerosol related parameters for solar irradiance forecasts in WRF-solar

Model parameters are a major source of uncertainty in numerical weather prediction. Recently, the Weather Research and Forecasting model with Solar extensions (WRF-Solar) has been upgraded by enhancing the treatment of sub-grid scale cloud and aerosols with augmentations of a sub-grid scale cloud scheme (CLD3) and an upgraded aerosol-aware Thompson-Eidhammer scheme (TE14). However, the value of model parameters associated with these parameterizations are assigned based on limited measurements or theoretical calculations. Calibrating the most sensitive parameters has the potential to improve solar irradiance predictions. Here, we adopted a multiobjective surrogate-based optimization (SBO) framework to calibrate nine parameters used in CLD3 and TE14 that lead to the largest sensitivity in simulated irradiance. The normalized mean-absolute-error (NMAE) of global horizontal irradiance (GHI) and direct normal irradiance (DNI) are minimized by calibrating WRF-Solar over two regions including the Southern Great Plains (SGP) and Central California, in order to focus on parameter calibration under cloudy conditions with different aerosol loading. The results show that generalized linear model (GLM)-based surrogate models approximate physical models well, particularly when the third order and three-way interaction terms are considered. The SBO framework efficiently searches the parameter space for optimal solutions with less computational costs than directly calibrating the physical model. We first calibrate CLD3 parameters over the less-polluted SGP region. Optimized CLD3 parameters alone result in NMAE reduction by 14% for the site-mean and up to 33% for individual cases over the SGP region. With further calibration of TE14 parameters over the Central California during active fire periods, the optimized parameters lead to over 20% reductions of NMAE. Our investigation reveals, however, that optimizing TE14 has a limited impact on irradiance simulations under less-polluted conditions in the SGP.

14 SOLAR ENERGY↗
DOE OSTI2020

Making coal relevant for small scale applications: Modular gasification for syngas/engine CHP applications in challenging environments

Small-scale coal gasification technology, coupled to a reciprocating engine generator, has the potential for making coal a cost competitive resource for meeting the flexible energy needs and resiliency requirements of utilities across the United States. To maintain grid stability and reliability, electrical generation must be regulated to match the load at the proper voltage and frequency. With the expansion of intermittent sources into the grid, such as wind and solar, frequency and voltage regulation become increasingly important and challenging. This work presents the results of a Front End Engineering Design (FEED) effort to detail the engineering and preliminary economics of a small-scale, air blown, fixed-bed gasification process, operating at near-atmospheric pressure, with gas cleanup to provide syngas and pyrolysis liquid fuels for use in reciprocating engine generators for combined heat and power at the University of Alaska Fairbanks. A very detailed assessment of capital and operating costs allows the evaluation of a levelized cost of electricity of 208.06 $/MWh, which can be reduced through the selling of by-products (steam for heating purposes and pyrolysis liquids to be used for power generation in an existing diesel engine). Here, a combined sensitivity analysis, based on the Monte Carlo approach, has been carried out to evaluate the effects of the uncertainties (capital and operating costs and plant annual availability) on the LCOE.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

THERMAL MODELING OF HANFORD CESIUM AND STRONTIUM CANISTERS DURING SIMULATED LOADING

A computational fluid dynamics (CFD) model was built to simulate planned testing of heater assemblies within a canister and overpack for the Hanford Lead Canister (HLC) project. The HLC is a canister storage system that will contain heaters to simulate the decay heat of nuclear material and provide the canister storage system with environmental conditions equivalent to the operating conditions on a dry storage pad. The HLC will be equipped with long-term data collection and monitoring systems to provide an early warning of corrosion, pitting, cracking, or other signs of canister degradation that might threaten the integrity of the containment boundary over the potentially long term of dry storage. An important part of the HLC development is to make pretest numerical predictions for the behavior of the heated canister during the simulated radiolytic decay heat testing, which simulates the dry storage system during loading operations. The simulated radiolytic decay heat test is planned for mid-2024 in a configuration that includes the heater assembly, overpack, and canister, but with the lids removed to allow loading cesium and strontium capsules into the canister. One of the goals of the test is to evaluate the thermal behavior of the canister and overpack assembly in the ambient air of the test facility, which will provide data critical to validating the thermal models and understanding how the HLC will perform as a system once deployed. To best approximate real-world conditions, the CFD model includes the full air volume of the mock-up truck bay the heated canister test will be performed in, enabling detailed investigation of how the heated canister affects airflow around it. Rigorous pre-deployment testing of the complete HLC cask and canister system is intended to be completed before the HLC is deployed in the 2028 timeframe. This study presents the pre-test temperature predictions of the simulated radiolytic decay heat test. A description of the heater assembly, canister, and overpack system is presented. The model was developed with the commercial CFD software STAR-CCM+. An uncertainty analysis was run with the CFD model to determine the uncertainty in the temperature predictions and provide a range over which the predicted temperatures are expected to vary. The uncertainty analysis was preformed by coupling STAR-CCM+ with the software Dakota, which provides advanced parametric analyses, including quantification of margins and uncertainty with computational models. This work is expected to provide insight into SNF canister behavior.

Carpenter-Graffy, Dina E.↗
DOE OSTI2022

Resolving glacial isostatic adjustment (GIA) in response to modern and future ice loss at marine grounding lines in West Antarctica

Accurate glacial isostatic adjustment (GIA) modelling in the cryosphere is required for interpreting satellite, geophysical and geological records and for assessing the feedbacks of Earth deformation and sea-level change on marine ice-sheet grounding lines. GIA modelling in areas of active ice loss in West Antarctica is particularly challenging because the ice is underlain by laterally varying mantle viscosities that are up to several orders of magnitude lower than the global average, leading to a faster and more localised response of the solid Earth to ongoing and future ice-sheet retreat and necessitating GIA models that incorporate 3-D viscoelastic Earth structure. Improvements to GIA models allow for computation of the viscoelastic response of the Earth to surface ice loading at sub-kilometre resolution, and ice-sheet models and observational products now provide the inputs to GIA models at comparably unprecedented detail. However, the resolution required to accurately capture GIA in models remains poorly understood, and high-resolution calculations come at heavy computational expense. We adopt a 3-D GIA model with a range of Earth structure models based on recent seismic tomography and geodetic data to perform a comprehensive analysis of the influence of grid resolution on predictions of GIA in the Amundsen Sea Embayment (ASE) in West Antarctica. Through idealised sensitivity testing down to sub-kilometre resolution with spatially isolated ice loading changes, we find that a grid resolution of ~$\frac{1}{3}$ of the radius of the load or higher is required to accurately capture the elastic response of the Earth. However, when we consider more realistic, spatially coherent ice loss scenarios based on modern observational records and future ice-sheet model projections and adopt a viscoelastic Earth, we find that predicted deformation and sea-level change along the grounding line converge to within 5 % with grid resolutions of 7.5 km or higher, and to within 2% for grid resolutions of 3.75 km and higher, even when the input ice model is on a 1 km grid. Furthermore, we show that low mantle viscosities beneath the ASE lead to viscous deformation that contributes to the instrumental record on decadal timescales and equals or dominates over elastic effects by the end of the 21st century. Our findings suggest that for the range of resolutions of 1.9–15 km that we considered, the error due to adopting a coarser grid in this region is negligible compared to the effect of neglecting viscous effects and the uncertainty in the adopted mantle viscosity structure.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Bi-Level Adaptive Storage Expansion Strategy for Microgrids Using Deep Reinforcement Learning

Battery energy storage (BES) is a versatile resource for the secure and economic operation of microgrids (MGs). Prevailing stochastic optimization-based approaches for BES expansion planning for MGs are computationally complicated. This work proposes a data-driven bi-level multi-period BES expansion planning framework to determine the siting, sizing, and timing of BES installations. The proposed planning framework unifies deep reinforcement learning (DRL) and linear programming, thereby decoupling the determinations for the integer and continuous decision variables in two time scales, respectively. In the upper level, a rainbow DRL agent with quantile regression is trained to provide dynamic planning policies to accommodate stochastic renewable energy resources (RESs), load, and battery price changes efficiently. Further, the lower level computes the optimal operation of MGs with frequency constraints to hedge the islanding contingency. The two levels communicate with one another by exchanging storage configuration and operating expenses in order to accomplish the shared goal of minimizing investment and operation costs. Comparative case studies on an MG are carried out to demonstrate the superiority of the proposed DRL-based solution to the mixed-integer linear programming counterpart on efficiency, scalability, and adaptability.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Timber Harvesting Economics in Early Commercial Thinning and Group Selection Prescriptions

Forest management decisions that consider both short- and long-term objectives are becoming more challenging and costly, due to market uncertainty, labor shortages, and the evolving forest industry. Under such unpredictable circumstances, a better understanding of the production economics of various silvicultural prescriptions is crucial in implementing sound forest management activities. The objectives of this study were to: a) estimate stump to landing productivity and cost for Early Commercial Thinning (ECT) and Group Selection (GS) prescriptions; and b) compare the harvest system cost ($/PMH) using machine rate and PATH (Planning and Analysis in Timber Harvesting) methods. Detailed time-motion studies were conducted in the spruce-fir forest of Northern Maine, USA, during the summer of 2019. The results, as expected, indicated the harvesting costs were 23% higher for ECT (21.08 $/m3) compared to GS (16.32 $/m3). The average productivity of the harvester and forwarder were reduced by 15 and 37%, respectively, for ECT compared to GS. The estimated harvesting system cost ($/PMH) was 4% higher as indicated by PATH compared to machine rate method. Higher operational productivity in GS can be attributed to the larger volume of logs, and shorter travel loaded distance than ECT. These findings can help forest managers make more informed decisions while implementing appropriate silvicultural prescriptions under varying conditions based on management objectives and market conditions.

Forestry↗