Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximate computing”

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 1,315 records · Page 73

Clifford Neural Operators on Atmospheric Data Influenced Partial Differential Equations

Mathematical representations of the atmosphere are key to forecasting and research tasks across Earth science. Numerically solving the underlying partial differential equations(PDEs) of the atmosphere, however, can be difficult and computationally expensive with numerous trade-offs between computing efficiency and accuracy. Utilizing neural net-works to learn approximations of the PDE solutions from the data can help us model complex phenomena more efficiently than traditional numerical schemes. Here, we have applied Clifford algebra-based neural operators for predicting atmospheric variables. Clifford Fourier neural operators are used with two different backbone architectures, ResNet and UNet, on custom data of U10, V10 and surface pressure as well as U500, V500 and Z500. Clifford Fourier neural operators, coupled with ResNet and UNet architectures, are applied to a key reanalysis dataset. Model performance is initially strong, but we observe increasing errors, resulting in the model becoming highly unstable.

Sujit Roy↗

Space Shuttle Main Engine performance analysis

For a number of years, NASA has relied primarily upon periodically updated versions of Rocketdyne's power balance model (PBM) to provide space shuttle main engine (SSME) steady-state performance prediction. A recent computational study indicated that PBM predictions do not satisfy fundamental energy conservation principles. More recently, SSME test results provided by the Technology Test Bed (TTB) program have indicated significant discrepancies between PBM flow and temperature predictions and TTB observations. Results of these investigations have diminished confidence in the predictions provided by PBM, and motivated the development of new computational tools for supporting SSME performance analysis. A multivariate least squares regression algorithm was developed and implemented during this effort in order to efficiently characterize TTB data. This procedure, called the 'gains model,' was used to approximate the variation of SSME performance parameters such as flow rate, pressure, temperature, speed, and assorted hardware characteristics in terms of six assumed independent influences. These six influences were engine power level, mixture ratio, fuel inlet pressure and temperature, and oxidizer inlet pressure and temperature. A BFGS optimization algorithm provided the base procedure for determining regression coefficients for both linear and full quadratic approximations of parameter variation. Statistical information relative to data deviation from regression derived relations was also computed. A new strategy for integrating test data with theoretical performance prediction was also investigated. The current integration procedure employed by PBM treats test data as pristine and adjusts hardware characteristics in a heuristic manner to achieve engine balance. Within PBM, this integration procedure is called 'data reduction.' By contrast, the new data integration procedure, termed 'reconciliation,' uses mathematical optimization techniques, and requires both measurement and balance uncertainty estimates. The reconciler attempts to select operational parameters that minimize the difference between theoretical prediction and observation. Selected values are further constrained to fall within measurement uncertainty limits and to satisfy fundamental physical relations (mass conservation, energy conservation, pressure drop relations, etc.) within uncertainty estimates for all SSME subsystems. The parameter selection problem described above is a traditional nonlinear programming problem. The reconciler employs a mixed penalty method to determine optimum values of SSME operating parameters associated with this problem formulation.

Santi, L. Michael↗

On three dimensional object recognition and pose-determination: An abstraction based approach

A method of computing reliable Gaussian and mean curvature sign-map descriptors from the polynomial approximation of surfaces was demonstrated. Such descriptors which are invariant under perspective variation are suitable for hypothesis generation. A means for determining the pose of constructed geometric forms whose algebraic surface descriptors are nonlinear in terms of their orienting parameters was developed. This was done by means of linear functions which are capable of approximating nonlinear forms and determining their parameters. It was shown that biquadratic surfaces are suitable companion linear forms for cylindrical approximation and parameter estimation. The estimates provided the initial parametric approximations necessary for a nonlinear regression stage to fine tune the estimates by fitting the actual nonlinear form to the data. A hypothesis-based split-merge algorithm for extraction and pose determination of cylinders and planes which merge smoothly into other surfaces was developed. It was shown that all split-merge algorithms are hypothesis-based. A finite-state algorithm for the extraction of the boundaries of run-length regions was developed. The computation takes advantage of the run list topology and boundary direction constraints implicit in the run-length encoding.

Quek, Kok How Francis↗

Superthermal electrons and Bernstein waves in Jupiter's inner magnetosphere

A theoretical model for generation of banded electrostatic emissions by low density, superthermal electrons is developed for application to Jupiter's magnetosphere. The model employs a power law form for the energy dependence and a loss cone pitch angle distribution of the superthermals to drive convective instability of Bernstein modes. A direct correspondence between spectral features of the 3/2 band and resonant superthermal electrons is found. The concept of a critical flux of resonant electrons able to provide 10 e-foldings of electric field amplification yields an explicit relation in terms of the background thermal electron pressure. This result is used to construct a theoretical/empirical model of thermal electron density and temperature from 6-20 Jupiter radii in the Jovian magnetosphere which suggests that the electron temperature is less than the ion temperature which is approximately equal to 10 times the electron temperature in this region. Finally, wave ray paths are computed for propagation in the magnetic equator and in the magnetic meridional plane of a dipole magnetic field. These ray paths suggest that intense wave activity is tightly confined to a small latitudinal extent, less than + or - approximately 4 deg, about the magnetic equator.

Barbosa, D. D.↗

An improved computational approach for multilevel optimum design

A penalty-function algorithm employing Newton's method with approximate second derivatives (Haftka and Starnes, 1980) is developed for two-level hierarchical design optimization problems. The difficulties posed by discontinuous behavior in typical multilevel problems are explained and illustrated for the case of a three-bar truss; the algorithm is formulated; and its advantages are demonstrated in the problem of a portal framework having three beams (described by six cross-section parameters), subjected to two loading conditions, and to be constructed in six different materials for comparison. The final design parameters are listed in a table.

Haftka, R. T.↗

High order filtering methods for approximating hyperbolic systems of conservation laws

The essentially nonoscillatory (ENO) schemes, while potentially useful in the computation of discontinuous solutions of hyperbolic conservation-law systems, are computationally costly relative to simple central-difference methods. A filtering technique is presented which employs central differencing of arbitrarily high-order accuracy except where a local test detects the presence of spurious oscillations and calls upon the full ENO apparatus to remove them. A factor-of-three speedup is thus obtained over the full-ENO method for a wide range of problems, with high-order accuracy in regions of smooth flow.

Lafon, F.↗

A numerical algorithm for optimal feedback gains in high dimensional LQR problems

A hybrid method for computing the feedback gains in linear quadratic regulator problems is proposed. The method, which combines the use of a Chandrasekhar type system with an iteration of the Newton-Kleinman form with variable acceleration parameter Smith schemes, is formulated so as to efficiently compute directly the feedback gains rather than solutions of an associated Riccati equation. The hybrid method is particularly appropriate when used with large dimensional systems such as those arising in approximating infinite dimensional (distributed parameter) control systems (e.g., those governed by delay-differential and partial differential equations). Computational advantage of the proposed algorithm over the standard eigenvector (Potter, Laub-Schur) based techniques are discussed and numerical evidence of the efficacy of our ideas presented.

Banks, H. T.↗

Finite element dynamic analysis on CDC STAR-100 computer

Computational algorithms are presented for the finite element dynamic analysis of structures on the CDC STAR-100 computer. The spatial behavior is described using higher-order finite elements. The temporal behavior is approximated by using either the central difference explicit scheme or Newmark's implicit scheme. In each case the analysis is broken up into a number of basic macro-operations. Discussion is focused on the organization of the computation and the mode of storage of different arrays to take advantage of the STAR pipeline capability. The potential of the proposed algorithms is discussed and CPU times are given for performing the different macro-operations for a shell modeled by higher order composite shallow shell elements having 80 degrees of freedom.

Noor, A. K.↗

Performance of FFT methods in local gravity field modelling

Fast Fourier transform (FFT) methods provide a fast and efficient means of processing large amounts of gravity or geoid data in local gravity field modelling. The FFT methods, however, has a number of theoretical and practical limitations, especially the use of flat-earth approximation, and the requirements for gridded data. In spite of this the method often yields excellent results in practice when compared to other more rigorous (and computationally expensive) methods, such as least-squares collocation. The good performance of the FFT methods illustrate that the theoretical approximations are offset by the capability of taking into account more data in larger areas, especially important for geoid predictions. For best results good data gridding algorithms are essential. In practice truncated collocation approaches may be used. For large areas at high latitudes the gridding must be done using suitable map projections such as UTM, to avoid trivial errors caused by the meridian convergence. The FFT methods are compared to ground truth data in New Mexico (xi, eta from delta g), Scandinavia (N from delta g, the geoid fits to 15 cm over 2000 km), and areas of the Atlantic (delta g from satellite altimetry using Wiener filtering). In all cases the FFT methods yields results comparable or superior to other methods.

Forsberg, Rene↗

Solution of complex nonlinear problems by a generalized application of the method of base and comparison solutions with applications to aerodynamics problems

A theory for obtaining approximate solutions to nonlinear problems whose exact solutions require the use of large computational procedures is described. The technique represents in some respects a generalization of the method of base and comparison solutions for flows depending on a parameter. For the generalized problem, the input variable is no longer a parameter but a function that is incremented over its entire domain. After performing calculations for a base configuration and a small number of variations of it, solutions for a large class of configurations can be obtained by forming linear combinations of the solution increments. For a restricted class of problems, approximate solutions can be obtained for general variations of a base configuration by using a function-space derivative estimate obtained from a base solution and a single variation.

Barger, R. L.↗

Computational solution of the defect stream-function equation for nonequilibrium turbulent boundary layers

The derivation of the accurate, second-order, almost linear, approximate equation governing the defect stream function for nonequilibrium compressible turbulent boundary layers is reviewed. The similarity of this equation to the heat conduction equation is exploited in the development of an unconditionally stable, tridiagonal computational method which is second-order accurate in the marching direction and fourth-order accurate in the surface-normal direction. Results compare well with experimental data. Nonlinear effects are shown to be small. This two-dimensional method is simple and has been implemented on a programmable calculator.

Barnwell, Richard W.↗

Ground-based time-guidance algorithm for control of airplanes in a time-metered air traffic control environment: A piloted simulation study

The rapidly increasing costs of flight operations and the requirement for increased fuel conservation have made it necessary to develop more efficient ways to operate airplanes and to control air traffic for arrivals and departures to the terminal area. One concept of controlling arrival traffic through time metering has been jointly studied and evaluated by NASA and ONERA/CERT in piloted simulation tests. From time errors attained at checkpoints, airspeed and heading commands issued by air traffic control were computed by a time-guidance algorithm for the pilot to follow that would cause the airplane to cross a metering fix at a preassigned time. These tests resulted in the simulated airplane crossing a metering fix with a mean time error of 1.0 sec and a standard deviation of 16.7 sec when the time-metering algorithm was used. With mismodeled winds representing the unknown in wind-aloft forecasts and modeling form, the mean time error attained when crossing the metering fix was increased and the standard deviation remained approximately the same. The subject pilots reported that the airspeed and heading commands computed in the guidance concept were easy to follow and did not increase their work load above normal levels.

Knox, C. E.↗

An efficient massively parallel Euler solver for unstructured grids

A data parallel mesh-vertex upwind finite-volume scheme for solving the Euler equations on triangular unstructured meshes is described. A novel vertex-based partitioning of the problem is introduced which minimizes the computation and communication costs associated with distributing the computation to the processors of a massively parallel computer. Finally, the performance of this unstructured computation on 8K processors of the Connection Machine CM-2 is compared with one processor of a Cray-YMP. The experiments show that 8K processors of the CM-2 achieve approximately 70 percent of the performance of one processor of the Cray-YMP on the unstructured mesh computations described here.

Hammond, Steven W.↗

Mesoscale roughness of Venus

The global distribution of multi-kilometer (approx. 9 km) length scale 'roughness' (hereafter mesoscale roughness or MR) on Venus can be estimated from the Magellan global altimetry dataset (GxDR) and then compared with MR data derived for Earth from 5' ETOP5 data and for Mars (from USGS Mars DTM dataset). The mesoscale roughness parameter (MR) represents the RMS variance in meters of the actual planetary surface topography relative to the best fitting tangent plane defined on the basis of a 3x3 pixel sliding window. The best-fit plane was computed using a least-squares solution which minimizes delta H, the sum of the squares of the differences between the 9 local elevation values (H(sub i)), and the elevation of best-fit plane at the same grid location. Using the best-fit plane and delta H, we have computed the RMS 'roughness' var(delta R), where this parameter is always minimized on the basis of its calculation using least squares. We have called this 'ruggedness' parameter the Mesoscale Roughness (MR) because it is directly related to the high-frequency variance of topography after mesoscale slopes and tilts (i.e., for Venus, the baseline over which MR is computed (dx) is approx. 8.8 km and dx for Earth is approx. 9.3 km) are removed. As such, MR represents the degree to which a planetary surface is more rugged than approximately 10 km scale facets or tilts. It should not be confused with the radar 'RMS Roughness' parameter computed at 0.1 to 10 m length scales on the basis of the Magellan radar altimeter echo. We will use our MR parameter to investigate the global ruggedness properties of Venus as they relate to geological provinces and in comparison with the spatial pattern of MR for Earth and Mars.

Garvin, J. B.↗

An Implicit Upwind Algorithm for Computing Turbulent Flows on Unstructured Grids

An implicit, Navier-Stokes solution algorithm is presented for the computation of turbulent flow on unstructured grids. The inviscid fluxes are computed using an upwind algorithm and the solution is advanced in time using a backward-Euler time-stepping scheme. At each time step, the linear system of equations is approximately solved with a point-implicit relaxation scheme. This methodology provides a viable and robust algorithm for computing turbulent flows on unstructured meshes. Results are shown for subsonic flow over a NACA 0012 airfoil and for transonic flow over a RAE 2822 airfoil exhibiting a strong upper-surface shock. In addition, results are shown for 3 element and 4 element airfoil configurations. For the calculations, two one equation turbulence models are utilized. For the NACA 0012 airfoil, a pressure distribution and force data are compared with other computational results as well as with experiment. Comparisons of computed pressure distributions and velocity profiles with experimental data are shown for the RAE airfoil and for the 3 element configuration. For the 4 element case, comparisons of surface pressure distributions with experiment are made. In general, the agreement between the computations and the experiment is good.

Anerson, W. Kyle↗

Stellar model chromospheres. III - Arcturus /K2 III/

Models are constructed for the upper photosphere and chromosphere of Arcturus based on the H, K, and IR triplet lines of Ca II and the h and k lines of Mg II. The chromosphere model is derived from complete redistribution solutions for a five-level Ca II ion and a two-level Mg II ion. A photospheric model is derived from the Ca II wings using first the 'traditional' complete-redistribution limit and then the more realistic partial-redistribution approximation. The temperature and mass column densities for the temperature-minimum region and the chromosphere-transition region boundary are computed, and the pressure in the transition region and corona are estimated. It is found that the ratio of minimum temperature to effective temperature is approximately 0.77 for Arcturus, Procyon, and the sun, and that mass tends to increase at the temperature minimum with decreasing gravity. The pressure is found to be about 1 percent of the solar value, and the surface brightness of the Arcturus transition region and coronal spectrum is estimated to be much less than for the sun. The partial-redistribution calculation for the Ca II K line indicates that the emission width is at least partially determined by damping rather than Doppler broadening, suggesting a reexamination of previous explanations for the Wilson-Bappu effect.

Ayres, T. R.↗

A numerical method for computing three dimensional viscous supersonic flow fields about slender bodies

A numerical method is presented for the calculation of steady, three-dimensional, viscous, compressible flow fields about slender bodies at angle of attack and at supersonic speeds. Approximations are introduced in modeling the flow in the longitudinal direction. Accordingly, the flow fields calculated with the program were computed with a model that permits viscous crossflow together with inviscid axial flow. An analysis of the errors introduced by such a treatment is presented. Numerical calculations were made and compared with experimental results for an ogive-cylinder and an airplane fuselage configuration. Generally, good agreement with experiment was obtained. However, boundary layer separation and body vortex positions differed from experimental locations on the ogive-cylinder, and the shock induced by the fuselage canopy was predicted at a slightly different location.

Walitt, L.↗