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 235 records · Page 13

Nickel-coated Aluminum Particles: A Promising Fuel for Mars Missions

Combustion of metals in carbon dioxide is a promising source of energy for propulsion on Mars. This approach is based on the ability of some metals (e.g. Mg, Al) to burn in CO2 atmosphere and suggests use of the Martian carbon dioxide as an oxidizer in jet or rocket engines. Analysis shows that CO2/metal propulsion will reduce significantly the mass of propellant transported from Earth for long-range mobility on Mars and sample return missions. Recent calculations for the near-term missions indicate that a 200-kg ballistic hopper with CO2/metal rocket engines and a CO2 acquisition unit can perform 10-15 flights on Mars with the total range of 10-15 km, i.e. fulfill the exploration program typically assigned for a rover. Magnesium is currently recognized as a candidate fuel for such engines owing to easy ignition and fast burning in CO2. Aluminum may be more advantageous if a method for reducing its ignition temperature is found. Coating it by nickel is one such method. It is known that a thin nickel layer of nickel on the surface of aluminum particles can prevent their agglomeration and simultaneously facilitate their ignition, thus increasing the efficiency of aluminized propellants.

Shafirovich, Evgeny↗

Nickel-Coated Aluminum Particles: A Promising Fuel for Mars Missions

Combustion of metals in carbon dioxide is a promising source of energy for propulsion on Mars. This approach is based on the ability of some metals (e.g. Mg, Al) to burn in CO2 atmosphere and suggests use of the Martian carbon dioxide as an oxidizer in jet or rocket engines. Analysis shows that CO2/metal propulsion will reduce significantly the mass of propellant transported from Earth for long-range mobility on Mars and sample return missions. Recent calculations for the near-term missions indicate that a 200-kg ballistic hopper with CO2/metal rocket engines and a CO2 acquisition unit can perform 10-15 flights on Mars with the total range of 10-15 km, i.e. fulfill the exploration program typically assigned for a rover. Magnesium is currently recognized as a candidate fuel for such engines owing to easy ignition and fast burning in CO2. Aluminum may be more advantageous if a method for reducing its ignition temperature is found. Coating it by nickel is one such method. It is known that a thin nickel layer of nickel on the surface of aluminum particles can prevent their agglomeration and simultaneously facilitate their ignition, thus increasing the efficiency of aluminized propellants. Combustion of single Ni-coated Al particles in different gas environments (O2, CO2, air) was studied using electrodynamic levitation and laser ignition. It was shown that the combustion mechanisms depend on the ambient atmosphere. Combustion in CO2 is characterized by the smaller size and lower brightness of flame than in O2, and by phenomena such as micro-flashes and fragment ejection. The size and brightness of flame gradually decrease as the particle burns.

Shafirovich, Evgeny↗

Interactions Between Mineral Dust, Climate, and Ocean Ecosystems

Over the past decade, technological improvements in the chemical and physical characterization of dust have provided insights into a number of phenomena that were previously unknown or poorly understood. In addition, models are now incorporating a wider range of physical processes, which will allow us to better quantify the climatic and ecological impacts of dust. For example, some models include the effect of dust on oceanic photosynthesis and thus on atmospheric CO 2 (Friedlingstein et al. 2006). The impact of long-range dust transport, with its multiple forcings and feedbacks, is a relatively new and complex area of research, where input from several disciplines is needed. So far, many of these effects have only been parameterized in models in very simple terms. For example, the representation of dust sources remains a major uncertainty in dust modeling and estimates of the global mass of airborne dust. This is a problem where Earth scientists could make an important contribution, by working with climate scientists to determine the type of environments in which easily erodible soil particles might have accumulated over time. Geologists could also help to identify the predominant mineralogical composition of dust sources, which is crucial for calculating the radiative and chemical effects of dust but is currently known for only a few regions. Understanding how climate and geological processes control source extent and characterizing the mineral content of airborne dust are two of the fascinating challenges in future dust research.

Gasso, Santiago↗

Fast Plasma Investigation for MMS: Simulation of the Burst Triggering System

The Magnetospheric Multiscale (MMS) mission will study small-scale reconnection structures and their rapid motions from closely spaced platforms using instruments capable of high angular, energy, and time resolution measurements. To meet these requirements, the Fast Plasma Instrument (FPI) consists of eight (8) identical half top-hat electron sensors and eight (8) identical ion sensors and an Instrument Data Processing Unit (IDPU). The sensors (electron or ion) are grouped into pairs whose 6 degree x 180 degree fields-of-view (FOV) are set 90 degrees apart. Each sensor is equipped with electrostatic aperture steering to allow the sensor to scan a 45 degree x 180 degree fan about the its nominal viewing (0 deflection) direction. Each pair of sensors, known as the Dual Electron Spectrometer (DES) and the Dual Ion Spectrometer (DIS), occupies a quadrant on the MMS spacecraft and the combination of the eight electron/ion sensors, employing aperture steering, image the full-sky every 30-ms (electrons) and 150-ms (ions), respectively. To probe the diffusion regions of reconnection, the highest temporal/spatial resolution mode of FPI results in the DES complement of a given spacecraft generating 6.5-Mb (raised dot) per second of electron data while the DIS generates 1.1-Mb (raised dot) per second of ion data yielding an FPI total data rate of 6.6-Mb (raised dot) per second. The FPI electron/ion data is collected by the IDPU then transmitted to the Central Data Instrument Processor (CIDP) on the spacecraft for science interest ranking. Only data sequences that contain the greatest amount of temporal/spatial structure will be intelligently down-linked by the spacecraft. This requires a data ranking process known as the burst trigger system. The burst trigger system uses pseudo physical quantities to approximate the local plasma environments. As each pseudo quantity will have a different value, a set of two scaling factors is employed for each pseudo term. These pseudo quantities are then combined at the instrument, spacecraft, and observatory level leading to a final ranking of data based on expected scientific interest. Here, we present simulations of the fixed point burst trigger system for the FPI. A variety of data sets based on previous mission data as well as analytical formulations are tested. Comparisons of floating point calculations versus the fixed point hardware simulation are shown. Analysis of the potential sources of error from overflows, quantization, etc. are examined and mitigation methods are presented. Finally a series of calibration curves are presented, showing the expected error in pseudo quantities based solely on the scale parameters chosen and the expected data range. We conclude with a presentation of the current base-lined FPI burst trigger approach.

Barrie, A. C.↗

On the effects of the ocean on atmospheric CFC-11 lifetimes and emissions

The ocean is a reservoir for CFC-11, a major ozone-depleting chemical. Anthropogenic production of CFC-11 dramatically decreased in the 1990s under the Montreal Protocol, which stipulated a global phase out of production by 2010. However, studies raise questions about current overall emission levels and indicate unexpected increases of CFC-11 emissions of about 10 Gg/yr after 2013 (based upon measured atmospheric concentrations and an assumed atmospheric lifetime). These findings heighten the need to understand processes that could affect the CFC-11 lifetime, including ocean fluxes. We evaluate how ocean uptake and release through 2300 affects CFC-11 lifetimes, emission estimates, and the long-term return of CFC-11 from the ocean reservoir. We show that ocean uptake yields a shorter total lifetime and larger inferred emission of atmospheric CFC-11 from 1930 to 2075 compared to estimates using only atmospheric processes. Ocean flux changes over time result in small but not completely negligible effects on the calculated unexpected emissions change (decreasing it by 0.4 ± 0.3 Gg/yr). Moreover, it is expected that the ocean will eventually become a source of CFC-11, increasing its total lifetime thereafter. Ocean outgassing should produce detectable increases in global atmospheric CFC-11 abundances by the mid-2100s, with emission of around 0.5 Gg/yr; this should not be confused with illicit production at that time. An illustrative model projection suggests that climate change is expected to make the ocean a weaker reservoir for CFC-11, advancing the detectable change in the global atmospheric mixing ratio by about 5 yr.

Anthropogenic production of CFC-11↗

Forecasting of Clear Air Turbulence using a Diagnostic Richardson Number Tendency formulation

The results of several case studies of Clear Air Turbulence (CAT) using the Diagnostic Richardson Number Tendency (DRT) formulation are highlighted. The performance of this technique in resolving regions of documented CAT encounters is encouraging. Its operational adaptability seems particularly attractive in that the input data can be supplied by the currently operational Rawinsonde system. A CAT index is calculated deterministically, sensing synoptic-scale changes in static stability and vertical wind shear conducive for supporting meso-scale CAT layers. This index reveals volumes of the troposphere which act as 'source regions' of Kelvin-Helmholtz instabilities. The results suggest that these regions are particularly efficient with respect to the synoptic/meso-scale energetical coupling necessary for supporting significantly turbulent layers. The output highlights specific regions of the atmosphere which can be interpreted operationally in terms of CAT-encounter probabilities.

Keller, J. L.↗

New Data Source for Studying and Modelling the Topside Ionosphere

The existing uncertainties about density profiles in the topside ionosphere, i.e., in the height regime from hmF2 to approx. 2000 km, requires the search for new data sources. Millions of ionograms had been recorded by the ISIS and Alouette satellites in the sixties and seventies, that never were analyzed in terms of electron density profiles. In recent years an effort started to digitize the analog recordings to prepare the ionograms for computerized analysis. This paper shows how the digital ionograms are processed and the electron density profiles (from satellite orbit altitude, 1400 km for ISIS-2, down to the F peak) are calculated. The most difficult part of the task is the automatic scaling of the echo traces in the ISIS ionograms. Unlike the ionograms from modern ionosondes, the ISIS ionograms do not identify the wave polarization of the different echo traces, so physical logic must be applied to identify the ordinary ()) and extraordinary (X) traces, and this is not always successful. Characteristic resonance features seen in the topside ionograms occur at the gyro and plasma frequencies. An elaborate scheme was developed to identify these resonance frequencies in order to determine the local plasma and gyrofrequencies. This information helps in the identification of the O and X traces, and it provides the starting density of the electron density profile. The inversion of the echo traces into electron density profiles uses the same modified Chebyshev polynomial fitting technique that is successfully applied in the ground-based Digisonde network. The automatic topside ionogram scaler with true height algorithm TOPIST is successfully scaling approx. 70% of the ionograms. An 'editing process' is available to manually scale the more difficult ionograms. The home page for the ISIS project is at http://nssdc.gsfc.nasa.gov/space/isis/isis-status.html. It provides access to as of January 2001, 3000,000 digitized ISIS ionogram data and to related software. A search page lets users select data location, time, and a host of other search criteria. The automated processing of the ISIS ionograms will begin later this year and the electron density profiles will be made available from the project home page. The ISIS data restoration efforts are supported through NASA's Applied Systems and Information Research Program.

Huang, Xue-Qin↗

SBUV version 8.6 Retrieval Algorithm: Error Analysis and Validation Technique

SBUV version 8.6 algorithm was used to reprocess data from the Back Scattered Ultra Violet (BUV), the Solar Back Scattered Ultra Violet (SBUV) and a number of SBUV/2 instruments, which 'span a 41-year period from 1970 to 2011 (except a 5-year gap in the 1970s)[see Bhartia et al, 2012]. In the new version Daumont et al. [1992] ozone cross section were used, and new ozone [McPeters et ai, 2007] and cloud climatologies Doiner and Bhartia, 1995] were implemented. The algorithm uses the Optimum Estimation technique [Rodgers, 2000] to retrieve ozone profiles as ozone layer (partial column, DU) on 21 pressure layers. The corresponding total ozone values are calculated by summing ozone columns at individual layers. The algorithm is optimized to accurately retrieve monthly zonal mean (mzm) profiles rather than an individual profile, since it uses monthly zonal mean ozone climatology as the A Priori. Thus, the SBUV version 8.6 ozone dataset is better suited for long-term trend analysis and monitoring ozone changes rather than for studying short-term ozone variability. Here we discuss some characteristics of the SBUV algorithm and sources of error in the SBUV profile and total ozone retrievals. For the first time the Averaging Kernels, smoothing errors and weighting functions (or Jacobians) are included in the SBUV metadata. The Averaging Kernels (AK) represent the sensitivity of the retrieved profile to the true state and contain valuable information about the retrieval algorithm, such as Vertical Resolution, Degrees of Freedom for Signals (DFS) and Retrieval Efficiency [Rodgers, 2000]. Analysis of AK for mzm ozone profiles shows that the total number of DFS for ozone profiles varies from 4.4 to 5.5 out of 6-9 wavelengths used for retrieval. The number of wavelengths in turn depends on solar zenith angles. Between 25 and 0.5 hPa, where SBUV vertical resolution is the highest, DFS for individual layers are about 0.5.

Kramarova, N. A.↗

Propagation and Breaking at High Altitudes of Gravity Waves Excited by Tropospheric Forcing

An anelastic approximation is used with a time-variable coordinate transformation to formulate a two-dimensional numerical model that describes the evolution of gravity waves. The model is solved using a semi-Lagrangian method with monotone (nonoscillatory) interpolation of all advected fields. The time-variable transformation is used to generate disturbances at the lower boundary that approximate the effect of a traveling line of thunderstorms (a squall line) or of flow over a broad topographic obstacle. The vertical propagation and breaking of the gravity wave field (under conditions typical of summer solstice) is illustrated for each of these cases. It is shown that the wave field at high altitudes is dominated by a single horizontal wavelength; which is not always related simply to the horizontal dimension of the source. The morphology of wave breaking depends on the horizontal wavelength; for sufficiently short waves, breaking involves roughly one half of the wavelength. In common with other studies, it is found that the breaking waves undergo "self-acceleration," such that the zonal-mean intrinsic frequency remains approximately constant in spite of large changes in the background wind. It is also shown that many of the features obtained in the calculations can be understood in terms of linear wave theory. In particular, linear theory provides insights into the wavelength of the waves that break at high altitudes, the onset and evolution of breaking. the horizontal extent of the breaking region and its position relative to the forcing, and the minimum and maximum altitudes where breaking occurs. Wave breaking ceases at the altitude where the background dissipation rate (which in our model is a proxy for molecular diffusion) becomes greater than the rate of dissipation due to wave breaking, This altitude, in effect, the model turbopause, is shown to depend on a relatively small number of parameters that characterize the waves and the background state.

Prusa, Joseph M.↗

Theoretical quasar emission-line profiles. II - Effects of source function vibration and total thickness

By calculating the actual limb brightening functions for the emission-line clouds, rather than using the approximation that all observed radiation originates from the same part of the growth curve, the present consideration of radiation transfer effects in the broad emission line regions of quasars and Seyfert 1 galaxies allows the study of clouds whose total optical thickness is intermediate between the limiting cases of an optically thin Doppler core and Voigt wing emissions. The ability to predict minor differences in observed line profiles in a given object in terms of radiative transfer effects leads to the proposal of a test which may be able to distinguish between radial and orbital geometries. The test involves the difference in width between emission lines from the same family having different optical depth scales.

Puetter, R. C.↗

Neutrino tomography - Tevatron mapping versus the neutrino sky

The feasibility of neutrino tomography of the earth's interior is discussed, taking the 80-GeV W-boson mass determined by Arnison (1983) and Banner (1983) into account. The opacity of earth zones is calculated on the basis of the preliminary reference earth model of Dziewonski and Anderson (1981), and the results are presented in tables and graphs. Proposed tomography schemes are evaluated in terms of the well-posedness of the inverse-Radon-transform problems involved, the neutrino generators and detectors required, and practical and economic factors. The ill-posed schemes are shown to be infeasible; the well-posed schemes (using Tevatrons or the neutrino sky as sources) are considered feasible but impractical.

Wilson, T. L.↗

Causal compensated perturbations in cosmology

A theoretical framework is developed to calculate linear perturbations in the gravitational and matter fields which arise causally in response to the presence of stiff matter sources in a FRW cosmology. It is shown that, in order to satisfy energy and momentum conservation, the gravitational fields of the source must be compensated by perturbations in the matter and gravitational fields, and the role of such compensation in containing the initial inhomogeneities in their subsequent evolution is discussed. A complete formal solution is derived in terms of Green functions for the perturbations produced by an arbitrary source in a flat universe containing cold dark matter. Approximate Green function solutions are derived for the late-time density perturbations and late-time gravitational waves in a universe containing a radiation fluid. A cosmological energy-momentum pseudotensor is defined to clarify the nature of energy and momentum conservation in the expanding universe.

Veeraraghavan, Shoba↗

Comparison of theories for gravity wave induced fluctuations in airglow emissions

A comparison is undertaken of theories for the gravity wave induced fluctuations in the intensity of airglow emissions and the associated temperature of the source region. The comparison is made in terms of Krassovsky's ratio eta(sub E) for a vertically extended emission region (eta(sub E) is the ratio of the vertically integrated normalized intensity perturbation to the vertically integrated intensity-weighted temperature perturbation). It is shown that the formulas for eta(sub E) in the works by Tarasick and Hines (1990) and Schubert et al. (1991) are in agreement for the case of an inviscid atmosphere. The calculation of eta(sub E) using the theory of Tarasick and Hines (1990) requires determination of their function chi; we show that chi is simply related to the 'single-level' Krassovsky's ratio eta of Schubert et al. (1991). The general relationship between chi and eta is applied to a simple chemical-dynamical model of the O2 atmospheric airglow and the altitude dependence of these quantities is evaluated for nonsteady state chemistry. Though the Tarasick and Hines (1990) formula for eta(sub E) does not explicitly depend on the scale heights of the minor constituents involved in airglow chemistry, eta(sub E) implicitly depends upon these scale heights through its dependences on chemical production and loss contained in chi. We demonstrate this dependence of eta(sub E) for the OH nightglow on atomic oxygen scale height by direct numerical evaluation of eta(sub E) in this case the dependence originates in the chemical production of perturbed ozone.

Walterscheid, R. L.↗

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↗

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↗