Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “calculator”

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 1,441 records · Page 80

A nonproprietary, nonsecret program for calculating Stirling cryocoolers

A design program for an integrated Stirling cycle cryocooler was written on an IBM-PC computer. The program is easy to use and shows the trends and itemizes the losses. The calculated results were compared with some measured performance values. The program predicts somewhat optimistic performance and needs to be calibrated more with experimental measurements. Adding a multiplier to the friction factor can bring the calculated rsults in line with the limited test results so far available. The program is offered as a good framework on which to build a truly useful design program for all types of cryocoolers.

Martini, W. R.↗

The calculation of theoretical chromospheric models and the interpretation of solar spectra from rockets and spacecraft

Calculated results based on two chromospheric flare models F1 and F2 of Machado, et al., (1980) are presented. Two additional models are included: F1*, which has enhanced temperatures relative to the weak-flare model F1 in the upper photosphere and low chromosphere, and F3 which has enhanced temperatures relative to the strong flare model F2 in the upper chromosphere. Each model is specified by means of a given variation of the temperature as a function of column mass. The corresponding variation of particle density and the geometrical height scale are determined by assuming hydrostatic equilibrium. The coupled equations of statistical equilibrium is solved as is radiative transfer for H, H-, He I-II, C I-IV, Si I-II, Mg I-II, Fe, Al, O I-II, Na, and Ca II. The overall absorption and emission of radiation by lines throughout the spectrum is determined by means of a reduced set of opacities sampled from a compilation of over 10 to the 7th power individual lines. That the white flight flare continuum may arise by extreme chromospheric overheating as well as by an enhancement of the minimum temperature region is also shown. The radiative cooling rate calculations for our brightest flare model suggest that chromospheric overheating provides enhanced radiation that could cause significant heating deep in the flare atmosphere.

Avrett, E. H.↗

Calculating far-field radiated sound pressure levels from NASTRAN output

FAFRAP is a computer program which calculates far field radiated sound pressure levels from quantities computed by a NASTRAN direct frequency response analysis of an arbitrarily shaped structure. Fluid loading on the structure can be computed directly by NASTRAN or an added-mass approximation to fluid loading on the structure can be used. Output from FAFRAP includes tables of radiated sound pressure levels and several types of graphic output. FAFRAP results for monopole and dipole sources compare closely with an explicit calculation of the radiated sound pressure level for those sources.

Lipman, R. R.↗

Optical model calculations of heavy-ion target fragmentation

The fragmentation of target nuclei by relativistic protons and heavy ions is described within the context of a simple abrasion-ablation-final-state interaction model. Abrasion is described by a quantum mechanical formalism utilizing an optical model potential approximation. Nuclear charge distributions of the excited prefragments are calculated by both a hypergeometric distribution and a method based upon the zero-point oscillations of the giant dipole resonance. Excitation energies are estimated from the excess surface energy resulting from the abrasion process and the additional energy deposited by frictional spectator interactions of the abraded nucleons. The ablation probabilities are obtained from the EVA-3 computer program. Isotope production cross sections for the spallation of copper targets by relativistic protons and for the fragmenting of carbon targets by relativistic carbon, neon, and iron projectiles are calculated and compared with available experimental data.

Townsend, L. W.↗

Characterization of real gas properties for space shuttle main engine fuel turbine and performance calculations

Real thermodynamic and transport properties of hydrogen, steam, the SSME mixture, and air are developed. The SSME mixture properties are needed for the analysis of the space shuttle main engine fuel turbine. The mixture conditions for the gases, except air, are presented graphically over a temperature range from 800 to 1200 K, and a pressure range from 1 to 500 atm. Air properties are given over a temperature range of 320 to 500 K, which are within the bounds of the thermodynamics programs used, in order to provide mixture data which is more easily checked (than H2/H2O). The real gas property variation of the SSME mixture is quantified. Polynomial expressions, needed for future computer analysis, for viscosity, Prandtl number, and thermal conductivity are given for the H2/H2O SSME fuel turbine mixture at a pressure of 305 atm over a range of temperatures from 950 to 1140 K. These conditions are representative of the SSME turbine operation. Performance calculations are presented for the space shuttle main engine (SSME) fuel turbine. The calculations use the air equivalent concept. Progress towards obtaining the capability to evaluate the performance of the SSME fuel turbine, with the H2/H2O mixture, is described.

Harloff, G. J.↗

Calculations of axisymmetric vortex sheet roll-up using a panel and a filament model

A method for calculating the self-induced motion of a vortex sheet using discrete vortex elements is presented. Vortex panels and vortex filaments are used to simulate two-dimensional and axisymmetric vortex sheet roll-up. A straight forward application using vortex elements to simulate the motion of a disk of vorticity with an elliptic circulation distribution yields unsatisfactroy results where the vortex elements move in a chaotic manner. The difficulty is assumed to be due to the inability of a finite number of discrete vortex elements to model the singularity at the sheet edge and due to large velocity calculation errors which result from uneven sheet stretching. A model of the inner portion of the spiral is introduced to eliminate the difficulty with the sheet edge singularity. The model replaces the outermost portion of the sheet with a single vortex of equivalent circulation and a number of higher order terms which account for the asymmetry of the spiral. The resulting discrete vortex model is applied to both two-dimensional and axisymmetric sheets. The two-dimensional roll-up is compared to the solution for a semi-infinite sheet with good results.

Kantelis, J. P.↗

Integral processing in beyond-Hartree-Fock calculations

The increasing rate at which improvements in processing capacity outstrip improvements in input/output performance of large computers has led to recent attempts to bypass generation of a disk-based integral file. The direct self-consistent field (SCF) method of Almlof and co-workers represents a very successful implementation of this approach. This paper is concerned with the extension of this general approach to configuration interaction (CI) and multiconfiguration-self-consistent field (MCSCF) calculations. After a discussion of the particular types of molecular orbital (MO) integrals for which -- at least for most current generation machines -- disk-based storage seems unavoidable, it is shown how all the necessary integrals can be obtained as matrix elements of Coulomb and exchange operators that can be calculated using a direct approach. Computational implementations of such a scheme are discussed.

Taylor, P. R.↗

Short wavelength laser calculations for electron pumping in Be I and B I isoelectronic sequences (Z greater than or equal to 18 and less than or equal to 36)

There has been considerable interest in creating large laser gains in the X-ray region from atomic transitions in ne-like ions formed in high-temperature plasmas. Although the Ne-like ions habe received the most attention, it is also possible to produce population inversions in ions of other isoelectronic sequences. Accordingly, it is shown that certain levels of the 2s3p and 2s(2)3p configurations in the Be I and B I isoelectronic sequences can be inverted by collisional pumping relative to levels of the 2s3s and 2s(2)3s configurations. Level populations and gain coefficients were calculated for the most promising transitions, and opacities were calculated for 2s2p-2s3s and 2s(2)2p-2s(2)3s transitions.

Feldman, U.↗

Definition and verification of a complex aircraft for aerodynamic calculations

Techniques are reviewed which are of value in CAD/CAM CFD studies of the geometries of new fighter aircraft. In order to refine the computations of the flows to take advantage of the computing power available from supercomputers, it is often necessary to interpolate the geometry of the mesh selected for the numerical analysis of the aircraft shape. Interpolating the geometry permits a higher level of detail in calculations of the flow past specific regions of a design. A microprocessor-based mathematics engine is described for fast image manipulation and rotation to verify that the interpolated geometry will correspond to the design geometry in order to ensure that the flow calculations will remain valid through the interpolation. Applications of the image manipulation system to verify geometrical representations with wire-frame and shaded-surface images are described.

Edwards, T. A.↗

Calculations of two and three-dimensional transonic cascade flow fields using the Navier-Stokes equations

The multidimensional, ensemble-averaged, compressible, time-dependent Navier-Stokes equations have been used to study the turbulent flow field in two and three-dimensional turbine cascades. The viscous regions of the flow were resolved and non-slip boundary conditions were utilized on solid surfaces. The calculations were performed in a constructive 'O'-type grid which allows representation of the blade rounded trailing edge. Converged solutions were obtained in relatively few time steps (about 80-150) and comparisons for both surface pressure and heat transfer showed good agreement with data. The three-dimensional turbine cascade calculation showed many of the expected flow-field features.

Weinberg, B. C.↗

Calculated occultation profiles of Io and the hot spots

Occultations of Io by other Galilean satellites in 1985 provide a means to locate volcanic hot spots and to model their temperatures. The expected time variations in the integral reflected and emitted radiation of the occultations are computed as a function of wavelength (visual to 8.7 microns). The best current ephemerides were used to calculate the geometry of each event as viewed from earth. Visual reflectances were modeled from global mosaics of Io. Thermal emission from the hot spots was calculated from Voyager 1 IRIS observations and, for regions unobserved by IRIS, from a model based on the distribution of low-albedo features. The occultations may help determine (1) the location and temperature distribution of Loki; (2) the source(s) of excess emission in the region from long 50 deg to 200 deg and (3) the distribution of small, high-temperature sources.

Mcewen, A. S.↗

Calculations of microwave brightness temperature of rough soil surfaces

A model for simulating the brightness temperatures of soils with rough surfaces is developed. The surface emissivity of the soil media is obtained by the integration of the bistatic scattering coefficients for rough surfaces. The roughness of a soil surface is characterized by two parameters, the surface height standard deviation sigma and its horizontal correlation length l. The model calculations are compared to the measured angular variations of the polarized brightness temperatures at both 1.4 GHz and 5 GHz frequencies. A nonlinear least-squares fitting method is used to obtain the values of delta and l that best characterize the surface roughness. The effect of shadowing is incorporated by introducing a function S(theta), which represents the probability that a point on a rough surface is not shadowed by other parts of the surface. The model results for the horizontal polarization are in excellent agreement with the data. However, for the vertical polarization, some discrepancies exist between the calculations and data, particularly at the 1.4 GHz frequency. Possible causes of the discrepancy are discussed.

Mo, T.↗

Monochromatic calculations of atmospheric radiative transfer due to molecular line absorption

Sensitivity studies related to the effects of line cutoff, spectral resolution, and temperature and pressure interpolations in radiative transfer have been performed so that a data set of absorption coefficients for water vapor, CO2, and O3 may be created efficiently. Results show that computations of absorption coefficients are affected only slightly by cutting a line off at a wave number 190 times the Lorentz half width from the center, or equivalently, cutting off 0.33 percent of the line intensity from the wings. To achieve a relative cooling rate error smaller than 2 percent, it is sufficient to precompute the absorption coefficient at three temperatures (210, 250, and 290 K) and 19 pressures with Delta (log 10 p) = 0.2. The absorption coefficient at other conditions can be interpolated linearly with pressure and exponentially with a quadratic in temperature. For the spectral resolution the absorption coefficients can be adequately computed at 0.01, 0.002, 0.005, and 0.025/cm intervals in the thermal water vapor, the CO2 and O3 bands, and the solar water vapor bands, respectively, which limits the error to only a few percent in the cooling and heating rates. Using the precomputed absorption coefficients, repeated monochromatic calculations of atmospheric heating/cooling rates for radiation model developments and for comparison with less detailed calculations are no longer difficult.

Chou, M.-D.↗

Continuum calculations of continental deformation in transcurrent environments

A thin viscous sheet approximation is used to investigate continental deformation near a strike-slip boundary. The vertically averaged velocity field is calculated for a medium characterized by a power law rheology with stress exponent n. Driving stresses include those applied along boundaries of the sheet and those arising from buoyancy forces related to lateral differences in crustal thickness. Exact and approximate analytic solutions for a region with a sinusoidal strike-slip boundary condition are compared with solutions for more geologically relevant boundary conditions obtained using a finite element technique. The across-strike length scale of the deformation is approximately 1/4pi x sq rt n times the dominant wavelength of the imposed strike-slip boundary condition for both the analytic and the numerical solutions; this result is consistent with length scales observed in continental regions of large-scale transcurrent faulting. An approximate, linear relationship between displacement and rotation is found that depends only on the deformation length scale and the rheology. Calculated displacements, finite rotations, and distribution of crustal thicknesses are consistent with those observed in the region of the Pacific-North America plate boundary in California.

Sonder, L. J.↗

Unsteady transonic flow calculations for wing-fuselage configurations

The development of the XTRAN3S wing fuselage capability for predicting transonic unsteady aerodynamic loads for aeroplastic applications is described. The modifications of the XTRAN3S code, and XTRAN3S fuselage modeling based on an alternating-direction implicit finite-difference algorithm for solving a modified transonic small-disturbance equation, are analyzed. The XTRAN3S transonic code is applied to the steady and unsteady calculations for three wing-fuselage configurations and a transport-type wing-fuselage model. The effects of fuselage aerodynamic interference on the transonic unsteady air load and flutter characteristics of wings are investigated. Good correlation is detected for the calculated and experimental pressure distributions of the models confirming the validity of the code.

Batina, J. T.↗

Practical method for highly accurate large-scale surface calculations

An accurate and efficient film linearized muffin-tin orbital (FLMTO) technique for surface electronic-structure calculations is presented which uses only 60-70 basis functions, as opposed to the 300 functions used in the linear augmented plane-wave method. Calculations for three different (3d and 4d) transition-metal films resulted in high quality results for five-layer slabs of Cu(001), Fe(001), and Ru(001), in addition to good results for the work functions and projected density of states. By retaining the LMTO small basis size, computer time and memory are reduced, making practical the study of systems with a larger number of atoms in the two-dimensional unit cell.

Fernando, G. W.↗

Fast approach for calculating film thicknesses and pressures in elastohydrodynamically lubricated contacts at high loads

The film thicknesses and pressures in elastohydrodynamically lubricated contacts have been calculated for a line contact by using an improved version of Okamura's approach. The new approach allows for lubricant compressibility, the use of Roelands' viscosity, a general mesh (nonconstant step), and accurate calculations of the elastic deformation. The new approach is described, and the effects on film thickness, pressure, and pressure spike of each of the improvements are discussed. Successful runs have been obtained at high pressure (to 4.8 GPa) with low CPU times.

Houpert, L. G.↗

Thermal-capillary analysis of small-scale floating zones Steady-state calculations

Galerkin finite element analysis of a thermal-capillary model of the floating zone crystal growth process is used to predict the dependence of molten zone shape on operating conditions for the growth of small silicon boules. The model accounts for conduction-dominated heat transport in the melt, feed rod and growing crystal and for radiation between these phases, the ambient and a heater. Surface tension acting on the shape of the melt/gas meniscus counteracts gravity to set the shape of the molten zone. The maximum diameter of the growing crystal is set by the dewetting of the melt from the feed rod when the crystal radius is large. Calculations with small Bond number show the increased zone lengths possible for growth in a microgravity environment. The sensitivity of the method to the shape and intensity of the applied heating distribution is demonstrated. The calculations are compared with experimental observations.

Duranceau, J. L.↗