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 379 records · Page 21

Monte Carlo calculations of elementary particle properties

The object of this project is to calculate the masses of the elementary particles. This ambitious goal apparently is not possible using analytic methods or known approximation methods. However, it is probable that the power of a modern super computer will make at least part of the low lying mass spectrum accessible through direct numerical computation. Initial attempts by several groups at calculating this spectrum on small lattices of space time points have been very promising. Using new methods and super computers considerable progress has been made towards evaluating the mass spectrum on comparatively large lattices. Only more time and faster machines with increased storage will allow calculations of systems with guaranteed minimal boundary effects. The ideas that currently go into this calculation are outlined.

Guralnik, G. S.↗

Transonic small disturbance calculations including entropy corrections

Murman's fully conservative mixed type finite-difference operators are first modified. A special sonic point operator with an iterative damping term is introduced which helps the convergence and does not affect the spatial conservative differences. Reliable calculations with second order supersonic schemes are obtained using two sonic operators, the regular sonic point operator followed by a first order supersonic scheme. Also, shock point operator is shown to be equivalent to fitting a locally normal shock terminating the supersonic region. The potential calculations are then modified to account for the non-isentropic jump conditions using a simple shock fitting procedure based on Prandtl relation. The entropy increase across the shock is calculated in terms of the Mach number upstream of the shock and the effect of the generated vorticity is estimated via Crocco relation. Different examples are calculated and extensions to the full potential equation are discussed.

Hafez, M.↗

Comparison of measured and calculated low latitude ionospheric properties

Measurements of ionospheric parameters above Arecibo, Puerto Rico, have been compared with a computer simulation for a variety of conditions. Agreement was found between the measured and calculated electron concentration during geomagnetically quiet conditions. Comparisons for more active conditions indicate a significant upward flow of ionization during the mid-afternoon. Calculated electron temperatures were found to be consistently lower than measured temperatures during the daytime. Calculated values of NmF2 and h(max) agreed with measured results except during the post-midnight period. Calculated values of the ion flux indicate a 24-hour net flow of ionization from the northern to the Southern Hemisphere amounting to 12 percent and 6 percent of the equilibrium flux tube content above 1000 km for the winter solstice and equinox cases, respectively.

Chandler, M. O.↗

Navier-Stokes calculations for unsteady three-dimensional vortical flows in unbounded domains

Finite-difference Navier-Stokes calculations for unsteady, three-dimensional, incompressible, viscous flows induced by initial vorticity distributions are presented and discussed in this paper. The initial vorticity distributions are assumed to be embedded in a flow field of infinite extent that is quiescent at infinity. These vorticity distributions are typical of vortex rings and other closed vortical tubes or structures. Such structures are important elements in fluid flows such as jets, atmospheric convection and the far-field wakes of aircraft; studies of their interaction may aid in an understanding of complex fluid flows. The calculations employ a method recently proposed by Ting to approximate the infinite-domain boundary value problem with a finite boundary computational domain, and this method is shown to yield accurate three-dimensional results for reasonable expenditures of computer time. Because of the efficiency of the boundary condition technique and the resulting Navier-Stokes code, a 16-bit minicomputer with virtual memory was capable of performing the calculations for the unsteady motion of two obliquely colliding vortex rings. The results of these calculations are presented in the paper.

Chamberlain, J. P.↗

Numerical calculation of subsonic jets in crossflow with reduced numerical diffusion

A series of calculations are reported for two, subsonic jet in crossflow geometries. The parametric variation examined are the lateral spacing of a row of jets. The first series of calculations corresponds to a widely space jet geometry, S/D = 4, and the second series corresponds to closely spaced jets, S/D = 2. The calculations are done with alternate differencing schemes to illustrate the impact of numerical diffusion. The calculated jet trajectories agreed well with experimental data in the widely spaced jet geometry, but not in the closely spaced geometry.

Claus, R. W.↗

Application of axisymmetric analogue for calculating heating in three-dimensional flows

A rapid, approximate method has been developed for calculating the heating rates on three-dimensional vehicles such as the Space Shuttle Orbiter and other advanced reentry configurations. The method is based on the axisymmetric analogue for three-dimensional boundary layers. It uses information obtained from a three-dimensional inviscid flowfield solution, such as HALIS, to calculate inviscid surface streamlines along which approximate heating rates are calculated independent of what happens along other streamlines. Three-dimensional effects are included through the metric coefficient that describes the divergence or convergence of streamlines. Boundary-layer edge properties are obtained from the inviscid flowfield solution by interpolating in the inviscid flowfield at a distance equal to the boundary-layer thickness away from the wall. This accounts, approximately, for the variable boundary-layer edge entropy. Using this method, heating calculations can be made along a typical streamline in a few seconds. This method has been used to accurately predict heating rates for simple shapes such as a spherically blunted cone and more complex shapes such as the Shuttle Orbiter for a variety of wind-tunnel and flight conditions. A unique feature of the method is its ability to accurately predict heating rates on the Shuttle Orbiter wing.

Hamilton, H. H.↗

Fast numerical calculations of EHD sliding traction forces Application to rolling bearings

Based on analytical calculations assuming isothermal elastohydrodynamic (EHD) lubrication conditions, and on curve-fitting of a thermal correction factor, a new formula is proposed to calculate the sliding traction force developed in a concentrated contact, assuming a nonlinear, viscoelastic lubricant and undirectional sliding. When sliding occurs in the rolling and transverse directions (bidimensional sliding), a simplified numerical method is outlined by which the components of the local shear stress are quickly calculated. The latter method is applied to the calculation of the shear stress and temperature distribution in the ball-raceway contacts of an angular-contact ball bearing. Because of the viscoelastic behavior of the lubricant, a nonzero lateral traction force is obtained.

Houpert, L.↗

Numerical calculation of subsonic jets in crossflow with reduced numerical diffusion

A series of calculations are reported for two, subsonic jet in crossflow geometries. The parametric variation examined are the lateral spacing of a row of jets. The first series of calculations corresponds to a widely space jet geometry, S/D = 4, and the second series corresponds to closely spaced jets, S/D = 2. The calculations are done with alternate differencing schemes to illustrate the impact of numerical diffusion. The calculated jet trajectories agreed well with experimental data in the widely spaced jet geometry, but not in the closely spaced geometry.

Claus, R. W.↗

Computer program for calculating flow parameters and power requirements for cryogenic wind tunnels

A computer program has been written that performs the flow parameter calculations for cryogenic wind tunnels which use nitrogen as a test gas. The flow parameters calculated include static pressure, static temperature, compressibility factor, ratio of specific heats, dynamic viscosity, total and static density, velocity, dynamic pressure, mass-flow rate, and Reynolds number. Simplifying assumptions have been made so that the calculations of Reynolds number, as well as the other flow parameters can be made on relatively small desktop digital computers. The program, which also includes various power calculations, has been developed to the point where it has become a very useful tool for the users and possible future designers of fan-driven continuous-flow cryogenic wind tunnels.

Dress, D. A.↗

A novel method of calculating far-field patterns of large aperture antennas

A method is described for calculation of the radiation pattern of large aperture antennas. A piece-wise linear approximation of the aperture field using overlapping pyramidal basis functions allows the radiation pattern of an aperture antenna to be calculated as though it were a two-dimensional array. The calculation of radiation pattern data versus theta and phi, suitable for 3-D or contour plot algorithms, is achieved by locating the array in the yz-plane and performing a summation over the aperture field data sampled on a square grid. A FORTRAN subroutine is provided for performing radiation pattern calculations. Numerical results are included to demonstrate the accuracy and convergence of the method. These numerical results indicate that typical accuracies of + or - 0.1 dB for Directivity, + or - dB for the 1st Sidelobe Level, and + - 2dB for the 2nd Sidelobe Level can be obtained with an aperture grid of 45x45 points and requires approximately 0.02 seconds CPU time per far-field data point on a VAX 11/750 with a floating point accelerator.

Bailey, M. C.↗

Net thrust calculation sensitivity of an afterburning turbofan engine to variations in input parameters

The calculated value of net thrust of an aircraft powered by a General Electric F404-GE-400 afterburning turbofan engine was evaluated for its sensitivity to various input parameters. The effects of a 1.0-percent change in each input parameter on the calculated value of net thrust with two calculation methods are compared. This paper presents the results of these comparisons and also gives the estimated accuracy of the overall net thrust calculation as determined from the influence coefficients and estimated parameter measurement accuracies.

Hughes, D. L.↗

Numerical algorithms for transonic inviscid flow calculations

It is pointed out that the supercritical wing is one of the most important features of modern transonic aerodynamics. The design of its shock free airfoil section depends on potential flow calculations. The present paper is concerned with the development of inviscid flow simulation methods based on potential formulations, taking into account also the problem of nonuniqueness of the potential solution. Nonisentropic and nonisoenergetic models are considered, and an alternative approach using the stream function is discussed. Attention is given to transonic small disturbance calculations, calculations based on the full potential equation, iterative methods, wave drag calculations, and an alternative form of Euler equations.

Hafez, M. M.↗

Calculation of three-dimensional boundary layers on rotating turbine blades

An assessment has been made of the applicability of a three dimensional boundary layer analysis to the calculation of heat transfer and streamline flow patterns on the surfaces of both stationary and rotating turbine passages. In support of this effort, an analysis has been developed to calculate a general nonorthogonal surface coordinate system for arbitrary three dimensional surfaces and also to calculate the boundary layer edge conditions for compressible flow using the surface Euler equations and experimental pressure distributions. Using available experimental data to calibrate the method, calculations are presented for the endwall, and suction surfaces of a stationary cascade and for the pressure surface of a rotating turbine blade. The results strongly indicate that the three dimensional boundary layer analysis can give good predictions of the flow field and heat transfer on the pressure, suction, and endwall surfaces in a gas turbine passage.

Anderson, O. L.↗

A comparison of three algebraic stress closures for combustor flow calculations

A comparison is made of the performance of two locally nonequilibrium and one equilibrium algebraic stress closures in calculating combustor flows. Effects of four different pressure-strain models on these closure models are also analyzed. The results show that the pressure-strain models have a much greater influence on the calculated mean velocity and turbulence field than the algebraic stress closures, and that the best mean strain model for the pressure-strain terms is that proposed by Launder, Reece and Rodi (1975). However, the equilibrium algebraic stress closure with the Rotta return-to-isotropy model (1951) for the pressure-strain terms gives as good a correlation with measurements as when the Launder et al. mean strain model is included in the pressure-strain model. Finally, comparison of the calculations with the standard k-epsilon closure results show that the algebraic stress closures are better suited for simple turbulent flow calculations.

Nikjooy, M.↗

Calculations of electric currents in Europa

Electrical currents should flow in the Galilean satellite, Europa, because it is located in Jupiter's corotating magnetosphere. The possible magnitudes of these currents are calculated by assuming that Europa is a differentiated body consisting of an outer H2O layer and a silicate core. Two types of models are considered here: one in which the water is completely frozen and a second in which there is an intermediate liquid layer. For the transverse electric mode (eddy currents), the calculated current density in a liquid layer is approximately 10 to the -5/Am. For the transverse magnetic mode (unipolar generator), the calculated current density in the liquid is severely constrained by the ice layer to a range of only 10 to the -10 to -11th power/ Am, for a total H2O thickness of 100 km, provided that neither layer is less than 4 km thick. The current density is less for a completely frozen H2O layer. If transient cracks were to appear in the ice layer, thereby exposing liquid, the calculated current density could rise to a range of 10 to the -6 to 10 to the -5/Am, depending on layer thicknesses, which would require an exposed area of 10 to the -9 to 10 to the -8 of the Europa surface. The corresponding total current of 2.3x10 to the 5th power A could in 1 yr. electrolyze 7x10 to the 5th power kg of water (and more if the cells were in series), and thereby store up to 10 the 8th power J of energy, but it is not clear how electrolysis can take place in the absence of suitable electrodes. Electrical heating would be significant only if the ice-layer thickness were on the order of 1 m, such as might occur if an exposed liquid surface were to freeze over; the heating under this condition could hinder the thickening of the ice layer.

Colburn, D. S.↗

A VLSI Processor Engine to Support Calculation of the Vector Wave Equation

The development of a VLSI processor engine tailored to support computation of the near-field vector wave equation (VWE) is presented. Calculation of the near field is essential for evaluation of the electromagnetic fields generated by a complex radiating structure. The processor engine presented is optimized such that it takes advantage of the extensive parallelism inherent in the computational algorithm. The utilization of very large scale integration (VLSI) circuit technology and a customized reduced instruction set computer (RISC) architecture support development of a system of these processor engines, each of which is capable of performing the VWE algorithm calculation. By taking maximum advantage of the parallelism present in the calculation and dividing the problem among an array of these processor engines, the calculation can be computed in a realistic time frame.

J W Degroat↗

Mixing rates calculated from potential vorticity

Horizontal mixing coefficients K(yy) and K(yz) were computed using potential vorticity and gradients calculated from NMC data for geopotential heights and temperatures. It was found that the mid-latitude to polar mixing coefficients depend on the use of nongeostrophic winds, which yield K(yy) values which are mostly positive, whereas geostophic winds yield K(yy) values which were largely negative near the polar region. A comparison of the K(yy) values calculated from the potential vorticity fluxes and gradients with the K(yy) values calculated using the momentum balance between the flux of potential vorticity and the residual circulation has shown that the residual-circulation-balanced K(yy) values in the upper stratosphere were much larger than the values calculated by the present method. It is suggested that this difference is the result of breaking gravity waves in the upper stratosphere and lower mesosphere.

Newman, Paul A.↗

Noise produced by turbulent flow into a rotor: Theory manual for noise calculation

An analysis is presented for the calculation of noise produced by turbulent flow into a helicopter rotor. The method is based on the analysis of Amiet for the sound produced by an airfoil moving in rectilinear motion through a turbulent flow field. The rectilinear motion results are used in a quasi-steady manner to calculate the instantaneous spectrum of the rotor noise at any given rotor position; the overall spectrum is then found by averaging the instantaneous spectrum over all rotor azimuth angles. Account is taken of the fact that the rotor spends different amounts of retarded time at different rotor positions. Blade to blade correlation is included in the analysis, leading to harmonics of blade passing frequency. The spectrum of the turbulence entering the rotor is calculated by applying rapid distortion theory to an isotropic turbulence spectrum, assuming that the turbulence is stretched as it is pulled into the rotor. The inputs to the program are obtained from the atmospheric turbulence model and mean flow distortion calculation, described in another volume of this set of reports. The analytical basis is provided for a module which was incorporated in NASA's ROTONET helicopter noise prediction program.

Amiet, R. K.↗