Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical integration”

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 901 records · Page 50

Investigation of turbulent shear flows by numerical simulation

A data base (Kim et al., 1986) generated by direct numerical integration of the Navier-Stokes equations for unsteady three-dimensional turbulent channel flow at low Reynolds numbers is applied to investigate the physical processes acting in wall-bounded shear flows. The techniques used to construct the data base are outlined; the modeling of turbulence, streaks and streamwise vortices, horseshoe and inverted-horseshoe vortices, and vortex rings is described; the conditional sampling techniques employed are discussed; and typical results are presented graphically. An organized vorticity-field structure is revealed and found to be associated with the bursting process and a local peak in the wall pressure, indicating a possible approach to reducing viscous drag.

Kim, J.↗

A piloted simulator evaluation of a ground-based 4D descent advisor algorithm

A ground-based, four-dimensional (4D) descent-advisor algorithm is under development at NASA Ames Research Center. The algorithm combines detailed aerodynamic, propulsive, and atmospheric models with an efficient numerical integration scheme to generate 4D descent advisories. This paper investigates the ability of the 4D descent advisor algorithm to provide adequate control of arrival time for aircraft not equipped with on-board 4D guidance systems. A piloted simulation was conducted to determine the precision with which the descent advisor could predict the 4D trajectories of typical straight-in descents flown by airline pilots under different wind conditions. The effects of errors in the estimation of wind and initial aircraft weight were also studied. A description of the descent advisor as well as the results of the simulation studies are presented.

Green, Steven M.↗

Wind measurements for non-uniform wind fields from spaceborne scatterometers

Radar backscattering coefficient measurements by spaceborne scatterometers are presently simulated for the case of nonuniform wind fields, by means of a detailed numerical integration of the radar equation. The winds thus estimated are then compared with a nominal field which is defined as the average wind vector over the wind cell. The simulation results obtained for the NASA scatterometer are presented for cases of random wind fields whose spectra are consistent with the Seasat scatterometer sea surface wind spectrum. When the nonuniformity is small, system noise dominates the wind error; wind error degradation is therefore small for both perfect and imperfect coregistration cases. When it is relatively large, however, the wind error degradation persistently increases for both perfect and imperfect coregistrations.

Chi, Chong-Yung↗

Estimating Geophysical Parameters From Gravity Data

ORBSIM program developed for accurate extraction of parameters of geophysical models from Doppler-radio-tracking data acquired from orbiting planetary spacecraft. Model of proposed planetary structure used in numerical integration along simulated trajectories of spacecraft around primary body. Written in FORTRAN 77.

Sjogren, William L.↗

Flap-lag-torsion aeroelastic stability of a circulation control rotor in forward flight

The aeroelastic stability of a circulation control rotor blade undergoing three degrees of motion (flap, lag, and torsion) is investigated in forward flight. Quasi-steady strip theory is used to evaluate the aerodynamics forces; and the airfoil characteristics are from data tables. The propulsive and the auxiliary power trims are calculated from vehicle and rotor equilibrium equations through the numerical integration of element forces in azimuth as well as in radial directions. The nonlinear time dependent periodic blade response is calculated using an iterative procedure based on Floquet theory. The periodic perturbation equations are solved for stability using Floquet transition matrix theory. The effects of several parameters on blade stability are examined, including advance ratio, collective pitch, thrust level, shaft tilt, structural stiffnesses variation, and propulsive and auxiliary power trims.

Chopra, Inderjit↗

Trajectory characteristics and heating of hypervelocity projectiles having large ballistic coefficients

A simple, approximate equation describing the velocity-density relationship (or velocity-altitude) has been derived from the flight of large ballistic coefficient projectiles launched at high speeds. The calculations obtained by using the approximate equation compared well with results for numerical integrations of the exact equations of motion. The flightpath equation was used to parametrically calculate maximum body decelerations and stagnation pressures for initial velocities from 2 to 6 km/s. Expressions were derived for the stagnation-point convective heating rates and total heat loads. The stagnation-point heating was parametrically calculated for a nonablating wall and an ablating carbon surface. Although the heating rates were very high, the pulse decayed quickly. The total nose-region heat shield weight was conservatively estimated to be only about 1 percent of the body mass.

Tauber, Michael E.↗

Computer architecture for efficient algorithmic executions in real-time systems: New technology for avionics systems and advanced space vehicles

Improvements and advances in the development of computer architecture now provide innovative technology for the recasting of traditional sequential solutions into high-performance, low-cost, parallel system to increase system performance. Research conducted in development of specialized computer architecture for the algorithmic execution of an avionics system, guidance and control problem in real time is described. A comprehensive treatment of both the hardware and software structures of a customized computer which performs real-time computation of guidance commands with updated estimates of target motion and time-to-go is presented. An optimal, real-time allocation algorithm was developed which maps the algorithmic tasks onto the processing elements. This allocation is based on the critical path analysis. The final stage is the design and development of the hardware structures suitable for the efficient execution of the allocated task graph. The processing element is designed for rapid execution of the allocated tasks. Fault tolerance is a key feature of the overall architecture. Parallel numerical integration techniques, tasks definitions, and allocation algorithms are discussed. The parallel implementation is analytically verified and the experimental results are presented. The design of the data-driven computer architecture, customized for the execution of the particular algorithm, is discussed.

Carroll, Chester C.↗

The viscoplasticity theory based on overstress applied to the modeling of a nickel base superalloy at 815 C

Short term strain rate change, creep and relaxation tests were performed in an MTS computer controlled servohydraulic testing machine. Aging and recovery were found to be insignificant for test times not exceeding 30 hrs. The material functions and constants of the theory were identified from results of strain rate change tests. Numerical integration of the theory for relaxation and creep tests showed good predictive capabilities of the viscoplasticity theory based on overstress.

Krempl, E.↗

An error floor in tone calibrated transmission

Use of a low level pilot tone has been shown to eliminate the error floor in fading channels. This paper demonstrates that non-idealities in the receiver's pilot tone filter cause reappearance of the error floor. It also presents the bit error rate (BER) in closed form, in contrast to the multidimensional numerical integration of previous work.

Cavers, James K.↗

A study of the dynamics of the Intertropical Convergence Zone (ITCZ) in a symmetric atmosphere-ocean model

A numerical model of the circulation of a coupled axisymmetric atmosphere-ocean system was constructed to investigate the physical factors governing the location and intensity of the Intertropical Convergence Zone (ITCZ) over oceans and over land. The results of several numerical integrations are presented to illustrate the interaction of the individual atmospheric and oceanic circulations. It is shown that the ITCA cannot be located at the equator because the atmosphere-ocean system is unstable for lateral displacements of the ITCA from an equilibrium position at the equator.

Charney, J. G.↗

High resolution absorption cross-sections and band oscillator strengths of the Schumann-Runge bands of oxygen at 79 K

Cross sections of O2 at 79 K have been obtained from photoabsorption measurements at various pressures throughout the wavelength region 179.3-198.0 nm with a 6.65-m photoelectric scanning spectrometer equipped with a 2400-lines/mm grating and having an instrumental width (FWHM) of 0.0013 nm. The measured absorption cross sections of the Schumann-Runge bands (12,0) through (2,0) are independent of the instrumental width. The measured cross-sections are presented graphically here and are available at wavenumber intervals of about 0.1/cm as numerical compilations stored on magnetic tape from the National Space Science Data Center, NASA/Goddard. Band oscillator strengths of these bands have been determined by direct numerical integration of the measured cross sections.

Yoshino, K.↗

Temporal and spatial adaptive algorithm for reacting flows

A numerical integration scheme for quasi-one-dimensional unsteady flows with finite-rate chemistry is developed and demonstrated. The governing Euler and species-conservation equations are derived; the integration method and the spatial and temporal grid embedding techniques are explained; and results for sample problems involving stream-tube flow with one dissociating gas, shock-tube flow with one dissociating gas, and diverging-channel flow with multiple reactions are presented in extensive graphs and briefly characterized. Accuracy comparable to that of globally fine grid solutions is obtained with significant savings in CPU time.

Pervaiz, Mehtab M.↗

Particle paths and phase plane for time-dependent similarity solutions of the one-dimensional Vlasov-Maxwell equations

The phase trajectories of particles in a plasma described by the one-dimensional Vlasov-Maxwell equations are determined qualitatively, analyzing exact general similarity solutions for the cases of temporally damped and growing (sinusoidal or localized) electric fields. The results of numerical integration in both untransformed and Lie-group point-transformed coordinates are presented in extensive graphs and characterized in detail. The implications of the present analysis for the stability of BGK equilibria are explored, and the existence of nonlinear solutions arbitrarily close to and significantly different from the BGK solutions is demonstrated.

Roberts, Dana Aaron↗

Outcomes of planetary close encounters - A systematic comparison of methodologies

Several methods for estimating the outcomes of close planetary encounters are compared on the basis of the numerical integration of a range of encounter types. An attempt is made to lay the foundation for the development of predictive rules concerning the encounter outcomes applicable to the refinement of the statistical mechanics that apply to planet-formation and similar problems concerning planetary swarms. Attention is given to Oepik's (1976) formulation of the two-body approximation, whose predicted motion differs from the correct three-body behavior.

Greenberg, Richard↗

The rate parameters for coupled vibration-dissociation in a generalized SSH approximation

A theoretical study of vibrational excitations and dissociations of nitrogen undergoing a nonequilibrium relaxation process upon heating and cooling is reported. The rate coefficients for collisional induced vibrational transitions and transitions from a bound vibrational state into a dissociative state have been calculated using an extension of the theory originally proposed by Schwarz (SSH) et al. (1952). High-lying vibrational states and dissociative states were explicitly included but rotational energy transfer was neglected. The transition probabilities calculated from the SSH theory were fed into the master equation, which was integrated numerically to determine the population distribution of the vibrational states as well as bulk thermodynamic properties. The results show that: (1) the transition rates have a minimum near the middle of the bound vibrational levels, causing a bottleneck in the vibrational relaxation and dissociation rates; (2) high vibrational states are always in equilibrium with the dissociative state; (3) for the heating case, only the low vibrational states relax according to the Landau-Teller theory; (4) for the cooling case, vibrational relaxation cannot be described by a rate equation; (5) Park's (1985, 1988) two-temperature model is approximately valid; and (6) the average vibrational energy removed in dissociation is about 30 percent of the dissociation energy.

Sharma, Surendra P.↗

Developments in the simulation of a geopotential research mission

An essential element of any satellite system that will be used to recover information about high degree and order terms in the geopotential model of the earth is one or more low altitude (about 160 km) satellites equipped with a drag compensation mechanism. To study the effects of various error sources and to test the new theoretical and numerical techniques that will be developed for such a mission, two simulated scenarios have been used with a geopotential model complete to degree and order 360 and Encke's method to numerically integrate the equations of motion. The two scenarios include the low-low dual satellite system with integrated, one-way Doppler measurements and the single satellite system with gradiometer measurements. The simulations include a reference orbit which is assumed to be available from conventional tracking systems.

Schutz, B. E.↗

The nonlinear development of Goertler vortices in growing boundary layers

The development of Goertler vortices in boundary layers over curved walls in the nonlinear regime is investigated. The growth of the boundary layer makes a parallel flow analysis impossible except in the high wavenumber regime so in general the instability equations must be integrated numerically. Here the spanwise dependence of the basic flow is described using Fourier series expansion while the normal and streamwise variations are taken into account using finite differences. The calculations suggest that a given disturbance imposed at some position along the wall will eventually reach a local equilibrium state essentially independent of the initial conditions. In fact, the equilibrium state reached is qualitatively similar to the large amplitude high wave-number solution described asymptotically by Hall (1982). In general, it is found that the nonlinear interactions are dominated by a mean field type of interaction between the mean flow and the fundamental. Thus, even though higher harmonics of the fundamental are necessarily generated, most of the disturbance energy is confined to the mean flow correction and the fundamental. A major result of the calculations is finding that the downstream velocity field develops a strongly inflectional character as the flow moves downstream. The latter result suggests that the major effect of Goertler vortices on boundary layers of practical importance might be to make them highly receptive to rapidly growing Rayleigh modes of instability.

Hall, Philip↗

Comprehensive Silicon-Solar-Cell Program

Comprehensive Silicon Solar Cell Computer Modeling (SICELL) program simulates silicon solar cell. Predicts device parameters as efficiency, voltage-vs.-current characteristic fill factor, and temperature coefficients of parameters. Technique used similar to numerical-integration methods, but commonality described by use of mesh-point field. Validation studies show accuracies of simulations range from 0.08 percent to 3.6 percent for 27 experimental data points over temperature range of 300K to 421K. Results obtained by use of 10 mesh points in n- and p-type regions and for two iterations. SICELL runs interactively on VAX computer under VMS and written in VAX/VMS FORTRAN 77.

Lamorte, Michael F.↗