Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Source Term Analysis”

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 127 records · Page 7

Determination of space vehicle contamination

The contamination and degradation of externally exposed sensitive experimental and functional surfaces of the Skylab by the self-induced atmosphere surrounding the vehicle are studied. It is shown that the major contamination sources of long-term deposits on Skylab are primarily outgassing from nonmetallic coatings and rocket engine exhaust products. Analysis showed that the characteristic decay time of the outgassing source rate may be as long as 4100 hr, and that essentially 20% of RCS deposits that stick at temperatures of 0 C will remain after sublimating at a rate whose characteristic decay time is 72 hr.

Rantanen, R. O.↗

Mechanistic Source Term Considerations for Advanced Non-LWRs (Revision 1)

This report is a functional review of the radionuclide containment strategies of fluoride-salt-cooled high temperature reactor (FHR), molten salt reactor (MSR) and high temperature gas reactor (HTGR) systems. This analysis serves as a starting point for further, more in-depth analyses geared towards identifying phenomenological gaps that still exist, hindering the creation of a mechanistic source term for these reactor types. As background information to this review, an overview of how a mechanistic source term is created and used for consequence assessment necessary for licensing is provided. How a mechanistic source term is used within the Licensing Modernization Project (LMP) is also provided. Lastly, the characteristics of non-LWR mechanistic source terms are examined. This report does not assess the viability of any software system for use with advanced reactor designs, but instead covers system function requirements. Future work within the Nuclear Energy Advanced Modeling and Simulations (NEAMS) program will address such gaps. This document is an update of SAND 2020-6730. An additional chapter is included as well as edits to original content.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Analysis and Synthesis of Tonal Aircraft Noise Sources

Fixed and rotary wing aircraft operations can have a significant impact on communities in proximity to airports. Simulation of predicted aircraft flyover noise, paired with listening tests, is useful to noise reduction efforts since it allows direct annoyance evaluation of aircraft or operations currently in the design phase. This paper describes efforts to improve the realism of synthesized source noise by including short term fluctuations, specifically for inlet-radiated tones resulting from the fan stage of turbomachinery. It details analysis performed on an existing set of recorded turbofan data to isolate inlet-radiated tonal fan noise, then extract and model short term tonal fluctuations using the analytic signal. Methodologies for synthesizing time-variant tonal and broadband turbofan noise sources using measured fluctuations are also described. Finally, subjective listening test results are discussed which indicate that time-variant synthesized source noise is perceived to be very similar to recordings.

Allen, Matthew P.↗

Global Premature Mortality By Dust and Pollution PM 2.5 Estimated From Aerosol Reanalysis of the Modern-Era Retrospective Analysis for Research and Applications, Version 2

This study quantifies global premature deaths attributable to long-term exposure of ambient PM 2.5 , or PM 2.5 -attributable mortality, by dust and pollution sources. We used NASA’s Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2) aerosol reanalysis product for PM 2.5 and the cause-specific relative risk (RR) from the integrated exposure-response (IER) model to estimate global PM2.5-attributable mortality for five causes of deaths, namely ischaemic heart disease (IHD), cerebrovascular disease (CEV) or stroke, lung cancer (LC), chronic obstructive pulmonary disease (COPD), and acute lower respiratory infection (ALRI). The estimated yearly global PM 2.5 -attributable mortality in 2019 amounts to 2.89 (1.38–4.48) millions, which is composed of 1.19 (0.73–1.84) million from IHD, 1.01 (0.35–1.55) million from CEV, 0.29 (0.11–0.48) million from COPD, 0.23 (0.14–0.33) million from ALRI, and 0.17 (0.04–0.28) million from LC (the numbers in parentheses represent the estimated mortality range due corresponding to RR spread at the 95% confidence interval). The mortality counts vary with geopolitical regions substantially, with the highest number of deaths occurring in Asia. China and India account for 40% and 23% of the global PM 2.5 -attributable deaths, respectively. In terms of sources of PM 2.5 , about 22% of the global all-cause PM 2.5 -attributable deaths are caused by desert dust. The largest dust attribution is 37% for ALRI. The relative contributions of dust and pollution sources vary with the causes of deaths and geographical regions. Enforcing air pollution regulations to transfer areas from PM 2.5 nonattainment to PM2.5 attainment can have great health benefits. Being attainable with the United States air quality standard (AQS) of 15 μg/m 3 globally would have avoided nearly 40% or 1.2 million premature deaths. The most recent update of PM 2.5 guideline from 10 to 5 μg/m 3 by the World Health Organization (WHO) would potentially save additional one million lives. Our study highlights the importance of distinguishing aerodynamic size from geometric size in accurately assessing the global health burden of PM 2.5 and particularly for dust. A use of geometric size in diagnosing dust PM 2.5 from the model simulation, a common approach in current health burden assessment, could overestimate the PM 2.5 level in the dust belt by 40–170%, leading to an overestimate of global all-cause mortality by 1 million or 32%.

PM2.5↗

A high order Cartesian grid, finite volume method for elliptic interface problems

We present a higher-order finite volume method for solving elliptic PDEs with jump conditions on interfaces embedded in a 2D Cartesian grid. Second, fourth, and sixth order accuracy is demonstrated on a variety of tests including problems with high-contrast and spatially varying coefficients, large discontinuities in the source term, and complex interface geometries. We include a generalized truncation error analysis based on cell-centered Taylor series expansions, which then define stencils in terms of local discrete solution data and geometric information. In the process, we develop a simple method based on Green's theorem for computing exact geometric moments directly from an implicit function definition of the embedded interface. This approach produces stencils with a simple bilinear representation, where spatially-varying coefficients and jump conditions can be easily included and finite volume conservation can be enforced.

97 MATHEMATICS AND COMPUTING↗

A mathematical framework for ejecta cloud dynamics with application to source models and piezoelectric mass measurements

We present a mathematical framework for describing the dynamical evolution of an ejecta cloud generated by a generic ejecta source model. We consider a piezoelectric sensor fielded in the path of an ejecta cloud, for experimental configurations in which the ejecta are created at a singly shocked planar surface and fly ballistically through vacuum to the stationary sensor. To do so, we introduce the concept of a time- and velocity-dependent ejecta “areal mass function.” We derive expressions for the analytic (“true”) accumulated ejecta areal mass at the sensor and the measured (“inferred”) value obtained via the standard method for analyzing piezoelectric voltages. In this way, we derive an exact expression and upper bound for the error imposed upon a piezoelectric ejecta mass measurement (in a perfect system) by the assumption of instantaneous creation, which is commonly required for momentum diagnostic analyses. This error term is zero for truly instantaneous source models; otherwise, the standard piezoelectric analysis is guaranteed to overestimate the true mass. When combined with a piezoelectric dataset, this framework provides a unique solution for the ejecta particle velocity distribution, subject to the assumptions inherent in the data analysis. The framework also leads to strong boundary conditions that any ejecta source model must satisfy in order to be consistent with apparently global properties of piezoelectric measurements from a wide range of experiments. We demonstrate this methodology by applying it to the Richtmyer–Meshkov instability+self-similar velocity distribution ejecta source model currently under development at Los Alamos National Laboratory.

97 MATHEMATICS AND COMPUTING↗

As-Run Physics Analysis for the EPRI Zirconium Growth C Capsule

This engineering calculation and analysis report documents the as-run nuclear-heating term, source-term, and radiation damage (measured in displacements per atom [dpa]) accumulation for the Electric Power Research Institute zirconium growth experiment's C capsule (EPRI-ZG-C). These data have been calculated using MCNP 6.1 and SCALE 6.2.3. These data can be used for: the shipping of the EPRI-ZG-C capsule, the acceptance of the capsule at the Materials and Fuels Complex, and programmatic scientific needs.

36 MATERIALS SCIENCE↗

Linear ion trap based atomic frequency standard

In order to develop a trapped ion-based fieldable frequency standard with stability 1 x 10 to the -13th/sq rt tau for averaging times tau greater than 10,000 s, a hybrid RF/DC linear ion trap was developed which permits storage of large numbers of ions with reduced susceptibility to the second-order Doppler effect caused by the RF confining fields. The authors have confined Hg-199(+) ions in this trap and have measured very high Q transitions with good SNRs. In preliminary measurements they obtained stabilities of 1.6 x 10 to the -13th/sq rt tau (tau between 50 and 800 s) with a 160-mHz wide atomic resonance linewidth and a signal-to-noise ratio of 40 for each measurement cycle. Atomic resonance lines as narrow as 30 mHz on the 40.5-GHz clock transition have been measured with no appreciable reduction in the ion signal. A stability of 7 x 10 to the -14th/sq rt tau is made possible by the signal-to-noise and line Q of this measured transition. Analysis of fundamental sources of frequency instability indicates that a long-term stability of 2 x 10 to the -16th is feasible for this device with existing technology for tau = 10 to the 6th s or more.

Prestage, John D.↗

An Overview of Atmospheric Features Over the Western North Atlantic Ocean and North American East Coast – Part 1: Analysis of Aerosols, Gases, and Wet Deposition Chemistry

The Western North Atlantic Ocean (WNAO) and adjoining East Coast of North America are of great importance for atmospheric research and have been extensively studied for several decades. This broad region exhibits complex meteorological features and a wide range of conditions associated with gas and particulate species from many sources regionally and other continents. As Part 1 of a 2‐part paper series, this work characterizes quantities associated with atmospheric chemistry, including gases, aerosols, and wet deposition, by analyzing available satellite observations, ground‐based data, model simulations, and reanalysis products. Part 2 provides insight into the atmospheric circulation, boundary layer variability, three‐dimensional cloud structure, properties, and precipitation over the WNAO domain. Key results include spatial and seasonal differences in composition along the North American East Coast and over the WNAO associated with varying sources of smoke and dust and meteorological drivers such as temperature, moisture, and precipitation. Spatial and seasonal variations of tropospheric carbon monoxide and ozone highlight different pathways towards the accumulation of these species in the troposphere. Spatial distributions of speciated aerosol optical depth and vertical profiles of aerosol mass mixing ratios show a clear seasonal cycle highlighting the influence of different sources in addition to the impact of intercontinental transport. Analysis of long‐term climate model simulations of aerosol species and satellite observations of carbon monoxide confirm that there has been a significant decline in recent decades among anthropogenic constituents owing to regulatory activities.

ACTIVATE↗

Assessing High Burnup U-19Pu-10Zr Fuel Performance against Historical and Modeled Behavior

Advancing the deployment of sodium-cooled fast reactors (SFRs) requires thorough testing of metallic fuel pins under accident conditions to establish safe operational limits of high burnup fuel. To conduct transient testing, a comprehensive understanding of steady-state fuel behavior obtained through both experimental characterization and accurate predictive capabilities is needed. This study comparatively assesses the steady-state irradiation performance of two high burnup U-19Pu-10Zr fuel pins, DP-36 and DP-40, irradiated under prototypic fast reactor conditions in preparation for planned safety testing at the Transient Reactor Test Facility. Since DP-40 was designated for use in the test and DP-36 serves as its sibling pin, non-destructive, engineering-scale post-irradiation examinations (PIE) were conducted on both pins while destructive examinations were performed exclusively on DP-36. The results were then assessed against historical performance data from similar fuel pins irradiated in the Experimental Breeder Reactor-II. Additionally, the steady-state irradiation of each pin was modeled using the BISON fuel performance code to assess the accuracy of current modeling capabilities in predicting the baseline irradiation behavior. Non-destructive examinations included neutron radiography to measure fuel column elongation, gamma scanning to verify pin integrity and fission product migration, and profilometry to assess dimensional changes. Benchmarking against existing PIE data revealed consistent patterns in axial fuel column growth and cladding diametral strain, though both pins exhibited longer low-density “fluff” structures, which can have implications for core reactivity and source term calculations. Destructive examinations on DP-36 included fission gas release analysis and sectioning for optical microscopy, which showed more complex constituent redistribution patterns than the traditionally accepted 3-ring model. The axial evolution of fractional areas and porosities of each of the redistributed zones were quantified and presented. Modeling comparisons showed agreement in fractional fission gas release but consistently overestimated axial and radial swelling and disagreed with measured axial porosity patterns. These conservative overpredictions suggested that the pins would appear closer to failure or operational limits at the start of transient tests, potentially leading to higher strain accumulation during the transient. While conservative estimates provide safety margins, they can negatively impact fuel economics. A review of the swelling models identified areas for improvement in the gaseous swelling, solid swelling, and fuel hot-pressing models when applied to ternary fuel. The results of this study highlight the critical importance of conducting pre-test characterization on both test and sibling pins to accurately capture steady-state fuel behavior, providing a precise baseline for post-test evaluations and essential inputs for transient modeling of the planned experiments. The analysis also revealed significant data gaps that require further investigation to enhance the understanding and prediction of fuel swelling and pore dynamics. Collecting comprehensive data across different irradiation conditions, burnup levels, and fuel compositions are essential for refining existing models and developing mechanistic models for both binary and ternary metallic fuels, ultimately improving the integration of modeling and experimental approaches in accident testing.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Development of a BlueCRAB/MELCOR Framework for Supporting Realistic Mechanistic Source Term Calculations in Microreactors

Efforts are currently underway to deploy microreactor modeling and simulation tools to better support vendors and regulatory authorities in submitting and reviewing licensing applications. In particular, the Nuclear Regulatory Commission is expected to rely on the Comprehensive Reactor Analysis Bundle (BlueCRAB) multiphysics toolset in performing design- and beyond-design-basis accident analyses. In addition, the Nuclear Regulatory Commission has been using the MELCOR code to estimate mechanistic source terms during accidents. As MELCOR relies on isotopic inventory and reactor temperature/power evolution profiles during accident conditions—all of which can theoretically be obtained from BlueCRAB—the ultimate goal of this activity is to establish a common BlueCRAB-MELCOR framework. However, prior to the present research, BlueCRAB had never been used to calculate such quantities of interest at the full-core level. While there are many Monte Carlo (MC) codes capable of computing such quantities of interest, they are unable to readily account for multiphysics feedback. BlueCRAB allows for the coupling of different physics codes together to perform multiphysics-informed calculations. Therefore, the purpose of this fiscal year 2023 work is to investigate the feasibility and challenges of performing such calculations within BlueCRAB so as to generate the data that MELCOR relies on. To demonstrate the methodology, the proposed workflow was applied to a prototypical heat pipe-cooled microreactor model. To predict isotopic concentrations (taking into account the ac- cumulation of fission products during operation), the necessary microscopic cross sections were generated via OpenMC and tabulated with respect to temperature and burnup. Next, a recently developed capability in Griffin (the reactor physics application in BlueCRAB) was used to convert the OpenMC output format into the ISOXML format used by Griffin. A multiphysics microscopic depletion calculation that involved performing a coupled full-core, heterogeneous neutron trans- port and thermal calculation at each depletion step was conducted to deplete the core to end of life (EOL) conditions so as to provide both isotopics and the initial condition for the transient calculation. Following a brief null-transient to verify that the initial condition had been properly restarted and was indeed in thermal equilibrium, a heat pipe failure transient was simulated. Thus, the entire workflow of using BlueCRAB to generate MELCOR inputs, from cross-section generation to producing isotopic inventory and power/temperature evolution profiles during transients, is demonstrated. This report also details the identified gaps in the workflow and how they were (for the most part) addressed. Future work should focus on directly including MEL- COR into the workflow by performing a MELCOR calculation using the BlueCRAB-generated input data. In addition, the heat pipe reactor design should be improved so as to reflect more prototypical burnup characteristics at EOL.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Task III: Development of an Effective Computational Methodology for Body Force Representation of High-speed Rotor 37

A framework for an effective computational methodology for characterizing the stability and the impact of distortion in high-speed multi-stage compressor is being developed. The methodology consists of using a few isolated-blade row Navier-Stokes solutions for each blade row to construct a body force database. The purpose of the body force database is to replace each blade row in a multi-stage compressor by a body force distribution to produce same pressure rise and flow turning. To do this, each body force database is generated in such a way that it can respond to the changes in local flow conditions. Once the database is generated, no hrther Navier-Stokes computations are necessary. The process is repeated for every blade row in the multi-stage compressor. The body forces are then embedded as source terms in an Euler solver. The method is developed to have the capability to compute the performance in a flow that has radial as well as circumferential non-uniformity with a length scale larger than a blade pitch; thus it can potentially be used to characterize the stability of a compressor under design. It is these two latter features as well as the accompanying procedure to obtain the body force representation that distinguish the present methodology from the streamline curvature method. The overall computational procedures have been developed. A dimensional analysis was carried out to determine the local flow conditions for parameterizing the magnitudes of the local body force representation of blade rows. An Euler solver was modified to embed the body forces as source terms. The results from the dimensional analysis show that the body forces can be parameterized in terms of the two relative flow angles, the relative Mach number, and the Reynolds number. For flow in a high-speed transonic blade row, they can be parameterized in terms of the local relative Mach number alone.

Tan, Choon-Sooi↗

Application of a computational model for vortex generators in subsonic internal flows

A model for the analysis of vortex generators in a fully viscous subsonic internal flow is evaluated. A vorticity source term is used in a modified form of the Parabolized Navier-Stokes equations to model the shed vortex. Computed results are compared with idealized flow vortex paths, and with experimental data for vortex generators embedded in a thick turbulent boundary layer. The analysis is also compared with experimental data for a separated diffusing S-duct and for a diffusing S-duct with vortex generators. Quantitative comparisons are shown for the latter three cases. Emphasis is placed on verifying the ability of the model to predict global distortions in the flow field.

Kunik, W. G.↗

Application of a computational model for vortex generators in subsonic internal flows

A model for the analysis of vortex generators in a fully viscous subsonic internal flow is evaluated. A vorticity source term is used in a modified form of the Parabolized Navier-Stokes equations to model the shed vortex. Computed results are compared with idealized flow vortex paths, and with experimental data for vortex generators embedded in a thick turbulent boundary layer. The analysis is also compared with experimental data for a separated diffusing S-duct and for a diffusing S-duct with vortex generators. Quantitative comparisons are shown for the latter three cases. Emphasis is placed on verifying the ability of the model to predict global distortions in the flow field.

Kunik, W. G.↗

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↗

Toward a high-fidelity tritium transport modeling for retention and permeation experiments

We report Tritium Migration Analysis Program version 8 (TMAP8), the latest version of TMAP, was developed within the framework of the Multiphysics Object-Oriented Simulation Environment (MOOSE). Created at Idaho National Laboratory (INL), MOOSE is an open-source, dimension-agnostic, fully coupled, and fully implicit multiphysics platform featuring massively parallel computation capabilities. Using TMAP8, tritium transport in a divertor monoblock was analyzed to elucidate the effects of pulsed operation (up to fifty 1,600 s plasma discharge and cool-down cycles) on the tritium in-vessel inventory source term and ex-vessel release term (i.e., tritium retention and permeation) for safety analysis. With its built-in Message Passing Interface capability, TMAP8 can, in under 2 h, simulate tritium transport in three different layered materials (i.e., tungsten, copper, and copper-chromium-zirconium alloy) in 2D geometry, using a single device/computer with 10 cores. The MOOSE-based TMAP8 code can leverage other MOOSE tools developed under the Nuclear Energy Advanced Modeling and Simulation program to perform tritium and thermal transport in complex geometries and multiphysics environments. And via its massively parallel computation, MOOSE will enable the fusion pilot plant designers to conduct high-fidelity multiphysics modeling for the design of the divertor and blanket systems as well as for the safety analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Study, optimization, and design of a laser heat engine

Laser heat engine concepts, proposed for satellite applications, are analyzed to determine which engine concept best meets the requirements of high efficiency (50 percent or better), continuous operation in space using near-term technology. The analysis of laser heat engines includes the thermodynamic cycles, engine design, laser power sources, collector/concentrator optics, receiving windows, absorbers, working fluids, electricity generation, and heat rejection. Specific engine concepts, optimized according to thermal efficiency, are rated by their technological availability and scaling to higher powers. A near-term experimental demonstration of the laser heat engine concept appears feasible utilizing an Otto cycle powered by CO2 laser radiation coupled into the engine through a diamond window. Higher cycle temperatures, higher efficiencies, and scalability to larger sizes appear to be achievable from a laser heat engine design based on the Brayton cycle and powered by a CO laser.

Taussig, R. T.↗

Sensitivity analysis of small modular reactor mPower with ATHLET

Small modular reactor (SMR) possesses the characteristics of less radioactive source term, simpler structure, modular design and features like passive safety as well as low risk of proliferation. Developed and designed by Babcock and Wilcox, the mPower reactor is one of pronouncing small modular reactors that adopts the integrated design of the pressurized water reactor, steam generator and a built-in internal control rod drive mechanism. In this paper, the code ATHLET (Analysis of fossil-hydraulics of Leaks and Transients) is used to model the mPower reactor. At the same time, the sensitivity analysis of the core structure and geometric parameters of the pipeline loop, as well as the steady-state simulation and post-exothermic simulation are carried out depending on the given scenarios. In addition, the sensitivity analysis of pressure, temperature and mass flow in the secondary loop is performed. The results of this study will mitigate some of the frequent misunderstanding in the field and will help to better understand the design of SMR and provide basic research data for further research. However, there are limitations in the current research. Considering that the mPower is still in the conceptual design, the established model is only based on the limited public data while detailed parameters are still in myth, which leads to certain level of uncertainties in the calculation. Besides, the research is mainly on the influence of sensitivity analysis under steady state conditions. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗