Search NASA⌕ Search

SEARCH · Search NASA

Results for “Matrix Theory”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

On the Stability of Power Transmission Systems Under Persistent Inverter Attacks: A Bi-Linear Matrix Approach

We investigate the stability and robustness properties of a power transmission system under persistent deceiving attacks on inverter-interfaced energy resources. The attacks can corrupt the damping coefficients in the inverters' controllers and measurements of the frequency at the points of coupling. Leveraging tools from hybrid dynamical systems theory, we characterize a broad family of persistent (and not necessarily periodic) attacks acting on the inverters, under which the stability properties of the transmission system can be shown to not be compromised. To address potentially conservative conditions identified through conventional bounding techniques, sufficient conditions on the average activation time of the attacks are identified via Lyapunov theory, as well as the formulation and solution of a class of bilinear matrix inequalities (BMI). The results are obtained for constant and slowly time-varying loads via input-to-state stability (ISS) tools. Numerical simulations on the IEEE 39-bus test system are also presented.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1994

Micromechanical Modeling of the Thermal Expansion of Graphite/copper Composites with Nonuniform Microstructure

Two micromechanical models were developed to investigate the thermal expansion of graphite/copper (Gr/Cu) composites. The models incorporate the effects of temperature-dependent material properties, matrix inelasticity, initial residual stresses due to processing history, and nonuniform fiber distribution. The first model is based on the multiple concentric cylinder geometry, with each cylinder treated as a two-phase composite with a characteristic fiber volume fractions. By altering the fiber volume fraction of the individual cylinders, unidirectional composites with radially nonuniform fiber distributions can be investigated using this model. The second model is based on the inelastic lamination theory. By varying the fiber content in the individual laminae, composites with nonuniform fiber distribution in the thickness direction can be investigated. In both models, the properties of the individual regions (cylinders or laminae) are calculated using the method of cells micromechanical model. Classical incremental plasticity theory is used to model the inelastic response of the copper matrix at the microlevel. The models were used to characterize the effects of nonuniform fiber distribution on the thermal expansion of Gr/Cu. These effects were compared to the effects of matrix plasticity, choice of stress-free temperature, and slight fiber misalignment. It was found that the radially nonuniform fiber distribution has little effect on the thermal expansion of Gr/Cu but could become significant for composites with large fiber-matrix transverse CTE and Young's modulus mismatch. The effect of nonuniform fiber distribution in the through-thickness direction of a laminate was more significant, but only approached that of the stress-free temperature for the most extreme cases that include large amounts of bending. Subsequent comparison with experimental thermal expansion data indicated the need for more accurate characterization of the graphite fiber thermomechanical properties. Correlation with cyclic data revealed the presence of a mechanism not considered in the developed models. The predicted response did, however, exhibit ratcheting behavior that has been observed experimentally in Gr/Cu. Finally, simulation of the actual fiber distribution of particular specimens had little effect on the predicted thermal expansion.

Bednarcyk, Brett A.↗
NASA NTRS1988

An economical semi-analytical orbit theory for micro-computer applications

An economical algorithm is presented for predicting the position of a satellite perturbed by drag and zonal harmonics J sub 2 through J sub 4. Simplicity being of the essence, drag is modeled as a secular decay rate in the semi-axis (retarded motion); with the zonal perturbations modeled from a modified version of the Brouwers formulas. The algorithm is developed as: an alternative on-board orbit predictor; a back up propagator requiring low energy consumption; or a ground based propagator for microcomputer applications (e.g., at the foot of an antenna). An O(J sub 2) secular retarded state partial matrix (matrizant) is also given to employ with state estimation. The theory was implemented in BASIC on an inexpensive microcomputer, the program occupying under 8K bytes of memory. Simulated trajectory data and real tracking data are employed to illustrate the theory's ability to accurately accommodate oblateness and drag effects.

Gordon, R. A.↗
NASA NTRS1974

Micromechanics strength theories

The available micromechanics theories for predicting the uniaxial strengths of unidirectional composites are described. The theories presuppose that the required properties of the constituents of the composites are known and that the uniaxial loads are coincident with the principal axes of the unidirectional composite. The uniaxial strengths considered are those required to resist either in-plane loads, flexural (bending), or horizontal-beam shear. The discussion deals primarily with fiber-nonmetallic matrix composites where all the fibers are parallel and lie in the same plane. However, the micromechanics theories presented here are transferable to fiber-metallic matrix composites so long as the underlying assumptions are not violated. Some of the concepts described herein may also be applicable to particulate composites.

Chamis, C. C.↗
DOE OSTI2026

Glauber-theory analysis of nuclear reactions on a 12 C target with variational Monte Carlo wave functions

The application of Glauber theory has been playing an increasingly important role with the study of unstable or exotic nuclei. Its adaptation to medium and high-energy nucleus-nucleus collisions is severely limited because one has to evaluate the matrix elements of multiple-scattering operators. The extraction of physical observables has been done using ‘approximate’ Glauber theory whose validity is hard to evaluate. Here, we perform a full calculation of the matrix elements using Monte Carlo integration and analyze the elastic differential cross sections and the total reaction cross sections for p+¹²C, ⁴,⁶He+¹²C, and ¹²C+¹²C collisions. We use the variational Monte Carlo wave functions for ⁴,⁶He and ¹²C obtained by using realistic two- and three-nucleon potentials. We demonstrate the performance of the Glauber-theory calculations by comparing with available experimental data. We further discuss the accuracy of the conventional approximate methods in the light of the cumulant expansion for Glauber’s phase-shift function.

Horiuchi, W. [Osaka Metropolitan University (Japan↗
NASA NTRS1973

Polarization phenomena in multiphoton ionization of atoms

The theory of multiphoton ionization for an atomic system of arbitrary complexity is developed using a density matrix formalism. An expression is obtained which determines the differential N-photon ionization cross section as a function of the polarization states of the target atom and the incident radiation. The parameters which characterize the photoelectron angular distribution are related to the general reduced matrix elements for the N-photon transition. Two-photon ionization of unpolarized atoms is treated as an illustration of the use of the theory. The dependence of the multiphoton ionization cross section on the polarization state of the incident radiation, which has been observed in two- and three-photon ionization of Cs, is accounted for by the theory. Finally, the photoelectron spin polarization produced by the multiphoton ionization of unpolarized atoms, like the analogous polarization resulting from single-photon ionization, is found to depend on the circular polarization of the incident radiation.

Jacobs, V. L.↗
NASA NTRS1973

Polarization phenomena in multiphoton ionization of atoms.

The theory of multiphoton ionization for an atomic system of arbitrary complexity is developed using a density matrix formalism. An expression is obtained which determines the differential N-photon ionization cross section as a function of the polarization states of the target atom and the incident radiation. The parameters which characterize the photo-electron angular distribution are related to the general reduced matrix elements for the N-photon transition. Two-photon ionization of unpolarized atoms is treated as an illustration of the use of the theory. The dependence of the multiphoton ionization cross section on the polarization state of the incident radiation, which has been observed in two- and three-photon ionization of Cs, is accounted for by the theory. Finally, the photoelectron spin polarization produced by the multiphoton ionization of unpolarized atoms, like the analogous polarization resulting from single-photon ionization, is found to depend on the circular polarization of the incident radiation.

Jacobs, V. L.↗
NASA NTRS1990

A Mechanical Model for Elastic Fiber Microbuckling

A two-dimensional mechanical model is presented to predict the compressive strength of unidirectional fiber composites using technical beam theory and classical elasticity. First, a single fiber resting on a matrix half-plane is considered. Next, a more elaborate analysis of a uniformly laminated, unidirectional fiber composite half-plane is presented. The model configuration incorporates a free edge which introduces a buckling mode that originates at the free edge and decays into the inferior of the half-plane. It is demonstrated that for composites of low volume fraction (less than 0.3), this decay mode furnishes values of buckling strain that are below the values predicted by the Rosen (1965) model. At a higher volume fraction the buckling mode corresponds to a half wavelength that is in violation of the usual assumptions of beam theory. Causes for deviations of the model prediction froi?i existing experimental results are discussed.

Waas, M. A.↗
NASA NTRS1990

A mechanical model for elastic fiber microbuckling

A two-dimensional mechanical model is presented to predict the compressive strength of unidirectional fiber composites using technical beam theory and classical elasticity. First, a single fiber resting on a matrix half-plane is considered. Next, a more elaborate analysis of a uniformly laminated, unidirectional fiber composite half-plane is presented. The model configuration incorporates a free edge which introduces a buckling mode that originates at the free edge and decays into the interior of the half-plane. It is demonstrated that for composites of low volume fraction (less than 0.3), this decay mode furnishes values of buckling strain that are below the values predicted by the Rosen (1965) model. At a higher volume fraction the buckling mode corresponds to a half-wavelength that is in violation of the usual assumptions of beam theory. Causes for deviations of the model prediction from existing experimental results are discussed.

Waas, A. M.↗
DOE OSTI2024

Position-space renormalization schemes for four-quark operators in HQET

Abstract X-space schemes are gauge-invariant, regulator-independent renormalization schemes that are defined by requiring position-space correlation functions of gauge-invariant operators to be equal to their noninteracting values at particular kinematic points. These schemes can be used to nonperturbatively renormalize composite operators in Lattice Quantum Chromodynamics (LQCD), and by computing matching coefficients between theX-space scheme and$$ \overline{\textrm{MS}} $$ MS ¯ in the dimensionally-regulated continuum, matrix elements calculated with LQCD can be converted to$$ \overline{\textrm{MS}} $$ MS ¯ -renormalized matrix elements. UsingX-space schemes for Heavy Quark Effective Theory (HQET) operators has the additional benefit that appropriate ratios of position-space correlation functions cancel the power-divergent static-quark self-energy of Lattice HQET nonperturbatively. This work presents theO(α S ) matching coefficients betweenX-space renormalized four-quark flavor-nonsinglet HQET operators relevant for the lifetimes of charm- and bottom-hadrons, and four-quark HQET operators relevant for mixing between neutral mesons containing a heavy quark, such asB−$$ \overline{B} $$ B ¯ mixing.

Physics↗
DOE OSTI2025

Strong Correlation DMRG and DFT

This project developed new ways to improve computer simulations of materials where electrons interact strongly with each other, a challenge for today’s most widely used method, density functional theory (DFT). We used an exact numerical method, the density matrix renormalization group (DMRG), to create highly accurate reference results for simple model systems, and used these to test DFT, prove when it will converge, and even train machine-learned functionals. We also invented new kinds of localized basis functions (“gausslets” and “multi-sliced gausslets”) and a “sliced-basis” approach that make high-accuracy simulations faster and more practical. These methods were applied to extended hydrogen systems, enabling the direct derivation of accurate low-energy models from first-principles calculations. We also introduced a new formalism, Conditional-Probability DFT, which could bypass traditional approximations. The tools and results from this work, including open-source software releases, will help scientists design and understand complex quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1976

Light scattering by tenuous particles - A generalization of the Rayleigh-Gans-Rocard approach

We consider scattering by arbitrarily shaped particles that satisfy two conditions: (1) that the polarizability of the particle relative to the ambient medium be small compared to 1 and (2) that the phase shift introduced by the particle be less than 2. We solve the integro-differential equation proposed by Shifrin by using the method of successive iterations and then applying a Fourier transform. For the second iteration, results are presented that accurately describe scattering by a broad class of particles. The phase function and other elements of the scattering matrix are shown to be in excellent agreement with Mie theory for spherical scatterers.

Acquista, C.↗
NASA NTRS1977

Efficient numerical treatment of periodic systems with application to stability problems

Two efficient numerical methods for dealing with the stability of linear periodic systems are presented. Both methods combine the use of multivariable Floquet-Liapunov theory with an efficient numerical scheme for computing the transition matrix at the end of one period. The numerical properties of these methods are illustrated by applying them to the simple parametric excitation problem of a fixed end column. The practical value of these methods is shown by applying them to some helicopter rotor blade aeroelastic and structural dynamics problems. It is concluded that these methods are numerically efficient, general and practical for dealing with the stability of large periodic systems.

Friedmann, P.↗
NASA NTRS1982

On the inversion of block tridiagonals without storage constraints

A strategy was developed to permit trade-offs between the number of floating point operations required and the storage requirements for the solution of certain difference problems, such as block tridiagonal systems of equations. This is done by recomputing some intermediate results instead of storing them. Reducing the storage to the square root of the current requirement roughly doubles the number of computations. Reducing the storage more than this tends to make the number of computations prohibitively large. In theory, if m is the order of each sub-matrix in the block tridiagonal matrix, one can solve any linear system with only 5m(2) + 1 temporary storage cells. In many cases m is a constant and quite small. For example, in solving a factored form of the three-dimensional Navier-Stokes equations, the size m of the block tridiagonals is 5. This method lends itself to efficient use on computers with parallel processing or vector processing architectures. On these computers the larger number of floating point operations is more than offset by the decrease in I/O and the increased percentage of vector operations made possible by this algorithm.

Merriam, M. L.↗
NASA NTRS1985

Aircraft cabin noise prediction and optimization

Theoretical and experimental studies were conducted to determine the noise transmission into acoustic enclosures ranging from simple rectangular box models to full scale light aircraft in flight. The structural models include simple, stiffened, curved stiffened, and orthotropic panels and double wall windows. The theoretical solutions were obtained by model analysis. Transfer matrix and finite element procedures were utilized. Good agreement between theory and experiment has been achieved. An efficient acoustic add-on treatment was developed for interior noise control in a twin engine light aircraft.

Vaicaitis, R.↗
NASA NTRS1985

On the factorization of block-tridiagonals without storage constraints

In many programs solving difference equations, problem size is restricted by the number of available memory cells. A strategy has been developed to permit trade-offs between the number of floating point operations required and storage requirements for the solution of certain problems such as block tridiagonal systems of equations. This is done by recomputing some intermediate results instead of storing them. Reducing the storage to the square root of the current requirement will roughly double the number of computations. In theory, if m is the order of each sub-matrix in the block tridiagonal matrix, one can solve any linear system with only 5 sq m + 1 temporary storage cells. This method lends itself to efficient use on computers with parallel processing or vector processing architectures. On these computers the larger number of floating point operations is more than offset by the decrease in I/O and the increased percentage of vector operations made possible by this algorithm.

Merriam, M. L.↗
NASA NTRS2002

Three Dimensional CFD Analysis of the GTX Combustor

The annular combustor geometry of a combined-cycle engine has been analyzed with three-dimensional computational fluid dynamics. Both subsonic combustion and supersonic combustion flowfields have been simulated. The subsonic combustion analysis was executed in conjunction with a direct-connect test rig. Two cold-flow and one hot-flow results are presented. The simulations compare favorably with the test data for the two cold flow calculations; the hot-flow data was not yet available. The hot-flow simulation indicates that the conventional ejector-ramjet cycle would not provide adequate mixing at the conditions tested. The supersonic combustion ramjet flowfield was simulated with frozen chemistry model. A five-parameter test matrix was specified, according to statistical design-of-experiments theory. Twenty-seven separate simulations were used to assemble surrogate models for combustor mixing efficiency and total pressure recovery. ScramJet injector design parameters (injector angle, location, and fuel split) as well as mission variables (total fuel massflow and freestream Mach number) were included in the analysis. A promising injector design has been identified that provides good mixing characteristics with low total pressure losses. The surrogate models can be used to develop performance maps of different injector designs. Several complex three-way variable interactions appear within the dataset that are not adequately resolved with the current statistical analysis.

Steffen, C. J., Jr.↗