SEARCH · Engineering Papers
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.
Time transformations of the equation dt = cr to the n ds, where s is a variable called the intermediate anomaly, are known to reduce global error in the solution of gravitational systems obtained by numerical integration. Attention is given to the Sundman time transformation, and its relation to equations of Keplerian elliptical motion.
The perturbed motion of a rigid body about its center of mass, is formulated in terms of the six elements: l, the magnitude of the angular momentum vector; h, the total energy; delta and epsilon, two linear functions of the independent variable; and psi(1) and theta (1), two Euler angles that orientate the inertial frame with respect to the unperturbed solution. Solutions from the element formulation and the original Euler equations are numerically compared using shuttle-type data. For applied torques smaller than a given magnitude, the element formulation produced the following results: (1) larger step sizes in the numerical integration of the differential equations, resulting in an overall computational time-saving, and (2) more significant figures of accuracy in the computation of the variables describing the state of the rigid body.
Numerical integration of the gravitational N-body problem has been carried out for a variety of protoplanetary clusters in the range N = 100 to 200. Particles are assumed to coagulate at collisions irrespective of relative velocity and mass ratio of the particles. It is shown graphically how the dispersed N-bodies accumulate to a single planet through mutual collisions. The velocity distribution and size distribution of bodies are also investigated as functions of time in the accretion process. Accretion rates of planets are found to be dependent strongly on the initial number density distribution, the initial size distribution, and the initial velocity distribution of bodies. Formation of satellites of about 10% in the planet mass is common to most cases in the present study. A substantial mass of bodies also escapes from the cluster. Many satellites and escapers formed during the accretion process of planets may be source materials of heavy bombardment in the early history of planets.
A model is proposed for single close encounters between two small masses, m1 and m2, which orbit a much larger mass, M. The main new feature of the model is the assumption of conic motion of the center of mass of m1 and m2 in the gravitational field of M. Comparisons of the model with the three-body equations of motion indicate that the model is a useful approximation for m1, m2 not exceeding 0.00001 M. The model is therefore applicable for encounters between bodies of the order of an earth mass or smaller in the presence of the sun. Comparisons are also made of outcomes obtained by the model with outcomes of numerical integration for a large variety of close encounters. The above comparisons reveal that for many purposes the model is an adequate approximation for encounters with an eccentricity of the hyperbolic orbit of m1 about m2 not less than 4.
A third-order solution is developed for the motions of artificial satellites moving in the gravitational field of the earth, whose potential includes the second-, third-, and fourth-order zonal harmonics. Third-order periodic perturbations with fourth-order secular perturbations are derived by Hori's perturbations method. All quantities are expanded into power series of the eccentricity, but the solution is obtained so as to be closed with respect to the inclination. A comparison with the results of numerical integration of the equations of motion indicates that the solution can predict the position of a close-earth satellite with a small eccentricity with an accuracy of better than 1 cm over 1 month.
The precise numerical integration of Cowell's equations of satellite motion is frequently performed with an independent variable s defined by an equation of the form dt = cr to the n-th power ds, where t represents time, r the radial distance from the center of attraction, c is a constant, and n is a parameter. This has been primarily motivated by the 'uniformizing' effects of such a transformation resulting in desirable 'analytic' stepsize control for elliptical orbits. This report discusses the 'proper' choice of the parameter n defining the independent variable s for various types of orbits and perturbation models, and develops a criterion for its selection.
The process by which odd nitrogen species (atomic nitrogen and nitric oxide) are formed during atmospheric entry of meteoroids is analyzed theoretically. An ablating meteoroid is assumed to be a point source of mass with a continuum regime evolving in its wake. The amounts of odd nitrogen species, produced by high-temperature reactions of air in the continuum wake, are calculated by numerical integration of chemical rate equations. Flow properties are assumed to be uniform across the wake, and 29 reactions involving five neutral species and five singly ionized species are considered, as well as vibrational and electron temperature nonequilibrium phenomena. The results, when they are summed over the observed mass, velocity, and entry-angle distribution of meteoroids, provide odd-nitrogen-species annual global production rates as functions of altitude. The peak production of nitric oxide is found to occur at an altitude of about 85 km; atomic nitrogen production peaks at about 95 km. The total annual rate for nitric oxide is 40 million kg; for atomic nitrogen it is 170 million kg.
A rapid orbit generator based on a first-order application of the Generalized Method of Averaging has been developed for the Research and Development (R&D) version of the Goddard Trajectory Determination System (GTDS). The evaluation of the averaged equations of motion can use both numerically averaged and recursively evaluated, analytically averaged perturbation models. These equations are numerically integrated to obtain the secular and long-period motion. Factors affecting efficient orbit prediction are discussed and guidelines are presented for treatment of each major perturbation. Guidelines for obtaining initial mean elements compatible with the theory are presented. An overview of the orbit generator is presented and comparisons with high precision methods are given.
The extensive computation time required for the numerical integration of the differential equations obtained in studies of the concentrations of pollutants emitted by gas turbine combustors, can be reduced significantly by assuming the heat releasing hydrocarbon reactions to be in local equilibrium. In determining the CO and NO concentrations produced in spray combustion, it is, therefore, tempting to use the same local equilibrium assumption in order to simplify the kinetic calculations. An investigation of the validity of the local equilibrium assumption is conducted by using a simple analytical model, and then by actually carrying out the kinetic and local equilibrium calculations for typical case.
A system of differential equations is integrated numerically in space and time over several different matrices in an effort to simulate the atmospheric wave structures which organize severe local storms. Preliminary results with case studies indicate that dynamical fields produced by the numerical simulations can be translated into very fine scale space and time zones where severe storm forecast indices can be developed. These fine scale indices are now available in real-time when run on Langley's STAR 100 computer system.
The maintenance of the quasi-stationary waves obtained through numerically integrating a 2-level quasi-geostrophic spectral model on a beta-plane is investigated. An idealized topography which has only wavenumber n in the zonal direction and the first mode in the meridional direction is used to force the quasi-stationary waves. It is shown that the topographical forcing is not necessarily the mechanism for maintaining the quasi-stationary waves.
Recent work on analytical satellite perturbation theory has involved the completion of a revision to 4th order for zonal harmonics, the addition of a treatment for ocean tides, an extension of the treatment for the noninertial reference system, and the completion of a theory for direct solar-radiation pressure and earth-albedo pressure. Combined with a theory for tesseral-harmonics, lunisolar, and body-tide perturbations, these formulations provide a comprehensive orbit-computation program. Detailed comparisons with numerical integration and observations are presented to assess the accuracy of each theoretical development.
Indirect mapping ocean tides by means of land and island-based tidal gravity measurements is presented. The inverse scheme of linear programming is used for indirect mapping of ocean tides. Open ocean tides were measured by the numerical integration of Laplace's tidal equations.
Aerodynamic experiments were performed on an oscillating NACA 0012 airfoil utilizing a tunnel-spanning wing in both unswept and 30 degree swept configurations. The airfoil was tested in steady state and in oscillatory pitch about the quarter chord. The unsteady aerodynamic loading was measured using pressure transducers along the chord. Numerical integrations of the unsteady pressure transducer responses were used to compute the normal force, chord force, and moment components of the induced loading. The effects of sweep on the induced aerodynamic load response was examined. For the range of parameters tested, it was found that sweeping the airfoil tends to delay the onset of dynamic stall. Sweeping was also found to reduce the magnitude of the unsteady load variation about the mean response. It was determined that at mean incidence angles greater than 9 degrees, sweep tends to reduce the stability margin of the NACA 0012 airfoil; however, for all cases tested, the airfoil was found to be stable in pure pitch. Turbulent eddies were found to convect downstream above the upper surface and generate forward-moving acoustic waves at the trailing edge which move upstream along the lower surface.
A computer model to predict the frost formation process on a flat plate was developed for application to most environmental conditions under which frost occurs. The model was analytically based on a generalized frost thermal conductivity expression, on frost density and thickness rate equations, and on modified heat and mass transfer coefficients designed to fit the available experimental data. The broad experimental ranges reflected by the extremes in ambient humidities, wall temperatures, and convective flow properties in the various publications which were examined served to severely test the flexibility of the model. An efficient numerical integration scheme was developed to solve for the frost surface temperature, density, and thickness under the changing environmental conditions. The comparison of results with experimental data was very encouraging.
The three-dimensional time dependent equations of motion were numerically integrated for fully-developed turbulent channel flow. A large scale flow field was obtained directly from the solution of these equations, and small scale field motions were simulated through an eddy viscosity model. The calculations were carried out on the ILLIAC 4 computer. The computed flow patterns show that the wall layer consists of coherent structures of low speed and high speed streaks alternating in the spanwise direction. These structures were absent in the regions away from the wall. Hot spots, small localized regions of very large turbulent shear stress, were frequently observed. The profiles of the pressure velocity-gradient correlations show a significant transfer of energy from the normal to the spanwise component of turbulent kinetic energy in the immediate neighborhood of the wall ('the splatting effect').
An analysis of the dumbbell gravity gradiometer concept for measuring short wavelength variations in the earth's gravity gradient is presented. Variations in the gradient are recorded by measuring tension variations in a vertically stabilized satellite consisting of heavy masses connected by a long wire or rod. Tension noise arises from the excitation of various mechanical oscillations of the system. The principal noise sources that were identified are fluctuations in atmospheric drag heating and drag force resulting from density variations and winds. Approximate analytical expressions are presented for the tension noise as a function of the system design parameters for various possible configurations. Computer simulations using numerical integration were performed to study the tension noise for several sample cases. Three designs consistent with Shuttle launch capabilities are discussed.
A revised orbit and extended ephemeris are calculated for Leda, the thirteenth satellite of Jupiter, by means of numerical integration. Thirty-one observations between September 1974 and December 1977 satisfy the new orbit to 0.9 arcsec on the average.