Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inverse 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 199 records · Page 11

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

Efficient multimode Wigner tomography

Abstract Advancements in quantum system lifetimes and control have enabled the creation of increasingly complex quantum states, such as those on multiple bosonic cavity modes. When characterizing these states, traditional tomography scales exponentially with the number of modes in both computational and experimental measurement requirement, which becomes prohibitive as the system size increases. Here, we implement a state reconstruction method whose sampling requirement instead scales polynomially with system size, and thus mode number, for states that can be represented within such a polynomial subspace. We demonstrate this improved scaling with Wigner tomography of multimode entangled W states of up to 4 modes on a 3D circuit quantum electrodynamics (cQED) system. This approach performs similarly in efficiency to existing matrix inversion methods for 2 modes, and demonstrates a noticeable improvement for 3 and 4 modes, with even greater theoretical gains at higher mode numbers.

Science & Technology - Other Topics↗

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↗

Novel estimation of aerosol processes with particle size distribution measurements: a case study with the TOMAS algorithm v1.0.0

Atmospheric aerosol microphysical processes are a significant source of uncertainty in predicting climate change. Specifically, aerosol nucleation, emissions, and growth rates, which are simulated in chemical transport models to predict the particle size distribution, are not understood well. However, long-term size distribution measurements made at several ground-based sites across Europe implicitly contain information about the processes that created those size distributions. This work aims to extract that information by developing and applying an inverse technique to constrain aerosol emissions as well as nucleation and growth rates based on hourly size distribution measurements. We developed an inverse method based upon process control theory into an online estimation technique to scale aerosol nucleation, emissions, and growth so that the model–measurement bias in three measured aerosol properties exponentially decays. The properties, which are calculated from the measured and predicted size distributions, used to constrain aerosol nucleation, emission, and growth rates are the number of particles with a diameter between 3 and 6 nm, the number with a diameter greater than 10 nm, and the total dry volume of aerosol ( N 3–6 , N 10 , V dry ), respectively. In this paper, we focus on developing and applying the estimation methodology in a zero-dimensional “box” model as a proof of concept before applying it to a three-dimensional simulation in subsequent work. The methodology is first tested on a dataset of synthetic and perfect measurements that span diverse environments in which the true particle emissions, growth, and nucleation rates are known. The inverse technique accurately estimates the aerosol microphysical process rates with an average and maximum error of 2 % and 13 %, respectively. Next, we investigate the effect that measurement noise has on the estimated rates. The method is robust to typical instrument noise in the aerosol properties as there is a negligible increase in the bias of the estimated process rates. Finally, the methodology is applied to long-term datasets of in situ size distribution measurements in western Europe from May 2006 through June 2007. At Melpitz, Germany, and Hyytiälä, Finland, the average diurnal profiles of estimated 3 nm particle formation rates are reasonable, having peaks near noon local time with average peak values of 1 and 0.15 cm –3 s –1 , respectively. The normalized absolute error in estimated N 3–6 , N 10 , and V dry at three European measurement sites is less than 15 %, showing that the estimation framework developed here has potential to decrease model–measurement bias while constraining uncertain aerosol microphysical processes.

54 ENVIRONMENTAL SCIENCES↗

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

Variational encoder geostatistical analysis (VEGAS) with an application to large scale riverine bathymetry

Estimation of riverbed profiles, also known as bathymetry, plays a vital role in many applications, such as safe and efficient inland navigation, prediction of bank erosion, land subsidence, and flood risk management. The high cost and complex logistics of direct bathymetry surveys, i.e, depth imaging, have encouraged the use of indirect measurements such as surface flow velocities. However, estimating high-resolution bathymetry from indirect measurements is an inverse problem that can be computationally challenging. Here, we propose a reduced-order model (ROM) based approach that utilizes a variational autoencoder (VAE), a type of deep neural network with a narrow layer in the middle, to compress bathymetry and flow velocity information and accelerate bathymetry inverse problems from flow velocity measurements. In our application, the shallow-water equations (SWE) with appropriate boundary conditions (BCs), e.g., the discharge and/or the free surface elevation, constitute the forward problem, to predict flow velocity. Then, ROMs of the SWEs are constructed on a nonlinear manifold of low dimensionality through a variational encoder and the bathymetry inversion problem is derived on the low-dimensional latent space in a Hierarchical Bayesian setting. Further, the reformulation allows variational inference with a small number (e.g., $\mathscr{O}$ (100) of ROM runs and efficient uncertainty quantification. We have tested our inversion approach on a one-mile reach of the Savannah River, GA, USA. Once the neural network is trained (offline stage), the proposed technique can perform the inversion operation orders of magnitude faster than traditional inversion methods that are commonly based on linear projections, such as principal component analysis (PCA), or the principal component geostatistical approach (PCGA). Furthermore, tests show that the algorithm can estimate the bathymetry with good accuracy even with sparse flow velocity measurements.

54 ENVIRONMENTAL SCIENCES↗

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↗

Machine learning inversion from scattering for mechanically driven polymers

A machine learning inversion method is developed for analyzing scattering functions of mechanically driven polymers and extracting the corresponding feature parameters, which include energy parameters and conformation variables. The polymer is modeled as a chain of fixed-length bonds constrained by bending energy, and it is subject to external forces such as stretching and shear. We generate a data set consisting of random combinations of energy parameters, including bending modulus, stretching and shear force, along with Monte Carlo-calculated scattering functions and conformation variables such as end-to-end distance, radius of gyration and off-diagonal component of the gyration tensor. The effects of the energy parameters on the polymer are captured by the scattering function, and principal component analysis ensures the feasibility of the machine learning inversion. Finally, we train a Gaussian process regressor using part of the data set as a training set and validate the trained regressor for inversion using the rest of the data. The regressor successfully extracts the feature parameters.

Gaussian process regressors↗

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↗

Antimony nanobelt asymmetric membranes for sodium ion battery

In this study, composite asymmetric membranes containing antimony (Sb) nanobelts are prepared via a straightforward phase inversion method in combination with post-pyrolysis treatment. Sb nanobelt asymmetric membranes demonstrate improved cyclability and specific capacity as the alloy anode of sodium ion battery compared to Sb nanobelt thin films without asymmetric porous structure. The unique structure can effectively accommodate the large volume expansion of Sb-based alloy anodes, prohibit the loss of fractured active materials, and aid in the formation of stable artificial solid electrolyte interphases as evidenced by an outstanding capacity retention of ~98% in 130 cycles at 60 mA g -1 . A specific capacity of ~600 mAh g -1 is obtained at 15 mA g -1 (1/40C). When the current density is increased to 240 mA g -1 , ~80% capacity can be maintained (~480 mAh g -1 ). The relations among phase inversion conditions, structures, compositions, and resultant electrochemical properties are revealed through comprehensive characterization.

25 ENERGY STORAGE↗

Microseismic Monitoring at the Farnsworth CO 2 -EOR Field

The Farnsworth Unit in northern Texas is a field site for studying geologic carbon storage during enhanced oil recovery (EOR) using CO 2 . Microseismic monitoring is essential for risk assessment by detecting fluid leakage and fractures. We analyzed borehole microseismic data acquired during CO 2 injection and migration, including data denoising, event detection, event location, magnitude estimation, moment tensor inversion, and stress field inversion. We detected and located two shallow clusters, which occurred during increasing injection pressure. The two shallow clusters were also featured by large b values and tensile cracking moment tensors that are obtained based on a newly developed moment tensor inversion method using single-borehole data. The inverted stress fields at the two clusters showed large deviations from the regional stress field. The results provide evidence for microseismic responses to CO 2 /fluid injection and migration.

02 PETROLEUM↗