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 865 records · Page 48

Adaptive Flight Control Design with Optimal Control Modification on an F-18 Aircraft Model

In the presence of large uncertainties, a control system needs to be able to adapt rapidly to regain performance. Fast adaptation is referred to as the implementation of adaptive control with a large adaptive gain to reduce the tracking error rapidly; however, a large adaptive gain can lead to high-frequency oscillations which can adversely affect the robustness of an adaptive control law. A new adaptive control modification is presented that can achieve robust adaptation with a large adaptive gain without incurring high-frequency oscillations as with the standard model-reference adaptive control. The modification is based on the minimization of the Y2 norm of the tracking error, which is formulated as an optimal control problem. The optimality condition is used to derive the modification using the gradient method. The optimal control modification results in a stable adaptation and allows a large adaptive gain to be used for better tracking while providing sufficient robustness. A damping term (v) is added in the modification to increase damping as needed. Simulations were conducted on a damaged F-18 aircraft (McDonnell Douglas, now The Boeing Company, Chicago, Illinois) with both the standard baseline dynamic inversion controller and the adaptive optimal control modification technique. The results demonstrate the effectiveness of the proposed modification in tracking a reference model.

Burken, John J.↗

Optimal Aerosol Parameterization for Remote Sensing Retrievals

We have developed a new algorithm for the retrieval of aerosol and gases from SAGE It1 solar transmission measurements. This algorithm improves upon the NASA operational algorithm in several key aspects, including solving the problem non-linearly and incorporating a new methodology for separating the contribution of aerosols and gases. In order to extract aerosol information we have built a huge database of aerosol models for both stratospheric and tropospheric aerosols, and polar stratospheric cloud particles. This set of models allows us to calculate a vast range of possible extinction spectra for aerosols. and from these, derive a set of eigenvectors which then provide the basis set used in our inversion algorithm. Our aerosol algorithm and retrievals are described in several articles (listed in References Section) published under this grant. In particular they allow us to analyze the spectral properties of aerosols and PSCs and ultimately derive their microphysical properties. We have found some considerable differences between our spectra and the ones derived from the SAGE III operational algorithm. These are interesting as they provide an independent check on the validity of published aerosol data and, in particular, on their associated uncertainties. In order to understand these differences, we are assembling independent aerosol data from other sources with which to make comparisons. We have carried out extensive comparisons of our ozone retrievals with both SAGE III and independent lidar, ozonesonde, and satellite measurements (Polyakov et al., 2004). These show very good agreement throughout the stratosphere and help to quantify differences which can be attributed to natural variation in ozone versus that produced by algorithmic differences. In the mid - upper stratosphere, agreement with independent data was generally within 5 - 20%. but in the lower stratosphere the differences were considerably larger. We believe that a large proportion of this discrepancy in the lower stratosphere is attributable to natural variation, and is also seen in comparisons between lidar and ozonesonde measurements. NO2 profiles obtained with our algorithm were compared to those obtained through the SAGE III operational algorithm and exhibited differences of 20 - 40%. Our retrieved profiles agree with the HALOE NO2 measurements significantly better than those of the operational retrieval. In other work (described below), we are extending our aerosol retrievals into the infrared regime and plan to perform retrievals from combined uv-visible-infrared spectra. This work will allow us to use the spectra to derive the size and composition of aerosols, and we plan to employ our algorithms in the analysis of PSC spectra. We are presently also developing a limb-scattering algorithm to retrieve aerosol data from limb measurements of solar scattered radiation.

Newchurch, Michael J.↗

Estimates of trapped radiation encountered on low-thrust trajectories through the Van Allen belts

Estimates were made of the number of trapped protons and electrons encountered by vehicles on low-thrust trajectories through the Van Allen belts. The estimates serve as a first step in assessing whether these radiations present a problem to on-board sensitive components and payload. The integrated proton spectra and electron spectra are presented for the case of a trajectory described by a vehicle with a constant-thrust acceleration A sub c equal to 0.001 meter/sq sec. This value of acceleration corresponds to a trip time of about 54 days from low earth orbit to synchronous orbit. It is shown that the time spent in the belts and hence the radiation encountered vary nearly inversely with the value of thrust acceleration. Thus, the integrated spectral values presented for the case of A sub c = 0.001 meter/sq sec can be generalized for any other value of thrust acceleration by multiplying them by the factor 0.001/A sub c.

Karp, I. M.↗

The elastic layer with a cylindrical hole subjected to a nonuniform axisymmetric radial displacement.

A problem in the linear theory of elasticity is considered wherein a layer with a circular cylindrical hole is subjected to a nonuniform axisymmetric radial displacement. The solution utilizes Navier's equations of elasticity which are solved by means of extended Hankel transforms. A special case in which the radial displacement is a linear function of the axial coordinate is presented. Numerical results are given in graphical form for the case when hole radius and layer thickness are equal. The inversion integrals were evaluated numerically using Longman's technique for computing infinite integrals of oscillatory functions.

Grissom, D. S.↗

Mass, heat and nutrient fluxes in the Atlantic Ocean determined by inverse methods

Inverse methods are applied to historical hydrographic data to address two aspects of the general circulation of the Atlantic Ocean. The method allows conservation statements for mass and other properties, along with a variety of other constraints, to be combined in a dynamically consistent way to estimate the absolute velocity field and associated property transports. The method was first used to examine the exchange of mass and heat between the South Atlantic and the neighboring ocean basins. The second problem addressed concerns the circulation and property fluxes across the 24 and 36 deg N in the subtropical North Atlantic. Conservation statements are considered for the nutrients as well as mass, and the nutrients are found to contribute significant information independent of temperature and salinity.

Rintoul, Stephen Rich↗

A spatial operator algebra for manipulator modeling and control

A recently developed spatial operator algebra, useful for modeling, control, and trajectory design of manipulators is discussed. The elements of this algebra are linear operators whose domain and range spaces consist of forces, moments, velocities, and accelerations. The effect of these operators is equivalent to a spatial recursion along the span of a manipulator. Inversion of operators can be efficiently obtained via techniques of recursive filtering and smoothing. The operator algebra provides a high level framework for describing the dynamic and kinematic behavior of a manipulator and control and trajectory design algorithms. The interpretation of expressions within the algebraic framework leads to enhanced conceptual and physical understanding of manipulator dynamics and kinematics. Furthermore, implementable recursive algorithms can be immediately derived from the abstract operator expressions by inspection. Thus, the transition from an abstract problem formulation and solution to the detailed mechanizaton of specific algorithms is greatly simplified. The analytical formulation of the operator algebra, as well as its implementation in the Ada programming language are discussed.

Rodriguez, G.↗

Trajectory optimization for the National Aerospace Plane

The objective of this second phase research is to investigate the optimal ascent trajectory for the National Aerospace Plane (NASP) from runway take-off to orbital insertion and address the unique problems associated with the hypersonic flight trajectory optimization. The trajectory optimization problem for an aerospace plane is a highly challenging problem because of the complexity involved. Previous work has been successful in obtaining sub-optimal trajectories by using energy-state approximation and time-scale decomposition techniques. But it is known that the energy-state approximation is not valid in certain portions of the trajectory. This research aims at employing full dynamics of the aerospace plane and emphasizing direct trajectory optimization methods. The major accomplishments of this research include the first-time development of an inverse dynamics approach in trajectory optimization which enables us to generate optimal trajectories for the aerospace plane efficiently and reliably, and general analytical solutions to constrained hypersonic trajectories that has wide application in trajectory optimization as well as in guidance and flight dynamics. Optimal trajectories in abort landing and ascent augmented with rocket propulsion and thrust vectoring control were also investigated. Motivated by this study, a new global trajectory optimization tool using continuous simulated annealing and a nonlinear predictive feedback guidance law have been under investigation and some promising results have been obtained, which may well lead to more significant development and application in the near future.

Lu, Ping↗

Towards a theory of tropical/midlatitude mass exchange from the earth's surface through the stratosphere

The main focus of this work is to understand the dynamics of mass exchange between the tropics and the midlatitudes and to determine any links between tropospheric exchange and that in the stratosphere. We have approached this problem from two different perspectives. The first is aimed towards understanding the troposphere's role in inducing lower stratospheric tropical/midlatitude exchange. For this project we focus on observational analysis of the lower stratosphere to assess the key regions of transport in/out of the tropics and to what extent this transport is driven by tropospheric processes. The second approach is to determine the extent to which stratospheric processes influence the troposphere. In this project we are performing potential vorticity (PV) inversions to assess the winds induced near the tropopause when the stratospheric polar vortex is displaced equatorward. These are each discussed in more detail in the subsections below. Also, we have organized a session on Tropical/Midlatitude Interaction and Transport at the Fall AGU where we will be showing our latest results.

Hartley, Dana↗

Unsteady aerodynamic modeling and active aeroelastic control

Unsteady aerodynamic modeling techniques are developed and applied to the study of active control of elastic vehicles. The problem of active control of a supercritical flutter mode poses a definite design goal stability, and is treated in detail. The transfer functions relating the arbitrary airfoil motions to the airloads are derived from the Laplace transforms of the linearized airload expressions for incompressible two dimensional flow. The transfer function relating the motions to the circulatory part of these loads is recognized as the Theodorsen function extended to complex values of reduced frequency, and is termed the generalized Theodorsen function. Inversion of the Laplace transforms yields exact transient airloads and airfoil motions. Exact root loci of aeroelastic modes are calculated, providing quantitative information regarding subcritical and supercritical flutter conditions.

Edwards, J. W.↗

Cluster Analysis of Spectroscopic Line Profiles and EUV Emission in RMHD Simulations and Observations of the Solar Atmosphere

Spatially-resolved observations from the IRIS, SDO/AIA, and other space mission and ground-based telescopes, coupled with realistic 3D RMHD simulations, are a powerful tool for analysis of processes in the solar atmosphere. To better understand the dynamical and thermodynamic properties in the simulation data and their connection to observations, it is essential to determine similarities in the behaviors of the synthesized and observed emission. However, the complexity of observational data and physical processes makes comparison of observations and modeling results difficult. In this work, we show the initial results of application of K-Means clustering (unsupervised machine learning) algorithm to two different problems: 1) recognition of the typical spectroscopic line profiles observed by IRIS during solar flares and their typical dynamic behavior; 2) recognition of shocks and heating events in synthetic AIA emission data obtained from StellarBox quiet-Sun simulations. The average silhouette width technique for the KMeans algorithm is utilized in different ways to obtain optimal numbers of clusters. We discuss application of the emission clustering to visualizations of the computational volume, understanding its evolutionary trends and behavior patterns, and inversion (reconstruction) of physical properties of the solar atmosphere from synthesizes emission data.

Sadykov, Viacheslav↗

Mixed finite-difference scheme for free vibration analysis of noncircular cylinders

A mixed finite-difference scheme is presented for the free-vibration analysis of simply supported closed noncircular cylindrical shells. The problem is formulated in terms of eight first-order differential equations in the circumferential coordinate which possess a symmetric coefficient matrix and are free of the derivatives of the elastic and geometric characteristics of the shell. In the finite-difference discretization, two interlacing grids are used for the different fundamental unknowns in such a way as to avoid averaging in the difference-quotient expressions used for the first derivative. The resulting finite-difference equations are symmetric. The inverse-power method is used for obtaining the eigenvalues and eigenvectors.

Noor, A. K.↗

A New, More Physically Based Algorithm, for Retrieving Aerosol Properties over Land from MODIS

The MOD Imaging Spectrometer (MODIS) has been successfully retrieving aerosol properties, beginning in early 2000 from Terra and from mid 2002 from Aqua. Over land, the retrieval algorithm makes use of three MODIS channels, in the blue, red and infrared wavelengths. As part of the validation exercises, retrieved spectral aerosol optical thickness (AOT) has been compared via scatterplots against spectral AOT measured by the global Aerosol Robotic NETwork (AERONET). On one hand, global and long term validation looks promising, with two-thirds (average plus and minus one standard deviation) of all points falling between published expected error bars. On the other hand, regression of these points shows a positive y-offset and a slope less than 1.0. For individual regions, such as along the U.S. East Coast, the offset and slope are even worse. Here, we introduce an overhaul of the algorithm for retrieving aerosol properties over land. Some well-known weaknesses in the current aerosol retrieval from MODIS include: a) rigid assumptions about the underlying surface reflectance, b) limited aerosol models to choose from, c) simplified (scalar) radiative transfer (RT) calculations used to simulate satellite observations, and d) assumption that aerosol is transparent in the infrared channel. The new algorithm attempts to address all four problems: a) The new algorithm will include surface type information, instead of fixed ratios of the reflectance in the visible channels to the mid-IR reflectance. b) It will include updated aerosol optical properties to reflect the growing aerosol retrieved from eight-plus years of AERONE". operation. c) The effects of polarization will be including using vector RT calculations. d) Most importantly, the new algorithm does not assume that aerosol is transparent in the infrared channel. It will be an inversion of reflectance observed in the three channels (blue, red, and infrared), rather than iterative single channel retrievals. Thus, this new formulation of the MODIS aerosol retrieval over land includes more physically based surface, aerosol and radiative transfer with fewer potentially erroneous assumptions.

Levy, Robert C.↗

The Spectral Properties of Gamma-ray Bursts: a Review of Recent Development

Developments in the spectroscopy of gamma ray bursts (GRB) are reviewed. The general question of the validity of the spectral results, particularly with regard to features in the spectrum, is discussed. Confirmations of these spectral features are summarized. Results from the KONUS experiments on Venera 13 and 14 are reviewed. The status of models of the continuum spectrum is summarized. A number of different radiation mechanisms appear capable of fitting the data. These include thermal bremsstrahlung, thermal synchrotron and inverse Compton. Rapid variability of the spectra shape on time scales 76] 0.25 sec. was reported. The characteristic energy of the spectrum was observed to vary over nearly an order of magnitude during individual events. A strong correlation between spectral hardness and luminosity was found. Low-energy (50 keV) absorption features and high-energy (400 keV) emission features continue to appear in GRB spectra. Understanding the origin of these lines in the context of the existing continuum models remains a difficult problem.

Teegarden, B. J.↗

An entropy method for induced drag minimization

A fundamentally new approach to the aircraft minimum induced drag problem is presented. The method, a 'viscous lifting line', is based on the minimum entropy production principle and does not require the planar wake assumption. An approximate, closed form solution is obtained for several wing configurations including a comparison of wing extension, winglets, and in-plane wing sweep, with and without a constraint on wing-root bending moment. Like the classical lifting-line theory, this theory predicts that induced drag is proportional to the square of the lift coefficient and inversely proportioinal to the wing aspect ratio. Unlike the classical theory, it predicts that induced drag is Reynolds number dependent and that the optimum spanwise circulation distribution is non-elliptic.

Greene, George C.↗

Gamma guidance of trajectories for coplanar, aeroassisted orbital transfer

The optimization and guidance of trajectories for coplaner, aeroassisted orbital transfer (AOT) from high Earth orbit (HEO) to low Earth orbit (LEO) are examined. In particular, HEO can be a geosynchronous Earth orbit (GEO). It is assumed that the initial and final orbits are circular, that the gravitational field is central and is governed by the inverse square law, and that at most three impulses are employed: one at HEO exit, one at atmospheric exit, and one at LEO entry. It is also assumed that, during the atmospheric pass, the trajectory is controlled via the lift coefficient. The presence of upper and lower bounds on the lift coefficient is considered. First, optimal trajectories are computed by minimizing the total velocity impulse (hence, the propellant consumption) required for AOT transfer. The sequential gradient-restoration algorithm (SGRA) is used for optimal control problems. The optimal trajectory is shown to include two branches: a relatively short descending flight branch (branch 1) and a long ascending flight branch (branch 2). Next, attention is focused on guidance trajectories capable of approximating the optimal trajectories in real time, while retaining the essential characteristics of simplicity, ease of implementation, and reliability. For the atmospheric pass, a feedback control scheme is employed and the lift coefficient is adjusted according to a two-stage gamma guidance law. Further improvements are possible via a modified gamma guidance which is more stable with respect to dispersion effects arising from navigation errors, variations of the atmospheric density, and uncertainties in the aerodynamic coefficients than gamma guidance trajectory. A byproduct of the studies on dispersion effects is the following design concept. For coplaner aeroassisted orbital transfer, the lift-range-to-weight ratio appears to play a more important role than the lift-to-drag ratio. This is because the lift-range-to-weight ratio controls mainly the minimum altitude (hence, the peak heating rate) of the guidance trajectory; on the other hand, the lift-to-drag ratio controls mainly the duration of the atmospheric pass of the guidance trajectory.

Miele, A.↗

The investigation of advanced remote sensing techniques for the measurement of aerosol characteristics

Advanced remote sensing techniques and inversion methods for the measurement of characteristics of aerosol and gaseous species in the atmosphere were investigated. Of particular interest were the physical and chemical properties of aerosols, such as their size distribution, number concentration, and complex refractive index, and the vertical distribution of these properties on a local as well as global scale. Remote sensing techniques for monitoring of tropospheric aerosols were developed as well as satellite monitoring of upper tropospheric and stratospheric aerosols. Computer programs were developed for solving multiple scattering and radiative transfer problems, as well as inversion/retrieval problems. A necessary aspect of these efforts was to develop models of aerosol properties.

Deepak, A.↗

On the Retrieval of Lightning Radio Sources from Time-of-Arrival Data

We examine the problem of retrieving three-dimensional lightning locations from radio frequency Time-Of-Arrival (TOA) measurements. Arbitrary antenna locations are considered. By judiciously differencing measurements that are related to the location of the antennas and their excitation times, the problem is converted from the initial spherical nonlinear form to a system of linear equations. In the linear formalism, the source location and time-of-occurrence is viewed geometrically as an intersection of hyperplanes in the four-dimensional Minkowski space (x,y,z,t). The linear equations are solved to obtain explicit analytic expressions for the location and time variables. Retrieval errors are not interpreted with conventional Geometrical Dilution of Precision (GDOP) arguments as discussed by Holmes and Reedy (1951), but with more recent inversion analyses considered by Twomey (1977). Measurement errors are propagated analytically so that the specific effect of these errors on the solution is clarified. The sensitivity of the solution on the number of antennas used, antenna network geometry, source position, and measurement differencing schemes are discussed in terms of the eigenvalues of the linear system.

Koshak, William J.↗

Current Induced Forces on Adsorbates at Conducting Carbon Nanotubes

The problem addressed is force on an impurity due to electric current. In order to determine the size and direction of the force, as well as its dependence on the voltage, LCAO (Linear Combination of Atomic Orbitals) and the Green's Function approach were implemented. The resultant induced mean force was directed outwards in addition to a smaller component in the electron flow direction. At low bias, force is larger for larger radius, and crosses over the opposite trend when increasing the bias. The crossover current is approximately inversely proportional to the radius.

Mingo, N.↗