Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model inversion”

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 919 records · Page 51

Radiative-hydrodynamical simulations of accretion disk coronae

An algorithm has been developed that will carry out detailed, two-dimensional, fully self-consistent radiative-hydrodynamical simulations of accretion disk coronae in X-ray-luminous compact sources. The calculation reported here, for an accreting neutron star radiating at 0.5 times the Eddington luminosity, reveals several striking features. (1) The corona is comprised of two main regions - an inner highly dynamic portion whose vertical structure varies cyclically on a dynamical time scale, and an outer, more stable zone in which the evaporated plasma rises to form a 'sheath' that gradually merges into a wind at large radii. (2) The coronal structure shows a density inversion that contrasts sharply with the Gaussian profiles of earlier hydrostatic models. (3) Interestingly, flow velocities as high as a few billion cm/s are not uncommon in portions of the corona. The relevance of this point to the large measured width of the emission features seen in low-mass X-ray binaries and the Galactic black hole candidates is discussed.

Melia, Fulvio↗

Reliable automatic plan selection for visual robotic positioning

Reliability analysis techniques developed for robotics are applied to a six-DOF (degree-of-freedom) visual positioning system. The task is to measure an object's pose and then grip the object reliably subject to time and accuracy constraints. Errors due to pixel truncation in the stadimetric vision system, forward kinematic and inverse kinematic uncertainties, and joint position and velocity noise are stochastically modeled. Simulation results for PUMA 500 kinematic and dynamic characteristics are presented.

Mcinroy, John E.↗

Ozone vertical distribution retrieval from ground-based high resolution infrared solar spectra

A practical procedure for the retrieval of ozone vertical profiles from ground-based high resolution Fourier transform infrared solar spectra has been developed. The analysis is based on a multilayer line-by-line forward model and a semi-empirical version of the optimal estimation inversion method of Rodgers. The 1002.6-1003.2 cm(exp -1) spectral interval has been selected for the analysis on the basis of synthetic spectrum calculations. This interval contains numerous ozone lines covering a range of intensities and providing retrieval sensitivity from ground level to about 35 km. Characterization of the method and an error analysis have been performed. For a spectral resolution of 0.05-0.01 cm(exp -1) and a signal-to-noise ratio greater than or equal to 100 the retrieval is stable with a vertical resolution of approximately 5 km attainable near the surface degrading to approximately 10 km in the stratosphere. Synthetic spectra studies show that the a priori profile and weak constraints selected for the retrievals introduce no significant biases for a wide range of ozone profiles.

Pougatchev, N. S.↗

Global Analysis, Interpretation, and Modelling: First Science Conference

Topics considered include: Biomass of termites and their emissions of methane and carbon dioxide - A global database; Carbon isotope discrimination during photosynthesis and the isotope ratio of respired CO2 in boreal forest ecosystems; Estimation of methane emission from rice paddies in mainland China; Climate and nitrogen controls on the geography and timescales of terrestrial biogeochemical cycling; Potential role of vegetation feedback in the climate sensitivity of high-latitude regions - A case study at 6000 years B.P.; Interannual variation of carbon exchange fluxes in terrestrial ecosystems; and Variations in modeled atmospheric transport of carbon dioxide and the consequences for CO2 inversions.

Sahagian, Dork↗

Spectroscopy of Solid State Laser Materials

We retrieved the vertical distribution of ozone from a series 0.005-0.013/cm resolution infrared solar spectra recorded with the McMath Fourier Transform spectrometer at the Kitt Peak National Solar Observatory. The analysis is based on a multi-layer line-by-line forward model and a semi-empirical version of the optimal estimation inversion method by Rodgers. The 1002.6-1003.2/cm spectral interval has been selected for the analysis on the basis of synthetic spectrum calculations. The characterization and error analysis of the method have been performed. It was shown that for the Kitt Peak spectral resolution and typical signal-to-noise ratio (greater than or equal to 100) the retrieval is stable, with the vertical resolution of approximately 5 km attainable near the surface degrading to approximately 10 km in the stratosphere. Spectra recorded from 1980 through 1993 have been analyzed. The retrieved total ozone and vertical profiles have been compared with total ozone mapping spectrometer (TOMS) satellite total columns for the location and dates of the Kitt Peak Measurements and about 100 ozone ozonesoundings and Brewer total column measurements from Palestine, Texas, from 1979 to 1985. The total ozone measurements agree to +/- 2%. The retrieved profiles reproduce the seasonally averaged variations with altitude, including the ozone spring maximum and fall minimum measured by Palestine sondes, but up to 15% differences in the absolute values are obtained.

Buoncristiani, A. M.↗

The 2010 May Flaring Episode of Cygnus X-3 in Radio, X-Rays, and gamma-Rays

In 2009, Cygnus X-3 (Cyg X-3) became the first microquasar to be detected in the GeV gamma-ray regime, via the satellites Fermi and AGILE. The addition of this new band to the observational toolbox holds promise for building a more detailed understanding of the relativistic jets of this and other systems. We present a rich dataset of radio, hard and soft X-ray, and gamma-ray observations of Cyg X-3 made during a flaring episode in 2010 May. We detect a approx.3-d softening and recovery of the X-ray emission, followed almost immediately by a approx.1-Jy radio flare at 15 GHz, followed by a 4.3sigma gamma-ray flare (E > 100 MeV) approx.1.5 d later. The radio sampling is sparse, but we use archival data to argue that it is unlikely the gamma-ray flare was followed by any significant unobserved radio flares. In this case, the sequencing of the observed events is difficult to explain in a model in which the gamma-ray emission is due to inverse Compton scattering of the companion star's radiation field. Our observations suggest that other mechanisms may also be responsible for gamma-ray emission from Cyg X-3.

Williams, Peter K. G.↗

Multi-Angle Implementation of Atmospheric Correction for MODIS (MAIAC). Part 3: Atmospheric Correction

This paper describes the atmospheric correction (AC) component of the Multi-Angle Implementation of Atmospheric Correction algorithm (MAIAC) which introduces a new way to compute parameters of the Ross-Thick Li-Sparse (RTLS) Bi-directional reflectance distribution function (BRDF), spectral surface albedo and bidirectional reflectance factors (BRF) from satellite measurements obtained by the Moderate Resolution Imaging Spectroradiometer (MODIS). MAIAC uses a time series and spatial analysis for cloud detection, aerosol retrievals and atmospheric correction. It implements a moving window of up to 16 days of MODIS data gridded to 1 km resolution in a selected projection. The RTLS parameters are computed directly by fitting the cloud-free MODIS top of atmosphere (TOA) reflectance data stored in the processing queue. The RTLS retrieval is applied when the land surface is stable or changes slowly. In case of rapid or large magnitude change (as for instance caused by disturbance), MAIAC follows the MODIS operational BRDF/albedo algorithm and uses a scaling approach where the BRDF shape is assumed stable but its magnitude is adjusted based on the latest single measurement. To assess the stability of the surface, MAIAC features a change detection algorithm which analyzes relative change of reflectance in the Red and NIR bands during the accumulation period. To adjust for the reflectance variability with the sun-observer geometry and allow comparison among different days (view geometries), the BRFs are normalized to the fixed view geometry using the RTLS model. An empirical analysis of MODIS data suggests that the RTLS inversion remains robust when the relative change of geometry-normalized reflectance stays below 15%. This first of two papers introduces the algorithm, a second, companion paper illustrates its potential by analyzing MODIS data over a tropical rainforest and assessing errors and uncertainties of MAIAC compared to conventional MODIS products.

Lyapustin, A.↗

Two Distinct Roles of Atlantic SSTs in ENSO Variability: North Tropical Atlantic SST and Atlantic Nino

Two distinct roles of the Atlantic sea surface temperatures (SSTs), namely, the North Tropical Atlantic (NTA) SST and the Atlantic Nino, on the El Nino-Southern Oscillation (ENSO) variability are investigated using the observational data from 1980 to 2010 and coupled model experiments. It appears that the NTA SST and the Atlantic Nino can be used as two independent predictors for predicting the development of ENSO events in the following season. Furthermore, they are likely to be linked to different types of El Nino events. Specifically, the NTA SST cooling during February, March, and April contributes to the central Pacific warming at the subsequent winter season, while the negative Atlantic Nino event during June, July, and August contributes to enhancing the eastern Pacific warming. The coupled model experiments support these results. With the aid of a lagged inverse relationship, the statistical forecast using two Atlantic indices can successfully predict various ENSO indices.

Atlantic Nino↗

Carbon Monitoring System Flux Estimation and Attribution: Impact of ACOS-GOSAT X(CO2) Sampling on the Inference of Terrestrial Biospheric Sources and Sinks

Using an Observing System Simulation Experiment (OSSE), we investigate the impact of JAXA Greenhouse gases Observing SATellite 'IBUKI' (GOSAT) sampling on the estimation of terrestrial biospheric flux with the NASA Carbon Monitoring System Flux (CMS-Flux) estimation and attribution strategy. The simulated observations in the OSSE use the actual column carbon dioxide (X(CO2)) b2.9 retrieval sensitivity and quality control for the year 2010 processed through the Atmospheric CO2 Observations from Space algorithm. CMS-Flux is a variational inversion system that uses the GEOS-Chem forward and adjoint model forced by a suite of observationally constrained fluxes from ocean, land and anthropogenic models. We investigate the impact of GOSAT sampling on flux estimation in two aspects: 1) random error uncertainty reduction and 2) the global and regional bias in posterior flux resulted from the spatiotemporally biased GOSAT sampling. Based on Monte Carlo calculations, we find that global average flux uncertainty reduction ranges from 25% in September to 60% in July. When aggregated to the 11 land regions designated by the phase 3 of the Atmospheric Tracer Transport Model Intercomparison Project, the annual mean uncertainty reduction ranges from 10% over North American boreal to 38% over South American temperate, which is driven by observational coverage and the magnitude of prior flux uncertainty. The uncertainty reduction over the South American tropical region is 30%, even with sparse observation coverage. We show that this reduction results from the large prior flux uncertainty and the impact of non-local observations. Given the assumed prior error statistics, the degree of freedom for signal is approx.1132 for 1-yr of the 74 055 GOSAT X(CO2) observations, which indicates that GOSAT provides approx.1132 independent pieces of information about surface fluxes. We quantify the impact of GOSAT's spatiotemporally sampling on the posterior flux, and find that a 0.7 gigatons of carbon bias in the global annual posterior flux resulted from the seasonally and diurnally biased sampling when using a diagonal prior flux error covariance.

biased sampling↗

Retrieval of Aerosol Microphysical Properties from AERONET Photo-Polarimetric Measurements. 2: A New Research Algorithm and Case Demonstration

A new research algorithm is presented here as the second part of a two-part study to retrieve aerosol microphysical properties from the multispectral and multiangular photopolarimetric measurements taken by Aerosol Robotic Network's (AERONET's) new-generation Sun photometer. The algorithm uses an advanced UNified and Linearized Vector Radiative Transfer Model and incorporates a statistical optimization approach.While the new algorithmhas heritage from AERONET operational inversion algorithm in constraining a priori and retrieval smoothness, it has two new features. First, the new algorithmretrieves the effective radius, effective variance, and total volume of aerosols associated with a continuous bimodal particle size distribution (PSD) function, while the AERONET operational algorithm retrieves aerosol volume over 22 size bins. Second, our algorithm retrieves complex refractive indices for both fine and coarsemodes,while the AERONET operational algorithm assumes a size-independent aerosol refractive index. Mode-resolved refractive indices can improve the estimate of the single-scattering albedo (SSA) for each aerosol mode and thus facilitate the validation of satellite products and chemistry transport models. We applied the algorithm to a suite of real cases over Beijing_RADI site and found that our retrievals are overall consistent with AERONET operational inversions but can offer mode-resolved refractive index and SSA with acceptable accuracy for the aerosol composed by spherical particles. Along with the retrieval using both radiance and polarization, we also performed radiance-only retrieval to demonstrate the improvements by adding polarization in the inversion. Contrast analysis indicates that with polarization, retrieval error can be reduced by over 50% in PSD parameters, 10-30% in the refractive index, and 10-40% in SSA, which is consistent with theoretical analysis presented in the companion paper of this two-part study.

aerosol retrieval algorithm↗

Assessing Uncertainties in Gridded Emissions: A Case Study for Fossil Fuel Carbon Dioxide (FFCO2) Emission Data

Fossil fuel carbon dioxide (CO2) emissions (FFCO2) are the largest input to the global carbon cycle on a decadal time scale. Because total emissions are assumed to be reasonably well constrained by fuel statistics, FFCO2 often serves as a reference in order to deduce carbon uptake by poorly understood terrestrial and ocean sinks. Conventional atmospheric CO2 flux inversions solve for spatially explicit regional sources and sinks and estimate land and ocean fluxes by subtracting FFCO2. Thus, errors in FFCO2 can propagate into the final inferred flux estimates. Gridded emissions are often based on disaggregation of emissions estimated at national or regional level. Although national and regional total FFCO2 are well known, gridded emission fields are subject to additional uncertainties due to the emission disaggregation. Assessing such uncertainties is often challenging because of the lack of physical measurements for evaluation. We first review difficulties in assessing uncertainties associated with gridded FFCO2 emission data and present several approaches for evaluation of such uncertainties at multiple scales. Given known limitations, inter-emission data differences are often used as a proxy for the uncertainty. The popular approach allows us to characterize differences in emissions, but does not allow us to fully quantify emission disaggregation biases. Our work aims to vicariously evaluate FFCO2 emission data using atmospheric models and measurements. We show a global simulation experiment where uncertainty estimates are propagated as an atmospheric tracer (uncertainty tracer) alongside CO2 in NASA's GEOS model and discuss implications of FFCO2 uncertainties in the context of flux inversions. We also demonstrate the use of high resolution urban CO2 simulations as a tool for objectively evaluating FFCO2 data over intense emission regions. Though this study focuses on FFCO2 emission data, the outcome of this study could also help improve the knowledge of similar gridded emissions data for non-CO2 compounds with similar emission characteristics.

FFCO↗

Constraining Venus and Mars Internal Structures from the Gravity Response to Atmospheric Loading

The gravity fields of celestial bodies that possess an atmosphere are periodically perturbed by the redistribution of fluid mass associated with the atmospheric dynamics. A component of this perturbation is due to the gravitational response of the body to the deformation of its surface induced by the atmospheric pressure loading. The magnitude of this effect depends on the relation between the loading and the response in terms of geopotential variations measured by the load Love numbers. In this work, we simulate and analyze the gravity field generated by the atmospheres of Venus and Mars by accounting for different models of their internal structure. By precisely characterizing the phenomena that drive the mass transportation in the atmosphere through general circulation models, we determine the effect of the interior structure on the response to the atmospheric loading. An accurate estimation of the time-varying gravity field, which measures the atmospheric contribution, may provide significant constraints on the interior structure through the measurement of the load Love numbers. A combined determination of tidal and load Love numbers would enhance our knowledge of the interior of planetary bodies, providing further geophysical constraints in the inversion of internal structure models.The gravity fields of celestial bodies that possess an atmosphere are periodically perturbed by the redistribution of fluid mass associated with the atmospheric dynamics. A component of this perturbation is due to the gravitational response of the body to the deformation of its surface induced by the atmospheric pressure loading. The magnitude of this effect depends on the relation between the loading and the response in terms of geopotential variations measured by the load Love numbers. In this work, we simulate and analyze the gravity field generated by the atmospheres of Venus and Mars by accounting for different models of their internal structure. By precisely characterizing the phenomena that drive the mass transportation in the atmosphere through general circulation models, we determine the effect of the interior structure on the response to the atmospheric loading. An accurate estimation of the time-varying gravity field, which measures the atmospheric contribution, may provide significant constraints on the interior structure through the measurement of the load Love numbers. A combined determination of tidal and load Love numbers would enhance our knowledge of the interior of planetary bodies, providing further geophysical constraints in the inversion of internal structure models.

Mars↗

July 16, 1987 revisited: Lessons for modelers

A cross pattern was flown by the NCAR Electra in the restricted air space around San Nicolas Island (SNI). There was available a LANDSAT image coincident with the mission. The variability of cloud thickness along the flight track is discussed. The evidence points to variability in the water vapor content above the inversion as the controlling factor. This compounds the difficulty of parameterizing these clouds in General Circulation Models (GCMs). The importance of specifying correctly the upper boundary condition in cloud models is stressed. Particularly, in GCMs, calculating the humidity in the lower layers of the model is crucial for successful parameterization of marine stratocumulus clouds.

Hanson, Howard P.↗

INVESTIGATING THE FEASIBILITY AND STABILITY FOR MODELING ACOUSTIC WAVE SCATTERING USING A TIME-DOMAIN BOUNDARY INTEGRAL EQUATION WITH IMPEDANCE BOUNDARY CONDITION

Reducing aircraft noise is a major objective in the field of computational aeroacoustics. When designing next generation quiet and environmentally friendly aircraft, it is important to be able to accurately and efficiently predict the acoustic scattering by an aircraft body from a given noise source. Acoustic liners are an effective tool for aircraft noise reduction and are characterized by a frequency-dependent impedance. Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic wave scattering by geometric bodies treated with acoustic liners. This work considers using either an impedance or an admittance (inverse of impedance) boundary condition to allow for acoustic scattering problems to be modeled with geometries consisting of both un-lined and lined surfaces. Three acoustic liner models are discussed: the Extended Helmholtz Resonator Model, the Three-Parameter Impedance Model, and the Broadband Impedance Model. In both the Helmholtz and Three-Parameter models, liner impedance is specified at a given frequency, whereas the Broadband model allows for the investigation of multiple frequencies simultaneously. The impedance and admittance boundary conditions for acoustic liners are derived for each model and coupled with a time-domain boundary integral equation. The scattering solution is obtained iteratively using a boundary element method with constant spatial and third-order temporal basis functions. Time-domain boundary integral equations are unfortunately prone to numerical instabilities due to resonant frequencies resulting from non-trivial solutions in the interior domain. When reformulated with the Burton-Miller method, the instabilities are eliminated. Using a Burton-Miller reformulation, the stability of the boundary element method assuming a liner boundary condition is assessed using eigenvalue analysis. The stability of each liner model is discussed, and it is shown that the Three-Parameter and Broadband models are sufficient for modeling an acoustic liner on the surface of scattering bodies. The Helmholtz model demonstrates strict limitations for stability, whereas the Three-Parameter and Broadband models are stable for most cases. Also included in this work is an assessment of the spatial accuracy of the time-domain boundary element method with respect to the surface element basis functions, as well as a performance study of the numerical algorithm.

computational mathematics↗

Ocean circulation modeling by use of radar altimeter data

The project will investigate the use of radar altimetry (RA) data in the determination of the ocean circulation models. RA data will be used to verify prognostic experiments of the steady state and seasonal cycle of large-scale circulation models and the statistical steady state of eddy-resolving models. The data will serve as initial and update conditions in data assimilation experiments and as constraints in inverse calculations. The aim of the project is a better understanding of ocean physics, the determination and mapping of ocean currents, and a contribution to the establishment of ocean circulation models for climate studies. The goal of the project is to use satellite radar altimetry data for improving our knowledge of ocean circulation both in a descriptive sense and through the physics that govern the circulation state. The basic tool is a series of ocean circulation models. Depending on the model, different techniques will be applied to incorporate the RA data.

Olbers, Dirk↗

MoSST DAS: The First Working Geomagnetic Data Assimilation System

The Earth possesses an internal magnetic field (geomagnetic field) generated by convection in the outer core (geodynamo). Previous efforts have been focused along two distinct paths: (1) numerical geodynamo modeling to understand the origin of the geomagnetic field, and the mechanisms of geomagnetic secular variations (SV); and (2) geomagnetic field modeling to map the spatial/temporal variations of the field from geomagnetic data, and to derive core properties, e.g. inversion of core flow near the core-mantle boundary (CMB). Geomagnetic data assimilation is a new approach emerged over the past 5 years: surface observations are assimilated with geodynamo models for better understanding of the core dynamical state, and accurately prediction of SV. In collaboration with several geomagnetic research groups, we have developed the first working geomagnetic data assimilation system, Modular, Scalable, Self-consistent, and Three-dimensional (MoSST) DAS, that includes the MoSST numerical dynamo model; 7000 years of geomagnetic field maps from several field models utilizing satellite and ground observatory data, historical magnetic records and archeo/paleo magnetic data; and an ensemble based optimal interpolation (01) assimilation algorithm. With this system, we have demonstrated clearly that the assimilated core dynamical state is substantially different from those of pure geodynamo simulations. Ensemble assimilation runs also show the convergence of the assimilated solutions inside the core, suggesting that the simulation state is pulled closer to the truth via data assimilation. The forecasts from this system are also very accurate: the 5-year forecast of the geomagnetic field agrees very well with the observations; and the 5-year secular variation forecast is more accurate than the IGRF SV forecast models in the past. Using geomagnetic records up to 2009, we have made an SV forecast for the period from 2010-2015, and is a candidate SV model for IGRF-11.

Kuang, Weijia↗

Synchrotron and inverse Compton emission from expanding sources in jets - Application to SS 433

The formalism for an expanding source is extended to include continuous injection of electrons and the production of X-rays by inverse Compton scattering of external photons and synchrotron-self-Compton scattering. The model is applied to the expanding nonthermal sources entrained in the jets of SS 433. These sources are found capable of producing both the flaring radio component and the X-ray fluxes only if the relativistic electrons and the magnetic field are far out of equipartition. If the expanding source in the jets is in equipartition, the X-rays may originate from a stationary source within the binary system, which is self-absorbed at radio frequencies and produces X-rays by inverse Compton scattering of optical photons from the binary.

Band, David L.↗

Multivariate space - time analysis of PRE-STORM precipitation

This paper presents the methodologies and results of the multivariate modeling and two-dimensional spectral and correlation analysis of PRE-STORM rainfall gauge data. Estimated parameters of the models for the specific spatial averages clearly indicate the eastward and southeastward wave propagation of rainfall fluctuations. A relationship between the coefficients of the diffusion equation and the parameters of the stochastic model of rainfall fluctuations is derived that leads directly to the exclusive use of rainfall data to estimate advection speed (about 12 m/s) as well as other coefficients of the diffusion equation of the corresponding fields. The statistical methodology developed here can be used for confirmation of physical models by comparison of the corresponding second-moment statistics of the observed and simulated data, for generating multiple samples of any size, for solving the inverse problem of the hydrodynamic equations, and for application in some other areas of meteorological and climatological data analysis and modeling.

Polyak, Ilya↗