Search NASA⌕ Search

SEARCH · Search NASA

Results for “Newton-Raphson”

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 91 records · Page 5

A two-level trajectory decomposition algorithm featuring optimal intermediate target selection

A decomposition algorithm is presented which optimizes complex missions by partitioning the trajectory into natural segments such as ascent or entry. Each segment defines a full-rank targeting subproblem. These are solved sequentially using the Newton-Raphson algorithm. The master problem, representing the complete mission, is to determine subproblem targets and master-problem controls that optimize the mission objective subject to intersegment constraints. The gradient projection algorithm solves this problem using derivatives obtained analytically from finite-difference subproblem sensitivities. Thus, the mission is optimized by coordinating the solution of tractible subproblems. Computational results for a synchronous equatorial mission are included.

Petersen, F. M.↗

A finite element thermal analysis procedure for several temperature-dependent parameters

A finite-element thermal analysis procedure for elements with several temperature-dependent thermal parameters is presented. The procedure, based on an application of the Newton-Raphson iteration technique, is formulated by resolving element matrices into component matrices, one component for each thermal parameter. Component conductance matrices are evaluated by assuming constant thermal parameters within an element and are computed once per unit thermal parameter. Significant savings in computer time result from the unit thermal parameter concept. The solution procedure applied to a convectively cooled structure with significantly varying thermal parameters converged in four iterations.

Thornton, E. A.↗

TAP 2: A finite element program for thermal analysis of convectively cooled structures

A finite element computer program (TAP 2) for steady-state and transient thermal analyses of convectively cooled structures is presented. The program has a finite element library of six elements: two conduction/convection elements to model heat transfer in a solid, two convection elements to model heat transfer in a fluid, and two integrated conduction/convection elements to represent combined heat transfer in tubular and plate/fin fluid passages. Nonlinear thermal analysis due to temperature-dependent thermal parameters is performed using the Newton-Raphson iteration method. Transient analyses are performed using an implicit Crank-Nicolson time integration scheme with consistent or lumped capacitance matrices as an option. Program output includes nodal temperatures and element heat fluxes. Pressure drops in fluid passages may be computed as an option. User instructions and sample problems are presented in appendixes.

Thornton, E. A.↗

Testing of turbulent seals for rotodynamic coefficients

A test program developed for dynamic testing of straight and convergent-tapered seals, with the capability of separately determining both direct and cross-coupled stiffness, damping, and added mass coefficients is described. The test apparatus causes the seal journal to execute small-eccentricity centered circular orbits within its bearings. Dynamic measurements are made and recorded of the seal-displacement-vector components, and of the pressure field. The pressure field is integrated to yield seal reaction force components. The displacement and force vector components are analyzed via a generalized Newton-Raphson procedure to yield the desired seal dynamic coefficients. Representative test data are provided and discussed.

Childs, D. W.↗

A general algorithm for the solution of Kepler's equation for elliptic orbits

An efficient algorithm is presented for the solution of Kepler's equation f(E)=E-M-e sin E=0, where e is the eccentricity, M the mean anomaly and E the eccentric anomaly. This algorithm is based on simple initial approximations that are cubics in M, and an iterative scheme that is a slight generalization of the Newton-Raphson method. Extensive testing of this algorithm has been performed on the UNIVAC 1108 computer. Solutions for 20,000 pairs of values of e and M show that for single precision, 42.0% of the cases require one iteration, 57.8% two and 0.2% three. For double precision one additional iteration is required.

Ng, E. W.↗

NLSCIDNT user's guide maximum likehood parameter identification computer program with nonlinear rotorcraft model

A nonlinear, maximum likelihood, parameter identification computer program (NLSCIDNT) is described which evaluates rotorcraft stability and control coefficients from flight test data. The optimal estimates of the parameters (stability and control coefficients) are determined (identified) by minimizing the negative log likelihood cost function. The minimization technique is the Levenberg-Marquardt method, which behaves like the steepest descent method when it is far from the minimum and behaves like the modified Newton-Raphson method when it is nearer the minimum. Twenty-one states and 40 measurement variables are modeled, and any subset may be selected. States which are not integrated may be fixed at an input value, or time history data may be substituted for the state in the equations of motion. Any aerodynamic coefficient may be expressed as a nonlinear polynomial function of selected 'expansion variables'.

Source record↗

Recent developments in quasi-Newton methods for structural analysis and synthesis

Unlike the Newton-Raphson method, quasi-Newton methods by virture of the updates and step length control procedures are globally convergent and hence better suited for the solution of nonlinear problems of structural analysis and synthesis. Extension of quasi-Newton algorithms to large scale problems has led to the development of sparse update algorithms and to economical strategies for evaluating sparse Hessians. Ill-conditioning problems have led to the development of self-scaled variable metric and conjugate gradient algorithms, as well as the use of the singular perturbation theory. This paper emphasizes the effectiveness of such quasi-Newton algorithms for nonlinear structural analysis and synthesis.

Kamat, M. P.↗

HADY-I, a FORTRAN program for the compressible stability analysis of three-dimensional boundary layers

A computer program HADY-I for calculating the linear incompressible or compressible stability characteristics of the laminar boundary layer on swept and tapered wings is described. The eigenvalue problem and its adjoint arising from the linearized disturbance equations with the appropriate boundary conditions are solved numerically using a combination of Newton-Raphson interative scheme and a variable step size integrator based on the Runge-Kutta-Fehlburh fifth-order formulas. The integrator is used in conjunction with a modified Gram-Schmidt orthonormalization procedure. The computer program HADY-I calculates the growth rates of crossflow or streamwise Tollmien-Schlichting instabilities. It also calculates the group velocities of these disturbances. It is restricted to parallel stability calculations, where the boundary layer (meanflow) is assumed to be parallel. The meanflow solution is an input to the program.

El-Hady, N. M.↗

The SPAR thermal analyzer: Present and future

The SPAR thermal analyzer, a system of finite-element processors for performing steady-state and transient thermal analyses, is described. The processors communicate with each other through the SPAR random access data base. As each processor is executed, all pertinent source data is extracted from the data base and results are stored in the data base. Steady state temperature distributions are determined by a direct solution method for linear problems and a modified Newton-Raphson method for nonlinear problems. An explicit and several implicit methods are available for the solution of transient heat transfer problems. Finite element plotting capability is available for model checkout and verification.

Marlowe, M. B.↗

On nonlinear finite element analysis in single-, multi- and parallel-processors

Numerical solution of nonlinear equilibrium problems of structures by means of Newton-Raphson type iterations is reviewed. Each step of the iteration is shown to correspond to the solution of a linear problem, therefore the feasibility of the finite element method for nonlinear analysis is established. Organization and flow of data for various types of digital computers, such as single-processor/single-level memory, single-processor/two-level-memory, vector-processor/two-level-memory, and parallel-processors, with and without sub-structuring (i.e. partitioning) are given. The effect of the relative costs of computation, memory and data transfer on substructuring is shown. The idea of assigning comparable size substructures to parallel processors is exploited. Under Cholesky type factorization schemes, the efficiency of parallel processing is shown to decrease due to the occasional shared data, just as that due to the shared facilities.

Utku, S.↗

Self-adaptive closed constrained solution algorithms for nonlinear conduction

Self-adaptive solution algorithms are developed for nonlinear heat conduction problems encountered in analyzing materials for use in high temperature or cryogenic conditions. The nonlinear effects are noted to occur due to convection and radiation effects, as well as temperature-dependent properties of the materials. Incremental successive substitution (ISS) and Newton-Raphson (NR) procedures are treated as extrapolation schemes which have solution projections bounded by a hyperline with an externally applied thermal load vector arising from internal heat generation and boundary conditions. Closed constraints are formulated which improve the efficiency and stability of the procedures by employing closed ellipsoidal surfaces to control the size of successive iterations. Governing equations are defined for nonlinear finite element models, and comparisons are made of results using the the new method and the ISS and NR schemes for epoxy, PVC, and CuGe.

Padovan, J.↗

Three-dimensional flow measurements in a vaneless radial turbine scroll

The flow behavior in a vaneless radial turbine scroll was examined experimentally. The data was obtained using the slant sensor technique of hot film anemometry. This method used the unsymmetric heat transfer characteristics of a constant temperature hot film sensor to detect the flow direction and magnitude. This was achieved by obtaining a velocity vector measurement at three sensor positions with respect to the flow. The true magnitude and direction of the velocity vector was then found using these values and a Newton-Raphson numerical technique. The through flow and secondary flow velocity components are measured at various points in three scroll sections.

Tabakoff, W.↗

Efficient computerized model for dynamic analysis of energy conversion systems

In searching for the optimum parameters that minimize the total life cycle cost of an energy conversion system, various combinations of components are examined and the resulting system performance and associated economics are studied. The systems performance and economics simulation computer program (SPECS) was developed to fill this need. The program simulates the fluid flow, thermal, and electrical characteristics of a system of components on a quasi-steady state basis for a variety of energy conversion systems. A unique approach is used in which the set of characteristic equations is solved by the Newton-Raphson technique. This approach eliminates the tedious iterative loops which are found in comparable programs such as TRNSYS or SOLTES-1. Several efficient features were also incorporated such as the centralized control and energy management scheme, and analogous treatment of energy flow in electrical and mechanical components, and the modeling of components of similar fundamental characteristics using generic subroutines. Initial tests indicate that this model can be used effectively with a relatively small number of time steps and low computer cost.

Hughes, R. D.↗

On the solution of elastic-plastic static and dynamic postbuckling collapse of general structure

Many investigations have considered structural collapse from strictly the transient point of view. While such an approach is ideally correct, certain difficulties have to be overcome in its implementation. The present investigation is concerned with the development of self-adaptive algorithms which make it possible to conduct the analysis of both static elastic and elastic-plastic postbuckling, as well as static loading to the onset of buckling followed by subsequent dynamic postbuckling. The approach employed to solve the static portion of loading is to extend the constrained Incremental Newton-Raphson (INR) algorithm by incorporating elastic-plastic constitutive characterizations. Large deformation moderate strain theory is adopted to establish the overall strategy. Attention is given to governing field equations, aspects of algorithmic development, and numerical experiments conducted to illustrate the efficiency and stability of the developed schemes.

Padovan, J.↗

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

Longitudinal stability and control derivatives obtained from flight data of a PA-30 aircraft

In order to obtain reliable and accurate values of the stability and control derivatives, the Dryden Fligh Research Center (DFRC) developed a technique for extracting the derivatives from flight data. This technique is implemented by a set of FORTRAN computer programs that is based on a modified maximum likelihood estimator that uses the Newton-Raphson algorithm to perform the required minimization of the derivatives. Data was obtained with a PA-30, light twin-engine general aviation aircraft in zero, half, and full flap configuration in level unaccelerated flight with the landing gear retracted. The derivatives were plotted as functions of angle of attack using various graphical arrangements to show variations of wind tunnel and flight determined values at zero flap settings. Also, data was displayed to show the effects of flap deflection and thrust variation on the longitudinal stability derivatives. The angle of attack and angle of sideslip were measured. The dynamic pressure, velocity, and altitude were calculated, using a FORTRAN computer program, from the static and dynamic pressures. The control deflections of the stabilator, ailerons, and rudder also were recorded along with left throttle position, engine rpm, and manifold pressure.

Turley, D. R.↗