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 55 records · Page 3

Current Induced Forces on Adsorbates at Conducting Carbon Nanotubes

The problem addressed is force on an impurity due to electric current. In order to determine the size and direction of the force, as well as its dependence on the voltage, LCAO (Linear Combination of Atomic Orbitals) and the Green's Function approach were implemented. The resultant induced mean force was directed outwards in addition to a smaller component in the electron flow direction. At low bias, force is larger for larger radius, and crosses over the opposite trend when increasing the bias. The crossover current is approximately inversely proportional to the radius.

Mingo, N.↗

Parallel Preconditioning for CFD Problems on the CM-5

Up to today, preconditioning methods on massively parallel systems have faced a major difficulty. The most successful preconditioning methods in terms of accelerating the convergence of the iterative solver such as incomplete LU factorizations are notoriously difficult to implement on parallel machines for two reasons: (1) the actual computation of the preconditioner is not very floating-point intensive, but requires a large amount of unstructured communication, and (2) the application of the preconditioning matrix in the iteration phase (i.e. triangular solves) are difficult to parallelize because of the recursive nature of the computation. Here we present a new approach to preconditioning for very large, sparse, unsymmetric, linear systems, which avoids both difficulties. We explicitly compute an approximate inverse to our original matrix. This new preconditioning matrix can be applied most efficiently for iterative methods on massively parallel machines, since the preconditioning phase involves only a matrix-vector multiplication, with possibly a dense matrix. Furthermore the actual computation of the preconditioning matrix has natural parallelism. For a problem of size n, the preconditioning matrix can be computed by solving n independent small least squares problems. The algorithm and its implementation on the Connection Machine CM-5 are discussed in detail and supported by extensive timings obtained from real problem data.

Simon, Horst D.↗

A Portable MPI Implementation of the SPAI Preconditioner in ISIS++

A parallel MPI implementation of the Sparse Approximate Inverse (SPAI) preconditioner is described. SPAI has proven to be a highly effective preconditioner, and is inherently parallel because it computes columns (or rows) of the preconditioning matrix independently. However, there are several problems that must be addressed for an efficient MPI implementation: load balance, latency hiding, and the need for one-sided communication. The effectiveness, efficiency, and scaling behavior of our implementation will be shown for different platforms.

Barnard, Stephen T.↗

Recent Progress in Parallel Schur Complement Preconditioning for Computational Fluid

We consider preconditioning methods for nonself-adjoint advective-diffusive systems based on a nonoverlapping Schur complement procedure for arbitrary triangulated domains. The triangulation is first partitioned using the METIS multi-level $k$-way partitioning code. This partitioning of the triangulation induces a natural 2x2 partitioning of the demoralization matrix. By considering various inverse approximations of the 2x2 system we have developed a family of robust preconditioning techniques. The performance of these approximations will be discussed and numerous examples shown to illustrate the efficiency of the technique.

Barth, Tim↗

Progress in Parallel Schur Complement Preconditioning for Computational Fluid Dynamics

We consider preconditioning methods for nonself-adjoint advective-diffusive systems based on a non-overlapping Schur complement procedure for arbitrary triangulated domains. The ultimate goal of this research is to develop scalable preconditioning algorithms for fluid flow discretizations on parallel computing architectures. In our implementation of the Schur complement preconditioning technique, the triangulation is first partitioned into a number of subdomains using the METIS multi-level k-way partitioning code. This partitioning induces a natural 2X2 partitioning of the p.d.e. discretization matrix. By considering various inverse approximations of the 2X2 system, we have developed a family of robust preconditioning techniques. A computer code based on these ideas has been developed and tested on the IBM SP2 and the SGI Power Challenge array using MPI message passing protocol. A number of example CFD calculations will be presented to illustrate and assess various Schur complement approximations.

Barth, Timothy J.↗

Plasma Propulsion Research at NASA Marshall Space Flight Center

The Propulsion Research Center at NASA Marshall Space Flight Center is pursuing a range of research efforts aimed at identifying and developing new technologies for primary spacecraft propulsion. Efficient high-power electric propulsion (Ep) thrusters are a particular area of emphasis; these would enable the relatively rapid transit of large payloads about the solar system for unmanned or manned science and exploration. Such a mission would make heavy demands on the propulsion system, which may be required to run reliably for several years at a specific impulse approaching 10,OOO s with an efficiency of turning electrical power into jet power of at least 70%. The transit time to a destination scales approximately inversely with the cube root of the specific power, which is the ratio of jet power to power-plant mass. Consequently, reducing a trip time by half requires roughly an eight-fold increase in specific power. Given a renewed NASA commitment to space nuclear power, developing efficient EP thrusters with high jet power (> 100 kW) would seem to provide the most direct means of significantly increasing the specific power and hence reducing trip times. In particular, electromagnetic devices, with their high inherent thrust densities, should be better suited to high power applications than thrusters which depend exclusively on electrostatic forces for propellant acceleration.

Sheehy, Jeffrey A.↗

Bubble Eliminator Based on Centrifugal Flow

The fluid bubble eliminator (FBE) is a device that removes gas bubbles from a flowing liquid. The FBE contains no moving parts and does not require any power input beyond that needed to pump the liquid. In the FBE, the buoyant force for separating the gas from the liquid is provided by a radial pressure gradient associated with a centrifugal flow of the liquid and any entrained bubbles. A device based on a similar principle is described in Centrifugal Adsorption Cartridge System (MSC- 22863), which appears on page 48 of this issue. The FBE was originally intended for use in filtering bubbles out of a liquid flowing relatively slowly in a bioreactor system in microgravity. Versions that operate in normal Earth gravitation at greater flow speeds may also be feasible. The FBE (see figure) is constructed as a cartridge that includes two concentric cylinders with flanges at the ends. The outer cylinder is an impermeable housing; the inner cylinder comprises a gas-permeable, liquid-impermeable membrane covering a perforated inner tube. Multiple spiral disks that collectively constitute a spiral ramp are mounted in the space between the inner and outer cylinders. The liquid enters the FBE through an end flange, flows in the annular space between the cylinders, and leaves through the opposite end flange. The spiral disks channel the liquid into a spiral flow, the circumferential component of which gives rise to the desired centrifugal effect. The resulting radial pressure gradient forces the bubbles radially inward; that is, toward the inner cylinder. At the inner cylinder, the gas-permeable, liquid-impermeable membrane allows the bubbles to enter the perforated inner tube while keeping the liquid in the space between the inner and outer cylinders. The gas thus collected can be vented via an endflange connection to the inner tube. The centripetal acceleration (and thus the radial pressure gradient) is approximately proportional to the square of the flow speed and approximately inversely proportional to an effective radius of the annular space. For a given FBE geometry, one could increase the maximum rate at which gas could be removed by increasing the rate of flow to obtain more centripetal acceleration. In experiments and calculations oriented toward the original microgravitational application, centripetal accelerations between 0.001 and 0.012 g [where g normal Earth gravitation (.9.8 m/s2)] were considered. For operation in normal Earth gravitation, it would likely be necessary to choose the FBE geometry and the rate of flow to obtain centripetal acceleration comparable to or greater than g.

Gonda, Steve R.↗

Adaptive Sampling of Spatiotemporal Phenomena with Optimization Criteria

This work was designed to find a way to optimally (or near optimally) sample spatiotemporal phenomena based on limited sensing capability, and to create a model that can be run to estimate uncertainties, as well as to estimate covariances. The goal was to maximize (or minimize) some function of the overall uncertainty. The uncertainties and covariances were modeled presuming a parametric distribution, and then the model was used to approximate the overall information gain, and consequently, the objective function from each potential sense. These candidate sensings were then crosschecked against operation costs and feasibility. Consequently, an operations plan was derived that combined both operational constraints/costs and sensing gain. Probabilistic modeling was used to perform an approximate inversion of the model, which enabled calculation of sensing gains, and subsequent combination with operational costs. This incorporation of operations models to assess cost and feasibility for specific classes of vehicles is unique.

Chien, Steve A.↗

Conditions in the WR 140 Wind-collision Region Revealed by the 1.083-um He I Line Profile

We present spectroscopy of the P Cygni profile of the 1.083-μm He i line in the WC7 + O5 colliding-wind binary (CWB) WR140 (HD 193793), observed in 2008, before its periastron passage in 2009, and in 2016–17, spanning the subsequent periastron passage. Both absorption and emission components showed strong variations. The variation of the absorption component as the O5 star was occulted by the wind-collision region (WCR) sets a tight constraint on its geometry. While the sightline to the O5 star traversed the WCR, the strength and breadth of the absorption component varied significantly on time-scales of days. An emission sub-peak was observed on all our profiles. The variation of its radial velocity with orbital phase was shown to be consistent with formation in the WCR as it swung round the stars in their orbit. Modelling the profile gave a measure of the extent of the sub-peak forming region. In the phase range 0.93–0.99, the flux in the sub-peak increased steadily, approximately inversely proportionally to the stellar separation, indicating that the shocked gas in the WCR where the line was formed was adiabatic. After periastron, the sub-peak flux was anomalously strong and varied rapidly, suggesting formation in clumps down-stream in the WCR. For most of the time, its flux exceeded the 2–10-keV X-ray emission, showing it to be a significant coolant of the shocked wind.

Peredur M Williams↗

Conditions in the WR 140 Wind-Collision Region Revealed by the 1.083-mu m Hei Line Profile

We present spectroscopy of the P Cygni profile of the 1.083-μm He I line in the WC7+O5 colliding-wind binary (CWB) WR 140 (HD 193793), observed in 2008, before its periastron passage in 2009, and in 2016–2017, spanning the subsequent periastron passage. Both absorption and emission components showed strong variations. The variation of the absorption component as the O5 star was occulted by the wind-collision region (WCR) sets a tight constraint on its geometry. While the sightline to the O5 star traversed the WCR, the strength and breadth of the absorption component varied significantly on time-scales of days. An emission subpeak was observed on all our profiles. The variation of its radial velocity with orbital phase was shown to be consistent with formation in the WCR as it swung round the stars in their orbit. Modelling the profile gives a measure of the extent of the subpeak-forming region. In the phase range 0.93–0.99, the flux in the subpeak increased steadily, approximately inversely proportionally to the stellar separation, indicating that the shocked gas in the WCR where the line was formed was adiabatic. After periastron, the subpeak flux was anomalously strong and varied rapidly, suggesting formation in clumps downstream in the WCR. For most of the time, its flux exceeded the 2–10-keV X-ray emission, showing it to be a significant coolant of the shocked wind.

Peredur M Williams↗

Approximation methods for inverse problems involving the vibration of beams with tip bodies

Two cubic spline based approximation schemes for the estimation of structural parameters associated with the transverse vibration of flexible beams with tip appendages are outlined. The identification problem is formulated as a least squares fit to data subject to the system dynamics which are given by a hybrid system of coupled ordinary and partial differential equations. The first approximation scheme is based upon an abstract semigroup formulation of the state equation while a weak/variational form is the basis for the second. Cubic spline based subspaces together with a Rayleigh-Ritz-Galerkin approach were used to construct sequences of easily solved finite dimensional approximating identification problems. Convergence results are briefly discussed and a numerical example demonstrating the feasibility of the schemes and exhibiting their relative performance for purposes of comparison is provided.

Rosen, I. G.↗

Accelerated Training for Large Feedforward Neural Networks

In this paper we introduce a new training algorithm, the scaled variable metric (SVM) method. Our approach attempts to increase the convergence rate of the modified variable metric method. It is also combined with the RBackprop algorithm, which computes the product of the matrix of second derivatives (Hessian) with an arbitrary vector. The RBackprop method allows us to avoid computationally expensive, direct line searches. In addition, it can be utilized in the new, 'predictive' updating technique of the inverse Hessian approximation. We have used directional slope testing to adjust the step size and found that this strategy works exceptionally well in conjunction with the Rbackprop algorithm. Some supplementary, but nevertheless important enhancements to the basic training scheme such as improved setting of a scaling factor for the variable metric update and computationally more efficient procedure for updating the inverse Hessian approximation are presented as well. We summarize by comparing the SVM method with four first- and second- order optimization algorithms including a very effective implementation of the Levenberg-Marquardt method. Our tests indicate promising computational speed gains of the new training technique, particularly for large feedforward networks, i.e., for problems where the training process may be the most laborious.

Stepniewski, Slawomir W.↗

Inversion of the anomalous diffraction approximation for variable complex index of refraction near unity

The Fymat analytic inversion method for retrieving a particle-area distribution function from anomalous diffraction multispectral extinction data and total area is generalized to the case of a variable complex refractive index m(lambda) near unity depending on spectral wavelength lambda. Inversion tests are presented for a water-haze aerosol model. An upper-phase shift limit of 5 pi/2 retrieved an accurate peak area distribution profile. Analytical corrections using both the total number and area improved the inversion.

Smith, C. B.↗

Approximation Of Multi-Valued Inverse Functions Using Clustering And Sugeno Fuzzy Inference

Finding the inverse of a continuous function can be challenging and computationally expensive when the inverse function is multi-valued. Difficulties may be compounded when the function itself is difficult to evaluate. We show that we can use fuzzy-logic approximators such as Sugeno inference systems to compute the inverse on-line. To do so, a fuzzy clustering algorithm can be used in conjunction with a discriminating function to split the function data into branches for the different values of the forward function. These data sets are then fed into a recursive least-squares learning algorithm that finds the proper coefficients of the Sugeno approximators; each Sugeno approximator finds one value of the inverse function. Discussions about the accuracy of the approximation will be included.

Walden, Maria A.↗

Steady induction effects in geomagnetism. Part 1B: Geomagnetic estimation of steady surficial core motions: A non-linear inverse problem

The problem of estimating a steady fluid velocity field near the top of Earth's core which induces the secular variation (SV) indicated by models of the observed geomagnetic field is examined in the source-free mantle/frozen-flux core (SFI/VFFC) approximation. This inverse problem is non-linear because solutions of the forward problem are deterministically chaotic. The SFM/FFC approximation is inexact, and neither the models nor the observations they represent are either complete or perfect. A method is developed for solving the non-linear inverse motional induction problem posed by the hypothesis of (piecewise, statistically) steady core surface flow and the supposition of a complete initial geomagnetic condition. The method features iterative solution of the weighted, linearized least-squares problem and admits optional biases favoring surficially geostrophic flow and/or spatially simple flow. Two types of weights are advanced radial field weights for fitting the evolution of the broad-scale portion of the radial field component near Earth's surface implied by the models, and generalized weights for fitting the evolution of the broad-scale portion of the scalar potential specified by the models.

Voorhies, Coerte V.↗

A Theoretical and Computational Revisit of Conversions Between Whitham’s F-Function and Equivalent Area

This paper proves mathematically that the integral transforms between Whitham’s F-function and equivalent area are the inverse transforms of each other if and only if the slope of the equivalent area at the origin is zero. This mathematical fact contradicts the accepted unconditional inverse relation between Whitham’s F-function and equivalent area in the sonic boom research literature. Piecewise linear approximations of an F-function and of the second derivative of an equivalent area are used to derive numerical formulas for conversions between Whitham’s F-function and equivalent area. Numerical results are included to show convergence of the numerical conversions as the maximum length of the segments for piecewise linear approximations goes to zero. These numerical conversions are approximately the inverse transforms of each other when the second derivative of an equivalent area is continuous and the slope of the equivalent area at the origin is zero.

Whitham's F-function↗

A systematic experimental investigation of significant parameters affecting model tire hydroplaning

The results of a comprehensive parametric study of model and small pneumatic tires operating on a wet surface are presented. Hydroplaning inception (spin down) and rolling restoration (spin up) are discussed. Conclusions indicate that hydroplaning inception occurs at a speed significantly higher than the rolling restoration speed. Hydroplaning speed increases considerably with tread depth, surface roughness and tire inflation pressure of footprint pressure, and only moderately with increased load. Water film thickness affects spin down speed only slightly. Spin down speed varies inversely as approximately the one-sixth power of film thickness. Empirical equations relating tire inflation pressure, normal load, tire diameter and water film thickness have been generated for various tire tread and surface configurations.

Wray, G. A.↗