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 469 records · Page 26

On the theoretical limitations in estimating thickness of a plate-like structure from a full-field single-tone response Lamb wave measurement

A persistent question in wavenumber analysis in the estimation of thickness from a steady-state wave field is identifying the theoretical sensitivity of the system. Here, a well-known trade-off between spatial frequency/wavenumber resolution and thickness resolution exists. The current work presents a calculation of the Cramer-Rao Lower bound (CRLB), specifically as applied to thickness estimates, for a 2-dimensional multi-mode waveform in Gaussian noise. Cases of near-field and far-field excitation are considered, and transducer position with respect to the scan area is also varied. Additionally, we consider the CRLB in a plate with multiple sources, simulated as sources placed on the boundary of the plate. We conclude by presenting the CRLB values in terms of frequency for various thicknesses, and by presenting optimal excitation frequencies for a nominal thickness, based on the CRLB.

47 OTHER INSTRUMENTATION↗

Nuclear Technology Selection and Preliminary Integration at St. Charles and Carrollton Chemical Plants

Nuclear energy is identified as a valuable candidate for decarbonizing industrial facilities as a non-emitting baseload energy source capable of dispatching energy as needed for an integrated facility. As industrial facilities are unique, analyzing what nuclear technology is the best candidate should be evaluated on a case-by-case basis. For two identified facilities of the Dow Chemical Company, an ideal mix of nuclear energy resources is calculated, and integration strategies are put forward. This research was funded through a Gateway for Accelerated Innovation in Nuclear (GAIN) Nuclear Energy Voucher, specifically seeking near-term engineering solutions for the technical issue of introducing nuclear generated energy into the Dow chemical facilities. The research is split into four main phases: background on reactors and chemical facilities, preliminary engineering analysis determining how reactor technologies may be integrated with specific facilities, site characterization and discrete hazards analysis, and finally an optimization of reactor-facility selection through the use of Framework for Optimization of ResourCes and Economics (FORCE) tools.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

LU-SGS implicit algorithm for three-dimensional incompressible Navier-Stokes equations with source term

A numerical method is developed for solving the incompressible Navier-Stokes equations using the concept of pseudocompressibility. A lower-upper symmetric-Gauss-Seidel implicit scheme is developed for three-dimensional incompressible viscous flow computations. The present algorithm offers additional advantages when solving the flow equations with source terms. Complete vectorizability of the algorithm on oblique planes of sweep in three-dimensions is accomplished in a new flow solver, INS3D-LU code. Spatial differencing is a second-order accurate semi-discrete finite-volume method augmented by a third-order accurate numerical dissipation model which is based on spectral-radii. Comparison of numerical solutions for a curved duct with experimental data shows good agreement. The method is applied to calculate the inducer flow of the Space Shuttle Main Engine turbopump.

Yoon, Seokkwan↗

Comparison of low-latitude ion and neutral zonal drifts using DE 2 data

We have used data from the ion drift meter and the wind and temperature spectrometer on the DE 2 spacecraft to make statistical comparisons of the zonal ion and neutral drifts at dip latitudes (DLAT) in the +/- 35 deg range over all local times. Fourier analysis indicates that the superrotation and the diurnal components of both flows are strongly peaked at the dip equator, with the superrotation term becoming negative for the absolute value DLAT greater than or equal 20 deg. One interesting feature is the presence of a period (2200-0500 solar local time) in the 300-400 km altitude region near the dip equator where the ion drift is more strongly eastward than the neutral flow. This would seem to indicate the presence of an electric field source of greater strength than the F region dynamo elsewhere along the geomagnetic field line. Model calculations indicate that a possible mechanism for this source lies in the vertical shear in the zonal neutral wind in the 100-200 km altitude region.

Coley, W. R.↗

Deterministic Stress Modeling of Hot Gas Segregation in a Turbine

Simulation of unsteady viscous turbomachinery flowfields is presently impractical as a design tool due to the long run times required. Designers rely predominantly on steady-state simulations, but these simulations do not account for some of the important unsteady flow physics. Unsteady flow effects can be modeled as source terms in the steady flow equations. These source terms, referred to as Lumped Deterministic Stresses (LDS), can be used to drive steady flow solution procedures to reproduce the time-average of an unsteady flow solution. The goal of this work is to investigate the feasibility of using inviscid lumped deterministic stresses to model unsteady combustion hot streak migration effects on the turbine blade tip and outer air seal heat loads using a steady computational approach. The LDS model is obtained from an unsteady inviscid calculation. The LDS model is then used with a steady viscous computation to simulate the time-averaged viscous solution. Both two-dimensional and three-dimensional applications are examined. The inviscid LDS model produces good results for the two-dimensional case and requires less than 10% of the CPU time of the unsteady viscous run. For the three-dimensional case, the LDS model does a good job of reproducing the time-averaged viscous temperature migration and separation as well as heat load on the outer air seal at a CPU cost that is 25% of that of an unsteady viscous computation.

Busby, Judy↗

Python Based Plume Dynamics Estimation Tool (PyPDET) Rapid Plume Strike Analysis for RPOD Maneuvers in Deep Space Operations

I worked as a NASA Intern during the Summer 2023 term in the DS-00 division under the supervision of my mentor, Dr. Jonathan Pitt. Our goal was to build on our previous work from 2022 to develop a plume strike estimation tool using a prescribed physics methodology and model plume impingement effects while considering the dynamics of a rendezvous, operations, proximity, and docking (RPOD) maneuver. This tool supports previously configured CFD-DSMC calculations by allowing for rapid analysis of initial designs using a low-fidelity source flow model. Engineers can then use the high-fidelity CFD-DSMC tool to consolidate results as they work towards finalizing a design. This year’s project was focused on developing a software application that other engineers would be using in their analysis. Thus, the user’s experience was considered in the development of this application. Proper documentation, testability, and modularity of the codebase was our priority. For example, the project included auto documentation procedures to start building towards a User Manual, while also including dedicated demonstration cases for more explicit communication of functionality. Also, this project included a framework for testing the source code for future developments. Additionally, care was taken to develop the code using an Object-Oriented Programming approach. Thus, allowing for a modular extensibility of functionality in anticipation of future developments. The core work of this project was developing an algorithm that would transform the visiting vehicle and associated thruster data according to the kinematics described in the jet firing history. It would then calculate the estimated plume strikes on two of the target vehicles and write data accordingly into a VTK file. Summer work is to conclude by developing and presenting a PowerPoint slide deck at the intern exit briefing on August 11 th , 2023. Once the model for simple plume strike calculations is developed and tested there are several avenues to explore to continue development of this tool. These are also discussed in this report.

Plume Impingement↗

First STAX detector installation at the National Institute for Radioelements (IRE)

The Source Term Analysis of Xenon (STAX) project has been installing stack detectors at medical isotope production facilities to measure radioxenon emissions to investigate the effect of radioxenon releases on nuclear explosion monitoring. This paper outlines the installation of the first STAX detection system at the National Institute for Radioelements (IRE) in Fleurus, Belgium which has been operating for over three years and transferring collected data to the STAX repository. Information about the equipment installed, the data flow established, and calculations for determination of radioxenon releases from the facility are presented. Further. data quality was investigated to confirm values reported by STAX automated data processing and in a comparison of collected STAX data with data collected by IRE for regulatory reporting.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Thermodynamics of the internal combustion engine

The ubiquitous automobile has played both the role of hero and villain during the course of its history. At the present time it is being cast as the villain with increasing frequency because it has been indicted as a major source of air pollution. It is just this characterization that is responsible for our renewed interest in the automobile. When we speak of the automobile as a source of air pollution, we really mean the heart of the automobile, the internal combustion engine. The prototype of the modern internal combustion engine was first successfully operated by Nicholas August Otto in 1876. In spite of its long history, the internal combustion engine is still poorly understood in terms of the details of the physical processes that take place during its operations. The internal combustion engine operates on what is known as the Otto cycle, and Otto cycle calculations have been with us for many years. The emphasis in earlier calculations was usually on the performance rather than on the pollution aspects. Performance characteristics can often be estimated reasonably accurately even with a considerable degree of simplification. We are interested in the pollution aspects of the internal combustion engine, and therefore the calculations must be performed with more care. Our extensive experience with equilibrium calculations places us in a position to carry out such careful calculations. We have just recently begun such an analysis, and its objective is the prediction of exhaust gas properties. I will describe our progress by first reviewing the Otto cycle, and then I will illustrate the effect of various assumptions on the calculated results.

Frank J. Zeleznik↗

Estimates of Power Plant NOx Emissions and Lifetimes from OMI NO2 Satellite Retrievals

Isolated power plants with well characterized emissions serve as an ideal test case of methods to estimate emissions using satellite data. In this study we evaluate the Exponentially-Modified Gaussian (EMG) method and the box model method based on mass balance for estimating known NOx emissions from satellite retrievals made by the Ozone Monitoring Instrument (OMI). We consider 29 power plants in the USA which have large NOx plumes that do not overlap with other sources and which have emissions data from the Continuous Emission Monitoring System (CEMS). This enables us to identify constraints required by the methods, such as which wind data to use and how to calculate background values. We found that the lifetimes estimated by the methods are too short to be representative of the chemical lifetime. Instead, we introduce a separate lifetime parameter to account for the discrepancy between estimates using real data and those that theory would predict. In terms of emissions, the EMG method required averages from multiple years to give accurate results, whereas the box model method gave accurate results for individual ozone seasons.

Emissions inventory↗

Issues in Absolute Spectral Radiometric Calibration: Intercomparison of Eight Sources

The application of atmospheric models to AVIRIS and other spectral imaging data to derive surface reflectance requires that the sensor output be calibrated to absolute radiance. Uncertainties in absolute calibration are to be expected, and claims of 92% accuracy have been published. Measurements of accurate surface albedos and cloud absorption to be used in radiative balance calculations depend critically on knowing the absolute spectral-radiometric response of the sensor. The Earth Observing System project is implementing a rigorous program of absolute radiometric calibration for all optical sensors. Since a number of imaging instruments that provide output in terms of absolute radiance are calibrated at different sites, it is important to determine the errors that can be expected among calibration sites. Another question exists about the errors in the absolute knowledge of the exoatmospheric spectral solar irradiance.

Goetz, Alexander F. H.↗

A projection-based analytical Jacobian framework for chemical kinetics applications

A major challenge in simulating complex combustion systems with detailed chemical kinetic models is the cost of integrating the chemical source terms, often done using stiff ODE solvers that require frequent Jacobian evaluations. Using analytically derived Jacobian matrices instead of divided-difference-based numerical Jacobian approximations can significantly reduce the associated computational cost. However, ambiguities arise in the formulation of analytical Jacobians because the chemical state of the system, or state vector, can be expressed in multiple ways, involving variables that are typically not independent from one another. Here, in this work, the consequences of those ambiguities on practical calculations are characterized in detail, and a generalized, projection-based framework is proposed as a mitigation strategy. Performances are assessed in a series of test cases involving a variety of configurations and numerical solution approaches. Results show that with proper treatment, commonly used analytical Jacobian formulations can be considered as equivalent for practical purposes, thereby alleviating concerns that the state vector chosen to express the governing equations and corresponding analytical Jacobian may significantly impact the accuracy of the simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Consequence analyses of sabotage-induced radiological releases in sodium-cooled fast microreactors

Analysis of three sodium-cooled fast microreactors (SFMs) with thermal powers of 10, 30, and 50 MWt showed that smaller reactors result in lower radiological consequences during a postulated sabotage-induced event because of their reduced core inventory. All SFMs used U-10Zr metal fuel enriched to 15 wt% high-assay low-enriched uranium and operated until their respective effective multiplication factor (k eff ) reduced to less than 1 or until the end of their operational lifespan. Sabotage scenarios were simulated at this point, when the fuel inventory within the core contains the highest-level of radioactivity. Radionuclide core inventories were calculated using the SCALE code at shutdown and 3 days post-shutdown. Dose consequence analyses were performed for three sabotage scenarios using the RASCAL tool. As microreactor developers plan for minimal on-site or complete off-site emergency response, it remains essential to evaluate their physical protection needs and potential hazards, including assessing postulated sabotage-induced events that could become more relevant. SFM licensees should identify a credible worst-case, major accident, estimate release source terms, and perform dose consequence analyses to evaluate site-specific physical protection measures. In conclusion, this recommendation supports a risk-informed, performance-based approach, aligning with applicable regulatory requirements, i.e., 10 CFR Parts 100 and 53 rulemaking in the United States.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Voyager 2 plasma ion observations in the magnetosphere of Uranus

Positive ion measurements in the magnetosphere of Uranus have been made by the Voyager 2 plasma science experiment. The paper presents an overview of the entire data set and a detailed analysis of the observations from the inner magnetosphere which complements and extends results reported elsewhere. Densities and temperatures are obtained from an analysis which incorporates details of the instrumental response. These results are then used to calculate flux tube particle and energy content to support the hypothesis that the plasma transport is controlled by a solar wind-driven magnetospheric convection system. Variations in the flux tube content suggest both a local source of plasma, produced from the neutral hydrogen corona of Uranus, and a nonlocal source, convected inwared and heated by adiabatic compression. In each case a proton composition is inferred. Sharp boundaries in the high-energy (approximately 1 keV) plasma population are interpreted in terms of the spatial extent of the magnetospheric convection, with significant shielding of the convection electric field. The convection theory is also used in a simulation of the low-energy (approximately 10 eV) ion component using the neutral hydrogen source, resulting in distribution functions which qualitatively agree with the observations.

Selesnick, Richard S.↗

The construction and use of divergence free vector expansions for incompressible fluid flow calculations

For incompressible fluids the law of mass conservation reduces to a constraint on the velocity vector, namely that it be divergence free. This constraint has long been a source of great difficulty to the numericist seeking to discretize the Navier-Stokes and Euler equations. A spectral method is discussed which overcomes this difficulty. Its efficacy is demonstrated on some simple problems. The velocity is approximated by a finite sum of divergence free vectors, each of which satisfies the same boundary conditions as the velocity. Projecting the governing equation onto the space of inviscid vector fields eliminates the pressure term and produces a set of ordinary differential equations that must be solved for the coefficents in the velocity. The pressure can then be recovered if it is needed.

Mhuiris, N. M. G.↗

Natural noise above 50 MHZ from terrestrial and extraterrestrial sources

This paper offers a brief overview of natural radio noise for frequencies above 50 MHz in terms of brightness temperature as observed from two vantage points. The first is from an Earth station located at 40 degrees north latitude and observing at elevation angles from 0 to 90 degrees with an ideal antenna. The second is a satellite in geostationary orbit communicating with the Earth. Earth station noise at VHF and UHF is dominated by galactic and solar noise. Emission from the atmosphere, gases and hydrometeors, are dominant at EHF and SHF. Radiative transfer theory is invoked in the calculation of brightness temperature from the atmosphere. The situation is not vastly different from geostationary orbit if communications is with the Earth. Emission from the land and sea, even under idealized conditions, enters significantly. Land is a much more effective emitter than sea water, but at frequencies above 30 GHz the differential becomes much less due to the increasing significance of atmospheric emission.

Smith, E. K.↗

Shedding light on the 239Pu fission source term with new high-precision experiments and advanced fission modeling

In the last decade, there has been a renaissance of fission research resulting in new high-precision experiments and advanced fission modeling. For instance, the Chi-Nu and CEA teams supplied, for the first time, the 239 Pu prompt fission neutron spectrum (PFNS) for broad ranges of incident and outgoing neutron energies. The CEA team also measured 239 Pu average prompt neutron multiplicities, ν ̄ p , with lower statistical uncertainties and a technique significantly different than the one used in the past. The NIFFTE collaboration provided 239 Pu( n ,f)/ 235 U( n ,f) cross section shape ratios with uncertainties below 1% utilizing a novel detector type. Advanced fission event generators were developed, among them CGMF, FIFRELIN, FREYA, and GEF, which calculate post-scission fission observables in a correlated manner. These new experimental data and more consistent fission models change the evaluated PFNS, ν ̄ p , and ( n ,f) cross sections, some only modestly, compared to ENDF/B-VIII.0. In turn, the individual new nuclear data distinctly change simulated effective neutron multiplication factors of fast critical assemblies, but their combined impact is small, while affecting the prediction of LLNL pulsed sphere neutron leakage spectra and reaction rates only within experimental uncertainties. Also, the parameters obtained from fitting to ν ̄ p reproduce various post-scission fission observables within the uncertainties of experimental data. This indicates that new differential experiments and consistent fission modeling reduce compensating errors present in ENDF/B-VIII.0.

fission cross section↗

Integrated Bosch Process System Models for In-Situ Oxygen and Carbon Production

In-Situ Resource Utilization (ISRU) technology is a vital component to NASA’s mission of a sustainable presence on the Moon and Mars. Local resources can be leveraged to reduce resupply frequency and mass. Elements of the Bosch process, combined with the carbothermal reduction process, can produce oxygen on the lunar surface with minimal consumables. The Bosch process can also produce oxygen on the Martian surface by using the CO 2 -rich environment. Between both systems, adsorption pump, solar thermal energy, carbon formation reactor, and water recovery subsystems are modeled and integrated to create a functional model in MATLAB software. The model is used to simulate performance of the system and reduce mass, power, and volume requirements. This integrated system model provides a tool to scale ISRU technologies for oxygen and carbon production. The MATLAB model is created by developing a system of independent subsystem models that are solved for their quasi-steady state values which can be integrated with respect to time to determine the change in current states. A flexible time stepping method is used to ensure a high level of accuracy during periods of rapid change while still making use of a simple explicit integration method. The flexible time step is calculated for each independent subsystem and the minimum value from those is used as the overall time step. A flexible time step is calculated by dividing a resolution value, or the maximum change per time step, by the variables current rate of change. The maximum value from all points in space is used for subsystem models that contain multiple values. The process is done for every variable that is being monitored in each subsystem and the global minimum is used as that iteration’s timestep. Several assumptions used in the MATLAB model for fluid flow dynamics, such as 1-D gas flow through the sorption pump, are supported by modeling in Ansys Fluent software. The Lunar oxygen production system is outlined in Fig. 1. The carbothermal reduction subsystem uses solar energy to heat a mixture of lunar regolith and carbon powder to produce carbon monoxide. To begin, the carbon monoxide feeds to the modified Bosch subsystem along with hydrogen gas. The reactants then enter the carbon formation reactor where water and carbon powder are produced. Solar thermal energy is used to add energy to the reactor, but waste heat from the carbothermal process is another potential heat source. The water is collected and electrolyzed to produce hydrogen which reenters the Bosch subsystem, and the oxygen is stored for downstream use. The carbon powder is collected and feeds back into the carbothermal subsystem. The Martian oxygen production system uses the full Bosch process and is outlined in Fig 2. A CO 2 adsorption pump thermally cycles to scrub and pressurize CO 2 from the environment. Along with an initial supply of hydrogen, the reactants enter the Reverse Water Gas Shift Reactor (RWGSR) which produces carbon monoxide and water. Carbon monoxide and unreacted hydrogen enter the carbon formation reactor to produce water and carbon powder. The water is collected from both reactors and electrolyzed to reintroduce hydrogen and store oxygen for propellant production or life support. Carbon is removed from the carbon formation reactor and stored. The adsorption pump utilizes rapid cycle temperature swings within a stack of zeolite coated surfaces. The subsystem model solves 1-D quasi-steady conservation laws of the quasi-steady form, shown in Eq. 1, for the gas stream and heat exchange liquid to predict performance parameters such as breakthrough capacity and optimum cycle time. The source term S is used to capture interactions between the fluid flows and the sorbent. A quasi-steady-state scheme is used where no time derivatives appear in the governing equations, except for those in the source terms. This results in an autonomous system, where ∂F/∂x = ƒ(F). The fluxes F are provided at the inlet, and an explicit method is used to solve for the spatial distribution of F. The heat and mass flows to the sorbent are then extracted from the source terms. These flows are numerically integrated to produce a 1-D solution for the system’s state as a function of both time and space. The body of the adsorption pump is separated into two semi-independent models: the heat exchanger fluid flow and gas flow through the zeolite coated surfaces. Both models are solved using the above-described method to find a 1-D solution as a function of space and interact only once a timestep is taken. The interaction point is the sorbent through which all heat transfer between the two models must occur. Sorbent mass adsorption is calculated using the Lagergren model, shown in Eq. 2, where the transfer coefficient, λ D , is found by solving a system of nondimensionalized equations derived by using the heat and mass transfer analogy for transport phenomena. Using Grade 544 Type 13X zeolite as the sorbent material, the equilibrium concentration, θ eq , is calculated using the k-site Langmuir isotherm and fit parameters. Additionally, the enthalpy of adsorption used in the model is computed by interpolation of available data [1]. The subsystem model was validated using the Rapid Cycle Temperature Swing Adsorption (RC-TSA) pump. The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC). The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC).

In situ Resource Utilization↗

SCALE Code System

The SCALE Code System is a widely-used modeling and simulation suite for nuclear safety analysis and design that is developed, maintained, tested, and managed by the Reactor and Nuclear Systems Division (RNSD) of Oak Ridge National Laboratory (ORNL). SCALE provides a comprehensive, verified and validated, user-friendly tool set for criticality safety, reactor and lattice physics, radiation shielding, spent fuel and radioactive source term characterization, and sensitivity and uncertainty analysis. Since 1980, regulators, licensees, and research institutions around the world have used SCALE for safety analysis and design. SCALE provides an integrated framework with dozens of computational modules including three deterministic and three Monte Carlo radiation transport solvers that are selected based on the desired solution strategy. SCALE includes current nuclear data libraries and problem-dependent processing tools for continuous-energy (CE) and multigroup (MG) neutronics and coupled neutron-gamma calculations, as well as activation, depletion, and decay calculations. SCALE includes unique capabilities for automated variance reduction for shielding calculations, as well as sensitivity and uncertainty analysis. SCALE’s graphical user interfaces assist with accurate system modeling, visualization of nuclear data, and convenient access to desired results.

97 MATHEMATICS AND COMPUTING↗