SEARCH · Search NASA

Results for “Arithmetic”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Data reduction programs for a laser radar system

The listing and description of software routines which were used to analyze the analog data obtained from LIDAR - system are given. All routines are written in FORTRAN - IV on a HP - 1000/F minicomputer which serves as the heart of the data acquisition system for the LIDAR program. This particular system has 128 kilobytes of highspeed memory and is equipped with a Vector Instruction Set (VIS) firmware package, which is used in all the routines, to handle quick execution of different long loops. The system handles floating point arithmetic in hardware in order to enhance the speed of execution. This computer is a 2177 C/F series version of HP - 1000 RTE-IVB data acquisition computer system which is designed for real time data capture/analysis and disk/tape mass storage environment.

Badavi, F. F.↗

NASA NTRS1984

Mathematical algorithms to maximize performance in numerical weather prediction

Numerical weather prediction models, which involve the solution of non-linear partial differential equations at points on an extensive three dimensional grid, are ideally suited for processing on vector machines. It was logical therefore that the new global forecast model to be implemented at the Meteorological Office should be written in vector code for the CYBER 205. In order to achieve full efficiency and to reduce storage requirements the model used 32-bit arithmetic which was found to provide high enough precision. Unfortunately, however, the trigonometrical and logarithmic functions provided by CDC could only handle 64-bit vectors and, although written in efficient scalar code, did not take advantage of the special facilities of a vector processor. It was therefore necessary to rewrite the functions in vector code to handle both 32 and 64-bit vectors. There was also no half-precision compiler available for the Cyber 205 at that time and so the functions, like the model, had to make extensive use of the special call syntax. This made the code more difficult to write but it allowed much greater flexibility in that it became possible to access the exponent of a floating-point number independently of its coefficient. A description is given of the technique and the results which were achieved are summarized.

Foreman, A.↗

NASA NTRS1984

A highly optimized vectorized code for Monte Carlo simulations of SU(3) lattice gauge theories

New methods are introduced for improving the performance of the vectorized Monte Carlo SU(3) lattice gauge theory algorithm using the CDC CYBER 205. Structure, algorithm and programming considerations are discussed. The performance achieved for a 16(4) lattice on a 2-pipe system may be phrased in terms of the link update time or overall MFLOPS rates. For 32-bit arithmetic, it is 36.3 microsecond/link for 8 hits per iteration (40.9 microsecond for 10 hits) or 101.5 MFLOPS.

Barkai, D.↗

NASA NTRS1984

Navier-Stokes Simulation of Homogeneous Turbulence on the CYBER 205

A computer code which solves the Navier-Stokes equations for three dimensional, time-dependent, homogenous turbulence has been written for the CYBER 205. The code has options for both 64-bit and 32-bit arithmetic. With 32-bit computation, mesh sizes up to 64 (3) are contained within core of a 2 million 64-bit word memory. Computer speed timing runs were made for various vector lengths up to 6144. With this code, speeds a little over 100 Mflops have been achieved on a 2-pipe CYBER 205. Several problems encountered in the coding are discussed.

Wu, C. T.↗

NASA NTRS1984

Implementation of the Sun Position Calculation in the PDC-1 Control Microprocessor

The several computational approaches to providing the local azimuth and elevation angles of the Sun as a function of local time and then the utilization of the most appropriate method in the PDC-1 microprocessor are presented. The full algorithm, the FORTRAN form, is felt to be very useful in any kind or size of computer. It was used in the PDC-1 unit to generate efficient code for the microprocessor with its floating point arithmetic chip. The balance of the presentation consists of a brief discussion of the tracking requirements for PPDC-1, the planetary motion equations from the first to the final version, and the local azimuth-elevation geometry.

Stallkamp, J. A.↗

NASA NTRS1984

An efficient approximate factorization implicit scheme for the equations of gasdynamics

An efficient implicit finite-difference algorithm for the gas dynamic equations utilizing matrix reduction techniques is presented. A significant reduction in arithmetic operations is achieved while maintaining the same favorable stability characteristics and generality found in the Beam and Warming approximate factorization algorithm. Steady-state solutions to the conservative Euler equations in generalized coordinates are obtained for transonic flows about a NACA 0012 airfoil. The theoretical extension of the matrix reduction technique to the full Navier-Stokes equations in Cartesian coordinates is presented in detail. Linear stability, using a Fourier stability analysis, is demonstrated and discussed for the one-dimensional Euler equations. It is shown that the method offers advantages over the conventional Beam and Warming scheme and can retrofit existing Beam and Warming codes with minimal effort.

Barth, T. J.↗

NASA NTRS1983

Infrared astronomical imaging using a microcomputer data system

A small, portable data-acquisition and display system with a unique coprocessor was designed and built to collect data from infrared array detectors for astronomical and laboratory use. The system consists of a microcomputer board, two 12-bit analog-to-digital data-acquisition boards, a color-graphics driver board, and two coprocessor boards. The coprocessors are arithmetic logic units designed to achieve the high computation rates necessary for real-time infrared background cancellation. The system can collect data from infrared detectors of 1 to 256 pixel elements with a computation rate of 768 kHz. The software for the system was designed to collect data in the system's background, while performing other tasks in the system's foreground. The data for the detector arrays is stored in block fashion to facilitate variable image size.

Stafford, P. S.↗

NASA NTRS1984

An implicit form for the Osher upwind scheme

Conservative upwind schemes for the Euler equations, such as the Osher scheme, accurately resolve flow discontinuities and correctly model the physics of the problem. However, these schemes require many more arithmetic operations per integration step than simple central-difference schemes and hence result in large computing times. An implicit version of the first-order- and second-order-accurate Osher schemes in two spatial dimensions and generalized coordinates is developed in this study. Because implicit schemes permit the use of large integration steps, in many cases they require fewer integration steps to reach steady-state (especially in calculations on grids with widely varying mesh-cell sizes). The implicit scheme developed in this study accelerated convergence speeds by almost an order of magnitude in the problems considered. Test cases include quasi-one-dimensional nozzle flow and supersonic flow past a cylinder.

Rai, M. M.↗

NASA NTRS1984

An investigation into the probabilistic combination of quasi-static and random accelerations

The development of design load factors for aerospace and aircraft components and experiment support structures, which are subject to a simultaneous vehicle dynamic vibration (quasi-static) and acoustically generated random vibration, require the selection of a combination methodology. Typically, the procedure is to define the quasi-static and the random generated response separately, and arithmetically add or root sum square to get combined accelerations. Since the combination of a probabilistic and a deterministic function yield a probabilistic function, a viable alternate approach would be to determine the characteristics of the combined acceleration probability density function and select an appropriate percentile level for the combined acceleration. The following paper develops this mechanism and provides graphical data to select combined accelerations for most popular percentile levels.

Schock, R. W.↗

NASA NTRS1984

Modeling algorithm execution time on processor arrays

An approach to modelling the execution time of algorithms on parallel arrays is presented. This time is expressed as a function of the number of processors and system parameters. The resulting model has been applied to a parallel implementation of the conjugate-gradient algorithm on NASA's FEM. Results of experiments performed to compare the model predictions against actual behavior show that the floating-point arithmetic, communication, and synchronization components of the parallel algorithm execution time were correctly modelled. The results also show that the overhead caused by the interaction of the system software and the actual parallel hardware must be reflected in the model parameters. The model has been used to predict the performance of the conjugate gradient algorithm on a given problem as the number of processors and machine characteristics varied.

Adams, L. M.↗

NASA NTRS1985

Compensating for Range Migration in Synthetic-Aperture Radar

Correlation algorithm sharpens radar imagery. Azimuthal correlation of point-target response performed in frequency domain by fast Fourier transforms relatively efficient in arithmetic operations. Frequencydomain azimuthal response has simple form that enables correction for range migration made simultaneously and over wide area, thereby providing many picture elements.

Jin, M. Y.↗

NASA NTRS1984

Flight tests of the Digital Integrated Automatic Landing System (DIALS)

The design, development, implementation and flight tests of the Digital Integrated Automatic Landing System (DIALS) are discussed. The system was implemented and flight tested on the Transport Systems Research Vehicle (TSRV), a Boeing 737-100. The design uses modern optimal control methods. The direct digital design obtained uses a 10 Hz rate for the sampling of sensors and the control commands. The basic structure of the control law consists of a steady state Kalman filter followed by a control gain matrix. The sensor information used includes Microwave Landing System (MLS) position, attitude, calibrated airspeed, and body accelerations. The phases of the final approach considered are localized and steep glideslope capture (which may be performed simultaneously or independently), localizer and glideslope track, crab/decrab, and flare to touchdown. The system can capture, track, and flare from conventional, as well as steep, glideslopes ranging from 2.5 deg to 5.5 deg. All of the modes of the control law including the Kalman filters were implemented on the TSRV flight computers which use fixed point arithmetic with 16 bit words. The implementation considerations are described as well as an analysis of the flight test results.

Halyo, N.↗

NASA NTRS1985

Effect of counterface material type and its topography on the tribological properties of polyimide composites

Graphite fiber reinforced polyimide composite pins were slid against seven different counterfaces to determine the effect of material type on the tribological properties of polymer composites. In addition, the effect of sliding a new pin on a pre-established transfer film was investigated. The results indicated that almost a five order of magnitude difference in composite wear rate can occur just by varying the counterface material. An attempt to make all surfaces as smooth as possible was made, but due to differences in material composition this was not possible and a range of surface roughnesses were obtained. The results indicate that the smoother the surface, the lower the composite wear rate; but that small protrusions (not discernible with arithmetic surface roughness measurements) can markedly increase wear rates. A pre-established transfer film improved both run in and steady state wear rates.

Fusaro, R. L.↗

NASA NTRS1985

The BLAZE language: A parallel language for scientific programming

A Pascal-like scientific programming language, Blaze, is described. Blaze contains array arithmetic, forall loops, and APL-style accumulation operators, which allow natural expression of fine grained parallelism. It also employs an applicative or functional procedure invocation mechanism, which makes it easy for compilers to extract coarse grained parallelism using machine specific program restructuring. Thus Blaze should allow one to achieve highly parallel execution on multiprocessor architectures, while still providing the user with onceptually sequential control flow. A central goal in the design of Blaze is portability across a broad range of parallel architectures. The multiple levels of parallelism present in Blaze code, in principle, allow a compiler to extract the types of parallelism appropriate for the given architecture while neglecting the remainder. The features of Blaze are described and shows how this language would be used in typical scientific programming.

Mehrotra, P.↗

NASA NTRS1985

SPAR improved structural-fluid dynamic analysis capability

The results of a study whose objective was to improve the operation of the SPAR computer code by improving efficiency, user features, and documentation is presented. Additional capability was added to the SPAR arithmetic utility system, including trigonometric functions, numerical integration, interpolation, and matrix combinations. Improvements were made in the EIG processor. A processor was created to compute and store principal stresses in table-format data sets. An additional capability was developed and incorporated into the plot processor which permits plotting directly from table-format data sets. Documentation of all these features is provided in the form of updates to the SPAR users manual.

Pearson, M. L.↗

NASA NTRS1984

Perturbation analysis of the limit cycle of the free van der Pol equation

A power series expansion in the damping parameter, epsilon, of the limit cycle of the free van der Pol equation is constructed and analyzed. Coefficients in the expansion are computed in exact rational arithmetic using the symbolic manipulation system MACSYMA and using a FORTRAN program. The series is analyzed using Padeapproximants. The convergence of the series for the maximum amplitude of the limit cycle is limited by two pair of complex conjugate singularities in the complex epsilon-plane. A new expansion parameter is introduced which maps these singularities to infinity and leads to a new expansion for the amplitude which converges for all real values of epsilon. Amplitudes computed from this transformed series agree very well with reported numerical and asymptotic results. For the limit cycle itself, convergence of the series expansion is limited by three pair of complex conjugate branch point singularities. Two pair remain fixed throughout the cycle, and correspond to the singularities found in the maximum amplitude series, while the third pair moves in the epsilon-plane as a function of t from one of the fixed pairs to the other. The limit cycle series is transformed using a new expansion parameter, which leads to a new series that converges for larger values of epsilon.

Dadfar, M. B.↗

NASA NTRS1983

The middle and high latitude winter ionosphere at the Ariel 4 satellite altitude

The electron (0+) density variations over the northern and southern winter high latitude ionosphere are comprehensively analyzed using the technique of Brinton et al. (1978). Two-hour Magnetic Local Time (MLT) arithmetic means of electron densities are studied in terms of invariant magnetic latitude and in terms of magnetic activity as classified by the three-hour planetary magnetic activity index. It is found that the southern hemisphere densities are significantly lower than those in the northern hemisphere. Further, the maximum electron densities observed in the northern hemisphere are located in a MLT range symmetrical about the 14-02 MLT meridian, whereas in the southern hemisphere the maxima are observed about the noon midnight magnetic meridian. A deep localized ionization hole on the nightside of the polar cap is not observed although the polar cavity is apparent.

Tulunay, Y. K.↗

NASA NTRS1985

A fast efficient implicit scheme for the gasdynamic equations using a matrix reduction technique

An efficient implicit finite-difference algorithm for the gasdynamic equations utilizing matrix reduction techniques is presented. A significant reduction in arithmetic operations is achieved without loss of the stability characteristics generality found in the Beam and Warming approximate factorization algorithm. Steady-state solutions to the conservative Euler equations in generalized coordinates are obtained for transonic flows and used to show that the method offers computational advantages over the conventional Beam and Warming scheme. Existing Beam and Warming codes can be retrofit with minimal effort. The theoretical extension of the matrix reduction technique to the full Navier-Stokes equations in Cartesian coordinates is presented in detail. Linear stability, using a Fourier stability analysis, is demonstrated and discussed for the one-dimensional Euler equations.

Barth, T. J.↗