Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inverse problem”

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 901 records · Page 50

Large-scale kinetic energy spectra from Eulerian analysis of EOLE wind data

A data set of 56,000 winds determined from the horizontal displacements of EOLE balloons at the 200 mb level in the Southern Hemisphere during the period October 1971-February 1972 is utilized for the computation of planetary- and synoptic-scale kinetic energy space spectra. However, the random distribution of measurements in space and time presents some problems for the spectral analysis. Two different approaches are used, i.e., a harmonic analysis of daily wind values at equi-distant points obtained by space-time interpolation of the data, and a correlation method using the direct measurements. Both methods give similar results for small wavenumbers, but the second is more accurate for higher wavenumbers (k above or equal to 10). The spectra show a maximum at wavenumbers 5 and 6 due to baroclinic instability and then decrease for high wavenumbers up to wavenumber 35 (which is the limit of the analysis), according to the inverse power law k to the negative p, with p close to 3.

Desbois, M.↗

Study of genetic direct search algorithms for function optimization

The results are presented of a study to determine the performance of genetic direct search algorithms in solving function optimization problems arising in the optimal and adaptive control areas. The findings indicate that: (1) genetic algorithms can outperform standard algorithms in multimodal and/or noisy optimization situations, but suffer from lack of gradient exploitation facilities when gradient information can be utilized to guide the search. (2) For large populations, or low dimensional function spaces, mutation is a sufficient operator. However for small populations or high dimensional functions, crossover applied in about equal frequency with mutation is an optimum combination. (3) Complexity, in terms of storage space and running time, is significantly increased when population size is increased or the inversion operator, or the second level adaptation routine is added to the basic structure.

Zeigler, B. P.↗

Potential flow around two-dimensional airfoils using a singular integral method

The problem of potential flow around two-dimensional airfoils is solved by using a new singular integral method. The potential flow equations for incompressible potential flow are written in a singular integral equation. The equation is solved at N collocation points on the airfoil surface. A unique feature of this method is that the airfoil geometry is specified as an independent variable in the exact integral equation. Compared to other numerical methods, the present calculation procedure is much simpler and gives remarkable accuracy for many body shapes. An advantage of the present method is that it allows the inverse design calculation and the results are extremely accurate.

Nguyen, Yves↗

Three-dimensional unsteady lifting surface theory in the subsonic range

The methods of the unsteady lifting surface theory are surveyed. Linearized Euler's equations are simplified by means of a Galileo-Lorentz transformation and a Laplace transformation so that the time and the compressibility of the fluid are limited to two constants. The solutions to this simplified problem are represented as integrals with a differential nucleus; these results in tolerance conditions, for which any exact solution must suffice. It is shown that none of the existing three-dimensional lifting surface theories in subsonic range satisfy these conditions. An oscillating elliptic lifting surface which satisfies the tolerance conditions is calculated through the use of Lame's functions. Numerical examples are calculated for the borderline cases of infinitely stretched elliptic lifting surfaces and of circular lifting surfaces. Out of the harmonic solutions any such temporal changes of the down current are calculated through the use of an inverse Laplace transformation.

Kuessner, H. G.↗

Dynamic Forms: Application to Aircraft Guidance - Part 2

The paper describes a method for guiding a dynamic system through a given set of points. The paradigm is a fully automatic aircraft subject to air traffic control (ATC). The ATC provides a sequence of waypoints through which the aircraft trajectory must pass. The waypoints typically specify time, position, and velocity. The guidance problem is to synthesize a system state trajectory that satisfies both the ATC and aircraft constraints. Complications arise because the controlled process is multidimensional, multiaxis, nonlinear, highly coupled, and the state space is not flat. In addition, there is a multitude of operating modes, which may number in the hundreds. Each such mode defines a distinct state space model of the process by specifying the state space coordinatization, the partition of the controls into active controls and configuration controls, and the output map. Furthermore, mode transitions are required to be smooth. The proposed guidance algorithm is based on the inversion of the pure feedback approximation, followed by correction for the effects of zero dynamics. The paper describes the structure and major modules of the algorithm, and the performance is illustrated by several example aircraft maneuvers.

Meyer, George↗

Optimal two-impulse transfer between specified terminal states of Keplerian orbits

The integrals of motion for an inverse-square force field and a necessary condition for optimality are used in a transfer orbit co-ordinate system to formulate the solution of the optimal two-impulse transfer between fixed position and velocity vectors on Keplerian orbits. In this co-ordinate system, the equations reveal two asymptotes that are useful in analysing the solution. It is shown that there are only two real extremals (both minimums) which are separated by an asymptote. An approximate analytic solution is also obtained by ignoring a quadratic term whose coefficient is approximately zero for a large class of orbit transfer problems.

Jezewski, D. J.↗

A Bayesian approach to microwave precipitation profile retrieval

A multichannel passive microwave precipitation retrieval algorithm is developed. Bayes theorem is used to combine statistical information from numerical cloud models with forward radiative transfer modeling. A multivariate lognormal prior probability distribution contains the covariance information about hydrometeor distribution that resolves the nonuniqueness inherent in the inversion process. Hydrometeor profiles are retrieved by maximizing the posterior probability density for each vector of observations. The hydrometeor profile retrieval method is tested with data from the Advanced Microwave Precipitation Radiometer (10, 19, 37, and 85 GHz) of convection over ocean and land in Florida. The CP-2 multiparameter radar data are used to verify the retrieved profiles. The results show that the method can retrieve approximate hydrometeor profiles, with larger errors over land than water. There is considerably greater accuracy in the retrieval of integrated hydrometeor contents than of profiles. Many of the retrieval errors are traced to problems with the cloud model microphysical information, and future improvements to the algorithm are suggested.

Evans, K. Franklin↗

Charge-exchange plasma generated by an ion thruster

The use of high voltage solar arrays greatly reduces or eliminates power processing requirements in space electric propulsion systems. This use also requires substantial areas of solar array to be at high positive potential relative to space and most of the spacecraft. The charge exchange plasma conducts electrons from the ion beam to such positive surfaces, and thereby electrically load the high voltage solar array. To evaluate this problem, the charge-exchange plasma generated by an ion beam was investigated experimentally. Based upon the experimental data, a simple model was derived for the charge-exchange plasma. This model is conservative in the sense that both the electron/ion density and the electron current density should be equal to, or less than, the preducted value for all directions in the hemisphere upstream of the ion beam direction. Increasing the distance between a positive potential surface (such as a high voltage solar array) and the thruster is the simplest way to control interactions. Both densities and currents vary as the inverse square of this distance.

Kaufman, H. R.↗

Optimal low-thrust takeoff from an orbit about an oblate planet

Future space missions to the outer planets may depend upon the use of low-thrust propulsion systems. As these planets are decidedly oblate, the question of the effect of that oblateness on a low-thrust trajectory is of some interest. In this paper the problem of optimal energy increase is attacked under the assumption that the coefficients for the second zonal harmonic, and the nondimensional thrust acceleration are the same order of magnitude. By means of a two-variable asymptotic expansion technique, a near optimal control program is generated and the first-order uniformly valid approximation for the corresponding trajectory is obtained. Tangential thrust is shown to be a good near-optimal thrust program even in the presence of oblateness effects. The optimal control program is found to be oscillatory and quite similar to the optimal control for energy increase in an inverse square gravitational field.

Jacobson, R. A.↗

Guidance of Nonlinear Systems

The paper describes a method for guiding a dynamic system through a given set of points. The paradigm is a fully automatic aircraft subject to air traffic control (ATC). The ATC provides a sequence of way points through which the aircraft trajectory must pass. The way points typically specify time, position, and velocity. The guidance problem is to synthesize a system state trajectory which satisfies both the ATC and aircraft constraints. Complications arise because the controlled process is multi-dimensional, multi-axis, nonlinear, highly coupled, and the state space is not flat. In addition, there is a multitude of possible operating modes, which may number in the hundreds. Each such mode defines a distinct state space model of the process by specifying the state space coordination, the partition of the controls into active controls and configuration controls, and the output map. Furthermore, mode transitions must be smooth. The guidance algorithm is based on the inversion of the pure feedback approximations, which is followed by iterative corrections for the effects of zero dynamics. The paper describes the structure and modules of the algorithm, and the performance is illustrated by several example aircraft maneuvers.

Meyer, George↗

Guidance of Nonlinear Systems

The paper describes a method for guiding a dynamic system through a given set of points. The paradigm is a fully automatic aircraft subject to air traffic control (ATC). The ATC provides a sequence of way points through which the aircraft trajectory must pass. The way points typically specify time, position, and velocity. The guidance problem is to synthesize a system state trajectory which satisfies both the ATC and aircraft constraints. Complications arise because the controlled process is multi-dimensional, multi-axis, nonlinear, highly coupled, and the state space is not flat. In addition, there is a multitude of possible operating modes, which may number in the hundreds. Each such mode defines a distinct state space model of the process by specifying the state space coordinatization, the partition of the controls into active controls and configuration controls, and the output map. Furthermore, mode transitions must be smooth. The guidance algorithm is based on the inversion of the pure feedback approximations, which is followed by iterative corrections for the effects of zero dynamics. The paper describes the structure and modules of the algorithm, and the performance is illustrated by several example aircraft maneuvers.

Meyer, George↗

A comprehensive comparison between experiment and prediction for a transonic turbulent separated flow

Attempts to predict surface pressure distributions on lifting surfaces have been relatively unsuccessful in the transonic regime when the shock wave is of sufficient strength to produce an extensive region of turbulent separated flow. For these conditions, the viscous flow behavior must be accurately described even to obtain reasonable predictions of surface pressure. The present paper addresses this problem. Detailed comparisons between prediction and experiment are made for a transonic, turbulent boundary-layer separation (freestream Mach number = 0.875) for which the turbulent flow properties (including the turbulent Reynolds stress) had been measured by the laser velocimeter technique from upstream of the separated region through reattachment. The flow was generated on an axisymmetric 'bump' model designed to simulate the flow on an airfoil at transonic conditions. The numerical methods used in the comparisons include the solution of the time-dependent, mass-averaged Navier-Stokes equations, and the solution of the compressible boundary-layer equations by the inverse method. Solutions were obtained for the well established Cebeci-Smith algebraic turbulence model and the more recently developed Wilcox-Rubesin two-equation turbulence model.

Johnson, D. A.↗

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

Recursive mass matrix factorization and inversion: An operator approach to open- and closed-chain multibody dynamics

This report advances a linear operator approach for analyzing the dynamics of systems of joint-connected rigid bodies.It is established that the mass matrix M for such a system can be factored as M=(I+H phi L)D(I+H phi L) sup T. This yields an immediate inversion M sup -1=(I-H psi L) sup T D sup -1 (I-H psi L), where H and phi are given by known link geometric parameters, and L, psi and D are obtained recursively by a spatial discrete-step Kalman filter and by the corresponding Riccati equation associated with this filter. The factors (I+H phi L) and (I-H psi L) are lower triangular matrices which are inverses of each other, and D is a diagonal matrix. This factorization and inversion of the mass matrix leads to recursive algortihms for forward dynamics based on spatially recursive filtering and smoothing. The primary motivation for advancing the operator approach is to provide a better means to formulate, analyze and understand spatial recursions in multibody dynamics. This is achieved because the linear operator notation allows manipulation of the equations of motion using a very high-level analytical framework (a spatial operator algebra) that is easy to understand and use. Detailed lower-level recursive algorithms can readily be obtained for inspection from the expressions involving spatial operators. The report consists of two main sections. In Part 1, the problem of serial chain manipulators is analyzed and solved. Extensions to a closed-chain system formed by multiple manipulators moving a common task object are contained in Part 2. To retain ease of exposition in the report, only these two types of multibody systems are considered. However, the same methods can be easily applied to arbitrary multibody systems formed by a collection of joint-connected regid bodies.

Rodriguez, G.↗

Analytical and phenomenological studies of rotating turbulence

A framework, which combines mathematical analysis, closure theory, and phenomenological treatment, is developed to study the spectral transfer process and reduction of dimensionality in turbulent flows that are subject to rotation. First, we outline a mathematical procedure that is particularly appropriate for problems with two disparate time scales. The approach which is based on the Green's method leads to the Poincare velocity variables and the Poincare transformation when applied to rotating turbulence. The effects of the rotation are now reflected in the modifications to the convolution of a nonlinear term. The Poincare transformed equations are used to obtain a time-dependent analog of the Taylor-Proudman theorem valid in the asymptotic limit when the non-dimensional parameter mu is identical to Omega(t) approaches infinity (Omega is the rotation rate and t is the time). The 'split' of the energy transfer in both direct and inverse directions is established. Secondly, we apply the Eddy-Damped-Quasinormal-Markovian (EDQNM) closure to the Poincare transformed Euler/Navier-Stokes equations. This closure leads to expressions for the spectral energy transfer. In particular, an unique triple velocity decorrelation time is derived with an explicit dependence on the rotation rate. This provides an important input for applying the phenomenological treatment of Zhou. In order to characterize the relative strength of rotation, another non-dimensional number, a spectral Rossby number, which is defined as the ratio of rotation and turbulence time scales, is introduced. Finally, the energy spectrum and the spectral eddy viscosity are deduced.

Mahalov, Alex↗

Spaceborne Radar Measurements of Rainfall Vertical Velocity

This paper studies the performance of a spaceborne precipitation radar in measuring vertical Doppler velocity of rainfall. As far as a downward pointing precipitation radar is concerned, one of the major problems affecting Doppler measurement at the nadir direction arises from the Non-Uniform Beam-Filling effect (NUBF). That is, when significant variation in rain rate is present within the radar IFOV (Instrument Field of View) in the along track direction. the Doppler shift caused by the radial component of the horizontal speed of the satellite is weighted differently among the portions of IFOV. The effects of this non-uniform weighting may dominate any other contribution. Under this condition, shape, average value and width of the Doppler spectrum may not be directly correlated with the vertical velocity of the precipitating particles. However, by using an inversion technique which over-samples the radar measurements in the along track direction, we show that the shift due to NUBF can be evaluated, and that the NUBF induced errors on average fall speed can be reduced.

Im, Eastwood↗

Coronal Diagnostics of Intermediate Activity Star XI Boo A

The analysis of Xi Boo A proved difficult to adapt to our line-by-line approach because of the strong wings of the RGS instrumental profile, as has been detailed in earlier reports. While progress was also delayed because of problems in using SAS v4, we succeeded in the past year or so to bring the analysis to conclusion. Abundances have been derived using both EPIC and RGS data, confirming earlier EUVE findings of a mild solar-like FIP effect, though with some evidence of a turn-up in abundances of elements with higher FIP. Plasma densities appear normal for a moderately active stellar corona. Xi Boo A nicely bridges the gap between the very active stars and stars like the Sun, and it indeed does appear that these are the stars in which the solar-like FIP effects begins to change to the "inverse FIP" type of effect seen in the very active stars. Probing this divide was the main goal of the proposal. These results are in the process of being prepared for publication, though we have not decided the target journal as yet.

Drake, Jeremy↗

A Study of the Effect of Increased Size and Speed of Pursuit Airplanes on the Aileron Balancing Problem

The present trend is toward faster and larger pursuit airplanes. Because both speed and size increase the aileron control forces, the design of ailerons for manual operation is becoming increasingly difficult. In order to obtain a clearer picture of the future problem of balancing ailerons, and inspection has been made of the effects of airplane size and speed on the control forces. Computations were made of the aileron control forces required to meet specified rolling conditions for plain ailerons on wings with spans from 40 to 80 feet and for speeds up to 500 miles per hour. The rolling conditions were specified by two alternative criterions. One was the rolling criterion pb/2V of reference 1. For reasons, which will be discussed later, a value of 0.09 rather than the recommended value of 0.07 was assigned to this criterion. For the criterion pb/2V, the required value of the rolling velocity p varies inversely with the airplane span b. There is some question as to whether the rolling velocity of a pursuit airplane can be permitted to decrease simply because its size is increased. For the second criterion, therefore, the rolling velocity is independent of span (p/V is a constant). The value assigned to this criterion was so chosen that for a wing of 40-foot span the value of pb/2V would be 0.09. The computations neglected compressibility effects. Available experimental data and the results of tests given in reference 2 indicate that the effect of compressibility is to increase the control force. Recent flight tests have indicated that, with certain types of aileron, serious compressibility effects may cause discontinuity at speeds of approximately 400 miles per hour in the aileron control force curves.

Soule, H. A.↗