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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

Ion heating in a plasma focus

Ion acceleration and heating in a plasma focus were investigated by the numerical integration of the three-dimensional equations of motion. The electric and magnetic fields given were derived from experimental data. The results obtained show that during the collapse phase of focus formation, ions are efficiently heated to temperatures of several keV. During the phase of rapid current reduction, ions are accelerated to large velocities in the axial direction. The results obtained with the model are in general agreement with experimental results.

Hohl, F.↗
NASA NTRS1974

Laser Doppler velocimeter system simulation for sensing aircraft wake vortices

A hydrodynamic model of aircraft vortex wakes in an irregular wind shear field near the ground is developed and used as a basis for modeling the characteristics of a laser Doppler detection and vortex location system. The trailing vortex sheet and the wind shear are represented by discrete free vortices distributed over a two-dimensional grid. The time dependent hydrodynamic equations are solved by direct numerical integration in the Boussinesq approximation. The ground boundary is simulated by images, and fast Fourier Transform techniques are used to evaluate the vorticity stream function. The atmospheric turbulence was simulated by constructing specific realizations at time equal to zero, assuming that Kolmogoroff's law applies, and that the dissipation rate is constant throughout the flow field. The response of a simulated laser Doppler velocimeter is analyzed by simulating the signal return from the flow field as sensed by a simulation of the optical/electronic system.

Thomson, J. A. L.↗
NASA NTRS1973

Dynamic response and stability of a gas-lubricated Rayleigh-step pad

The quasi-static, pressure characteristics of a gas-lubricated thrust bearing with shrouded, Rayleigh-step pads are determined for a time-varying film thickness. The axial response of the thrust bearing to an axial forcing function or an axial rotor disturbance is investigated by treating the gas film as a spring having nonlinear restoring and damping forces. These forces are related to the film thickness by a power relation. The nonlinear equation of motion in the axial mode is solved by the Ritz-Galerkin method as well as the direct, numerical integration. Results of the nonlinear response by both methods are compared with the response based on the linearized equation. Further, the gas-film instability of an infinitely wide Rayleigh step thrust pad is determined by solving the transient Reynolds equation coupled with the equation of the motion of the pad. Results show that the Rayleigh-step geometry is very stable for bearing number A up to 50. The stability threshold is shown to exist only for ultrahigh values of Lambda equal to or greater than 100, where the stability can be achieved by making the mass heavier than the critical mass.

Cheng, C.↗
NASA NTRS1974

On the theory of thermal runaway in the hydrogen envelope of a white dwarf

An evolutionary sequence was calculated for a white dwarf on which a hydrogen rich envelope is assumed to increase with time. The stellar models were computed by integrating numerically the system of stellar structure equations by means of the Henyey method. The accretion of matter was assumed to be quasi-static. Near the point of the chemical discontinuity, a temperature maximum is shown to arise, and hydrogen ignites; a thin shell energy source is formed which is thermally unstable. The resulting thermal runaway was investigated numerically and by means of a simple approximation based on the energy balance in the hydrogen burning shell source.

Redkoborodyy, Y. N.↗
NASA NTRS1975

Some procedures for computerized electronic data processing of absorption measurements from artificial earth satellites

The processing of data obtained from solar absorption radiation measurements is discussed. The position of the satellite was obtained by numerical integration of the differential equations of motion using initial conditions. The position of the sun was calculated as a function of time, and the tangential elevation was determined approximately from the positions of the satellite and the sun. The coefficients of an approximation formula and of a data smoothing process were determined, and the inversion of an Abel integral equation is solved analytically.

Trinkkeller, B.↗
NASA NTRS1975

Analytical method for calculation of navigational data for the position of a satellite

A method is described for calculating the position of a satellite at the instants when measurements are made on board. The initial conditions used were the mean orbital elements of the satellite and their time derivatives in one orbit. The results of the calculation are compared with those obtained by numerical integration, and it is found that results are identical at the beginning of an orbit, but change as the orbit progresses. The advantages and disadvantages of the analytical method are presented.

Lala, P.↗
NASA NTRS1975

The NASTRAN SRB slapdown water impact analysis

A computer analysis was made of the hydroelastic interaction between a shell structure representing the space shuttle solid rocket motor case and an incompressible fluid during the slapdown phase of water impact. The large motions and hydroelastic response of the system were obtained by numerical integration of the combined hydrodynamics and structural equations of motion and are analyzed for a selected interval after the initial nozzle impact. The computerization of the slapdown hydroelastic capability was incorporated into the general purpose NASTRAN computer code. The development of the theoretical basis, a guide to the program's usage, results of correlation and parameter studies, and a detailed description of the computer code are also included.

Herting, D. N.↗
NASA NTRS1975

Prediction and verification of creep behavior in metallic materials and components for the space shuttle thermal protection system. Volume 2: Phase 2 subsize panel cyclic creep predictions

A method for predicting permanent cyclic creep deflections in stiffened panel structures was developed. The resulting computer program may be applied to either the time-hardening or strain-hardening theories of creep accumulation. Iterative techniques were used to determine structural rotations, creep strains, and stresses as a function of time. Deflections were determined by numerical integration of structural rotations along the panel length. The analytical approach was developed for analyzing thin-gage entry vehicle metallic-thermal-protection system panels subjected to cyclic bending loads at high temperatures, but may be applied to any panel subjected to bending loads. Predicted panel creep deflections were compared with results from cyclic tests of subsize corrugation and rib-stiffened panels. Empirical equations were developed for each material based on correlation with tensile cyclic creep data and both the subsize panels and tensile specimens were fabricated from the same sheet material. For Vol. 1, see N75-21431.

Cramer, B. A.↗
NASA NTRS1975

On upstream blocking in a viscous diffusive stratified flow

The effect of diffusion of specie upon the flow about a transverse flat plate moving horizontally in a viscous stratified medium is considered. Asymptotic expansions are used to define a parameter regime where a viscous-diffusive-buoyancy balance is dominant. The solution, expressed in terms of an inverse Fourier transform, is numerically integrated. The results show that, as in the non-diffusive problem, a region of closed streamlines exists ahead of the body. However, unlike the case where diffusion is neglected, the density field within this recirculating region is uniquely determined and found to be statically stable. It is also found that varying the relative amount of diffusion affects not only the density distribution, but the velocity profile as well, indicating a strong coupling between the vorticity and specie equation.

Koop, C. G.↗
NASA NTRS1975

ISOFINEL: Isoparametric finite element code for elastic analysis of two-dimensional bodies

A formulation is presented for the development of a finite element program for the elastic analysis of two-dimensional bodies using the eight-node isoparametric quadrilateral. The program solves for both plane stress and plane strain problems. The finite element formulation based on the isoparametric displacement functions is presented. The program structure is given in the form of flow diagrams with descriptions of the numerical procedure used to obtain the element stiffness matrix, and the solution method employed to solve for nodal displacements. Three numerical examples (a plate under uniaxial tension, a plate under pure shear, and a beam under pure bending) are presented to illustrate the capability and limitations of the element implementation. The first problem is solved exactly by the element, as predicted by the form of its displacement functions. In the other two problems the accuracy of the solution is highly dependent upon the slenderness of the element, the number of elements in the map, and the numerical integration scheme used to build the element stiffness matrix.

Marino, C.↗
NASA NTRS1975

Nonlinear seismic analysis of a reactor structure impact between core components

The seismic analysis of the FFTF-PIOTA (Fast Flux Test Facility-Postirradiation Open Test Assembly), subjected to a horizontal DBE (Design Base Earthquake) is presented. The PIOTA is the first in a set of open test assemblies to be designed for the FFTF. Employing the direct method of transient analysis, the governing differential equations describing the motion of the system are set up directly and are implicitly integrated numerically in time. A simple lumped-nass beam model of the FFTF which includes small clearances between core components is used as a "driver" for a fine mesh model of the PIOTA. The nonlinear forces due to the impact of the core components and their effect on the PIOTA are computed.

Hill, R. G.↗
NASA NTRS1975

Evaluation of geopotential and luni-solar perturbations by a recursive algorithm

The disturbing functions due to the geopotential and Luni-solar attractions are linear and bilinear forms in spherical harmonics. Making use of recurrence relations for the solid spherical harmonics and their derivatives, recurrence formulas are obtained for high degree terms as function of lower degree for any term of those disturbing functions and their derivative with respect to any element. The equations obtained are effective when a numerical integration of the equations of motion is appropriate. In analytical theories, they provide a fast way of obtaining high degree terms starting from initial very simple functions.

Giacaglia, G. E. O.↗
NASA NTRS1975

Programs for computing abscissas and weights for classical and nonclassical Gaussian quadrature formulas

Computer programs for computing Gaussian quadrature abscissas and weights are described. For the classical case the programs use Laguerre iteration to compute abscissas as zeros of orthogonal polynomials. The polynomials are evaluated from known recursion coefficients. The nonclassical case is handled similarly except that the recursion coefficients are computed by numerical integration. A sample problem, with input and output, is presented to illustrate the use of the programs. It computes the quadrature abscissas and weights associated with the weight function over the interval (0,1) for quadrature orders from 16 to 96 in increments of 8.

Desmarais, R. N.↗
NASA NTRS1975

Modern control concepts in hydrology

Two approaches to an identification problem in hydrology are presented, based upon concepts from modern control and estimation theory. The first approach treats the identification of unknown parameters in a hydrologic system subject to noisy inputs as an adaptive linear stochastic control problem; the second approach alters the model equation to account for the random part in the inputs, and then uses a nonlinear estimation scheme to estimate the unknown parameters. Both approaches use state-space concepts. The identification schemes are sequential and adaptive and can handle either time-invariant or time-dependent parameters. They are used to identify parameters in the Prasad model of rainfall-runoff. The results obtained are encouraging and confirm the results from two previous studies; the first using numerical integration of the model equation along with a trial-and-error procedure, and the second using a quasi-linearization technique. The proposed approaches offer a systematic way of analyzing the rainfall-runoff process when the input data are imbedded in noise.

Duong, N.↗
NASA NTRS1975

Short-periodic variations and second-order numerical averaging

The principal disadvantage of the method of numerical averaging is that it provides only the average time history of the orbital elements and yields no information about the high-frequency or short-periodic variations that occur inside the averaging interval. This paper contains a description of a technique for recovering the short-periodic variations by minor modifications to the averaging process so as to permit the construction of a Fourier series for the osculating elements. The availability of this series permits the extension of the averaging technique to higher order and allows us to account for short-periodic coupling of the perturbations. Comparisons of the results with numerically integrated solutions are presented for three distinct orbit prediction problems.

Lutzky, D.↗
NASA NTRS1975

Creep analysis capability for stiffened panel structures

A method of analysis developed for predicting permanent cyclic creep deflections in stiffened panel structures is presented. The resulting computer program (TPSC, for Thermal Protection System Creep) applies either the time hardening or strain hardening theories of creep accumulation, using iterative techniques to determine structural rotations, creep strains, and stresses as a function of time. Deflections are determined by numerical integration of structural rotations along the panel length. The analytical approach was developed for purposes of analyzing thin gage entry vehicle metallic thermal protection system panels subjected to cyclic bending loads at high temperatures. The method is applicable to any panel subjected to bending loads. Predicted panel creep deflections are compared with results from cyclic tests of subsize and full size corrugation and rib stiffened panels. Empirical equations used in the analysis were developed for each material based on correlation with tensile cyclic creep data.

Cramer, B. A.↗
NASA NTRS1975

Time elements

Time elements are introduced for use with Sundman time transformations for satellite equations of motion. Two time elements are presented, one providing maximum accuracy when the exponent of the governing equation equals 1, the other when the exponent equals 2. Time elements and time transformations reduce local truncation error and Liapunov (in-track) instability, and provide analytical step size control. Numerical results show accuracy improvements of more than one order of magnitude when time elements are employed with time transformations in the numerical integration of the satellite equations, compared with using time transformations alone.

Nacozy, P. E.↗
NASA NTRS1975

Simulation of multicorrelated random processes using the FFT algorithm

A technique for the digital simulation of multicorrelated Gaussian random processes is described. This technique is based upon generating discrete frequency functions which correspond to the Fourier transform of the desired random processes, and then using the fast Fourier transform (FFT) algorithm to obtain the actual random processes. The main advantage of this method of simulation over other methods is computation time; it appears to be more than an order of magnitude faster than present methods of simulation. One of the main uses of multicorrelated simulated random processes is in solving nonlinear random vibration problems by numerical integration of the governing differential equations. The response of a nonlinear string to a distributed noise input is presented as an example.

Wittig, L. E.↗