SEARCH · Search NASA

Results for “Symmetric positive definite”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS1999

Object-Oriented Design for Sparse Direct Solvers

We discuss the object-oriented design of a software package for solving sparse, symmetric systems of equations (positive definite and indefinite) by direct methods. At the highest layers, we decouple data structure classes from algorithmic classes for flexibility. We describe the important structural and algorithmic classes in our design, and discuss the trade-offs we made for high performance. The kernels at the lower layers were optimized by hand. Our results show no performance loss from our object-oriented design, while providing flexibility, case of use, and extensibility over solvers using procedural design.

Dobrian, Florin↗

NASA NTRS1973

On stability of discrete composite systems

Conditions are developed under which exponential stability of a composite discrete system is implied by exponential stability of its subsystems and the nature of their interactions. Stability of the system is determined by testing positive definiteness property of a real symmetric matrix the dimension of which is equal to the number of subsystems.

Grujic, L. T.↗

NASA NTRS1973

On stability of discrete composite systems.

Grujic, L. T.↗

NASA NTRS1985

Application of Ritz vectors for dynamic analysis of large structures

The use of an orthogonal set of specially selected Ritz vectors is shown to be very effective in reducing the cost of dynamic analysis by modal superposition. Several mechanical structures are examined, and the Ritz vector approach is compared to the classical eigenvector approach on the basis of cost, accuracy and elapsed analysis (throughput) time. Mathematical proof of the completeness of orthogonal Ritz vectors is provided for the case of a positive definite mass matrix and a symmetric stiffness matrix.

Arnold, R. R.↗

NASA NTRS2011

Efficient Kriging Algorithms

More efficient versions of an interpolation method, called kriging, have been introduced in order to reduce its traditionally high computational cost. Written in C++, these approaches were tested on both synthetic and real data. Kriging is a best unbiased linear estimator and suitable for interpolation of scattered data points. Kriging has long been used in the geostatistic and mining communities, but is now being researched for use in the image fusion of remotely sensed data. This allows a combination of data from various locations to be used to fill in any missing data from any single location. To arrive at the faster algorithms, sparse SYMMLQ iterative solver, covariance tapering, Fast Multipole Methods (FMM), and nearest neighbor searching techniques were used. These implementations were used when the coefficient matrix in the linear system is symmetric, but not necessarily positive-definite.

Memarsadeghi, Nargess↗

NASA NTRS1986

Preconditioned conjugate residual methods for the solution of spectral equations

Conjugate residual methods for the solution of spectral equations are described. An inexact finite-difference operator is introduced as a preconditioner in the iterative procedures. Application of these techniques is limited to problems for which the symmetric part of the coefficient matrix is positive definite. Although the spectral equation is a very ill-conditioned and full matrix problem, the computational effort of the present iterative methods for solving such a system is comparable to that for the sparse matrix equations obtained from the application of either finite-difference or finite-element methods to the same problems. Numerical experiments are shown for a self-adjoint elliptic partial differential equation with Dirichlet boundary conditions, and comparison with other solution procedures for spectral equations is presented.

Wong, Y. S.↗

DOE OSTI2022

Solving generalized eigenvalue problems by ordinary differential equations on a quantum computer

Many eigenvalue problems arising in practice are often of the generalized form A x = λ B x . One particularly important case is symmetric, namely A , B are Hermitian and B is positive definite. The standard algorithm for solving this class of eigenvalue problems is to reduce them to Hermitian eigenvalue problems. For a quantum computer, quantum phase estimation is a useful technique to solve Hermitian eigenvalue problems. In this work, we propose a new quantum algorithm for symmetric generalized eigenvalue problems using ordinary differential equations. The algorithm has lower complexity than the standard one based on quantum phase estimation. Moreover, it works for a wider case than symmetric: B is invertible, B - 1 A is diagonalizable and all the eigenvalues are real.

Science & Technology - Other Topics↗

NASA NTRS1973

Eigenproblem solution by a combined Sturm sequence and inverse iteration technique.

Description of an efficient and numerically stable algorithm, along with a complete listing of the associated computer program, developed for the accurate computation of specified roots and associated vectors of the eigenvalue problem Aq = lambda Bq with band symmetric A and B, B being also positive-definite. The desired roots are first isolated by the Sturm sequence procedure; then a special variant of the inverse iteration technique is applied for the individual determination of each root along with its vector. The algorithm fully exploits the banded form of relevant matrices, and the associated program written in FORTRAN V for the JPL UNIVAC 1108 computer proves to be most significantly economical in comparison to similar existing procedures. The program may be conveniently utilized for the efficient solution of practical engineering problems, involving free vibration and buckling analysis of structures. Results of such analyses are presented for representative structures.

Gupta, K. K.↗

NASA NTRS1989

Factorization Of Positive Definite, Banded Hermitian Matrices

Report discusses application of Cholesky factorization algorithm to positive definite, banded Hermitian matrices. Begins by extending Cholesky factorization algorithm to cover uniformly-partitioned, banded, positive definite matrices of rank n that is real symmetric or Hermitian. Then two stratagems given for use of algorithm in concurrent-processing system in which N less than it has to be to enable factorization of matrix in as few serial steps as possible and where uniformly high efficiency expected from all processing elements. One of major purposes of this and related studies to maximize speedup and efficiency in system of concurrent-data-processing elements.

Salama, Moktar A.↗

NASA NTRS1986

Concurrent Cholesky factorization of positive definite banded Hermitian matrices

First, the Cholesky factorization is extended to cover uniformly partitioned banded positive definite matrices of rank n which may be real symmetric or Hermitian. Then, two stratagems are given for the use of the algorithm in concurrent machines where the number of processing elements is less than required to factor the matrix in as few serial steps as possible, and where uniformly high efficiency is expected from all processing elements. Expressions are given for the efficiency factor e appearing in the speed-up expression q = eN, and these are specialized for the N node hypercube machine as a function of partition size s, the number N of processing elements of the hypercube machine, and the cost mu of interelement transmission relative to computation. It is shown that the efficiency factor e is inversely proportional to mu/s, and that e is almost independent of N when N is large and mu/s = 0. The task is completed in n/s serial steps with no limit on n. The half bandwidth b of the matrix is 2 Ns.

Utku, S.↗

NASA NTRS1996

Schwarz Methods: To Symmetrize or not to Symmetrize

A preconditioning theory for Schwarz methods is presented. The theory establishes sufficient conditions for multiplicative and additive Schwarz algorithms to yield self-adjoint positive definite preconditioners. It allows for the analysis and use of non-variational and non-convergent linear methods as preconditioners for conjugate gradient methods, and it is applied to domain decomposition and multigrid. This paper illustrates why symmetrizing may be a bad idea for linear methods. Numerical examples are presented for a test problem.

Holst, Michael↗

NASA NTRS1999

Flight Test of an Adaptive Configuration Optimization System for Transport Aircraft

A NASA Dryden Flight Research Center program explores the practical application of real-time adaptive configuration optimization for enhanced transport performance on an L-1011 aircraft. This approach is based on calculation of incremental drag from forced-response, symmetric, outboard aileron maneuvers. In real-time operation, the symmetric outboard aileron deflection is directly optimized, and the horizontal stabilator and angle of attack are indirectly optimized. A flight experiment has been conducted from an onboard research engineering test station, and flight research results are presented herein. The optimization system has demonstrated the capability of determining the minimum drag configuration of the aircraft in real time. The drag-minimization algorithm is capable of identifying drag to approximately a one-drag-count level. Optimizing the symmetric outboard aileron position realizes a drag reduction of 2-3 drag counts (approximately 1 percent). Algorithm analysis of maneuvers indicate that two-sided raised-cosine maneuvers improve definition of the symmetric outboard aileron drag effect, thereby improving analysis results and consistency. Ramp maneuvers provide a more even distribution of data collection as a function of excitation deflection than raised-cosine maneuvers provide. A commercial operational system would require airdata calculations and normal output of current inertial navigation systems; engine pressure ratio measurements would be optional.

Gilyard, Glenn B.↗

NASA NTRS1992

Design of dissipative low-authority controllers using an eigensystem assignment technique

A novel method for the design of dissipative, low-authority controllers has been developed. The method uses a sequential approach along with eigensystem assignment to compute rate and position gain matrices that assign a number of closed-loop poles of the system to desired locations. Because the feedback gain matrices are symmetric and nonnegative definite, the closed-loop stability is always guaranteed regardless of the model order or parameter inaccuracies. The resulting (nominal) closed-loop system can have specified damping ratios for m modes, which makes the plant amenable to high-authority controller design, using methods such as LQG/LTR or H-infinity. A numerical example is worked out for a flexible structure in order to demonstrate the proposed technique.

Maghami, P. G.↗

DOE OSTI2021

Two-nucleon 𝑆-wave interactions at the SU(3) flavor-symmetric point with 𝑚 𝑢⁢𝑑 ≈ 𝑚$^{phys}_𝑠$: A first lattice QCD calculation with the stochastic Laplacian Heaviside method

We report on the first application of the stochastic Laplacian Heaviside method for computing multiparticle interactions with lattice QCD to the two-nucleon system. Like the Laplacian Heaviside method, this method allows for the construction of interpolating operators which can be used to construct a set of positive-definite two-nucleon correlation functions, unlike nearly all other applications of lattice QCD to two nucleons in the literature. It also allows for a variational analysis in which optimal linear combinations of the interpolating operators are formed that couple predominantly to the eigenstates of the system. Utilizing such methods has become of paramount importance to help resolve the discrepancy in the literature on whether two nucleons in either isospin channel form a bound state at pion masses heavier than physical, with the discrepancy persisting even in the SU(3)-flavor-symmetric point with all quark masses near the physical strange quark mass. This is the first in a series of papers aimed at resolving this discrepancy. In the present work, we employ the stochastic Laplacian Heaviside method without a hexaquark operator in the basis at a lattice spacing of 𝑎 ≈0.086 fm, lattice volume of 𝐿 = 48⁢𝑎 ≈ 4.1 fm and pion mass 𝑚 𝜋 ≈ 714 MeV. With this setup, the observed spectrum of two-nucleon energy levels strongly disfavors the presence of a bound state in either the deuteron or dineutron channel.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1962

Empire. a study of early manned interplanetary missions final report, may 26 - nov. 25, 1962

This report summarizes the investigations and results of the EMPIRE Study Program undertaken by Aeronutronic Division of Ford Motor Company for the Future Projects Office, Marshall Space Flight Center, under Contract NAS8-5025. The dual planet flyby missions of the Crocco and Symmetric trajectory classes are discussed. The Crocco mission with an August 1971 launch window requires an interplanetary injection velocity increment of i0.i km/sec, has a return velocity of 13.5 km/sec, and takes approximately 400 days. The Symmetric mission with a July 1970 launch window has an injection velocity increment of 5.3 km/sec, a return velocity of 15.8 km/sec, and takes approximately 630 days. Additional results of the trajectory studies and abort trajectories are reported. The guidance and navigation subsystem, midcourse corrections, and planetary approach corrections are discussed. A detailed analysis of the reentry phase of EMPIRE includes consideration of an Apollo-type, a Drag Brake, and a lifting-type reentry vehicle to return the six-man crew at mission completion or in an aborted condition. The High L/D reentry vehicle is used in the missions considered. The various technological areas required for design criteria are developed and several spacecraft designs are considered. The all chemical propulsion Crocco system is discarded due to weight, complexity, and cost. The nuclear injected Crocco is treated in a similar manner. The lower energy injection for the Symmetric Mission leads to the feasibility of a nuclear injected vehicle with an Earth orbit weight of about 180,000 kilograms (400,000 pounds) before interplanetary transit. In addition, two chemical symmetric vehicles are treated. Conservative radiation exposures are derived, for the 630 day mission, of less than 200 REM and a polyethelene radiation shelter is designed. Scientific aspects of the missions are discussed. Mission Success Probabilities are presented for the various missions considered and for Saturn C-5, Nova, and Super-Nova Earth launch vehicles in light of possible development. The need for acceleration of nuclear rocket engine developments and auxiliary power developments is indicated. Definition of a larger nuclear engine of the order of 200,000 pounds thrust and about 800 seconds burning time or 50,000 pound thrust and 3600 seconds burning time is indicated for the Symmetric Mission in 1970 (energy requirements are higher in 1972 and for later launch due to the less favorable position of Mars)_ Immediate development of this advanced nuclear propulsion capability is recommended. A Development Plan and Funding Schedule is given for the 1970 launch window pinpointing the critical development areas and indicating a total program cost of $12.6 billion independent of other programmed R&D costs. In conclusion, technological feasibility for an early manned dual planet Mars-Venus flyby is believed to be demonstrated in this study. Several areas of accelerated development and experimental confirmation of theory are pinpointed. The necessary funding and development of Nova or orbital operations capability with Saturn C-5's is required. The 1970 launch window appears to offer the least expensive Symmetric Mission for several years into the 1980's.

F. P. Dixon↗

DOE OSTI2021

Intrinsic Anomalous Hall Conductivity in a Nonuniform Electric Field

We study how the intrinsic anomalous Hall conductivity is modified in two-dimensional crystals with broken time-reversal symmetry due to weak inhomogeneity of the applied electric field. Focusing on a clean noninteracting two-band system without band crossings, we derive the general expression for the Hall conductivity at small finite wave vector q to order q 2 , which governs the Hall response to the second gradient of the electric field. Using the Kubo formula, we show that the answer can be expressed through the Berry curvature, Fubini-Study quantum metric, and the rank-3 symmetric tensor which is related to the quantum geometric connection and physically corresponds to the gauge-invariant part of the third cumulant of the position operator. We further compare our results with the predictions made within the semiclassical approach. By deriving the semiclassical equations of motion, we reproduce the result obtained from the Kubo formula in some limits. We also find, however, that the conventional semiclassical description in terms of the definite position and momentum of the electron is not fully consistent because of singular terms originating from the Heisenberg uncertainty principle. We thus present a clear example of a case when the semiclassical approach inherently suffers from the uncertainty principle, implying that it should be applied to systems in nonuniform fields with extra care.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2010

Establishing a Robotic, LEO-to-GEO Satellite Servicing Infrastructure as an Economic Foundation for Exploration

The strategy for accomplishing civilian exploration goals and objectives is in the process of a fundamental shift towards a potential new approach called Flexible Path. This paper suggests that a government-industry or public-private partnership in the commercial development of low Earth orbit to geostationary orbit (LEO-to-GEO (LTG)) space, following or in parallel with the commercialization of Earth-to-LEO and International Space Station (ISS) operations, could serve as a necessary, logical step that can be incorporated into the flexible path approach. A LTG satellite-servicing infrastructure and architecture concept is discussed within this new strategic context. The concept consists of a space harbor that serves as a transport facility for a fleet of specialized, fully- or semi-autonomous robotic servicing spacecraft. The baseline, conceptual system architecture is composed of a space harbor equipped with specialized servicer spacecraft; a satellite command, communication, and control system; a parts station; a fuel station or depot; and a fuel/parts replenishment transport. The commercial servicer fleet would consist of several types of spacecraft, each designed with specialized robotic manipulation subsystems to provide services such as refueling, upgrade, repair, inspection, relocation, and removal. The space harbor is conceptualized as an ISS-type, octagonal truss structure equipped with radiation tolerant subsystems. This space harbor would be primarily capable of serving as an operational platform for various commercially owned and operated servicer spacecraft positioned and docked symmetrically on four of the eight sides. Several aspects of this concept are discussed, such as: system-level feasibility in terms of ISS-truss-type infrastructure and subsystems emplacement and maintenance between LEO and GEO; infrastructure components assembly in LEO, derived from ISS assembly experience, and transfer to various higher orbital locations; the evolving Earth-to-orbit (ETO) capability to deliver humans and cargo to LEO for assembly purposes; system architectural definition, optimal orbital parameters, mass estimations, delta velocity ( V) estimations, power and propulsion options, and assessments of various critical technologies. Large-scale, robotic, LTG satellite servicing is considered as an essential economic pre-condition and next parallel or sequential step on the road toward exploration beyond LEO. Such a step might produce the necessary pre-requisite economic value that can be used by future decision makers to justify further investment in exploration beyond LEO.

Horsham, Gary A. P.↗

DOE OSTI2023

On de Sitter spacetime and string theory

We review various aspects of de Sitter spacetime in string theory: its status as an Effective Field Theory spacetime solution, its relation to the vacuum energy problem in string theory, its (global) holographic definition in terms of two entangled and noncanonical conformal field theories as well as a realization of a realistic de Sitter universe endowed with the observed visible matter and the necessary dark sector in order to reproduce the realistic cosmological structure. In particular, based on the new insight regarding the cosmological constant problem in string theory, we argue that in a doubled, [Formula: see text]-duality-symmetric, phase-space-like and noncommutative generalized-geometric formulation, string theory can naturally lead to a small and positive cosmological constant that is radiatively stable and technically natural. Such a formulation is fundamentally based on a quantum spacetime, but in an effective spacetime description of this general formulation of string theory, the curvature of the dual spacetime is the cosmological constant of the observed spacetime, while the size of the dual spacetime is the gravitational constant of the same observed spacetime. Also, the three scales associated with intrinsic noncommutativity of string theory, the cosmological constant scale, the Planck scale as well as the Higgs scale, can be arranged to satisfy various seesaw-like formulae. Along the way, we show that these new features of string theory can be implemented in a particular deformation of cosmic-string-like models.

Astronomy & Astrophysics↗