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.

SourceResearch connectionFrom yearThrough year

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

At least 649 records · Page 36

All sources NASA NTRS DOE OSTI

NASA NTRS1982

A Numerical simulation of transition in plane channel flow

A numerical simulation of the final stages of transition to turbulence in plane channel flow at a Reynolds number of 7500 is described. Three dimensional, incompressible Navier-Stokes equations are numerically integrated to obtain the time evolution of two and three dimensional finite amplitude disturbances. Computations are performed on the CYBER-203 vector processor for a 32 by 33 by 32 grid. Solutions indicate the existence of structures similar to those observed in the laboratory and which are characteristic of various stages of transition that lead to final breakdown. Details of the resulting flow field after breakdown indicate the evolution of streak-like formations found in turbulent flows. Although the flow field does approach a steady state (turbulent channel flow), implementation of subgrid-scale terms are necessary to obtain proper turbulent statistics.

Goglia, G.↗

NASA NTRS1981

An autonomous Nimbus Doppler positioning system

A single-pass Doppler positioning system was developed at the GSFC in support of the Nimbus-6 mission. The system was designed to satisfy the following requirements: compatibility with the PDP 11/70 minicomputer; single-pass recovery of up to 250 positions with a two-sigma accuracy of 5 km; high reliability with a minimum of human interaction; and near real-time responsiveness. The system consists of a numerical integrator which includes only the J2 term (earth flattening) of the earth's gravity field, an editing routine, a first guess algorithm, a least squares position recovery program, and a routine for generating a 95% confidence circle based on a computer covariance matrix. An analysis and quantification of the following error sources that limit system accuracy are presented: orbit error, radial aliasing, computational difficulties with overhead passes, oscillator drift, and ambiguity problems. Results of system reliability tests and sample outputs are included.

Garza-Robles, R.↗

NASA NTRS1981

The long-term motion of comet Halley

The orbital motion of comet Halley is numerically integrated back to 1404 BC. Starting with an orbit based on the 1759, 1682, and 1607 observations of the comet, the integration was run back in time with full planetary perturbations and nongravitational forces taken into account at each 0.5 day time-step. Small empirical corrections were made to the computed perihelion passage time in 837 and to the osculating orbital eccentricity in 800. In nine cases, the perihelion passage times calculated by Kiang (1971) from Chinese observations have been redetermined, and osculating orbital elements are given at each apparition from 1910 back to 1404 BC.

Yeomans, D. K.↗

NASA NTRS1981

Temporal development of composition, spectra, and anisotropies during upstream particle events

An analysis is conducted of the temporal development of particle intensity, composition, and anisotropies during diffuse particle events. The observations were made with an ultra-low energy charge analyzer (ULECA) sensor. ULECA basically consists of an electrostatic deflection system and an array of solid state detectors. The simultaneous measurement of energy per charge and total energy makes it possible to determine the charge state of a particle. Spectra are derived by an iterative technique involving the numerical integration of particle spectra with the energy response functions of the individual detectors. The presented data are explained in terms of temporal and spatial structures in the upstream region. The temporal structures during the event onset are consistent with a time-dependent first-order Fermi acceleration model, as discussed by Scholer et al. (1980).

Ipavich, F. M.↗

NASA NTRS1982

The origin of the Kirkwood gaps - A mapping for asteroidal motion near the 3/1 commensurability

A mapping of the phase space onto itself with the same low-order resonance structure as the 3/1 commensurability in the planar-elliptic restricted three-body problem is obtained. This mapping is about 1,000 times faster than the usual method of numerically integrating the averaged equations of motion. It exhibits some surprising behavior that might provide a key to the origin of the Kirkwood gaps. It is noted that a test asteroid placed in the gap may evolve for a million years with low eccentricity (less than 0.05) and then suddenly jump to large eccentricity (greater than 0.3), becoming a Mars crosser. The removal of the asteroid by a close encounter with Mars would then be possible. As a first test of this hypothesis, a distribution of 300 test asteroids in the area of the 3/1 commensurability was evolved for two million years. When the Mars crossers are removed, the distribution of initial conditions reveals a gap at the location of the 3/1 Kirkwood gap.

Wisdom, J.↗

NASA NTRS1982

The tailward magnetopause field beyond 10 RE

The geomagnetic field in the near-earth magnetotail is made up of three major contributions related to the dipole field, the magnetotail current system, and the magnetopause current system exclusive of the return currents in the magnetotail current system. The present investigation is concerned with the third contribution. The calculation of the magnetospheric field is discussed, taking into account the calculation of the field inside the magnetopause on the basis of the numerical integration of the Biot Savart integral. The representation of the field for distances greater than 10 earth radii on the dark side of the magnetosphere is considered. The field internal to the magnetopause is expressed as the negative gradient of a scalar potential which is expressed as a sum of solutions to Laplace's equation in cylindrical coordinates. The origin of coordinates is the earth.

Beard, D. B.↗

NASA NTRS1982

Nonlinear analysis of rotor-bearing systems using component mode synthesis

The method of component mode synthesis is developed to determine the forced response of nonlinear, multishaft, rotor-bearing systems. The formulation allows for simulation of system response due to blade loss, distributed unbalance, base shock, maneuver loads, and specified fixed frame forces. The motion of each rotating component of the system is described by superposing constraint modes associated with boundary coordinates and constrained precessional modes associated with internal coordinates. The precessional modes are truncated for each component and the reduced component equations are assembled with the nonlinear supports and interconnections to form a set of nonlinear system equations of reduced order. These equations are then numerically integrated to obtain the system response. A computer program, which is presently restricted to single shaft systems, has been written and results are presented for transient system response associated with blade loss dynamics with squeeze film dampers, and with interference rubs.

Nelson, H. D.↗

NASA NTRS1982

Recent advances in electromagnetic synthesis and analysis of dual-shaped reflector antennas

Dual-shaped reflectors have been used for many years. Thus, these reflectors have been used as high gain antennas on Voyagers 1 and 2. The objectives of the geometrical optics (GO) dual shaped synthesis are considered. Concerning the synthesis of dual shaped reflectors, it has been shown for circular symmetric reflectors that an exact GO solution can be found to the problem of transforming, by two reflections, any feed pattern into any aperture distribution. This problem involves solving two simultaneous nonlinear ordinary differential equations. The same approach for offset geometry leads to two simultaneous nonlinear partial differential equations. It is shown that these equations could also be integrated numerically, except that in general these equations are not total and therefore, in general, they do not have a 'smooth' solution. It is further shown that the offset partial differentials often very nearly form a total differential in many cases of practical importance.

Galindo-Israel, V.↗

NASA NTRS1982

Symmetric baroclinic instability of a Hadley cell

A symmetric baroclinic instability is examined in terms of a Boussinesq fluid contained between two horizontal plates to determine the effects of the Ekman and thermal layers. Governing equations are written for a rotating reference frame, taking into account the Rossby, Ekman, and Prandtl numbers. Equations are defined for the perturbation functions, treated as an eigenvalue problem, and a numerical integration of the full eighth order differential system is performed by a shooting technique. An instability is found to occur in the Hadley cell containing both Ekman and thermal boundary layers when the Richardson number is close to unity. If the Prandtl number is fixed the critical Richardson number decreases with an increasing Ekman number until the Ekman number reaches a certain value, at which time the fluid is stable.

Antar, B. N.↗

NASA NTRS1984

Improvements in Vibration-Analysis Technique

Two reports present details of vibration-analysis technique. Technique scheme for accurate and efficient numerical integration of coupled equations of motion of vehicle and its payload(s). Applicable to large stationary structures, land vehicles, airplanes, ships, or wherever vibration analysis done repeatedly after each of many design changes.

Engels, R. C.↗

NASA NTRS1982

Tone-excited jet: Theory and experiments

A detailed study to understand the phenomenon of broadband jet-noise amplification produced by upstream discrete-tone sound excitation has been carried out. This has been achieved by simultaneous acquisition of the acoustic, mean velocity, turbulence intensities, and instability-wave pressure data. A 5.08 cm diameter jet has been tested for this purpose under static and also flight-simulation conditions. An open-jet wind tunnel has been used to simulate the flight effects. Limited data on heated jets have also been obtained. To improve the physical understanding of the flow modifications brought about by the upstream discrete-tone excitation, ensemble-averaged schlieren photographs of the jets have also been taken. Parallel to the experimental study, a mathematical model of the processes that lead to broadband-noise amplification by upstream tones has been developed. Excitation of large-scale turbulence by upstream tones is first calculated. A model to predict the changes in small-scale turbulence is then developed. By numerically integrating the resultant set of equations, the enhanced small-scale turbulence distribution in a jet under various excitation conditions is obtained. The resulting changes in small-scale turbulence have been attributed to broadband amplification of jet noise. Excellent agreement has been found between the theory and the experiments. It has also shown that the relative velocity effects are the same for the excited and the unexcited jets.

Ahuja, K. K.↗

NASA NTRS1982

Structural dynamics payload loads estimates

Methods for the prediction of loads on large space structures are discussed. Existing approaches to the problem of loads calculation are surveyed. A full scale version of an alternate numerical integration technique to solve the response part of a load cycle is presented, and a set of short cut versions of the algorithm developed. The implementation of these techniques using the software package developed is discussed.

Engels, R. C.↗

NASA NTRS1982

A wake traverse technique for use in a 2 dimensional transonic flexible walled test section

Reported two dimensional validation data from the Transonic Self-Streamlining Wind Tunnel (TSWT) concerns model lift. The models tested provided data on their pressure distributions. This information was numerically integrated over the model surface to determine lift, pressure drag and pitching moment. However, the pressure drag is only a small component of the total drag at nominal angles of attack and cannot be used to assess the quality of flow simulation. An intrusive technique for obtaining information on the total drag of a model in TSWT is described. The technique adopted is the wake traverse method. The associated tunnel hardware and control and data reduction software are outlined and some experimental results are presented for discussion.

Wolf, S. W. D.↗

NASA NTRS1982

A mechanism for plasma waves at the harmonics of the plasma frequency foreshock boundary

A bump-on-tail unstable reduced velocity distribution, constructed from data obtained at the upstream boundary of the electron foreshock by the GSFC electron spectrometer experiment on the ISEE-1 satellite, is used as the initial plasma state for a numerical integration of the 1D-Vlasov-Maxwell system of equations. The integration is carried through the growth of the instability, beyond its saturation, and well into the stabilized plasma regime. A power spectrum computed for the electric field of the stabilized plasma is dominated by a narrow peak at the Bohm-Gross frequency of the unstable field mode but also contains significant power at the harmonics of the Bohm-Gross frequency. The harmonic power is in sharp peaks which are split into closely spaced doublets. The fundamental peak at the Bohm-Gross frequency is split into a closely spaced triplet. The mechanism for excitation of the second harmonic is shown to be second order wave-wave coupling.

Klimas, A. J.↗

NASA NTRS1982

Study of tethered satellite active attitude control

Existing software was adapted for the study of tethered subsatellite rotational dynamics, an analytic solution for a stable configuration of a tethered subsatellite was developed, the analytic and numerical integrator (computer) solutions for this "test case' was compared in a two mass tether model program (DUMBEL), the existing multiple mass tether model (SKYHOOK) was modified to include subsatellite rotational dynamics, the analytic "test case,' was verified, and the use of the SKYHOOK rotational dynamics capability with a computer run showing the effect of a single off axis thruster on the behavior of the subsatellite was demonstrated. Subroutines for specific attitude control systems are developed and applied to the study of the behavior of the tethered subsatellite under realistic on orbit conditions. The effect of all tether "inputs,' including pendular oscillations, air drag, and electrodynamic interactions, on the dynamic behavior of the tether are included.

Colombo, G.↗

NASA NTRS1981

Unsteady boundary-layer injection

The boundary-layer equations for two-dimensional incompressible flow are integrated numerically for the flow over a flat plate and a Howarth body. Injection is introduced either impulsively or periodically along a narrow strip. Results indicate that injection perpendicular to the wall is transmitted instantly across the boundary layer and has little effect on the velocity profile parallel to the wall. The effect is a little more noticeable for flows with adverse pressure gradients. Injection parallel to the wall results in fuller velocity profiles. Parallel and oscillatory injection appears to influence the mean. The amplitude of oscillation decreases with distance from the injection strip but further downstream it increases again in a manner reminiscent of an unstable process.

Telionis, D. P.↗

NASA NTRS1982

Transmission of sound across a vortex layer enclosing a cylindrical column of jet

An approximate solution to the problem of transmission of sound across a cylindrical vortex was obtained. Results are considerably different from the plane vortex sheet case because of the added role played by the curvature of the jet. In comparison with the plane case, the specularly transmitted waves are more complex and require some numerical integration. Resonance waves are identically predicted for M 2, but there is also a wave field whose modified effect appears to extend the region of resonance just as the instability waves cover a region in space and time. The instability waves are predicted to exist for all Mach numbers but vanish for wavelengths that are large compared to the jet radius. The region of propagation is similarly wavelength dependent.

Luh, R.↗

NASA NTRS1983

Experimental results of a deflected thrust V/STOL nozzle research program

Four deflected thrust nozzle concepts, designed to operate at the low pressure ratio typical of high bypass-ratio turbofan engines for medium speed (subsonic) V/STOL aircraft, were studied. Maps of overall performance characteristics and exit velocity distributions are used to highlight similarities and differences between the four concepts. Analytically determined secondary flows at the exit of a 90 deg circular pipe bend are compared with the experimental results from the more complex three dimensional geometries. The relative impact of total-pressure losses and secondary flows on nozzle thrust coefficient is addressed by numerical integration of exit velocity measurements.

Burstadt, P. L.↗