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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Spin-lattice relaxation and the calculation of gain, pump power, and noise temperature in ruby

The use of a quantitative analysis of the dominant source of relaxation in ruby spin systems to make predictions of key maser amplifier parameters is described. The spin-lattice Hamiltonian which describes the interaction of the electron spins with the thermal vibrations of the surrounding lattice is obtained from the literature. Taking into account the vibrational anisotropy of ruby, Fermi's rule is used to calculate the spin transition rates between the maser energy levels. The spin population rate equations are solved for the spin transition relaxation times, and a comparison with previous calculations is made. Predictions of ruby gain, inversion ratio, and noise temperature as a function of physical temperature are made for 8.4-GHz and 32-GHz maser pumping schemes. The theory predicts that ruby oriented at 90 deg will have approximately 50 percent higher gain in dB and slightly lower noise temperature than a 54.7-deg ruby at 32 GHz (assuming pump saturation). A specific calculation relating pump power to inversion ratio is given for a single channel of the 32-GHz reflected wave maser.

Lyons, J. R.↗
DOE OSTI2020

Inverse energy transfer in decaying, three-dimensional, non-helical magnetic turbulence due to magnetic reconnection

ABSTRACT It has been recently shown numerically that there exists an inverse transfer of magnetic energy in decaying, non-helical, magnetically dominated, magnetohydrodynamic turbulence in three dimensions (3D). We suggest that magnetic reconnection is the underlying physical mechanism responsible for this inverse transfer. In the two-dimensional (2D) case, the inverse transfer is easily inferred to be due to smaller magnetic islands merging to form larger ones via reconnection. We find that the scaling behaviour is similar between the 2D and 3D cases, i.e. the magnetic energy evolves as t−1, and the magnetic power spectrum follows a slope of k−2. We show that on normalizing time by the magnetic reconnection time-scale, the evolution curves of the magnetic field in systems with different Lundquist numbers collapse on to one another. Furthermore, transfer function plots show signatures of magnetic reconnection driving the inverse transfer. We also discuss the conserved quantities in the system and show that the behaviour of these quantities is similar between the 2D and 3D simulations, thus making the case that the dynamics in 3D could be approximately explained by what we understand in 2D. Lastly, we also conduct simulations where the magnetic field is subdominant to the flow. Here, too, we find an inverse transfer of magnetic energy in 3D. In these simulations, the magnetic energy evolves as t−1.4 and, interestingly, a dynamo effect is observed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A persistent adjoint method with dynamic time-scaling and an application to mass action kinetics

In this article, we consider an optimization problem where the objective function is evaluated at the fixed-point of a contraction mapping parameterized by a control variable, and optimization takes place over this control variable. Since the derivative of the fixed-point with respect to the parameter can usually not be evaluated exactly, an adjoint dynamical system can be used to estimate gradients. Using this estimation procedure, the optimization algorithm alternates between derivative estimation and an approximate gradient descent step. We analyze a variant of this approach involving dynamic time-scaling, where after each parameter update the adjoint system is iterated until a convergence threshold is passed. Here, we prove that, under certain conditions, the algorithm can find approximate stationary points of the objective function. We demonstrate the approach in the settings of an inverse problem in chemical kinetics, and learning in attractor networks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Using the Shallow Strain Tensor to Characterize Deep Geologic Reservoirs

Abstract Storing and recovering water, carbon, and heat from geologic reservoirs is central to managing resources in a changing climate. We tested the hypothesis that the strain tensor caused by injecting or producing fluids can be measured at shallow depths and interpreted to advance understanding of underlying deep aquifers or reservoirs. Geodetic‐grade strainmeters were deployed at 30 m depth overlying the Bartlesville Formation, a 500‐m‐deep sandstone near Tulsa, OK. The strainmeters are 220 m east of injection well 9A completed in a permeable lens at the base of the Bartlesville Formation. Water was injected into well 9A at approximately 1.0 L/s during four tests that ranged in duration from a few hours to a few weeks. The horizontal strain increased (tension) and the circumferential strain was a few times larger than the radial strain. The vertical strain decreased (compression) during injection. Strain rates were approximately 100 nε/day during the first few hours, but the rates decreased and were approximately 10 nε/day during most of the tests. Four independent methods of poroelastic simulation and inversion predict reservoir properties and geometries that are similar to each other and consistent with independent information about the reservoir. All strain interpretations predict that a boundary to the permeable lens occurs beneath the vicinity of the strainmeters, which is consistent with core data from the site. The boundary of the permeable lens is located by matching the vertical, radial and circumferential strains, which demonstrates the value of measuring the strain tensor.

Murdoch, Lawrence C.↗
NASA NTRS1976

Formation of hard X-ray tails

X-ray observations of Cygnus X-1 have shown a power-law spectrum. A natural way that this spectrum might arise is discussed. Inverse Compton scattering of a soft photon flux in a very hot accretion disk, kTe at approximately 100 keV is a possible cause. The slope of the spectrum may vary if source conditions vary.

Eardley, D. M.↗
NASA NTRS1976

An inverse problem design method for branched and unbranched axially symmetrical ducts

This paper concerns the potential flow design of axially symmetrical ducts of both circular and annular cross section with or without wall suction or blowing slots. The objective of the work was to develop a method by which such ducts could be designed with directly prescribed wall pressure variation. Previous axially symmetrical design methods applied only to circular cross sectional ducts and required that the pressure distribution be prescribed along the duct centerline and not along the duct wall. The present method uses an inverse problem approach which extends the method of Stanitz to the axially symmetrical case, and an approximation is used to account for the stagnation point in branched duct designs. Two examples of successful designs of diffusers with suction slots are presented.

Nelson, C. D.↗
NASA NTRS1977

Implications of Navier-Stokes turbulence theory for plasma turbulence

Plasma turbulence is considered within the framework of the Navier-Stokes turbulence theory. Two-dimensional turbulence is discussed, noting inverse cascades, and compared to three-dimensional turbulence. MHD turbulence is described with reference to applications of the Navier-Stokes theory and the possibility of inverse magnetic cascades. Turbulence in Vlasov plasmas is outlined on the basis of the direct interaction approximation developed for Navier-Stokes fluids (Kraichnan, 1958-1959).

Montgomery, D.↗
NASA NTRS1986

Connections between conventional and singular-value-based multi-variable flight control system design techniques

It is proposed that frequency-domain multivariable robustness techniques, when combined with classical multivariable procedures, can offer an additional means of evaluating FCS designs. A lateral-directional FCS for an advanced fighter is used as an example. Robustness to unstructured aircraft-input uncertainties is assessed using purely numerical singular-value procedures. Literal approximations for the singular values of the open-loop plant and controller and for the inverse return difference are shown to provide a means of decomposing and diagnosing robustness problems that are insoluble via purely numerical methods.

Mcruer, D. T.↗
NASA NTRS1997

Mid Infrared Hydrogen Recombination Line Emission from the Maser Star MWC 349A

We have detected and spectrally resolved the mid-IR hydrogen recombination lines H6(alpha)(12.372 micrometers), H7(alpha)(19.062 micrometers), H7(beta)(l1.309 micrometers) and H8(gamma)(12.385 micrometers) from the star MWC349A. This object has strong hydrogen maser emission (reported in the millimeter and submillimeter hydrogen recombination lines from H36(alpha) to H21(alpha)) and laser emission (reported in the H15(alpha), H12(alpha) and H10(alpha) lines). The lasers/masers are thought to arise predominantly in a Keplerian disk around the star. The mid-IR lines do not show evident signs of lasing, and can be well modeled as arising from the strong stellar wind, with a component arising from a quasi-static atmosphere around the disk, similar to what is hypothesized for the near IR (less than or equal to 4 micrometers) recombination lines. Since populations inversions in the levels producing these mid-IR transitions are expected at densities up to approximately 10(exp 11)/cu cm, these results imply either that the disk does not contain high-density ionized gas over long enough path lengths to produce a gain approximately 1, and/or that any laser emission from such regions is small compared to the spontaneous background emission from the rest of the source as observed with a large beam. The results reinforce the interpretation of the far-IR lines as true lasers.

Smith, Howard A.↗
NASA NTRS1988

Flight test results for the Daedalus and Light Eagle human powered aircraft

The results of the flight test program of the Daedalus and Light Eagle human powered aircraft in the winter of 1987/88 are given. The results from experiments exploring the Light Eagle's rigid body and structural dynamics are presented. The interactions of these dynamics with the autopilot design are investigated. Estimates of the power required to fly the Daedalus aircraft are detailed. The system of sensors, signal conditioning boards, and data acquisition equipment used to record the flight data is also described. In order to investigate the dynamics of the aircraft, flight test maneuvers were developed to yield maximum data quality from the point of view of estimating lateral and longitudinal stability derivatives. From this data, structural flexibility and unsteady aerodynamics have been modeled in an ad hoc manner and are used to augment the equations of motion with flexibility effects. Results of maneuvers that were flown are compared with the predictions from the flexibility model. To extend the ad hoc flexibility model, a fully flexible aeroelastic model has been developed. The model is unusual in the approximate equality of many structural natural frequencies and the importance of unsteady aerodynamic effects. the Gossamer Albatross. It is hypothesized that this inverse ground effect is caused by turbulence in the Earth's boundary layer. The diameters of the largest boundary layer eddies (which represent most of the turbulent kinetic energy) are proportional to altitude; thus, closer to the ground, the energy in the boundary layer becomes concentrated in eddies of smaller and smaller diameter. Eventually the eddies become sufficiently small (approximately 0.5 cm) that they trip the laminar boundary layer on the wing. As a result, a greater percentage of the wing area is covered with turbulent flow. Consequently the aircraft's drag and the pow er required both increase as the aircraft flies closer to the ground. The results of the flight test program are examined critically, especially with respect to future applications. Maneuvers that allow the observation of stability derivatives for flexible aircraft are detailed. Considerations for the design of autopilots for future human power aircraft and high-altitude RPV's are discussed. Techniques useful in estimating the power required to fly aircraft of very high lift to drag ratio are described.

Sullivan, R. Bryan↗
NASA NTRS2007

X-Ray, UV, and Optical Observations of Supernova 2006bp with Swift: Detection of Early X-Ray Emission

We present results on the X-ray and optical/UV emission from the Type IIP supernova (SN) 2006bp and the interaction of the SW shock with its environment, obtained with the X-Ray Telescope (XRT) and UV/Optical Telescope (UVOT) on-board the Swift observatory. SN 2006bp is detected in X-rays at a 4.5 sigmalevel of significance in the merged XRT data from days 1 to 12 after the explosion. If the (0.2-10 keV band) X-ray luminosity of L(sub 0.2-10) = (1.8 plus or minus 0.4) x l0(exp 39 ergs s(exp -1) is caused by interaction of the SN shock with circumstellar material (CSM), deposited by a stellar wind from the progenitor's companion star, a mass-loss rate of M is approximately 2x10(exp -6) solar mass yr(exp -1) (v(sub w)/10 km s(exp -l) is inferred. The mass-loss rate is one of the lowest ever recorded for a core-collapse SN and consistent with the non-detection in the radio with the VLA on days 2, 9, and 11 after the explosion. The Swift data further show a fading of the X-ray emission starting around day 12 after the explosion. In combination with a follow-up XMM-Newton observation obtained on day 21 after the explosion, an X-ray rate of decline Lx, varies as t(exp -n) with index n = 1.2 plus or minus 0.6 is inferred. Since no other SN has been detected in X-rays prior to the optical peak and since Type IIP SNe have an extended 'plateau' phase in the optical, we discuss the scenario that the X-rays might be due to inverse Compton scattering of photospheric optical photons off relativistic electrons produced in circumstellar shocks. However, due to the high required value of the Lorentz factor (approximately 10-100), inconsistent with the ejecta velocity inferred from optical line widths, we conclude that Inverse Compton scattering is an unlikely explanation for the observed X-ray emission. The fast evolution of the optical/ultraviolet (1900-5500A) spectral energy distribution and the spectral changes observed with Swift reveal the onset of metal line-blanketing and cooling of the expanding photosphere during the first few weeks after the outburst.

Immler, S.↗
DOE OSTI2024

Advancing attenuation estimation through integration of the Hessian in multiparameter viscoacoustic full-waveform inversion

Accurate seismic attenuation models of subsurface structures not only enhance subsequent migration processes by improving fidelity, resolution, and facilitating amplitude-compliant angle gather generation but also provide valuable constraints on subsurface physical properties. Leveraging full-wavefield information, multiparameter viscoacoustic full-waveform inversion ( Q-FWI) simultaneously estimates seismic velocity and attenuation ( Q) models. However, a major challenge in Q-FWI is the contamination of crosstalk artifacts, where inaccuracies in the velocity model are mistakenly mapped to the inverted attenuation model. While incorporating the Hessian is expected to mitigate these artifacts, the explicit implementation is prohibitively expensive due to its formidable computational cost. In this study, we formulate and develop a Q-FWI algorithm via the Newton-conjugate gradient (CG) framework, where the search direction at each iteration is determined through an internal CG loop. In particular, the Hessian is integrated into each CG step in a matrix-free fashion using the second-order adjoint-state method. We find through synthetic experiments that our Newton-CG Q-FWI significantly mitigates crosstalk artifacts compared with the limited-memory Broyden-Fletcher-Goldfarb-Shanno method and the CG method, albeit with a notable computational cost. In the discussion of several key implementation details, we also determine the significance of the approximate Gauss-Newton Hessian, the second-order adjoint-state method, and the two-stage inversion strategy.

Geochemistry & Geophysics↗
NASA NTRS1987

Inferring spectral reflectances of plant elements by simple inversion of bidirectional reflectance measurements

Inverting previously developed explicit expressions for a vertical architecture, bidirectional reflectances measured over corn viewing from the solar quadrant at azimuths near the principal plane are used to determine the spectral reflectances of plant elements. The leaf reflectance values extracted in three visible bands at viewing zenith angles of 70 deg, 60 deg, and 45 deg agree closely with laboratory-measured reflectances of corn leaves. At viewing zenith angle of 30 deg, the inversion breaks down, inasmuch as the inferred plant element reflectances are too high. Satisfactory results are also achieved when the same approach is applied to bidirectional reflectances measured over potted balsam firs, but when applied to soybeans reflectances, the procedure yields unreasonably high leaf reflectances. The failure in this case is attributed to the nonvertical architecture of the soybean canopy; however, for this canopy, inversion based on horizontal architecture is possible. The bidirectional reflectances measured from the solar quadrant, at viewing angles appreciably far from 'hot spot' viewing, approximately equal in magnitude half of the leaf reflectance. The 0.5 ratio is predicted by a previous analysis of opaque horizontal Lambertian facets, as the asymptotic value for a dense canopy at any viewing angle. For soybeans, this ratio applies very closely at 15 deg viewing zenith angle. The results suggest that inversion based on simple architecture, applying explicit expressions, might be of value, either in itself or as a preliminary step before inversion applying complex models.

Otterman, J.↗
NASA NTRS1979

Elastic lithosphere thickness on the moon from mare tectonic features - A formal inversion

The thickness (T) of the lunar elastic lithosphere at the time (3.6 to 3.8 billion years ago) of the earliest preserved basalt flows in circular mare basins can be estimated by inverting the observed locations of extensional tectonic features in and surrounding the maria. In performing the inversion, the lithosphere is modeled as an elastic shell with a liquid interior, and the basalt load for each mare is approximated by a set of concentric cylinders. To permit solving the forward problem of placing radial limits on the positions of the rilles around a given mare, an additional parameter F, the ratio of the radial stress at the radial limits to the maximum radial stress, is introduced. T and F are chosen to give the best weighted-squares fit of the radial limits to the observations, and are used as the initial values in a linearized matrix inversion to check the resolution and estimate errors. The application of the procedure to three maria with prominent extensional features, Humorum, Orientale, and Serenitatis, gives values of T from about 40 + or - 10 to 50 + or - 10 km, and in each case the linearized matrix equation has an exact inverse.

Comer, R. P.↗
NASA NTRS1988

Computational fluid dynamics and aerothermodynamics

Approximations applicable to the radiating, reacting, and conducting stagnation region of a hypervelocity vehicle were incorporated into a method for rapidly obtaining approximate solutions. This solution utilizes a coordinate system based upon the origin of the radiative losses and includes in a phenomenologically correct manner the effects of chemical and thermal nonequilibrium, and nonequilibrium, nongray radiative transfer. Results were presented which demonstrate the usefulness of the method and indicate which radiation parameters require further study and definition. Excellent comparisons were obtained with published results for the Fire2 data. An axisymmetric nonequilibrium inverse method was modified and extended and used to investigate and compare various vibration dissociation chemistry coupling models and radiative heat transfer approximations. The similarities, differences, and consequences of using these models in the Aero-assist Orbital Transfer Vehicles flight regime will be discussed.

Carlson, Leland A.↗
DOE OSTI2025

Hierarchical Bayesian Inverse Problems: A High-Dimensional Statistics Viewpoint

This paper analyzes hierarchical Bayesian inverse problems using techniques from highdimensional statistics. Furthermore, our analysis leverages a property of hierarchical Bayesian regularizers that we call approximate decomposability to obtain non-asymptotic bounds on the reconstruction error attained by maximum a posteriori estimators. The new theory explains how hierarchical Bayesian models that exploit sparsity, group sparsity, and sparse representations of the unknown parameter can achieve accurate reconstructions in high-dimensional settings.

MAP estimation↗
DOE OSTI2022

A low-rank power iteration scheme for neutron transport criticality problems

Computing effective eigenvalues for neutron transport often requires a fine numerical resolution. Here, the main challenge of such computations is the high memory effort of classical solvers, which limits the accuracy of chosen discretizations. In this work, we derive a method for the computation of effective eigenvalues when the underlying solution has a low-rank structure. This is accomplished by utilizing dynamical low-rank approximation (DLRA), which is an efficient strategy to derive time evolution equations for low-rank solution representations. The main idea is to interpret the iterates of the classical inverse power iteration as pseudo-time steps and apply the DLRA concepts in this framework. In our numerical experiment, we demonstrate that our method significantly reduces memory requirements while achieving the desired accuracy. Analytic investigations show that the proposed iteration scheme inherits the convergence speed of the inverse power iteration, at least for a simplified setting.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Coefficient-to-Basis Network: a fine-tunable operator learning framework for inverse problems with adaptive discretizations and theoretical guarantees

We propose a Coefficient-to-Basis Network (C2BNet), a novel framework for solving inverse problems within the operator learning paradigm. C2BNet efficiently adapts to different discretizations through fine-tuning, using a pre-trained model to significantly reduce computational cost while maintaining high accuracy. Unlike traditional approaches that require retraining from scratch for new discretizations, our method enables seamless adaptation without sacrificing predictive performance. Furthermore, we establish theoretical approximation and generalization error bounds for C2BNet by exploiting low-dimensional structures in the underlying datasets. Our analysis demonstrates that C2BNet adapts to low-dimensional structures without relying on explicit encoding mechanisms, highlighting its robustness and efficiency. To validate our theoretical findings, we conducted extensive numerical experiments that showcase the superior performance of C2BNet on several inverse problems. The results confirm that C2BNet effectively balances computational efficiency and accuracy, making it a promising tool to solve inverse problems in scientific computing and engineering applications.

97 MATHEMATICS AND COMPUTING↗