The product form of inverses of sparse matrices and graph theory.
Geometric interpretation of product form of inverse applied to sparse matrices in linear programming, using graph theory
LINEAR PROGRAMMING↗
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.
Fundamental studies in geodynamics
Progress in modeling instantaneous plate kinematics is reviewed, with emphasis on recently developed models of present day plate motions derived by the systematic inversion of globally distributed data sets. Rivera plate motions, the Caribbean South American boundary, Indian plate deformation, Pacific-North America, seismicity and subduction processes, and the study of slow earthquakes and free oscillations are discussed.
Anderson, D. L.↗
Emission mechanisms of the March 5th 1979 gamma ray transient
This paper reviews the spectral and temporal constraints on the emission mechanisms of the initial burst and subsequent pulses of the March 5, 1979 gamma-transient. Interpretation of the hard tail above 300 keV as due to inverse comptonization, when combined with the synchrotron-pair annihilation model of Ramaty et al. (1981), makes it possible to derive from first principles the intrinsic luminosity of the source. This is in general agreement with that expected from N49 in LMC. The pulsed emissions have characteristics difficult to explain in terms of conventional models.
Liang, E. P. T.↗
An analysis of the role of drift waves in equatorial spread F
An account is given of results of rocket measurements of the wave number spectrum of equatorial spread F irregularities, with emphasis on wavelengths less than 100 m. The measurements were made from two sounding rockets launched from Peru as part of Project Condor. The Condor density fluctuation spectra display a break at a wavelength near 100 m, identical to that found in the PLUMEX experiment (Kelley et al., 1982). The Condor data also confirm a subrange in which the density and the wave potential obey the Boltzmann relation - a strong indication of the presence of low-frequency electrostatic waves with finite wavelength parallel to the magnetic field, perhaps low-frequency drift waves as proposed by Kelley et al. The Condor data are also consistent with the previous conjecture that drift waves only exist above 300 km altitude. To investigate the difference in spectra observed over two altitude ranges, the data must be fitted to a form for the power spectrum taken from Keskinen and Ossakow (1981). The fitted spectrum, along with empirically determined growth and dissipation rates, is used to calculate the energy pumped into the spectrum at long wavelengths as well as the energy dissipated at shorter wavelengths. It is found that the energy is balanced by classical collisional effects in the low-altitude case, but energy balance in the high-altitude case requires an enhanced dissipation of about 500 times that due to classical diffusion. The model is consistent with, but does not uniquely imply, an inverse cascade of drift wave turbulence in equatorial spread F.
Labelle, J.↗
Estimation of canopy parameters for inhomogeneous vegetation canopies from reflectance data. III - TRIM: A model for radiative transfer in heterogeneous three-dimensional canopies
A model for radiative transfer in heterogeneous three-dimensional canopies such as those found in forests is proposed. Its use in estimating important biophysical variables such as leaf area index and canopy architecture from bidirectional canopy reflectance data is discussed. The model and its use in estimating canopy parameters through its inversion are validated with measured canopy reflectance data for corn canopies.
Goel, Narendra S.↗
An observational study of stratocumulus entrainment and thermodynamics
The stratocumulus (SC) entrainment and thermodynamics are studied using data from the Dynamics and Chemistry of Marine Stratocumulus (DYCOMS) experiment. The rate of entrainment of air from the free troposphere into the cloud-topped PBL is estimated using a technique based on the measurement of ozone flux and mean distribution. The average measured value of entrainment rate was found to be 3.0 mm/s with a range of 1.0 to 5.0 mm/s for cloudy cases. Thermodynamic budgets are constructed for eight DYCOMS cases. It was found that the divergence of the solar radiative flux is an important component of the boundary-layer energetics during midday, the factor which must be accounted for in mixed-layer models. The net longwave radiative flux profiles show good agreement with theoretical models, but an unambiguous partitioning of the flux divergence between inversion and mixed layers could not be established.
Kawa, S. R.↗
Visualizing characteristics of ocean data collected during the Shuttle Imaging Radar-B experiment
Topographic measurements of sea surface elevation collected by the Surface Contour Radar (SCR) during NASA's Shuttle Imaging Radar (SIR-B) experiment are plotted as three dimensional surface plots to observe wave height variance along the track of a P-3 aircraft. Ocean wave spectra were computed from rotating altimeter measurements acquired by the Radar Ocean Wave Spectrometer (ROWS). Fourier power spectra computed from SIR-B synthetic aperture radar (SAR) images of the ocean are compared to ROWS surface wave spectra. Fourier inversion of SAR spectra, after subtraction of spectral noise and modeling of wave height modulation, yields topography similar to direct measurements made by SCR. Visual perspectives on the SCR and SAR ocean data are compared. Threshold distinctions between surface elevation and texture modulations of SAR data are considered within the context of a dynamic statistical model of rough surface scattering. The result of these endeavors is insight as to the physical mechanism governing the imaging of ocean waves with SAR.
Tilley, David G.↗
Efficient control with an order (n) recursive inversion of the Jacobian for an n-link serial manipulator
A method is presented to order (n) recursively invert the Jacobian matrix for a serial n-link manipulator. By realizing that the Jacobian relationship involving accelerations is a special case for the forward dynamics problem, it is possible to obtain the Jacobian inverse by directly applying Rodriguez's (1987) recursive forward dynamics algorithm. In this technique, the n-link robot equation is formulated as a spatially recursive algorithm in the form of a filtering and smoothing problem. To compute the Jacobian inverse via this algorithm, the n-link manipulator is modeled with only a unit end-point mass and zero applied torques.
Meldrum, D. R.↗
Coronal holes and solar wind streams during the sunspot cycle
Complementary synoptic observations of the Sun and interplanetary space have been obtained nearly continuously for more than two sunspot cycles and have led to new ideas about the origin of the solar wind. These observations show an inverse correlation between wind speed at Earth and magnetic flux tube expansion in the corona, with fast wind originating from slowly diverging tubes and vice versa. Although this result is consistent with the Skylab-era concept that fast wind originates from the center of a large isolated coronal hole, it implies that the wind may be even faster at the facing edges of like-polarity holes where the flux-tubes converge as they begin their outward extension. Thus, very fast wind ought to originate from the high-latitude edges of the circumpolar holes soon after sunspot maximum and from the mid-latitude necks of the polar-hole lobes during the declining phase of the cycle. The observed inverse correlation may be understood physically in terms of a model in which Alfven waves boost the wind to high speed provided that the wave energy flux is distributed approximately uniformly at the coronal base.
Sheeley, N. R., Jr.↗
Constraints on electron acceleration in the Crab nebula
Using the radio through hard X-ray images of the Crab nebula to derive the spatial dependence of the electron spectrum and the magnetic field distribution from MHD flow models, we have rederived the high-energy spectrum of inverse Compton scattered gamma rays. We find agreement with the observed spectrum at TeV energies, but it is clear that the inverse Compton flux does not contribute significantly to the unpulsed nebular emission observed by COS-B from 50 to 500 MeV, which is consistent with a smooth continuation of the spectrum in hard X-rays. The emission at these energies must therefore be due to synchrotron radiation by electrons of at least PeV energies. It appears that the emission in the high-energy gamma-ray range, sensitive to the highest energy electrons in the nebula, can put interesting constraints on the acceleration mechanism.
Harding, A. K.↗
The 1995 scientific assessment of the atmospheric effects of stratospheric aircraft
This report provides a scientific assessment of our knowledge concerning the impact of proposed high-speed civil transport (HSCT) aircraft on the atmosphere. It comes at the end of Phase 1 of the Atmospheric Effects of Stratospheric Aircraft element of the NASA High-Speed Research Program. The fundamental problem with stratospheric flight is that pollutant residence times are long because the stratosphere is a region of permanent temperature inversion with stable stratification. Using improved two-dimensional assessment models and detailed fleet emissions scenarios, the assessment examines the possible impact of the range of effluents from aircraft. Emphasis is placed on the effects of NO(x) and H2O on the atmospheric ozone content. Measurements in the plume of an in-flight Concorde supersonic transport indicated a large number of small particles. These measurements, coupled with model sensitivity studies, point out the importance of obtaining a more detailed understanding of the fate of sulfur in the HSCT exhaust. Uncertainties in the current understanding of the processes important for determining the overall effects of HSCT's on the atmosphere are discussed and partially quantified. Research directions are identified to improve the quantification of uncertainties and to reduce their magnitude.
Stolarski, Richard S.↗
Structure of magnetic field in the solar wind
This work is concerned with empirical data on magnetic field in the solar wind in frame of a concept of dissipative solar wind, developed in papers (Solar Wind 7 Conf., Pergamon Press, 1992, 165 and 1992 STEP/5th COSPAR Coll. Pergamon Press, 1994, 117; 235; 803). Interplanetary magnetic fields should be classified with respect to their origin. It is very important for all the theoretical problems from the necessity to specify correctly boundary and initial conditions: the magnetic field must be sewed with its source. One should select the field, connected directly with the Sun (stretched out from it), and the field of moving electric currents. It occured central in discussion about the velocity of Alfven waves, probably warming up the solar wind, relative to the Sun, the magnetic field and solar wind plasma. The selection problem corresponds to an inverse problem and obviously has no single solution. The dissipative model of the solar wind introduce the slipping and leakage of plasma relative to magnetic field. There are no 'interplanetary current sheets' in it. But temporal fluctuations from the filamentation of electric currents play the key role. As a whole, the new concept requires the re-interpretation of main objects in the interplanetary magnetic field.
Chertkov, A. D.↗
Concentric Coplanar Capacitive Sensor System with Quantitative Model
A concentric coplanar capacitive sensor includes a charged central disc forming a first electrode, an outer annular ring coplanar with and outer to the charged central disc, the outer annular ring forming a second electrode, and a gap between the charged central disc and the outer annular ring. The first electrode and the second electrode may be attached to an insulative film. A method provides for determining transcapacitance between the first electrode and the second electrode and using the transcapacitance in a model that accounts for a dielectric test piece to determine inversely the properties of the dielectric test piece.
Bowler, Nicola↗
Observations of mode coupling in the solar corona and bipolar noise storms
High-spatial-resolution observations of the sun which reflect on the role of mode coupling in the solar corona, and a number of new observations are presented. It is shown that typically, polarization inversion is seen at 5 GHz in active region sources near the solar limb, but not at 1.5 GHz. Although this is apparently in contradiction to the simplest form of mode coupling theory, it remains consistent with current models for the active region emission. Microwave bursts show no strong evidence for polarization inversion. Bipolar noise storm continuum emission is discussed in some detail, utilizing recent VLA observations at 327 MHz. It is shown that bipolar sources are common at 327 MHz. Further, the trailing component of the bipole is frequently stronger than the leading component, in apparent conflict with the 'leading-spot' hypothesis. The observations indicate that, at 327 MHz, mode coupling is apparently strong at all mode-coupling layers in the solar corona. The 327 MHz observations require a much weaker magnetic field strength in the solar corona to explain this result than did earlier lower-frequency observations: maximum fields are 0.2 G. This is a much weaker field than is consistent with current coronal models.
White, S. M.↗
Damage Diagnosis in Semiconductive Materials Using Electrical Impedance Measurements
Recent aerospace industry trends have resulted in an increased demand for real-time, effective techniques for in-flight structural health monitoring. A promising technique for damage diagnosis uses electrical impedance measurements of semiconductive materials. By applying a small electrical current into a material specimen and measuring the corresponding voltages at various locations on the specimen, changes in the electrical characteristics due to the presence of damage can be assessed. An artificial neural network uses these changes in electrical properties to provide an inverse solution that estimates the location and magnitude of the damage. The advantage of the electrical impedance method over other damage diagnosis techniques is that it uses the material as the sensor. Simple voltage measurements can be used instead of discrete sensors, resulting in a reduction in weight and system complexity. This research effort extends previous work by employing finite element method models to improve accuracy of complex models with anisotropic conductivities and by enhancing the computational efficiency of the inverse techniques. The paper demonstrates a proof of concept of a damage diagnosis approach using electrical impedance methods and a neural network as an effective tool for in-flight diagnosis of structural damage to aircraft components.
Ross, Richard W.↗
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↗
On the Synergy Between Numerics and Subgrid Scale Modeling in LES of Stratified Flows: Grid Convergence of a Stratocumulus-Topped Boundary Layer
The effectiveness of a linear upwinding scalar advection scheme to suppress numerical dispersion errors near sharp inversions in large-eddy simulations of a nocturnal stratocumulus-topped boundary layer is assessed. Linear upwinding is a trade-off between non-dissipative and non-linear positive definite advection schemes. It is shown that linear upwinding does not negatively impact the model's grid convergence properties and a sharp inversion free of numerical artifacts is maintained. Even though mean profiles and turbulence fluxes show good grid convergence characteristics the liquid water amount varies significantly with grid resolution. The entrainment rate is identical for all resolutions and independent of the liquid water amount. For the present stratocumulus case, the impact of cloud-top radiative cooling is negligible and turbulence is largely driven by convection emanating from the surface.
Matheou, Georgios↗
Adding GPU Support to the Markov Chain Monte Carlo Code Catmip
In geophysics, we are confronted with many under-determined inverse problems. For example, all of our observations of earthquakes are made at the Earth’s surface. So, when we try to infer how slip during an earthquake evolves in space and time, we find that there are many potential slip histories that are consistent with our limited observations and our understanding of earthquake physics. One way to approach these problems is with Bayesian analysis which allows us to infer the ensemble of all potential slip models that satisfy the observations and our prior knowledge of earthquake physics. In Bayesian analysis, our prior knowledge is known as the prior probability density function or prior PDF, the fit to the data is known as the data likelihood, and the target PDF that satisfies both the prior PDF and data likelihood is known as the posterior PDF. However, simulating the posterior PDF typically requires using Markov Chain Monte Carlo (MCMC) to draw tens of billions of random realizations of earthquake slip models, which may not be computationally feasible. To make this and similar geophysical inversions computationally tractable, we developed the Cascading Adaptive Transitional Metropolis In Parallel (CATMIP) algorithm. CATMIP is an efficient parallel Markov Chain Monte Carlo (MCMC) sampler that is used for model fitting and uncertainty quantification in geophysics. Example use cases are earthquake rupture modeling, determining mineral composition on Mars, reconstructing the history of ocean salinity, and historical earthquake relocation. CATMIP employs many parallel instances of the Metropolis algorithm for sampling in a transitioning framework. Transitioning is a process in which a set of random samples at equilibrium with a known probability density function (PDF) are used as seeds for the Markov chains to sample successive target PDFs that incrementally move the distribution from the starting seeds to the final desired PDF that describes the relative plausibility of potential values for the model parameters. The algorithm is implemented as a Master-Worker model employing MPI for communication. The worker processes are loosely coupled with global parameters periodically optimized by the master process. This provides a very high amount of parallelism with little communication between updates. During the presentation we will discuss the history of the algorithm and elaborate the earthquake rupture modeling use case for the CATMIP package. Our first step toward GPU optimization was to optimize the code for the CPU. CPU profiling revealed that most of the compute time is spent in calls to level 2 BLAS routines and calls to GSL random number generators. We revised the algorithm to employ level 3 BLAS routines instead. In our presentation we will describe how this was accomplished. Adding GPU support to CATMIP consisted mostly of replacing the calls to GSL with calls to GPU vendor-provided library routines. A small number of loops were directly implemented in CUDA. In the presentation will provide implementation details. Finally, we will discuss methods for profiling and opportunities for further optimizing GPU execution. By creating a code with the flexibility to run on either a CPU or GPU architecture, CATMIP can be used on systems ranging from large CPU-based HPC environments to single servers with GPU acceleration and everything in between.
HECC↗