Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inversion methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

Experimental Flight Validation of the Prandtl 1933 Bell Spanload

This report describes the validation of the 1933 Prandtl bell spanload. This spanload is the minimum induced drag of a wing for a given structural weight with properties that eliminate adverse yaw. Aircraft using the Prandtl bell spanload were flown and investigated. The results of this research show that many previously held assumptions should be rethought, and the creation of aircraft using the Prandtl bell spanload will require considerable new techniques. Part of this work centered on the use of inverse methods. The usual first-step approach to a computational fluids problem is to create a geometry of the aircraft. Once this geometry exists, the computational fluids solution has been solved; however, the problem of creating the geometry still exists. A very small segment of the computational fluids world has concentrated on inverse solutions. This design approach begins with an end result - the computational fluids solution; from this end result, the geometry is sought. A more generic inverse tool was desired that would allow for the design of wings - specifically Prandtl bell spanload wings. Such a tool has been developed and allows for variable taper, aspect ratio, sweep, airfoils, and design-lift coefficients; output from the tool results in twist distribution of wings.

Albion H. Bowers↗

A class of airfoils having finite trailing edge pressure gradients

Some new developments relevant to the design of single-element airfoils using potential flow methods are presented. In particular, the ramifications of the unbounded trailing edge pressure gradients generally present in the potential flow solution for the flow over an airfoil are examined, and the conditions necessary to obtain a class of airfoils having finite trailing edge pressure gradients developed. The incorporation of these conditions into the inverse method of Eppler for the design of low-speed airfoils is discussed, and designs generated using the modified scheme are presented for consideration. A detailed viscous analysis of one of these airfoils demonstrates a significant reduction in the strong inviscid-viscid interactions generally present near the trailing edge. These reductions offer the possibility of improved airfoil performance, as well as the possibility of improved accuracy in the methods of airfoil design and analysis.

Ormsbee, A. I.↗

Development of direct-inverse 3-D methods for applied transonic aerodynamic wing design and analysis

Progress in the direct-inverse wing design method in curvilinear coordinates has been made. This includes the remedying of a spanwise oscillation problem and the assessment of grid skewness, viscous interaction, and the initial airfoil section on the final design. It was found that, in response to the spanwise oscillation problem that designing at every other spanwise station produced the best results for the cases presented, a smoothly varying grid is especially needed for the accurate design at the wing tip, the boundary layer displacement thicknesses must be included in a successful wing design, the design of high and medium aspect ratio wings is possible with this code, and the final airfoil section designed is fairly independent of the initial section.

Carlson, Leland A.↗

Real-Time Adaptive Drag Minimization Wind Tunnel Investigation of a Flexible Wing with Variable Camber Continuous Trailing Edge Flap System

This paper reports the results of a recently completed real-time adaptive drag minimization wind tunnel investigation of a highly flexible wing wind tunnel model equipped with the Variable Camber Continuous Trailing Flap (VCCTEF) technology at the University of Washington Aeronautical Laboratory (UWAL). The wind tunnel investigation is funded by NASA SBIR Phase II contract with Scientific Systems Company, Inc. (SSCI) and University of Washington (UW) as a subcontractor. The wind tunnel model is a sub-scale Common Research Model (CRM) wing constructed of foam core and fiberglass skin and is aeroelastically scaled to achieve a wing tip deflection of 10% of the wing semi-span which represents a typical wing tip deflection for a modern transport such as Boeing 787. The jig-shape twist of the CRM wing is optimized using a CART3D aero-structural model to achieve the minimum induced drag for the design cruise lift coefficient of 0.5. The wing is equipped with two chord wise cambered segments for each of the six span wise flap sections for a total of 12 individual flap segments that comprise the VCCTEF system. Each of the 12 flap segments is actively controlled by an electric servo-actuator. The real-time adaptive drag optimization strategy includes an on-board aerodynamic model identification, a model excitation, and a real-time drag optimization. The on-board aerodynamic model is constructed parametrically as a function of the angle of attack and flap positions to model the lift and drag coefficients of the wing. The lift coefficient models include a linear model and a second-order model. The drag coefficient models include a quadratic model and a higher-order up to 6th-order model to accurately model the drag coefficient at high angles of attack. The onboard aerodynamic model identification includes a recursive least-squares (RLS) algorithm and a batch least-squares (BLS) algorithm designed to estimate the model parameters. The model excitation method is designed to sample the input set that comprises the angle of attack and the flap positions. Three model excitation methods are developed: random excitation method, sweep method, and iterative angle-of-attack seeking method. The real-time drag optimization includes a generic algorithm developed by SSCI and several optimization methods developed by NASA which include a second-order gradient Newton-Raphson optimization method, an iterative gradient optimization method, a pseudo-inverse optimization method, an analytical optimization method, and an iterative refinement optimization method. The first wind tunnel test entry took place in September 2017. This test revealed major hardware issues and required further redesign of the flap servo mechanisms. The second test entry took place in April 2018. However, the test was not successful due to the issues with the onboard aerodynamic model identification RLS algorithm which incorrectly identified model parameters. This test also provides an experimental comparison study between the VCCTEF and a variable camber discrete trailing edge flap (VCDTEF) without the elastomer transition mechanisms. The experimental result confirms the benefit of the VCCTEF which produces lower drag by 5% than the VCDTEF. The third and final test entry took place in June 2018 after the issues with the RLS algorithm have been identified and corrected. Additional improvements were implemented. These include the BLS algorithm, the iterative angle-of-attack seeking method, the iterative gradient optimization method, and the pseudo-inverse optimization method. The test objectives were successfully demonstrated as the real-time drag optimization identifies several optimal solutions at off-design lift coefficients. The iterative gradient optimization method is found to achieve up to 4.7% drag reduction for the off-design lift coefficient of 0.7. The pseudo-inverse optimization method which does not require the drag coefficient model is found to be quite effective in reducing drag. Up to 9.4% drag reduction for the off-design lift coefficient of 0.7 is achieved with the pseudo-inverse optimization method. The wind tunnel investigation demonstrates the potential of real-time drag optimization technology. Several new capabilities are developed that could enable future adaptive wing technologies for flexible wings equipped with drag control devices such as the VCCTEF.

Drag Optimization↗

Remote sensing of ozone in the middle ultraviolet

Earlier studies of the remote sensing of the ozone distribution from satellites in the middle ultraviolet are extended, giving particular attention to the analytic model method. After a brief introduction, the historical background of this subject is presented. This is followed by a summary and review of early work based upon the analytic model. Then a linearized approach to the problem and a linearized equation for a density correction are discussed. These indicate some of the general limitations of the inversion method. Finally the parametrized model approach is reexamined by using both simulated data and some actual satellite data.

Green, A. E. S.↗

A multi-grid method for the computation of viscid/inviscid interactions on airfoils

An improved version of the 'GRUMFOIL' code has been developed for the computation of airfoil flows. The method employs a conservative difference scheme for the potential equation, integral methods for the boundary layer, and viscous coupling conditions that fully account for the wake and strong interaction effects at trailing edges. The improved version uses Jameson's 'MAD' scheme to accelerate convergence of the inviscid solution, an improved 2nd order artificial viscosity and far field 'TARE' correction to reduce spatial truncation errors and Carter's semi-inverse method for the viscous solution. Results are presented which demonstrate a factor of ten reduction in computing cost.

Melnik, R. E.↗

Synoptic Gulf Stream velocity profiles through simultaneous inversion of hydrographic and acoustic Doppler data

Data from a shipborne acoustic profiling device have been combined with conductivity, temperature, depth/O2 sections across the Gulf Stream to form estimates of the absolute flow fields. The procedure for the combination was a form of inverse method. The results suggest that at the time of the observations (June 1982) the net Gulf Stream transport off Hatteras was 107 + or - 11 Sv and that across a section near 72.5 W it had increased to 125 + or - 6 Sv. The transport of the deep western boundary current was 9 + or - 3 Sv. For comparison purposes an inversion was done using the hydrographic/O2 data alone as in previously published results and obtained qualitative agreement with the combined inversion. Inversion of the acoustic measurements alone, when corrected for instrument biases, leaves unacceptably large mass transport residuals in the deep water.

Joyce, T. M.↗

High Resolution Doppler Lidar

This Grant supported the development of an incoherent lidar system to measure winds and aerosols in the lower atmosphere. During this period the following activities occurred: (1) an active feedback system was developed to improve the laser frequency stability; (2) a detailed forward model of the instrument was developed to take into account many subtle effects, such as detector non-linearity; (3) a non-linear least squares inversion method was developed to recover the Doppler shift and aerosol backscatter without requiring assumptions about the molecular component of the signal; (4) a study was done of the effects of systematic errors due to multiple etalon misalignment. It was discovered that even for small offsets and high aerosol loadings, the wind determination can be biased by as much as 1 m/s. The forward model and inversion process were modified to account for this effect; and (5) the lidar measurements were validated using rawinsonde balloon measurements. The measurements were found to be in agreement within 1-2 m/s.

Source record↗

A case study of the sensitivity of forecast skill to data and data analysis techniques

A series of experiments have been conducted to examine the sensitivity of forecast skill to various data and data analysis techniques for the 0000 GMT case of January 21, 1979. These include the individual components of the FGGE observing system, the temperatures obtained with different satellite retrieval methods, and the method of vertical interpolation between the mandatory pressure analysis levels and the model sigma levels. It is found that NESS TIROS-N infrared retrievals seriously degrade a rawinsonde-only analysis over land, resulting in a poorer forecast over North America. Less degradation in the 72-hr forecast skill at sea level and some improvement at 500 mb is noted, relative to the control with TIROS-N retrievals produced with a physical inversion method which utilizes a 6-hr forecast first guess. NESS VTPR oceanic retrievals lead to an improved forecast over North America when added to the control.

Baker, W. E.↗

Satellite Optical Remote Sensing of Clouds and Aerosols: From Particle Single-Scattering and Gaseous Absorption Through Radiative Transfer to Retrieval Products

Clouds and aerosols are fundamental regulators of Earth’s radiation budget and climate system, influencing both solar and terrestrial radiation through scattering, absorption, and emission processes. Accurate characterization of their physical and radiative properties from space requires a rigorous understanding of particle single-scattering, gaseous absorption, and radiative transfer in the atmosphere, as well as reliable inversion methods. This review synthesizes the physical foundations and algorithmic implementations of satellite-based passive optical remote sensing of clouds and aerosols, spanning the ultraviolet to thermal infrared spectral range. Beginning with electromagnetic scattering theory and state-of-the-art methods for computing single-scattering by nonspherical particles and computationally efficient methods for accounting for atmospheric absorption, we discuss the radiative transfer framework underpinning cloud and aerosol retrievals. The connection between single-scattering and multiple-scattering is rigorously formulated. We then summarize operational and research-grade retrieval techniques, including cloud masking and thermodynamic phase determination, CO₂ slicing for cloud-top pressure, the Nakajima-King shortwave bi-spectral, and infrared split-window approaches for cloud optical thickness and effective particle size, inversion algorithms for determining aerosol properties from multi-spectral and/or multi-angle radiometric and polarimetric measurements, and active-passive sensing synergy. Examples of the global cloud and aerosol climatologies are illustrated using observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Multi-angle Imaging SpectroRadiometer (MISR). Furthermore, the unique strengths of active remote sensing techniques based on spaceborne lidar observations are briefly elaborated in the context of studying ice clouds composed of randomly and horizontally oriented ice crystals, which is a significant challenge for conventional passive remote sensing techniques. By connecting physical theory to practical retrievals, this review highlights both the maturity of current methodologies and the remaining challenges in reducing uncertainties in particle morphology, vertical structure, absorption, and aerosol-cloud interactions. Furthermore, the impact of artificial intelligence (AI) on atmospheric remote sensing is briefly addressed.

Aerosols↗

Assessing the impact of different satellite retrieval methods on forecast available potential energy

The effects of the inclusion of satellite temperature retrieval data, and of different satellite retrieval methods, on forecasts made with the NASA Goddard Laboratory for Atmospheres (GLA) fourth-order model were investigated using, as the parameter, the available potential energy (APE) in its isentropic form. Calculation of the APE were used to study the differences in the forecast sets both globally and in the Northern Hemisphere during 72-h forecast period. The analysis data sets used for the forecasts included one containing the NESDIS TIROS-N retrievals, the GLA retrievals using the physical inversion method, and a third, which did not contain satellite data, used as a control; two data sets, with and without satellite data, were used for verification. For all three data sets, the Northern Hemisphere values for the total APE showed an increase throughout the forecast period, mostly due to an increase in the zonal component, in contrast to the verification sets, which showed a steady level of total APE.

Whittaker, Linda M.↗

Trailing edge flow conditions as a factor in airfoil design

Some new developments relevant to the design of single-element airfoils using potential flow methods are presented. In particular, the role played by the non-dimensional trailing edge velocity in design is considered and the relationship between the specified value and the resulting airfoil geometry is explored. In addition, the ramifications of the unbounded trailing edge pressure gradients generally present in the potential flow solution of the flow over an airfoil are examined, and the conditions necessary to obtain a class of airfoils having finite trailing edge pressure gradients developed. The incorporation of these conditions into the inverse method of Eppler is presented and the modified scheme employed to generate a number of airfoils for consideration. The detailed viscous analysis of airfoils having finite trailing edge pressure gradients demonstrates a reduction in the strong inviscid-viscid interactions generally present near the trailing edge of an airfoil.

Ormsbee, A. I.↗

Time-Distance Helioseismology with the MDI Instrument: Initial Results

In time-distance helioseismology, the travel time of acoustic waves is measured between various points on the solar surface. To some approximation, the waves can be considered to follow ray paths that depend only on a mean solar model, with the curvature of the ray paths being caused by the increasing sound speed with depth below the surface. The travel time is effected by various inhomogeneities along the ray path, including flows, temperature inhomogeneities, and magnetic fields. By measuring a large number of times between different locations and using an inversion method, it is possible to construct 3-dimensional maps of the subsurface inhomogeneities. The SOI/MDI experiment on SOHO has several unique capabilities for time-distance helioseismology. The great stability of the images observed without benefit of an intervening atmosphere is quite striking. It his made it possible for us to detect the travel time fo separations of points as small as 2.4 Mm in the high-resolution mode of MDI (0.6 arc sec 1/pixel). This has enabled the detection of the supergranulation flow. Coupled with the inversion technique, we can now study the 3-dimensional evolution of the flows near the solar surface.

Duvall, T. L., Jr.↗

Retrievals of Profiles of Fine And Coarse Aerosols Using Lidar And Radiometric Space Measurements

In couple of years we expect the launch of the CALIPSO lidar spaceborne mission designed to observe aerosols and clouds. CALIPSO will collect profiles of the lidar attenuated backscattering coefficients in two spectral wavelengths (0.53 and 1.06 microns). Observations are provided along the track of the satellite around the globe from pole to pole. The attenuated backscattering coefficients are sensitive to the vertical distribution of aerosol particles, their shape and size. However the information is insufficient to be mapped into unique aerosol physical properties and vertical distribution. Infinite number of physical solutions can reconstruct the same two wavelength backscattered profile measured from space. CALIPSO will fly in formation with the Aqua satellite and the MODIS spectro-radiometer on board. Spectral radiances measured by MODIS in six channels between 0.55 and 2.13 microns simultaneously with the CALIPSO observations can constrain the solutions and resolve this ambiguity, albeit under some assumptions. In this paper we describe the inversion method and apply it to aircraft lidar and MODIS data collected over a dust storm off the coast of West Africa during the SHADE experiment. It is shown that the product of the single scattering albedo, omega, and the phase function, P, for backscattering can be retrieved from the synergism between measurements avoiding a priori hypotheses required for inverting lidar measurements alone. The resultant value of (omega)P(180 deg.) = 0.016/sr are significantly different from what is expected using Mie theory, but are in good agreement with recent results obtained from lidar observations of dust episodes. The inversion is robust in the presence of noise of 10% and 20% in the lidar signal in the 0.53 and 1.06 pm channels respectively. Calibration errors of the lidar of 5 to 10% can cause an error in optical thickness of 20 to 40% respectively in the tested cases. The lidar calibration errors cause degradation in the ability to fit the MODIS data. Therefore the MODIS measurements can be used to identify the calibration problem and correct for it. The CALIPSO-MODIS measurements of the profiles of fine and coarse aerosols, together with CALIPSO measurements of clouds vertical distribution, is expected to be critically important in understanding aerosol transport across continents and political boundaries, and to study aerosol-cloud interaction and its effect on precipitation and global forcing of climate.

Kaufman, Yoram↗

Thermal Analysis Methods and Basic Heat-Transfer Data for a Turbulent Heating Flight Experiment at Mach 20 (Reentry F)

A heat-transfer experiment was flight conducted on a 5 deg half-angle cone, 396.2 cm (13 ft) in length, as it entered the sensible atmosphere under laminar, transitional, and turbulent boundary-layer conditions at a free-stream Mach number of about 20. Accurate turbulent-heat-transfer data with natural transition were obtained for correlation with theories in regions of simultaneous high Mach number, Reynolds number, enthalpy, and total-to-wall temperature ratio. Temperatures were measured at four depths through the 15.24-mm-thick (0.600-in.) beryllium wall. Experimental heating rates at 20 stations on the cone were determined independently from the outermost temperature measurement and from the temperature measurement at the second depth by a single-thermocouple inverse method and also from the temperature histories at all four depths by an integral method. The thermal data analysis procedure, associated problems, and results are presented herein.

Howard, Floyd G.↗

Mass, heat and freshwater fluxes in the South Indian Ocean

Six hydrographic sections were used to examine the circulation and property fluxes in the South Indian Ocean from 10 to 32 deg S. The calculations were made by applying an inverse method to the data. In the interior of the South Indian Ocean, the geostrophic flow is generally northward. At 18 deg S, the northward interior mass flux is balanced by the southward Ekman mass flux at the surface, whereas at 32 deg S the northward interior mass flux is balanced by the southward mass flux of the Agulhas Current. There is a weak, southward mass flux of 6 x 10 to the 9th kg/s in the Mozambique Channel. The rate of water exchange between the Pacific Ocean and the Indian Ocean is dependent on the choice of the initial reference level used in the inverse calculation. The choice of 1500 m, the depth of the deep oxygen minimum, has led to a flux of water from the Pacific Ocean to the Indian Ocean at a rate of 6.6 x 10 to the 9th kg/s. Heat flux calculations indicate that the Indian Ocean is exporting heat to the rest of the world's oceans at a rate of -0.69 x 10 to the 15th W at 18 deg S and -0.25 x 10 to the 15th W at 32 deg S (negative values being southward).

Fu, Lee-Lueng↗

Satellite inventory of Minnesota forest resources

The methods and results of using Landsat Thematic Mapper (TM) data to classify and estimate the acreage of forest covertypes in northeastern Minnesota are described. Portions of six TM scenes covering five counties with a total area of 14,679 square miles were classified into six forest and five nonforest classes. The approach involved the integration of cluster sampling, image processing, and estimation. Using cluster sampling, 343 plots, each 88 acres in size, were photo interpreted and field mapped as a source of reference data for classifier training and calibration of the TM data classifications. Classification accuracies of up to 75 percent were achieved; most misclassification was between similar or related classes. An inverse method of calibration, based on the error rates obtained from the classifications of the cluster plots, was used to adjust the classification class proportions for classification errors. The resulting area estimates for total forest land in the five-county area were within 3 percent of the estimate made independently by the USDA Forest Service. Area estimates for conifer and hardwood forest types were within 0.8 and 6.0 percent respectively, of the Forest Service estimates. A trial of a second method of estimating the same classes as the Forest Service resulted in standard errors of 0.002 to 0.015. A study of the use of multidate TM data for change detection showed that forest canopy depletion, canopy increment, and no change could be identified with greater than 90 percent accuracy. The project results have been the basis for the Minnesota Department of Natural Resources and the Forest Service to define and begin to implement an annual system of forest inventory which utilizes Landsat TM data to detect changes in forest cover.

Bauer, Marvin E.↗

Simulation of the Microwave Emission of Multi-layered Snowpacks Using the Dense Media Radiative Transfer Theory: the DMRT-ML Model

DMRT-ML is a physically based numerical model designed to compute the thermal microwave emission of a given snowpack. Its main application is the simulation of brightness temperatures at frequencies in the range 1-200 GHz similar to those acquired routinely by spacebased microwave radiometers. The model is based on the Dense Media Radiative Transfer (DMRT) theory for the computation of the snow scattering and extinction coefficients and on the Discrete Ordinate Method (DISORT) to numerically solve the radiative transfer equation. The snowpack is modeled as a stack of multiple horizontal snow layers and an optional underlying interface representing the soil or the bottom ice. The model handles both dry and wet snow conditions. Such a general design allows the model to account for a wide range of snow conditions. Hitherto, the model has been used to simulate the thermal emission of the deep firn on ice sheets, shallow snowpacks overlying soil in Arctic and Alpine regions, and overlying ice on the large icesheet margins and glaciers. DMRT-ML has thus been validated in three very different conditions: Antarctica, Barnes Ice Cap (Canada) and Canadian tundra. It has been recently used in conjunction with inverse methods to retrieve snow grain size from remote sensing data. The model is written in Fortran90 and available to the snow remote sensing community as an open-source software. A convenient user interface is provided in Python.

snowpacks↗