Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “method of multipliers”

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 145 records · Page 8

Modified Interior Distance Functions (Theory and Methods)

In this paper we introduced and developed the theory of Modified Interior Distance Functions (MIDF's). The MIDF is a Classical Lagrangian (CL) for a constrained optimization problem which is equivalent to the initial one and can be obtained from the latter by monotone transformation both the objective function and constraints. In contrast to the Interior Distance Functions (IDF's), which played a fundamental role in Interior Point Methods (IPM's), the MIDF's are defined on an extended feasible set and along with center, have two extra tools, which control the computational process: the barrier parameter and the vector of Lagrange multipliers. The extra tools allow to attach to the MEDF's very important properties of Augmented Lagrangeans. One can consider the MIDFs as Interior Augmented Lagrangeans. It makes MIDF's similar in spirit to Modified Barrier Functions (MBF's), although there is a fundamental difference between them both in theory and methods. Based on MIDF's theory, Modified Center Methods (MCM's) have been developed and analyzed. The MCM's find an unconstrained minimizer in primal space and update the Lagrange multipliers, while both the center and the barrier parameter can be fixed or updated at each step. The MCM's convergence was investigated, and their rate of convergence was estimated. The extension of the feasible set and the special role of the Lagrange multipliers allow to develop MCM's, which produce, in case of nondegenerate constrained optimization, a primal and dual sequences that converge to the primal-dual solutions with linear rate, even when both the center and the barrier parameter are fixed. Moreover, every Lagrange multipliers update shrinks the distance to the primal dual solution by a factor 0 less than gamma less than 1 which can be made as small as one wants by choosing a fixed interior point as a 'center' and a fixed but large enough barrier parameter. The numericai realization of MCM leads to the Newton MCM (NMCM). The approximation for the primal minimizer one finds by Newton Method followed by the Lagrange multipliers update. Due to the MCM convergence, when both the center and the barrier parameter are fixed, the condition of the MDF Hessism and the neighborhood of the primal ninimizer where Newton method is 'well' defined remains stable. It contributes to both the complexity and the numerical stability of the NMCM.

Polyak, Roman A.↗

Analytic Method for Computing Instrument Pointing Jitter

A new method of calculating the root-mean-square (rms) pointing jitter of a scientific instrument (e.g., a camera, radar antenna, or telescope) is introduced based on a state-space concept. In comparison with the prior method of calculating the rms pointing jitter, the present method involves significantly less computation. The rms pointing jitter of an instrument (the square root of the jitter variance shown in the figure) is an important physical quantity which impacts the design of the instrument, its actuators, controls, sensory components, and sensor- output-sampling circuitry. Using the Sirlin, San Martin, and Lucke definition of pointing jitter, the prior method of computing the rms pointing jitter involves a frequency-domain integral of a rational polynomial multiplied by a transcendental weighting function, necessitating the use of numerical-integration techniques. In practice, numerical integration complicates the problem of calculating the rms pointing error. In contrast, the state-space method provides exact analytic expressions that can be evaluated without numerical integration.

Bayard, David↗

Discrete fourier transform (DFT) analysis for applications using iterative transform methods

According to various embodiments, a method is provided for determining aberration data for an optical system. The method comprises collecting a data signal, and generating a pre-transformation algorithm. The data is pre-transformed by multiplying the data with the pre-transformation algorithm. A discrete Fourier transform of the pre-transformed data is performed in an iterative loop. The method further comprises back-transforming the data to generate aberration data.

Dean, Bruce H.↗

Extension of Case's method to the two-region critical cylinder and its application as an analytic standard

The normal mode expansion technique is applied to the transformed monoenergetic integral transport equation to develop a solution for the rotationally invariant and axially infinite, critical, two-region cylinder with a finite outer reflector boundary. Isotropic scattering and identical neutron mean free paths in the core and reflector regions are assumed. The solution in terms of singular integral equations is obtained by applying a completeness theorem found for the singular eigenfunctions. Numerical results for a variety of core and reflector multiplying properties and reflector thickness are presented and compared with the results of other method. An example of this type of application is given in a study of approximations inherent in the neutronic design analysis of a small, fast-neutron-spectrum reactor concept proposed as a space power source.

Westfall, R. M.↗

Aerodynamic design and optimization in one shot

This paper describes an efficient numerical approach for the design and optimization of aerodynamic bodies. As in classical optimal control methods, the present approach introduces a cost function and a costate variable (Lagrange multiplier) in order to achieve a minimum. High efficiency is achieved by using a multigrid technique to solve for all the unknowns simultaneously, but restricting work on a design variable only to grids on which their changes produce nonsmooth perturbations. Thus, the effort required to evaluate design variables that have nonlocal effects on the solution is confined to the coarse grids. However, if a variable has a nonsmooth local effect on the solution in some neighborhood, it is relaxed in that neighborhood on finer grids. The cost of solving the optimal control problem is shown to be approximately two to three times the cost of the equivalent analysis problem. Examples are presented to illustrate the application of the method to aerodynamic design and constraint optimization.

Ta'asan, Shlomo↗

Relaxation methods in fluid mechanics

The present work considers the iterative solution of a coupled set of difference equations and examines methods that carry successive approximates to a state that is invariant with further iteration and independent of the initial guess. Methods are studied with regard to their efficiency and economy of computer resources. The basic principles of classical relaxation are set forth, with attention confined to linear elliptic equations. This discussion involves the evaluation of the spectral radius that is the magnitude of the eigenvalue with largest modulus. The subject of relaxation is then related to the study of ordinary differential equations and hyperbolic partial differential equations. Problems that occur when linearly dependent eigenvectors appear in the relaxation matrix are discussed, leading to multiply connected eigenvalues in the Jordan canonical form. Finally, a brief survey of relaxation methods used in aerodynamics is given.

Lomax, H.↗

A flow study in radial inflow turbine scroll-nozzle assembly

The present analysis describes the flow behavior in the combined scroll-nozzle assembly of a radial inflow turbine. This model was chosen to provide a better understanding of the mutual interaction effects of these two components on the flow. The finite element method is used in the solution of the flow field in this multiply connected domain. The mass flow rates in the different nozzle channels is not presumed constant, but is determined from the solution.

Hamed, A.↗

Scattering of waves from periodic surfaces

In order to study the scattering of waves from periodic surfaces, the basic grating equations are reviewed, and a general approach for scattering from impenetrable and penetrable media is formulated. Three analytical methods for scattering on a conducting sinusoidal surface for both TE and TM polarized waves are compared. In The Masel, Merrill, and Miller (MMM) method for quantum scattering of atoms (1975, 1976), the surface field expansions are expressed in terms of Fourier series. The Modified Physical Optics (MPO) method (DeSanto, 1975; Whitman and Schwering, 1977) uses surface field expansions consisting of a leading term proportional to the physical optics approximation multiplied by a Fourier series expansion. The Waterman's Plane Harmonics (WPH) approach (1975) makes use of basis functions which are downward plane harmonics evaluated on the surface. The MMM method proved to be the most efficient one in terms of rate and range of convergence. For dielectric media with periodic rough surfaces, an improved method is developed for calculating the reflected and transmitted powers, and the results are compared with experimental data obtained at optical frequencies.

Chuang, S.-L.↗

A Variational Method for the Optimization of Interplanetary Round-Trip Trajectories

The indirect method of the calculus of variations is used to optimize interplanetary round-trip trajectories for the case of a single, central, attracting body. The method of solution makes use of certain partial derivative properties of the Lagrangian multipliers associated with the Mayer formulation of the variational problem. This property of the multipliers allows the construction of mathematical expressions for certain other partial derivatives that must vanish when an optimum round trip has been found. These expressions are developed for the cases of propulsion systems using (1) fixed thrust and specific impulse or (2) variable thrust and constant exhaust jet power. Two numerical examples demonstrate how the analytical results may be applied to the solution of round-trip problems including (1) actual three-dimensional planetary positions and (2) planetocentric maneuvers.

VARIATIONAL CALCULUS↗

Modified-signed-digit optical computing by using fan-out elements

We present a novel three-step modified-signed-digit optical computing system based on the disrete correlation of space-coded input matrices that are multiply-imaged by using optical fan-out elements. The advantages of our method over other approaches include a larger throughput, more flexible operation, and an easier optical implementation.

Zhou, Shaomin↗

Tropospheric Delay Signal in SAR Interferogram and its Correction for Precise Surface Change Detection

Synthetic Aperture Radar (SAR) interferometry has become an important tool for measuring the surface deformation and mapping topography. The largest error source of the SAR interferometry measurements is differential atmospheric delay of water vapor. It reflects detailed distribution of water vapor in troposphere at data acquisition. We found phase difference associated with atmospheric waves and severe local atmospheric phenomena in interferograms. To distinguish phase difference associated with surface deformation from tropospheric effect, we need several SAR interferograms including the time period of the deformation. Averaging the interferograms is an effective way to reduce the tropospheric delay from horizontal inhomogeneity of the water vapor distribution. Apart from the tropospheric delay of the horizontal water vapor inhomogeneity, we often find the differential phase correlated to the topography (elevation) in interferograms, which might cause error in interpretation of surface deformation. This phase is due to the differential tropospheric delay caused by the topography and vertical change of water vapor between two images in different atmospheric condition. Theoretical calculation shows that the phase difference can be approximated by linear expression of the elevation. We applied a simple and effective correction method that the error is removed by subtracting the DEM (Digital Elevation Model) multiplied a coefficient.

Fujiwara, Satoshi↗

Performance of NOAA-20 VIIRS Solar Diffuser and Solar Diffuser Stability Monitor

Visible Infrared Imaging Radiometer Suite (VIIRS) radiometrically calibrates its reflective solar bands (RSBs) primarily through a sunlit onboard solar diffuser (SD). The sunlit SD provides a known radiance under the condition that the absolute product of the SD screen transmittance and the bidirectional reflectance distribution function (BRDF) along the SD-to-telescope direction is accurately known. The BRDF changes due to solar exposure. The change, referred to as the H-factor, is monitored by the onboard SD stability monitor (SDSM). The accuracy of the retrieved H-factor propagates to the retrieved F-factor which corrects the scene spectral radiance. High accuracy of the retrieved H-factor relies on high accuracies in the SDSM screen relative effective transmittance and the relative product of the SD screen effective transmittance and the BRDF at the mission start, and a high SDSM detector signal-to-noise ratio (SNR). This article briefly reviews the algorithms used for the NOAA-20 (N20) VIIRS RSB on-orbit radiometric calibration. Additionally, we show the performance of the N20 VIIRS SDSM, giving the SDSM detector SNRs and the SDSM detector gain temporal changes. We develop a model for the SNRs. The model shows that the decreased SNRs in time are due to the detector gain decreases. We also show the N20 VIIRS SD on-orbit performance, measured by the retrieved H-factor and the estimated standard deviation of its error. The H-factor for the telescope view is obtained from the H-factor for the SDSM view, multiplied by an H-factor angular dependence term. We use an innovative method to determine the angular dependence, using the dependence obtained for the Suomi National Polar-orbiting Partnership (SNPP) VIIRS.

VIIRS↗

Pseudo-compressibility methods for the incompressible flow equations

Preconditioning methods to accelerate convergence to a steady state for the incompressible fluid dynamics equations are considered. The analysis relies on the inviscid equations. The preconditioning consists of a matrix multiplying the time derivatives. Thus the steady state of the preconditioned system is the same as the steady state of the original system. The method is compared to other types of pseudo-compressibility. For finite difference methods preconditioning can change and improve the steady state solutions. An application to viscous flow around a cascade with a non-periodic mesh is presented.

Turkel, Eli↗

Structural Optimization Using the Newton Modified Barrier Method

The Newton Modified Barrier Method (NMBM) is applied to structural optimization problems with large a number of design variables and constraints. This nonlinear mathematical programming algorithm was based on the Modified Barrier Function (MBF) theory and the Newton method for unconstrained optimization. The distinctive feature of the NMBM method is the rate of convergence that is due to the fact that the design remains in the Newton area after each Lagrange multiplier update. This convergence characteristic is illustrated by application to structural problems with a varying number of design variables and constraints. The results are compared with those obtained by optimality criteria (OC) methods and by the ASTROS program.

Khot, N. S.↗