Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time integrator”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

Estimation of elastic aircraft parameters using the maximum likelihood method

The application of the maximum likelihood method to estimate the aerodynamic parameters of elastic flight vehicles in a symmetric flight condition is discussed. In this application, particular attention is directed toward the center of mass, elastic deformation, and sensor equations of motion. It is shown that the two major computational problems to be overcome are the inversion of large-sized matrices and the time-wise integration of a large number of linear, ordinary, differential equations.

Schwanz, R. C.↗
NASA NTRS1974

Calculation of turbulent shear stress in supersonic boundary layer flows

Turbulent shear stress distributions for supersonic boundary layer flows have been computed from experimental mean boundary layer data. The computations have been made by numerically integrating the time averaged continuity and streamwise momentum equations. Distributions have been obtained for flows upstream and downstream of shock-wave-boundary layer interactions and for both two-dimensional and axisymmetric flows. The computed results are compared with recently reported shear stress measurements which were obtained by hot wire anemometer and laser velocimeter techniques.

Sun, C. C.↗
NASA NTRS1973

Initiation and structure of axisymmetric eddies in a rotating stream

Axisymmetric flow of a rotating stream is examined numerically to determine conditions under which an isolated eddy will form on the axis of rotation. An explicit finite difference procedure is used to integrate the time dependent transport equations. Solutions provide details of the flow structure and are presented for a range of Reynolds numbers and swirl ratios. Calculated results are interpreted in terms of recent physical experiments insofar as is possible.

Kopecky, R. M.↗
NASA NTRS1974

Steady-state and dynamic analysis of a jet engine, gas lubricated shaft seal

Dynamic response of a gas-lubricated, jet-engine main shaft seal was analytically established as a function of collar misalignment and secondary seal friction. Response was obtained by a forward integration-in-time (time-transient) scheme, which traces a time history of seal motions in all its degrees of freedom. Results were summarized in the form of a seal tracking map which indicated regions of acceptable collar misalignments and secondary seal friction. Methodology, results and interpretations are comprehensively described.

Shapiro, W.↗
NASA NTRS1979

Evaluation of MOSTAS computer code for predicting dynamic loads in two-bladed wind turbines

Calculated dynamic blade loads are compared with measured loads over a range of yaw stiffnesses of the DOE/NASA Mod-0 wind turbine to evaluate the performance of two versions of the MOSTAS computer code. The first version uses a time-averaged coefficient approximation in conjunction with a multiblade coordinate transformation for two-bladed rotors to solve the equations of motion by standard eigenanalysis. The results obtained with this approximate analysis do not agree with dynamic blade load amplifications at or close to resonance conditions. The results of the second version, which accounts for periodic coefficients while solving the equations by a time history integration, compare well with the measured data.

Kaza, K. R. V.↗
NASA NTRS1979

Evaluation of improved algorithms for nonlinear transient structural heat transfer

The paper presents the evaluation of algorithms having potential for reducing computer time to integrate large systems of equations governing nonlinear transient heat transfer in complex structures. A series of test problems is being developed to test two algorithms. The Shuttle Orbiter frame model having 190 nodes was the subject of an analytical/experimental study. Plans are to perform calculations for the complete wing model having about 4000 nodes as well as a regional (plug) model having about 200 nodes.

Adelman, H. M.↗
NASA NTRS1979

Calculated global conditions for the unsteady Euler's flow equations

This paper considers the numerical solution of a fluid dynamical system with time-dependent oscillating boundaries. The problem detailed, a compressible and inviscid flow in a one-dimensional pipe with vibrating ends, is examined to develop techniques needed to solve the more complex unsteady flows in turbomachine cascades. The problem is shown to be mathematically ill-posed but by complementing the governing equations with energy conservation, a unique solution is obtained. An artificial damping term is introduced into the governing equations to improve the performance of the time-marching integration.

Kuzo, D. M.↗
NASA NTRS1981

A hybrid method for improving MCSCF convergence

It has been found that the convergence problems for many ill conditioned single-configuration SCF calculations arise from mixing among only a small number of orbitals. This orbital set includes the highest closed, the partially filled, and (possibly) a few of the lowest virtual orbitals. The fact that convergence problems can be, in very large measure, linked to a small orbital set is used to design a hybrid MCSCF procedure in which these orbitals are treated using a second-order MCSCF method, while other mixings are treated with a lower-order method which avoids the time consuming integral transformation. Tests on BeO show that the hybrid method yields convergence even when the simple lower-order treatment diverges. The method is expected to facilitate determination of MCSCF wave functions for large basis problems and for the construction of potential energy surfaces.

Bauschlicher, C. W., Jr.↗
NASA NTRS1981

Superflux I, II, and III experiment designs: Water sampling and analyses

Both airborne remote sensors and seagoing research vessels were used to study the effects of man's continual use of the Chesapeake Bay offshore environments. The major focus of the study was to: (1) advance the development and transfer of improved remote sensing systems and techniques for monitoring environmental quality and effects on living marine resources; (2) increase understanding of the influence of estuarine outwellings (plumes) on contiguous shelf ecosystems; and (3) provide a synoptic, integrated and timely data base for application to problems of marine resources and environmental quality.

Thomas, J. P.↗
NASA NTRS1982

Three-dimensional relativistic field-electron interaction in a multicavity high-power klystron. Part 2: Working Equations

A high power multicavity klystron amplifier was designed and a computation package containing all equations and procedures needed is presented. The rigorously derived three dimensional relativistic axisymmetric equations of motion are used to compute the bunched current and the induced RF gap voltage for all interaction cavities except the input and second cavities, where the linear space charge wave theory data are employed in order to reduce the computation time. Both distance step and time step integration methods are used to compute the Fourier coefficients of both the beam current and induced current.

Kosmahl, H. G.↗
NASA NTRS1982

An experimental investigation of pure tone generation by vortex shedding in a duct

An experimental investigation has been carried out for acoustic oscillations sustained by flow through a duct containing two baffles. Pure acoustic tones corresponding to longitudinal resonant modes of the duct are produced when certain flow and geometrical conditions are satisfied. The conditions are such as to ensure close coincidence between the frequency of vortex shedding from the forward baffle, and a natural frequency of the duct. Flow visualization has shown that under these conditions a stable vortex structure exists between the baffles, containing at all times an integral number of vortices.

Nomoto, H.↗
NASA NTRS1982

A numerical simulation of an atmospheric vortex street

A simplified mesoscale mixed-layer model is utilized to study the flow past an ocean island during a case of cold air outbreak over the Kuroshio Current. Using data taken during AMTEX '75, the governing equations are integrated in time to simulate the development of a Karman vortex street downstream of the island of Cheju-do. The surface layer, inversion layer, and overlying stable layer are all parametrized in the model. Very weak sinusoidal waviness but no vortex street was found in the case of the island. For a mixed layer barely deep enough to allow no part of the mountain on the island to protrude, there was more pronounced waviness but still no vortices. This result may not hold for boundary layers with poor vertical mixing. A numerical simulation of the vortex street yielded good results in comparison. The speed ratio, spacing ratio, and Strouhal number agreed rather closely with observations.

Ruscher, P. H.↗
NASA NTRS1984

Comparison of numerical techniques for integration of stiff ordinary differential equations arising in combustion chemistry

The efficiency and accuracy of several algorithms recently developed for the efficient numerical integration of stiff ordinary differential equations are compared. The methods examined include two general-purpose codes, EPISODE and LSODE, and three codes (CHEMEQ, CREK1D, and GCKP84) developed specifically to integrate chemical kinetic rate equations. The codes are applied to two test problems drawn from combustion kinetics. The comparisons show that LSODE is the fastest code currently available for the integration of combustion kinetic rate equations. An important finding is that an interactive solution of the algebraic energy conservation equation to compute the temperature does not result in significant errors. In addition, this method is more efficient than evaluating the temperature by integrating its time derivative. Significant reductions in computational work are realized by updating the rate constants (k = at(supra N) N exp(-E/RT) only when the temperature change exceeds an amount delta T that is problem dependent. An approximate expression for the automatic evaluation of delta T is derived and is shown to result in increased efficiency.

Radhakrishnan, K.↗
NASA NTRS1985

New integration techniques for chemical kinetic rate equations. 2: Accuracy comparison

A comparison of the accuracy of several techniques recently developed for solving stiff differential equations is presented. The techniques examined include two general purpose codes EEPISODE and LSODE developed for an arbitrary system of ordinary differential equations, and three specialized codes CHEMEQ, CREKID, and GCKP84 developed specifically to solve chemical kinetic rate equations. The accuracy comparisons are made by applying these solution procedures to two practical combustion kinetics problems. Both problems describe adiabatic, homogeneous, gas phase chemical reactions at constant pressure, and include all three combustion regimes: induction, heat release, and equilibration. The comparisons show that LSODE is the most efficient code - in the sense that it requires the least computational work to attain a specified accuracy level. An important finding is that an iterative solution of the algebraic enthalpy conservation equation for the temperature can be more accurate and efficient than computing the temperature by integrating its time derivative.

Radhakrishnan, K.↗
NASA NTRS1985

Gamma rays of energy or = 10(15) eV from Cyg X-3

The experimental data of extensive air showers observed at Akeno have been analyzed to detect the gamma ray signal from Cyg X-3. After muon poor air showers are selected, the correlation of data acquisition time with 4.8 hours X-ray period is studied, giving the data concentration near the phase 0.6, the time of X-ray maximum. The probability that uniform backgrounds create the distribution is 0.2%. The time averaged integral gamma ray flux is estimated as (1.1 + or - 0.4)x 10 to the -14th power cm(-2) sec(-1) for Eo 10 to the 15th power eV and (8.8 + or - 5.0)x 10 to the 14th power cm(-2) sec(-1) for Eo 6 x 10 to the 14th power eV.

Kifune, T.↗
NASA NTRS1985

Electromagnetic particle simulation codes

Electromagnetic particle simulations solve the full set of Maxwell's equations. They thus include the effects of self-consistent electric and magnetic fields, magnetic induction, and electromagnetic radiation. The algorithms for an electromagnetic code which works directly with the electric and magnetic fields are described. The fields and current are separated into transverse and longitudinal components. The transverse E and B fields are integrated in time using a leapfrog scheme applied to the Fourier components. The particle pushing is performed via the relativistic Lorentz force equation for the particle momentum. As an example, simulation results are presented for the electron cyclotron maser instability which illustrate the importance of relativistic effects on the wave-particle resonance condition and on wave dispersion.

Pritchett, P. L.↗
NASA NTRS1986

New integration techniques for chemical kinetic rate equations. II - Accuracy comparison

A comparison of the accuracy of several techniques recently developed for solving stiff differential equations is presented. The techniques examined include two general purpose codes EEPISODE and LSODE developed for an arbitrary system of ordinary differential equations, and three specialized codes CHEMEQ, CREKID, and GCKP84 developed specifically to solve chemical kinetic rate equations. The accuracy comparisons are made by applying these solution procedures to two practical combustion kinetics problems. Both problems describe adiabatic, homogeneous, gas phase chemical reactions at constant pressure, and include all three combustion regimes: induction heat release, and equilibration. The comparisons show that LSODE is the most efficient code - in the sense that it requires the least computational work to attain a specified accuracy level. An important finding is that an iterative solution of the algebraic enthalpy conservation equation for the temperature can be more accurate and efficient than computing the temperature by integrating its time derivative.

Radhakrishnan, K.↗
NASA NTRS1988

An explicit Runge-Kutta method for unsteady rotor/stator interaction

A quasi-three-dimensional rotor/stator analysis has been developed for blade-to-blade flows in turbomachinery. The analysis solves the unsteady Euler or thin-layer Navier-Stokes equations in a body-fitted coordinate system. It accounts for the effects of rotation, radius change, and stream-surface thickness. The Baldwin-Lomax eddy-viscosity model is used for turbulent flows. The equations are integrated in time using a four-stage Runge-Kutta scheme with a constant timestep. Results are shown for the first stage of the Space Shuttle Main Engine high pressure fuel turbopump. Euler and Navier-Stokes results are compared on the scaled single- and multi-passage machine. The method is relatively fast and the quasi-three-dimensional formulation is applicable to a wide range of turbomachinery geometries.

Jorgenson, Philip C. E.↗