Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical techniques”

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 865 records · Page 48

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 determine the areas in the aim-plane that would result in trajectories that impact the satellite and 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.↗

Numerical integration routines for near-earth operations

Two general purpose numerical integration schemes were built into the NASA-JSC computer system. The state-of-the-art of numerical integration, the particular integrators built into the JSC computer system, and the use of the new integration packages are described. Background information about numerical integration and the variable-order, variable-stepsize Adams numerical integration technique is discussed. Results concerning the PEACE parameter optimization program are given along with recommendations and conclusions.

Powers, W. F.↗

A study of numerical methods of solution of the equations of motion of a controlled satellite under the influence of gravity gradient torque

Numerical methods of integration of the equations of motion of a controlled satellite under the influence of gravity-gradient torque are considered. The results of computer experimentation using a number of Runge-Kutta, multi-step, and extrapolation methods for the numerical integration of this differential system are presented, and particularly efficient methods are noted. A large bibliography of numerical methods for initial value problems for ordinary differential equations is presented, and a compilation of Runge-Kutta and multistep formulas is given. Less common numerical integration techniques from the literature are noted for further consideration.

Thompson, J. F.↗

Round-off errors in cutting plane algorithms based on the revised simplex procedure

This report statistically analyzes computational round-off errors associated with the cutting plane approach to solving linear integer programming problems. Cutting plane methods require that the inverse of a sequence of matrices be computed. The problem basically reduces to one of minimizing round-off errors in the sequence of inverses. Two procedures for minimizing this problem are presented, and their influence on error accumulation is statistically analyzed. One procedure employs a very small tolerance factor to round computed values to zero. The other procedure is a numerical analysis technique for reinverting or improving the approximate inverse of a matrix. The results indicated that round-off accumulation can be effectively minimized by employing a tolerance factor which reflects the number of significant digits carried for each calculation and by applying the reinversion procedure once to each computed inverse. If 18 significant digits plus an exponent are carried for each variable during computations, then a tolerance value of 0.1 x 10 to the minus 12th power is reasonable.

Moore, J. E.↗

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.↗

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.↗

Survey of computational methods for three-dimensional supersonic inviscid flows with shocks

Hyperbolic systems of partial differential equations governing supersonic inviscid flows are discussed and analyzed. Finite-difference analogues for integrating these systems in the interior of fluid domains are described from two points of view: a differential form approach and an integral form approach. The algorithms presented are analyzed for stability and accuracy. The concept of time splitting is discussed and supplied to these methods to achieve increased numerical efficiency. Techniques for treating conditions at the boundaries of the fluid domain and shock-wave discontinuities at surfaces within the domain are described.

Maccormack, R. W.↗

Systems identification - reprise and projections

A state-of-the-arts review is given for the field of system identification. Progress in the field is traced from the early models of dynamic systems by Sir Isaac Newton up to the present day use of advanced techniques for numerous applications.

Taylor, L. W., Jr.↗

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.↗

Combustion instability with partial length acoustic liners

The development of analytical techniques and numerical methods for the prediction of the stability behavior of liquid propellant rocket combustors is discussed. The studies involved nonlinear wave propagation effects and nonlinear combustion zone, liner, and nozzle responses. A second primary extension of the basic theory was directed at the inclusion of the effects of injector face baffles. Attention was restricted to the linear (small amplitude) problem, finite Mach numbers, three dimensional oscillations, and distributed combustion effects.

Mitchell, C. E.↗

Suppression of nonlinear oscillations in combustors with partial length acoustic liners

An analytical model is formulated for a three-dimensional nonlinear stability problem in a rocket motor combustion chamber. The chamber is modeled as a right circular cylinder with a short (multi-orifice) nozzle, and an acoustic linear covering an arbitrary portion of the cylindrical periphery. The combustion is concentrated at the injector and the gas flow field is characterized by a mean Mach number. The unsteady combustion processes are formulated using the Crocco time lag model. The resulting equations are solved using a Green's function method combined with numerical evaluation techniques. The influence of acoustic liners on the nonlinear waveforms is predicted. Nonlinear stability limits and regions where triggering is possible are also predicted for both lined and unlined combustors in terms of the combustion parameters.

Espander, W. R.↗

Turbulent flow between a rotating disk and a parallel wall

A scalar effective viscosity method is used to calculate the incompressible flow in the gap between two infinite, parallel smooth disks, one of them rotating and the other stationary. The flow is calculated for two values of the gap Reynolds number omega(s squared)/nu, namely 296 and 2852, where omega is the disk angular velocity, s the gap width, and nu the kinematic viscosity of the fluid; and over a range of radial Reynolds number omega(r squared)/nu from 0 to 10 million, where r is the radius from the axis of rotation. A numerical shooting technique is used to match the two boundary-layer flows arising on the disk and wall. Gradual transition is assumed for radial Reynolds numbers from 160,000 to 250,000. The results are in good agreement with Daily and Nece's measurements of velocity profiles and disk friction drag for enclosed rotating disks, and at sufficiently large radius the local skin friction approaches that of turbulent Couette flow.

Cooper, P.↗

Finite element techniques applied to cracks interacting with selected singularities

The finite-element method for computing the extensional stress-intensity factor for cracks approaching selected singularities of varied geometry is described. Stress-intensity factors are generated using both displacement and J-integral techniques, and numerical results are compared to those obtained experimentally in a photoelastic investigation. The selected singularities considered are a colinear crack, a circular penetration, and a notched circular penetration. Results indicate that singularities greatly influence the crack-tip stress-intensity factor as the crack approaches the singularity. In addition, the degree of influence can be regulated by varying the overall geometry of the singularity. Local changes in singularity geometry have little effect on the stress-intensity factor for the cases investigated.

Conway, J. C.↗

STICAP: A linear circuit analysis program with stiff systems capability. Volume 1: Theory manual

STICAP (Stiff Circuit Analysis Program) is a FORTRAN 4 computer program written for the CDC-6400-6600 computer series and SCOPE 3.0 operating system. It provides the circuit analyst a tool for automatically computing the transient responses and frequency responses of large linear time invariant networks, both stiff and nonstiff (algorithms and numerical integration techniques are described). The circuit description and user's program input language is engineer-oriented, making simple the task of using the program. Engineering theories underlying STICAP are examined. A user's manual is included which explains user interaction with the program and gives results of typical circuit design applications. Also, the program structure from a systems programmer's viewpoint is depicted and flow charts and other software documentation are given.

Cooke, C. H.↗

Preheated, combustion-driven gasdynamic lasers

The effect of preheated nitrogen on the available laser specific energy is examined for two classes of combustion-driven gasdynamic lasers (GDLs). The analysis combines a nonequilibrium model of the relevant GDL processes with a simplified nondissociating combustion model handling a variety of candidate GDL fuels. For a more comprehensive evaluation of the capabilities of preheating, an optimization routine is added to the combined laser/combustion models. The effectiveness of preheating to enhance laser specific energy is demonstrated, although preheating may impose unacceptable hardware requirements. The application of numerical optimization techniques appears to offer a means to improve laser design.

Saunders, R. C., III↗

Laser beam scintillation beyond the turbulent atmosphere A numerical computation

The extended Huygens-Fresnel formulation for propagation through turbulence is used to examine scintillation of a finite laser beam. The method is demonstrated analytically for propagation beyond a weak Gaussian phase screen. A numerical integration technique is used to extend the results to a more realistic turbulence model. Results are compared with existing Gaussian beam propagation theory.

Bufton, J. L.↗

Solar generated quasi-biennial geomagnetic variation

The existence of highly correlated quasi-biennial variations in the geomagnetic field and in solar activity is demonstrated. The analysis uses a numerical filter technique applied to monthly averages of the geomagnetic horizontal component and of the Zurich relative sunspot number. Striking correlations are found between the quasi-biennial geomagnetic variations determined from several magnetic observatories located at widely different longitudes, indicating a worldwide nature of the obtained variation. The correlation coefficient between the filtered Dst index and the filtered relative sunspot number is found to be -0.79 at confidence level greater than 99% with a time-lag of 4 months, with solar activity preceding the Dst variation. The correlation between the unfiltered data of Dst and of the sunspot number is also high with a similar time-lag. Such a timelag has not been discussed in the literature, and a further study is required to establish the mode of sun-earth relationship that gives this time delay.

Sugiura, M.↗

Numerical optimization of gasdynamic lasers

A sophisticated numerical optimization procedure has been coupled to a detailed theoretical model of the carbon dioxide gasdynamic laser (GDL) for the purpose of systematically evaluating the performance of these systems. The integrated model was used to predict the performance of a wide variety of combustion-driven gasdynamic lasers and was then used to maximize the theoretical available specific energy of these lasers. Up to six design variables and three nonlinear constraints were used in these optimizations, although the method is not limited to this number. Results are presented to show the performance of typical lasers before and after optimization. Improvements approaching a factor of two in the available specific energy are found, suggesting that numerical optimization techniques are valuable in improving the design of gasdynamic lasers.

Saunders, R. C., III↗