SEARCH · Engineering Papers
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.
A simple modification of Gilbert's formula to account for slight lateral heterogeneity of the Earth leads to a convenient formula to calculate synthetic long period seismograms. Partial derivatives are easily calculated, thus the formula is suitable for direct inversion of seismograms for lateral heterogeneity of the Earth.
A new method of reliability prediction for complex systems is defined. Calculation of both upper and lower bounds are involved, and a procedure for combining the two to yield an approximately true prediction value is presented. Both mission success and crew safety predictions can be calculated, and success probabilities can be obtained for individual mission phases or subsystems. Primary consideration is given to evaluating cases involving zero or one failure per subsystem, and the results of these evaluations are then used for analyzing multiple failure cases. Extensive development is provided for the overall mission success and crew safety equations for both the upper and lower bounds.
A chemical reaction scheme comprising gas phase as well as grain surface reactions is used for the numerical calculation of the chemical composition of mantles accreting on interstellar grains. The composition of the growing grain mantle is determined on the basis of the relative accretion rates of the gas phase molecules and diffusion-controlled surface reactions. Results show that grain mantles generally consist of the molecules H2O, H2CO, N2, O2, CO, CO2, H2O2, and NH3. A detailed examination of the formation of H2 on grain surfaces leads to the conclusion that it proceeds through hydrogen abstraction from such molecules as H2CO, H2S, N2H2, and N2H4. The IR characteristics of the calculated grain mantles are discussed with emphasis on the observed ice band at the 3250/cm wavelength.
A method for calculating the longitudinal aerodynamic coefficients and the pressure distributions on a body at reasonably high angles of attack is presented. The body is represented by a combination of source elements and vortex-lattice elements, including separation of the vortices at increasing angles of attack. The method is self-consistent in that the body and the separated vortex wake are treated as an integrated interacting system. The location of the separation line can be included as an arbitrary input from experimental data or can be evaluated approximately by a pressure-dependent criterion. Calculated values of the aerodynamic coefficients and pressure distributions on cone-cylinder and ogive-cylinder bodies compare well, qualitatively and quantitatively, with experimental data, including simulation of the dependence on Reynolds number.
Procedures for calculating bound-bound and bound-continuum (free) radiative transition probabilities and radiative lifetimes are summarized. Calculations include rotational dependence and R-dependent electronic transition moments (no Franck-Condon or R-centroid approximation). Detailed comparisons of theoretical results with experimental measurements are made for bound-bound transitions in the A-X systems of LiH and Na2. New bound-free results are presented for LiH. New bound-free results and comparisons with very recent fluorescence experiments are presented for Na2.
The atomic data presented here and in Bhatia and Mason (1980) allow the calculation of theoretical intensity ratios for all the EUV, UV, and X-ray lines from Fe XX. Tabulations are presently given for the transitions between levels in the 2s2 2p3, 2s2 2p2 3s, and 2s2 2p2 3d configurations of Fe(19+), and electron collision strengths are calculated by means of the 'distorted wave' approximation. In addition to the theoretical X-ray line intensity ratios, new spectral line identifications from a solar flare are presented.
Calculations of electron impact collision strengths and spontaneous radiative decay rates are made for neon-like krypton (Kr XXVII) for the 2s2 2p6, 2s2 2p5 3s, 2s2 2p5 3p, and 2s2 2p5 3d configurations. From these atomic data, the level populations as a function of the electron density are calculated at two temperatures, 1 x 10 to the 7th K and 3 x 10 to the 7th K. An analysis of level populations reveals that a volume of krypton in which a significant number of the ions are in the Kr XXVII degree of ionization can produce a significant gain in transition between the 2s2 2p5 3s and 2s2 2p5 3p configurations. At an electron density of 1 x 10 to the 19th/cu cm the plasma length has to be of the order of 1 m; at a density of 1 x 10 to the 21st/cu cm the length is reduced to approximately 0.5 cm; and at an electron density of 1 x 10 to the 22nd/cu cm the length of the plasma is further reduced to approximately 1 mm.
Navier-Stokes calculations have been performed for a supercritical airfoil at a transonic design condition and at a subsonic condition. Wind-tunnel pressure-rail measurements were employed as boundary data in the calculations to account for wall-interference effects. A fine mesh was used so that most details of the flows were resolved, particular attention having been given to the trailing-edge region. Detailed comparisons are made with the experimental data. Good agreement was obtained on the airfoil except at the trailing edge where separation occurred. Flow details in the trailing-edge region are examined and differences are shown to be attributable to the turbulence model employed.
A method is given for determining the sensitivity of certain types of calculations to the uncertainties in the input physics or model parameters; this method is applicable to problems that involve solutions to coupled, ordinary differential equations. In particular the sensitivity of calculations of the thermal structure of neutron star envelopes to uncertainties in the opacity and equation of state is examined. It is found that the uncertainties in the relationship between the surface and interior temperatures of a neutron star are due almost entirely to the imprecision in the values of the conductive opacity in the region where the ions form a liquid; here the conductive opacity is, for the most part, due to the scattering of electrons from ions.
An unsteady inverse boundary-layer method is developed which can be used to calculate steady flows with separation. Two versions of Keller's box method with the Mechul function formulation developed by Cebeci (1976) are employed, depending on the complexity of the flow. The regular box is employed in regions of positive streamwise velocity component u, whereas the zigzag box is employed in regions where u becomes negative (t greater than 0). The regular box with the FLARE approximation is employed when t = 0 and u becomes negative in some region across the layer. Results of calculations show that the use of a time-dependent inverse boundary-layer method in which time is used as an iteration parameter provides a good approach in improving the accuracy of the solutions obtained from the FLARE approximation.
Tangents and normals to ellipse derived for situations with two correlated aerodynamic force or moment variables having normally (Gaussian) distributed uncertainty. Although transformations and vector calculations tedious when done by hand, equations in closed form and give rapid results when performed by computers or programable desk calculators. Probability extremes of ratios and magnitude obtained directly, without repeated table lookup or other complication.
This paper presents a technique for calculating the modal participation factors for base excitation problems using a DMAP alter to the NASTRAN real eigenvalue analysis Rigid Format. The DMAP program automates the generation of the seismic mass to add to the degrees of freedom representing the shaker input directions and calculates the modal participation factors. These are shown in the paper to be a good measure of the maximum acceleration expected at any point on the structure when the subsequent frequency response analysis is run.
Using the quasi-steady, full potential code, ROT22, pressures were calculated on straight and swept tip model helicopter rotor blades at advance ratios of 0.40 and 0.45, and into the transonic tip speed range. The calculated pressures were compared with values measured in the tip regions of the model blades. Good agreement was found over a wide range of azimuth angles when the shocks on the blade were not too strong. However, strong shocks persisted longer than predicted by ROT22 when the blade was in the second quadrant. Since the unsteady flow effects present at high advance ratios primarily affect shock waves, the underprediction of shock strengths is attributed to the simplifying, quasi-steady, assumption made in ROT22.
A modified local plasma model is established. The Gordon-Kim's molecular charged density model is employed to obtain a formula to evaluate the stopping power of many useful molecular systems. The stopping power of H2 and He gas was calculated for incident proton energy ranging from 100 keV to 2.5 MeV. The stopping power of O2, N2, and water vapor was also calculated for incident proton energy ranging from 40 keV. to 2.5 MeV. Good agreement with experimental data was obtained. A discussion of molecular effects leading to department from Bragg's rule is presented. The equipartition rule and the effect of nuclear momentum recoiling in stopping power are also discussed.
MERIDL is a program that calculates a meridional plane stream function solution, and TSONIC is a program that calculates a blade to blade stream function solution for turbomachine blade passages. Both programs are discussed, including input required and assumptions and limitations. Examples of use and references are included.
To identify the fastest algorithm currently available for the numerical integration of chemical kinetic rate equations, several algorithms were examined. Findings to date are summarized. The algorithms examined include two general-purpose codes EPISODE and LSODE and three special-purpose (for chemical kinetic calculations) codes CHEMEQ, CRK1D, and GCKP84. In addition, an explicit Runge-Kutta-Merson differential equation solver (IMSL Routine DASCRU) is used to illustrate the problems associated with integrating chemical kinetic rate equations by a classical method. Algorithms were applied to two test problems drawn from combustion kinetics. These problems included all three combustion regimes: induction, heat release and equilibration. Variations of the temperature and species mole fraction are given with time for test problems 1 and 2, respectively. Both test problems were integrated over a time interval of 1 ms in order to obtain near-equilibration of all species and temperature. Of the codes examined in this study, only CREK1D and GCDP84 were written explicitly for integrating exothermic, non-isothermal combustion rate equations. These therefore have built-in procedures for calculating the temperature.
A computer model which describes the solidification of a binary metal alloy in an insulated rectangular mold with a temperature gradient is presented. A numerical technique, applicable to a broad class of moving boundary problems, was implemented therein. The solidification model described is used to calculate the macrosegregation within the solidified casting by coupling the equations for liquid flow in the solid/liquid or mushy zone with the energy equation for heat flow throughout the ingot and thermal convection in the bulk liquid portion. The rate of development of the solid can be automatically calculated by the model. Numerical analysis of such solidification parameters as enthalpy and boundary layer flow is displayed. On-line user interface and software documentation are presented.
An algorithm for calculating the atmospheric transmittance in the 10 to 20 micro m spectral band from a known temperature and dewpoint profile, and then using this transmittance to estimate the surface (skin) temperature from a VISSR observation in the 11 micro m window is presented. Parameterizations are drawn from the literature for computing the molecular absorption due to the water vapor continuum, water vapor lines, and carbon dioxide lines. The FORTRAN code is documented for this application, and the sensitivity of the derived skin temperature to variations in the model's parameters is calculated. The VISSR calibration uncertainties are identified as the largest potential source of error.