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 505 records · Page 28

Model Evaluation of Methods for Estimating Surface Emissions and Chemical Lifetimes from Satellite Data

Column densities from satellite retrievals can provide valuable information for estimating emissions and chemical lifetimes objectively across the globe. To better understand the uncertainties associated with these estimates, we test four methods using simulated column densities from a point source: a box model approach, a 2D Gaussian fit, an Inverse Radius fit and an Exponentially-Modified Gaussian fit. The model results were simulated using the WRF and CAMx models for the year 2005, for a single point source outside Atlanta in Georgia, USA with specified emissions and three chemical scenarios: no chemical reactions, 12 h chemical lifetime and 1 h chemical lifetime. No other sources were included in the simulations. We find that the box model provides reliable estimates irrespective of plume speed and plume direction, if the plume speed and the chemical lifetime are known accurately. The 2D Gaussian fit was found to be sensitive to plume speed and direction, and requires omnidirectional dispersion in order to have a decent fit. However, the 2D Gaussian fit is only an approximate fit to the data, and the discrepancies mean that the results are dependent on the geographical domain used for the optimization. An Inverse Radius fit is introduced to correct this issue, which is found to provide improved emissions and lifetime estimates. The Exponentially-Modified Gaussian fit also gave improved estimates. It is however dependent on accurate plume rotation such that reported chemical lifetimes with this method could be significantly underestimated.

Emission inventory↗

Modeling Atmospheric CO2 Processes to Constrain the Missing Sink

We report on a NASA supported modeling effort to reduce uncertainty in carbon cycle processes that create the so-called missing sink of atmospheric CO2. Our overall objective is to improve characterization of CO2 source/sink processes globally with improved formulations for atmospheric transport, terrestrial uptake and release, biomass and fossil fuel burning, and observational data analysis. The motivation for this study follows from the perspective that progress in determining CO2 sources and sinks beyond the current state of the art will rely on utilization of more extensive and intensive CO2 and related observations including those from satellite remote sensing. The major components of this effort are: 1) Continued development of the chemistry and transport model using analyzed meteorological fields from the Goddard Global Modeling and Assimilation Office, with comparison to real time data in both forward and inverse modes; 2) An advanced biosphere model, constrained by remote sensing data, coupled to the global transport model to produce distributions of CO2 fluxes and concentrations that are consistent with actual meteorological variability; 3) Improved remote sensing estimates for biomass burning emission fluxes to better characterize interannual variability in the atmospheric CO2 budget and to better constrain the land use change source; 4) Evaluating the impact of temporally resolved fossil fuel emission distributions on atmospheric CO2 gradients and variability. 5) Testing the impact of existing and planned remote sensing data sources (e.g., AIRS, MODIS, OCO) on inference of CO2 sources and sinks, and use the model to help establish measurement requirements for future remote sensing instruments. The results will help to prepare for the use of OCO and other satellite data in a multi-disciplinary carbon data assimilation system for analysis and prediction of carbon cycle changes and carbodclimate interactions.

Kawa, S. R.↗

The Relaxation Matrix for Symmetric Tops with Inversion Symmetry: Line Mixing Effects in the V1 Band of NH3 - II

Line mixing effects have been calculated in the ν1 parallel band of self-broadened NH3. The theoretical approach is an extension of a semi-classical model to symmetric-top molecules with inversion symmetry developed in the companion paper [Q. Ma and C. Boulet, J. Chem. Phys. 144, 224303 (2016)]. This model takes into account line coupling effects and hence enables the calculation of the entire relaxation matrix. A detailed analysis of the various coupling mechanisms is carried out for Q and R inversion doublets. The model has been applied to the calculation of the shape of the Q branch and of some R manifolds for which an obvious signature of line mixing effects has been experimentally demonstrated. Comparisons with measurements show that the present formalism leads to an accurate prediction of the available experimental line shapes. Discrepancies between the experimental and theoretical sets of first order mixing parameters are discussed as well as some extensions of both theory and experiment.

Boulet, C.↗

Study of a Simulation Tool to Determine Achievable Control Dynamics and Control Power Requirements with Perfect Tracking

This paper contains a study of two methods for use in a generic nonlinear simulation tool that could be used to determine achievable control dynamics and control power requirements while performing perfect tracking maneuvers over the entire flight envelope. The two methods are NDI (nonlinear dynamic inversion) and the SOFFT(Stochastic Optimal Feedforward and Feedback Technology) feedforward control structure. Equivalent discrete and continuous SOFFT feedforward controllers have been developed. These equivalent forms clearly show that the closed-loop plant model loop is a plant inversion and is the same as the NDI formulation. The main difference is that the NDI formulation has a closed-loop controller structure whereas SOFFT uses an open-loop command model. Continuous, discrete, and hybrid controller structures have been developed and integrated into the formulation. Linear simulation results show that seven different configurations all give essentially the same response, with the NDI hybrid being slightly different. The SOFFT controller gave better tracking performance compared to the NDI controller when a nonlinear saturation element was added. Future plans include evaluation using a nonlinear simulation.

Ostroff, Aaron J.↗

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

Predictions of transonic separated flow with an eddy-viscosity/Reynolds-shear-stress closure model

In previous work, a new turbulence-closure model was specifically developed for two-dimensional turbulent boundary layers subjected to strong, adverse-pressure gradients and the attendant separation. This closure model was shown to perform well, but the inverse, boundary-layer calculations used in evaluating the model were not truly predictive since the mass-flux parameter distribution obtained from the experimental data had to be specified. The objective of this study was to evaluate this closure model for a series of strong, transonic, inviscid-viscous interactions with varying degrees of separation using a fully predictive calculation method based on the Reynolds-averaged, Navier-Stokes equations. Calculated results are compared with experimental results and are shown to be in excellent agreement, even for the interactions with massive separation. This closure model has the very favorable property of requiring little more computational effort than equilibrium, algebraic models.

Johnson, D. A.↗

An equivalent layer magnetization model for Australia based on Magsat data

An equivalent layer magnetization model for Australia and adjacent oceanic areas is presented. The model is obtained by linear inversion of Magsat anomaly data measured in the altitude range 325-550 km. The anomaly data set has been isolated from the raw data set by use of models of the core field and very long wavelength external fields, and is internally consistent. Certain major structural features of the Australian continent are geographically associated with magnetization anomalies. A first-order difference is seen between the Tasman Zone and the Precambrian cratonic areas: magnetization anomalies are much more subdued in the former, possibly reflecting a shallowing of the Curie isotherm within the crust. A profile of the vertical integral of magnetization is presented for a crustal section extending from the Gawler Block to the southeast coast. It is shown that the magnetization variations are probably due partly, but not wholly, to depth to Curie isotherm variations; gross magnetization variations among at least three distinct crustal units must be involved.

Mayhew, M. A.↗

Cosmic-ray-modified supernova remnant shocks

The evolution of cosmic-ray-modified SNRs during the early Sedov phase is investigated using the time-dependent, two-fluid model for diffusive shock acceleration. Consideration is given to the sensitivity of net acceleration efficiency to model assumptions regarding density structure of the external medium as well as time dependence in the diffusion coefficient and the cosmic-ray adiabatic index. Of the model assumptions explored, the greatest sensitivity was found to time variations in the specific heat ratio to cosmic rays. It was found that if tau is greater than or approximately equal to few times 100, dynamically significant cosmic-ray pressures are produced early in the Sedov inefficient cosmic-ray accelerators. For the SNR models under consideration, the total energy channeled into the cosmic rays can be of order 10 percent of the initial blast energy for both the uniform density and inverse-square density interstellar medium models.

Jones, T. W.↗

Bayesian inference in geomagnetism

The inverse problem in empirical geomagnetic modeling is investigated, with critical examination of recently published studies. Particular attention is given to the use of Bayesian inference (BI) to select the damping parameter lambda in the uniqueness portion of the inverse problem. The mathematical bases of BI and stochastic inversion are explored, with consideration of bound-softening problems and resolution in linear Gaussian BI. The problem of estimating the radial magnetic field B(r) at the earth core-mantle boundary from surface and satellite measurements is then analyzed in detail, with specific attention to the selection of lambda in the studies of Gubbins (1983) and Gubbins and Bloxham (1985). It is argued that the selection method is inappropriate and leads to lambda values much larger than those that would result if a reasonable bound on the heat flow at the CMB were assumed.

Backus, George E.↗

Landscape Evolution, A Comparison of Form and Process

This project's goals were to collect, analyze and interpret 3-dimensional physiographic data for understanding the processes responsible for landscape modification. The primary landforms to be studied were Neogene cinder cones in Arizona (San Francisco Volcanic Field (SFVF), Coconino and Kaibab National Forests, Arizona). We also obtained and are still analyzing digital topographic data for the Long Valley-White Mountains area of California, which display Quaternary normal fault scarps, as well as extensive evidence of degradation. The work resulted in a large database of measured rates of downslope transport of slope debris. It was hypothesized that the work would increase our understanding of process-response models of hillslope degradation, and of the effects of climate change and other parameters on degradation rates. In greater detail, our primary goal was to compare evolutionary sequences of hillslopes, as exemplified by the topography of landforms of the same type but of different ages, with measurements of the surficial processes active on the landforms. Assuming that other parameters, such as hillslope materials and vegetation are held constant, and that the effects of changing climate are negligible, then the sediment transport rates measured today on the landforms should be the same as those calculated from the inversion of landform topography by use of a diffusion-type model. However, if the effects of changing climate or other factors are not negligible, then the observed transport rates would differ from those which must be invoked to explain the current topography. We hypothesized in fact that because degradation on the event scale is highly transient and localized, we would find a wide divergence between modern, measured transport rates, and rates calculated by global landform inversion or modeling. Because of the length of time involved in collection of sufficient data on current degradation rates, we are still continuing to analyse and interpret the data. Completion of the work will increase our understanding of the potential effects of anthropogenic climate and surficial change on the Earth's solid surface, and possibly allow us to constrain paths of hill-slope evolution following anthropogenic modifications, as well as compare the short-term with the long-term rates of hillslope degradation.

Bursik, Marcus I.↗

(abstract) Asteroid 4769 Casalia: Intrepretation of Optical Lightcurves Using a Radar-Derived Shape Model

A 167-parameter, 3-D shape model of Castalia, obtained from inversion of delay-Doppler images, constrained the object's subradar latitude on the date of the observations by defining an annular sky region that contains Castalia's pole. We have used the shape model and nominal Hapke parameters for S asteroids to model Castalia lightcurves obtained at Table Mountain Observatory on August 23-25, 1989, at solar phase angles ranging from 58(deg) to 90(deg). This investigation defines the projected pole's azimuthal orientation around the radar line of sight and therefore yields a least-squares estimate of the orientation of Castalia's north or south pole; there is insufficient geometric leverage to distinguish the rotation sense or to estimate values of Hapke parameters. The model lightcurves match the TMO data surprisingly well, given the modest number of shape parameters in the model and the fact that part of the asteroid was not seen by the radar.

asteroid Castalia shape radar 3-D model optical li↗

Stellar occultations by turbulent planetary atmospheres. I - A heuristic scattering model. II - The Beta Scorpii events

Effects of atmospheric turbulence on stellar-occultation inversion procedures are investigated using a heuristic scattering model that is believed to reproduce the essential features of turbulence. A quantitative estimate is made of the size of the error in deducing the mean refractivity profile of a planetary atmosphere, taking into account constant as well as exponential scattering. It is shown that ordinary turbulence has no important effect on the average intensity profile in a stellar occultation but could have an important instantaneous effect. A critical examination of possible manifestations of turbulent scattering during occultations of Beta Sco by Jupiter indicates that all observed phenomena during these events can be understood in terms of scintillations produced by turbulence.

Hubbard, W. B.↗

Effective Tree Scattering at L-Band

For routine microwave Soil Moisture (SM) retrieval through vegetation, the tau-omega [1] model [zero-order Radiative Transfer (RT) solution] is attractive due to its simplicity and eases of inversion and implementation. It is the model used in baseline retrieval algorithms for several planned microwave space missions, such as ESA's Soil Moisture Ocean Salinity (SMOS) mission (launched November 2009) and NASA's Soil Moisture Active Passive (SMAP) mission (to be launched 2014/2015) [2 and 3]. These approaches are adapted for vegetated landscapes with effective vegetation parameters tau and omega by fitting experimental data or simulation outputs of a multiple scattering model [4-7]. The model has been validated over grasslands, agricultural crops, and generally light to moderate vegetation. As the density of vegetation increases, sensitivity to the underlying SM begins to degrade significantly and errors in the retrieved SM increase accordingly. The zero-order model also loses its validity when dense vegetation (i.e. forest, mature corn, etc.) includes scatterers, such as branches and trunks (or stalks in the case of corn), which are large with respect to the wavelength. The tau-omega model (when applied over moderately to densely vegetated landscapes) will need modification (in terms of form or effective parameterization) to enable accurate characterization of vegetation parameters with respect to specific tree types, anisotropic canopy structure, presence of leaves and/or understory. More scattering terms (at least up to first-order at L-band) should be included in the RT solutions for forest canopies [8]. Although not really suitable to forests, a zero-order tau-omega model might be applied to such vegetation canopies with large scatterers, but that equivalent or effective parameters would have to be used [4]. This requires that the effective values (vegetation opacity and single scattering albedo) need to be evaluated (compared) with theoretical definitions of these parameters. In a recent study [9], effective vegetation opacity of coniferous trees was compared with two independent estimates of the same parameter. First, a zero-order RT model was fitted to multiangular microwave emissivity data in a least-square sense to provide effective vegetation optical depth as done in spaceborne retrieval algorithms. Second, a ratio between radar backscatter measurements with a corner reflector under trees and in an open area was calculated to obtain measured tree propagation characteristics. Finally, the theoretical propagation constant was determined by forward scattering theorem using detailed measurements of size/angle distributions and dielectric constants of the tree constituents (trunk, branches, and needles). Results indicated that the effective attenuation values are smaller than but of similar magnitude to both the theoretical and measured values. This study will complement the previous work [9] and will focus on characterization of effective scattering albedo by assuming that effective vegetation opacity is same as theoretical opacity. The resultant effective albedo will not be the albedo of single forest canopy element anymore, but it becomes a global parameter, which depends on all the processes taking place within the canopy including multiple scattering as described.

Kurum, Mehmet↗

A Parametric Study of Erupting Flux Rope Rotation: Modeling the 'Cartwheel CME' on 9 April 2008

The rotation of erupting filaments in the solar corona is addressed through a parametric simulation study of unstable, rotating flux ropes in bipolar force-free initial equilibrium. The Lorentz force due to the external shear-field component and the relaxation of tension in the twisted field are the major contributors to the rotation in this model, while reconnection with the ambient field is of minor importance, due to the field's simple structure. In the low-beta corona, the rotation is not guided by the changing orientation of the vertical field component's polarity inversion line with height. The model yields strong initial rotations which saturate in the corona and differ qualitatively from the profile of rotation vs. height obtained in a recent simulation of an eruption without preexisting flux rope. Both major mechanisms writhe the flux rope axis, converting part of the initial twist helicity, and produce rotation profiles which, to a large part, are very similar within a range of shear-twist combinations. A difference lies in the tendency of twist-driven rotation to saturate at lower heights than shear-driven rotation. For parameters characteristic of the source regions of erupting filaments and coronal mass ejections, the shear field is found to be the dominant origin of rotations in the corona and to be required if the rotation reaches angles of order 90 degrees and higher; it dominates even if the twist exceeds the threshold of the helical kink instability. The contributions by shear and twist to the total rotation can be disentangled in the analysis of observations if the rotation and rise profiles are simultaneously compared with model calculations. The resulting twist estimate allows one to judge whether the helical kink instability occurred. This is demonstrated for the erupting prominence in the "Cartwheel CME" on 9 April 2008, which has shown a rotation of approximately 115 deg. up to a height of 1.5 Solar R above the photosphere. Out of a range of initial equilibria which include strongly kink-unstable (Phi = 5 pi), weakly kink-unstable (Phi = 3.5 pi), and kink-stable (Phi = 2.5 pi) configurations, only the evolution of the weakly kink-unstable flux rope matches the observations in their entirety.

initiation and propagation↗

Mirages and the nature of Pluto's atmosphere

We present model occultation lightcurves demonstrating that a strong thermal inversion layer at the base of Pluto's stratosphere can reproduce the minimum flux measured by the Kuiper Airborne Observatory (KAO) during the 1988 occultation of a star by Pluto. The inversion layer also forms the occultation equivalent of a mirage at a radius of 1198 km, which is capable of hiding tropospheres of significant depth. Pluto's surface lies below 1198 km, its radius depending on the depth of the troposphere. We begin by computing plausible temperature structures for Pluto's lower atmosphere, constrained by a calculation of the temperature of the atmosphere near the surface. We then trace rays from the occulted star through the model atmosphere, computing the resultant bending of the ray. Model light curves are obtained by summing the contribution of individual rays within the shadow of Pluto on Earth. We find that we can reproduce the KAO lightcurve using model atmospheres with a temperature inversion and no haze. We have explored models with tropospheres as deep as 40 km (implying a Pluto radius of 1158 km) that reproduce the suite of occultation data. Deeper tropospheres can be fitted to the data, but the mutual event radius of 1150 km probably provides a lower bound. If Pluto has a shallow or nonexistent troposphere, its density is consistent with formation in the solar nebula with modest water loss due to impact ejection. If the troposhere is relatively deep, implying a smaller radius and larger density, significant amounts of water loss are required.

Stansberry, J. A.↗

Boundary Layer Aerosol Composition over Sierra Nevada Mountains using 9.11- and 10.59-micron CW Lidars and Modeled Backscatter from Size Distribution Data

An inversion technique has been developed to determine volume fractions of an atmospheric aerosol composed primarily of ammonium sulfate and ammonium nitrate and water combined with fixed concentration of elemental and organic carbon. It is based on measured aerosol backscatter obtained with 9.11 - and 10.59-micron wavelength continuous wave CO2 lidars and modeled backscatter from aerosol size distribution data. The technique is demonstrated during a flight of the NASA DC-8 aircraft over the Sierra Nevada Mountain Range, California on 19 September, 1995. Volume fraction of each component and effective complex refractive index of the composite particle were determined assuming an internally mixed composite aerosol model. The volume fractions were also used to re-compute aerosol backscatter, providing good agreement with the lidar-measured data. The robustness of the technique for determining volume fractions was extended with a comparison of calculated 2.1,-micron backscatter from size distribution data with the measured lidar data converted to 2.1,-micron backscatter using an earlier derived algorithm, verifying the algorithm as well as the backscatter calculations.

Cutten, D. R.↗

Estimates of surface roughness derived from synthetic aperture radar (SAR) data

Radar remote sensing data provide a unique perspective of the earth's crust and the processes that have influenced its evolution. Physically based models are required, however, to relate the geophysical quantities being measured by the radar sensor to useful geologic information. In this study, synthetic aperture radar (SAR) data over the Cima volcanic field in the Mojave Desert of California are quantitatively connected with microtopography through inversion of a radar backscatter model. Changes in surface roughness inferred from the derived microtopography are modeled and found to be consistent with aeolian mantling as surfaces age. Estimated rates of aeolian deposition for the Cima area are compared to the Lunar Crater volcanic field in Nevada. Rates of deposition appear to be higher at Cima volcanic field, most likely because of its proximity to Soda Lake, the main source of the aeolian material.

Evans, Diane L.↗