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

Delivery-Risk-Aware Flexibility Scheduling and Dispatch for Aggregated Flexible Loads

Flexible loads like smart thermostats and water heaters can shift energy consumption and provide flexibility to the grid. However, this flexibility is dependent on occupant behavior and can lead to delivery risk, which causes utilities and grid operations to consider them as unreliable for purposes of grid operation. To date, they have not been well integrated into wholesale electricity markets or ancillary service offerings. With proper consideration of uncertainty and risk, these resources can be one of the most cost-effective sources of flexibility. This work uses stochastic optimization to quantify and bid flexibility from a fleet of flexible resources while considering their delivery risk.

DER↗
DOE OSTI2021

Development of Digital Twin Predictive Model for PWR Components: Updates on Multi Times Series Temperature Prediction Using Recurrent Neural Network, DMW Fatigue Tests, System Level Thermal-Mechanical-Stress Analysis

The long-term operation (LTO) of nuclear power plant (NPP) beyond their original design life of 40 years, can lead to more material damage associated with cyclic fatigue under thermal-mechanical loading cycles and associated long-term exposure of reactor material to the deleterious reactor-coolant environments. However, under this LTO condition the reactor components can still safely operate but may require more frequent Nondestructive Evaluation (NDE) of reactor components. Frequent NDE requirement may lead to frequent shutdown of the NPP. This in turn can lead to power outage and additional NDE-inspection-cost related economic loss. The economic loss can be minimized by reducing uncertainty in life estimation of safety-critical pressure boundary components and by implementing more digital approach such as by using upcoming digital-twin (DT) technology for predicting the structural states (e.g., time and location dependent inside/outside thickness temperature, stress, strain, plastic deformation, etc.) and associated fatigue life of a component in real time. Towards this goal Argonne National Laboratory (ANL) with the sponsorship of DOE Light Water Reactor Sustainability (LWRS) program is working on the development of a DT framework that can be used for real time environmental fatigue prediction of reactor components. The DT framework is based on limited experiment-data, Artificial-intelligence (AI) – Machine-Learning (ML) - Deep-Learning (DL) based techniques and Multiphysics-computational-mechanics such as finite element (FE) based modeling tools. Towards this overall goal, following are some of the major contributions made during the FY21: 1) Multiple 82/182 dissimilar metal weld (DMW) specimens (both solid-weld and joint-weld representing the actual reactor multi-metal nozzles) were fatigue tested. The resulting fatigue lives were compared to the NUREG-6909 based best-fit and design fatigue curves. Additionally, the results of 52/152 DMW fatigue specimens (which were recently tested at Republic of Korea under the sponsorship of International Nuclear Energy Research Initiative - INERI program) were compared to the NUREG-6909 based best-fit and design fatigue curves. From the comparison of 82/182 and 52/152 DMW test data with NUREG-6909 best-fit curve, most of the reported test data fall way away from the NUREG-6909 suggested best-fit or mean curve. The NUREG-6909 suggested best-fit curve is the best-fit curve of austenitic stainless steel and due to lack of enough data on Nickel-based welds, this is currently being used for predicting the life of Nickel-alloy-based welded components. However, the above observation may require higher scaling factor (e.g., ASME suggested factor of 20 on cycles rather than the current NUREG-6909 suggested factor of 12 on cycles) for scaling the austenitic-stainless-steel best-fit-curve for estimating the design or safe-life of a welded component. Accordingly, for example, if a DMW component experience a strain amplitude of 0.6% the PWR-water life of the component would be 52 cycles instead of 85 cycles. However, more DMW tests are required to further ascertain the above-mentioned observations. 2) A system level CAD and finite element model were developed which consists of reactor pressure vessel (RPV), part of steam generator (SG), part of pressurizer (PRZ), hot leg (HL), and surge line (SL). This is with detailed nozzle geometry and thermal-mechanical material properties of different metals to simulate realistic thermal-mechanical stress under connected system global thermal-mechanical boundary conditions. 3) Different system level heat transfer analyses were performed with estimation of relevant heat transfer coefficients. The resulting data were used in subsequent system level thermal-mechanical stress analysis and for generating spatial-temporal training and validation data for a system level digital-twin based temperature predictor. Transient heat transfer analyses were performed considering thermal boundary condition under design-basis (DB) loading and EDF (Électricité de France) data-based grid-load-following (EDF-GLF) loading cycles. 4) System level thermal-mechanical stress analysis was performed for identifying damage-prone hotspots and for future extension of the model for cyclic state prediction. From the system-level model simulation under DB loading cycle it is found that HL and the SL nozzle that connect to the HL can experience significant stress and strain and could be one of the weakest links in the overall reactor coolant system (RCS). 5) An AI/ML based DT model was developed for multi-time-series temperature prediction at any inside/outside thickness locations of PWR pressure boundary components. This is by using Recurrent-neural-network (RNN) and keras machine learning libraries. The RNN model was validated against two laboratory test-based data sets with one obtained through ANL’s in-air fatigue test system and other through PWR-water test loop. The experimentally validated DT model further validated against FE model results to predict thermal scarification related spatialtemporal temperatures at random locations of a component. The well validated DT model was then used for demonstrating spatial-temporal temperature prediction under 100+ years of reactor operation subjected to combined DB, EDF-GLF and randomized grid-load-following (RANDOMGLF) loading Cycles. The expert-elicitation DT model framework was developed assuming field/input/process measurements can be available from a few existing plant sensors and can readily be used by the NPP operators. The above temperature prediction model will feed to the next-step stress analysis model based on which the life of a component can be predicted in realtime, which is one of our future works.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Evaluation of the performance of event reconstruction algorithms in the JSNS 2 experiment using a 252 Cf calibration source

JSNS 2 investigates short-baseline neutrino oscillations using a 24-meter baseline and a 17-tonne Gd-loaded liquid scintillator target. Accurate event-reconstruction algorithms are crucial for analyzing experimental data. The algorithms undergo meticulous validation through calibration with a 252 Cf source. This paper outlines the methodology and evaluates the reconstruction performance, focusing on neutrino interactions up to approximately 50 MeV for sterile neutrino searches. Both 252 Cf and Michel electron events are studied to evaluate reconstruction accuracy. The analysis concludes that the uncertainty of the fiducial volume, with an appropriate correction, is much less than the requirement of JSNS 2 requirement (10%). Furthermore, the energy resolution is measured to be 3.3 ± 0.1% for the Michel electron endpoint and 4.3 ± 0.1% for the n-Gd peak in the central region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

An Investigation of Strategies for Sensing Window Attachment Position: Survey of Candidate Technologies, Evaluation, and Deployment Recommendations

While windows are desirable feature in buildings for views, daylighting, and aesthetics, it is important to recognize that windows are a major contributor to a building’s energy load. Interior window coverings, such as shades and blinds, can help improve the energy performance of windows without sacrificing the benefits; however, while appropriate use of window coverings provide the potential for energy savings, there is little consensus regarding how occupants typically use these window attachments, which sometimes puts into question how much energy they actually save in an occupied home. Automated shading control strategies can reduce the uncertainty of energy savings; however, users can overwrite pre-defined schedules or write their own shade control programs to reach their personal comfort and privacy goals. Because it is difficult to identify how occupants interact with window shades, it is difficult to identify the energy savings potential of blinds. To address this problem, this study was developed to identify a candidate sensor that can be deployed to sense blind position. The intended application for this sensor is deployment in energy efficiency field studies so that shade usage can be measured to better understand results and identify the energy impact of specific window attachment technologies.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

BeyondPlanck: V. Minimal ADC Corrections for Planck LFI

We describe the correction procedure for Analog-to-Digital Converter (ADC) differential non-linearities (DNL) adopted in the Bayesian end-to-end BEYONDPLANCK analysis framework. This method is nearly identical to that developed for the official Planck Low Frequency Instrument (LFI) Data Processing Center (DPC) analysis, and relies on the binned rms noise profile of each detector data stream. However, rather than building the correction profile directly from the raw rms profile, we first fit a Gaussian to each significant ADC-induced rms decrement, and then derive the corresponding correction model from this smooth model. The main advantage of this approach is that only samples which are significantly affected by ADC DNLs are corrected, as opposed to the DPC approach in which the correction is applied to all samples, filtering out signals not associated with ADC DNLs. The new corrections are only applied to data for which there is a clear detection of the non-linearities, and for which they perform at least comparably with the DPC corrections. Out of a total of 88 LFI data streams (sky and reference load for each of the 44 detectors) we apply the new minimal ADC corrections in 25 cases, and maintain the DPC corrections in 8 cases. All these corrections are applied to 44 or 70 GHz channels, while, as in previous analyses, none of the 30 GHz ADCs show significant evidence of non-linearity. By comparing the BEYONDPLANCK and DPC ADC correction methods, we estimate that the residual ADC uncertainty is about two orders of magnitude below the total noise of both the 44 and 70 GHz channels, and their impact on current cosmological parameter estimation is small. However, we also show that non-idealities in the ADC corrections can generate sharp stripes in the final frequency maps, and these could be important for future joint analyses with the Planck High Frequency Instrument (HFI), Wilkinson Microwave Anisotropy Probe (WMAP), or other datasets. We therefore conclude that, although the existing corrections are adequate for LFI-based cosmological parameter analysis, further work on LFI ADC corrections is still warranted.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Long-Range Transport of Biomass Burning Aerosols from Southern Africa: A Case Study Using Layered Atlantic Smoke Interactions with Clouds Observations

A case study of an incoming biomass burning aerosol plume at Ascension Island is analyzed for the peak of the 2017 fire season using satellites, reanalysis and in situ observations. Measurements from the Atmospheric Radiation Measurement Mobile Facility 1 reveal an abrupt change from relatively clean conditions (~70 parts per billion by volume of carbon monoxide) to a more polluted state (~150 parts per billion by volume of carbon monoxide). Corresponding changes in aerosol size reveal a broadening of size distributions toward larger optical diameters, consistent with the arrival of aged aerosols. Within a 24 h period, black carbon fraction increases ~500% from ~300 ng m e to ~1500 ng m 3 , while light absorption coefficients increase ~300%. Long-range transport of these aerosols is primarily confined between 2 and 5 km above sea level along the northwesterly trade winds. Our results show that the primary driver of increases in aerosol loading over Ascension Island is an intensification of the St. Helena high-pressure system (anticyclone) that leads to a weakening of trade winds and increases westward transport on its northern flank. A better understanding of the complex interactions between air quality, meteorology and long-range aerosol transport is important for future modeling studies focused on aerosol–cloud–radiation interactions over the open ocean and reducing its associated uncertainties.

aerosol size distribution↗
NASA NTRS2005

Aerocapture Performance Analysis of A Venus Exploration Mission

A performance analysis of a Discovery Class Venus Exploration Mission in which aerocapture is used to capture a spacecraft into a 300km polar orbit for a two year science mission has been conducted to quantify its performance. A preliminary performance assessment determined that a high heritage 70 sphere-cone rigid aeroshell with a 0.25 lift to drag ratio has adequate control authority to provide an entry flight path angle corridor large enough for the mission s aerocapture maneuver. A 114 kilograms per square meter ballistic coefficient reference vehicle was developed from the science requirements and the preliminary assessment s heating indicators and deceleration loads. Performance analyses were conducted for the reference vehicle and for sensitivity studies on vehicle ballistic coefficient and maximum bank rate. The performance analyses used a high fidelity flight simulation within a Monte Carlo executive to define the aerocapture heating environment and deceleration loads and to determine mission success statistics. The simulation utilized the Program to Optimize Simulated Trajectories (POST) that was modified to include Venus specific atmospheric and planet models, aerodynamic characteristics, and interplanetary trajectory models. In addition to Venus specific models, an autonomous guidance system, HYPAS, and a pseudo flight controller were incorporated in the simulation. The Monte Carlo analyses incorporated a reference set of approach trajectory delivery errors, aerodynamic uncertainties, and atmospheric density variations. The reference performance analysis determined the reference vehicle achieves 100% successful capture and has a 99.87% probability of attaining the science orbit with a 90 meters per second delta V budget for post aerocapture orbital adjustments. A ballistic coefficient trade study conducted with reference uncertainties determined that the 0.25 L/D vehicle can achieve 100% successful capture with a ballistic coefficient of 228 kilograms per square meter and that the increased ballistic coefficient increases post aerocapture V budget to 134 meters per second for a 99.87% probability of attaining the science orbit. A trade study on vehicle bank rate determined that the 0.25 L/D vehicle can achieve 100% successful capture when the maximum bank rate is decreased from 30 deg/s to 20 deg/s. The decreased bank rate increases post aerocapture delta V budget to 102 meters per second for a 99.87% probability of attaining the science orbit.

Starr, Brett R.↗
NASA NTRS2010

A Mesoscale Analysis of Column-Integrated Aerosol Properties in Northern India During the TIGERZ 2008 Pre-Monsoon Period and a Comparison to MODIS Retrievals

The Indo-Gangetic Plain (IGP) of the northern Indian subcontinent produces anthropogenic pollution from urban, industrial and rural combustion sources nearly continuously and is affected by convection-induced winds driving desert and alluvial dust into the atmosphere during the premonsoon period. Within the IGP, the NASA Aerosol Robotic Network (AERONET) project initiated the TIGERZ measurement campaign in May 2008 with an intensive operational period from May 1 to June 23, 2008. Mesoscale spatial variability of aerosol optical depth (AOD, tau) measurements at 500mn was assessed at sites around Kanpur, India, with averages ranging from 0.31 to 0.89 for spatial variability study (SVS) deployments. Sites located downwind from the city of Kanpur indicated slightly higher average aerosol optical depth (delta Tau(sub 500)=0.03-0.09). In addition, SVS AOD area-averages were compared to the long-tenn Kanpur AERONET site data: Four SVS area-averages were within +/- 1 cr of the climatological mean of the Kanpur site, while one SVS was within 2sigma below climatology. For a SVS case using AERONET inversions, the 440-870mn Angstrom exponent of approximately 0.38, the 440-870mn absorption Angstrom exponent (AAE) of 1.15-1.53, and the sphericity parameter near zero suggested the occurrence of large, strongly absorbing, non-spherical aerosols over Kanpur (e.g., mixed black carbon and dust) as well as stronger absorption downwind of Kanpur. Furthermore, the 3km and lOkm Terra and Aqua MODIS C005 aerosol retrieval algorithms at tau(sub 550) were compared to the TIGERZ data set. Although MODIS retrievals at higher quality levels were comparable to the MODIS retrieval uncertainty, the total number of MODIS matchups (N) were reduced with subsequent quality levels (N=25, QA>=0; N=9,QA>=l; N=6, QA>=2; N=1, QA=3) over Kanpur during the premonsoon primarily due to the semi-bright surface, complex aerosol mixture and cloud-contaminated pixels. The TIGERZ 2008 data set provided a unique opportunity to measure the spatial and temporal variations of aerosol loading in the IGP. The strong aerosol absorption derived from ground-based sun/sky radiometer measurements suggested the presence of a predominately black carbon and dust mixture during the pre-monsoon period. Consistent with the elevated heat-pump hypothesis, these absorbing aerosols found across Kanpur and the greater IGP region during the pre-monsoon period likely induced regional atmospheric warming, which lead to a more rapid advance of the southwest Asian monsoon and above normal precipitation over northern India in June 2008.

Giles, D. M.↗
DOE OSTI2022

OPF-Learn: An Open-Source Framework for Creating Representative AC Optimal Power Flow Datasets: Preprint

Increasing levels of renewable generation motivate a growing interest in data-driven approaches for AC optimal power flow (AC OPF) to manage uncertainty. However, a lack of disciplined dataset creation and benchmarking prohibits useful comparison between approaches in the literature. To instigate confidence, models must be able to reliably predict solutions across a wide range of operating conditions. This paper develops the OPF-Learn package for Julia and Python which uses a computationally efficient approach to create representative datasets that span a wide spectrum of the AC OPF feasible region. Load profiles are uniformly sampled from a convex set that contains the AC OPF feasible set. For each infeasible point found, the convex set is reduced using infeasibility certificates, found by utilizing properties of a relaxed formulation. The framework is shown to generate datasets which are more representative of the entire feasible space versus traditional techniques seen in the literature, improving machine learning model performance.

dataset↗
NASA NTRS1982

Study of aerodynamic technology for single-cruise-engine V/STOL fighter/attack aircraft

Conceptual designs and analyses were conducted on two V/STOL supersonic fighter/attack aircraft. These aircraft feature low footprint temperature and pressure thrust augmenting ejectors in the wings for vertical lift, combined with a low wing loading, low wave drag airframe for outstanding cruise and supersonic performance. Aerodynamic, propulsion, performance, and mass properties were determined and are presented for each aircraft. Aerodynamic and Aero/Propulsion characteristics having the most significant effect on the success of the up and away flight mode were identified, and the certainty with which they could be predicted was defined. A wind tunnel model and test program are recommended to resolve the identified uncertainties.

Mark, L.↗
NASA NTRS2010

Calibration for Thrust and Airflow Measurements in the CE-22 Advanced Nozzle Test Facility

CE-22 facility procedures and measurements for thrust and airflow calibration obtained with choked-flow ASME nozzles are presented. Six calibration nozzles are used at an inlet total pressure from 20 to 48 psia. Throat areas are from 9.9986 to 39.986 sq. in.. Throat Reynolds number varies from 1.8 to 7.9 million. Nozzle gross thrust coefficient (CFG) uncertainty is 0.25 to 0.75 percent, with smaller uncertainly generally for larger nozzles and higher inlet total pressure. Nozzle discharge coefficient (CDN) uncertainty is 0.15 percent or less for all the data. ASME nozzle calibrations need to be done before and after research model testing to achieve these uncertainties. In addition, facility capability in terms of nozzle pressure ratio (NPR) and nozzle airflow are determined. Nozzle pressure ratio of 50 or more is obtainable at 40 psia for throat areas between 20 and 30 sq. in.. Also presented are results for two of the ASME nozzles vectored at 10deg, a dead-weight check of the vertical (perpendicular to the jet axis) force measurement, a calibration of load cell forces for the effects of facility tank deflection with tank pressure, and the calibration of the metric-break labyrinth seal.

Werner, Roger A.↗
DOE OSTI2021

OC6 Phase Ib: Validation of the CFD predictions of difference-frequency wave excitation on a FOWT semisubmersible

During the previous OC5 11 project, state-of-the-art mid-fidelity engineering tools for floating wind systems were found to consistently underpredict the nonlinear, low-frequency responses of semisubmersible offshore wind platforms, leading to substantial errors in the structural loads because of low-frequency surge and pitch resonance. To examine this underprediction, a coordinated investigation with computational fluid dynamics (CFD) simulations and model-basin experiments was carried out. Both investigations involved a fixed and simplified OC5-DeepCwind semisubmersible in bichromatic waves. The wave excitations - especially the nonlinear, difference-frequency excitation - on the structure from the CFD simulations were compared to the experimental measurements for validation, with uncertainty analyses for both the experimental and the CFD results. Furthermore, the wave excitations on each column of the semisubmersible were measured separately in the experiment, allowing the validation of the CFD results to be done on a per-column basis. Overall, the CFD predictions of the difference-frequency excitations agree with the experimental measurements, suggesting the CFD solutions can be used as a reference for tuning and improving the engineering-level tools and can provide a means to better understand the underprediction at low frequencies.

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Support of the Development of EPP Plus SMT Design Method (FY2020 Test Results)

Experiments in support of the development of the integrated Elastic–Perfectly Plastic (EPP) analysis and Simplified Model Test (SMT) design methodology continued in FY 2020. The previously developed single-bar SMT (SBSMT) technique was successfully extended to tubular shaped Alloy 617 specimens for evaluation of the effect of primary pressure load on SMT creep-fatigue (CF) properties at 950°C. A series of pressurized SMT tests were performed with various elastic follow-up factors and strain ranges. The results show that the primary pressure load reduces the SMT CF cycles. The reduction of SMT CF life due to primary load was found to be dependent on strain ranges, elastic follow up and test temperatures. These results are crucial for verification of the EPP strain range analysis in the EPP+SMT evaluation procedure. Additionally, hold time effect on the SMT CF life cycles is a determining factor for generating the EPP+SMT CF design curves. Lack of test data in the high cycle and low strain range region results in major uncertainty in developing the design curves. In FY2020, experiments were designed, and testing started on Alloy 617 at 950°C to fill the gap in experimental data. The on-going CF and additional SMT CF with elastic follow-up at low strain ranges will generate the information needed to verify the methods for extrapolation to long hold time and finalize the EPP+SMT design method.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2004

Assessing Videogrammetry for Static Aeroelastic Testing of a Wind-Tunnel Model

The Videogrammetric Model Deformation (VMD) technique, developed at NASA Langley Research Center, was recently used to measure displacements and local surface angle changes on a static aeroelastic wind-tunnel model. The results were assessed for consistency, accuracy and usefulness. Vertical displacement measurements and surface angular deflections (derived from vertical displacements) taken at no-wind/no-load conditions were analyzed. For accuracy assessment, angular measurements were compared to those from a highly accurate accelerometer. Shewhart's Variables Control Charts were used in the assessment of consistency and uncertainty. Some bad data points were discovered, and it is shown that the measurement results at certain targets were more consistent than at other targets. Physical explanations for this lack of consistency have not been determined. However, overall the measurements were sufficiently accurate to be very useful in monitoring wind-tunnel model aeroelastic deformation and determining flexible stability and control derivatives. After a structural model component failed during a highly loaded condition, analysis of VMD data clearly indicated progressive structural deterioration as the wind-tunnel condition where failure occurred was approached. As a result, subsequent testing successfully incorporated near- real-time monitoring of VMD data in order to ensure structural integrity. The potential for higher levels of consistency and accuracy through the use of statistical quality control practices are discussed and recommended for future applications.

Spain, Charles V.↗
DOE OSTI2024

Online Model-Free DER Dispatch Via Adaptive Voltage Sensitivity Estimation and Chance Constrained Programming

This paper proposes an online data-driven distributed energy resource management system (DERMS) for distribution system optimal DER dispatch as well as voltage regulation. Here, the key innovation is to leverage the Local Sensitivity Factor (LSF) for transforming the DER control into a computationally efficient linear programming (LP) problem. By taking real-time measurements, the estimation of LSF eliminates the need for an accurate distribution system model as well as full nodal load information, which is difficult to achieve in practice. A robust recursive least squares method is also developed to ensure the robust estimation of LSF, which is initialized using reasonable values from model-derived LSFs. This allows the system to adapt to changing operational conditions effectively. A scenario-based, chance-constrained framework is further employed to ensure voltage remains within acceptable limits in the presence of measurement and estimation uncertainties. Test results on a real-world, 759-node distribution network located in western Colorado, U.S., validate the effectiveness and robustness of the proposed control approach and demonstrate its superior performance as compared to alternative methods.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1999

Measuring Greenland Ice Mass Variation With Gravity Recovery and the Climate Experiment Gravity and GPS

The response of the Greenland ice sheet to climate change could significantly alter sea level. The ice sheet was much thicker at the last glacial maximum. To gain insight into the global change process and the future trend, it is important to evaluate the ice mass variation as a function of time and space. The Gravity Recovery and Climate Experiment (GRACE) mission to fly in 2001 for 5 years will measure gravity changes associated with the current ice variation and the solid earth's response to past variations. Our objective is to assess the separability of different change sources, accuracy and resolution in the mass variation determination by the new gravity data and possible Global Positioning System (GPS) bedrock uplift measurements. We use a reference parameter state that follows a dynamic ice model for current mass variation and a variant of the Tushingham and Peltier ICE-3G deglaciation model for historical deglaciation. The current linear trend is also assumed to have started 5 kyr ago. The Earth model is fixed as preliminary reference Earth model (PREM) with four viscoelastic layers. A discrete Bayesian inverse algorithm is developed employing an isotropic Gaussian a priori covariance function over the ice sheet and time. We use data noise predicted by the University of Texas and JPL for major GRACE error sources. A 2 mm/yr uplift uncertainty is assumed for GPS occupation time of 5 years. We then carry out covariance analysis and inverse simulation using GRACE geoid coefficients up to degree 180 in conjunction with a number of GPS uplift rates. Present-day ice mass variation and historical deglaciation are solved simultaneously over 146 grids of roughly 110 km x 110 km and with 6 time increments of 3 kyr each, along with a common starting epoch of the current trend. For present-day ice thickness change, the covariance analysis using GRACE geoid data alone results in a root mean square (RMS) posterior root variance of 2.6 cm/yr, with fairly large a priori uncertainties in the parameters and a Gaussian correlation length of 350 km. Simulated inverse can successfully recover most features in the reference present-day change. The RMS difference between them over the grids is 2.8 cm/yr. The RMS difference becomes 1.1 cm/yr when both are averaged with a half Gaussian wavelength of 150 km. With a fixed Earth model, GRACE alone can separate the geoid signals due to past and current load fairly well. Shown are the reference geoid signatures of direct and elastic effects of the current trend, the viscoelastic effect of the same trend starting from 5 kyr ago, the Post Glacial Rebound (PGR), and the predicted GRACE geoid error. The difference between the reference and inverse modeled total viscoelastic signatures is also shown. Although past and current ice mass variations are allowed the same spatial scale, their geoid signals have different spatial patterns. GPS data can contribute to the ice mass determination as well. Additional information is contained in the original.

Wu, Xiao-Ping↗
NASA NTRS2002

Surface Deformation and Gravity Changes from Surface and Internal Loads

Air and space borne remote sensing have made it possible to monitor the mass and energy transport at various scales within the cryosphere-hydrosphere-atmosphere system. The recent surface mass balance (the rate of net gain of snow and ice at a geographic point) map for the Antarctic ice sheet is constructed by interpolating sparse in situ observations (about 1,800 points) with empirically calibrated satellite data of passive back emission of microwaves. The digital elevation model obtained from satellite radar altimetry is used to improve the delineation of the ice flow drainage basins. As important as these results are, the uncertainty remains up to about 2 mm/yr of eustatic sea level change with the net imbalance. In other words, we are still unable to determine even the sign of the contribution of the Antarctic ice sheet to contemporary sea level change. The problem is more likely with the discharge rather than accumulation.

Hager, Bradford H.↗
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↗