Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximate inverses”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

An approximation theory for the identification of nonlinear distributed parameter systems

An abstract approximation framework for the identification of nonlinear distributed parameter systems is developed. Inverse problems for nonlinear systems governed by strongly maximal monotone operators (satisfying a mild continuous dependence condition with respect to the unknown parameters to be identified) are treated. Convergence of Galerkin approximations and the corresponding solutions of finite dimensional approximating identification problems to a solution of the original finite dimensional identification problem is demonstrated using the theory of nonlinear evolution systems and a nonlinear analog of the Trotter-Kato appproximation result for semigroups of bounded linear operators. The nonlinear theory developed here is shown to subsume an existing linear theory as a special case. It is also shown to be applicable to a broad class of nonlinear elliptic operators and the corresponding nonlinear parabolic partial differential equations to which they lead. An application of the theory to a quasilinear model for heat conduction or mass transfer is discussed.

Banks, H. T.↗
DOE OSTI2021

Hydraulic Testing of the Lost Lake Aquifer Near Recovery Well RWM019

An aquifer pumping test was conducted on the Lost Lake Aquifer Zone (LLAZ) at the recently installed recovery well RWM019 in accordance with the approved test plan (Dixon, 2020). The objective of the testing was to determine baseline well performance parameters and aquifer hydraulic conductivity. This testing consisted of a step-drawdown test to determine well performance properties, a constant pumping rate aquifer test to determine aquifer hydraulic properties, and a post-test aquifer recovery monitoring period also used to estimate aquifer hydraulic properties. Well performance parameters determined included specific capacity, well efficiency, and head loss coefficients. The specific capacity of RWM019 was determined to be approximately 2.6 gpm/ft of drawdown based on the final step of the step-drawdown test. Well efficiency was inversely related to pumping rate and decreased from 87% to 81% over a pumping range of approximately 49 to 79 gpm. The aquifer head loss coefficient was determined to be 2.3 ft/ft 3 /min and the well loss coefficient was determined to be 0.05 min 2 /ft 5 .

42 ENGINEERING↗
NASA NTRS1986

An inverse method for subcritical flows

An exact method is presented for two-dimensional subsonic flow within the limitations of the tangent gas approximation. While Woods (1961) studied these equations and proposed iterative methods for the solution of both the iterative and inverse problems, the inverse method here presented is noniterative and exact. It is shown that the direct Euler solution over the designed airfoil is very close to the input speed distribution, and that the constraints necessitated by upstream condition and closure requirements are easily incorporated.

Daripa, P. K.↗
NASA NTRS1994

Improvements to TOVS retrievals over sea ice and applications to estimating Arctic energy fluxes

Modeling studies suggest that polar regions play a major role in modulating the Earth's climate and that they may be more sensitive than lower latitudes to climate change. Until recently, however, data from meteorological stations poleward of 70 degs have been sparse, and consequently, our understanding of air-sea-ice interaction processes is relatively poor. Satellite-borne sensors now offer a promising opportunity to observe polar regions and ultimately to improve parameterizations of energy transfer processes in climate models. This study focuses on the application of the TIROS-N operational vertical sounder (TOVS) to sea-ice-covered regions in the nonmelt season. TOVS radiances are processed with the improved initialization inversion ('3I') algorithm, providng estimates of layer-average temperature and moisture, cloud conditions, and surface characteristics at a horizontal resolution of approximately 100 km x 100 km. Although TOVS has flown continuously on polar-orbiting satellites since 1978, its potential has not been realized in high latitudes because the quality of retrievals is often significantly lower over sea ice and snow than over the surfaces. The recent availability of three Arctic data sets has provided an opportunity to validate TOVS retrievals: the first from the Coordinated Eastern Arctic Experiment (CEAREX) in winter 1988/1989, the second from the LeadEx field program in spring 1992, and the third from Russian drifting ice stations. Comparisons with these data reveal deficiencies in TOVS retrievals over sea ice during the cold season; e.g., ice surface temperature is often 5 to 15 K too warm, microwave emissivity is approximately 15% too low at large view angles, clear/cloudy scenes are sometimes misidentified, and low-level inversions are often not captured. In this study, methods to reduce these errors are investigated. Improvements to the ice surface temperature retrieval have reduced rms errors from approximately 7 K to 3 K; correction of microwaves (50 GHz) brightness temperatures for view angle dependence now allows the surface type (open water versus sea ice) to be determined; modifications to the clear/cloud tests have improved cloud detection over sea ice, especially where it is inhomogeneous; and improved surface temperature estimates have resulted in more successful diagnoses of low-level stratification. Applications of improved retrievals to studies of the Arctic energy budget are encouraging. Preliminary calculations of air-ice stress vectors, 10-m wind speed, and horizontal advection of heat and moisture are presented.

Francis, Jennifer A.↗
DOE OSTI2023

Applying Bayesian inference and deterministic anisotropy to retrieve the molecular structure ∣Ψ(R)∣2 distribution from gas-phase diffraction experiments

Abstract Currently, our general approach to retrieving molecular structures from ultrafast gas-phase diffraction heavily relies on complex ab initio electronic or vibrational excited state simulations to make conclusive interpretations. Without such simulations, inverting this measurement for the structural probability distribution is typically intractable. This creates a so-called inverse problem. Here we address this inverse problem by developing a broadly applicable method that approximates the molecular frame structure ∣Ψ( R , t )∣ 2 distribution independent of these complex simulations. We retrieve the vibronic ground state ∣Ψ( R )∣ 2 for both simulated stretched NO 2 and measured N 2 O. From measured N 2 O, we observe 40 mÅ coordinate-space resolution from 3.75 Å −1 reciprocal space range and poor signal-to-noise, a 50X improvement over traditional Fourier transform methods. In simulated NO 2 diffraction experiments, typical to high signal-to-noise levels predict 100–1000X resolution improvements, down to 0.1 mÅ. By directly measuring the width of ∣Ψ( R )∣ 2 , we open ultrafast gas-phase diffraction capabilities to measurements beyond current analysis approaches. This method has the potential to effectively turn gas-phase ultrafast diffraction into a discovery-oriented technique to probe systems that are prohibitively difficult to simulate.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

The concept of spin ice graphs and a field theory for their charges

Originally detected in rare earth pyrochlores, spin ice physics is now being artificially extended to a variety of geometries that control collective behavior and exotic properties, making graph theory their proper framework. We relate spin ice notions, such as ice rule, ice manifold, Coulomb phases, charges, and monopoles, to graph-theoretical notions, such as balance, in/out-degrees, and Euler paths. We then propose a field-theoretical treatment in which topological charges and monopoles are the degrees of freedom, while the binary spins are subsumed in an entropic interaction among charges. We show that for a spin ice on a graph in a Gaussian approximation, the kernel of the entropic interaction is the inverse of the graph Laplacian, and we compute screening functions from the graph spectra as Green operators for the screened Poisson problem on a graph. We then apply the treatment to star graphs, tournaments, cycles, and regular spin ice in different dimensions. Our aim is twofold: to set spin ice physics in a proper graph setting, where only topological rather than geometrical notions hold, and to invite graph theorists to contribute their powerful tools to the field of spin ice.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Machine learning-based real-time kinetic profile reconstruction in DIII-D

Abstract Kinetic equilibrium reconstruction plays a vital role in the physical analysis of plasma stability and control in fusion tokamaks. However, the traditional approach is subjective and prone to human biases. To address this, the consistent automatic kinetic equilibrium reconstruction (CAKE) method was introduced, providing objective results. Nonetheless, its offline nature limits its application in real-time plasma control systems (PCSs). To address this limitation, we present RTCAKENN, a machine learning model that approximates 7 CAKE-level output profiles, namely pressure, inverse q , toroidal current density, electron temperature and density, carbon ion impurity temperature and rotation profiles, using real-time available inputs. The deep neural network consists of an encoder layer, where the scalars and interdependent inputs such as plasma boundary coordinates and motional Stark effect data are encoded using multi-layer perceptrons (MLPs), while profile inputs are encoded by 1D convolutional layers. The encoded data is passed through a MLP for latent feature extraction, before being decoded in the decoding layers, which consist of upsampling and convolutional layers. RTCAKENN has been implemented in the DIII-D PCS and our model achieves accuracy comparable to CAKE and surpasses existing real-time alternatives. Through clever dropout training, RTCAKENN exhibits robustness and can operate even in the absence of Thomson scattering data or charge exchange recombination data. It executes in under 8 ms in the real-time environment, enabling future application in real-time control and analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1975

Doppler measurements of the ionosphere on the occasion of the Apollo-Soyuz test project. Part 2: Inversion of differential and rotating Doppler shifts

The preparation of the analytical approach and of the related software used in the inversion of the differential and rotating Doppler data obtained from the ionospheric experiment of the Apollo-Soyuz Test Project (ASTP) is discussed. These data were collected in space-to-space paths (between the ASTP Docking Module (DM) and the Apollo Command Service Module and in space-to-ground paths (between the DM and ground). The Doppler links operated at 162 and 324 MHz and have an accuracy better than 3 MHz over 10-sec integration time. The inversion approach was tested with dummy data obtained with a computer simulation. It was found that a measurement accuracy of 1 to 10% in the value of the horizontal electron density gradient at 221-km altitude can be achieved, in space-to-space paths. For space-to-ground paths near the orbital plane, possible effects of the horizontal gradients on the received differential Doppler shifts were identified. It was possible to reduce the gradient-associated errors in the inversion that leads to the columnar electron content by approximately one-half. Accuracies of 5 to 10% in columnar electron content are achievable, with this gradient-compensation technique.

Gay, R. H.↗
NASA NTRS1975

Forward marching procedure for separated boundary-layer flows

A forward-marching procedure for separated boundary-layer flows which permits the rapid and accurate solution of flows of limited extent is presented. The streamwise convection of vorticity in the reversed flow region is neglected, and this approximation is incorporated into a previously developed (Carter, 1974) inverse boundary-layer procedure. The equations are solved by the Crank-Nicolson finite-difference scheme in which column iteration is carried out at each streamwise station. Instabilities encountered in the column iterations are removed by introducing timelike terms in the finite-difference equations. This provides both unconditional diagonal dominance and a column iterative scheme, found to be stable using the von Neumann stability analysis.

Carter, J. E.↗
NASA NTRS1976

Aerosol size distributions - Remote determination from air-borne measurements of the solar aureole

Aerosol-size distributions were determined remotely using an airborne scanning photometer. This instrument consists of a sunshade, optics, and a large-dynamic-range photodetector-amplifier system that scans about the vertical, recording both the direct attenuated solar beam and the diffuse sky radiance in the almucantar of the sun at an effective wavelength of 0.54 micron. Inversion of the measured radiances determined aerosol radii from approximately 0.2 to 8.0 microns. Data were taken during flights over western Lake Superior and southwestern Wisconsin. Inverted aerosol-size distributions are compared with results obtained by other methods.

Twitty, J. T.↗
NASA NTRS1986

The possible importance of synchrotron/inverse Compton losses to explain fast mm-wave and hard X-ray emission of a solar event

The solar burst of 21 May 1984, presented a number of unique features. The time profile consisted of seven major structures (seconds), with a turnover frequency of greater than or approximately 90 GHz, well correlated in time to hard X-ray emission. Each structure consisted of multiple fast pulses (0.1 seconds), which were analyzed in detail. A proportionality between the repetition rate of the pulses and the burst fluxes at 90 GHz and greater than or approximately 100 keV hard X-rays, and an inverse proportionality between repetition rates and hard X-ray power law indices were found. A synchrotron/inverse Compton model was applied to explain the emission of the fast burst structures, which appear to be possible for the first three or four structures.

Correia, E.↗
NASA NTRS1989

A multigrid and upwind viscous flow solver on 3-D embedded and overlapped grids

A numerically efficient method is presented for solving the three-dimensional governing equations of the viscous compressible flow about complex configurations with topologically different components. The physical domain is decomposed into regions for which the grid generation is relatively simple and virtually with no significant restrictions. The Navier-Stokes equations are solved by an implicit, approximately factored, upwind, finite-volume scheme. The block inversions and the diagonalized scalar inversions of the coefficient matrices are modified to allow the holes created in the computational domain by the embedded and overlapped grids. The convergence is accelerated by a multigrid algorithm despite the existence of such holes. The solution for s supersonic flow past a blunt-nose-cylinder at high angle-of-attack is obtained using a C-O grid embedded in a global Cartesian grid.

Baysal, Oktay↗
NASA NTRS2011

The Structure of Titan's Atmosphere from Cassini Radio Occultations

We present results from the two radio occultations of the Cassini spacecraft by Titan in 2006, which probed mid-southern latitudes. Three of the ingress and egress soundings occurred within a narrow latitude range, 31.34 deg S near the surface, and the fourth at 52.8 deg S. Temperature - altitude profiles for all four occultation soundings are presented, and compared with the results of the Voyager 1 radio occultation (Lindal et al., 1983), the HASI instrument on the Huygens descent probe (Fulchignoni et al., 2005), and Cassini CIRS results (Flasar et al., 2005; Achterberg et al., 2008b). Sources of error in the retrieved temperature - altitude profiles are also discussed, and a major contribution is from spacecraft velocity errors in the reconstructed ephemeris. These can be reduced by using CIRS data at 300 km to make along-track adjustments of the spacecraft timing. The occultation soundings indicate that the temperatures just above the surface at 31-34 deg S are about 93 K, while that at 53 deg S is about 1 K colder. At the tropopause, the temperatures at the lower latitudes are all about 70 K, while the 53 deg S profile is again 1 K colder. The temperature lapse rate in the lowest 2 km for the two ingress (dawn) profiles at 31 and 33 deg S lie along a dry adiabat except within approximately 200m of the surface, where a small stable inversion occurs. This could be explained by turbulent mixing with low viscosity near the surface. The egress profile near 34 deg S shows a more complex structure in the lowest 2 km, while the egress profile at 53 deg S is more stable.

Schinder, Paul J.↗
NASA NTRS1999

Inversion of Gravity Fields From the Spacecraft Orbital Data Using an Adjoint Operator Approach

In perturbation approximation, the forward problem of orbital dynamics (equations with initial conditions) is linear with respect to variations of coordinates and/or velocities of the spacecraft and to corresponding variations of the gravity field in the models used. The linear operator adjoint to the linear operator of such forward problem turns out to be instrumental in inversion of differences between observed and predicted coordinates/velocities in terms of the updates of harmonics in the initial gravity field model. Based on this approach, the solution of resulting adjoint problem of orbital dynamics can be used to directly evaluate the matrix of partial derivatives of observable differences with respect to the gravity field harmonics. General discussion of the adjoint problem of orbital dynamics is given and an example of a mathematical formalism for the practical retrieval algorithm is presented.

Ustinov, E. A.↗
NASA NTRS2016

Inverse Energy Dispersion of Energetic Ions Observed in the Magnetosheath

We present a case study of energetic ions observed by the Energetic Particle Detector (EPD) on the Magnetospheric Multiscale spacecraft in the magnetosheath just outside the subsolar magnetopause that occurred at 1000 UT on 8 December 2015. As the magnetopause receded inward, the EPD observed a burst of energetic (approximately 50-1000 keV) proton, helium, and oxygen ions that exhibited an inverse dispersion, with the lowest energy ions appearing first. The prolonged interval of fast antisunward flow observed in the magnetosheath and transient increases in the H components of global ground magnetograms demonstrate that the burst appeared at a time when the magnetosphere was rapidly compressed. We attribute the inverse energy dispersion to the leakage along reconnected magnetic field lines of betatron-accelerated energetic ions in the magnetosheath, and a burst of reconnection has an extent of about 1.5 R(sub E) using combined Super Dual Auroral Radar Network radar and EPD observations.

Lee, S. H.↗
NASA NTRS1985

Least squares inversion of Stokes profiles in the presence of velocity gradients

The Auer, Heasley and House Stokes inversion procedure in use at High Altitude Observatory is based on the analytic solution of the equation of transfer for polarized light where the representation of the thermodynamic and magnetic structure of the atmosphere is assumed to have a high degree of invariance, namely, a Milne-Eddington (ME) structure with a constant magnetic field. In the presence of invariance breaking gradients the resultant Stokes profiles are represented only approximately, if at all, by analytic forms. The accuracy of the inversion parameters and their significance as measures of actual structure are explored for the ME and the Landman-Finn sunspot models under the effects of velocity gradients. The resulting field parameters are good to a few percent and prove to be insensitive to the errors committed by the use of a ME-representation, but the resulting ME parameters yield a less precise measure of thermal structure.

Skumanich, A.↗
NASA NTRS1975

Asymptotic form of the charge-exchange cross section in three-body rearrangement collisions

A three-body general-type rearrangement collision is considered in which the initial and final bound states are described by hydrogen-like wave functions. It is shown that the charge-exchange amplitude in the first Born approximation can be expanded at all incident energies in terms of the inverse powers of the principal quantum number (n). By expanding the exchange amplitude in this way, it is demonstrated conclusively that the cross section for capture into the s, p, and d states as well as for the sum over all the angular-momentum states is proportional to 1/n-cubed plus terms proportional to higher inverse odd powers of n. It is found that the low-lying levels cannot be scaled to the 1/n-cubed law irrespective of the value of the incident energy except in the case of capture into the s states in accordance with the Oppenheimer-Brinkman-Kramers approximation. Zeros and minima in the differential cross sections are given in the limit of high principal quantum number for electron capture by protons from atomic hydrogen and for positronium formation by proton-atomic hydrogen collisions.

Omidvar, K.↗
DOE OSTI2025

Closed-Form Approximation of the Total Variation Proximal Operator

Total variation (TV) is a widely used function for regularizing imaging inverse problems that is particularly appropriate for images whose underlying structure is piecewise constant. TV regularized optimization problems are typically solved using proximal methods, but the way in which they are applied is constrained by the absence of a closed-form expression for the proximal operator of the TV function. A closed-form approximation of the TV proximal operator has previously been proposed, but its accuracy was not theoretically explored in detail. Here, we address this gap by making several new theoretical contributions, proving that the approximation leads to a proximal operator of some convex function, it is equivalent to a gradient descent step on a smoothed version of TV, and that its error can be fully characterized and controlled with its scaling parameter. We experimentally validate our theoretical results on image denoising and sparse-view computed tomography (CT) image reconstruction.

97 MATHEMATICS AND COMPUTING↗