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.
By introducing variable section silencers in divergent and convergent sections or elbows of ventilation plants, the length of classical silencers can be reduced. Presented is an original calculation method for active noise silencers with a continuously variable section, as well as calculation relations worked out by this method for eight common practical cases.
Iterative finite-difference techniques are developed for integrating the boundary-layer equations, without approximation, through a region of reversed flow. The numerical procedures are used to calculate incompressible laminar separated flows and to investigate the conditions for regular behavior at the point of separation. Regular flows are shown to be characterized by an integrable saddle-type singularity that makes it difficult to obtain numerical solutions which pass continuously into the separated region. The singularity is removed and continuous solutions ensured by specifying the wall shear distribution and computing the pressure gradient as part of the solution. Calculated results are presented for several separated flows and the accuracy of the method is verified. A computer program listing and complete solution case are included.
A method, developed for making measurements in reverberation chambers in the steady state mode, is reported. Spatial distribution patterns of the square of the sound pressure vs type of exciting signal were produced theoretically and checked experimentally in a reverberation chamber. Results of measurements of spatial irregularities of the field with pure tones and with noise bands in the 200 to 2,000 Hz frequency range were used for calculating statistical parameters of the reverbation fields by computer. The calculation program included determination of density probability function, integral distribution pattern, mathematical expectation, normalized dispersion of the mean square of the pressure.
The method and basic elements of computer solutions for both potential flow and viscous flow calculations for engine inlets are described. The procedure is applicable to subsonic conventional, short-haul, and vertical takeoff aircraft engine nacelles operating in a compressible viscous flow. The calculated results compare well with measured surface-pressure distributions for a number of model inlets. The uses of the program in both the design and analysis of engine inlets are discussed, with several examples given for VTOL lift fans, acoustic splitters, and for STOL engine nacelles. Several test support applications are also given.
Using the optimized valence configurations (OVC) multiconfiguration self-consistent-field (MCSCF) method, the dipole moment function for the ground state of CO in the vicinity of the equilibrium internuclear distance has been calculated. The OVC MCSCF calculation results are compared with existing Hartree-Fock and configuration interaction treatments of this molecule at single points and also the dipole moment function deduced from experimental infrared intensities. A general prescription for constructing OVC wavefunctions for diatomic molecules is also presented.
Description of a numerical technique for solving axisymmetric, incompressible, turbulent swirling flow problems. Isothermal flow calculations are presented for a coaxial flow configuration of special interest. The calculation results are discussed in regard to their implications for the design of gas turbine combustors.
The theory for calculating the absorption of laser radiation by hydrogen is outlined for the temperatures and pressures of common laboratory plasmas. Nonhydrogenic corrections for determining the absorption by helium are also included. The coefficients for the absorption of He-Ne laser radiation at the wavelengths of 0.633, 1.15, and 3.39 microns in a H plasma is presented for temperatures in the range from 10,000 to 40,000 K and electron number densities in the range from 10 to the 15th power to 10 to the 18th power per cu cm. The total absorption of a H-He plasma calculated from this theory is compared with the measured absorption. The theoretical composition of the H-He absorption is analyzed with respect to the significant absorption processes, inverse bremsstrahlung, photoionization, resonance excitation, and photodetachment.
Calculational procedure evolved in process of designing inlets. Douglas axisymmetric potential flow program called EOD calculates incompressible potential flow about arbitrary bodies. Program SCIRCL generates input for EOD from inlet components. Program COMBYN takes basic solutions output by EOD and combines them into solutions of interest and applied compressibility correction.
Methods for calculating the transfer of infrared radiation in planetary atmospheres are summarized. Direct spectral integration is compared with high precision measurements, suggesting that calculational techniques are sufficient for molecular atmospheres. Infrared transfer is used to study the atmospheric energy budget and atmospheric structures.
The use of linear theory in calculating the twist distribution of a wing designed for cruise at supercritical speeds is justified on the basis of the transonic equivalence rule. A modified version of Multhopp's subsonic lifting-surface theory was used to calculate the twist distribution. The lifting-surface theory is compared with both slender wing theory and experimental results. A study was also made of the effect of wing sweep on the twist distribution required to maintain an elliptic span load at cruise conditions. The important parameters used in establishing this twist distribution are identified.
A numerical scheme using the method of characteristics to calculate the flow properties and pressures behind decaying shock waves for materials under hypervelocity impact is developed. Time-consuming double interpolation subroutines are replaced by a technique based on orthogonal polynomial least square surface fits. Typical calculated results are given and compared with the double interpolation results. The complete computer program is included.
Film density data that have been digitalized and recorded in binary coded decimal (BCD) format are converted into a number representing the film density difference between the unexposed film border and any point on a photograph by program CONVERT. The difference AVERAGE, the BCD number representing the unexposed film area and the BCD number of a specific data point is calculated and then multiplied by three constants to produce the desired film density difference. The angle off the principal axis of the lens is computed for each data point in order that a correction factor be applied to compensate for atmospheric backscattering differences and light intensity fall-off, inherent in the camera lens system in use. The program is capable of plotting the calculated density differences. The percentage of points in specific density ranges can also be computed.
A FORTRAN IV subprogram called GASP is discussed which calculates the thermodynamic and transport properties for 10 pure fluids: parahydrogen, helium, neon, methane, nitrogen, carbon monoxide, oxygen, fluorine, argon, and carbon dioxide. The pressure range is generally from 0.1 to 400 atmospheres (to 100 atm for helium and to 1000 atm for hydrogen). The temperature ranges are from the triple point to 300 K for neon; to 500 K for carbon monoxide, oxygen, and fluorine; to 600 K for methane and nitrogen; to 1000 K for argon and carbon dioxide; to 2000 K for hydrogen; and from 6 to 500 K for helium. GASP accepts any two of pressure, temperature and density as input conditions along with pressure, and either entropy or enthalpy. The properties available in any combination as output include temperature, density, pressure, entropy, enthalpy, specific heats, sonic velocity, viscosity, thermal conductivity, and surface tension. The subprogram design is modular so that the user can choose only those subroutines necessary to the calculations.
A method is described for calculating the position of a satellite at the instants when measurements are made on board. The initial conditions used were the mean orbital elements of the satellite and their time derivatives in one orbit. The results of the calculation are compared with those obtained by numerical integration, and it is found that results are identical at the beginning of an orbit, but change as the orbit progresses. The advantages and disadvantages of the analytical method are presented.
A description is given of the modification of a theory on the reflectance of particulate media so as to apply it to analysis of the infrared spectra obtained by the IRIS instrument on Mariner 9. With the aid of this theory and the optical constants of muscovite mica, quartz, andesite, anorthosite, diopside pyroxenite, and dunite, modeling calculations were made to refine previous estimates of the mineralogical composition of the Martian dust particles. These calculations suggest that a feldspar rich mixture is a very likely composition for the dust particles. The optical constants used for anorthosite and diopside pyroxenite were derived during this program from reflectance measurements. Those for the mica were derived from literature reflectance data. Finally, a computer program was written to invert the measured radiance data so as to obtain the absorption coefficient spectrum which should then be independent of the temperature profile and gaseous component effects.
A code was developed for the calculation of interplanetary solar proton fluences at 1 AU for the active years 1977 to 1983. The fluences are presented as functions of mission duration tau, energy threshold E, and confidence level Q. For a given combination of tau and Q, the routine determines whether ordinary or anomalously large events are to be considered, and in the latter case, the number of anomalously large events that are predicted by probabilistic theory for the specified mission duration. The code is described in detail. A listing and sample calculations are included.
A viscid-inviscid interaction procedure for the calculation of subsonic and transonic flow over a boattail was developed. This method couples a finite-difference inviscid analysis with an integral boundary-layer technique. Results indicate that the effect of the boundary layer is as important as an accurate inviscid method for this type of flow. Theoretical results from the solution of the full transonic-potential equation, including boundary layer effects, agree well with the experimental pressure distribution for a boattail. Use of the small disturbance transonic potential equation yielded results that did not agree well with the experimental results even when boundary-layer effects were included in the calculations.
A user's manual is presented for a computer program which calculates inviscid flow about lifting configurations in the free-stream Mach-number range from zero to low supersonic. Angles of attack of the order of the configuration thickness-length ratio and less can be calculated. An approximate formulation was used which accounts for shock waves, leading-edge separation and wind-tunnel wall effects.