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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
NASA NTRS1973

Prior to failure extension of flaws in a rate sensitive Tresca solid.

Theoretical analysis shows that the amount of slow growth occurring prior to fracture in a plane stress tension specimen subjected to a subcritical stress intensity level is affected by the ductility and rheological sensitivity of the material, the rate of loading, the initial crack size, and the geometrical configuration of the test. These parameters enter the governing nonlinear differential equation which describes the quasi-static crack extension in a rate-sensitive Tresca solid. Examples of numerical integration are given and employed to predict the maximum fracture toughness attained in a plane stress tension test.

Wnuk, M. P.↗
NASA NTRS1973

Orbit improvement from satellite imaging data obtainable from outer planet missions.

Equations of motion are established for a dynamical system in which a spacecraft flies close to and interacts with an outer planet and one or more of its satellites. For the computation of the state and mass partials needed in a simultaneous orbit correction of n interacting bodies, a notably compact set of variational equations is derived. The above system of differential equations is integrated numerically on a computer. Spacecraft-satellite direction measurements accurate to plus or minus 10 sec were simulated along three representative trajectories (Mariner/Jupiter/Saturn 1977 missions) approaching Io, Titan, and Iapetus to within 41000, 13000, and 7000 km, respectively. The paper concludes with a brief discussion of the need for future work on the orbits of the satellites of the outer planets.

Aksnes, K.↗
NASA NTRS1973

Acoustic wave propagation in an axisymmetric swirling jet.

An analysis has been developed to study the acoustic wave propagation in an axisymmetric swirling subsonic jet flow. The governing convected wave equation derived in the spherical coordinates includes mean shears, shear gradients and pressure gradients. The directivity patterns for various spinning and non-spinning modes due to the influence of the mean jet swirl were obtained by numerically integrating the governing wave equation. The mean flow field used in the computation was that obtained semiempirically for subsonic swirling turbulent jet and is completely specified once the degree of swirl is known. The dependence of sound directivity on jet Mach number, swirl ratio and frequency are discussed.

Yu, J. C.↗
NASA NTRS1973

Large amplitude flutter of a low aspect ratio panel at low supersonic speeds comparison of theory and experiment

Flutter boundaries, as well as flutter limit cycle amplitudes, frequencies and stresses were computed for a panel of length-width ratio 4.48 exposed to applied in-plane and transverse loads. The Mach number range was 1.1 to 1.4. The method used involved direct numerical integration of modal equations of motion derived from the nonlinear plate equations of von Karman, coupled with linearized potential flow aerodynamic theory. The flutter boundaries agreed reasonably well with experiment, except when the in-plane loading approached the buckling load. Structural damping had to be introduced, to produce frequencies comparable to the experimental values. Attempts to compute panel deflections or stress at a given point met with limited success. There is some evidence, however, that deflection and stress maxima can be estimated with somewhat greater accuracy.

Ventres, C. S.↗
NASA NTRS1973

Three isoparametric solid elements for NASTRAN

Linear, quadratic, and cubic isoparametric hexahedral solid elements have been added to the element library of NASTRAN. These elements are available for static, dynamic, buckling, and heat-transfer analyses. Because the isoparametric element matrices are generated by direct numerical integration over the volume of the element, variations in material properties, temperatures, and stresses within the elements are represented in the computations. In order to compare the accuracy of the new elements, three similar models of a slender cantilever were developed, one for each element. All elements performed well. As expected, however, the linear element model yielded excellent results only when shear behavior predominated. In contrast, the results obtained from the quadratic and cubic element models were excellent in both shear and bending.

Johnson, S. E.↗
NASA NTRS1973

The early history of the lunar inclination

The effect of tidal friction on the inclination of the lunar orbit to the earth's equator for earth-moon distances of less than 10 earth radii is examined. The results obtained bear on a conclusion drawn by Gerstenkorn and others which has been raised as a fatal objection to the fission hypothesis of lunar origin, namely, that the present nonzero inclination of the moon's orbit to the ecliptic implies a steep inclination of the moon's orbit to the earth's equatorial plane in the early history of the earth-moon system. This conclusion is shown to be valid only for particular rheological models of the earth. The earth is assumed to behave like a highly viscous fluid in response to tides raised in it by the moon. The moon is assumed to be tideless and in a circular orbit about the earth. The equations of tidal friction are integrated numerically to give inclination of the lunar orbit as a function of earth-moon distance.

Rubincam, D. P.↗
NASA NTRS1974

Dynamic analysis of a system of hinge-connected rigid bodies with nonrigid appendages

Equations of motion are derived for use in simulating a spacecraft or other complex electromechanical system amenable to idealization as a set of hinge-connected rigid bodies of tree topology, with rigid axisymmetric rotors and nonrigid appendages attached to each rigid body in the set. In conjunction with a previously published report on finite-element appendage vibration equations, this report provides a complete minimum-dimension formulation suitable for generic programming for digital computer numerical integration.

Likins, P. W.↗
NASA NTRS1974

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 conform with results from two previous studies; the first using numerical integration of the model equation along with a trial-and-error procedure, and the second, by 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 NTRS1973

Asymptotic solution to the tangential low thrust energy increase trajectory

A approximate analytic solution to the low thrust constant acceleration energy increase trajectory using a tangential steering program is presented. The solution is based on a two-variable asymptotic expansion of the equations of motion formulated in regularized variables. The high accuracy of the second order solution derived is demonstrated by comparing it with numerically integrated trajectories.

Schwenzfeger, K. J.↗
NASA NTRS1974

A recurrence matrix method for the analysis of longitudinal and torsional vibrations in non-uniform multibranch beams with variable boundary conditions

An approximate method for calculating the longitudinal and torsional natural frequencies and associated modal data of a beamlike, variable cross section multibranch structure is presented. The procedure described is the numerical integration of the first order differential equations that characterize the beam element in longitudinal motion and that satisfy the appropriate boundary conditions.

Davis, R. B.↗
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 NTRS1974

Investigation of highly efficient satellite solution methods

Methods for analyzing the stability of satellites are discussed. The subjects considered are: (1) time elements, (2) stabilization by external energy corrections, and (3) long term global solutions for the synchronous satellite. A set of canonical two-body elements referred to as Delaunay-similar elements is presented. In contrast to the classical Delaunay theory which has time as the independent variable, the D-S theory uses an independent variable which is a generalized true anomaly. The numerical integration of the canonical perturbation equations of these elements is developed.

Nacozy, P. E.↗
NASA NTRS1971

Numerical solution of stiff systems of ordinary differential equations with applications to electronic circuits

Systems of ordinary differential equations in which the magnitudes of the eigenvalues (or time constants) vary greatly are commonly called stiff. Such systems of equations arise in nuclear reactor kinetics, the flow of chemically reacting gas, dynamics, control theory, circuit analysis and other fields. The research reported develops an A-stable numerical integration technique for solving stiff systems of ordinary differential equations. The method, which is called the generalized trapezoidal rule, is a modification of the trapezoidal rule. However, the method is computationally more efficient than the trapezoidal rule when the solution of the almost-discontinuous segments is being calculated.

Rosenbaum, J. S.↗
NASA NTRS1974

Contribution to the theory of tidal oscillations of an elastic earth. External tidal potential

The differential equations of the tidal oscillations of the earth were established under the assumption that the interior of the earth is laterally inhomogeneous. The theory was developed using vectorial and dyadic symbolism to shorten the exposition and to reduce the differential equations to a symmetric form convenient for programming and for numerical integration. The formation of tidal buldges on the surfaces of discontinuity and the changes in the internal density produce small periodic variations in the exterior geopotential which are reflected in the motion of artificial satellites. The analoques of Love elastic parameters in the expansion of exterior tidal potential reflect the asymmetric and inhomogeneous structure of the interior of the earth.

Musen, P.↗
NASA NTRS1974

A study to determine the usefulness of interval analysis in solving problems in celestial mechanics

This investigation was undertaken to determine the usefulness of interval analysis to numerical integration and matrix inversion techniques and to combine these results to determine the value of interval analysis in bounding computational errors in the two-body problem. Conclusions were that interval analysis may be worthwhile in certain small scale isolated problems, but its usefulness in any large scale problem is doubtful.

Walling, D.↗
NASA NTRS1973

Low-thrust trajectory analysis for the geosynchronous mission

Methodology employed in development of a computer program designed to analyze optimal low-thrust trajectories is described, and application of the program to a Solar Electric Propulsion Stage (SEPS) geosynchronous mission is discussed. To avoid the zero inclination and eccentricity singularities which plague many small-force perturbation techniques, a special set of state variables (equinoctial) is used. Adjoint equations are derived for the minimum time problem and are also free from the singularities. Solutions to the state and adjoint equations are obtained by both orbit averaging and precision numerical integration; an evaluation of these approaches is made.

Jasper, T. P.↗
NASA NTRS1973

Quasi-static extension of a tensile crack contained in a viscoelastic-plastic solid

Final stretch criterion of failure is applied to the problem of quasi-static extension of a crack embedded in an elastic-plastic or viscoelastic-plastic matrix. The slow growth under subcritical conditions in a rate-sensitive Tresca solid is shown to be a superposition of creep rupture and McClintock's ductile growth. This type of growth occurs at subcritical magnitude of the imposed K-factor and can be accounted for only through a recognition of inelastic properties of solids. In the subcritical range there is no unique value for K sub c independent of geometrical configuration and flaw size. Not only the produced states of stress and strain are dependent on the loading path, but also the material resistance to fracture turns out to be a function of the history of loading that precedes catastrophic failure. A nonlinear integro-differential equation of motion is derived for a crack progressing through a viscoelastic medium with some limited ability to plastic flow. Examples of numerical integration are given incorporating both monotonic and cyclic loading programs.

Wnuk, M. P.↗