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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Three-dimensional computations of transverse hydrogen jet combustion in a supersonic airstream

A computational fluid dynamics (CFD) code is being developed to compute the mixing and combustion of hydrogen fuel in the turbulent flow fields of supersonic combustion ramjets (scramjet). The code solves the three-dimensional Reynolds time-averaged complete Navier-Stokes equations including transport equations for a four species, two reaction, global finite rate chemistry model. The code was applied to the case of transverse injection of hydrogen from a sonic circular orifice into a supersonic airstream. The equations were numerically integrated using MacCormack's explicit method, and the algebraic eddy viscosity model of Baldwin-Lomax was used to model the turbulence. In the species transport and energy equations, diffusion coefficients based on Fick's Law and an assumption of unit Lewis number were applied. Computed features of the three-dimensional flow field are depicted by static pressure, static temperature, mass fraction of species, and velocity vectors. For engineering interest, mixing and combustion parameters were examined to assess the effect of injector diameter, injected fuel pressure, fuel-air ratio, and spacing of fuel injectors. The objective of the present paper is to demonstrate the capability of the present three-dimensional spatially elliptic, CFD code for turbulent, reacting flow. Application of the code to specific supersonic combustion configurations is planned.

Uenishi, K.↗
NASA NTRS1986

Numerical planetary and lunar ephemerides - Present status, precision and accuracies

Features of the emphemeris creation process are described with attention given to the equations of motion, the numerical integration, and the least-squares fitting process. Observational data are presented and ephemeride accuracies are estimated. It is believed that radio measurements, VLBI, occultations, and the Space Telescope and Hipparcos will improve ephemerides in the near future. Limitations to accuracy are considered as well as relativity features. The export procedure, by which an outside user may obtain and use the JPL ephemerides, is discussed.

Standish, E. Myles, Jr.↗
NASA NTRS1987

Stability analysis of two coupled Lorenz lasers and the coupling-induced periodic to chaotic transition

A linear stability analysis of two Lorenz lasers coupled by their electric fields has been performed, and it is shown that the bad cavity condition becomes a function of coupling and that a good cavity instability may occur if the injected fields are inverted before injection. In addition, it is shown that the symmetrically coupled Lorenz system is isomorphic to the original Lorenz system with new parameters. The stability analysis also predicts a lowering of the second laser threshold with coupling for both the chaotic and self-pulsing regimes. Numerical integration of the equations is in agreement with these predictions and has revealed a coupling induced transition from self-pulsing to chaotic behavior. The classification of the behavior of the coupled system in the parameter space of the coupling constants has been investigated and shows that the results of symmetric coupling allow enough of a margin for an experimental test of the theory. This would allow experimentalists to observe the actual Lorenz instability at excitations as low as 4-5 times above threshold.

Lawandy, N. M.↗
NASA NTRS1987

An approximate semi-analytical method for prediction of interlaminar shear stresses in an arbitrarily laminated thick plate

An approximate semianalytical method for determination of interlaminar shear stress distribution through the thickness of an arbitrarily laminated thick plate has been presented. The method is based on the assumptions of transverse inextensibility and layerwise constant shear angle theory (LCST) and utilizes an assumed quadratic displacement potential energy based finite element method (FEM). Centroid of the triangular surface has been proved from a rigorous mathematical point of view (Aubin-Nitsche theory), to be the point of exceptional accuracy for the interlaminar shear stresses. Numerical results indicate close agreement with the available three-dimensional elasticity theory solutions. A comparison between the present theory and that due to an assumed stress hybrid FEM suggest that the (normal) traction-free-edge condition is not satisfied in the latter approach. Furthermore, the present paper is the first to present the results for interlaminar shear stresses in a two-layer thick square plate of balanced unsymmetric angle-ply construction. A comparison with the recently proposed Equilibrium Method (EM) indicates the superiority of the present method, because the latter assures faster convergence as well as simultaneous vanishing of the transverse shear stresses on both of the exposed surfaces of the laminate. Superiority of the present method over the EM, in the case of a symmetric laminate, is limited to faster convergence alone. It has also been demonstrated that the combination of the present method and the reduced (quadratic order) numerical integration scheme yields convergence of the interlaminar shear stresses almost as rapidly as that of the nodal displacements, in the case of a thin plate.

Chaudhuri, Reaz A.↗
NASA NTRS1987

Investigation of turbulent shear flows by numerical simulation

A data base (Kim et al., 1986) generated by direct numerical integration of the Navier-Stokes equations for unsteady three-dimensional turbulent channel flow at low Reynolds numbers is applied to investigate the physical processes acting in wall-bounded shear flows. The techniques used to construct the data base are outlined; the modeling of turbulence, streaks and streamwise vortices, horseshoe and inverted-horseshoe vortices, and vortex rings is described; the conditional sampling techniques employed are discussed; and typical results are presented graphically. An organized vorticity-field structure is revealed and found to be associated with the bursting process and a local peak in the wall pressure, indicating a possible approach to reducing viscous drag.

Kim, J.↗
NASA NTRS1987

A piloted simulator evaluation of a ground-based 4D descent advisor algorithm

A ground-based, four-dimensional (4D) descent-advisor algorithm is under development at NASA Ames Research Center. The algorithm combines detailed aerodynamic, propulsive, and atmospheric models with an efficient numerical integration scheme to generate 4D descent advisories. This paper investigates the ability of the 4D descent advisor algorithm to provide adequate control of arrival time for aircraft not equipped with on-board 4D guidance systems. A piloted simulation was conducted to determine the precision with which the descent advisor could predict the 4D trajectories of typical straight-in descents flown by airline pilots under different wind conditions. The effects of errors in the estimation of wind and initial aircraft weight were also studied. A description of the descent advisor as well as the results of the simulation studies are presented.

Green, Steven M.↗
NASA NTRS1987

Wind measurements for non-uniform wind fields from spaceborne scatterometers

Radar backscattering coefficient measurements by spaceborne scatterometers are presently simulated for the case of nonuniform wind fields, by means of a detailed numerical integration of the radar equation. The winds thus estimated are then compared with a nominal field which is defined as the average wind vector over the wind cell. The simulation results obtained for the NASA scatterometer are presented for cases of random wind fields whose spectra are consistent with the Seasat scatterometer sea surface wind spectrum. When the nonuniformity is small, system noise dominates the wind error; wind error degradation is therefore small for both perfect and imperfect coregistration cases. When it is relatively large, however, the wind error degradation persistently increases for both perfect and imperfect coregistrations.

Chi, Chong-Yung↗
NASA NTRS1988

Estimating Geophysical Parameters From Gravity Data

ORBSIM program developed for accurate extraction of parameters of geophysical models from Doppler-radio-tracking data acquired from orbiting planetary spacecraft. Model of proposed planetary structure used in numerical integration along simulated trajectories of spacecraft around primary body. Written in FORTRAN 77.

Sjogren, William L.↗
NASA NTRS1987

Flap-lag-torsion aeroelastic stability of a circulation control rotor in forward flight

The aeroelastic stability of a circulation control rotor blade undergoing three degrees of motion (flap, lag, and torsion) is investigated in forward flight. Quasi-steady strip theory is used to evaluate the aerodynamics forces; and the airfoil characteristics are from data tables. The propulsive and the auxiliary power trims are calculated from vehicle and rotor equilibrium equations through the numerical integration of element forces in azimuth as well as in radial directions. The nonlinear time dependent periodic blade response is calculated using an iterative procedure based on Floquet theory. The periodic perturbation equations are solved for stability using Floquet transition matrix theory. The effects of several parameters on blade stability are examined, including advance ratio, collective pitch, thrust level, shaft tilt, structural stiffnesses variation, and propulsive and auxiliary power trims.

Chopra, Inderjit↗
NASA NTRS1986

Trajectory characteristics and heating of hypervelocity projectiles having large ballistic coefficients

A simple, approximate equation describing the velocity-density relationship (or velocity-altitude) has been derived from the flight of large ballistic coefficient projectiles launched at high speeds. The calculations obtained by using the approximate equation compared well with results for numerical integrations of the exact equations of motion. The flightpath equation was used to parametrically calculate maximum body decelerations and stagnation pressures for initial velocities from 2 to 6 km/s. Expressions were derived for the stagnation-point convective heating rates and total heat loads. The stagnation-point heating was parametrically calculated for a nonablating wall and an ablating carbon surface. Although the heating rates were very high, the pulse decayed quickly. The total nose-region heat shield weight was conservatively estimated to be only about 1 percent of the body mass.

Tauber, Michael E.↗
NASA NTRS1987

Computer architecture for efficient algorithmic executions in real-time systems: New technology for avionics systems and advanced space vehicles

Improvements and advances in the development of computer architecture now provide innovative technology for the recasting of traditional sequential solutions into high-performance, low-cost, parallel system to increase system performance. Research conducted in development of specialized computer architecture for the algorithmic execution of an avionics system, guidance and control problem in real time is described. A comprehensive treatment of both the hardware and software structures of a customized computer which performs real-time computation of guidance commands with updated estimates of target motion and time-to-go is presented. An optimal, real-time allocation algorithm was developed which maps the algorithmic tasks onto the processing elements. This allocation is based on the critical path analysis. The final stage is the design and development of the hardware structures suitable for the efficient execution of the allocated task graph. The processing element is designed for rapid execution of the allocated tasks. Fault tolerance is a key feature of the overall architecture. Parallel numerical integration techniques, tasks definitions, and allocation algorithms are discussed. The parallel implementation is analytically verified and the experimental results are presented. The design of the data-driven computer architecture, customized for the execution of the particular algorithm, is discussed.

Carroll, Chester C.↗
NASA NTRS1988

The viscoplasticity theory based on overstress applied to the modeling of a nickel base superalloy at 815 C

Short term strain rate change, creep and relaxation tests were performed in an MTS computer controlled servohydraulic testing machine. Aging and recovery were found to be insignificant for test times not exceeding 30 hrs. The material functions and constants of the theory were identified from results of strain rate change tests. Numerical integration of the theory for relaxation and creep tests showed good predictive capabilities of the viscoplasticity theory based on overstress.

Krempl, E.↗
NASA NTRS1988

An error floor in tone calibrated transmission

Use of a low level pilot tone has been shown to eliminate the error floor in fading channels. This paper demonstrates that non-idealities in the receiver's pilot tone filter cause reappearance of the error floor. It also presents the bit error rate (BER) in closed form, in contrast to the multidimensional numerical integration of previous work.

Cavers, James K.↗
NASA NTRS1988

A study of the dynamics of the Intertropical Convergence Zone (ITCZ) in a symmetric atmosphere-ocean model

A numerical model of the circulation of a coupled axisymmetric atmosphere-ocean system was constructed to investigate the physical factors governing the location and intensity of the Intertropical Convergence Zone (ITCZ) over oceans and over land. The results of several numerical integrations are presented to illustrate the interaction of the individual atmospheric and oceanic circulations. It is shown that the ITCA cannot be located at the equator because the atmosphere-ocean system is unstable for lateral displacements of the ITCA from an equilibrium position at the equator.

Charney, J. G.↗
NASA NTRS1987

High resolution absorption cross-sections and band oscillator strengths of the Schumann-Runge bands of oxygen at 79 K

Cross sections of O2 at 79 K have been obtained from photoabsorption measurements at various pressures throughout the wavelength region 179.3-198.0 nm with a 6.65-m photoelectric scanning spectrometer equipped with a 2400-lines/mm grating and having an instrumental width (FWHM) of 0.0013 nm. The measured absorption cross sections of the Schumann-Runge bands (12,0) through (2,0) are independent of the instrumental width. The measured cross-sections are presented graphically here and are available at wavenumber intervals of about 0.1/cm as numerical compilations stored on magnetic tape from the National Space Science Data Center, NASA/Goddard. Band oscillator strengths of these bands have been determined by direct numerical integration of the measured cross sections.

Yoshino, K.↗
NASA NTRS1988

Temporal and spatial adaptive algorithm for reacting flows

A numerical integration scheme for quasi-one-dimensional unsteady flows with finite-rate chemistry is developed and demonstrated. The governing Euler and species-conservation equations are derived; the integration method and the spatial and temporal grid embedding techniques are explained; and results for sample problems involving stream-tube flow with one dissociating gas, shock-tube flow with one dissociating gas, and diverging-channel flow with multiple reactions are presented in extensive graphs and briefly characterized. Accuracy comparable to that of globally fine grid solutions is obtained with significant savings in CPU time.

Pervaiz, Mehtab M.↗
NASA NTRS1987

Particle paths and phase plane for time-dependent similarity solutions of the one-dimensional Vlasov-Maxwell equations

The phase trajectories of particles in a plasma described by the one-dimensional Vlasov-Maxwell equations are determined qualitatively, analyzing exact general similarity solutions for the cases of temporally damped and growing (sinusoidal or localized) electric fields. The results of numerical integration in both untransformed and Lie-group point-transformed coordinates are presented in extensive graphs and characterized in detail. The implications of the present analysis for the stability of BGK equilibria are explored, and the existence of nonlinear solutions arbitrarily close to and significantly different from the BGK solutions is demonstrated.

Roberts, Dana Aaron↗
NASA NTRS1988

Outcomes of planetary close encounters - A systematic comparison of methodologies

Several methods for estimating the outcomes of close planetary encounters are compared on the basis of the numerical integration of a range of encounter types. An attempt is made to lay the foundation for the development of predictive rules concerning the encounter outcomes applicable to the refinement of the statistical mechanics that apply to planet-formation and similar problems concerning planetary swarms. Attention is given to Oepik's (1976) formulation of the two-body approximation, whose predicted motion differs from the correct three-body behavior.

Greenberg, Richard↗