Search NASA⌕ Search

SEARCH · Search NASA

Results for “Source term calculation”

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 595 records · Page 33

Suite of Comprehensive Rapid Analysis Tools for Estimation (SOCRATES) - 19076

Long-term monitoring and remedy implementation require ongoing assessment and analytics, much of which relies on intensive and time-consuming data integration and engineering, data manipulation, and 'big-data' analytics. Many of these analyses, such as groundwater flow direction, contaminant rate attenuation, and non-parametric statistics, are spatial-temporal in nature and involve coalescing disparate data sets from a variety of sources with measurement frequencies from milliseconds to years, as well as total data volumes often exceeding the petabyte scale. Web-based analytical tools are advantageous for ongoing real-time or near-real-time ingestion of data from various sources for general data access and visualization, and also allow for complex, data rich calculations by leveraging expansive networks and virtually unlimited compute power. A web-based framework offers a natural architecture for such endeavors, as access to compute nodes through auto-scaling and dynamic frameworks enables speedy data manipulation and unlimited access to any data set that is made available via an http request. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Investigating old-growth ponderosa pine physiology using tree-rings, δ 13 C, δ 18 O, and a process-based model

In dealing with predicted changes in environmental conditions outside those experienced today, forest managers and researchers rely on process-based models to inform physiological processes and predict future forest growth responses. The carbon and oxygen isotope ratios of tree-ring cellulose (d13Ccell, d18Ocell) reveal long-term, integrated physiological responses to environmental conditions. We incorporated a submodel of d18Ocell into the widely used Physiological Principles in Predicting Growth (3-PG) model for the first time, to complement a recently added d13Ccell submodel. We parameterized the model using previously reported stand characteristics and long-term trajectories of tree-ring growth, d13Ccell, and d18Ocell collected from the Metolius AmeriFlux site in central Oregon (upland trees). We then applied the parameterized model to a nearby set of riparian trees to investigate the physiological drivers of differences in observed basal area increment (BAI) and d13Ccell trajectories between upland and riparian trees. The model showed that greater available soil water and maximum canopy conductance likely explain the greater observed BAI and lower d13Ccell of riparian trees. Unexpectedly, both observed and simulated d18Ocell trajectories did not differ between the upland and riparian trees, likely due to similar d18O of source water isotope composition. The d18Ocell submodel with a Peclet effect improved model estimates of d18Ocell because its calculation utilizes 3-PG growth and allocation processes. Because simulated stand-level transpiration (E) is used in the d18O submodel, aspects of leaf-level anatomy such as the effective path length for transport of water from the xylem to the sites of evaporation could be estimated.

54 ENVIRONMENTAL SCIENCES↗

Effects of solar radiation on the orbits of small particles

A modification of the Robertson (1937) equations of particle motion in the presence of solar radiation is developed which allows for partial reflection of sunlight as a result of rapid and varying particle rotations caused by interaction with the solar wind. The coefficients and forces in earlier forms of the equations are compared with those in the present equations, and secular rates of change of particle orbital elements are determined. Orbital dimensions are calculated in terms of time, probable sizes and densities of meteoric and cometary particles are estimated, and times of infall to the sun are computed for a particle moving in an almost circular orbit and a particle moving in an elliptical orbit of high eccentricity. Changes in orbital elements are also determined for particles from a long-period sun-grazing comet. The results show that the time of infall to the sun from a highly eccentric orbit is substantially shorter than from a circular orbit with a radius equal to the mean distance in the eccentric orbit. The possibility is considered that the free orbital kinetic energy of particles drawn into the sun may be the energy source for the solar corona.

Lyttleton, R. A.↗

A comparison of observed and numerically predicted eddy kinetic energy budgets for a developing extratropical cyclone

The eddy kinetic energy budget is calculated for a 48-hour forecast of an intense occluding winter cyclone associated with a strong well-developed jet stream. The model output consists of the initialized (1200 GMT January 9, 1975) and the 12, 24, 36, and 48 hour forecast fields from the Drexel/NCAR Limited Area Mesoscale Prediction System (LAMPS) model. The LAMPS forecast compares well with observations for the first 24 hours, but then overdevelops the low-level cyclone while inadequately developing the upper-air wave and jet. Eddy kinetic energy was found to be concentrated in the upper-troposphere with maxima flanking the primary trough. The increases in kinetic energy were found to be due to an excess of the primary source term of kinetic energy content, which is the horizontal flux of eddy kinetic energy over the primary sinks, and the generation and dissipation of eddy kinetic energy.

Dare, P. M.↗

The root-zone soil moisture spectrum in a mediterranean ecosystem

Storage of water within soil pores of the root zone introduce memory effects in the dynamics of soil moisture that are considerably longer than the integral timescale of many atmospheric processes. Thus, hydro-climatic states can be “sustained” through land-surface heat and water vapor fluxes primarily because they can “feed off” on this long-term soil moisture memory. Root-zone soil moisture memory is only but one feature characterizing the spectrum of soil moisture dynamics, which is analyzed here using a combination of long-term measurements and models. In particular, the spectrum of root-zone soil moisture content in a Mediterranean ecosystem is examined using 14-years of half-hourly measurements. A distinguishing hydro-climatic feature in such ecosystems is that sources (mainly rainfall) and sinks (mainly evapotranspiration) of soil moisture are roughly out of phase with each other. For over 4 decades of time scales and 7 decades of energy, the canonical shape of the measured soil moisture spectrum is shown to be approximately Lorentzian determined by the soil moisture variance and its memory but with two exceptions: the occurrences of a peak at diurnal-to-daily time scales and a weaker peak at near annual time scales. Model calculations and spectral analysis demonstrate that diurnal and seasonal variations in hydroclimate forcing responsible for variability in evapotranspiration had minor impact on the normalized shape of the soil moisture spectrum. However, their impact was captured by adjustments in the temporal variance. These findings indicate that precipitation and not evapotranspiration variability dominates the multi-scaling properties of soil moisture variability consistent with prior climate model simulations. Furthermore, the soil moisture memory inferred by the annual peak of soil moisture (340 d) is consistent with climate model simulations, while the memory evaluated from the loss function of a linearized mass balance approach leads to a smaller value (50 d), highlighting the effect of weak non-stationarity on soil moisture variability.

54 ENVIRONMENTAL SCIENCES↗

Decision support for United States—Canada energy integration is impaired by fragmentary environmental and electricity system modeling capacity

The renewable energy transition is leading to increased electricity trade between the United States and Canada, with Canadian hydropower providing firm lower-carbon power and buffering variability of wind and solar generation in the U.S. However, long-term power purchase agreements and transborder transmission projects are controversial, with two of four proposed transmission lines between Quebec, Canada and the northeast U.S. cancelled since 2018. Here, we argue that controversies are exacerbated by a lack of open-source data and tools to understand tradeoffs of new hydropower generation and transmission infrastructure in comparison to alternatives. This gap includes impacts that incremental transmission and generation projects have on the economics of the entire system, for example, how new transmission projects affect exports to existing markets or incentivize new generation. We identify priority areas for data synthesis and model development, such as integrating linked hydropower and hydrologic interactions in energy system models and openly releasing (by utilities) or back-calculating (by researchers) hydropower generation and operational parameters. Publicly available environmental (e.g. streamflow, precipitation) and techno-economic (e.g. costs, reservoir size,) data can be used to parameterize freely usable and extensible models. Existing models have been calibrated with operational data from Canadian utilities that are not publicly available, limiting the range of scientific and commercial questions these tools have been used to answer and the range of parties that have been involved. Studies conducted using highly resolved, national-scale public data exist in other countries, notably, the United States, and demonstrate how greater transparency and extensibility can drive industry action. Improved data availability in Canada could facilitate approaches that (1) increase participation in decarbonization planning by a broader range of actors; (2) allow independent characterizations of environmental, health, and economic outcomes of interest to the public; and (3) identify decarbonization pathways consistent with community values.

13 HYDRO ENERGY↗

Methodology for sensitivity analysis, approximate analysis, and design optimization in CFD for multidisciplinary applications

The straightforward automatic-differentiation and the hand-differentiated incremental iterative methods are interwoven to produce a hybrid scheme that captures some of the strengths of each strategy. With this compromise, discrete aerodynamic sensitivity derivatives are calculated with the efficient incremental iterative solution algorithm of the original flow code. Moreover, the principal advantage of automatic differentiation is retained (i.e., all complicated source code for the derivative calculations is constructed quickly with accuracy). The basic equations for second-order sensitivity derivatives are presented; four methods are compared. Each scheme requires that large systems are solved first for the first-order derivatives and, in all but one method, for the first-order adjoint variables. Of these latter three schemes, two require no solutions of large systems thereafter. For the other two for which additional systems are solved, the equations and solution procedures are analogous to those for the first order derivatives. From a practical viewpoint, implementation of the second-order methods is feasible only with software tools such as automatic differentiation, because of the extreme complexity and large number of terms. First- and second-order sensitivities are calculated accurately for two airfoil problems, including a turbulent flow example; both geometric-shape and flow-condition design variables are considered. Several methods are tested; results are compared on the basis of accuracy, computational time, and computer memory. For first-order derivatives, the hybrid incremental iterative scheme obtained with automatic differentiation is competitive with the best hand-differentiated method; for six independent variables, it is at least two to four times faster than central finite differences and requires only 60 percent more memory than the original code; the performance is expected to improve further in the future.

Taylor, Arthur C., III↗

Iodine Effects on Medical Isotope Producer Delay & Guard Bed Temperature and Performance in Support of NA-23 Program

In support of NA-23’s objective to accelerate the establishment of a reliable commercial United States (U.S.) domestic supply of 99Mo produced without the use of highly enriched uranium (HEU), the Pacific Northwest National Laboratory (PNNL) was tasked with modeling the performance of activated carbon delay and guard beds that may be typical for non-HEU 99Mo production processes. The specific tasks accomplished by PNNL included calculations of long-term effects of iodine on delay bed performance and development of a finite element simulation to better understand delay bed performance under various scenarios specifically with radioxenon and radioiodine decay heating of the delay beds. The main issue of interest in this study is understanding the effect of iodine on the abatement performance of the beds if iodine is not sufficiently removed from gases prior to reaching the delay beds. Finite element simulations of delay bed performance were developed using COMSOL Multiphysics® Software. These simulations are based on partial differential equations that describe xenon and iodine adsorption, heat transfer and decay of radioxenon and radioiodine isotopes. Heat transfer in the activated carbon beds is simulated with a heat source tied to the decay of radioactive xenon and iodine species. These simulations were used to investigate delay bed temperature and hold up time using radioxenon and radioiodine activities and flowrates expected from a non-HEU production process. A bounding case scenario was simulated, where all the iodine was released into the delay beds due to insufficient guarding of the delay beds. In addition, preliminary investigations of the guard beds were conducted.

07 ISOTOPE AND RADIATION SOURCES↗

Deciphering Degradation: Machine Learning on Real-World Performance Data

The solar investment community is in need of up-to-date and data-derived system degradation rates to use in financial models that calculate the risk and expectation of energy revenue. To date, financial models choose degradation rates from past studies that are limited and may not be representative of solar projects in development today. As the solar industry expands in all sectors (residential, commercial, and utility), there is a growing amount of energy generation data available to analyze. With the help of the open source degradation analysis code (RdTools) it is now possible to derive degradation on an ongoing basis and continuously provide up-to-date durability statistics to the investment community. Ongoing and accurate statistics of fielded photovoltaic systems allows financial stakeholders to constrain energy revenue projections, lower financial risk, and increase the bankability of solar. This presentation will provide a mid-term progress update on a 2 year SETO-funded degradation project that involves applying RdTools analysis on kWh Analytics’ industry database. With results from the RdTools analysis in hand, machine learning methods are employed to quantify how much variance of system performance degradation can be explained by predictor variables such as environmental factors, equipment bill of materials characteristics, and system design information. The goal of this presentation is to promote the usage of RdTools on an ongoing to data owners and to solicit feedback from researchers and stakeholders on how to maximize impact of the SETO-funded kWh Analytics degradation project.

14 SOLAR ENERGY↗

Calculating n-point charge correlations in evolving systems

Background: In dynamic systems such as heavy-ion collisions, charge susceptibilities and local charge correlations change with time. These changes are accompanied by nonlocal correlations, which spread diffusively with time and are constrained by local charge conservation. Such correlations have been measured at the Relativistic Heavy-Ion Collider (RHIC) and have been linked to the chemical evolution of the matter and possible phase separation or critical phenomena. The evolution of such correlations has been modeled for two-point correlators superimposed onto hydrodynamics, but not for three- or four-point correlators. These higher-order correlators represent an essential basis for calculating the kurtosis and skewness of charge distributions measured at RHIC. For short-range correlations, such as those between charges on the same particle, and the associated charge-balance correlations a theoretical formalism was lacking for handling higher-order correlations. In particular, one needed to understand how higher-order charge fluctuations was split onto lower numbers of particles. For example, correlations of order three have contributions with charge all on one particle, split onto two particle, or split onto three. By assuming local chemical equilibrium of short-range correlations, the sources of the various correlators were found to be uniquely determined. The relevant theoretical foundation is presented here, including a diagrammatic technique to correctly account for all the possible terms to higher-order correlations. Here, we found a consistent theoretical treatment and its applicability, viability, and tractability are discussed. The formalism derived here enables realistic and quantitative modeling of the three- and four-point charge correlations necessary for understanding measurements of the kurtosis and skewness of charge distributions at RHIC. These observables have been promoted as signals of phase transitions or measures of chemical evolution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Interseismic Strain Accumulation between the Colorado Plateau and the Eastern California Shear Zone: Implications for the Seismic Hazard near Las Vegas, Nevada

ABSTRACT Crustal deformation in the central Basin and Range between the Colorado plateau and the Eastern California Shear Zone is active but slow, making it a challenge to assess how strain is distributed and crustal motion transferred. However, knowledge of strain rates is very important, particularly for addressing the seismic hazard for both the Las Vegas urban area and the site of the proposed Yucca Mountain nuclear waste repository, in southern Nevada. Global Positioning System (GPS) data provide important constraints, particularly now that the GPS network in the area has substantially expanded in recent years. However, because deformation is slow, it is important to mitigate any transient tectonic and nontectonic signals to obtain the most accurate long-term interseismic motion and robust estimation of strain rates. We use data from all GPS stations in the region including both long-running continuous and semicontinuous stations. We model and remove postseismic displacements at these stations using source parameters for 41 events, dating back to the 1700 Cascadia megathrust earthquake, which contribute significantly to the deformation field within the central Basin and Range. We also remove correlated noise from the time series with the common-mode component imaging technique. We find that removal of both the postseismic transients and common-mode noise substantially reduces the uncertainties and spatial variation in the velocities. We find east–west extension across the Las Vegas Valley of 0.5–0.6 mm/yr. The interseismic strain rate field, calculated with the final velocities, reveals higher strain rates through southern Nevada than in previous studies, with rates within Las Vegas Valley of 8.5±2.4×10−9 yr−1. Our results also confirm shear along the Pahranagat shear zone, but the estimated amplitude is strongly affected by postseismic relaxation.

Geochemistry & Geophysics↗

Commercial Building Technology Demonstration (Final Scientific/Technical Report)

The background for this project was a unique and innovative concept that had been proposed and demonstrated by enVerid Systems, reducing the energy consumption in HVAC systems in commercial and public buildings by up to 50% at peak loads. The concept involves selective scrubbing of gas contaminants from indoor air as a practical substitute for outside air ventilation, the latter being the conventional practice today. This scrubbing enables recycling of most of the indoor air and significantly decreasing the amount of outside air, resulting in a reduced HVAC load and commensurate energy savings for the facility. Key to our approach is low cost scrubbing of all gas contaminants from indoor air using novel, efficient and regenerable sorbent materials. The total solution, labeled HLR (HVAC Load Reduction) is implemented by means of compact modules that can be retrofitted onto virtually any building, and utilizes novel sorbents that have been developed and brought to scale production by enVerid. Although the technology had been deployed and demonstrated, the US commercial building and HVAC markets, famously conservative with respect to innovation and new technologies, had yet to embrace the idea or the product. The overarching goal of the project was to help overcome market hesitancy and enable widespread adoption, by creating a critical mass of success stories across different regions and building types, all under oversight and verification of the DOE and NREL, to ensure credibility and accuracy. The project, overall, was a resounding success. Several installations were deployed and tested comprehensively for the impact on both energy consumption and indoor air quality, and results were not only excellent, but quantitatively in line with expectations. A total number of eight (8) installations were attempted, with 2 aborted due to building issues. Of the six successful installations, however, complete M&V was performed by NREL on 2 (Miami and NYC) and partial analysis on a third site (ArcBest), and these three, along with their results, were therefore presented from a financial and reporting standpoint. Although the number of case studies fell short of the original target of 10, we nevertheless had excellent demonstrations in vastly different climate zones: Southeast, Northeast and Midwest, and in very different building types. This included a modern high-rise office building in New York City and a University Campus health center in Miami, and a mid-rise high-density office in Arkansas. The cooling load reduction was shown to be significant, averaging 37% during the final cooling season in the Miami site. In the NYC site, the cooling energy savings were generally less, at around 12%, as expected. While this is a lower number for average/total savings, two comments are important to note: (i) the outside airflow (OA) reduction was not consistent during the data collection period, at least in part understating the attainable savings through incorrect measurement of the “high energy” reference state; and (ii) the peak impact in NYC is almost certainly as high as Miami, with critical ramifications for peak load reduction and equipment downsizing ability. Importantly, the measured energy savings in all locations were related to the outdoor climate conditions and the actual achievable reduction in outside air, very much in line with models and calculations; this leaves us with increased confidence in our understanding of the interaction between the technology and the building, and our ability to project the impact of this innovation in the future, through broad adoption. Simply put: the energy savings and load reduction impact of the HLR solution has been confirmed quantitatively. The results also confirmed the air quality objective. This may be even more significant, in that air quality is more complex and harder to measure and quantify, and the performance of the HLR as an air scrubber is the linchpin of the entire solution. In all sites, indoor air quality met or exceeded the required parameters, in terms of achieving the target (low) levels of indoor pollutants. Furthermore, the concentrations of outdoor sourced pollutants were markedly lower, as expected, mostly due to the reduction in outside air intake although in part also due to the scrubbing. The success was not immediate nor easy, and early challenges with the systems, especially (but not only) around controls and communication, caused delays in commissioning and undermined some of our M&V schedules. The good news is no fundamental issues were encountered and the number of such glitches was declining quickly by the end of the project. Other challenges had to do with coordinating the control of the building OA dampers, which the HLR needs to have in order to deliver the load reduction, but which the facility manager controls regularly and independently. This is not necessarily a problem for normal HLR operation but interferes with accurate M&V, where switching OA settings back and forth on a precise timetable is essential for clean data. The project provided valuable experience in retrofitting buildings that had been originally designed to operate with less recirculation and more outside air, and some key learnings related to the cyclical regeneration of the sorbents in the operational environment of the building. Some of these learnings gave already been incorporated into ongoing improvements in the design of the HLR module as well as the software that is used to control its operation. The cost and time of retrofit installations is much better understood; so is the test and balance procedure of the building that determines the available range of outside air reduction, which in turn determines the amount of load reduced. We had also implemented at least one annual cycle of sorbent cartridge replacement in the field, after which we tested the used sorbent to assess the rate of degradation of sorbent potency – an important factor affecting the long-term economics of HLR deployments. A notable shortfall of the project, other than the fact that the number of sites was below target, was that we had not (yet) completed testing to measure and verify (M&V) the winter energy savings. This is primarily due to late start on the NYC sites (Miami does not portend any heating savings associated with reduced outside air), combined with the issues of OA damper control. It is by no means suggested that such savings cannot be captured, only that the M&V with respect to winter in NYC has been incomplete. The value of this project and its interest to the general public are hard to overstate. The case studies created from Miami and especially from NYC have, without a doubt, opened up the interest of the broad HVAC industry to HLR in the past year. Major manufacturer reps have signed up to represent and promote the product, and dozens of top tier engineering consultants have introduced HLR as a base of design in upcoming construction and renovation projects – precisely the overarching goal of this project. These would not have happened without the highly visible installation in NYC and its ultimately satisfied customer. Furthermore, the DOE itself issued a report in late 2017 analyzing hundreds of candidate technologies for reducing energy consumption in HVAC; remarkably, the budding HLR solution ranked among the top 3(!); no doubt, its status as a viable and significant contender on this list is informed by the success of these demonstrations. And furthermore, in the wake of these demos, several major energy utilities, led by ConEdison in NYC, have endorsed the HLR solution and now routinely provide substantial cash rebates for customers choosing to incorporate the HLR in their buildings and reap the benefit of the reduced peak load as well as reduced cumulative energy demand. Last but not least, in January of 2019, as this report was being written, the HLR module won the AHR Expo Innovation Award for Green Buildings and, remarkably, the prestigious and only Product of the Year 2019, being selected by an independent professional panel of judges from AHR and ASHRAE. Given the confirmation of the enormous energy saving potential of HLR in widescale deployment, as well as ancillary but significant benefits in terms of equipment downsizing and indoor air quality, the benefit of this effort to the public speaks for itself. Accelerating adoption of HLR, without dependence on subsidies and driven by the private sector, is well under way. There is still much work to be done in terms of educating the HVAC and construction ecosystems and improving the product’s robustness and simplicity, to further facilitate the adoption of this technology and the capture of its potential value. But the American public will greatly benefit both directly and indirectly in the immediate future and beyond.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Preliminary Criticality and Radiation Shielding Analysis for the Storage and Transfer of MARVEL Reactor Spent Nuclear Fuel

This report documents the results of preliminary nuclear criticality and radiation shielding assessments during transfer and dry storage of irradiated Microreactor Applications Research Validation and Evaluation (MARVEL) reactor fuel at Idaho National Laboratory (INL). The assessments focus on transfer casks and storage canisters that are currently in use at INL, which may be compatible with the irradiated MARVEL reactor fuel. The criticality assessments were performed using the radiation transport code MCNP6 with 37 MARVEL reactor fuel elements in various configurations and scenarios. All transfer and storage configurations under dry conditions were below the assumed criticality safety limit of 0.93. Some storage and transfer configurations under wet conditions exceeded the criticality safety limit. This suggests that the appropriate administrative and engineering controls, in addition to reducing the number of MARVEL reactor fuel elements per transfer cask or storage canister, can be expected to ensure criticality safety under all scenarios. The radiation shielding assessments were performed by generating conservative neutron and photon source terms using the ORIGEN module in the SCALE suite of codes. These source spectra were used to estimate the dose equivalent rates using the radiation transport code MCNP6, both on contact and 1 m away from the fuel and transfer casks. The maximum estimated dose equivalent rate of 37 unshielded MARVEL reactor fuel elements on contact is approximately 42000 R/hr. The maximum estimated dose equivalent rates on contact to the ATR transfer cask, HFEF-5 transfer cask, and high load charger were approximately 233 mR/hr, 171 mR/hr, and 201 mR/hr, respectively. This suggests that with the appropriate administrative and engineering controls, all three transfer casks analyzed can be expected to provide sufficient radiation shielding to workers during transfer of irradiated MARVEL reactor fuel. These calculations are performed to support the planning and strategy for the MARVEL project and will demonstrate the technical viability of the different configurations discussed and help identify where engineered or administrative controls may be necessary. A complete criticality safety analysis and radiation shielding analysis, including validation and contingency and accident analysis must be completed by licensed and authorized personnel before any transfer or storage of MARVEL reactor nuclear fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Commercial, industrial, and institutional discount rate estimation for efficiency standards analysis: Sector-level data 1998–2023

Underlying each of the U.S. Department of Energy’s (DOE’s) federal appliance and equipment energy conservation standards are a set of complex analyses of the projected costs and benefits of regulation. Any new or amended standard must be designed to achieve significant additional energy conservation, provided that it is technologically feasible and economically justified (42 U.S.C. 6295(o)(2)(A)). DOE determines economic justification based on whether the benefits exceed the burdens, considering a variety of factors, including the economic impact of the standard on consumers of the product and the savings in lifetime operating cost compared to any increase in price or maintenance expenses (42 U.S.C. 6295(o)(2)(B)). As part of this determination, DOE conducts a life-cycle cost (LCC) analysis, which models the combined impact of appliance first cost and operating cost changes on a representative commercial building sample to identify the fraction of customers achieving LCC savings or incurring net cost at the considered efficiency levels. Thus, the commercial discount rate value(s) used to calculate the present value of energy cost savings within the LCC model implicitly plays a role in estimating the economic impact of potential standard levels. This report provides an in-depth discussion of the commercial discount rate estimation process relying on the Capital Asset Pricing Model (CAPM) to estimate a business’ cost of equity, and by adding a risk adjustment factor to the risk-free rate associated with long-term U.S. Treasury bonds to estimate their cost of debt. It is an update to previous reports on estimating commercial discount rates from firm-level and sector-level financial data (e.g., Fujita, 2021, 2016). Major topics covered in this report include the following: • Discount rate estimation methods and rationale • Data sources used and data limitations • Discount rate distributions for use in standards analysis • Discount rate estimation methods and distributions specific to the small business subgroup analysis.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Radiation force on a relativistic plasma and the Eddington limit

Calculations of the radiation force due to Thomson scattering on a relativistic plasma are presented and discussed in relation to certain astronomical objects which operate within a few orders of magnitude of their Eddington luminosity. The radiation force on a hot isotropic plasma is shown to exceed that on a cold plasma by a factor depending on the electron Lorentz factor, which can be substantial in a relativistic plasma. In such a plasma, radiative bulk acceleration is found to occur through the anisotropic loss of internal energy during the Thomson scattering, resulting in an effect termed a Compton rocket. The Compton rocket is shown to be a relevant acceleration mechanism in situations where a relativistic electron plasma lies in the vicinity of a luminous source operating near the classical Eddington limit, including compact galactic X-ray sources and objects associated with active galactic nuclei (quasars, blazars, and Seyfert nuclei).

Odell, S. L.↗

Theoretical Basis for Finite Difference Extrapolation of Sonic Boom Signatures

Calculation of sonic boom signatures for aircraft has traditionally followed the methods of Whitham' and Walkden. The wave disturbance generated by the vehicle is obtained by area rule linearized supersonic flow methods, which yield a locally axisymmetric asymptotic solution. This solution is acoustic in nature, i.e., first order in disturbance quantities, and corresponds to ray acoustics. Cumulative nonlinear distortion of the signature is incorporated by using this solution to adjust propagation speed to first order, thus yielding a solution second order in disturbance quantities. The effects of atmospheric gradients are treated by Blokhintzov's method of geometrical acoustics. Both nonlinear signature evolution and ray tracing are applied as if the pressure field very close to the vehicle were actually that given by the source term (the 'F-function') of the asymptotic linearized flow solution. The viewpoint is thus that the flow solution exists at a small radius near the vehicle, and may be treated as an input to an extrapolation procedure consisting of ray tracing and nonlinear aging. The F-function is often regarded as a representation of a near-field pressure signature, and it is common for computational implementations to treat it interchangeably with the pressure signature. There is a 'matching radius' between the source function and the subsequent propagation extrapolation. This viewpoint has been supported by wind tunnel tests of simple models, and very typically yields correct results for actual flight vehicles. The assumption that the F-function and near-field signature are interchangeable is generally not correct. The flowfield of a vehicle which is not axisymmetric contains crossflow components which are very significant at small radii and less so at larger distances. From an acoustical viewpoint, the crossflow is equivalent to source diffraction portions of the wave field. Use of the F-function as a near field signature effectively assumes that the diminution of the crossflow/diffraction component may be applied all at once at the matching radius noted above. This approximation, though not rigorously validated, is responsible for the usual correct far-field results. On the other hand, if an actual near-field signature (either from wind tunnel or CFD data) is used at a starting point rather than one based on th effective source distribution, the predicted far-field signature is generally wrong.

Plotkin, Kenneth J.↗

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↗

Photochemistry of N2(+) in the daytime F region

The photochemistry of N2(+) in the daytime F region continues to be a source of concern due primarily to the uncertain roles of the metastable O+(2D) ion and possible vibrational excitation. This investigation adopts a unique subset of data from the Atmospheric Explorer C satellite that spans distinct regions near the low-latitude F peak wherein either chemical reactions or electron processes alternatively control the N2(+) abundances. Concentrations of N2(+) calculated according to current theory and including recent laboratory data for O+(2D) losses via N2 are found to exceed ionospheric observations between 220 and 400 km by a factor near 2. Relevant characteristics of the basic molecular ion concentration measurements are described, along with an analysis for select orbits conducted according to current photochemical equilibrium theory. Possible simple modifications of this theory are discussed in terms of the constraints suggested by the nature of the observations.

Breig, E. L.↗