Engineering PapersSearch

SEARCH · Engineering Papers

Results for “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 55 records · Page 3

Methane on Mars: Thermodynamic Equilibrium and Photochemical Calculations

The detection of methane (CH4) in the atmosphere of Mars by Mars Express and Earth-based spectroscopy is very surprising, very puzzling, and very intriguing. On Earth, about 90% of atmospheric ozone is produced by living systems. A major question concerning methane on Mars is its origin - biological or geological. Thermodynamic equilibrium calculations indicated that methane cannot be produced by atmospheric chemical/photochemical reactions. Thermodynamic equilibrium calculations for three gases, methane, ammonia (NH3) and nitrous oxide (N2O) in the Earth s atmosphere are summarized in Table 1. The calculations indicate that these three gases should not exist in the Earth s atmosphere. Yet they do, with methane, ammonia and nitrous oxide enhanced 139, 50 and 12 orders of magnitude above their calculated thermodynamic equilibrium concentration due to the impact of life! Thermodynamic equilibrium calculations have been performed for the same three gases in the atmosphere of Mars based on the assumed composition of the Mars atmosphere shown in Table 2. The calculated thermodynamic equilibrium concentrations of the same three gases in the atmosphere of Mars is shown in Table 3. Clearly, based on thermodynamic equilibrium calculations, methane should not be present in the atmosphere of Mars, but it is in concentrations approaching 30 ppbv from three distinct regions on Mars.

Levine, J. S.

Spectral Efficiency and Tandem Performance Calculators [SWR-24-92]

The Spectral Efficiency and Tandem Performance Calculators allow users to test material types for use in tandem photovoltaic devices to better understand if a given material could work well. The Spectral Efficiency (SE) calculator takes a set of current vs voltage (I-V) and quantum efficiency (QE) data along with a spectrum input and plots the spectral efficiency of the cells represented by the datasets for the spectrum provided. The Tandem Performance calculator uses the calculated spectral efficiency of top and bottom cells and the transmission data of the top cell to calculate performance metrics for a tandem device in either a two or four terminal architecture. The SE and Tandem Performance calculators are based on the following publication: Yu, Z., Leilaeioun, M. & Holman, Z. Selecting tandem partners for silicon solar cells. Nat Energy 1, 16137 (2016).

Warren, Emily

R&D GREET Battery Carbon Footprint Calculator

The Battery Carbon Footprint (CF) Calculator was developed to help U.S. battery manufacturers meet the carbon footprint reporting requirements of the EU Battery Regulation (EU) 2023/1542. The calculator incorporates several major battery carbon footprint frameworks, including the Joint Research Centre's Rules for the Calculation of the Carbon Footprint of Electric Vehicle Batteries (CFB-EV), RECHARGE's Product Environmental Footprint Category Rules for High Specific Energy Rechargeable Batteries for Mobile Applications (PEFCR), the Catena-X Product Carbon Footprint Rulebook (CX-PCF Rules), Battery Pass's Battery Carbon Footprint: Rules for Calculating the Carbon Footprint of the "Distribution" and "End-of-Life and Recycling" Life Cycle Stages, the Global Battery Alliance's Greenhouse Gas Rulebook: Generic Rules, Version 2.1, and the Ministry of Economy, Trade and Industry's draft Carbon Footprint Calculation Method for Automotive Batteries. The tool pairs these frameworks with foreground data from Argonne's R&D GREET models and integrates user-supplied background data covering battery manufacturing and supply chain activities. By bringing multiple international methodologies together in a single platform, the calculator enables manufacturers to evaluate product carbon footprints, improve data consistency, and prepare for evolving regulatory compliance and global market reporting requirements.

Zhang, Jingyi

Equilipy: a python package for calculating phase equilibria

The CALPHAD (CALculation of PHAse Diagram) approach (Nigel Saunders & Miodownik, 1998) provides predictions for thermodynamically stable phases in multicomponent-multiphase materials across a wide range of temperatures. Consequently, the CALPHAD calculations became an essential tool in materials and process design (Luo, 2015). Such design tasks frequently require navigating a high-dimensional space due to multiple components involved in the system. This increasing complexity demands high-throughput CALPHAD calculations, especially in the rapidly evolving field of alloy design. In response to the need, we developed Equilipy an open-source Python package designed for calculating phase equilibria of multicomponent-multiphase systems. Equilipy is specifically tailored for high-throughput CALPHAD calculations, offering parallel computations across multiple processors and nodes with the given NPT input conditions namely elemental compositions (N), pressure (P), and temperature (T). Equilipy utilizes the program structure and Gibbs energy functions from the Fortran-based program, Thermochimica (Piro et al., 2013), with incorporating a new Gibbs energy minimization algorithm. This algorithm, originally developed by Capitani and Brown in 1987 (Capitani & Brown, 1987), has been revised and implemented to enhance the stability and performance of calculations. The Fortran codes are precompiled and interfaced with Python via F2PY, ensuring high computation speed. Benchmark tests shown in Figure 1 demonstrate that Equilipy’s computation speed is comparable to those of established commercial software, TC-Python and PanPython. This result highlights its efficiency and potential applications in various scientific and industrial fields.

97 MATHEMATICS AND COMPUTING

NASA-Lewis experiences with multigroup cross sections and shielding calculations

The nuclear reactor shield analysis procedures employed at NASA-Lewis are described. Emphasis is placed on the generation, use, and testing of multigroup cross section data. Although coupled neutron and gamma ray cross section sets are useful in two dimensional Sn transport calculations, much insight has been gained from examination of uncoupled calculations. These have led to experimental and analytic studies of areas deemed to be of first order importance to reactor shield calculations. A discussion is given of problems encountered in using multigroup cross sections in the resolved resonance energy range. The addition to ENDF files of calculated and/or measured neutron-energy-dependent capture gamma ray spectra for shielding calculations is questioned for the resonance region. Anomalies inherent in two dimensional Sn transport calculations which may overwhelm any cross section discrepancies are illustrated.

Lahti, G. P.

NASA-Lewis experiences with multigroup cross sections and shielding calculations.

Description of the nuclear reactor shield analysis procedures employed at NASA-Lewis Research Center. Emphasis is placed on the generation, use, and testing of multigroup cross-section data. Although coupled neutron and gamma ray cross-section sets are useful in two-dimensional Sn transport calculations, much insight has been gained from examination of uncoupled calculations. These have led to experimental and analytic studies of areas deemed to be of first-order importance to reactor shield calculations. A discussion is given of problems encountered in using multigroup cross sections in the resolved resonance energy range. The addition to ENDF files of calculated and/or measured neutron-energy-dependent capture gamma ray spectra for shielding calculations is questioned for the resonance region. Anomalies inherent in two-dimensional Sn transport calculations which may overwhelm any cross-section discrepancies are illustrated.

Lahti, G. P.

WETAIR: A computer code for calculating thermodynamic and transport properties of air-water mixtures

A computer program subroutine, WETAIR, was developed to calculate the thermodynamic and transport properties of air water mixtures. It determines the thermodynamic state from assigned values of temperature and density, pressure and density, temperature and pressure, pressure and entropy, or pressure and enthalpy. The WETAIR calculates the properties of dry air and water (steam) by interpolating to obtain values from property tables. Then it uses simple mixing laws to calculate the properties of air water mixtures. Properties of mixtures with water contents below 40 percent (by mass) can be calculated at temperatures from 273.2 to 1497 K and pressures to 450 MN/sq m. Dry air properties can be calculated at temperatures as low as 150 K. Water properties can be calculated at temperatures to 1747 K and pressures to 100 MN/sq m. The WETAIR is available in both SFTRAN and FORTRAN.

Fessler, T. E.

A discussion of theoretical ionization equilibrium calculations based on solar flare X-ray spectra

Several sets of ionization equilibrium calculations exist for use in interpreting X-ray and EUV spectra of astrophysical plasmas. In particular, the calculations of Jordan (1969, 1970), Jacobs et al. (1977, 1978) and Summers (1974) are well known. The temperatures of maximum fractional abundance calculated by Summers for the more highly ionized and heavier elements such as iron are about a factor of two higher than the temperatures calculated by Jordan and Jacobs et al. By the use of recently obtained X-ray spectra of solar flares, it is shown that the temperatures calculated by Summers (1974) for iron are incorrect. The temperatures calculated by Jordan or Jacobs et al. should be used until further improvements become available.

Feldman, U.

Calculation of the 3-D viscous flow at the endwall leading edge region of an axial annular turbine cascade

A three-dimensional viscous computer code (VANS/MD) was employed to calculate the turbulent flow field at the end wall leading edge region of a 20 inch axial annular turbine cascade. The initial boundary layer roll-up and formation of the end wall vortices were computed at the vane leading edge. The calculated flow field was found to be periodic with a frequency of approximately 1600 Hz. The calculated size of the separation region for the hub endwall vortex compared favorably with measured endwall oil traces. In an effort to determine the effects of the turbulence model on the calculated unsteadiness, a laminar calculation was made. The periodic nature of the calculated flow field persisted with the frequency essentially unchanged.

Walitt, L.

Calculations of cosmogenic nuclide production rates in the Earth's atmosphere and their inventories

The production rates of cosmogenic isotopes in the Earth's atmosphere and their resulting terrestrial abundances have been calculated, taking into account both geomagnetic and solar-modulatory effects. The local interstellar flux was assumed to be that of Garcia-Munoz, et al. Solar modulation was accounted for using the heliocentric potential model and expressed in terms of the Deep River neutron monitor count rates. The geomagnetic field was presented by vertical cutoffs calculated by Shea and Smart and the non-vertical cutoffs calculated using ANGRI. The local interstellar particle flux was first modulated using the heliocentric potential field. The modulated cosmic-ray fluxes reaching the earth's orbit then interacted with the geomagnetic field as though it were a high-pass filter. The interaction of the cosmic radiation with the Earth's atmosphere was calculated utilizing the Bolztmann transport equation. Spallation cross sections for isotope production were calculated using the formalism of Silberberg and Tsao and other cross sections were taken from standard sources. Inventories were calculated by accounting from the variation in solar modulation and geomagnetic field strength with time. Results for many isotope, including C-14, Be-7 and Be-10 are in generally good agreement with existing data. The C-14 inventory, for instance, amounts to 1.75/sq cm(e)/s, in excellent agreement with direct estimates.

Obrien, K.

An explicit finite-volume time-marching procedure for turbulent flow calculations

A method was developed which calculates two-dimensional, transonic, viscous flow in ducts. The finite-volume, time-marching formulation is used to obtain steady flow solutions of the Reynolds-averaged form of the Navier-Stokes equations. The entire calculation is performed in the physical domain. Control volumes are chosen so that smoothing of flow properties, typically required for stability, is not required. Different time steps are used in the different governing equations. A new pressure interpolation scheme is introduced which improves the shock capturing ability of the method. A multi-volume method for pressure changes in the boundary layer allows calculations which use very long and thin control volumes (length/height - 1000). The method is compared with two test cases. Essentially incompressible turbulent boundary layer flow in an adverse pressure gradient is calculated and the computed distributions of mean velocity and shear are in good agreement with the measurements. Transonic viscous flow in a converging diverging nozzle is calculated; the Mach number upstream of the shock is approximately 1.25. The agreement between the calculated and measured shock strength and total pressure losses is good.

Nicholson, Stephen

An explicit finite-volume time-marching procedure for turbulent flow calculations

A method was developed which calculates two-dimensional, transonic, viscous flow in ducts. The finite-volume, time-marching formulation is used to obtain steady flow solutions of the Reynolds-averaged form of the Navier-Stokes equations. The entire calculation is performed in the physical domain. Control volumes are chosen so that smoothing of flow properties, typically required for stability, is not required. Different time steps are used in the different governing equations. A new pressure interpolation scheme is introduced which improves the shock capturing ability of the method. A multi-volume method for pressure changes in the boundary layer allows calculations which use very long and thin control volumes (length/height - 1000). The method is compared with two test cases. Essentially incompressible turbulent boundary layer flow in an adverse pressure gradient is calculated and the computed distributions of mean velocity and shear are in good agreement with the measurements. Transonic viscous flow in a converging diverging nozzle is calculated; the Mach number upstream of the shock is approximately 1.25. The agreement between the calculated and measured shock strength and total pressure losses is good.

Nicholson, Stephen

Measurement effects on the calculation of in-flight thrust for an F404 turbofan engine

A study has been performed that investigates parameter measurement effects on calculated in-flight thrust for the General Electric F404-GE-400 afterburning turbofan engine which powered the X-29A forward-swept wing research aircraft. Net-thrust uncertainty and influence coefficients were calculated and are presented. Six flight conditions were analyzed at five engine power settings each. Results were obtained using the mass flow-temperature and area-pressure thrust calculation methods, both based on the commonly used gas generator technique. Thrust uncertainty was determined using a common procedure based on the use of measurement uncertainty and influence coefficients. The effects of data nonlinearity on the uncertainty calculation procedure were studied and results are presented. The advantages and disadvantages of using this particular uncertainty procedure are discussed. A brief description of the thrust-calculation technique along with the uncertainty calculation procedure is included.

Conners, Timothy R.

Measurement effects on the calculation of in-flight thrust for an F404 turbofan engine

A study was performed that investigates parameter measurement effects on calculated in-flight thrust for the General Electric F404-GE-400 afterburning turbofan engine which powered the X-29A forward-swept wing research aircraft. Net-thrust uncertainty and influence coefficients were calculated and are presented. Six flight conditions were analyzed at five engine power settings each. Results were obtained using the mass flow-temperature and area-pressure thrust calculation methods, both based on the commonly used gas generator technique. Thrust uncertainty was determined using a common procedure based on the use of measurement uncertainty and influence coefficients. The effects of data nonlinearity on the uncertainty calculation procedure were studied and results are presented. The advantages and disadvantages of using this particular uncertainty procedure are discussed. A brief description of the thrust-calculation technique along with the uncertainty calculation procedure is included.

Conners, Timothy R.

Evaluation of various thrust calculation techniques on an F404 engine

In support of performance testing of the X-29A aircraft at the NASA-Ames, various thrust calculation techniques were developed and evaluated for use on the F404-GE-400 engine. The engine was thrust calibrated at NASA-Lewis. Results from these tests were used to correct the manufacturer's in-flight thrust program to more accurately calculate thrust for the specific test engine. Data from these tests were also used to develop an independent, simplified thrust calculation technique for real-time thrust calculation. Comparisons were also made to thrust values predicted by the engine specification model. Results indicate uninstalled gross thrust accuracies on the order of 1 to 4 percent for the various in-flight thrust methods. The various thrust calculations are described and their usage, uncertainty, and measured accuracies are explained. In addition, the advantages of a real-time thrust algorithm for flight test use and the importance of an accurate thrust calculation to the aircraft performance analysis are described. Finally, actual data obtained from flight test are presented.

Ray, Ronald J.

Real-time aerodynamic heating and surface temperature calculations for hypersonic flight simulation

A real-time heating algorithm was derived and installed on the Ames Research Center Dryden Flight Research Facility real-time flight simulator. This program can calculate two- and three-dimensional stagnation point surface heating rates and surface temperatures. The two-dimensional calculations can be made with or without leading-edge sweep. In addition, upper and lower surface heating rates and surface temperatures for flat plates, wedges, and cones can be calculated. Laminar or turbulent heating can be calculated, with boundary-layer transition made a function of free-stream Reynolds number and free-stream Mach number. Real-time heating rates and surface temperatures calculated for a generic hypersonic vehicle are presented and compared with more exact values computed by a batch aeroheating program. As these comparisons show, the heating algorithm used on the flight simulator calculates surface heating rates and temperatures well within the accuracy required to evaluate flight profiles for acceptable heating trajectories.

Quinn, Robert D.

Calculations on the rate of the ion-molecule reaction between NH3(+) and H2

The rate coefficient for the ion-molecule reaction NH3(+) + H2 yields NH4(+) + H has been calculated as a function of temperature with the use of the statistical phase space approach. The potential surface and reaction complex and transition state parameters used in the calculation have been taken from ab initio quantum chemical calculations. The calculated rate coefficient has been found to mimic the unusual temperature dependence measured in the laboratory, in which the rate coefficient decreases with decreasing temperature until 50-100 K and then increases at still lower temperatures. Quantitative agreement between experimental and theoretical rate coefficients is satisfactory given the uncertainties in the ab initio results and in the dynamics calculations. The rate coefficient for the unusual three-body process NH3(+) + H2 + He yields NH4(+) + H + He has also been calculated as a function of temperature and the result found to agree well with a previous laboratory determination.

Herbst, Eric

Three dimensional model calculations of the global dispersion of high speed aircraft exhaust and implications for stratospheric ozone loss

Two-dimensional (zonally averaged) photochemical models are commonly used for calculations of ozone changes due to various perturbations. These include calculating the ozone change expected as a result of change in the lower stratospheric composition due to the exhaust of a fleet of supersonic aircraft flying in the lower stratosphere. However, zonal asymmetries are anticipated to be important to this sort of calculation. The aircraft are expected to be restricted from flying over land at supersonic speed due to sonic booms, thus the pollutant source will not be zonally symmetric. There is loss of pollutant through stratosphere/troposphere exchange, but these processes are spatially and temporally inhomogeneous. Asymmetry in the pollutant distribution contributes to the uncertainty in the ozone changes calculated with two dimensional models. Pollutant distributions for integrations of at least 1 year of continuous pollutant emissions along flight corridors are calculated using a three dimensional chemistry and transport model. These distributions indicate the importance of asymmetry in the pollutant distributions to evaluation of the impact of stratospheric aircraft on ozone. The implications of such pollutant asymmetries to assessment calculations are discussed, considering both homogeneous and heterogeneous reactions.

Douglass, Anne R.