Engineering Papers⌕ Search

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 181 records · Page 10

Calculation of the flow on a blunted cone at high angle of attack

A new technique for calculating the entire flow-field on spherically blunted circular cones at high angles of attack and high Reynolds numbers is described. The calculations are based on a single-layer system of three-dimensional parabolic equations which are approximations to the full steady Navier-Stokes equations. Initial conditions at the sphere-cone tangency plane are provided by using an inviscid time-dependent solution added to a viscous nonsimilar boundary layer solution. Calculated results are compared with experimental heat transfer and pressure data for a 15 deg half-angle cone with a 1-in. spherical nose at 15 deg angle of attack. The free-stream Mach number is 10.6, and the free-stream Reynolds number is 1,200,000 per foot. Excellent agreement between the calculated and experimental data for both pressure and heat transfer is obtained.

Lubard, S. C.↗

Calculation of autoionization of He and H/-/ using the projection-operator formalism

Improved calculations are reported for the first several autoionization states of the lower symmetries of He and H(-). Unshifted energies are calculated by diagonalizing QHQ using a Hylleraes basis with more terms than previously used; shifts, widths, and photoabsorption shape parameters are obtained with the additional use of exchange-approximate nonresonant continuum functions. Previous calculations of H(-) resonances are reviewed and slightly augmented by applying various nonresonant continua and correcting small errors. A comparison is made between the calculations and experimental results and is found to be very satisfactory except for the lowest 1P autoionization state of He, which is shown to need a more accurate experimental determination.

Bhatia, A. K.↗

Calculations of the radiative and dynamical state of the Venus atmosphere

Results are reported for accurate multiple-scattering calculations to determine the solar-energy deposition profile in the atmosphere of Venus. It is found that most of the absorbed energy is deposited in the main cloud-layer region, located at altitudes above 35 km, and that the ground receives approximately 3% of the energy absorbed in toto by Venus. Using these results, vertical temperature profiles are computed under conditions of pure radiative equilibrium and radiative-convective equilibrium. Since the latter results satisfactorily match the temperature structure determined from various spacecraft observations, it is inferred that the greenhouse effect can account for the high surface temperature. Aerosols make an important contribution to the infrared opacity in these calculations. Preliminary three-dimensional calculations of the general circulation of the atmosphere are discussed which incorporate the results of the radiative calculations.

Pollack, J. B.↗

Calculations of the Electron Energy Distribution Function in a Uranium Plasma by Analytic and Monte Carlo Techniques

Electron energy distribution functions were calculated in a U235 plasma at 1 atmosphere for various plasma temperatures and neutron fluxes. The distributions are assumed to be a summation of a high energy tail and a Maxwellian distribution. The sources of energetic electrons considered are the fission-fragment induced ionization of uranium and the electron induced ionization of uranium. The calculation of the high energy tail is reduced to an electron slowing down calculation, from the most energetic source to the energy where the electron is assumed to be incorporated into the Maxwellian distribution. The pertinent collisional processes are electron-electron scattering and electron induced ionization and excitation of uranium. Two distinct methods were employed in the calculation of the distributions. One method is based upon the assumption of continuous slowing and yields a distribution inversely proportional to the stopping power. An iteration scheme is utilized to include the secondary electron avalanche. In the other method, a governing equation is derived without assuming continuous electron slowing. This equation is solved by a Monte Carlo technique.

Bathke, C. G.↗

Calculation of subsonic and supersonic steady and unsteady aerodynamic forces using velocity potential aerodynamic elements

Expressions for calculation of subsonic and supersonic, steady and unsteady aerodynamic forces are derived, using the concept of aerodynamic elements applied to the downwash velocity potential method. Aerodynamic elements can be of arbitrary out of plane polygon shape, although numerical calculations are restricted to rectangular elements, and to the steady state case in the supersonic examples. It is suggested that the use of conforming, in place of rectangular elements, would give better results. Agreement with results for subsonic oscillating T tails is fair, but results do not converge as the number of collocation points is increased. This appears to be due to the form of expression used in the calculations. The methods derived are expected to facilitate automated flutter analysis on the computer. In particular, the aerodynamic element concept is consistent with finite element methods already used for structural analysis. The method is universal for the complete Mach number range, and, finally, the calculations can be arranged so that they do not have to be repeated completely for every reduced frequency.

Haviland, J. K.↗

A computer program to calculate the longitudinal aerodynamic characteristics of wing-flap configurations with externally blown flaps

A vortex lattice lifting-surface method is used to model the wing and multiple flaps. Each lifting surface may be of arbitrary planform having camber and twist, and the multiple-slotted trailing-edge flap system may consist of up to ten flaps with different spans and deflection angles. The engine wakes model consists of a series of closely spaced vortex rings with circular or elliptic cross sections. The rings are normal to a wake centerline which is free to move vertically and laterally to accommodate the local flow field beneath the wing and flaps. The two potential flow models are used in an iterative fashion to calculate the wing-flap loading distribution including the influence of the waves from up to two turbofan engines on the semispan. The method is limited to the condition where the flow and geometry of the configurations are symmetric about the vertical plane containing the wing root chord. The calculation procedure starts with arbitrarily positioned wake centerlines and the iterative calculation continues until the total configuration loading converges within a prescribed tolerance. Program results include total configuration forces and moments, individual lifting-surface load distributions, including pressure distributions, individual flap hinge moments, and flow field calculation at arbitrary field points.

Mendenhall, M. R.↗

Hybrid theory and calculation of e-N2 scattering

A theory of electron-molecule scattering is developed which is a synthesis of close-coupling and adiabatic-nuclei theories. Specifically, the theory is close-coupling with respect to vibrational degrees of freedom and adiabatic-nuclei with respect to rotation. It can be applied to any number of partial waves required; the remaining ones can be calculated purely in one or the other approximation. A theoretical criterion based on fixed-nuclei calculations is given which indicates those partial waves and energy domains requiring the various approximations. The theory allows all cross sections (pure rotational, vibrational, simultaneous vibration-rotation, differential, and total) to be calculated, and explicit formulas for all these cross sections are given. The theory is applied to low-energy e-N2 scattering. The fixed-nuclei results are such that the criterion shows clearly that vibrational close coupling is necessary, but only for the Pi sub g partial wave. It is found that the close-coupling calculation for this wave gives rise to the substructure as well as the gross structure of the 2.4-eV resonance and that vibrational excitation cross sections are about twice as large as previously inferred.

Chandra, N.↗

Twilight helium 10,830-A calculations and observations

The 10,830-A column emission rates from metastable He(2 3S) atoms in the upper atmosphere are calculated and compared with other calculations and observations. It is shown that the calculated emission rates agree well with the observed emission rates provided that the cross sections for loss of He(2 3S) by Penning ionization on atomic oxygen is near 5 A 2 and not 44 A 2 as previously measured. It is assumed that the local and conjugate photoelectron fluxes above 300 km are about double those calculated. The dawn/dusk intensity ratio observed in New Mexico, Brazil and Peru is accounted for by changes in the ambient thermal electron concentration in Brazil and Peru and by changes in plasmasphere opacity to conjugate point electrons for New Mexico and to some extent for Brazil. The rate at which He-4 would escape from the earth if sufficient energy were imparted to the He-4 atoms during the Penning reaction was recalculated and found to be inadequate by an order of magnitude to balance the terrestrial production.

Tinsley, B. A.↗

A calculation of Saturn's gravitational contraction history

A gravitational-contraction history is calculated for a homogeneous quasi-equilibrium Saturn model of solar composition without a rock-ice core and without allowance for the possible separation of hydrogen and helium. The calculations begin at the time when Saturn's radius was ten times its current value, and the subsequent gravitational contraction is followed for 4.5 billion years. The results obtained are given in terms of the path of the evolving model on the H-R diagram, the variation of central temperature with central density, the time variation of radius and internal luminosity, and the energy changes that take place during the evolution. These results are compared with those of similar calculations for Jupiter, and the internal structure of the evolutionary models is examined in detail. The Saturn calculations yield a radius at the current epoch that is 9% larger than the observed value and an excess luminosity comparable to that observed. It is noted that an inhomogeneous model containing a solar-mix envelope and a central rocky core has a radius equal to the observed value.

Pollack, J. B.↗

Calculated rate constants for the reaction ClO + O yields Cl + O2 between 220 and 1000 deg K

Classical trajectory calculations are presented for the reaction ClO + O yields Cl + O2, a reaction which is an important step in the chlorine-catalyzed destruction of ozone which is thought to occur in the 220 and 1000 K. The calculated rate constant is 4.36 x 10 to the minus 11th power exp (-191/T)cu cm molecule (-1)s(-1) and its value at 300 K is 2.3 plus or minus 10 to the 11th power cu cm molecule (-1)s(-1), about a factor of 2 lower than recent experimental data. The empirical potential energy surface used in the calculations was constructed to fit experimental data for ClO, O2 and ClOO molecules. Other important features of this potential surface, such as the barrier to reaction, were varied systematically and calculations were performed for a range of conditions to determine the best theoretical rate constants. Results demonstrate the utility of classical trajectory methods for determining activation energies and other kinetic data for important atmospheric reactions.

Jaffee, R. L.↗

Transonic flow over the NACA 64A006 with an oscillating flap-calculations based on the Euler equations

Exploratory calculations of transonic flows over the airfoil with a quarter-chord oscillating flap were made using a program which obtains approximate solutions to the Euler equations with an explicit, shock-capturing, finite-difference scheme. The calculations, essentially inviscid and for the airfoil at zero angle-of-attack in a free-stream, are at Mach numbers and reduced frequencies which were tested in experiments by Tijdeman. The oscillatory lifts from analogus calculations by various investigators generally agree with one another better than they agree with Tijdeman's data. Inclusion in the calculations of an approximate modelling of boundary conditions expected at slotted wind tunnel walls tends to shift some of the results closer to the experimental values.

Magnus, R. J.↗

Computer program to calculate three-dimensional boundary layer flows over wings with wall mass transfer

A system of computer programs for calculating three dimensional transonic flow over wings, including details of the three dimensional viscous boundary layer flow, was developed. The flow is calculated in two overlapping regions: an outer potential flow region, and a boundary layer region in which the first order, three dimensional boundary layer equations are numerically solved. A consistent matching of the two solutions is achieved iteratively, thus taking into account viscous-inviscid interaction. For the inviscid outer flow calculations, the Jameson-Caughey transonic wing program FLO 27 is used, and the boundary layer calculations are performed by a finite difference boundary layer prediction program. Interface programs provide communication between the two basic flow analysis programs. Computed results are presented for the NASA F8 research wing, both with and without distributed surface suction.

Mclean, J. D.↗

Effect of element density on the NASTRAN calculated mechanical and thermal stresses of a spar

A NASTRAN model of a spar was examined to determine the sensitivity of calculated axial thermal stresses and bending stresses to changes in element density of the model. The thermal stresses calculated with three different element densities resulted in drastically differing values. The position of the constraint also significantly affected the value of the calculated thermal stresses. Mechanical stresses calculated from an applied loading were insensitive to element density.

Jenkins, J. M.↗

SCF treatment of charge polarization effects in intermediate-energy electron scattering calculations with applications to N2

Converged rotational close coupling calculations of differential, integral, and momentum-transfer cross sections are obtained for seven model potentials for electron-N2 scattering at an impact energy of 30 eV. The model potentials involve a static potential calculated by the INDO/1s or INDOXI/1s method, and exchange potential calculated by the semiclassical exchange approximation from the INDO/1s or INDOXI/1s unperturbed electron density and a polarization potential. The polarization potentials used include the Buckley-Burke semiempirical one and various modifications of the INDOXI and INDO SCF (self-consistent-field) adiabatic polarization potentials. Without adjustment of parameters, good agreement is obtained with the angle dependence of the experimentally measured sum of the elastic and rotational excitation differential cross sections; however, the absolute value of the calculated cross sections is about 20-30% larger than the measured values in the best case, perhaps indicating that the model potentials are too strong or should have a nonzero imaginary part.

Onda, K.↗

High resolution atmospheric transmission calculations down to 28.7 km in the 200-243-nm spectral range

Decrease in stratospheric ozone absorption and increase in oxygen absorption with decreasing wavelength combine to produce a window of maximum atmospheric transmission near 210 nm. Since solar radiation in this spectral region dissociates molecular oxygen, the deep atmospheric penetration at this wavelength is of particular aeronomical interest. High resolution calculations of the transmittance down to 28.65 km were made for the 200-243-nm spectral range in this window region, in support of a stratospheric balloon flight from Fort Churchill in July 1974. The calculations were made by dividing the atmosphere into layers which were chosen so that each could be assumed homogeneous; optical depths were calculated separately for each of these layers and then summed to obtain the over-all transmittance of the atmosphere. Absorption by molecular oxygen (line and continuum) and by ozone was included, as well as extinction through Rayleigh scattering by air molecules. The calculated transmittances were combined with high altitude (above 100-km) rocket measurements of the sun-center spectrum and center-to-limb variations to give residual high resolution solar spectral flux for several altitudes and solar zenith angles.

Cann, M. W. P.↗

Extensions to the classical calculation of the effect of mutual shadowing in diffuse reflection

The classical method for accounting for the mutual shadowing among closely packed particles in multiple scattering calculations is extended in the following ways. (1) By modeling the particle distribution by a Poisson process with a varying density parameter, a 'Van der Waals' type approximation allows extension to a greater fractional volume density, D. In this case it is only required that D squared be much less than 1 instead of D being much less than 1. (2) In the case that the particle distribution is not uniform the classical calculation may be weighted by the pair correlation function of the distribution. (3) The use of the Markov chain formalism for radiative transfer allows inclusion of the effect of shadowing for two orders of scattering. For conditions such as might apply in Saturn's rings, the inclusion of this effect makes less than a 0.1% difference in the calculated phase curves, compared to previous calculations which have included shadowing only in the first scattering. The latter are thus shown to be quite accurate.

Esposito, L. W.↗

Comparison of calculated and measured model rotor loading and wake geometry

The calculated blade bound circulation and wake geometry are compared with measured results for a model helicopter rotor in hover and forward flight. Hover results are presented for rectangular tip and ogee tip planform blades. The correlation is quite good when the measured wake geometry characteristics are used in the analysis. Available prescribed wake geometry models are found to give fair predictions of the loading, but they do not produce a reasonable prediction of the induced power. Forward flight results are presented for twisted and untwisted blades. Fair correlation between measurements and calculations is found for the bound circulation distribution on the advancing side. The tip vortex geometry in the vicinity of the advancing blade in forward flight was predicted well by the free wake calculation used, although the wake geometry did not have a significant influence on the calculated loading and performance for the cases considered.

Johnson, W.↗

Calculation of Water Drop Trajectories to and About Arbitrary Three-Dimensional Bodies in Potential Airflow

Calculations can be performed for any atmospheric conditions and for all water drop sizes, from the smallest cloud droplet to large raindrops. Any subsonic, external, non-lifting flow can be accommodated; flow into, but not through, inlets also can be simulated. Experimental water drop drag relations are used in the water drop equations of motion and effects of gravity settling are included. Seven codes are described: (1) a code used to debug and plot body surface description data; (2) a code that processes the body surface data to yield the potential flow field; (3) a code that computes flow velocities at arrays of points in space; (4) a code that computes water drop trajectories from an array of points in space; (5) a code that computes water drop trajectories and fluxes to arbitrary target points; (6) a code that computes water drop trajectories tangent to the body; and (7) a code that produces stereo pair plots which include both the body and trajectories. Code descriptions include operating instructions, card inputs and printouts for example problems, and listing of the FORTRAN codes. Accuracy of the calculations is discussed, and trajectory calculation results are compared with prior calculations and with experimental data.

Norment, H. G.↗