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 757 records · Page 42

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

Phase change in liquid face seals. II - Isothermal and adiabatic bounds with real fluids

Analytical studies of phase change effects in parallel and tapered liquid face seals are presented. An isothermal and adiabatic model of low Reynolds number flow are considered by numerical integration of the descriptive equations for a real fluid, and its thermodynamic properties are calculated for each step, using a computer program for the steam tables or fluid thermodynamic properties. It was shown that for low leakage rate the isothermal model is more accurate and for high leakage rates the adiabatic model is more accurate; that both models yield the same conclusions regarding stability; and that the transient of collapse is described by the adiabatic model which predicts a catastrophic collapse and then either failure or explosive return to a larger film thickness value. Finally, it is shown that converging seals may become unstable and the mass leakage rate is reduced significantly below the all liquid value when boiling occurs.

Hughes, W. F.↗

Uncontrolled dynamics of the Skylab vehicle - Comparison of simulation, theory and flight data

When the final Skylab mission was completed in 1974, the onboard systems were used to maneuver the vehicle to a gravity gradient equilibrium attitude and were then shut down. Parked in this storage attitude, the vehicle was left undisturbed until March 1978 when contact was reestablished and the control systems were reactivated. During the intervening four years, the vehicle was subjected to subtle environmental torques which caused the vehicle to gradually progress from its initial stable state to a very dynamic condition just prior to the resumption of control. The paper deals with the dynamics of this process, including the analyses and simulations which were made to explain the motion and enable the planning of the reactivation process. The complementary roles of both numerical integration and closed-form analysis are demonstrated. Flight data obtained in the spring of 1978 are presented which support the analytical and simulation results.

Buchanan, H. J.↗

Flow over plates with suction through porous strips

This paper addresses the steady, incompressible, two-dimensional flow past a flat plate with suction through porous strips. Closed-form solutions for each flow quantity are developed in the context of linearized triple-deck theory using Fourier transforms. To demonstrate the validity of these closed-form solutions, we compare the wall shear stress and pressure coefficients and the streamwise velocity profiles from the linearized theory with those obtained by the numerical integration of both interacting and nonsimilar boundary-layer equations. The agreement between the linearized triple-deck and interacting boundary-layer equations is good; however, the nonsimilar boundary layers, which fail to account for upstream influence, are shown to be in poor agreement with both interacting boundary layers and the linearized triple deck. The linearized closed-form solutions will therefore be very useful in future stability calculations.

Nayfeh, A. H.↗

Long-time predictions in nonlinear dynamics

It is known that nonintegrable dynamical systems do not allow precise predictions concerning their behavior for arbitrary long times. The available series solutions are not uniformly convergent according to Poincare's theorem and numerical integrations lose their meaningfulness after the elapse of arbitrary long times. Two approaches are the use of existing global integrals and statistical methods. This paper presents a generalized method along the first approach. As examples long-time predictions in the classical gravitational satellite and planetary problems are treated.

Szebehely, V.↗