Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “source term estimation”

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 667 records · Page 37

Modeling the X-Ray Process, and X-Ray Flaw Size Parameter for POD Studies

Nondestructive evaluation (NDE) method reliability can be determined by a statistical flaw detection study called probability of detection (POD) study. In many instances the NDE flaw detectability is given as a flaw size such as crack length. The flaw is either a crack or behaving like a crack in terms of affecting the structural integrity of the material. An alternate approach is to use a more complex flaw size parameter. The X-ray flaw size parameter, given here, takes into account many setup and geometric factors. The flaw size parameter relates to X-ray image contrast and is intended to have a monotonic correlation with the POD. Some factors such as set-up parameters including X-ray energy, exposure, detector sensitivity, and material type that are not accounted for in the flaw size parameter may be accounted for in the technique calibration and controlled to meet certain quality requirements. The proposed flaw size parameter and the computer application described here give an alternate approach to conduct the POD studies. Results of the POD study can be applied to reliably detect small flaws through better assessment of effect of interaction between various geometric parameters on the flaw detectability. Moreover, a contrast simulation algorithm for a simple part-source-detector geometry using calibration data is also provided for the POD estimation.

Khoshti, Ajay↗

Modeling the X-ray Process, and X-ray Flaw Size Parameter for POD Studies

Nondestructive evaluation (NDE) method reliability can be determined by a statistical flaw detection study called probability of detection (POD) study. In many instances, the NDE flaw detectability is given as a flaw size such as crack length. The flaw is either a crack or behaving like a crack in terms of affecting the structural integrity of the material. An alternate approach is to use a more complex flaw size parameter. The X-ray flaw size parameter, given here, takes into account many setup and geometric factors. The flaw size parameter relates to X-ray image contrast and is intended to have a monotonic correlation with the POD. Some factors such as set-up parameters, including X-ray energy, exposure, detector sensitivity, and material type that are not accounted for in the flaw size parameter may be accounted for in the technique calibration and controlled to meet certain quality requirements. The proposed flaw size parameter and the computer application described here give an alternate approach to conduct the POD studies. Results of the POD study can be applied to reliably detect small flaws through better assessment of effect of interaction between various geometric parameters on the flaw detectability. Moreover, a contrast simulation algorithm for a simple part-source-detector geometry using calibration data is also provided for the POD estimation.

Koshti, Ajay M.↗

Mitigation Methods for Piperazine Oxidation in Post-Combustion Carbon Capture

Piperazine is a promising second-generation solvent for amine scrubbing in post-combustion CO2 capture. However, the oxidative degradation of PZ can cause environmental problems and economic loss. This work presents the effects of two mitigation methods: carbon treating and N2 sparging, on the PZ oxidation in long-term operations. The species and their respective quantities adsorbed by the activated carbon were tested in a bench-scale device. The carbon was then tested in the High Temperature Oxidation Reactor (HTOR), where the solvent oxidizes at a reasonably fast rate. Pilot plant campaigns were also performed at the UT Austin SRP and the National Carbon Capture Center (NCCC) and the effects were verified. When dissolved Fe is removed by carbon treating or other methods, available soluble Fe in the system dissolves and replaces the dissolved Fe. Therefore, all available Fe and ligands need to be removed for the mitigation to be effective. The sources of soluble Fe include fly ash and the corrosion of stainless steel, and the ligands are degradation products. 8 All PZ degraded solvents have two absorbance peaks at 320 nm and 538 nm. The 320 nm peak is caused by dissolved metals, especially Fe, complexed by degradation products. The 320 nm peak is related to the amine degradation level and can be used as a simple and efficient method to estimate the amine degradation rate. The pilot plant data suggest that NO2 can oxidize PZ significantly, possibly through radical reactions. 0.01 mmol/kg-hr absorption of NO2 increased the PZ oxidation rate from 1.2 mmol/kg-hr to 2.5 mmol/kg-hr.

Yuying, Wu↗

Territory and Population Attributes Affect Florida Scrub-Jay Fecundity in Fire Adapted Ecosystems

Fecundity, the number of young produced by a breeding pair during a breeding season, is a primary component in evolutionary and ecological theory and applications. Fecundity can be influenced by many environmental factors and requires long-term study due to the range of variation in ecosystem dynamics. Fecundity data often include a large proportion of zeros when many pairs fail to produce any young during a breeding season due to nest failure or when all young die independently after fledging. We conducted color banding and monthly censuses of Florida scrub-jays (Aphelocoma coerulescens) across 31 years, 15 populations, and 761 territories along central Florida’s Atlantic coast. We quantified how fecundity (juveniles/pair-year) was influenced by habitat quality, presence/absence of nonbreeders, population density, breeder experience, and rainfall, with a zero-inflated Bayesian hierarchical model including both a Bernoulli (e.g., brood success) and a Poisson (counts of young) submodel, and random effects for year, population, and territory. The results identified the importance of increasing “strong” quality habitat, which was a mid-successional state related to fire frequency and extent, because strong territories, and the proportion of strong territories in the overall population, influenced fecundity of breeding pairs. Populations subject to supplementary feeding also had greater fecundity. Territory size, population density, breeder experience, and rainfall surprisingly had no or small effects. Different mechanisms appeared to cause annual variation in fecundity, as estimates of random effects were not correlated between the success and count submodels. The increased fecundity for pairs with nonbreeders, compared to pairs without, identified empirical research needed to understand how the proportion of low-quality habitats influences population recovery and sustainability, because dispersal into low-quality habitats can drain nonbreeders from strong territories and decrease overall fecundity. We also describe how long term study resulted in reversals in our understanding because of complications involving habitat quality, sociobiology, and population density.

Long-term studies↗

The Global Sweep of Pollution: How You Think About a Problem Determines How You Will Respond

NASA is the preeminent source of data with synoptic coverage of the globe. NASA's budget for measuring attributes of the Earth is far larger than any other single organization. Satellites provide geophysical measurements such as radiance or emittance. These measurements can be used directly or can be used as input to numeric models, which in turn produce estimates of other parameters. In general the emphasis in NASA has been to use this capability to address basic science questions. The Agency is also required to apply these data and models to help improve lives and the economy. NASA's Applied Science Program has the mandate for directing and funding such applied work. The managers of the program have recognized the disciplines of air quality, pollution, atmospheric composition and public health are aspects of a single topic with local, regional, national and global implications. Significant effort is being made to bridge the disparate disciplines and deal with the multiple scales at NASA's Marshall Space Flight Center. A central portion of our strategic approach is to integrate our personnel with short-term (proposal) colleagues and with long-term, essentially permanent, partners. We consider such broad and inclusive viewpoints as important to making real, substantive differences in people's lives.

Rickman, D. L.↗

Estimate of Gamma Dose Rates from Arrays of Fermi-1 Blanket Elements During the MEDE Process

The Enrico Fermi Atomic Power Plant Unit 1 (Fermi-1) was a sodium-cooled fast breeder reactor located in Monroe County, Michigan. The reactor was powered by a core of enriched uranium metal alloy driver fuel, which was enveloped by an axial and radial blanket material consisting of depleted uranium metal alloyed with 3 wt.% molybdenum. There are 406 axial and 559 radial irradiated sodium-bonded Fermi-1 blanket assemblies in storage at INL, totaling 34 metric tons of heavy metal. Disposal of the Fermi-1 blanket material directly into a geological repository is prohibited due to the reactive characteristic of its bond sodium. A melt drain evaporate (MEDE) process can effectively remove bond sodium from Fermi-1 blanket material. Consequently, planning is underway to apply a MEDE process to treat the 34 metric tons of heavy metal of irradiated Fermi-1 blanket material. Given the irradiated Fermi-1 blanket material’s relatively low power history and 50+ years of decay, it can be treated in a shielded glovebox. To assess the requisite shielding, the dose rates of the Fermi-1 blanket elements in various process configurations are needed. The method to perform the dose rate calculations is to first generate an average source term and associated photon source spectra for the Fermi-1 blanket material in SCALE and then use the associated photon spectra to calculate dose rates using MCNP6.2 in seven representative geometries. Dose rates in rem/h were calculated on contact (1 cm from outer geometry surface), 30 cm away from outer geometry surface, and 1 m from outer geometry surface at axial heights spanning the length of the blanket material within the geometry. The maximum average dose rate for the single Fermi-1 radial blanket element is ~0.60 rem/h, for the Fermi-1 radial blanket assembly ~1.10 rem/h, for the single Fermi-1 axial blanket element ~0.5 rem/h, for the Fermi-1 axial blanket assembly ~1.01 rem/h, for the MEDE can ~ 0.98 rem/h, for the Fermi Storage Canister ~ 0.83 rem/h, and for the MEDE cans in the DOE Standard Canister ~0.51 rem/h. Based on these maximum average values, the bounding dose rates are assumed to be 3× the average dose rates calculated for the axial blanket material and 7.6× the average dose rates calculated for the radial blanket material, based on the distribution of Fermi-1 blanket material assembly burnup.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The interpretation of crustal dynamics data in terms of plate motions and regional deformation near plate boundaries

During our participation in the NASA Crustal Dynamics Project under NASA contract NAS-27339 and grant NAG5-814 for the period 1982-1991, we published or submitted for publication 30 research papers and 52 abstracts of presentations at scientific meetings. In addition, five M.I.T. Ph.D. students (Eric Bergman, Steven Bratt, Dan Davis, Jeanne Sauber, Anne Sheehan) were supported wholly or in part by this project during their thesis research. Highlights of our research progress during this period include the following: application of geodetic data to determine rates of strain in the Mojave block and in central California and to clarify the relation of such strain to the San Andreas fault and Pacific-North American plate motions; application of geodetic data to infer post seismic deformation associated with large earthquakes in the Imperial Valley, Hebgen Lake, Argentina, and Chile; determination of the state of stress in oceanic lithosphere from a systematic study of the centroid depths and source mechanisms of oceanic intraplate earthquakes; development of models for the state of stress in young oceanic regions arising from the differential cooling of the lithosphere; determination of the depth extent and rupture characteristics of oceanic transform earthquakes; improved determination of earthquake slip vectors in the Gulf of California, an important data set for the estimation of Pacific-North American plate motions; development of models for the state of stress and mechanics of fold-and-thrust belts and accretionary wedges; development of procedures to invert geoid height, residual bathymetry, and differential body wave travel time residuals for lateral variations in the characteristic temperature and bulk composition of the oceanic upper mantle; and initial GPS measurements of crustal deformation associated with the Imperial-Cerro Prieto fault system in southern California and northern Mexico. Full descriptions of the research conducted on these topics may be found in the Semi-Annual status Reports submitted regularly to NASA over the course of this project and in the publications listed.

Solomon, Sean C.↗

Ultrasonic measurement of temperature distributions in extreme environments: Electrical power plants testing in utility-scale steam generators

Thermal heterogeneities within energy conversion and storage, material processing, nuclear processes, aerospace, and military applications are often inaccessible to characterization by insertion sensors. When sensor deployment is possible, conventional pointwise temperature probes quickly degrade when inserted into harsh environments typical of such processes. We developed spatially-resolved ultrasonic thermometry to noninvasively measure the spatial distributions of thermal properties in such applications, even when sizable thermal gradients are present. Our method divides the path of ultrasonic propagation into segments bound by echogenic features, which create echoes in pulse-echo mode, encoding the information about interior temperature distributions. We use the acquired ultrasonic responses to estimate the internal temperature distributions by solving an inverse problem or concatenating segmental estimates. This work describes the implementation and industrial testing of the developed method at a coal-fired electrical power generation plant. We inserted an echogenically segmented Inconel 625 waveguide into the combustion zone of the utility-scale boiler and continuously acquired ultrasonic data while keeping sensitive components away from the damaging combustion environment. The accuracy of the time-dependent temperature distributions reconstructed from the ultrasonic measurements was comparable to that of thermocouples. The resiliency of ultrasonic thermometry to harsh combustion conditions was far superior to conventional insertion sensors. The measurements obtained during plant operation captured daily steam generation cycles in response to changing customer demand and intermittent contributions of renewable power sources to the power grid. These measurements have revealed new insights into the relationship between the dynamic power generation load and the conditions inside the steam generator. As a result, the successful industrial testing of spatially-resolved ultrasonic thermometry in solids indicates that the developed technology has matured to become an attractive alternative to conventional sensing in solving challenging problems of long-term thermal characterizations in extreme environments.

42 ENGINEERING↗

Reducing On-Board Computer Propagation Errors Due to Omitted Geopotential Terms by Judicious Selection of Uploaded State Vector

Several future, and some current missions, use an on-board computer (OBC) force model that is very limited. The OBC geopotential force model typically includes only the J(2), J(3), J(4), C(2,2) and S(2,2) terms to model non-spherical Earth gravitational effects. The Tropical Rainfall Measuring Mission (TRMM), Wide-field Infrared Explorer (WIRE), Transition Region and Coronal Explorer (TRACE), Submillimeter Wave Astronomy Satellite (SWAS), and X-ray Timing Explorer (XTE) all plan to use this geopotential force model on-board. The Solar, Anomalous, and Magnetospheric Particle Explorer (SAMPEX) is already flying this geopotential force model. Past analysis has shown that one of the leading sources of error in the OBC propagated ephemeris is the omission of the higher order geopotential terms. However, these same analyses have shown a wide range of accuracies for the OBC ephemerides. Analysis was performed using EUVE state vectors that showed the EUVE four day OBC propagated ephemerides varied in accuracy from 200 m. to 45 km. depending on the initial vector used to start the propagation. The vectors used in the study were from a single EUVE orbit at one minute intervals in the ephemeris. Since each vector propagated practically the same path as the others, the differences seen had to be due to differences in the inital state vector only. An algorithm was developed that will optimize the epoch of the uploaded state vector. Proper selection can reduce the previous errors of anywhere from 200 m. to 45 km. to generally less than one km. over four days of propagation. This would enable flight projects to minimize state vector uploads to the spacecraft. Additionally, this method is superior to other methods in that no additional orbit estimates need be done. The definitive ephemeris generated on the ground can be used as long as the proper epoch is chosen. This algorithm can be easily coded in software that would pick the epoch within a specified time range that would minimize the OBC propagation error. This techniques should greatly improve the accuracy of the OBC propagation on-board future spacecraft such as TRMM, WIRE, SWAS, and XTE without increasing complexity in the ground processing.

Greatorex, Scott↗

Automatically generated model for light alkene combustion

Light alkenes are common combustion intermediates for a variety of fuels. Therefore, understanding their oxidation and pyrolysis chemistry is key to building detailed mechanisms for heavier fuels. This work was focused on the development and evaluation of a detailed kinetic mechanism suitable for the combustion of light alkenes up to C 4 without the use of tuned parameters, instead, the parameter values come from first principles or direct measurements. The generated mechanism accurately estimates the laminar burning velocity (S u ) and ignition delay time (IDT) of light alkenes available in the literature, which represent fundamental combustion properties at a wide range of conditions. Because each parameter is thought to have a physically realistic value, not tuned to these measurements, the new model could be used as a sub-mechanism in models for other applications. The reaction network was generated with the open-source Reaction Mechanism Generator (RMG) software. Sensitivity analyses were performed under wide ranges of temperatures and pressures, allowing for the identification of the most impactful species and reactions. Based on these, a comprehensive thermochemistry database, including calculations on 550 molecules performed in this work at the CBS-QB3 level of theory, and a kinetic library, including theoretically-derived reaction rates retrieved from the literature, were built and used in the mechanism generation. The developed mechanism was compared against several existing detailed kinetic mechanisms for ethene, propene, 1-butene, 2-butene, and isobutene. Here, the newly generated model is the most accurate among the ones analyzed, in terms of fractional bias and normalized mean square error. Hence, this new model was used to analyze the chemistry of alkene combustion. Key rate coefficients were compared, to identify the cause of deviations between the models and possible areas for further improvements.

09 BIOMASS FUELS↗

Initial Validation of NDVI time seriesfrom AVHRR, VEGETATION, and MODIS

The paper will address Theme 7: Multi-sensor opportunities for VEGETATION. We present analysis of a long-term vegetation record derived from three moderate resolution sensors: AVHRR, VEGETATION, and MODIS. While empirically based manipulation can ensure agreement between the three data sets, there is a need to validate the series. This paper uses atmospherically corrected ETM+ data available over the EOS Land Validation Core Sites as an independent data set with which to compare the time series. We use ETM+ data from 15 globally distributed sites, 7 of which contain repeat coverage in time. These high-resolution data are compared to the values of each sensor by spatially aggregating the ETM+ to each specific sensors' spatial coverage. The aggregated ETM+ value provides a point estimate for a specific site on a specific date. The standard deviation of that point estimate is used to construct a confidence interval for that point estimate. The values from each moderate resolution sensor are then evaluated with respect to that confident interval. Result show that AVHRR, VEGETATION, and MODIS data can be combined to assess temporal uncertainties and address data continuity issues and that the atmospherically corrected ETM+ data provide an independent source with which to compare that record. The final product is a consistent time series climate record that links historical observations to current and future measurements.

Morisette, Jeffrey T.↗

Application of Synoptic Magnetograms for Prediction of Solar Activity Using Ensemble Kalman Filter

Solar activity predictions using the data assimilation approach have demonstrated great potential to build reliable long-term forecasts of solar activity. In particular, it has been shown that the Ensemble Kalman Filter (EnKF) method applied to a non-linear dynamo model is capable of predicting solar activity up to one sunspot cycle ahead in time, as well as estimating the properties of the next cycle a few years before it begins. These developments assume an empirical relationship between the mean toroidal magnetic field flux and the sunspot number. Estimated from the sunspot number series, variations of the toroidal field have been used to assimilate the data into the Parker-Kleeorin-Ruzmakin (PKR) dynamo model by applying the EnKF method. The dynamo model describes the evolution of the toroidal and poloidal components of the magnetic field and the magnetic helicity. Full-disk magnetograms provide more accurate and complete input data by constraining both the toroidal and poloidal global field components, but these data are available only for the last four solar cycles. In this presentation, using the available magnetogram data, we discuss development of the methodology and forecast quality criteria (including forecast uncertainties and sources of errors). We demonstrate the influence of limited time series observations on the accuracy of solar activity predictions. We present EnKF predictions of the upcoming Solar Cycle 25 based on both the sunspot number series and observed magnetic fields and discuss the uncertainties and potential of the data assimilation approach.

Kitiashvili, Irina N.↗

Water use in a young Pinus taeda bioenergy plantation: Effect of intensive management on stand evapotranspiration

Abstract The increasing demand for plant‐derived bioenergy is projected to expand tree plantations with intensive silviculture and improved tree genetics. These silvicultural practices result in faster stand development and canopy closure, which may also influence the systems' water dynamics. Here, we studied the evapotranspiration (ET) of a young (5 years old) intensively managed loblolly pine ( Pinus taeda ) stand and investigated the components of ET to determine its contribution to overall water use. We also compared ET with plantations that received less intensive management to determine whether our stand used more water. We used the eddy covariance method to estimate ecosystem‐level total ET (ET EC ), while plot‐level estimates of ET (ET P ) were obtained via soil lysimeters, sap flow sensors, and throughfall collectors, enabling measurement of the components of ET. Soil evaporation (Es) was the largest component of ET P (36%) over the course of the study, while transpiration and canopy interception accounted for 27% and 22%, respectively. Es decreased with stand development, while transpiration and canopy interception increased. Leaf area index (LAI) and precipitation were the most significant factors controlling ET and its components. Compared to previous studies in different sites that have similar age but lower LAI, our stand had higher water use. This high water use in the early stages of stand development was primarily due to high Es before the canopy was fully developed. While there are potential sources of uncertainty when comparing ET EC and the component fluxes in ET P , results from the two methods were not significantly different. This study had the advantage of using multiple methods to understand and verify the component processes that contribute to ET. Therefore, we recommend that multiple measurement techniques be used in the long‐term observation of ET, and in particular for the evaluation of the impact that intensively managed forests have on water resources in the southeastern United States.

54 ENVIRONMENTAL SCIENCES↗

Quantification of oil and gas methane emissions in the Delaware and Marcellus basins using a network of continuous tower-based measurements

Abstract. According to the United States Environmental Protection Agency (US EPA), emissions from oil and gas infrastructure contribute 30 % of all anthropogenic methane (CH4) emissions in the US. Studies in the last decade have shown emissions from this sector to be substantially larger than bottom-up assessments, including the EPA inventory, highlighting both the increased importance of methane emissions from the oil and gas sector in terms of their overall climatological impact and the need for independent monitoring of these emissions. In this study we present continuous monitoring of regional methane emissions from two oil and gas basins using tower-based observing networks. Continuous methane measurements were taken at four tower sites in the northeastern Marcellus basin from May 2015 through December 2016 and five tower sites in the Delaware basin in the western Permian from March 2020 through April 2022. These measurements, an atmospheric transport model, and prior emission fields are combined using an atmospheric inversion to estimate monthly methane emissions in the two regions. This study finds the mean overall emission rate from the Delaware basin during the measurement period to be 146–210 Mg CH4 h−1 (energy-normalized loss rate of 1.1 %–1.5 %, gas-normalized rate of 2.5 %–3.5 %). Strong temporal variability in the emissions was present, with the lowest emission rates occurring during the onset of the COVID-19 pandemic. Additionally, a synthetic model–data experiment performed using the Delaware tower network shows that the presence of intermittent sources is not a significant source of uncertainty in monthly quantification of the mean emission rate. In the Marcellus, this study finds the overall mean emission rate to be 19–28 Mg CH4 h−1 (gas-normalized loss rate of 0.30 %–0.45 %), with relative consistency in the emission rate over time. These totals align with aircraft top-down estimates from the same time periods. In both basins, the tower network was able to constrain monthly flux estimates within ±20 % uncertainty in the Delaware and ±24 % uncertainty in the Marcellus. The results from this study demonstrate the ability to monitor emissions continuously and detect changes in the emissions field, even in a basin with relatively low emissions and complex background conditions.

54 ENVIRONMENTAL SCIENCES↗

Quantifying the Aerosol Semi-Direct Effect in the NASA GEOS-5 AGCM

Aerosols such as black carbon, dust, and some organic carbon species both scatter and absorb incoming solar radiation. This direct aerosol radiative forcing (DARF) redistributes solar energy both by cooling the surface and warming the atmosphere. As a result, these aerosols affect atmospheric stability and cloud cover (the semi-direct effect, or SDE). Furthermore, in regions with persistent high loadings of absorbing aerosols (e.g. Asia), regional circulation patterns may be altered, potentially resulting in changes in precipitation patterns. Here we investigate aerosol-climate coupling using the NASA Goddard Earth Observing System model version 5 (GEOS-5) atmospheric general circulation model (AGCM), in which we have implemented an online version of the Goddard Chemistry, Aerosol, Radiation and Transport (GOCART) model. GOCART includes representations of the sources, sinks, and chemical transformation of externally mixed dust, sea salt, sulfate, and carbonaceous aerosols. We examine a series of free-running ensemble climate simulations of the present-day period (2000-2009) forced by observed sea surface temperatures to determine the impact of aerosols on the model climate. The SDE and response of each simulation is determined by differencing with respect to the control simulation (no aerosol forcing). In a free-running model, any estimate of the SDE includes changes in clouds due both to atmospheric heating from aerosols and changes in circulation. To try and quantify the SDE without these circulation changes we then examine the DARF and SDE in GEOS-5 with prescribed meteorological analyses introduced by the MERRA analysis. By doing so, we are able to examine changes in model clouds that occur on shorter scales (six hours). In the GEOS-5 data assimilation system (DAS), the analysis is defined as the best estimate of the atmospheric state at any given time, and it is determined by optimally combining a first-guess short-term GCM forecast with all available observations. The Incremental Analysis Update (IAU) is added to the model forecast tendencies to align them with the analysis every six hours, thus preventing longer timescale feedbacks due to the aerosol forcing. We calculate the SDE by comparing model runs with and without aerosols, and the difference in the IAU between these runs is a useful metric with which to evaluate the impact of the SDE on the model atmosphere and clouds. Decreasing the IAU indicates that the aerosol direct and semi-direct effects act to reduce the bias between the model and observations and vice versa.

Randles, Cynthia A.↗

Version 2 Ozone Monitoring Instrument SO2 Product (OMSO2 V2): New Anthropogenic SO2 Vertical Column Density Dataset

The Ozone Monitoring Instrument (OMI) has been providing global observations of SO2 pollution since 2004. Here we introduce the new anthropogenic SO2 vertical column density (VCD) dataset in the version 2 OMI SO2 product (OMSO2 V2). As with the previous version (OMSO2 V1.3), the new dataset is generated with an algorithm based on principal component analysis of OMI radiances, but features several updates. The most important among those is the use of expanded lookup tables and model a priori profiles to estimate SO2 Jacobians for individual OMI pixels, in order to better characterize pixel-to-pixel variations in SO2 sensitivity, including over snow and ice. Additionally, new data screening and spectral fitting schemes have been implemented to improve the quality of the spectral fit. As compared with the planetary boundary layer SO2 dataset in OMSO2 V1.3, the new dataset has substantially better data quality, especially over areas that are relatively clean or affected by the south Atlantic anomaly. The updated retrievals over snow/ice yield more realistic seasonal changes in SO2 at high latitudes and offer enhanced sensitivity to sources during wintertime. An error analysis has been conducted to assess uncertainties in SO2 VCDs from both the spectral fit and Jacobian calculations. The uncertainties from spectral fitting are reflected in SO2 slant column densities (SCDs) and largely depend on the signal-to-noise ratio of the measured radiances, as implied by the generally smaller SCD uncertainties over clouds or for smaller solar zenith angles. The SCD uncertainties for individual pixels are estimated to be~0.15-0.3 DU (Dobson Units) between ~40°S and ~40°N and to be~0.2-0.5 DU at higher latitudes. The uncertainties from the Jacobians are approximately ~50-100% over polluted areas, and primarily attributed to errors in SO2 a priori profiles and cloud pressures, as well as the lack of explicit treatment for aerosols. Finally, the daily mean and median SCDs over the presumably SO2-free equatorial East Pacific have increased by only~0.0035 DU and ~0.003 DU respectively over the entire 15-year OMI record; while the standard deviation of SCDs has grown by only~0.02 DU or ~10%. Such remarkable long-term stability makes the new dataset particularly suitable for detecting regional changes in SO2 pollution.

OMI, SO2, Remote Sensing↗

Economic Model for Estimation of GDP Losses in the MACCS Offsite Consequence Analysis Code

The MACCS (MELCOR Accident Consequence Code System) code is the U.S. Nuclear Regulatory Commission (NRC) tool used to perform probabilistic health and economic consequence assessments for atmospheric releases of radionuclides. It is also used by international organizations, both reactor owners and regulators. It is intended and most commonly used for hypothetical accidents that could potentially occur in the future rather than to evaluate past accidents or to provide emergency response during an ongoing accident. It is designed to support probabilistic risk and consequence analyses and is used by the NRC, U.S. nuclear licensees, the Department of Energy, and international vendors, licensees, and regulators. This report describes the modeling framework, implementation, verification, and benchmarking of a GDP-based model for economic losses that has recently been developed as an alternative to the original cost-based economic loss model in MACCS. The GDP-based model has its roots in a code developed by Sandia National Laboratories for the Department of Homeland Security to estimate short-term losses from natural and manmade accidents, called the Regional Economic Accounting analysis tool (REAcct). This model was adapted and modified for MACCS and is now called the Regional Disruption Economic Impact Model (RDEIM). It is based on input-output theory, which is widely used in economic modeling. It accounts for direct losses to a disrupted region affected by an accident, indirect losses to the national economy due to disruption of the supply chain, and induced losses from reduced spending by displaced workers. RDEIM differs from REAcct in its treatment and estimation of indirect loss multipliers, elimination of double counting associated with inter-industry trade in the affected area, and that it is designed to be used to estimate impacts for extended periods that can occur from a major nuclear reactor accident, such as the one that occurred at the Fukushima Daiichi site in Japan. Most input-output models do not account for economic adaptation and recovery, and in this regard RDEIM differs from its parent, REAcct, because it allows for a user-definable national recovery period. Implementation of a recovery period was one of several recommendations made by an independent peer review panel to ensure that RDEIM is state-of-practice. For this and several other reasons, RDEIM differs from REAcct. Both the original and the RDEIM economic loss models account for costs from evacuation and relocation, decontamination, depreciation, and condemnation. Where the original model accounts for an expected rate of return, based on the value of property, that is lost during interdiction, the RDEIM model instead accounts for losses of GDP based on the industrial sectors located within a county. The original model includes costs for disposal of crops and milk that the RDEIM model currently does not, but these costs tend to contribute insignificantly to the overall losses. This document discusses three verification exercises to demonstrate that the RDEIM model is implemented correctly in MACCS. It also describes a benchmark study at five nuclear power plants chosen to represent the spectrum of U.S. commercial sites. The benchmarks provide perspective on the expected differences between the RDEIM and the original cost-based economic loss models. The RDEIM model is shown to consistently predict larger losses than the original model, probably in part because it accounts for national losses by including indirect and induced losses; whereas, the original model only accounts for regional losses. Nonetheless, the RDEIM model predicts losses that are remarkably consistent with the original cost-based model, differing by 16% at most for the five sites combined with three source terms considered in this benchmark.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Analysis of near-field and far-field aerosol dispersion for microreactors

The current paper presents a simulation-based analysis of aerosol dispersion in the near-field and far-field of a generic, conceptual microreactor operating at pressures close to the ambient pressure; therefore, in the event of an accident that causes radionuclide leakage from the microreactor containment, the radionuclide particles are less likely to travel too far from the reactor, as opposed to conventional reactors. Accordingly, the presented work provides estimates of average and 95-percentile values of the relative effluent concentration. A parametric study is then performed to narrow down the parameters which affect the aerosol dispersion characteristics most significantly. Simulations were performed in the computer code ARCON96, and the parameters found to affect aerosol dispersion characteristics are the atmospheric stability class, and the distance between the release point and the receptor. It is recommended that the computer code RADTRAD be used to calculate actual dosage over distance, using the outputs from ARCON96 as inputs, along with reactor-specific core term inventories. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗