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.
Practical computational techniques developed. Report discusses practical considerations of calculation of three-dimensional radiative transfer. Extension of earlier theoretical developments described in "Calculating Atmospheric Effects in Satellite Imagery" NPO-16373.
Rotor/stator interactions analyzed with help of approximations. Analytical procedure assists designers of turbomachinery in calculating dynamic loads on blades and vanes. Procedure involves approximations simplifying calculation of phase angles and amplitudes of pressure fields throughout blade and vane rows at each disturbance frequency.
The spectroscopic parameters and separations between the three low-lying X 3 Pi u, A 3 Sigma g -, and a 1 Sigma g + states of Al2 are studied as a function of both the one-particle and n-particle basis set. Approximate correlation treatments are calibrated against full Cl calculations correlating the six valence electrons in a double-zeta plus two d-function basis set. Since the CASSCF/MRCI 3 Pi u to 3 Sigma g - separation is in excellent agreement wtih the FCI value, the MRCI calculations were carried out in an extended (20s13p6d4f)/(6s5p3d2f) gaussian basis. Including a small correction for relativistic effects, the best estimate is that 3 Sigma g - state lies 174/cm above the 3 Pi u ground state. The 1 Sigma g + state lies at least 2000/cm higher in energy. At the CPF level, inclusion of 2s and 2p correlation has little effect on D sub e, reduces T sub e by only 26/cm, and shortens the bond lengths by about 0.02 a sub o. Further strong support for a 3 Pi u ground state comes from the experimental absorption spectra, since both observed transitions can be convincingly assigned as 3 Pi u yields 3 Pi g. The (2) 3 Pi g state is observed to be sensitive to the level of correlation treatment, and to have its minimum shifted to shorter rho values, such that the strongest experimental absorption peak probably corresponds to the 0 yields 2 transition.
The method of calculating single-entry separation requirements for pairs of satellites is extended to include the interference on the top link as well as on the down link. Several heuristic models for analyzing the effects of shaped-beam antenna designs on required satellite separations are introduced and demonstrated with gain contour plots. The calculation of aggregate interference is extended to include the effects of up-link interference. The relationship between the single-entry C/I requirements, used in determining satellite separation constraints for various optimization procedures, and the aggregate C/I values of the resulting solutions is discussed.
A key step in the application of formal automated design techniques to structures under transient loading is the calculation of sensitivities of response quantities to the design parameters. This paper considers structures with general forms of damping acted on by general transient loading and addresses issues of computational errors and computational efficiency. The equations of motion are reduced using the traditional basis of vibration modes and then integrated using a highly accurate, explicit integration technique. A critical point constraint formulation is used to place constraints on the magnitude of each response quantity as a function of time. Three different techniques for calculating sensitivities of the critical point constraints are presented. The first two are based on the straightforward application of the forward and central difference operators, respectively. The third is based on explicit differentiation of the equations of motion. Condition errors, finite difference truncation errors, and modal convergence errors for the three techniques are compared by applying them to a simple five-span-beam problem. Sensitivity results are presented for two different transient loading conditions and for both damped and undamped cases.
This paper is concerned with the application of the acoustic analogy of Lighthill to the acoustic and aerodynamic problems associated with moving bodies. The Ffowcs Williams-Hawkings equation, which is an interpretation of the acoustic analogy for sound generation by moving bodies, manipulates the source terms into surface and volume sources. Quite often in practice the volume sources, or quadrupoles, are neglected for various reasons. Recently, Farassat, Long and others have attempted to use the FW-H equation with the quadrupole source and neglected to solve for the surface pressure on the body. The purpose of this paper is to examine the contribution of the quadrupole source to the acoustic pressure and body surface pressure for some problems for which the exact solution is known. The inviscid, incompressible, 2-D flow, calculated using the velocity potential, is used to calculate the individual contributions of the various surface and volume source terms in the FW-H equation. The relative importance of each of the sources is then assessed.
A 2-D steady-state Navier-Stokes solver has been upgraded to include the effects of frozen and equilibrium air chemistry for applications to high speed flight vehicles. To provide a computationally economical first order approximation to the high temperature physics, variable thermodynamic data is used for the chemically frozen mode to allow for a variation with temperature of the air specific heats and enthalpy. For calculations involving air in chemical equilibrium, a specially modified version of the NASA Lewis Chemical Equilibrium Code, CEC, is used to compute the chemical composition and resultant thermochemical properties. The upgraded solver is demonstrated by comparing results from calorically perfect (C sub p=constant), thermally perfect (frozen) and equilibrium air calculations for a variety of geometries, and flight Mach numbers.
Full-configuration-interaction (FCI) calculations have been performed for the Z 1Sigma(g)(+), a 3Pi(u), and b 3Sigma(g)(-) states of C2 and the X 3Sigma(g)(-) and A 3Pi(u) states of Si2 in valence double-zeta plus polarization one-particle basis sets. By carrying out the CASSCF multireference singles plus doubles CI (MRCI) procedure in a very large Gaussian basis set and correcting for errors with the FCI, the triplet-triplet separation in C2 is reproduced to within 47/cm. Analogous calculations on Si2 definitively determine the ground state to be 3Sigma(g)(-); the 3Pi(u) state is only 440 + or - 100/cm higher in energy. For the SiC molecule, which has never been observed spectroscopically but is believed to be an important component in carbon stars and to be present in interstellar space, the A 3Sigma(-) state is nearly 4000/cm above the X 3Pi ground state.
An interactive inviscid core flow-boundary layer method is presented for the calculation of turbomachine channel flows. For this method, a one-dimensional inviscid core flow is assumed. The end-wall and blade surface boundary layers are calculated using an integral entrainment method. The boundary layers are assumed to be collateral and thus are two-dimensional. The boundary layer equations are written in a streamline coordinate system. The streamwise velocity profiles are approximated by power law profiles. Compressibility is accounted for in the streamwise direction but not in the normal direction. Equations are derived for the special cases of conical and two-dimensional rectangular diffusers. For these cases, the assumptions of a one-dimensional core flow and collateral boundary layers are valid. Results using the method are compared with experiment and good quantitative agreement is obtained.
The usual assumption by which chemical reaction rates are calculated in two-dimensional atmospheric models is by using a product of zonal means of rate coefficients and constituent concentrations rather than the rigorous zonal mean of the corresponding products. This assumption has been tested for the reactions O + NO2 yields NO + O2 and NO + O3 yields NO2 + O2 using mapped limb infrared monitor of the stratosphere data from the Nimbus 7 satellite and found to be quite satisfactory for winter 1979 at 60 deg N in the upper stratosphere. Relative differences between the two-dimensional averaged rate and the more rigorous rate, calculated from the full, longitudinally varying temperatures and mixing ratios, were small (usually below 5 percent) and exceeded 15 percent only during times of strong dynamical activity. At those times or locations where stratospheric circulation is primarily zonal, the two averages agreed to within a few percent.
After calculating ambient electron densities as a function of altitude, latitude and local time, a simplified expression is used to calculate F-region eastward plasma drifts given a zonal neutral wind model. The derived eastward plasma drift on a magnetic flux tube is examined as a function of the flux tube apex height. If the neutral wind is assumed to be independent of latitude the plasma drift maximizes along the flux tube which intercepts the F-region peak concentration at the Appleton anomaly. Above this altitude the velocity decreases to reflect a decrease in the flux tube integrated F-region Pedersen conductivity. For a latitude dependent wind the plasma drift tends to maximize along the flux tube which intercepts the F-peak at the dip equator. Above this altitude the drift decreases to reflect the latitude distribution of the wind.
Numerical calculations are presented for the high-pressure dependence of the linewidths of the D2 polarized Q branch of the Stokes Raman nu = 0-1 transition. A comparison is made with the experimental linewidths. The relationship between the line shape and S matrices is reviewed as well as the calculation of the latter for the D2-He system.
Calculations of the influence of atmospheric conditions on solar cell short-circuit current (Isc) are made using a recently developed computer model for solar spectral irradiance distribution. The results isolate the dependence of Isc on changes in the spectral irradiance distribution without the direct influence of the total irradiance level. The calculated direct normal irradiance and percent diffuse irradiance are given as a reference to indicate the expected irradiance levels. This method can be applied to the calibration of photovoltaic reference cells. Graphic examples are provided for amorphous silicon and monocrystalline silicon solar cells under direct normal and global normal solar irradiances.
The effect of streamline geometry and pressure distributions on surface heating rates is examined for slender, spherically blunted cones. The modifications to the approximate aeroheating code include a curve fit of pressures computed by an Euler solution over a range of Mach numbers and cone angles. The streamline geometry is then found using the surface pressures and inviscid surface properties. Previously, streamlines were determined using the inviscid properties at the edge of the boundary layer when accounting for the effects of entropy-layer swallowing. Streamline calculations are now based on inviscid surface conditions rather than boundary-layer edge properties. However, the heating rates are calculated using inviscid properties at the edge of the boundary layer. Resulting heating rates compare favorably with solutions from the viscous-shock-layer equations.
Data on monthly ozone abundance trends near the stratopause, observed by the Ultraviolet Spectrometer (UVS) on the SME and by the Solar Backscatter Ultraviolet Instrument (SBUV) on NIMBUS-7 are presented for June, September, and January of the years 1982-1986. Globally averaged trends determined from the SME data (-0.5 + or - 1.3 percent/yr) were found to fall within model calculations by Rusch and Clancy (1988); the SBUV trends, on the other hand, were found to exceed maximum predicted ozone decreases by a factor of 3 or more. Detailed comparison of the two data sets indicated that an absolute offset of 3 percent/yr accounts for much of the difference between the two trends; the offset is considered to be due to incomplete characterization of the SBUV calibration drift. Both the UVS and SBUV data exhibited similar seasonal and latitudinal variations in ozone trends, which were reproduced by photochemical model calculations that included latitude-dependent NMC temperature trends over the 1982-1986 period.
BUNVIS-RG computer program designed to calculate vibration frequencies or buckling loads of prestressed lattice structures used in outer space. For buckling and vibration problems, BUNVIS-RG calculates deadload axial forces caused in members by any combination of externally-applied static point forces and moments at nodes, axial preload or prestrain in members, and such acceleration loads as those due to gravity. BUNVIS-RG is FORTRAN 77 computer program implemented on CDC CYBER and VAX computer.
Cumulative-binomial computer program, CUMBIN, one of set of three programs, calculates cumulative binomial probability distributions for arbitrary inputs. CUMBIN, NEWTONP (NPO-17556), and CROSSER (NPO-17557), used independently of one another. Reliabilities and availabilities of k-out-of-n systems analyzed. Used by statisticians and users of statistical procedures, test planners, designers, and numerical analysts. Used for calculations of reliability and availability. Program written in C.
Noniterative integral-entrainment method yields good approximations. Method of approximate calculation of flow in channel of turbomachine based on interaction of viscous flow in boundary layers with inviscid flow in core layer. Faster and more robust than other approximate methods of same type. Suitable for use in preliminary calculations for design and for off-design operation of turbomachinery. Flows in conical diffuser channels represented by two-dimensional boundary-layer and one-dimensional core flows described by equations of new method.