Search NASA⌕ Search

SEARCH · Search NASA

Results for “Inverse problems”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Practical considerations in waverider applications

This paper proposes a practical perspective for waverider shapes with respect to their use in the commercial and military aerospace industry. Also, integration problems that must be considered when deriving these shapes for operational aircraft are discussed and general approaches to solving them are suggested. The inverse design methodology that is used to generate waverider geometries is discussed as a potentially powerful tool for the aerospace industry in the design of practical operational high-speed (Mach 4 to 25) aircraft. Further, it is suggested that unless this methodology proves to be a valuable tool in the near future, it will not likely find a customer in the aerospace industry and may become no more than an interesting academic study in aerodynamic optimization.

Stevens, Daniel R.↗
NASA NTRS1996

Surface and Atmospheric Parameter Retrieval From AVIRIS Data: The Importance of Non-Linear Effects

AVIRIS data represent a new and important approach for the retrieval of atmospheric and surface parameters from optical remote sensing data. Not only as a test for future space systems, but also as an operational airborne remote sensing system, the development of algorithms to retrieve information from AVIRIS data is an important step to these new approaches and capabilities. Many things have been learned since AVIRIS became operational, and the successive technical improvements in the hardware and the more sophisticated calibration techniques employed have increased the quality of the data to the point of almost meeting optimum user requirements. However, the potential capabilities of imaging spectrometry over the standard multispectral techniques have still not been fully demonstrated. Reasons for this are the technical difficulties in handling the data, the critical aspect of calibration for advanced retrieval methods, and the lack of proper models with which to invert the measured AVIRIS radiances in all the spectral channels. To achieve the potential of imaging spectrometry, these issues must be addressed. In this paper, an algorithm to retrieve information about both atmospheric and surface parameters from AVIRIS data, by using model inversion techniques, is described. Emphasis is put on the derivation of the model itself as well as proper inversion techniques, robust to noise in the data and an inadequate ability of the model to describe natural variability in the data. The problem of non-linear effects is addressed, as it has been demonstrated to be a major source of error in the numerical values retrieved by more simple, linear-based approaches. Non-linear effects are especially critical for the retrieval of surface parameters where both scattering and absorption effects are coupled, as well as in the cases of significant multiple-scattering contributions. However, sophisticated modeling approaches can handle such non-linear effects, which are especially important over vegetated surfaces. All the data used in this study were acquired during the 1991 Multisensor Airborne Campaign (MAC-Europe), as part of the European Field Experiment on a Desertification-threatened Area (EFEDA), carried out in Spain in June-July 1991.

Green Robert O.↗
NASA NTRS2011

Obtaining the Grobner Initialization for the Ground Flash Fraction Retrieval Algorithm

At optical wavelengths and from the vantage point of space, the multiple scattering cloud medium obscures one's view and prevents one from easily determining what flashes strike the ground. However, recent investigations have made some progress examining the (easier, but still difficult) problem of estimating the ground flash fraction in a set of N flashes observed from space In the study by Koshak, a Bayesian inversion method was introduced for retrieving the fraction of ground flashes in a set of flashes observed from a (low earth orbiting or geostationary) satellite lightning imager. The method employed a constrained mixed exponential distribution model to describe the lightning optical measurements. To obtain the optimum model parameters, a scalar function of three variables (one of which is the ground flash fraction) was minimized by a numerical method. This method has formed the basis of a Ground Flash Fraction Retrieval Algorithm (GoFFRA) that is being tested as part of GOES-R GLM risk reduction.

Solakiewicz, R.↗
NASA NTRS2021

Improved Representation of the Global Dust Cycle Using Observational Constraints on Dust Properties and Abundance

Even though desert dust is the most abundant aerosol by mass in Earth’s atmosphere, atmospheric models 25 struggle to accurately represent its spatial and temporal distribution. These model errors are partially caused byfundamental difficulties in simulating dust emission in coarse-resolution models and in accurately representing dust microphysical properties. Here we mitigate these problems by developing a new methodology that yields an improved representation of the global dust cycle. We present an analytical framework that uses inverse modeling to integrate an ensemble of global model simulations with observational constraints on the dust size distribution, 30 extinction efficiency, and regional dust aerosol optical depth. We then compare the inverse model results against independent measurements of dust surface concentration and deposition flux and find that errors are reduced by approximately a factor of two relative to current model simulations of the Northern Hemisphere dust cycle. The inverse model results show smaller improvements in the less dusty Southern Hemisphere, most likely because both the model simulations and the observational constraints used in the inverse model are less accurate. On a global 35 basis, we find that the emission flux of dust with geometric diameter up to 20 μm (PM20) is approximately 5,000 Tg/year, which is greater than most models account for. This larger PM20 dust flux is needed to match observational constraints showing a large atmospheric loading of coarse dust. We obtain gridded data sets of dust emission, vertically integrated loading, dust aerosol optical depth, (surface) concentration, and wet and dry deposition fluxes that are resolved by season and particle size. As our results indicate that this data set is more accurate than current 40 model simulations and the MERRA-2 dust reanalysis product, it can be used to improve quantifications of dust impacts on the Earth system.

aerosols, dust↗
NASA NTRS2021

Improved representation of the global dust cycle using observational constraints on dust properties and abundance

Even though desert dust is the most abundant aerosol by mass in Earth's atmosphere, atmospheric models struggle to accurately represent its spatial and temporal distribution. These model errors are partially caused by fundamental difficulties in simulating dust emission in coarse-resolution models and in accurately representing dust microphysical properties. Here we mitigate these problems by developing a new methodology that yields an improved representation of the global dust cycle. We present an analytical framework that uses inverse modeling to integrate an ensemble of global model simulations with observational constraints on the dust size distribution, extinction efficiency, and regional dust aerosol optical depth. We then compare the inverse model results against independent measurements of dust surface concentration and deposition flux and find that errors are reduced by approximately a factor of 2 relative to current model simulations of the Northern Hemisphere dust cycle. The inverse model results show smaller improvements in the less dusty Southern Hemisphere, most likely because both the model simulations and the observational constraints used in the inverse model are less accurate. On a global basis, we find that the emission flux of dust with a geometric diameter up to 20 µm (PM20) is approximately 5000 Tg yr−1, which is greater than most models account for. This larger PM20 dust flux is needed to match observational constraints showing a large atmospheric loading of coarse dust. We obtain gridded datasets of dust emission, vertically integrated loading, dust aerosol optical depth, (surface) concentration, and wet and dry deposition fluxes that are resolved by season and particle size. As our results indicate that this dataset is more accurate than current model simulations and the MERRA-2 dust reanalysis product, it can be used to improve quantifications of dust impacts on the Earth system.

Dust cycle↗
NASA NTRS1987

The finite-difference time-domain (FD-TD) method for electromagnetic scattering and interaction problems

The formulation and recent applications of the finite-difference time-domain (FD-TD) method for the numerical modeling of electromagnetic scattering and interaction problems are considered. It is shown that improvements in FD-TD modeling concepts and software implementation often make it a preferable choice for structures which cannot be easily treated by conventional integral equations and asymptotic approaches. Recent FD-TD modeling validations in research areas including coupling to wires and wire bundles in free space and cavities, scattering from surfaces in relativistic motion, inverse scattering, and radiation condition theory, are reviewed. Finally, the advantages and disadvantages of FD-TD, and guidelines concerning when FD-TD should and should not be used in high-frequency electromagnetic modeling problems, are summarized.

Taflove, A.↗
NASA NTRS1999

Geodynamic Effects of Ocean Tides: Progress and Problems

Satellite altimetry, particularly Topex/Poseidon, has markedly improved our knowledge of global tides, thereby allowing significant progress on some longstanding problems in geodynamics. This paper reviews some of that progress. Emphasis is given to global-scale problems, particularly those falling within the mandate of the new IERS Special Bureau for Tides: angular momentum, gravitational field, geocenter motion. For this discussion I use primarily the new ocean tide solutions GOT99.2, CSR4.0, and TPXO.4 (for which G. Egbert has computed inverse-theoretic error estimates), and I concentrate on new results in angular momentum and gravity and their solid-earth implications. One example is a new estimate of the effective tidal Q at the M_2 frequency, based on combining these ocean models with tidal estimates from satellite laser ranging. Three especially intractable problems are also addressed: (1) determining long-period tides in the Arctic [large unknown effect on the inertia tensor, particularly for Mf]; (2) determining the global psi_l tide [large unknown effect on interpretations of gravimetry for the near-diurnal free wobble]; and (3) determining radiational tides [large unknown temporal variations at important frequencies]. Problems (2) and (3) are related.

Richard, Ray↗
NASA NTRS1982

On the inversion of block tridiagonals without storage constraints

A strategy was developed to permit trade-offs between the number of floating point operations required and the storage requirements for the solution of certain difference problems, such as block tridiagonal systems of equations. This is done by recomputing some intermediate results instead of storing them. Reducing the storage to the square root of the current requirement roughly doubles the number of computations. Reducing the storage more than this tends to make the number of computations prohibitively large. In theory, if m is the order of each sub-matrix in the block tridiagonal matrix, one can solve any linear system with only 5m(2) + 1 temporary storage cells. In many cases m is a constant and quite small. For example, in solving a factored form of the three-dimensional Navier-Stokes equations, the size m of the block tridiagonals is 5. This method lends itself to efficient use on computers with parallel processing or vector processing architectures. On these computers the larger number of floating point operations is more than offset by the decrease in I/O and the increased percentage of vector operations made possible by this algorithm.

Merriam, M. L.↗
NASA NTRS1972

A review of nonstatistical techniques for the estimation of vertical atmospheric structure from remote infrared measurements

Mathematical techniques used in atmospheric profile inversion are reviewed, with the discussion limited to methods that do not require extensive information on the statistical properties of the profiles. The basic problem of the calculation of vertical profiles from noise contaminated infrared spectral measurements is examined. Three specific examples of methods of solution are given: the iterative minimum estimation method, a second linear iterative method, and the algorithm of Chahine. Satellite data are used in an empirical comparison of the methods. Several factors are considered, including the stability of solutions, convergence behavior of the iterative methods, and dependence of the solutions on the first guess, a factor that is particularly important in the treatment of data from other planets.

Conrath, B. J.↗
NASA NTRS1981

Atmospheric structure of the outer planets from thermal emission data

Methods for determining atmospheric structure exploit the opacities provided by the collision induced H2 dipole and the nu4 fundamental of CH4. In addition to earth-based observations, useful measurements of thermal emission from Jupiter and Saturn have been or soon will be made by several spacecraft, with results cross-checked with independent radio occultaion results. For Uranus and Neptune, only a limited set of whole-disk earth-based data exists. All the outer planets show evidence for stratospheric temperature inversions; temperature minima range from about 105 K for Jupiter and 87 K for Saturn, to roughly 55 K for Uranus and Neptune. Remaining problems may be resolved by better quantitative understanding of gas and aerosol absorption and scattering properties, chemical composition, and non-LTE source functions. Ultimately, temperature structure results must be supplemented by quantitative energy equilibrium models which will allow some meaning to be given to the relationships between such characteristics as temperature, clouds, incident solar and planetary radiation, and chemical composition.

Orton, G. S.↗
NASA NTRS1985

Optimal Sensor Locations for Structural Identification

The optimum sensor location problem, OSLP, may be thought of in terms of the set of systems, S, the class of input time functions, I, and the identification algorithm (estimator) used, E. Thus, for a given time history of input, the technique of determining the OSL requires, in general, the solution of the optimization and the identification problems simultaneously. A technique which uncouples the two problems is introduced. This is done by means of the concept of an efficient estimator for which the covariance of the parameter estimates is inversely proportional to the Fisher Information Matrix.

Udwadia, F. E.↗
NASA NTRS1988

Random field estimation approach to multibody dynamics

The use of recursive techniques similar to random field models to factor and invert the inertia matrix of a multibody system is discussed. An equivalence is established between the composite multibody system inertia matrix and the covariance of the output of a described linear system model. Conditional mean estimation and sequential estimation problems are solved along with problems of filtering and smoothing. Formulas are developed to compute the covariance of several relevant quantities. The foregoing results are used to obtain the inverse of the composite multibody system inertia in closed form.

Rodriguez, G.↗
NASA NTRS1941

On the calculation of flow past an infinite screen of thin airfoils

This report deals with the flow past an infinite screen of thin airfoil (two-dimensional problem). The vortex distribution across the profile is established with appropriate expansion in series and the velocity distribution lift, moment, and profile shape deduced. Inversely, the distribution is deduced from the vorticity. The method is the extension of the Birnbaum-Glauert method for the isolated wing.

Pistolesi, E↗
NASA NTRS1993

Noncolocated Structural Vibration Suppression Using Zero Annihilation Periodic Control

The Zero Annihilation Periodic (ZAP) controller is applied to the problem of vibration control of a noncolocated flexible structure. It is shown that even though the transfer function is nonminimum-phase, a plant inverse controller can be designed which elicits a deadbeat closed-loop response.

periodic control noncolocated systems non-minimum ↗
NASA NTRS1972

Population inversion calculations using near resonant charge exchange as a pumping mechanism

Near resonance charge exchange between ions of a large ionization potential gas such as helium or neon and vapors of metals such as zinc, cadmium, selenium, or tellurium has produced laser action in the metal ion gas. The possibility of obtaining population inversions in near resonant charge exchange systems (Xe-Ca, Xe-Mg, Xe-Sr, Xe-Ba, Ar-Mg, N-Ca) was investigated. The analysis is an initial value problem that utilizes rate equations for the densities of relevant levels of the laser gas (Ca, Ba, Mg, or Sr) and an electron energy equation. Electron excitation rates are calculated using the Bohr-Thomson approximation for the cross section. Approximations to experimental values of the electron ionization cross section and the ion-atom charge exchange cross section are used. Preliminary results have been obtained for the Ca-Xe system and show that it is possible to obtain gains greater than 10 to the 14th power/m with inversion times up to 8x10 to the minus 7th power second. A possible charge exchange laser system using a MPD arc plasma accelerator is also described.

Chubb, D. L.↗
NASA NTRS2002

Containerless Ripple Turbulence

One of the longest standing unsolved problems in physics relates to the behavior of fluids that are driven far from equilibrium such as occurs when they become turbulent due to fast flow through a grid or tidal motions. In turbulent flows the distribution of vortex energy as a function of the inverse length scale [or wavenumber 'k'] of motion is proportional to 1/k(sup 5/3) which is the celebrated law of Kolmogorov. Although this law gives a good description of the average motion, fluctuations around the average are huge. This stands in contrast with thermally activated motion where large fluctuations around thermal equilibrium are highly unfavorable. The problem of turbulence is the problem of understanding why large fluctuations are so prevalent which is also called the problem of 'intermittency'. Turbulence is a remarkable problem in that its solution sits simultaneously at the forefront of physics, mathematics, engineering and computer science. A recent conference [March 2002] on 'Statistical Hydrodynamics' organized by the Los Alamos Laboratory Center for Nonlinear Studies brought together researchers in all of these fields. Although turbulence is generally thought to be described by the Navier-Stokes Equations of fluid mechanics the solution as well as its existence has eluded researchers for over 100 years. In fact proof of the existence of such a solution qualifies for a 1 M$ millennium prize. As part of our NASA funded research we have proposed building a bridge between vortex turbulence and wave turbulence. The latter occurs when high amplitude waves of various wavelengths are allowed to mutually interact in a fluid. In particular we have proposed measuring the interaction of ripples [capillary waves] that run around on the surface of a fluid sphere suspended in a microgravity environment. The problem of ripple turbulence poses similar mathematical challenges to the problem of vortex turbulence. The waves can have a high amplitude and a strong nonlinear interaction. Furthermore, the steady state distribution of energy again follows a Kolmogorov scaling law; in this case the ripple energy is distributed according to 1/k (sup 7/4). Again, in parallel with vortex turbulence ripple turbulence exhibits intermittency. The problem of ripple turbulence presents an experimental opportunity to generate data in a controlled, benchmarked system. In particular the surface of a sphere is an ideal environment to study ripple turbulence. Waves run around the sphere and interact with each other, and the effect of walls is eliminated. In microgravity this state can be realized for over 2 decades of frequency. Wave turbulence is a physically relevant problem in its own right. It has been studied on the surface of liquid hydrogen and its application to Alfven waves in space is a source of debate. Of course, application of wave turbulence perspectives to ocean waves has been a major success. The experiment which we plan to run in microgravity is conceptually straightforward. Ripples are excited on the surface of a spherical drop of fluid and then their amplitude is recorded with appropriate photography. A key challenge is posed by the need to stably position a 10cm diameter sphere of water in microgravity. Two methods are being developed. Orbitec is using controlled puffs of air from at least 6 independent directions to provided the positioning force. This approach has actually succeeded to position and stabilize a 4cm sphere during a KC 135 segment. Guigne International is using the radiation pressure of high frequency sound. These transducers have been organized into a device in the shape of a dodecahedron. This apparatus 'SPACE DRUMS' has already been approved for use for combustion synthesis experiments on the International Space Station. A key opportunity presented by the ripple turbulence data is its use in driving the development of codes to simulate its properties.

Putterman, Seth↗
NASA NTRS1983

Integrals of motion for the classical two-body problem with drag

Integrals of motion for the two-body problem with drag are obtained by operating on the second-order vector differential equation describing the motion. The force field consists of an inverse-square gravitational attraction and a drag force proportional to the velocity vector and inversely proportional to the square of the distance to the attracting center. The developed integrals are the analogs of the Keplerian scalar energy, the vector angular momentum, and the Laplace vector.

Jezewski, D. J.↗