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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
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.↗
NASA NTRS1975

Magnetic attitude control system for dual-spin satellites

A closed-loop control law is developed for a dual-spin satellite control system which utilizes the interaction of the geomagnetic field with the satellite dipole parallel to the spin axis. The control law consists of the linear combination of the pitch axis component of the rate of change of the geomagnetic field and the product of the roll angle and roll axis component of the geomagnetic field. Application of the method of multiple time scales yields approximate solutions for the feedback gains in terms of the system parameters. Approximate solutions are also obtained for the response of the system to disturbance torques. A comparison of the approximate solutions and numerical solutions obtained by numerical integration of the exact equations of motion is then given.

Alfriend, K. T.↗
NASA NTRS1976

Proton tissue dose

Program calculates dosage averaged over five major segments of blood-forming organ, treating human body geometry in detail but assuming isotropicity of incident primary particles. Approximate form of transport theory is used, incorporating nuclear star effects. Two numerical integrations are used to evaluate intermediate equation and then dosage equation.

Wilson, J. W.↗
NASA NTRS1975

Stabilization by modification of the Lagrangian

In order to reduce the error growth during a numerical integration, a method of stabilization of the differential equations of the Keplerian motion is offered. It is characterized by the use of the eccentric anomaly as an independent variable in such a way that the time transformation is given by a generalized Lagrange formalism. The control terms in the equations of motion obtained by this modified Lagrangian give immediately a completely Liapunov-stable set of differential equations. In contrast to other publications, here the equation of time integration is modified by a control term which leads to an integral which defined the time element for the perturbed Keplerian motion.

Baumgarte, J. W.↗
NASA NTRS1975

Mean colors of Cepheids

Color observations of Cepheids can be interpreted by a color effective temperature relation to give different T sub e values for each kind of mean color. Here, theoretical colors from numerical integrations of Cepheid pulsations are used to determine the proper method of taking the color mean in order to get, by the color T sub e relation, the correct nonpulsating T sub e.

Cox, A. N.↗