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 775 records · Page 43

All sources NASA NTRS DOE OSTI

NASA NTRS1980

Development and implementation of Shuttle/IUS proximity operations flight design software

The High Fidelity Relative Motion Program (HFRMP), a trajectory/attitude numerical integration program, was developed and implemented on the MPAD HP-9825A desk top computer systems. A solar and a lunar ephemeris is included in the HFRMP along with models of the oblate Earth, a rotating atmosphere, the orbiter's OMS/RCS/DAP system, orbiter vents, rotor dynamics, and upper stage propulsion systems. Although designed primarily for the analysis of proximity operations, it is useful in other areas such as attitude/stability analysis, propulsive consumables estimation, and trajector perturbation studies. An unique identification was assigned to each of the various configurations of the HFRMP that were developed to test new techniques and algorithms are briefly described. These include the HFRMP Versions 03H, 03M, 03T, 03U, and 05D. Development of orbiter/upper stage separation techniques including flight design support for the TDRS-A and Galileo deployment flights and design of standard maneuver sequences is discussed. Also, the development and implementation of the Euler angle conversion program is briefly addressed.

Wilson, S. W.↗

NASA NTRS1981

The prediction of normal force and rolling moment coefficients for a spinning wing

Nonlinear airfoil section data for angles of attack from 0 to 180 deg were used in a small computer code to numerically integrate the section normal force coefficients along the span as a function of the local velocity and angle of attack resulting from the combined spinning and descending motion. A correction was developed to account for the radial pressure gradient in the separated, rotating flow region above the wing. This correction was found to be necessary in order to obtain agreement, both in form and magnitude, with rotary balance test data.

Mccormick, B. W.↗

NASA NTRS1981

Acoustic properties of turbofan inlets

The finite element codes were improved using Hermitian elements and numerical integration of element relations. The question of real variable versus complex variable formulation was resolved and an integrated civil engineering system was implemented on the Georgia Tech Cyber 70/74. Efforts are underway to restructure the program to obtain the most efficient use of array storage.

Zinn, B. T.↗

NASA NTRS1981

Investigation of electrodynamic stabilization and control of long orbiting tethers

An algorithm for using electric currents to control pendular oscillations induced by various perturbing forces on the Skyhook wire is considered. Transverse and vertical forces on the tether; tether instability modes and causes during retrieval by space shuttle; simple and spherical pendulum motion and vector damping; and current generation and control are discussed. A computer program for numerical integration of the in-plane and out-of-plane displacements of the tether vs time was developed for heuristic study. Some techniques for controlling instabilities during payload retrieval and methods for employing the tether for launching satellites from the space shuttle are considered. Derivations and analyses of a general nature used in all of the areas studied are included.

Arnold, D. A.↗

NASA NTRS1979

Calculation of the flow field in supersonic inlets using a bicharacteristics method with shock wave fitting

The results of weakly viscous flow analysis are presented. The flow field, including molecular transport, is computed with the aid of a bicharacteristic method. The bow shock wave and the internal shock wave are computed with the aid of a three-dimensional shock wave fitting procedure. Characteristic equations are presented, and numerical integration procedure is discussed. Here, an inverse marching scheme is employed in which the solution is obtained on space-like planes of constant x and on space curves defined by the intersections of the internal shock wave with the solid boundaries. The distance between solution planes is arrived at by the Courant-Friedrichs-Lewy stability criterion.

Vadyak, J.↗

NASA NTRS1981

Rapid inversion of limb radiance data using an emissivity growth approximation

The time-consuming nature of limb relaxation-type inversion algorithms is due primarily to the numerous integrations over an absorption band to obtain forward radiance values with which to compare measured values. A new method has been devised for the quick and accurate (0.5% error) calculation of single gas broadband (approximately 100 per cm) limb radiance. The method uses a precalculated data base consisting of homogeneous path emissivity vs mass path data for a wide range of temperature and pressure. A 50-km altitude range, 1-km resolution, constituent inversion employing this method requires under 1 sec of computational time when run on modern computer hardware. The method does not rely upon a priori statistical knowledge.

Gordley, L. L.↗

NASA NTRS1981

Periodic solutions about the collinear Lagrangian solution in the general problem of three bodies

The article describes the solutions near Lagrange's circular collinear configuration in the planar problem of three bodies with three finite masses. The article begins with a detailed review of the properties of Lagrange's collinear solution. Lagrange's quintic equation is derived and several expressions are given for the angular velocity of the rotating frame. The equations of motion are then linearized near the circular collinear solution, and the characteristic equation is also derived in detail. The different types of roots and their corresponding solutions are discussed. The special case of two equal outer masses receives special attention, as well as the special case of two small outer masses. Finally, the fundamental family of periodic solutions is extended by numerical integration all the way up to and past a binary collision orbit. The stability and the bifurcations of this family are briefly enumerated.

Broucke, R.↗

NASA NTRS1981

Visible light scatter measurements of the Advanced X-ray Astronomical Facility /AXAF/ mirror samples

NASA is studying the properties of mirror surfaces for X-ray telescopes, the data of which will be used to develop the telescope system for the Advanced X-ray Astronomical Facility. Visible light scatter measurements, using a computer controlled scanner, are made of various mirror samples to determine surface roughness. Total diffuse scatter is calculated using numerical integration techniques and used to estimate the rms surface roughness. The data measurements are then compared with X-ray scatter measurements of the same samples. A summary of the data generated is presented, along with graphs showing changes in scatter on samples before and after cleaning. Results show that very smooth surfaces can be polished on the common substrate materials (from 2 to 10 Angstroms), and nickel appears to give the lowest visible light scatter.

Griner, D. B.↗

NASA NTRS1980

Representation and material charging response of geoplasma environments

The sensitivity of the charging response to the representation of the measured environments and material properties are discussed. Single and double Maxwellian representations are compared with direct numerical integration of the observed spectra. The effect of anisotropic incident flux distribution is modeled. In addition, the effect of the high energy radiation upon bulk conductivity and hence differential charging is examined.

Stannard, P. R.↗

NASA NTRS1981

Modelling of polymer photodegradation for solar cell modules

A computer program developed to model and calculate by numerical integration the varying concentrations of chemical species formed during photooxidation of a polymeric material over time, using as input data a choice set of elementary reactions, corresponding rate constants and a convenient set of starting conditions is evaluated. Attempts were made to validate the proposed mechanism by experimentally monitoring the photooxidation products of small liquid alkane which are useful starting models for ethylene segments of polymers like EVA. The model system proved in appropriate for the intended purposes. Another validation model is recommended.

Somersall, A. C.↗

NASA NTRS1980

An investigation of dynamic-analysis methods for variable-geometry structures

Selected space structure configurations were reviewed in order to define dynamic analysis problems associated with variable geometry. The dynamics of a beam being constructed from a flexible base and the relocation of the completed beam by rotating the remote manipulator system about the shoulder joint were selected. Equations of motion were formulated in physical coordinates for both of these problems, and FORTRAN programs were developed to generate solutions by numerically integrating the equations. These solutions served as a standard of comparison to gauge the accuracy of approximate solution techniques that were developed and studied. Good control was achieved in both problems. Unstable control system coupling with the system flexibility did not occur. An approximate method was developed for each problem to enable the analyst to investigate variable geometry effects during a short time span using standard fixed geometry programs such as NASTRAN. The average angle and average length techniques are discussed.

Austin, F.↗

NASA NTRS1981

Computer prediction of dual reflector antenna radiation properties

A program for calculating radiation patterns for reflector antennas with either smooth analytic surfaces or with surfaces composed of a number of panels. Techniques based on the geometrical optics (GO) approach were used in tracing rays over the following regions: from a feed antenna to the first reflector surface (subreflector); from this reflector to a larger reflector surface (main reflector); and from the main reflector to a mathematical plane (aperture plane) in front of the main reflector. The equations of GO were also used to calculate the reflected field components for each ray making use of the feed radiation pattern and the parameters defining the surfaces of the two reflectors. These resulting fields form an aperture distribution which is integrated numerically to compute the radiation pattern for a specified set of angles.

Christodoulou, C.↗

NASA NTRS1981

Quasi-Langrangian models of nascent thermals

The motions in and around an isolated thermal were studied and rising motion in the core, and sinking motion on the outside were found; while the circulation resembled that of a vortex ring. In an entity cloud model, cloudy thermal is tracked, in a Lagrangian fashion, as a discrete entity; the field of motion in and around the thermal is not explicitly simulated. Field of motion cloud models, in which the equations of motion are numerically integrated on an Eulerian grid were developed. It is shown that the great potential of a hybrid cloud model can combine the simplicity of the entity models with the generality and flexibility of the field-of-motion models. A key problem to be overcome in the development of a hybrid model is the formulation of a mathematical framework within which the cloud dynamics can be represented.

Rambaldi, S.↗

NASA NTRS1982

A prescribed wake rotor inflow and flow field prediction analysis, user's manual and technical approach

A user's manual is provided which includes the technical approach for the Prescribed Wake Rotor Inflow and Flow Field Prediction Analysis. The analysis is used to provide the rotor wake induced velocities at the rotor blades for use in blade airloads and response analyses and to provide induced velocities at arbitrary field points such as at a tail surface. This analysis calculates the distribution of rotor wake induced velocities based on a prescribed wake model. Section operating conditions are prescribed from blade motion and controls determined by a separate blade response analysis. The analysis represents each blade by a segmented lifting line, and the rotor wake by discrete segmented trailing vortex filaments. Blade loading and circulation distributions are calculated based on blade element strip theory including the local induced velocity predicted by the numerical integration of the Biot-Savart Law applied to the vortex wake model.

Egolf, T. A.↗

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.↗