Search NASA⌕ Search

SEARCH · Search NASA

Results for “STEEPEST DESCENT METHOD”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

A theory for the radiation of magnetohydrodynamic surface waves and body waves into the solar corona

The Green's function for the slab coronal hole is obtained explicitly. The Fourier integral representation for the radiated field inside and outside the coronal hole waveguide is obtained. The radiated field outside the coronal hole is calculated using the method of steepest descents. It is shown that the radiated field can be written as the sum of two contributions: (1) a contribution from the integral along the steepest descent path and (2) a contribution from all the poles of the integrand between the path of the original integral and the steepest descent path. The free oscillations of the waveguide can be associated with the pole contributions in the steepest descent representation for the Green's function. These pole contributions are essentially generalized surface waves with a maximum amplitude near the interface which separates the plasma inside the coronal hole from the surrounding background corona. The path contribution to the integral is essentially the power radiated in body waves.

Davila, Joseph M.↗
NASA NTRS1990

Asymptotic boundary conditions for dissipative waves: General theory

An outstanding issue in the computational analysis of time dependent problems is the imposition of appropriate radiation boundary conditions at artificial boundaries. Accurate conditions are developed which are based on the asymptotic analysis of wave propagation over long ranges. Employing the method of steepest descents, dominant wave groups are identified and simple approximations to the dispersion relation are considered in order to derive local boundary operators. The existence of a small number of dominant wave groups may be expected for systems with dissipation. Estimates of the error as a function of domain size are derived under general hypotheses, leading to convergence results. Some practical aspects of the numerical construction of the asymptotic boundary operators are also discussed.

Hagstrom, Thomas↗
NASA NTRS1991

Asymptotic boundary conditions for dissipative waves - General theory

An outstanding issue in computational analysis of time dependent problems is the imposition of appropriate radiation boundary conditions at artificial boundaries. Accurate conditions are developed which are based on the asymptotic analysis of wave propagation over long ranges. Employing the method of steepest descents, dominant wave groups are identified and simple approximations to the dispersion relation are considered in order to derive local boundary operators. The existence of a small number of dominant wave groups may be expected for systems with dissipation. Estimates of the error as a function of domain size are derived under general hypotheses, leading to convergence results. Some practical aspects of the numerical construction of the asymptotic boundary operators are also discussed.

Hagstrom, Thomas↗
NASA NTRS1994

A linear method for analyzing lightning field changes

A constrained, least-squares method for analyzing multiple-station measurements of lightning field changes (delta Es) is introduced. Previous methods have attempted to fit the spatial pattern of lightning delta Es using nonlinear models, such as a point charge (Q) or a point dipole (P) model. With the linear method, the delta Es are described not by models but by a general volume charge distribution that is deposited on a large (40 x 40 x 20 cu km) Cartesian grid above the measuring network. A linear system of equations is used to relate the measured delta Es to the charges that are deposited at each grid point. With this approach, the information content of the measurements can be quantified by an eigenanalysis of the covariance matrix of the linear system. Constraints can be used to reduce the infinity of possible solutions to the linear system and also to reduce systematic biases that can be introduced by the method of solution. It is shown that a Landweber iterative method, derived from the general method of steepest descent, can be used to solve the linear system and that the resulting volume charge distributions are generally consistent with computer-simulated charge sources, when these sources are over the measuring network. The Landweber iteration has also provided solutions for natural lightning events that are consistent with Q- and P-model results.

Koshak, William J.↗
NASA NTRS2000

The nu Andromedae System: Models and Stability

Radial velocity observations of the F8 V star nu Andromedae taken at Lick and at Whipple Observatories have revealed evidence of three periodicities in the line-of-sight velocity of the star. These periodicities have been interpreted as evidence for at least three low-mass companions (LMCs) revolving around nu Andromedae. The mass and orbital parameters inferred for these companions raise questions about the dynamical stability of the system. We report here results from our independent analysis of the published radial velocity data, as well as new unpublished data taken at Lick Observatory. Our results confirm the finding of three periods in the data. Our best fits to the data, on the assumption that these periods arise from the gravitational perturbations of companions in Keplerian orbits, are also generally in agreement but with some differences from the earlier findings. We find that the available data do not constrain well the orbital eccentricity of the middle companion in a three-companion model of the data. We also find that in order for our best-fit model to the Lick data to be dynamically stable over the lifetime of the star (approximately 2 billion years), the system must have a mean inclination to the plane of the sky greater than 13 deg. The corresponding minimum inclination for the best fit to the Whipple data set is 19 deg. These values imply that the maximum mass for the outer companion can be no greater than about 20 Jupiter masses. Our analysis of the stability of the putative systems also places constraints on the relative inclinations of the orbital planes of the companions. We comment on global versus local (i.e., method of steepest descent) means of finding best-fit orbits from radial velocity data sets.

Stepinski, Tomasz F.↗
NASA NTRS2010

Pitch Guidance Optimization for the Orion Abort Flight Tests

The National Aeronautics and Space Administration created the Constellation program to develop the next generation of manned space vehicles and launch vehicles. The Orion abort system is initiated in the event of an unsafe condition during launch. The system has a controller gains schedule that can be tuned to reduce the attitude errors between the simulated Orion abort trajectories and the guidance trajectory. A program was created that uses the method of steepest descent to tune the pitch gains schedule by an automated procedure. The gains schedule optimization was applied to three potential abort scenarios; each scenario tested using the optimized gains schedule resulted in reduced attitude errors when compared to the Orion production gains schedule.

Stillwater, Ryan Allanque↗
NASA NTRS1975

Application of multivariable search techniques to the optimization of airfoils in a low speed nonlinear inviscid flow field

Multivariable search techniques are applied to a particular class of airfoil optimization problems. These are the maximization of lift and the minimization of disturbance pressure magnitude in an inviscid nonlinear flow field. A variety of multivariable search techniques contained in an existing nonlinear optimization code, AESOP, are applied to this design problem. These techniques include elementary single parameter perturbation methods, organized search such as steepest-descent, quadratic, and Davidon methods, randomized procedures, and a generalized search acceleration technique. Airfoil design variables are seven in number and define perturbations to the profile of an existing NACA airfoil. The relative efficiency of the techniques are compared. It is shown that elementary one parameter at a time and random techniques compare favorably with organized searches in the class of problems considered. It is also shown that significant reductions in disturbance pressure magnitude can be made while retaining reasonable lift coefficient values at low free stream Mach numbers.

Hague, D. S.↗
NASA NTRS1973

Coordinated adaptive washout for motion simulators.

This paper introduces a new method of providing motion cues to a moving base six-degree-of-freedom flight simulator utilizing nonlinear filters. Coordinated adaptive filters, used to coordinate translational and rotational motion, are derived based on the method of continuous steepest descent, and the basic concept of the digital controllers used for the uncoordinated heave and yaw cues is also presented. The coordinated adaptive washout method is illustrated by an application in a six-degree-of-freedom fixed-base environment.

Parrish, R. V.↗
NASA NTRS2003

Dissociative Ionization and Product Distributions of Benzene and Pyridine by Electron Impact

We report a theoretical study of the dissociative ionization (DI) and product distributions of benzene (C6H6) and pyridine (C5H5N) from their low-lying ionization channels. Our approach makes use of the fact that electronic motion is much faster than nuclear motion allowing DI to be treated as a two-step process. The first step is the electron-impact ionization resulting in an ion with the same nuclear geometry as the neutral molecule. In the second step, the nuclei relax from the initial geometry and undergo unimolecular dissociation. For the ionization process we use the improved binary-encounter dipole (iBED) model [W.M. Huo, Phys. Rev. A64,042719-I (2001)]. For the unimolecular dissociation, we use multiconfigurational self-consistent field (MCSCF) methods to determine the steepest descent pathways to the possible product channels. More accurate methods are then used to obtain better energetics of the paths which are used to determine unimolecular dissociation probabilities and product distributions. Our analysis of the dissociation products and the thresholds of their productions for benzene are compared with the recent dissociative photoionization meausurements of benzene by Feng et al. [R. Feng, G. Cooper, C.E. Brion, J. Electron Spectrosc. Relat. Phenom. 123,211 (2002)] and the dissociative photoionization measurements of pyridine by Tixier et al. [S. Tixier, G. Cooper, R. Feng, C.E. Brion, J. Electron Spectrosc. Relat. Phenom. 123,185 (2002)] using dipole (e,e+ion) coincidence spectroscopy.

Dateo, Christopher E.↗
NASA NTRS1976

Visual/motion simulation of CTOL flare and touchdown comparing data obtained from two model board display systems

Acquisition of a modern terrain model board display system for 'out-the-window' scene presentation has allowed for the repetition of an evaluation study of combined visual/motion cues for CTOL flare and touchdown control that was originally carried out on a 1965 vintage landing display system. The motion drive system used in both studies was a nonlinear coordinated adaptive washout based on the method of continuous steepest descent optimization. Comparisons of the results of the studies present flare and touchdown data, as well as dynamic response data, from the two visual landing systems. The new visual system also allowed for the repetition, under visual conditions, of a subjective comparison of a linear and a nonlinear motion washout method done previously under instrument conditions.

Parrish, R. V.↗
NASA NTRS1991

New displacement-based methods for optimal truss topology design

Two alternate methods for maximum stiffness truss topology design are presented. The ground structure approach is used, and the problem is formulated in terms of displacements and bar areas. This large, nonconvex optimization problem can be solved by a simultaneous analysis and design approach. Alternatively, an equivalent, unconstrained, and convex problem in the displacements only can be formulated, and this problem can be solved by a nonsmooth, steepest descent algorithm. In both methods, the explicit solving of the equilibrium equations and the assembly of the global stiffness matrix are circumvented. A large number of examples have been studied, showing the attractive features of topology design as well as exposing interesting features of optimal topologies.

Bendsoe, Martin P.↗
NASA NTRS1976

Solution of the minimum time-to-climb problem by matched asymptotic expansions

Application of singular perturbation techniques to trajectory optimization problems of flight mechanics is discussed. The method of matched asymptotic expansions is used to obtain an approximate solution to the aircraft minimum time-to-climb problem. Outer, boundary-layer, and composite solutions are obtained to zeroth and first orders. A stability criterion is derived for the zeroth-order boundary-layer solutions (the theory requires a form of boundary-layer stability). A numerical example is considered for which it is shown that the stability criterion is satisfied and a useful numerical solution is obtained. The zeroth-order solution proves to be a poor approximation, but the first-order solution gives a good approximation for both the trajectory and the minimum time-to-climb. The computational cost of the singular perturbation solution is considerably less than that of a steepest descent solution. Thus singular perturbation methods appear to be promising for the solution of optimal control problems.

Ardema, M. D.↗
NASA NTRS1976

Gradient optimization and nonlinear control

The book represents an introduction to computation in control by an iterative, gradient, numerical method, where linearity is not assumed. The general language and approach used are those of elementary functional analysis. The particular gradient method that is emphasized and used is conjugate gradient descent, a well known method exhibiting quadratic convergence while requiring very little more computation than simple steepest descent. Constraints are not dealt with directly, but rather the approach is to introduce them as penalty terms in the criterion. General conjugate gradient descent methods are developed and applied to problems in control.

Hasdorff, L.↗
NASA NTRS1997

A Nonlinear Programming Perspective on Sensitivity Calculations for Systems Governed by State Equations

This paper discusses the calculation of sensitivities. or derivatives, for optimization problems involving systems governed by differential equations and other state relations. The subject is examined from the point of view of nonlinear programming, beginning with the analytical structure of the first and second derivatives associated with such problems and the relation of these derivatives to implicit differentiation and equality constrained optimization. We also outline an error analysis of the analytical formulae and compare the results with similar results for finite-difference estimates of derivatives. We then attend to an investigation of the nature of the adjoint method and the adjoint equations and their relation to directions of steepest descent. We illustrate the points discussed with an optimization problem in which the variables are the coefficients in a differential operator.

Lewis, Robert Michael↗
NASA NTRS2018

Modeling the Proton Radiation Belt With Van Allen Probes Relativistic Electron-Proton Telescope Data

An empirical model of the proton radiation belt is constructed from data taken during 2013-2017 by the Relativistic Electron-Proton Telescopes on the Van Allen Probes satellites. The model intensity is a function of time, kinetic energy in the range 18-600 megaelectronvolts, equatorial pitch angle, and L shell of proton guiding centers. Data are selected, on the basis of energy deposits in each of the nine silicon detectors, to reduce background caused by hard proton energy spectra at low L. Instrument response functions are computed by Monte Carlo integration, using simulated proton paths through a simplified structural model, to account for energy loss in shielding material for protons outside the nominal field of view. Overlap of energy channels, their wide angular response, and changing satellite orientation require the model dependencies on all three independent variables be determined simultaneously. This is done by least squares minimization with a customized steepest descent algorithm. Model uncertainty accounts for statistical data error and systematic error in the simulated instrument response. A proton energy spectrum is also computed from data taken during the 8 January 2014 solar event, to illustrate methods for the simpler case of an isotropic and homogeneous model distribution. Radiation belt and solar proton results are compared to intensities computed with a simplified, on-axis response that can provide a good approximation under limited circumstances.

Selesnick, R. S.↗
NASA NTRS1973

High speed flow past wings

The analytical solution to the transonic small perturbation equation which describes steady compressible flow past finite wings at subsonic speeds can be expressed as a nonlinear integral equation with the perturbation velocity potential as the unknown function. This known formulation is substituted by a system of nonlinear algebraic equations to which various methods are applicable for its solution. Due to the presence of mathematical discontinuities in the flow solutions, however, a main computational difficulty was to ensure uniqueness of the solutions when local velocities on the wing exceeded the speed of sound. For continuous solutions this was achieved by embedding the algebraic system in an one-parameter operator homotopy in order to apply the method of parametric differentiation. The solution to the initial system of equations appears then as a solution to a Cauchy problem where the initial condition is related to the accompanying incompressible flow solution. In using this technique, however, a continuous dependence of the solution development on the initial data is lost when the solution reaches the minimum bifurcation point. A steepest descent iteration technique was therefore, added to the computational scheme for the calculation of discontinuous flow solutions. Results for purely subsonic flows and supersonic flows with and without compression shocks are given and compared with other available theoretical solutions.

Norstrud, H.↗
NASA NTRS2004

Low-thrust Orbit Transfers Using Candidate Lyapunov Functions with a Mechanism for Coasting

We consider low-thrust orbit transfers around a central body, where specified changes are sought in orbit elements except true anomaly. The desired changes in the remaining five elements can be arbitrarily large. Candidate Lyapunov functions are created based on analytic expressions for maximum rates of change of the orbit elements and the desired changes in the elements. These functions may be thought of as proximity quotients because they provide R measure of the proximity to the target orbit. The direction of thrust needed for steepest descent to the target orbit is also available analytically. The thrust is shutoff if the effectivity of the thrust at the current location on the osculating orbit is below some threshhold value. Thus, the equations of motion can be numerically integrated to obtain quickly and simply a transfer to the target orbit. A series of transfers can be easily computed to assess the trade-off between propellant mass and flight time. Preliminary comparisons to optimal solutions show that the method, while sub-optimal, performs well.

Candidate Lyapunoc functions↗