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 685 records · Page 38

Second-order small-disturbance solutions for hypersonic flow over power-law bodies

Similarity solutions were found which give the adiabatic flow of an ideal gas about two-dimensional and axisymmetric power-law bodies at infinite Mach number to second order in the body slenderness parameter. The flow variables were expressed as a sum of zero-order and perturbation similarity functions for which the axial variations in the flow equations separated out. The resulting similarity equations were integrated numerically. The solutions, which are universal functions, are presented in graphic and tabular form. To avoid a singularity in the calculations, the results are limited to body power-law exponents greater than about 0.85 for the two-dimensional case and 0.75 for the axisymmetric case. Because of the entropy layer induced by the nose bluntness (for power-law bodies other than cones and wedges), only the pressure function is valid at the body surface. The similarity results give excellent agreement with the exact solutions for inviscid flow over wedges and cones having half-angles up to about 20 deg. They give good agreement with experimental shock-wave shapes and surface-pressure distributions for 3/4-power axisymmetric bodies, considering that Mach number and boundary-layer displacement effects are not included in the theory.

Townsend, J. C.↗

STICAP: A linear circuit analysis program with stiff systems capability. Volume 1: Theory manual

STICAP (Stiff Circuit Analysis Program) is a FORTRAN 4 computer program written for the CDC-6400-6600 computer series and SCOPE 3.0 operating system. It provides the circuit analyst a tool for automatically computing the transient responses and frequency responses of large linear time invariant networks, both stiff and nonstiff (algorithms and numerical integration techniques are described). The circuit description and user's program input language is engineer-oriented, making simple the task of using the program. Engineering theories underlying STICAP are examined. A user's manual is included which explains user interaction with the program and gives results of typical circuit design applications. Also, the program structure from a systems programmer's viewpoint is depicted and flow charts and other software documentation are given.

Cooke, C. H.↗

A two parameter survey of periodic orbits in the restricted problem of three bodies

Within the context of the restricted problem of three bodies the effects caused by varying the mass ratio of the primaries and the eccentricity of their orbits, upon periodic orbits of the infinitesimal mass which are numerical continuations of circular orbits in the ordinary problem of two bodies are shown. A recursive power series technique is used to numerically integrate the equations of motion as well as the first variational equations in order to generate a two parameter family of perodic orbits and identify the linear stability characteristics. Seven such families are investigated with equally spaced mass ratios from 0.0 to 1.0 and eccentricities of the orbits of the primaries in a range 0.0 to 0.6. Stable orbits are associated with large distances of the infinitesimal mass from the perturbing primary, nearly circular motion of the primaries, and small mass ratios of the primaries. Unstable orbits for the infinitesimal mass are associated with small distances from the perturbing primary, highly elliptic orbits of the primaries and large mass ratios.

Shelus, P. J.↗

Lunar motion analysis and laser data management

Work completed in lunar motion analysis and laser data management during the period July 1, 1971 - September 30, 1975 was reported. In this context, analysis refers to theoretical or numerical studies involving real or potential applications of such observations to improvement of the physical model, and data management refers to the process by which observed photon events are turned into observations and are made available to potential users. The data analysis work included: (1) bringing to operational status of computer programs for the numerical integration of the lunar orbit motion and for the application of lunar laser time delays for the improvement of the parameters of the physical model, (2) program improvement and program integrity, (3) three-dimensional ephemeris, and (4) miscellaneous independent studies. The data management work included: (1) data identification, (2) observatory interfaces, and (3) data distribution.

Mulholland, J. D.↗

Earth-moon system: Dynamics and parameter estimation

A theoretical development of the equations of motion governing the earth-moon system is presented. The earth and moon were treated as finite rigid bodies and a mutual potential was utilized. The sun and remaining planets were treated as particles. Relativistic, non-rigid, and dissipative effects were not included. The translational and rotational motion of the earth and moon were derived in a fully coupled set of equations. Euler parameters were used to model the rotational motions. The mathematical model is intended for use with data analysis software to estimate physical parameters of the earth-moon system using primarily LURE type data. Two program listings are included. Program ANEAMO computes the translational/rotational motion of the earth and moon from analytical solutions. Program RIGEM numerically integrates the fully coupled motions as described above.

Breedlove, W. J., Jr.↗

Combustion of hydrogen injected into a supersonic airstream (a guide to the HISS computer program)

A computer program based on a finite-difference, implicit numerical integration scheme is described for the prediction of hydrogen injected into a supersonic airstream at an angle ranging from normal to parallel to the airstream main flow direction. Results of calculations for flow and thermal property distributions were compared with 'cold flow data' taken by NASA/Langley and show excellent correlation. Typical results for equilibrium combustion are presented and exhibit qualitatively plausible behavior. Computer time required for a given case is approximately one minute on a CDC 7600. A discussion of the assumption of parabolic flow in the injection region is given which demonstrates that improvement in calculation in this region could be obtained by a partially-parabolic procedure which has been developed. It is concluded that the technique described provides an efficient and reliable means for analyzing hydrogen injection into supersonic airstreams and the subsequent combustion.

Dyer, D. F.↗

Wide angle lenses with aspheric correcting surfaces

Refracting elements having aspheric correcting surfaces near the center of curvature are analyzed. As in the case of the Schmidt reflector, such systems can have a wide aperture and a wide field of view; in addition, they are free from obstructing surfaces. Being uncorrected for dispersion, however, the refracting systems are restricted to nearly monochromatic radiation. Typical forms of the correcting surfaces have been determined, both by third order theory and by numerical integration of exact equations.

Jones, R. T.↗

On the Tesseral-Harmonics Resonance Problem in Artificial-Satellite Theory, Part 2

Equations were derived for the perturbations on an artificial satellite when the motion of the satellite is commensurable with that of the earth. This was done by first selecting the tesseral harmonics that contribute the most to the perturbations and then by applying Hori's method by use of Lie series. Here, are introduced some modifications to the perturbations, which now result in better agreement with numerical integration.

Romanowicz, B. A.↗

Satellite-tracking and earth-dynamics research programs

Observations and research progress of the Smithsonian Astrophysical Observatory are reported. Satellite tracking networks (ground stations) are discussed and equipment (Baker-Nunn cameras) used to observe the satellites is described. The improvement of the accuracy of a laser ranging system of the ground stations is discussed. Also, research efforts in satellite geodesy (tides, gravity anomalies, plate tectonics) is discussed. The use of data processing for geophysical data is examined, and a data base for the Earth and Ocean Physics Applications Program is proposed. Analytical models of the earth's motion (computerized simulation) are described and the computation (numerical integration and algorithms) of satellite orbits affected by the earth's albedo, using computer techniques, is also considered. Research efforts in the study of the atmosphere are examined (the effect of drag on satellite motion), and models of the atmosphere based on satellite data are described.

Source record↗

Sonic boom research

A computer program for CDC 6600 is developed for the nonlinear sonic boom analysis including the asymmetric effect of lift near the vertical plane of symmetry. The program is written in FORTRAN 4 language. This program carries out the numerical integration of the nonlinear governing equations from the input data at a finite distance from the airplane configuration at a flight altitude to yield the pressure signitude at ground. The required input data and the format for the output are described. A complete program listing and a sample calculation are given.

Zakkay, V.↗

LaRC design analysis report for National Transonic Facility for 9% nickel tunnel shell. Volume 1: Finite difference analysis of cone/cylinder junction

The results of extensive computer (finite element, finite difference and numerical integration), thermal, fatigue, and special analyses of critical portions of a large pressurized, cryogenic wind tunnel (National Transonic Facility) are presented. The computer models, loading and boundary conditions are described. Graphic capability was used to display model geometry, section properties, and stress results. A stress criteria is presented for evaluation of the results of the analyses. Thermal analyses were performed for major critical and typical areas. Fatigue analyses of the entire tunnel circuit are presented.

Ramsey, J. W., Jr.↗

LaRC design analysis report for National Transonic Facility for 304 stainless steel tunnel shell. Volume 1S: Finite difference analysis of cone/cylinder junction

The results of extensive computer (finite element, finite difference and numerical integration), thermal, fatigue, and special analyses of critical portions of a large pressurized, cryogenic wind tunnel (National Transonic Facility) are presented. The computer models, loading and boundary conditions are described. Graphic capability was used to display model geometry, section properties, and stress results. A stress criteria is presented for evaluation of the results of the analyses. Thermal analyses were performed for major critical and typical areas. Fatigue analyses of the entire tunnel circuit are presented.

Ramsey, J. W., Jr.↗

Lunar and solar perturbations on the orbit of a geosynchronous satellite

The luni-solar effects cause a large amplitude, long-period perturbation of the orbital plane. Canonical differential equations associated with this motion contain a singular point, and are expanded about this point to third order. A solution is given that is valid for long times and is not restricted to small inclinations. Higher order terms are investigated. Expressions are given for the inclination (I) and node as a function of time. Comparisons with a numerically integrated solution show a disagreement in inclination of only .03 deg after 11 years.

Graf, O. F., Jr.↗

On the alleged collisional origin of the Kirkwood Gaps

This paper examines two proposed mechanisms whereby asteroidal collisions and close approaches may have given rise to the Kirkwood Gaps. The first hypothesis is that asteroids in near-resonant orbits have markedly increased collision probabilities and so are preferentially destroyed, or suffer decay in population density, within the resonance zones. A simple order-of-magnitude analysis shows that this hypothesis is untenable since it leads to conclusions which are either unrealistic or not in accord with present understanding of asteroidal physics. The second hypothesis is the Brouwer-Jefferys theory that collisions would smooth an asteroidal distribution function, as a function of Jacobi constant, thus forming resonance gaps. This hypothesis is examined by direct numerical integration of 50 asteroid orbits near the 2:1 resonance, with collisions simulated by random variables. No tendency to form a gap was observed.

Heppenheimer, T. A.↗

Open-loop response of a burning liquid monopropellant

A direct comparison was obtained between the theoretical and experimental response of a burning liquid monopropellant surface to imposed pressure oscillations. The comparison extended beyond the quasi-steady frequency regime into the range where the combustion process interacts with the thermal wave in the liquid phase. The theory considered both gas and liquid-phase transient effects, involving a numerically integrated, first-order perturbation analysis. The experiment employed hydrazine as the test fuel, burning as a strand. With kinetic parameters fixed by steady burning rate measurements, the model gave good predictions of liquid-phase temperature distributions during steady combustion as well as the amplitude and phase of the oscillatory component of the burning rate in the presence of imposed pressure oscillations. The combustion response tends to peak in two frequency ranges; one corresponding to interaction with liquid-phase transient effects, which could be verified both theoretically and experimentally; the second corresponding to gas-phase transient effects, which could only be examined theoretically. Both ranges provide conditions where the amplifying power is sufficient to drive combustion instability.

Allison, C. B.↗

Dynamics of spin stabilized spacecraft during deployment of telescoping appendages

The dynamics of a spin stabilized spacecraft with telescoping appendages are treated both analytically and numerically. Two types of telescoping appendages are considered: (a) where an end mass is mounted at the end of an (assumed) massless boom; and (b) where the appendage is assumed to consist of a uniformly distributed homogeneous mass throughout its length. For the telescoping system Eulerian equations of motion are developed. During all deployment sequences it is assumed that the transverse component of angular momentum is much smaller than the component along the major spin axis. Closed form analytical solutions for the time response of the transverse components of angular velocities are obtained when the spacecraft hub has a nearly spherical mass distribution. For the more general case, a series solution is obtained and this solution is limited by its radius of convergence. The comparison of the different approximate analytical methods with numerical integration results are studied and it is observed that the oscillatory nature of the responses of the transverse angular velocity components reduces rapidly with faster extension rates.

Sellappan, R.↗

Thermospheric winds and exospheric temperatures from incoherent scatter radar measurements in four seasons

The time-dependent equations of neutral air motion are solved subject to three constraints: two of them are the usual upper and lower boundary conditions and the third is the value of the wind-induced ion drift at any given height. Using incoherent radar data, this procedure leads to a fast, direct numerical integration of the two coupled differential equations describing the horizontal wind components and yields time dependent wind profiles and meridional exospheric neutral temperature gradients. The diurnal behavior of the neutral wind system and of the exospheric temperature is presented for two solstice and two equinox days. The data used were obtained by the St. Santin and the Millstone Hill incoherent scatter radars. The derived geographic distributions of the exospheric temperatures are compared with those predicted by the OGO-6 empirical thermospheric model.

Antoniadis, D. A.↗