Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “residual minimization”

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 37 records · Page 2

On polynomial preconditioning for indefinite Hermitian matrices

The minimal residual method is studied combined with polynomial preconditioning for solving large linear systems (Ax = b) with indefinite Hermitian coefficient matrices (A). The standard approach for choosing the polynomial preconditioners leads to preconditioned systems which are positive definite. Here, a different strategy is studied which leaves the preconditioned coefficient matrix indefinite. More precisely, the polynomial preconditioner is designed to cluster the positive, resp. negative eigenvalues of A around 1, resp. around some negative constant. In particular, it is shown that such indefinite polynomial preconditioners can be obtained as the optimal solutions of a certain two parameter family of Chebyshev approximation problems. Some basic results are established for these approximation problems and a Remez type algorithm is sketched for their numerical solution. The problem of selecting the parameters such that the resulting indefinite polynomial preconditioners speeds up the convergence of minimal residual method optimally is also addressed. An approach is proposed based on the concept of asymptotic convergence factors. Finally, some numerical examples of indefinite polynomial preconditioners are given.

Freund, Roland W.↗

Electroform/Plasma-Spray Laminates for X-Ray Optics

Electroform/plasma-spray laminates have shown promise as lightweight, strong, low-thermal-expansion components for xray optics. The basic idea is to exploit both (1) the well-established art of fabrication of optical components by replication and (2) plasma spraying as a means of reinforcing a thin replica optic with one or more backing layer(s) having tailorable thermomechanical properties. In x-ray optics as in other applications, replication reduces the time and cost of fabrication because grinding and polishing can be limited to a few thick masters, from which many lightweight replicas can thereafter be made. The first step in the fabrication of a component of the type in question is to make a replica optic by electroforming a thin layer of nickel on a master. Through proper control of the electroforming process conditions, it is possible to minimize residual stress and, hence, to minimize distortion in the replica. Next, a powder comprising ceramic particles coated with a metal compatible with the electroformed nickel is plasma-sprayed onto the backside of the nickel replica. Then through several repetitions and variations of the preceding steps or perhaps a small compressive stress, alternating layers of electroformed nickel and plasma-sprayed metal-coated ceramic powder are deposited. The thicknesses of the layers and the composition of the metal-coated ceramic powder are chosen to optimize the strength, areal mass density, and toughness of the finished component. An important benefit of using both electroforming and plasma spraying is the possibility of balancing stresses to a minimum level, which could be zero or perhaps a small net compressive stress designed to enhance the function of the component in its intended application.

Ulmer, Melville P.↗

Optimal Output Trajectory Redesign for Invertible Systems

Given a desired output trajectory, inversion-based techniques find input-state trajectories required to exactly track the output. These inversion-based techniques have been successfully applied to the endpoint tracking control of multijoint flexible manipulators and to aircraft control. The specified output trajectory uniquely determines the required input and state trajectories that are found through inversion. These input-state trajectories exactly track the desired output; however, they might not meet acceptable performance requirements. For example, during slewing maneuvers of flexible structures, the structural deformations, which depend on the required state trajectories, may be unacceptably large. Further, the required inputs might cause actuator saturation during an exact tracking maneuver, for example, in the flight control of conventional takeoff and landing aircraft. In such situations, a compromise is desired between the tracking requirement and other goals such as reduction of internal vibrations and prevention of actuator saturation; the desired output trajectory needs to redesigned. Here, we pose the trajectory redesign problem as an optimization of a general quadratic cost function and solve it in the context of linear systems. The solution is obtained as an off-line prefilter of the desired output trajectory. An advantage of our technique is that the prefilter is independent of the particular trajectory. The prefilter can therefore be precomputed, which is a major advantage over other optimization approaches. Previous works have addressed the issue of preshaping inputs to minimize residual and in-maneuver vibrations for flexible structures; Since the command preshaping is computed off-line. Further minimization of optimal quadratic cost functions has also been previously use to preshape command inputs for disturbance rejection. All of these approaches are applicable when the inputs to the system are known a priori. Typically, outputs (not inputs) are specified in tracking problems, and hence the input trajectories have to be computed. The inputs to the system are however, difficult to determine for non-minimum phase systems like flexible structures. One approach to solve this problem is to (1) choose a tracking controller (the desired output trajectory is now an input to the closed-loop system and (2) redesign this input to the closed-loop system. Thus we effectively perform output redesign. These redesigns are however, dependent on the choice of the tracking controllers. Thus the controller optimization and trajectory redesign problems become coupled; this coupled optimization is still an open problem. In contrast, we decouple the trajectory redesign problem from the choice of feedback-based tracking controller. It is noted that our approach remains valid when a particular tracking controller is chosen. In addition, the formulation of our problem not only allows for the minimization of residual vibration as in available techniques but also allows for the optimal reduction fo vibrations during the maneuver, e.g., the altitude control of flexible spacecraft. We begin by formulating the optimal output trajectory redesign problem and then solve it in the context of general linear systems. This theory is then applied to an example flexible structure, and simulation results are provided.

Devasia, S.↗

Attitude-Independent Magnetometer Calibration for Spin-Stabilized Spacecraft

The paper describes a three-step estimator to calibrate a Three-Axis Magnetometer (TAM) using TAM and slit Sun or star sensor measurements. In the first step, the Calibration Utility forms a loss function from the residuals of the magnitude of the geomagnetic field. This loss function is minimized with respect to biases, scale factors, and nonorthogonality corrections. The second step minimizes residuals of the projection of the geomagnetic field onto the spin axis under the assumption that spacecraft nutation has been suppressed by a nutation damper. Minimization is done with respect to various directions of the body spin axis in the TAM frame. The direction of the spin axis in the inertial coordinate system required for the residual computation is assumed to be unchanged with time. It is either determined independently using other sensors or included in the estimation parameters. In both cases all estimation parameters can be found using simple analytical formulas derived in the paper. The last step is to minimize a third loss function formed by residuals of the dot product between the geomagnetic field and Sun or star vector with respect to the misalignment angle about the body spin axis. The method is illustrated by calibrating TAM for the Fast Auroral Snapshot Explorer (FAST) using in-flight TAM and Sun sensor data. The estimated parameters include magnetic biases, scale factors, and misalignment angles of the spin axis in the TAM frame. Estimation of the misalignment angle about the spin axis was inconclusive since (at least for the selected time interval) the Sun vector was about 15 degrees from the direction of the spin axis; as a result residuals of the dot product between the geomagnetic field and Sun vectors were to a large extent minimized as a by-product of the second step.

Natanson, Gregory↗

On conjugate gradient type methods and polynomial preconditioners for a class of complex non-Hermitian matrices

Conjugate gradient type methods are considered for the solution of large linear systems Ax = b with complex coefficient matrices of the type A = T + i(sigma)I where T is Hermitian and sigma, a real scalar. Three different conjugate gradient type approaches with iterates defined by a minimal residual property, a Galerkin type condition, and an Euclidian error minimization, respectively, are investigated. In particular, numerically stable implementations based on the ideas behind Paige and Saunder's SYMMLQ and MINRES for real symmetric matrices are proposed. Error bounds for all three methods are derived. It is shown how the special shift structure of A can be preserved by using polynomial preconditioning. Results on the optimal choice of the polynomial preconditioner are given. Also, some numerical experiments for matrices arising from finite difference approximations to the complex Helmholtz equation are reported.

Freund, Roland↗

Convergence acceleration of iterative solutions of Euler equations for transonic flow computations

In this paper, two acceleration techniques for Euler calculations are investigated. The first technique is an extrapolation procedure based on the Power Method; it is applicable when the iterative matrix has dominant eigenvalues. Both real and complex conjugate roots are allowed. The second technique is a generalization of the Minimal Residual Method, where the extrapolation step consists of a weighted combination of the corrections at different iteration levels and the weights are chosen to minimize the Lz norm of the residual. Numerical results, using Jameson's Runge-Kutta Multigrid Code, are presented. The extra computational work to apply either technique is negligible and the extra storage is not a problem on current supercomputers.

Hafez, M.↗

An extension of command shaping methods for controlling residual vibration using frequency sampling

The authors present an extension to the impulse shaping technique for commanding machines to move with reduced residual vibration. The extension, called frequency sampling, is a method for generating constraints that are used to obtain shaping sequences which minimize residual vibration in systems such as robots whose resonant frequencies change during motion. The authors present a review of impulse shaping methods, a development of the proposed extension, and a comparison of results of tests conducted on a simple model of the space shuttle robot arm. Frequency shaping provides a method for minimizing the impulse sequence duration required to give the desired insensitivity.

Singer, Neil C.↗

Stress In A Fiber During Curing Of Surrounding Matrix Resin

Experiments run to determine variation in tensile stress in single fiber during curing of matrix resin surrounding fiber. Study part of effort to understand physical mechanisms affecting residual stresses in matrix/fiber composites, with view toward optimizing curing cycles (in particular, optimizing temperature-vs.-time schedules of final cooldowns to ambient temperature) to minimize residual stresses. Results signify primary mechanisms affecting residual stress in fibers are thermal expansion and contraction and cure shrinkage of matrix material.

Bowles, Kenneth J.↗

The alignment-distribution graph

Implementing a data-parallel language such as Fortran 90 on a distributed-memory parallel computer requires distributing aggregate data objects (such as arrays) among the memory modules attached to the processors. The mapping of objects to the machine determines the amount of residual communication needed to bring operands of parallel operations into alignment with each other. We present a program representation called the alignment distribution graph that makes these communication requirements explicit. We describe the details of the representation, show how to model communication cost in this framework, and outline several algorithms for determining object mappings that approximately minimize residual communication.

Chatterjee, Siddhartha↗

The alignment-distribution graph

Implementing a data-parallel language such as Fortran 90 on a distributed-memory parallel computer requires distributing aggregate data objects (such as arrays) among the memory modules attached to the processors. The mapping of objects to the machine determines the amount of residual communication needed to bring operands of parallel operations into alignment with each other. We present a program representation called the alignment-distribution graph that makes these communication requirements explicit. We describe the details of the representation, show how to model communication cost in this framework, and outline several algorithms for determining object mappings that approximately minimize residual communication.

Chatterjee, Siddhartha↗

What is the minimal vestibular function required for compensation?

Living with an uncompensated, abnormal vestibular system requires oppressive modification of life style and often prevents return to work and activities of daily living. Patients with vestibular abnormalities were studied to determine the minimal residual vestibular function required to achieve compensation. Three groups of patients with (a) complete unilateral loss of vestibular function with normal horizontal canal-vestibulo-ocular (HCVOR) function in the opposite ear, (b) complete unilateral loss with abnormal HCVOR function in the opposite ear, and (c) bilateral reduction of vestibular function from aminoglycoside toxicity underwent vestibuloocular (VOR), optokinetic (OKN), visual-VOR (VVOR), and computerized dynamic posturography (CDP) tests before and after therapeutic procedures. Results suggest that a minimal VOR response amplitude must be present for compensation of VVOR function to occur. The roles of VOR and OKN phase shifts in vestibular compensation are more complicated and require further study. Compensation of vestibulospinal function does not necessarily accompany VOR or VVOR compensation. Ascending and descending vestibular compensatory mechanisms may involve different spatial sensory inputs. Results of these studies have important implications for the diagnosis and treatment of patients with vestibular disorders, including selection and monitoring of patients for therapeutic regimens such as vestibular nerve section and streptomycin therapy.

Non-NASA Center↗

Empirical orbit determination using Apollo 14 data

An empirical orbit determination method is shown to yield highly accurate navigation results when applied to lunar orbit tracking data. Regressions and predictions of free flight Apollo 14 tracking data exhibit minimal residual growth, and the solution orbital elements behave in a very consistent manner. Solutions from data acquired during propulsive maneuvers result in degraded predictions. The residual patterns from free flight processing are shown to be consistent from pass to pass and are correlated with lunar topographic features.

Ferrari, A. J.↗

Development of a gas pressure bonded four-pole alternator rotor

Methods were developed for fabrication of a solid four pole alternator rotor by hot isostatic pressure welding. The rotor blanks welded in this program had complex geometrical mating interfaces and were of considerable bulk, being approximately 3-1/2 inches (0.089 meters) in diameter and 14 inches (0.356 meters) long. Magnetic end pieces were machined from AlSl 4340 steel, while the non-magnetic central section was of Inconel 718. Excellent welds were produced which were shown to be responsive to post weld heat treatments which substantially improved joint strength. Prior to welding the rotors, test specimens of complex geometry were welded to demonstrate that complex surfaces with intentional mechanical misfit could be readily joined using HIP welding. This preliminary work demonstrated not only that interface compliance is achieved during welding but that welding pressure is developed in these thick sections sufficient to produce sound joints. Integral weld-heat treatment cycles were developed that permitted the attainment of magnetic properties while minimizing residual stress associated with the allotropic transformation of 4340 steel.

Lessmann, G. G.↗

An active thermal compensator for closed-cycle helium refrigerators

A technique was developed for reducing the amplitude of the temperature oscillation in He closed-cyle refrigerators. The device uses a semiconductor diode as a heating element to actively supply a small oscillating input of heat at a point between the laser and the cold-tip to cancel the heat oscillations due to the refrigerator. It was found that the heater diode could drive the temperature of the heat sink more effectively, i.e., with lower current and therefore less heat, if the heat sink was insulated slightly from the rest of the mount. A sine-wave generator was used to drive the programmable supply which provided the offset current to the heater diode. By matching the frequency and phase of the oscillator to that of the refrigerator cycle, and by adjusting the amplitude of the oscillator signal, the temperature fluctuations at the laser could be minimized. Residual fluctuations were about 0.003K peak-to-peak, at an operating temperature of 9.5K.

Jennings, D. E.↗

Three-Zone Catalyst Resists Sulfur Poisoning

Three-zone catalyst bed uses different types of nickel catalysts to convert sulfur-containing hydrocarbon fuels to hydrogen and carbon monoxide. Zones designed to achieve conversion with minimal residue of unconverted hydrocarbon, no soot and mimimal sulfur contamination.

Voecks, G. E.↗

Calculations of transonic potential flows by a parameter free procedure

A parameter-free computational procedure for transonic potential flows around an airfoil is presented. The method is based on an inexact, Newton-like iterative scheme, and a preconditioned minimal residual algorithm is used to solve the linear system of equations arising at each iterate. This iterative procedure is formulated so that the difficulties associated with the approximate factorization technique, namely the requirement of iteration parameters and the boundary conditions for the intermediate variable are eliminated. Numerical experiments are given for potential flows around NACA 0012 airfoil with different Mach numbers and different angles of attack. These results are then compared with those obtained by the approximate factorization techniques.

Wong, Y. S.↗

Some iterative schemes for transonic potential flows

The minimal residual (MR) method for the numerical solution of transonic potential flows is closely related to the conjugate gradient method, which has found widespread use in the solution of large sparse, symmetric, and positive-definite linear equations. The primary advantage of the MR method is its applicability to both symmetric and nonsymmetric matrices.

Wong, Y. S.↗

GMRES acceleration of computational fluid dynamics codes

The generalized minimal residual algorithm (GMRES) is a conjugate-gradient like method that applies directly to nonsymmetric linear systems of equations. In this paper, GMRES is modified to handle nonlinear equations characteristic of computational fluid dynamics. Attention is devoted to the concept of preconditioning and the role it plays in assuring rapid convergence. A formulation is developed that allows GMRES to be preconditioned by the solution procedures already built into existing computer codes. Examples are provided that demonstrate the ability of GMRES to greatly improve the robustness and rate of convergence of current state-of-the-art fluid dynamics codes. Theoretical aspects of GMRES are presented that explain why it works. Finally, the advantage GMRES enjoys over related methods such as conjugate gradients are discussed.

Wigton, L. B.↗