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 613 records · Page 34

Shuttle entry trajectory reconstruction using inflight accelerometer and gyro measurements

An error analysis has been made of a Shuttle postflight entry trajectory reconstruction process to obtain trajectory state estimation errors and to assess the impact of these errors on Shuttle aerodynamic force coefficient extraction. In this analysis, the entry trajectory is assumed to be reconstructed via numerical integration of onboard accelerometer and gyro measurements and constrained to satisfy ground-based radio tracking. The trajectory state estimation errors are calculated using a Kalman-Schmidt sequential filter assuming various measurement error models and combinations of ground-based tracking. The resultant trajectory estimation errors are analyzed in a simplified perturbation process to establish the accuracy to which postflight aerodynamic force coefficients can be determined. Results are presented which show that the principal error sources affecting the trajectory reconstruction and thus the force coefficient extraction, assuming perfect atmospheric density knowledge, are the accelerometer and gyro resolution, acceleration-sensitive gyro drifts, and the alignment uncertainties associated with integration on the Shuttle.

Compton, H. R.↗

Viscous flow analysis in mixed flow rotors

A method for analyzing the viscous flow through turbomachine rotors is presented. The field analysis is based upon the solution of the full Navier-Stokes equations over the rotor blade-to-blade stream channels. An Alternating-Direction-Implicit method is employed to carry out the necessary numerical integration of the elliptic governing equations. The flow analysis may be applied to various types of turbomachine rotors. Preliminarily, only the case of laminar flows are considered in this paper. The flow characteristics within the rotors of a mixed flow turbine and a radial bladed compressor are investigated over a wide range of operating conditions. Excellent results are obtained when compared with existing experimental data. The method of this analysis is quite general and can deal with a wide range of applications. Possible modification of the present study to deal with turbulent flow cases are also identified.

Khalil, I.↗

Present obliquity oscillations of Mars - Fourth-order accuracy in orbital e and I

A long period analysis of solar system orbital evolution, correct to fourth order in orbital eccentricities and inclinations (Bretagnon, 1974), and an improved value of the planet's moment of inertia (Reasenberg, 1977) have been incorporated in a recalculation of the obliquity oscillations of Mars. A linearized solution predicts a maximum oscillation amplitude of 13.6 deg centered on a long-term average value of 24.4 deg. A numerical integration of the obliquity for the past 10,000,000 years is also presented. Epochs of minimal oscillation like the present occur at intervals of the order of 4 m.y.

Ward, W. R.↗

Some simple properties of stellar pulsation modes

Except for the lowest harmonics, small-amplitude stellar pulsation modes possess many simple properties whose evaluation does not require numerical integration of the fourth-order equations of motion. All antinodes tend to have the same total kinetic energy except for those lying near physical or geometric boundaries. However, when kinetic energy per unit volume is considered, order-of-magnitude enhancements are seen in antinodes lying near the center of the star, and factor-of-2 enhancements occur near the polar axis. The nodes are distributed very regularly along the radius. They follow an exponential law in g-regions, and their separation is proportional to the sound travel time in p-regions. A simple graphical procedure is described for surveying the oscillation frequencies of a new stellar model. A precise condition is derived giving the division of energy between radial and angular motion. Another condition gives the fractional contribution to the velocity field of its two sources, the divergence and the curl. Certain simplifying results of weak coupling among the linear modes are briefly described.

Wolff, C. L.↗

An Improved Theoretical Aerodynamic Derivatives Computer Program for Sounding Rockets

The paper outlines a Theoretical Aerodynamic Derivatives (TAD) computer program for computing the aerodynamics of sounding rockets. TAD outputs include normal force, pitching moment and rolling moment coefficient derivatives as well as center-of-pressure locations as a function of the flight Mach number. TAD is applicable to slender finned axisymmetric vehicles at small angles of attack in subsonic and supersonic flows. TAD improvement efforts include extending Mach number regions of applicability, improving accuracy, and replacement of some numerical integration algorithms with closed-form integrations. Key equations used in TAD are summarized and typical TAD outputs are illustrated for a second-stage Tomahawk configuration.

Barrowman, J. S.↗

Diagrammatic theory of transition of pendulum like systems

Orbit-orbit and spin-orbit gravitational resonances are analyzed using the model of a rigid pendulum subject to both a time-dependent periodic torque and a constant applied torque. First, a descriptive model of passage through resonance is developed from an examination of the polynomial equation that determines the extremes of the momentum variable. From this study, a probability estimate for capture into libration is derived. Second, a lowest order solution is constructed and compared with the solution obtained from numerical integration. The steps necessary to systematically improve this solution are also discussed. Finally, the effect of a dissipative term in the pendulum equation is analyzed.

Yoder, C. F.↗

Nongravitational perturbations of long-period comets

The model of Marsden et al. (1973) is used to investigate the effect of nongravitational forces resulting from water ice sublimation on the orbits of long-period comets. The orbits of hypothetical comets with perihelion distances of 0.005 to 4 AU are integrated numerically along initially parabolic trajectories through one perihelion passage for nongravitational forces of 10 millionths to 200 millionths the solar attraction at 1 AU and lag angles of 0, 5, and 10 deg. The results indicate that the nongravitational perturbations are possibly of equal or greater importance than planetary perturbations for comets with small perihelia and that nongravitational forces could have brought sun-grazing comets to their current orbits in two or three returns. It is concluded that nongravitational forces are likely to be found only for those long-period comets that have large nongravitational accelerations to begin with and are bright enough to be discovered early and tracked over a very long arc of their orbits.

Weissman, P. R.↗

Analytical method for the effects of the asteroid belt on planetary orbits

Analytic expressions are derived for the perturbation of planetary orbits due to a thick constant-density asteroid belt. The derivations include extensions and adaptations of Plakhov's (1968) analytic expressions for the perturbations in five of the orbital elements for closed orbits around Saturn's rings. The equations of Plakhov are modified to include the effect of ring thickness, and additional equations are derived for the perturbations in the sixth orbital element, the mean anomaly. The gravitational potential and orbital perturbations are derived for the asteroid belt with and without thickness, and for a hoop approximation to the belt. The procedures are also applicable to Saturn's rings and the newly discovered rings of Uranus. The effects of the asteroid belt thickness on the gravitational potential coefficients and the orbital motions are demonstrated. Comparisons between the Mars orbital perturbations obtained by using the analytic expressions and those obtained by numerical integration are discussed. The effects of the asteroid belt on earth-based ranging to Mars are also demonstrated.

Mayo, A. P.↗

Galileo Jupiter encounter and satellite tour trajectory design

The trajectory design process for the Galilean satellite probe, expected to be launched in January 1982, is examined and its interplanetary path, emphasizing arrival dates and the Mars flyby, is presented. The Galileo spacecraft is of dual spin design, combining features of Pioneer, a spinning craft, and Voyager which was 3-axis stabilized. The design of the initial orbit from Orbiter deflection through Jupiter to the beginning of the satellite tour is discussed as well as the use of an Io gravity assist flyby to achieve Jupiter orbit insertion, and the effects of solar perturbation on the initial orbit. Detail is given to the entire satellite tour trajectory design process, starting from the requirements and constraints placed on the tour through the conic design to the final numerically integrated trajectories.

Diehl, R. E.↗

Decoupling control of a long flexible beam in orbit

The paper presents a method of control for large flexible systems using state variable feedback, with a long flexible beam given as an example. These feedback gains are selected: (1) based on the decoupling of the original coordinates and to obtain proper damping and (2) by applying the linear regulator problem to the individual modal coordinates separately. It is shown that the linear control law thus obtained are then evaluated by numerical integration of the non-linear system equations. Also included are results showing the effects (control spillover) on the uncontrolled modes when the number of controllers is less than the number of modes, and the effects of inaccurate knowledge of the control influence coefficients which lead to errors in the calculated feedback gains.

Reddy, A. S. S. R.↗

Long-term risk analysis associated with nuclear waste disposal in space

An assessment and verification of previous analytic results on the long term risk of earth reentry for hazardous payloads is presented. The two areas were studied: (1) stability of nominal, near-circular storage orbits in the regions between Venus and earth and between earth and Mars, and (2) probability of earth reentry for off-nominal planet-crossing orbits resulting from deployment system failures. In the first area, numerical integrations of the equations of motion are compared with stability predications based on secular perturbation theory. The agreement is good in terms of the heliocentric distances covered and the general behavior of the orbital history, although certain near-resonance situations can lead to difficulty. In the second area, a Monte Carlo simulation of orbital evolution is used and the results compared with Opik's analytic theory of planetary encounters and collision statistics, with data verified to within a close order-of-magnitude.

Friedlander, A. L.↗

Application of a microprocessor controlled lidar to tropospheric ozone measurements

A microprocessor controlled lidar system under construction at the NASA Goddard Space Flight Center is described and the problems in making space-based measurements of tropospheric ozone are considered. The differential absorption lidar using a dual wavelength, pulsed CO2 laser and direct detection receiver can significantly improve the existing global data base on tropospheric ozone burden. Sensitivity to tropospheric ozone can be obtained in the spaceborne version of lidar by selecting laser lines located in the wings of the target zone lines. Simulation studies using various laser line pairs in the P-branch of the CO2 9.4 micron band show that the ozone burden retrieval may be weighted to particular altitude regions. These simulation studies are based on numerical integration of differences in absorption coefficient at the two selected laser wavelengths using the AFGL absorption line parameter compilations, U.S. Standard Atmosphere ozone profile, and laser software. Simulation, data collection and reduction are performed by a microprocessor subsystem of the CO2 lidar.

Stewart, R. W.↗

Emissivity correction for interpreting thermal radiation from a terrestrial surface

A general method of accounting for emissivity in making temperature determinations of graybody surfaces from radiometric data is presented. The method differs from previous treatments in that a simple blackbody calibration and graphical approach is used rather than numerical integrations which require detailed knowledge of an instrument's spectral characteristics. Also, errors caused by approximating instrumental response with the Stephan-Boltzman law rather than with an appropriately weighted Planck integral are examined. In the 8-14 micron wavelength interval, it is shown that errors are at most on the order of 3 C for the extremes of the earth's temperature and emissivity. For more practical limits, however, errors are less than 0.5 C.

Sutherland, R. A.↗

Shell theory automated for rotational structures

Package of numerical integration programs static, buckling, vibration, and plastic analysis on thin sheels of revolution. Sheels may be subjected to distributed loads, concentrated line loads, and thermal strain. Outputs include stresses, displacement, plastic strains, and vibration and buckling results. Program aids design of aircraft bodies, spacecraft, submarines, and storage tanks. Written in FORTRAN IV for batch execution, program has been implemented on UNIVAC 1108.

Key, J.↗

Acid rain: Microphysical model

A microphysical model was used to simulate the case of a ground cloud without dilution by entrainment and without precipitation. The numerical integration techniques of the model are presented. The droplet size spectra versus time and the droplet molalities for each value of time are discussed.

Dingle, A. N.↗

A time variable model of Earth's albedo

A time variable model of Earth's albedo was prepared for use in climate studies and as an aid to the interpretation of satellite Earth radiation budget data. The features of the model include: a 10 deg latitude 10 deg longitude grid for numerical integration, surface albedo specified at 1 month intervals, calculation of zenith angle effect for surface albedo and of the additional effect of the atmosphere on the albedo. Percent cloud cover is specified for 29 different climatological cloud type regions at 8 times of the day for 12 months of the year. Cloud albedos were specified for each of the cloud climatological types. Diurnal and monthly variations of this model are described and results are compared with a model which is based on satellite measurements. A computer program was also written for use in studying the sampling effects in satellite radiation budget measurements. An example of the results of calculations with this program are compared with a previous study of the sampling effects. This program for satellite orbit characteristics is to be combined with the time-variable albedo model for further study of the sampling problem.

Bartman, F. L.↗

Finite-element numerical modeling of atmospheric turbulent boundary layer

A dynamic turbulent boundary-layer model in the neutral atmosphere is constructed, using a dynamic turbulent equation of the eddy viscosity coefficient for momentum derived from the relationship among the turbulent dissipation rate, the turbulent kinetic energy and the eddy viscosity coefficient, with aid of the turbulent second-order closure scheme. A finite-element technique was used for the numerical integration. In preliminary results, the behavior of the neutral planetary boundary layer agrees well with the available data and with the existing elaborate turbulent models, using a finite-difference scheme. The proposed dynamic formulation of the eddy viscosity coefficient for momentum is particularly attractive and can provide a viable alternative approach to study atmospheric turbulence, diffusion and air pollution.

Lee, H. N.↗

Multiclass Bayes error estimation by a feature space sampling technique

A general Gaussian M-class N-feature classification problem is defined. An algorithm is developed that requires the class statistics as its only input and computes the minimum probability of error through use of a combined analytical and numerical integration over a sequence simplifying transformations of the feature space. The results are compared with those obtained by conventional techniques applied to a 2-class 4-feature discrimination problem with results previously reported and 4-class 4-feature multispectral scanner Landsat data classified by training and testing of the available data.

Mobasseri, B. G.↗