SEARCH · Engineering Papers

Results for “numerical stability”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS1982

Numerical methods and computers used in elastohydrodynamic lubrication

Some of the methods of obtaining approximate numerical solutions to boundary value problems that arise in elastohydrodynamic lubrication are reviewed. The highlights of four general approaches (direct, inverse, quasi-inverse, and Newton-Raphson) are sketched. Advantages and disadvantages of these approaches are presented along with a flow chart showing some of the details of each. The basic question of numerical stability of the elastohydrodynamic lubrication solutions, especially in the pressure spike region, is considered. Computers used to solve this important class of lubrication problems are briefly described, with emphasis on supercomputers.

Hamrock, B. J.↗

NASA NTRS1984

An investigation of the transient thermal analysis of spur gears

A finite element computer program is developed for evaluating the transient behavior of surface temperature in high performance spur gears. The time dimension is implemented using two and three point finite difference schemes. The different schemes are provided for the purpose of numerical stability and convergence studies. A detailed explanation of the gear cooling process leading to the establishment of a modified Blok model is also included. Other conventional models for approximating the heat transfer coefficients are available for comparison. Preliminary results are given showing snap shots of gear temperature contours at the initial stages of tooth engagement.

El-Bayoumy, L. E.↗

NASA NTRS1983

A numerical procedure for predicting creep and delayed failures in laminated composites

A numerical procedure is described for predicting the viscoelastic response of general laminates. A nonlinear compliance model is used to predict the creep response of the individual laminae. A biaxial delayed failure model predicts ply failure. The numerical procedure, based on lamination theory, increases by increments through time to predict creep compliance and delayed failures in laminates. Numerical stability problems and experimental verification are discussed. Although the program has been quite successful in predicting creep of general laminates, the assumptions associated with lamination theory have resulted in erroneous bounds on the predicted material response. Delayed failure predictions have been conservative. Several improvements are suggested to increase the accuracy of the procedure.

Dillard, D. A.↗

NASA NTRS1986

Transonic airfoil analysis and design in nonuniform flow

A nonuniform transonic airfoil code is developed for applications in analysis, inverse design and direct optimization involving an airfoil immersed in propfan slipstream. Problems concerning the numerical stability, convergence, divergence and solution oscillations are discussed. The code is validated by comparing with some known results in incompressible flow. A parametric investigation indicates that the airfoil lift-drag ratio can be increased by decreasing the thickness ratio. A better performance can be achieved if the airfoil is located below the slipstream center. Airfoil characteristics designed by the inverse method and a direct optimization are compared. The airfoil designed with the method of direct optimization exhibits better characteristics and achieves a gain of 22 percent in lift-drag ratio with a reduction of 4 percent in thickness.

Chang, J. F.↗

NASA NTRS1986

Parametric effects of CFL number and artificial smoothing on numerical solutions using implicit approximate factorization algorithm

An implicit approximate factorization algorithm is employed to quantify the parametric effects of Courant number and artificial smoothing on numerical solutions of the unsteady 3-D Euler equations for a windmilling propeller (low speed) flow field. The results show that propeller global or performance chracteristics vary strongly with Courant number and artificial dissipation parameters, though the variation is such less severe at high Courant numbers. Candidate sets of Courant number and dissipation parameters could result in parameter-dependent solutions. Parameter-independent numerical solutions can be obtained if low values of the dissipation parameter-time step ratio are used in the computations. Furthermore, it is realized that too much artificial damping can degrade numerical stability. Finally, it is demonstrated that highly resolved meshes may, in some cases, delay convergence, thereby suggesting some optimum cell size for a given flow solution. It is suspected that improper boundary treatment may account for the cell size constraint.

Daso, E. O.↗

NASA NTRS1985

An investigation of the transient thermal analysis of spur gears

El-Bayoumy, L. E.↗

NASA NTRS1987

Lanczos modes for reduced-order control of flexible structures

Lanczos mode models represent low-frequency forced response better than do normal mode models and can be developed for both continuous and finite element structural representations. It was recommended that Lanczos mode models for systems with multiple input and/or rigid body modes should be developed; numerical stability of the Lanczos algorithm should be assessed; and control system designs employing the Lanczos mode models should be attempted.

Craig, Roy R., Jr.↗

NASA NTRS1986

Application of multiple objective optimization techniques to finite element model tuning

This report examines tuning a finite element model using vector optimization techniques. Structural models using finite element theory often need to be adjusted so they can accurately simulate the real structure. The goal is to tune the model such that it will reproduce data derived from structural tests. First, the performance indices are extremized using multiple objective optimization theory, producing a set of possible solutions. Next, the solutions are rank ordered according to a decision maker's preferences to select the best answer. The tuning process was applied to a T-38 horizontal stabilizer. Numerous weighted solutions contained a best static deformation model, a best frequency model and three intermediate combinations of these two models. This automated procedure proved to be a versatile method capable of producing solutions for many types of tuning problems.

Devore, Charles R.↗

NASA NTRS1986

Unsteady separation characteristics of airfoils operating under dynamic stall conditions

Unsteady viscous/inviscid interaction phenomena are investigated for airfoils operating under dynamic stall conditions, using coupling procedures between a time-dependent inviscid panel method and two-dimensional unsteady boundary layer codes. Two strategies are pursued: a coupling of the inviscid panel method with the boundary layer code, and a strong coupling of the inviscid panel method with the boundary layer code. Attention is given to the main features of the unsteady time-marching panel method and boundary layer codes, as well as to numerical stability and the phenomenon of unsteady separation.

Geissler, Wolfgang↗

NASA NTRS1987

The method of projected characteristics for the evolution of magnetic arches

A numerical method of solving fully nonlinear MHD equation is described. In particular, the formulation based on the newly developed method of projected characteristics (Nakagawa, 1981) suitable to study the evolution of magnetic arches due to motions of their foot-points is presented. The final formulation is given in the form of difference equations; therefore, the analysis of numerical stability is also presented. Further, the most important derivation of physically self-consistent, time-dependent boundary conditions (i.e. the evolving boundary equations) is given in detail, and some results obtained with such boundary equations are reported.

Nakagawa, Y.↗

NASA NTRS1987

Parallel pivoting combined with parallel reduction

Parallel algorithms for triangularization of large, sparse, and unsymmetric matrices are presented. The method combines the parallel reduction with a new parallel pivoting technique, control over generations of fill-ins and a check for numerical stability, all done in parallel with the work being distributed over the active processes. The parallel technique uses the compatibility relation between pivots to identify parallel pivot candidates and uses the Markowitz number of pivots to minimize fill-in. This technique is not a preordering of the sparse matrix and is applied dynamically as the decomposition proceeds.

Alaghband, Gita↗

NASA NTRS1988

Method to predict external store carriage characteristics at transonic speeds

Development of a computational method for prediction of external store carriage characteristics at transonic speeds is described. The geometric flexibility required for treatment of pylon-mounted stores is achieved by computing finite difference solutions on a five-level embedded grid arrangement. A completely automated grid generation procedure facilitates applications. Store modeling capability consists of bodies of revolution with multiple fore and aft fins. A body-conforming grid improves the accuracy of the computed store body flow field. A nonlinear relaxation scheme developed specifically for modified transonic small disturbance flow equations enhances the method's numerical stability and accuracy. As a result, treatment of lower aspect ratio, more highly swept and tapered wings is possible. A limited supersonic freestream capability is also provided. Pressure, load distribution, and force/moment correlations show good agreement with experimental data for several test cases. A detailed computer program description for the Transonic Store Carriage Loads Prediction (TSCLP) Code is included.

Rosen, Bruce S.↗

NASA NTRS1988

External store carriage loads prediction at transonic speeds

A computational method for prediction of external store carriage loads at transonic speeds is described. The geometric flexibility required for treatment of isolated and underwing, pylon-mounted stores is achieved by computing solutions on a five-level embedded grid arrangement. A completely automated grid generation procedure facilitates applications. Store modeling capability consists of bodies of revolution with multiple fore and aft fins. A body-conforming grid improves the accuracy of the computed store body flow field. A nonlinear finite difference relaxation scheme, developed specifically for modified transonic small disturbance flow equations, enhances numerical stability and accuracy. As a result, more accurate treatment of low aspect ratio, highly swept and tapered wing planforms is possible. A limited supersonic freestream capability is also provided. Pressure, load distribution, force and moment correlations show good agreement for several test cases.

Rosen, Bruce S.↗

NASA NTRS1988

Galilean satellite ephemeris improvement using Galileo tour encounter information

Accurate navigation of the satellite tour portion of the Galileo mission requires an accurate ephemeris of the Galilean satellites. The ephemeris is updated using radiometric and optical tracking data acquired during the satellite tour. The improved accuracy of the satellite ephemeris leads to improved targeting accuracy at subsequent encounters. The Galileo mission will benefit from improved targeting accuracy through reduced propellant costs and improved pointing accuracy. The predicted error in the updated ephemeris can be less than approximations inherent in the analytical theory used for the ephemeris, so an alternate numerical representation is applied. This alternate description shows promise but also raises questions of numerical stability.

Murrow, D. W.↗

NASA NTRS1989

Extensions and improvements on XTRAN3S

Improvements to the XTRAN3S computer program are summarized. Work on this code, for steady and unsteady aerodynamic and aeroelastic analysis in the transonic flow regime has concentrated on the following areas: (1) Maintenance of the XTRAN3S code, including correction of errors, enhancement of operational capability, and installation on the Cray X-MP system; (2) Extension of the vectorization concepts in XTRAN3S to include additional areas of the code for improved execution speed; (3) Modification of the XTRAN3S algorithm for improved numerical stability for swept, tapered wing cases and improved computational efficiency; and (4) Extension of the wing-only version of XTRAN3S to include pylon and nacelle or external store capability.

Borland, C. J.↗

NASA NTRS1989

Method to predict external store carriage loads at transonic speeds

A computational method for prediction of external store carriage loads at transonic speeds is described. The geometric flexibility required for treatment of isolated and underwing, pylon mounted stores is achieved by computing solutions on a five level embedded grid arrangement. A completely automated grid generation procedure facilitates applications. Store modeling capability consists of bodies of revolution with multiple fore and aft fins. A body conforming grid improves the accuracy of the computed store body flow field. A nonlinear finite difference relaxation scheme, developed specifically for modified transonic small disturbance flow equations, enhances numerical stability and accuracy. As a result, more accurate treatment of low aspect ratio, highly swept and tapered wing planforms is possible. A limited supersonic freestream capability is also provided. Pressure, load distribution, force and moment correlation show good agreement for several test cases.

Rosen, Bruce S.↗

NASA NTRS1989

Parallel pivoting combined with parallel reduction and fill-in control

Parallel algorithms for triangularization of large, sparse, and unsymmetric matrices are presented. The method combines the parallel reduction with a new parallel pivoting technique, control over generation of fill-ins and check for numerical stability, all done in parallel with the work being distributed over the active processes. The parallel pivoting technique uses the compatibility relation between pivots to identify parallel pivot candidates and uses the Markowitz number of pivots to minimize fill-in. This technique is not a preordering of the sparse matrix and is applied dynamically as the decomposition proceeds.

Alaghband, Gita↗

NASA NTRS1989

A generalized gradient algorithm for dynamic optimization

A gradient algorithm is developed that determines optimal trajectories with path equality constraints and terminal constraints. A generalized gradient is formed which improves both the performance index and the path equality constraints simultaneously. The algorithm is extended to treat terminal constraints by using Bryson's impulse response technique. The main features of this algorithm are its numerical stability and smooth convergence near the optimum.

Zhao, Yiyuan↗