Search NASA⌕ Search

SEARCH · Search NASA

Results for “Vectorized algorithm”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Acceleration of convergence of vector sequences

A general approach to the construction of convergence acceleration methods for vector sequence is proposed. Using this approach, one can generate some known methods, such as the minimal polynomial extrapolation, the reduced rank extrapolation, and the topological epsilon algorithm, and also some new ones. Some of the new methods are easier to implement than the known methods and are observed to have similar numerical properties. The convergence analysis of these new methods is carried out, and it is shown that they are especially suitable for accelerating the convergence of vector sequences that are obtained when one solves linear systems of equations iterative. A stability analysis is also given, and numerical examples are provided. The convergence and stability properties of the topological epsilon algorithm are likewise given.

Sidi, A.↗
DOE OSTI2025

Ground state energy and magnetization curve of a frustrated magnetic system from real-time evolution on a digital quantum processor

Models of interacting many-body quantum systems that may realize new exotic phases of matter, notably quantum spin liquids, are challenging to study using even state-of-the-art classical methods such as tensor network simulations. Quantum computing provides a promising route for overcoming these difficulties to find ground states, dynamics, and more. In this paper, we argue that recently developed hybrid quantum-classical algorithms based on real-time evolution are promising methods for solving a particularly important model in the search for spin liquids, the antiferromagnetic Heisenberg model on the two-dimensional kagome lattice. We show how to construct efficient quantum circuits to implement time evolution for the model and to evaluate key observables on the quantum computer, and we argue that the method has favorable scaling with increasing system size. We then restrict to a 12-spin star plaquette from the kagome lattice and a related 8-spin system, and we give an empirical demonstration on these small systems that the hybrid algorithms can efficiently find the ground state energy and the magnetization curve. For these demonstrations, we use four levels of approximation: exact state vectors, exact state vectors with statistical noise from sampling, noisy classical emulators, and (for the 8-spin system only) real quantum hardware, specifically the Quantinuum H1-1 processor; for the noisy simulations and hardware demonstration, we also employ error mitigation strategies based on the symmetries of the Hamiltonian. Our results strongly suggest that these hybrid algorithms present a promising direction for studying quantum spin liquids and more generally for resolving important unsolved problems in condensed matter theory and beyond.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Ab-initio simulation of spin-vibronic spectra of methoxy radical

Despite the fact that experimental and theoretical work on the spectrum of methoxy has stretched from the microwave to the ultraviolet and proceeded for nearly 50 years, parts of the spectrum have remained a challenge to simulate theoretically and make reliable line-by-line assignments. The spectral complexity arises because the radical has a non-zero electron spin and significant vibronic coupling between the two elec- tronic components of the ground state due to the presence of a conical intersection. This work describes a completely ab initio effort to understand and assign the spin- vibronic levels of the X 2E state from 0 to above 3000 cm−1, a region that includes the fundamental transitions of the C-H symmetric and asymmetric stretches that have not previously been identified uniquely. A potential energy surface for methoxy was calculated at the EOM-CCSDT/ANO1 level of theory. Subsequently this potential energy surface was fit to a quartic power series expansion of all nine vibrational nor- mal coordinates (as determined at the minimum of the conical intersection) by the use of a machine-learning-based algorithm. After the addition of spin-orbit coupling, the spin-vibronic problem was solved using both the Krylov-Schur and Lanczos algorithms with the SOCJT3 software to converge eigenvalues up to 3500 cm−1 and their eigen- vectors. The latter were used, in conjunction with the calculated dipole moment and its derivatives (calculated using finite differences at EOM-CCSDT/ANO1 level), to determine spectral intensities for the spin-vibronic spectra. The calculated transition frequencies and intensities were used to simulate and assign the observed transitions of the spin-vibronic spectra of the radical. The credibility of the assignments and their significance is discussed in detail.

Sharma, Ketan [University of Florida, Gainesville,↗
DOE OSTI2025

Chiral rank-$k$ truncations for the multigrid preconditioner of Wilson fermions in lattice QCD

We present a modification to the setup algorithm for the multigrid preconditioner of Wilson fermions in lattice QCD. A larger number of test vectors than that used in conventional multigrid is generated by the smoother. This set of test vectors is then truncated by a singular value decomposition on the chiral components of the test vectors, which are subsequently used to form the prolongation and restriction matrices of the multigrid hierarchy. This modification is demonstrated to improve the convergence of linear equations on an anisotropic lattice with 𝑚𝜋 ≈ 280 MeV from the Hadron Spectrum Collaboration and an isotropic lattice with 𝑚𝜋 ≈ 220 MeV from the MILC Collaboration. The lattice volume dependence of the method is also examined.

Whyte, Travis [Jülich Supercomputing Center, Jülic↗
NASA NTRS1978

Updating of ADAPT predictions of sunspot activity

An eigenvector analysis procedure was used to analyze and develop algorithms for predicting sunspot numbers. The predictors in these algorithms consist of sunspot numbers from the preceding two solar cycles and magnetic index data from the preceding cycle. Predictions are presented for cycles 21 and 22. The sunspot activity for cycle 21 is predicted to remain below 100 until early 1980 when it will rapidly reach a peak of approximately 120. The two sigma accuracy on these estimates is approximately 20 sunspot numbers in the region of the peak and 10 sunspot numbers early and late in the cycle. Algorithms were also developed for predicting the period of future sunspot cycles using the same predictor vector.

Hunter, H. C.↗
NASA NTRS1982

Survey of new vector computers: The CRAY 1S from CRAY research; the CYBER 205 from CDC and the parallel computer from ICL - architecture and programming

Problems which can arise with vector and parallel computers are discussed in a user oriented context. Emphasis is placed on the algorithms used and the programming techniques adopted. Three recently developed supercomputers are examined and typical application examples are given in CRAY FORTRAN, CYBER 205 FORTRAN and DAP (distributed array processor) FORTRAN. The systems performance is compared. The addition of parts of two N x N arrays is considered. The influence of the architecture on the algorithms and programming language is demonstrated. Numerical analysis of magnetohydrodynamic differential equations by an explicit difference method is illustrated, showing very good results for all three systems. The prognosis for supercomputer development is assessed.

Gentzsch, W.↗
NASA NTRS1985

Recent Developments in Nonlinear Normal Mode Initialization

The importance of a balanced initial condition upon GLAS GCM forecasts and assimilation cycles was assessed. An effort to combine previous work on normal mode initialization at GLA is underway to develop an initialization process for the production version of the GLAS 4th order GCM. The major aspects of this work fall into two parts: vectorization of the linear projector code and the insertion of the mode projector and Machenhauer iteration algorithm into the full GLAS GCM. Memory and paging constraints place restrictions on the number of horizontal modes stored for initialization purposes, and on the manner in which they are stored. Only the first five vertical structures of the gravity modes are used. Differing phase and normalization conventions provided many elusive coding errors. A Machenhauer nonlinear normal mode initialization technique is used. This method entails the insertion of a modified version of the mode projector into the full GCM, and the modification of the GCM to allow for iterative calls to the projector.

Bloom, S. C.↗
NASA NTRS1989

A fully-coupled implicit method for thermo-chemical nonequilibrium air at sub-orbital flight speeds

A CFD technique is described in which the finite-rate chemistry in thermal and chemical nonequilibrium air is fully and implicitly coupled with the fluid motion. Developed for use in the suborbital hypersonic flight speed range, the method accounts for nonequilibrium vibrational and electronic excitation and dissociation, but not ionization. The steady-state solution to the resulting system of equations is obtained by using a lower-upper factorization and symmetric Gauss-Seidel sweeping technique through Newton iteration. Inversion of the left-hand-side matrices is replaced by scalar multiplications through the use of the diagonal dominance algorithm. The code, named CENS2H (Compressible-Euler-Navier-Stokes Two-Dimensional Hypersonic), is fully vectorized and requires about 8.8 x 10 to the -5th sec per node point per iteration using a Cray X-MP computer. Converged solutions are obtained after about 2400 iterations. Sample calculations are made for a circular cylinder and a 10 percent airfoil at 5 deg angle of attack. The calculated cylinder flow field agrees with that obtained experimentally. The code predicts a 10 percent change in lift, drag, and pitching moment for the airfoil due to the thermochemical phenomena.

Park, Chul↗
NASA NTRS1990

Transonic Navier-Stokes solutions about a complex high-speed accelerator configuration

Three-dimensional transonic viscous flow computations are presented for a generic high-speed accelerator model which includes wing, body, fillets, and a no-flow through engine nacelle. Solutions are obtained from an algorithm for the compressible Navier-Stokes equations which incorporates an upwind-biased, flux-vector-splitting approach along with longitudinally-patched grids. Results are presented for fully turbulent flow assumptions and include correlations with wind tunnel data. A good quantitative agreement for the forebody surface pressure distribution is achieved between computations and the available wind-tunnel measurements at freestream Mach number equal to 0.9. Furthermore, it is demonstrated that the flow is stagnating around the boattail region due to separation from the aft-engine cowl lip.

Ghaffari, Farhad↗
NASA NTRS1993

Development of iterative techniques for the solution of unsteady compressible viscous flows

During the past two decades, there has been significant progress in the field of numerical simulation of unsteady compressible viscous flows. At present, a variety of solution techniques exist such as the transonic small disturbance analyses (TSD), transonic full potential equation-based methods, unsteady Euler solvers, and unsteady Navier-Stokes solvers. These advances have been made possible by developments in three areas: (1) improved numerical algorithms; (2) automation of body-fitted grid generation schemes; and (3) advanced computer architectures with vector processing and massively parallel processing features. In this work, the GMRES scheme has been considered as a candidate for acceleration of a Newton iteration time marching scheme for unsteady 2-D and 3-D compressible viscous flow calculation; from preliminary calculations, this will provide up to a 65 percent reduction in the computer time requirements over the existing class of explicit and implicit time marching schemes. The proposed method has ben tested on structured grids, but is flexible enough for extension to unstructured grids. The described scheme has been tested only on the current generation of vector processor architecture of the Cray Y/MP class, but should be suitable for adaptation to massively parallel machines.

Hixon, Duane↗
NASA NTRS1992

The explicit computation of integration algorithms and first integrals for ordinary differential equations with polynomials coefficients using trees

This note is concerned with the explicit symbolic computation of expressions involving differential operators and their actions on functions. The derivation of specialized numerical algorithms, the explicit symbolic computation of integrals of motion, and the explicit computation of normal forms for nonlinear systems all require such computations. More precisely, if R = k(x(sub 1),...,x(sub N)), where k = R or C, F denotes a differential operator with coefficients from R, and g member of R, we describe data structures and algorithms for efficiently computing g. The basic idea is to impose a multiplicative structure on the vector space with basis the set of finite rooted trees and whose nodes are labeled with the coefficients of the differential operators. Cancellations of two trees with r + 1 nodes translates into cancellation of O(N(exp r)) expressions involving the coefficient functions and their derivatives.

Crouch, P. E.↗
NASA NTRS1997

Propulsion Flight Research at NASA Dryden From 1967 to 1997

From 1967 to 1997, pioneering propulsion flight research activities have been conceived and conducted at the NASA Dryden Flight Research Center. Many of these programs have been flown jointly with the United States Department of Defense, industry, or the Federal Aviation Administration. Propulsion research has been conducted on the XB-70, F-111 A, F-111E, YF-12, JetStar, B-720, MD-11, F-15, F- 104, Highly Maneuverable Aircraft Technology, F-14, F/A-18, SR-71, and the hypersonic X-15 airplanes. Research studies have included inlet dynamics and control, in-flight thrust computation, integrated propulsion controls, inlet and boattail drag, wind tunnel-to-flight comparisons, digital engine controls, advanced engine control optimization algorithms, acoustics, antimisting kerosene, in-flight lift and drag, throttle response criteria, and thrust-vectoring vanes. A computer-controlled thrust system has been developed to land the F-15 and MD-11 airplanes without using any of the normal flight controls. An F-15 airplane has flown tests of axisymmetric thrust-vectoring nozzles. A linear aerospike rocket experiment has been developed and tested on the SR-71 airplane. This paper discusses some of the more unique flight programs, the results, lessons learned, and their impact on current technology.

Burcham, Frank W., Jr.↗
NASA NTRS1998

An Operational Concept for Flying FMS Trajectories in Center and TRACON Airspace

Current Flight Management Systems (FMS) do a good job of constructing and flying an optimal trajectory for a single aircraft. Unfortunately, flight crews are often unable to fly these FMS routes during arrivals at busy airports. The Center TRACON Automation System (CTAS) has been designed to aid Center and TRACON (Terminal Radar Approach Control) controllers in assigning runways, sequencing and vectoring all classes of aircraft. CTAS bases its advisories on trajectory predictions for arriving aircraft using algorithms very similar to those in airborne FMS systems. This paper presents near and far term operational concepts for how a ground ATM (air traffic management) automation system like CTAS could work more effectively with the airborne automation in FMS equipped aircraft. The concepts for a more compatible air-ground system include: 1) a common route databases for both CTAS and FMS; 2) datalink to downlink information on aircraft weight, final approach speed and trajectory intent and to uplink wind information; 3) new FMS functions to allow flight crews to easily update their FMS trajectory to match the trajectory suggested by the ground automation with voice clearances, and 4) in the far term, datalink to downlink user preferred trajectories and to uplink trajectory clearances in the terminal area. The paper analyses some of the human factors issues that may result in allowing aircraft to fly FMS routes during enroute descent and in the terminal area. A series of linked human in the loop flight deck and air traffic control simulations and a field test with the NASA 757 are being conducted at NASA's Ames and Langley Research Centers to address these issues and to evaluate the operational feasibility of these approaches to more efficient flight and increased airport throughput.

Palmer, Everett↗
NASA NTRS1991

Transonic Navier-Stokes Solutions about a Complex High-Speed Accelerator Configuration

Three dimensional transonic viscous flow computations are presented for a generic high-speed accelerator model which includes wing, body, filets,and a no-flow through engine nacelle. solutions are obtained from an algorithm for the compressible Navier-Stokes equations which incorporated an upwind-biased, flux-vector-splitting approach along with longitudinally-patched grids. Results are presented for fully turbulent flow assumptions and include correlations with wind tunnel data. A good quantitative agreement for the forebody surface pressure distribution is achieved between computations and the available wind-tunnel measurements at M(sub infinity) = 0.9. Furthermore, it is demonstrated that the flow is stagnating around the boattail region due to separation from the aft-engine cowl lip.

Ghaffari, F.↗
NASA NTRS2008

Sensing a Changing Chemical Mixture Using an Electronic Nose

A method of using an electronic nose to detect an airborne mixture of known chemical compounds and measure the temporally varying concentrations of the individual compounds is undergoing development. In a typical intended application, the method would be used to monitor the air in an inhabited space (e.g., the interior of a building) for the release of solvents, toxic fumes, and other compounds that are regarded as contaminants. At the present state of development, the method affords a capability for identifying and quantitating one or two compounds that are members of a set of some number (typically of the order of a dozen) known compounds. In principle, the method could be extended to enable monitoring of more than two compounds. An electronic nose consists of an array of sensors, typically made from polymer carbon composites, the electrical resistances of which change upon exposure to a variety of chemicals. By design, each sensor is unique in its responses to these chemicals: some or all of the sensitivities of a given sensor to the various vapors differ from the corresponding sensitivities of other sensors. In general, the responses of the sensors are nonlinear functions of the concentrations of the chemicals. Hence, mathematically, the monitoring problem is to solve the set of time-dependent nonlinear equations for the sensor responses to obtain the time dependent concentrations of individual compounds. In the present developmental method, successive approximations of the solution are generated by a learning algorithm based on independent-component analysis (ICA) an established information theoretic approach for transforming a vector of observed interdependent signals into a set of signals that are as nearly statistically independent as possible.

Duong, Tuan↗
NASA NTRS1991

An Explicit Upwind Algorithm for Solving the Parabolized Navier-Stokes Equations

An explicit, upwind algorithm was developed for the direct (noniterative) integration of the 3-D Parabolized Navier-Stokes (PNS) equations in a generalized coordinate system. The new algorithm uses upwind approximations of the numerical fluxes for the pressure and convection terms obtained by combining flux difference splittings (FDS) formed from the solution of an approximate Riemann (RP). The approximate RP is solved using an extension of the method developed by Roe for steady supersonic flow of an ideal gas. Roe's method is extended for use with the 3-D PNS equations expressed in generalized coordinates and to include Vigneron's technique of splitting the streamwise pressure gradient. The difficulty associated with applying Roe's scheme in the subsonic region is overcome. The second-order upwind differencing of the flux derivatives are obtained by adding FDS to either an original forward or backward differencing of the flux derivative. This approach is used to modify an explicit MacCormack differencing scheme into an upwind differencing scheme. The second order upwind flux approximations, applied with flux limiters, provide a method for numerically capturing shocks without the need for additional artificial damping terms which require adjustment by the user. In addition, a cubic equation is derived for determining Vegneron's pressure splitting coefficient using the updated streamwise flux vector. Decoding the streamwise flux vector with the updated value of Vigneron's pressure splitting improves the stability of the scheme. The new algorithm is applied to 2-D and 3-D supersonic and hypersonic laminar flow test cases. Results are presented for the experimental studies of Holden and of Tracy. In addition, a flow field solution is presented for a generic hypersonic aircraft at a Mach number of 24.5 and angle of attack of 1 degree. The computed results compare well to both experimental data and numerical results from other algorithms. Computational times required for the upwind PNS code are approximately equal to an explicit PNS MacCormack's code and existing implicit PNS solvers.

Korte, John J.↗
NASA NTRS1991

Choice of velocity variables for complex flow computation

The issue of adopting the velocity components as dependent velocity variables for the Navier-Stokes flow computations is investigated. The viewpoint advocated is that a numerical algorithm should preferably honor both the physical conservation law in differential form and the geometric conservation law in discrete form. With the use of Cartesian velocity vector, the momentum equations in curvilinear coordinates can retain the full conservation-law form and satisfy the physical conservation laws. With the curvilinear velocity components, source terms appear in differential equations and hence the full conservation law form can not be retained. In discrete expressions, algorithms based on the Cartesian components can satisfy the geometric conservation-law form for convection terms but not for viscous terms; those based on the curvilinear components, on the other hand, cannot satisfy the geometric conservation-law form for either convection or viscous terms. Several flow solutions for domain with 90 and 360 degree turnings are presented to illustrate the issues of using the Cartesian velocity components and the staggered grid arrangement.

Shyy, W.↗