Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Newton optimization”

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 109 records · Page 6

An historical survey of computational methods in optimal control.

Review of some of the salient theoretical developments in the specific area of optimal control algorithms. The first algorithms for optimal control were aimed at unconstrained problems and were derived by using first- and second-variation methods of the calculus of variations. These methods have subsequently been recognized as gradient, Newton-Raphson, or Gauss-Newton methods in function space. A much more recent addition to the arsenal of unconstrained optimal control algorithms are several variations of conjugate-gradient methods. At first, constrained optimal control problems could only be solved by exterior penalty function methods. Later algorithms specifically designed for constrained problems have appeared. Among these are methods for solving the unconstrained linear quadratic regulator problem, as well as certain constrained minimum-time and minimum-energy problems. Differential-dynamic programming was developed from dynamic programming considerations. The conditional-gradient method, the gradient-projection method, and a couple of feasible directions methods were obtained as extensions or adaptations of related algorithms for finite-dimensional problems. Finally, the so-called epsilon-methods combine the Ritz method with penalty function techniques.

Polak, E.↗

A Preconditioning Method for Shape Optimization Governed by the Euler Equations

We consider a classical aerodynamic shape optimization problem subject to the compressible Euler flow equations. The gradient of the cost functional with respect to the shape variables is derived with the adjoint method at the continuous level. The Hessian (second order derivative of the cost functional with respect to the shape variables) is approximated also at the continuous level, as first introduced by Arian and Ta'asan (1996). The approximation of the Hessian is used to approximate the Newton step which is essential to accelerate the numerical solution of the optimization problem. The design space is discretized in the maximum dimension, i.e., the location of each point on the intersection of the computational mesh with the airfoil is taken to be an independent design variable. We give numerical examples for 86 design variables in two different flow speeds and achieve an order of magnitude reduction in the cost functional at a computational effort of a full solution of the analysis partial differential equation (PDE).

Arian, Eyal↗

Tables Of Gaussian-Type Orbital Basis Functions

NASA technical memorandum contains tables of estimated Hartree-Fock wave functions for atoms lithium through neon and potassium through krypton. Sets contain optimized Gaussian-type orbital exponents and coefficients, and near Hartree-Fock quality. Orbital exponents optimized by minimizing restricted Hartree-Fock energy via scaled Newton-Raphson scheme in which Hessian evaluated numerically by use of analytically determined gradients.

Partridge, Harry↗

The Newton Modified Barrier Method for QP Problems

The Modified Barrier Functions (MBF) have elements of both Classical Lagrangians (CL) and Classical Barrier Functions (CBF). The MBF methods find an unconstrained minimizer of some smooth barrier function in primal space and then update the Lagrange multipliers, while the barrier parameter either remains fixed or can be updated at each step. The numerical realization of the MBF method leads to the Newton MBF method, where the primal minimizer is found by using Newton's method. This minimizer is then used to update the Lagrange multipliers. In this paper, we examine the Newton MBF method for the Quadratic Programming (QP) problem. It will be shown that under standard second-order optimality conditions, there is a ball around the primal solution and a cut cone in the dual space such that for a set of Lagrange multipliers in this cut cone, the method converges quadratically to the primal minimizer from any point in the aforementioned ball, and continues to do so after each Lagrange multiplier update. The Lagrange multipliers remain within the cut cone and converge linearly to their optimal values. Any point in this ball will be called a "hot start". Starting at such a "hot start", at most Omicron(1n 1n epsilon(exp -1)) Newton steps are sufficient to perform the primal minimization which is necessary for the Lagrange multiplier update. Here, epsilon > 0 is the desired accuracy. Because of the linear convergence of the Lagrange multipliers, this means that only Omicron(1n epsilon(exp -1))omicron(ln 1n epsilon(exp-1)) Newton steps are required to reach an epsilon-approximation to the solution from any "hot start". In order to reach the "hot start", one has to perform Omicron(square root(m) 1n C) Newton steps, where m characterizes the size of the problem and C > 0 is the condition number of the QP problem. This condition number will be characterized explicitly in terms of key parameters of the QP problem, which in turn depend on the input data and the size of the problem.

Melman, A.↗

The GPM Combined Algorithm

In this paper, the operational Global Precipitation Measurement (GPM) mission combined radar-radiometer algorithm is thoroughly described. The operational combined algorithm is designed to reduce uncertainties in GPM Core Observatory precipitation estimates by effectively integrating complementary information from the GPM Dual-Frequency Precipitation Radar (DPR) and the GPM Microwave Imager (GMI) into an optimal, physically consistent precipitation product. Although similar in many respects to previously developed combined algorithms, the GPM combined algorithm has several unique features that are specifically designed to meet the GPM objectives of deriving, based on GPM Core Observatory information, accurate and physically consistent precipitation estimates from multiple spaceborne instruments, and ancillary environmental data from reanalyses. The algorithm features an optimal estimation framework based on a statistical formulation of the Gauss-Newton method, a parameterization for the nonuniform distribution of precipitation within the radar fields of view, a methodology to detect and account for multiple scattering in Ka-band DPR observations, and a statistical deconvolution technique that allows for an efficient sequential incorporation of radiometer information into DPR precipitation retrievals.

Grecu, Mircea↗

Modeling Spectralon's Bidirectional Reflectance for In-Flight Calibration of Earth-Orbiting Sensors

This paper describes the results of a model inversion which was conducted using measured Spectralon BRDF data at several illumination angles. We can describe the bidirectional reflectance of the Spectralon and in paticular its slight anisotrophy with a physically-based model. The root mean square difference between the model and the observations is currently of the order of the noise in the data (+/-1%). It is expected that meaningful corrections to account for the anisotrophy of the Spectralon panels will be achieved in the next few years and will improve the calibration of the MISR instrument.

EOS↗

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

A decomposition algorithm is presented that 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 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.↗

Peak Seeking Control for Reduced Fuel Consumption with Preliminary Flight Test Results

The Environmentally Responsible Aviation project seeks to accomplish the simultaneous reduction of fuel burn, noise, and emissions. A project at NASA Dryden Flight Research Center is contributing to ERAs goals by exploring the practical application of real-time trim configuration optimization for enhanced performance and reduced fuel consumption. This peak-seeking control approach is based on Newton-Raphson algorithm using a time-varying Kalman filter to estimate the gradient of the performance function. In real-time operation, deflection of symmetric ailerons, trailing-edge flaps, and leading-edge flaps of a modified F-18 are directly optimized, and the horizontal stabilators and angle of attack are indirectly optimized. Preliminary results from three research flights are presented herein. The optimization system found a trim configuration that required approximately 3.5% less fuel flow than the baseline trim at the given flight condition. The algorithm consistently rediscovered the solution from several initial conditions. These preliminary results show the algorithm has good performance and is expected to show similar results at other flight conditions and aircraft configurations.

Brown, Nelson↗

Peak-Seeking Optimization of Spanwise Lift Distribution for Wings in Formation Flight

A method is presented for the optimization of the lift distribution across the wing of an aircraft in formation flight. The usual elliptical distribution is no longer optimal for the trailing wing in the formation due to the asymmetric nature of the encountered flow field. Control surfaces along the trailing edge of the wing can be configured to obtain a non-elliptical profile that is more optimal in terms of minimum drag. Due to the difficult-to-predict nature of formation flight aerodynamics, a Newton-Raphson peak-seeking controller is used to identify in real time the best aileron and flap deployment scheme for minimum total drag. Simulation results show that the peak-seeking controller correctly identifies an optimal trim configuration that provides additional drag savings above those achieved with conventional anti-symmetric aileron trim.

Hanson, Curtis E.↗

MASSIM, the Milli-Arc-Second Structure Imager

The MASSIM (Milli-Arc-Second Structure Imager) mission will use a set of achromatic diffractive-refractive Fresnel lenses to achieve imaging in the X-ray band with unprecedented angular resolution. It has been proposed for study within the context of NASA's "Astrophysics Strategic Mission Concept Studies" program. Lenses on an optics spacecraft will focus 5-11 keV X-rays onto detectors on a second spacecraft flying in formation 1000 km away. It will have a point-source sensitivity comparable with that of the current generation of major X-ray observatories (Chandra, XMM-Newton) but an angular resolution some three orders of magnitude better. MASSIM is optimized for the study of jets and other phenomena that occur in the immediate vicinity of black holes and neutron stars. It can also be used for studying other phenomena on the milli-arc-second scale, such as those involving proto-stars, the surfaces and surroundings of nearby active stars and interacting winds.

Skinner, Gerry↗

The Milli-Arc-Second Structure Imager, MASSIM: A New Concept for a High Angular Resolution X-ray Telescope

MASSIM, the Milli-Arc-Second Structure Imager, is a mission that has been proposed for study within the context of NASA's "Astrophysics Strategic Mission Concept Studies" program. It uses a set of achromatic diffractive-refractive Fresnel lenses on an optics spacecraft to focus 5-11 keV X-rays onto detectors on a second spacecraft flying in formation 1000 km away. It will have a point-source sensitivity comparable with that of the current generation of major X-ray observatories (Chandra, XMM-Newton) but an angular resolution some three orders of magnitude better. MASSIM is optimized for the study of jets and other phenomena that occur in the immediate vicinity of black holes and neutron stars. It can also be used for studying other astrophysical phenomena on the milli-arc-second scale, such as those involving proto-stars, the surfaces and surroundings of nearby active stars and interacting winds. After introducing the principle of diffractive imaging in the x-ray/gamma-ray regime, the MASSIM mission concept and baseline design will be described along with a discussion of the options and trade-offs within the X-ray optics design.

Skinner, Gerry↗

Optimization of flexible wing structures subject to strength and induced drag constraints

An optimization procedure for designing wing structures subject to stress, strain, and drag constraints is presented. The optimization method utilizes an extended penalty function formulation for converting the constrained problem into a series of unconstrained ones. Newton's method is used to solve the unconstrained problems. An iterative analysis procedure is used to obtain the displacements of the wing structure including the effects of load redistribution due to the flexibility of the structure. The induced drag is calculated from the lift distribution. Approximate expressions for the constraints used during major portions of the optimization process enhance the efficiency of the procedure. A typical fighter wing is used to demonstrate the procedure. Aluminum and composite material designs are obtained. The tradeoff between weight savings and drag reduction is investigated.

Haftka, R. T.↗

Peak-Seeking Optimization of Spanwise Lift Distribution for Wings in Formation Flight

A method is presented for the in-flight optimization of the lift distribution across the wing for minimum drag of an aircraft in formation flight. The usual elliptical distribution that is optimal for a given wing with a given span is no longer optimal for the trailing wing in a formation due to the asymmetric nature of the encountered flow field. Control surfaces along the trailing edge of the wing can be configured to obtain a non-elliptical profile that is more optimal in terms of minimum combined induced and profile drag. Due to the difficult-to-predict nature of formation flight aerodynamics, a Newton-Raphson peak-seeking controller is used to identify in real time the best aileron and flap deployment scheme for minimum total drag. Simulation results show that the peak-seeking controller correctly identifies an optimal trim configuration that provides additional drag savings above those achieved with conventional anti-symmetric aileron trim.

Hanson, Curtis E.↗

Optimization of self-acting herringbone journal bearing for maximum stability

Groove parameters were determined to maximize the stability of herringbone grooved journal bearings. Parameters optimized were groove depth, width, length, and angle. Optimization was performed using a small eccentricity, infinite groove analysis in conjunction with: (1) a previously developed Newton-Raphson procedure for bearings with the smooth member rotating or with the grooved member rotating at low compressibility numbers, and (2) a newly-developed vector technique for bearings with the grooved member rotating at high compressibility numbers. The design curves enable one to choose the optimum bearing for a wide range of operating conditions. These include: (1) compressibility numbers from 0 (incompressible) to 80, (2) length to diameter ratios from 1/4 to 2, and (3) smooth or grooved member rotating. Compared to bearings optimized to maximize load capacity, bearings optimized for stability: (1) allow a thousandfold increase in bearing-supported mass in some cases before onset of instability (the most dramatic increase are for bearings with small L/D operating at high compressibility numbers), and (2) lose no more than 77-percent of their load capacity in any case studied. Stability is much greater when the grooved member rotates.

Fleming, D. P.↗

Optimization of self-acting herringbone journal bearings for maximum stability

Groove parameters were determined to maximize the stability of herringbone grooved journal bearings. Parameters optimized were groove depth, width, length, and angle. Optimization was performed using a small eccentricity, infinite groove analysis in conjunction with: (1) a previously developed Newton-Raphson procedure for bearings with the smooth member rotating or with the grooved member rotating at low compressibility numbers; (2) a newly-developed vector technique for bearings with the grooved member rotating at high compressibility numbers. The design curves in this report enable one to choose the optimum bearing for a wide range of operating conditions. These include: (1) compressibility numbers from 0 to 80; (2) length to diameter ratios from 1/4 to 2; and (3) smooth or grooved member rotating. Compared to bearings optimized to maximize load capacity, bearings optimized for stability allow a thousandfold increase in bearing-supported mass in some cases before onset of instability. Stability is much greater when the grooved member rotates.

Fleming, D. P.↗

Optimization of self-acting herringbone-grooved journal bearings for maximum stability

Groove parameters were determined to maximize the stability of herringbone-grooved journal bearings. Parameters optimized were groove depth, width, length, and angle. Optimization was performed by using a small-eccentricity, infinite-groove analysis in conjunction with a previously developed Newton-Raphson procedure for bearings with the smooth member rotating or with the grooved member rotating at low compressibility numbers, and a newly developed vector technique for bearings with the grooved member rotating at high compressibility numbers. The design curves in this report enable one to choose the optimum bearing for a wide range of operating conditions. Compared with bearings optimized to maximize load capacity, bearings optimized for stability allow a thousandfold increase in bearing-supported mass in some cases before onset of instability, and lose no more than 77 percent of their load capacity in any case studied. Stability is much greater when the grooved member rotates.

Fleming, D. P.↗