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 739 records · Page 41

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Chaos and nonlinear dynamics of single-particle orbits in a magnetotaillike magnetic field

The properties of charged-particle motion in Hamiltonian dynamics are studied in a magnetotaillike magnetic field configuration. It is shown by numerical integration of the equation of motion that the system is generally nonintegrable and that the particle motion can be classified into three distinct types of orbits: bounded integrable orbits, unbounded stochastic orbits, and unbounded transient orbits. It is also shown that different regions of the phase space exhibit qualitatively different responses to external influences. The concept of 'differential memory' in single-particle distributions is proposed. Physical implications for the dynamical properties of the magnetotail plasmas and the possible generation of non-Maxwellian features in the distribution functions are discussed.

Chen, J.↗

NASA NTRS1986

Shell correction for the stopping power of K electrons

In view of the inapplicability of the asymptotic expressions for the stopping number available in the literature at high energies, an alternative approach is taken to compute the shell correction to the stopping number of K electrons. Anholt's formula for the K-shell ionization has been used to calculate the excitation function for longitudinal interaction and numerical integration over energy has been carried out to evaluate the shell correction. Comparison with other theoretical calculations is made. It is proposed that, with the inclusion of relativistic effects, an asymptotic expansion of the stopping number with a leading-term logarithmic in the energy of the incident particle would be more meaningful and might enable one to extract the relativistic contribution to the shell correction from it.

Leung, P. T.↗

NASA NTRS1985

Reassessment of the theory of stimulated Raman scattering

A modification of the standard theory of stimulated Raman scattering (SRS) first proposed by Sparks (1974, 1975) is analyzed and shown to incorporate a possibly important physical effect; however, its original formulation is incorrect. The analysis is based on an exact numerical integration of the coupled equations of the modified theory, the results of which are compared with both the conventional theory of SRS and with one set of experimental data. A reformulation of the modified theory is suggested that leads to a gain which is in somewhat better agreement with the data than is the conventional theory.

Fralick, G. C.↗

NASA NTRS1986

Three-dimensional computation of mixing of transverse injector in a ducted supersonic airstream

Numerical solutions of the three-dimensional mass-averaged elliptic Navier-Stokes equations, including species transport, are obtained for nonreacting, turbulent, mixing flow fields for the case of transverse sonic injection of a secondary gas into a supersonic airstream through a circular orifice injector. Results are presented for flow through a constant area duct and through a duct with a rearward-facing step upstream of the injector. The equations are numerically integrated using MacCormack's explicit method and turbulence is included using the Baldwin-Lomax algebraic eddy viscosity model. In the species transport and energy equations, diffusion coefficients based on Fick's law and an assumption of unit Lewis number are applied. The computations were performed on a CDC-VPS-32 (extended version of Cyber-205) using a grid consisting of approximately 200,000 points. The computed results are compared with experimentally observed penetration and spreading boundaries for an injected gas at two dynamic pressure ratios. Three-dimensional flow field structures are dipicted in terms of static pressure, mass fractions of species and velocity vectors.

Uenishi, K.↗

NASA NTRS1986

Semi-analytical theory for long-term rotational motion of uncontrolled satellites

Analytical solutions/integrals of the free rotational motion of uncontrolled satellites which contain spinning rotors are employed in conjunction with a suitable analytical/numerical averaging approach to predict long-term attitude motion. This semi-analytical method of averaging is basically similar to the one used in purely analytical 'general perturbation' theories. However, in view of the complexities involved in analytical treatment of the problem, its application is primarily numerical in character. An example is presented illustrating the use of the theory to predict the attitude motion of an axisymmetric satellite containing symmetrical constant-speed rotors under gravity-gradient torque. The use of the theory is considerably more efficient than straightforward numerical integration of the complete set of equations of motion.

Lahr, B. S.↗

NASA NTRS1986

High-accuracy computation of Delta V magnitude probability densities - Preliminary remarks

This paper describes an algorithm for the high accuracy computation of some statistical quantities of the magnitude of a random trajectory correction maneuver (TCM). The trajectory correction velocity increment Delta V is assumed to be a three-component random vector with each component being a normally distributed random scalar having a possibly nonzero mean. Knowledge of the statitiscal properties of the magnitude of a random TCM is important in the planning and execution of maneuver strategies for deep-space missions such as Galileo. The current algorithm involves the numerical integration of a set of differential equations. This approach allows the computation of density functions for specific Delta V magnitude distributions to high accuracy without first having to generate large numbers of random samples. Possible applications of the algorithm to maneuver planning, planetary quarantine evaluation, and guidance success probability calculations are described.

Chadwick, C.↗

NASA NTRS1987

Calculating Sonic-Boom Propagation

Nonlinear effects included, enabling more-realistic modeling. Modified Method of Characteristics Sonic Boom Extrapolation Program (MMOC) is computer program for sonic-boom propagation that includes shock coalescence and incorporates effects of asymmetries due to volume and lift. Numerically integrates nonlinear governing equations using data on initial data line approximately one body length from aircraft and yields sonic-boom pressure at ground as function of time or of position at given time. MMOC written in FORTRAN IV.

Darden, C. M.↗

NASA NTRS1987

Parallel Algorithm Solves Coupled Differential Equations

Numerical methods adapted to concurrent processing. Algorithm solves set of coupled partial differential equations by numerical integration. Adapted to run on hypercube computer, algorithm separates problem into smaller problems solved concurrently. Increase in computing speed with concurrent processing over that achievable with conventional sequential processing appreciable, especially for large problems.

Hayashi, A.↗

NASA NTRS1986

An all-purpose computer program for lidar simulations and reduction of lidar data

A general computer program was written for the dual purpose of simulating lidar returm signals and extracting H2O vapor and temperature (T) profiles from actual lidar returns. This program will be used to retrieve H2O vapor and T profiles from data obtained with the University of Maryland ground-based lidar systems. The program is designed to simulate and to analyze lidar returns from a monostatic system for two cases: when the output laser pulse is monochromatic, and when the pulse has a known wavelength spectrum. The simulation part of the program has several model H2O vapor, T, and aerosol profiles which can be combined to form different model atmospheres. The simulation program also permits the testing of different methods of numerical filtering to suppress noise, and testing the accuracy of different methods of numerical integration.

Braun, C.↗

NASA NTRS1986

The nonlinear development of Gortler vortices in growing boundary layers

The development of Gortler vortices in boundary layers over curved walls in the nonlinear regime is investigated. The growth of the boundary layer makes a parallel flow analysis impossible except in the high wavenumber regime so in general the instability equations must be integrated numerically. Here the spanwise dependence of the basic flow is described using Fourier series expansion whilst the normal and streamwise variations are taken into account using finite differences. The calculations suggest that a given disturbance imposed at some position along the wall will eventually reach a local equilibrium state essentially independent of the initial conditions. In fact, the equilibrium state reached is qualitatively similar to the large amplitude high wave-number solution described asymptotically by Hall (1982). In general, it is found that the nonlinear interactions are dominated by a mean field type of interaction between the mean flow and the fundamental. Thus, even though higher harmonics of the fundamental are necessarily generated, most of the disturbance energy is confined to the mean flow correction and the fundamental. A major result of the calculations is finding that the downstream velocity field develops a strongly inflection character as the flow moves downstream. The latter result suggests that the major effect of Gortler vortices on boundary layers of practical importance might be to make them highly receptive to rapidly growing Rayleigh modes of instability.

Hall, Philip↗

NASA NTRS1987

Transient response for interaction of two dynamic bodies

During the launch sequence of any space vehicle complicated boundary interactions occur between the vehicle and the launch stand. At the start of the sequence large forces exist between the two; contact is then broken in a short but finite time which depends on the release mechanism. The resulting vehicle response produces loads which are very high and often form the design case. It is known that the treatment of the launch pad as a second dynamic body is significant for an accurate prediction of launch response. A technique was developed for obtaining loads generated by the launch transient with the effect of pad dynamics included. The method solves uncoupled vehicle and pad equations of motion. The use of uncoupled models allows the simulation of vehicle launch in a single computer run. Modal formulation allows a closed-form solution to be written, eliminating any need for a numerical integration algorithm. When the vehicle is on the pad the uncoupled pad and vehicle equations have to be modified to account for the constraints they impose on each other. This necessitates the use of an iterative procedure to converge to a solution, using Lagrange multipliers to apply the required constraints. As the vehicle lifts off the pad the coupling between the vehicle and the pad is eliminated point by point until the vehicle flies free. Results obtained by this method were shown to be in good agreement with observed loads and other analysis methods. The resulting computer program is general, and was used without modification to solve a variety of contact problems.

Prabhakar, A.↗

NASA NTRS1987

A scaling law for accretion zone sizes

Current theories of runaway planetary accretion require small random velocities of the accreted particles. Two body gravitational accretion cross sections which ignore tidal perturbations of the Sun are not valid for the slow encounters which occur at low relative velocities. Wetherill and Cox have studied accretion cross sections for rocky protoplanets orbiting at 1 AU. Using analytic methods based on Hill's lunar theory, one can scale these results for protoplanets that occupy the same fraction of their Hill sphere as does a rocky body at 1 AU. Generalization to bodies of different sizes is achieved here by numerical integrations of the three-body problem. Starting at initial positions far from the accreting body, test particles are allowed to encounter the body once, and the cross section is computed. A power law is found relating the cross section to the radius of the accreting body (of fixed mass).

Greenzweig, Yuval↗

NASA NTRS1987

Research in computational fluid dynamics

The numerical integration of quasi-one-dimensional unsteady flow problems which involve finite rate chemistry are discussed, and are expressed in terms of conservative form Euler and species conservation equations. Hypersonic viscous calculations for delta wing geometries is also examined. The conical Navier-Stokes equations model was selected in order to investigate the effects of viscous-inviscid interations. The more complete three-dimensional model is beyond the available computing resources. The flux vector splitting method with van Leer's MUSCL differencing is being used. Preliminary results were computed for several conditions.

Murman, Earll M.↗

NASA NTRS1986

The global evolution of the primordial solar nebula

Complete radial, time-dependent calculations of the structure and evolution of the primordial solar nebula during the viscous diffusion stage are presented. The viscous stress is derived from analytic one-zone models of the vertical nebular structure based on detailed grain opacities. Comparisons with full numerical integrations indicate that the effective viscous alpha parameter is about 0.01. The evolution time of a minimum mass nebula is one-million yr or less. The flow pattern of fluid elements in the disk is examined and the implications the results have on the theory of the formation of the solar system are discussed.

Ruden, S. P.↗

NASA NTRS1985

Buckling of simply supported plate assemblies subject to shear loading

An analysis is presented of the shear buckling load of an arbitrary simply-supported prismatic plate assembly or stiffened panel configuration. The method utilizes stiffness matrices developed by Wittrick (1968). For the restraint of simple support along edges perpendicular to the longitudinal direction, various wavelengths of response are combined using an energy method that involves only a boundary integral. Numerical results for plate and panel configurations show the accuracy and convergence of the method.

Anderson, M. S.↗

NASA NTRS1986

A general theory for the Uranian satellites

A general analytical theory of the five main satellites of Uranus, including the secular and short period terms hereafter denoted by GUST, is presented. A comparison is made with an internal numerical integration with nominal masses of Veillet (1983). The precision of the theory goes from about 10 km for Miranda to 100 km for Oberon. The short period terms in the motions of Titania and Oberon are larger than 500 km. They should make possible the determination of the masses of the outer satellites through the optical data of Voyager encounter.

Laskar, J.↗

NASA NTRS1986

Formation and evolution of a circumterrestrial disk Constraints on the origin of the moon in geocentric orbit

A data base of about 25,000 numerically integrated trajectories of earth-encountering planetesimals is used to study the angular momentum problem of forming the moon out of material captured into a circumterrestrial disk from heliocentric orbits. Mass-orbital element distributions of incoming planetesimals are combined with this data base to calculate, as a function of distance from earth, the net geocentric specific angular momentum of disk-encountering material on heliocentric orbits. Results suggest that a permanent disk population is not possible.

Herbert, Floyd↗

NASA NTRS1987

Dynamic analysis of satellites with deployable hinged appendages

The nonlinear equations of motion determining the planar dynamical behavior of an orbiting satellite deploying both one and two rigid appendages have been formulated using Lagrange's equations. The analysis accounts for large angle rotations, Coriolis effects, and the gravitational gradient, and the resulting coupled governing equations are integrated numerically. The analysis is applied to the Space Shuttle based deployment of rigid truss-like members, and results show that spacecraft inertia parameters, appendage mass and length, deployment velocity, and initial conditions all influence the system response. It is found that the resulting librational movement is related to the size of the deployment payload, and that gravitational forces lead to vehicle stabilization.

Oakes, Kevin F.↗