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,585 records · Page 88

Recursive calculation of geared robot manipulator dynamics

A recursive formulation is presented for the calculation of the inverse and forward dynamics of rigid robot manipulators with gear systems on each joint. The complete effects of the gear ratios and the gyroscopic effects of the spinning motor/gear are included in the recursive formulation. The forward dynamics solution recursively calculates the joint accelerations when given motor torques, and the number of computations grows linearly with the number of links.

Murphy, Steve H.↗

Impeller response calculation due to complex pressure loading

An analysis technique is described to calculate the harmonic response of the first-stage impeller of the Space Shuttle Main Engine high pressure fuel turbopump. The excitation is a complex pressure loading at various locations on the impeller blades. The pressure loading was predicted using computational fluid dynamic techniques and was given as a Fourier series at 48 different locations on each of the three impeller blade types. The analysis consisted of mapping the pressures onto a three-dimensional finite element model, then converting these pressure loads to complex nodal force vectors at specified equal time intervals. The resulting vectors were then converted to modal force time histories that were then transformed into the frequency domain where the frequency response was calculated directly. This technique resulted in an improvement to the previously used direct integration technique and in a substantial analysis cost reduction.

Wellstein, Carl↗

A simple parameterization of ozone infrared absorption for atmospheric heating rate calculations

A simple parameterization of ozone absorption in the 9.6-micron region which is suitable for two- and three-dimensional stratospheric and tropospheric models is presented. The band is divided into two parts, a brand center region and a band wing region, grouping together regions for which the temperature dependence of absorption is similar. Each of the two regions is modeled with a function having the form of the Goody random model, with pressure and temperature dependent band parameters chosen by empirically fitting line-by-line equivalent widths for pressures between 0.25 and 1000 mbar and ozone absorber amounts between 1.0 x 10 to the -7th and 1.0 cm atm. The model has been applied to calculations of atmospheric heating rates using an absorber amount weighted mean pressure and temperature along the inhomogeneous paths necessary for flux computations. In the stratosphere, maximum errors in the heating rates relative to line-by-line calculations are 0.1 K/d, or 5 percent of the peak cooling at the stratopause. In the troposphere the errors are at most 0.005 K/d.

Rosenfield, Joan E.↗

Infrared cooling rate calculations in operational general circulation models - Comparisons with benchmark computations

The performance of several parameterized models is described with respect to numerical prediction and climate research at GFDL, NCAR, and GISS. The radiation codes of the models were compared to benchmark calculations and other codes for the intercomparison of radiation codes in climate models (ICRCCM). Cooling rates and fluxes calculated from the models are examined in terms of their application to established general circulation models (GCMs) from the three research institutions. The newest radiation parameterization techniques show the most significant agreement with the benchmark line-by-line (LBL) results. The LBL cooling rates correspond to cooling rate profiles from the models, but the parameterization of the water vapor continuum demonstrates uncertain results. These uncertainties affect the understanding of some lower tropospheric cooling, and therefore more accurate parameterization of the water vapor continuum, as well as the weaker absorption bands of CO2 and O3 is recommended.

Kiehl, J. T.↗

Calculations concerning the reaction C + H3(+) - CH(+) + H2

The ion-molecule reaction of C + H3(+) - CH(+) + H2 has not been studied in the laboratory but is thought to be important in the gas phase synthesis of organic molecules in dense interstellar clouds. This reaction was studied theoretically by performing quantum chemical ab initio calculations on the potential surface. It is found that there is no activation barrier to the reaction and that it proceeds smoothly to the first excited electronic state of CH(+). The rate coefficient as a function of temperature can then be estimated using the proper long-range potentials. The rate coefficient at 10 K is calculated to be 2.9 x 10 to the -9th cu cm/sec.

Talbi, D.↗

Euler calculations of unsteady transonic flow in cascades

In the present paper, Euler calculations of unsteady transonic flow in cascades are presented. A finite volume scheme is used to discretize the equations, which are implemented on a blade-fitted deformable mesh. The space-discretized equations are integrated forward in time using a multistage Runge-Kutta scheme. Adaptive dissipation terms of the type proposed by Jameson and Baker are added to capture shocks and to suppress nonphysical oscillations. Phase-shifted boundary conditions are used to reduce the computational domain to a single reference passage. No assumptions of small amplitudes or small flow deflections are made. Thus, the present code makes it possible to carry out aeroelastic calculations for cases where the shock strengths and oscillation amplitudes exceed the inherent limitations of potential flow codes.

Bendiksen, Oddvar O.↗

A row-by-row off-design performance calculation method for turbines

The turbine component of a gas turbine engine is frequently subjected to extreme operation conditions associated with significant changes in mass flow, turbine inlet temperature, pressure and rotational speed. These off-design operation conditions significantly affect the flow deflection within the turbine stage, which consists of individual stator and rotor rows. As a result, the stage parameters representing the velocity diagram will change and affect the efficiency and performance of the stage and, thus, the turbine. A row-by-row calculation method is presented for predicting the performance behavior of turbines under extreme off-design conditions. The method is applied to a multistage turbine for which the off-design performance is calculated and compared with the measurement.

Schobeiri, T.↗

Monte-Carlo calculations of Ge detector escape-peak efficiencies

The authors present calculated first and second escape peak efficiencies for a variety of n-type coaxial Ge detectors using the Monte Carlo method. The results are given for right-circular cylindrical geometries of size ranging from 3-cm diameter x 3-cm thick to 8-cm diameter x 8-cm thick. For completeness, the intrinsic interaction and photopeak efficiencies for energies ranging from 20 keV to 20 MeV are listed. The calculations have been tested against experimental measurement and found to be in agreement to within the quoted errors (typically 10 percent). It is found that the ratio of the second escape peak efficiency to the first escape peak efficiency, R, is approximately independent of energy for small detector volumes but has a weak power law energy dependence for large volumes. It is shown that R also varies as an inverse power law of the active volume.

Owens, Alan↗

Simulation of three-dimensional liquid sloshing flows using a strongly implicit calculation procedure

A coupled strongly implicit solution strategy for unsteady three-dimensional free surface flows has been developed based on an artificial compressibility formulation for the incompressible Navier-Stokes equations. A pseudotime term has been used in the continuity equation to permit time accurate calculations to be achieved. The scheme appears capable of tracking the free surface reasonably accurately inside a partially-filled spherical container undergoing a general rotating motion characteristic of that experienced by a spin-stabilized satellite. Five different free surface calculations have been presented. Some of the results exhibit an interesting Reynolds number dependent oscillatory behavior which is believed to be physical although no experimental results appears to be available for verification to date.

Chen, K.-H.↗

Two-dimensional model calculation of fluorine-containing reservoir species in the stratosphere

Two-dimensional model calculations have been carried out of the distributions of the fluorine-containing reservoir species HF, CF2O, and CFClO. HF constitutes the largest fluorine reservoir in the stratosphere, but CF2O also makes an important contribution to the inorganic fluorine budget. CFClO amounts are most important in the tropical lower stratosphere. HF amounts increase with altitude throughout the stratosphere, while those of CF2O and CFClO fall off above their mixing ratio peaks due to photolysis. The model is in good qualitative agreement with observed vertical profiles of HF and CF2O but tends to underestimate the total column of HF. The calculated CFClO distribution is in good agreement with the very limited data. The disagreement in the HF columns is likely due to small inaccuracies in the model's treatment of lower stratospheric photolysis of chlorofluorocarbons. The model results support the suggestion that CF2O may be heterogeneously converted to HF on the surface of polar stratospheric cloud particles. The model results also suggest that the quantum yield for photolysis of CF2O is near unity.

Kaye, Jack A.↗

Infrared measurements of HF and HCl total column abundances above Kitt Peak, 1977-1990 - Seasonal cycles, long-term increases, and comparisons with model calculations

The paper deals with hydrogen-chloride and hydrogen-fluoride total column measurements, their estimated long-term rates of increase, seasonal cycles, and variability, deduced from the analysis of a set of high-resolution infrared solar spectra recorded at Kitt Peak. The Kitt Peak observations and methods of analysis are described, the results are presented and compared with previously reported measurements and trends. The data is analyzed by using a multilayer nonlinear least-squares spectral fitting procedure and a consistent set of spectroscopic line parameters. Model-calculated hydrogen-chloride and hydrogen-fluoride total columns obtained with a two-dimensional model are discussed, and the model results are compared with the measured hydrogen-chloride and hydrogen-fluoride total columns, seasonal cycles, and trends. It is pointed out that the observed trends of both molecules are in satisfactory agreement with the model results calculated from emission histories and photooxidation rates for the source molecules.

Rinsland, C. P.↗

Sensitivity of calculated odd nitrogen distributions to the diabatic wind fields

This paper compares atmospheric total odd nitrogen and ozone computed with two different advective wind fields, one using climatological averages of ozone and temperature to obtain monthly averaged horizontal and vertical winds and the other using measurements from the LIMS instrument. Calculations using the former data show stronger poleward and downward motion in the winter season compared to those using the LIMS data. This leads to NO(y) mixing ratios in the lower stratosphere that are about 20 percent larger in the polar regions of both hemispheres and about 40 percent higher in the equatorial region for climatological transport fields compared to those derived from LIMS data. Consequently, the NO(y) distributions calculated with the LIMS advective field show worse agreement with the NO(y) values inferred from the LIMS measurements than similar results obtained with the climatological wind field.

Boughner, R. E.↗

Computational aspects of sensitivity calculations in linear transient structural analysis

The calculation of sensitivities of displacements, velocities, accelerations, and stresses in linear, structural, and transient response problems is studied. Several existing sensitivity calculation methods and two new methods are compared for three example problems. Approximation vectors such as vibration mode shapes are used to reduce the dimensionality of the finite model. To be computationally practical in large-order problems, the overall finite difference methods must use the approximation vectors from the original design in the analyses of the perturbed models. This was found to result in poor convergence of stress sensitivities in several cases. Two semianalytical techniques are developed to overcome this poor convergence. Both new methods result in very good convergence of the stress sensitivities; the computational cost is much less than would result if the vibration modes were recalculated and then used in an overall finite difference method.

Greene, W. H.↗

Impact cratering calculations

In the course of carrying out finite difference calculations, it was discovered that for large craters, a previously unrecognized type of crater (diameter) growth occurred which was called lip wave propagation. This type of growth is illustrated for an impact of a 1000 km (2a) silicate bolide at 12 km/sec (U) onto a silicate half-space at earth gravity (1 g). The von Misses crustal strength is 2.4 kbar. The motion at the crater lip associated with this wave type phenomena is up, outward, and then down, similar to the particle motion of a surface wave. It is shown that the crater diameter has grown d/a of approximately 25 to d/a of approximately 4 via lip propagation from Ut/a = 5.56 to 17.0 during the time when rebound occurs. A new code is being used to study partitioning of energy and momentum and cratering efficiency with self gravity for finite-sized objects rather than the previously discussed planetary half-space problems. These are important and fundamental subjects which can be addressed with smoothed particle hydrodynamic (SPH) codes. The SPH method was used to model various problems in astrophysics and planetary physics. The initial work demonstrates that the energy budget for normal and oblique impacts are distinctly different than earlier calculations for silicate projectile impact on a silicate half space. Motivated by the first striking radar images of Venus obtained by Magellan, the effect of the atmosphere on impact cratering was studied. In order the further quantify the processes of meteor break-up and trajectory scattering upon break-up, the reentry physics of meteors striking Venus' atmosphere versus that of the Earth were studied.

Ahrens, Thomas J.↗

Vegetation Canopies and Objects of Arbitrary Shapes: Computer Generation and Bidirectional Reflectance Calculations

A general computer graphics based model is presented for computer generation of objects of arbitrary shape and for calculating Bidirectional Reflectance Factor (BRF) and scattering from them, in the optical region. The computer generation uses a modified Lindemayer system (L system) approach. For rendering on a computer screen, the object is divided into polygons, and innovative computer graphics techniques are used to display the object and to calculate the scattering and reflectance from the object. The use of the technique is illustrated with scattering from canopies of simulated corn plants and from a snow covered mountain. The scattering is quantified using measures like BRF and albedo and by rendering the objects with brightness of each of the two facets of a polygon proportional to the amount of light scattered from the object in the viewer's direction.

Goel, Narenda S.↗

A loudness calculation procedure applied to shaped sonic booms

Described here is a procedure that can be used to calculate the loudness of sonic booms. The procedure is applied to a wide range of sonic booms, both classical N-waves and a variety of other shapes of booms. The loudness of N-waves is controlled by overpressure and the associated rise time. The loudness of shaped booms is highly dependent on the characteristics of the initial shock. A comparison of the calculated loudness values indicates that shaped booms may have significantly reduced loudness relative to N-waves having the same peak overpressure. This result implies that a supersonic transport designed to yield minimized sonic booms may be substantially more acceptable than an unconstrained design.

Shepherd, Kevin P.↗

Determination of Rotordynamic Coefficients for Labyrinth Seals and Application to Rotordynamic Design Calculations

In today's rotordynamic calculations, the input parameters for a finite element analysis (FEA) determine very much the reliability of eigenvalue and eigenmode predictions. While modeling of an elastic structure by means of beam elements etc. is relatively straightforward to perform and the input data for journal bearings are usually known exactly enough, the determination of stiffness and damping for labyrinth seals is still the subject of many investigations. Therefore, the rotordynamic influence of labyrinths is often not included in FEA for rotating machinery because of a lack of computer programs to calculate these parameters. This circumstance can give rise to severe vibration problems especially for high performance turbines or compressors, resulting in remarkable economic losses. The forces generated in labyrinths can be described for small motions around the seal center with a linearized force-motion relationship. Several years ago, we started with the development of computer codes for the determination of rotordynamic seal coefficients. Our different approaches to evaluate the dynamic fluid forces generated by turbulent, compressible seal flow are introduced.

Weiser, P.↗

Finite difference time domain calculation of transients in antennas with nonlinear loads

In this paper transient fields for antennas with more general geometries are calculated directly using Finite Difference Time Domain methods. In each FDTD cell which contains a nonlinear load, a nonlinear equation is solved at each time step. As a test case the transient current in a long dipole antenna with a nonlinear load excited by a pulsed plane wave is computed using this approach. The results agree well with both calculated and measured results previously published. The approach given here extends the applicability of the FDTD method to problems involving scattering from targets including nonlinear loads and materials, and to coupling between antennas containing nonlinear loads. It may also be extended to propagation through nonlinear materials.

Luebbers, Raymond J.↗