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

Zonal gravity harmonics from long satellite arcs by a seminumeric method.

A zonal geopotential is presented to degree 21 from evaluation of mean elements for 21 satellites, including two of low (less than 20 deg) inclination. Each satellite is represented by an arc of at least one apsidal rotation. The lengths range from 50 to 660 days. Differential correction of the initial elements in all of the arcs, together with radiation pressure and atmospheric drag coefficients, is accomplished simultaneously with the correction for the zonal harmonics. The satellite orbits and their variations are generated by numerical integration of the Lagrange equations for mean elements. Disturbances due to precession and nutation, atmospheric drag, radiation pressure, and lunisolar gravity are added at 1- to 8-day intervals in the integrated orbits.

Wagner, C. A.↗
NASA NTRS1973

Bayesian estimation of life parameters in the Weibull distribution.

Development of a Bayesian analysis of the scale and shape parameters in the Weibull distribution and the corresponding reliability function with respect to the usual life-testing procedures. For the scale parameter theta, Bayesian estimates of theta and reliability are obtained for the uniform, exponential, and inverted gamma prior probability densities. Bhattacharya's results (1967) for the one-parameter exponential life-testing distribution are reduced to a special case of these results. A fully Bayesian analysis of both the scale and shape parameters is developed by assuming independent prior distributions; since in the latter case, analytical tractability is not possible, Bayesian estimates are obtained through a conjunction of Monte Carlo simulation and numerical-integration techniques. In both cases, a computer simulation is carried out, and a comparison is made between the Bayesian and the corresponding minimum-variance unbiased, or maximum likelihood, estimates. As expected, the Bayesian estimates are superior.

Canavos, G. C.↗
NASA NTRS1973

Quarantine constraints as applied to satellites.

Plans for unmanned missions to planets beyond Mars in the 1970s include satellite encounters. Recently published observations of data for Titan, a satellite of Saturn, indicate that conditions may be hospitable for the growth of microorganisms. Therefore, the problem of satisfying possible quarantine constraints for outer planet satellites was investigated. This involved determining the probability of impacting a satellite of Jupiter or Saturn by a spacecraft for a planned satellite encounter during an outer planet mission. Mathematical procedures were formulated which (1) determine the areas in the aim-plane that would result in trajectories that impact the satellite and (2) provide a technique for numerically integrating the navigation error function over the impact area to obtain impact probabilities. The results indicate which of the planned spacecraft trajectory correction maneuvers are most critical in terms of satellite quarantine violation.

Hoffman, A. R.↗
NASA NTRS1973

Ion kinematics in a plasma focus.

The results of numerical integrations of three-dimensional equations of motion of ions subject to given electric and magnetic fields are presented. The fields represent those which may exist in the pinch phase of the plasma focus, although here they depend only on the radial coordinate. The ions initially have Maxwellian velocity distributions, and their trajectories are interpreted in terms of single-particle constants of the motion. Two models of the axial electric field Ez are considered. For strong Ez away from the axis, there is a cyclotron acceleration which leads to ion heating. For positive Ez on the axis, ions within a Larmor radius of the axis undergo very efficient acceleration; the results for this new model are in general agreement with experimental results.

Gary, S. P.↗
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.↗