SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,675 records · Page 93

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

A Machine Learning-Derived Atomistic Potential for Y2Si2O7

Incorporation of SiC/SiC ceramic matrix composite (CMC) hot section components into aircraft engines promises to increase efficiency and safety. However, SiC/SiC CMCs are subject to water vapor-induced oxidation and recession at the high temperatures of engine operation, and thus environmental barrier coatings (EBCs) are required to reduce this degradation and enable their widespread adoption. An understanding of EBCs failure mechanisms, including thermochemical and thermomechanical mechanisms, is essential as coating degradation leads to reduced CMC component service life. Computational modeling approaches can provide insight into EBC material properties important for coating design. However, density functional theory (DFT) is computationally expensive and atomistic potentials are lacking for materials of interest. In this work, we utilize a machine learning approach and DFT training data to parameterize atomistic potentials for two candidate EBC materials, Y2Si2O7 and Yb2Si2O7. These potentials enable near DFT-accurate calculations of thermodynamic and thermomechanical properties essential to EBC design.

Cameron J Bodenschatz↗

DOE OSTI2019

Scattering phase shifts and mixing angles for an arbitrary number of coupled channels on the lattice

We present a lattice method for determining scattering phase shifts and mixing angles for the case of an arbitrary number of coupled channels. Previous nuclear lattice effective field theory simulations were restricted to mixing of up to two partial waves for scattering of two spin-1/2 particles, which is insufficient for analyzing nucleon-nucleus or nucleus-nucleus scattering processes. In the proposed method, the phase shifts and mixing angles are extracted from the radial wave functions obtained by projecting the three-dimensional lattice Hamiltonian onto the partial wave basis. In this work, we use a spherical wall potential as a boundary condition along with a channel-mixing auxiliary potential to construct the full-rank S matrix. Our method can be applied to particles with any spin, but we focus here on scattering of two spin-1 bosons involving up to four coupled channels. For a considered test potential, the phase shifts and mixing angles extracted on the lattice are shown to agree with the ones calculated by solving the Schrödinger equation in the continuum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Communication Lower Bounds and Optimal Algorithms for Symmetric Matrix Computations

In this article, we focus on the communication costs of three symmetric matrix computations: (i) multiplying a matrix with its transpose, known as a symmetric rank-k update (SYRK) (ii) adding the result of the multiplication of a matrix with the transpose of another matrix and the transpose of that result, known as a symmetric rank-2k update (SYR2K) (iii) performing matrix multiplication with a symmetric input matrix (SYMM). All three computations appear in the Level 3 Basic Linear Algebra Subroutines (BLAS) and have wide use in applications involving symmetric matrices. We establish communication lower bounds for these kernels using sequential and distributed-memory parallel computational models, and we show that our bounds are tight by presenting communication-optimal algorithms for each setting. Our lower bound proofs rely on applying a geometric inequality for symmetric computations and analytically solving constrained nonlinear optimization problems. As a result, the symmetric matrix and its corresponding computations are accessed and performed according to a triangular block partitioning scheme in the optimal algorithms.

Al Daas, Hussam [Rutherford Appleton Laboratory, D↗

NASA NTRS1976

Optimal control of a low wing-loading STOL aircraft

Linear optimal quadratic control theory is applied to a low wing-loading STOL aircraft for ride quality and flight path following. Design criteria include minimum rms response to wind turbulence and desired transient response characteristics. Design techniques include proper choosing of design versus evaluation models, choosing appropriate performance index responses, and use of classical evaluation techniques. Results are obtained through a combination of frequency response shaping and gust observation. Effects of control rate and authority saturation are examined with a new rapid calculation of random input describing functions. Parameter sensitivity is also evaluated using a Liapunov type matrix equation.

Cunningham, T. B.↗

NASA NTRS1982

Effects of cobalt in nickel-base superalloys

A study has been carried out to assess the role of cobalt in Udimet 700, a representative nickel-base superalloy containing 17 percent or more cobalt. The study spans the spectrum of microstructural, microchemical, and mechanical behavior aspects which together form a basis for superalloy performance in jet engines. The results suggest that cobalt affects the solubility of elements in the gamma matrix, which leads to enhanced gamma-prime volume fraction and to the stabilization of MC-type carbides and sigma phase. However, these microstructural and microchemical changes are too slight to significantly affect the strength and ductile properties. Depending on the heat treatment, the creep and stress rupture resistance can be cobalt-sensitive. In the coarse-grained, fully solutioned and aged condition, all of the alloy's 17 percent Co can be replaced by nickel without decreasing the creep and stress rupture resistance. These findings are discussed with reference to existing theories and experimental data obtained by other workers.

Tien, J. K.↗

DOE OSTI2024

Revisiting Artifacts of Kohn–Sham Density Functionals for Biosimulation

We revisit the problem of unphysical charge density delocalization/fractionalization induced by the self-interaction error of common approximate Kohn–Sham (KS) density functional theory functionals on simulation of small to medium-sized proteins in a vacuum. Aside from producing unphysical electron densities and total energies, the vanishing of the HOMO–LUMO gap associated with the unphysical charge delocalization leads to an unphysical low-energy spectrum and catastrophic failure of most popular solvers for the KS self-consistent field (SCF) problem. We apply a robust quasi-Newton SCF solver to obtain solutions for some of these difficult cases. The anatomy of the charge delocalization is revealed by the natural deformation orbitals obtained from the density matrix difference between the Hartree–Fock and KS solutions; the charge delocalization not only can occur between charged fragments (such as in zwitterionic polypeptides) but also involves neutral fragments. The vanishing-gap phenomenon and troublesome SCF convergence are both attributed to the unphysical KS Fock operator eigenspectra of molecular fragments (e.g., amino acids or their side chains). Analysis of amino acid pairs suggests that the unphysical charge delocalization can be partially ameliorated by the use of some range-separated hybrid functionals but not by semilocal or standard hybrid functionals. Last, we demonstrate that solutions without the unphysical charge delocalization can be located even for semilocal KS functionals highly prone to such defects, but such solutions have non-Aufbau character and are unstable with respect to mixing of the non-overlapping “frontier” orbitals. Caution should be exercised when unexpectedly small (or vanishing) HOMO–LUMO gaps and atypical SCF convergence patterns (e.g., oscillatory) are observed in KS DFT simulations in any context (bio or otherwise).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Two-color polarization spectroscopy measurements of Zeeman state-to-state collision induced transitions of nitric oxide in binary gas mixtures

We investigated collision induced transitions in the (0, 0) band of the A2Σ+–X2Π electronic transition of nitric oxide (NO) using two-color polarization spectroscopy (TCPS). Two sets of TCPS spectra for 1% NO, diluted in different buffer gases at 295 K and 1 atm, were obtained with the pump beam tuned to the R11(11.5) and OP12(1.5) transitions. The buffer gases were He, Ar, and N2. The probe was scanned while the pump beam was tuned to the line center. Theoretical TCPS spectra, calculated by solving the density matrix formulation of the time-dependent Schrödinger wave equation, were compared with the experimental spectra. A collision model based on the modified exponential-gap law was used to model the rotational level-to-rotational level collision dynamics. A model for collisional transfer from an initial to a final Zeeman state was developed based on the difference in cosine of the rotational quantum number J projection angle with the z-axis for the two Zeeman states. Rotational energy transfer rates and Zeeman state collisional dynamics were varied to obtain good agreement between theory and experiment for the two different TCPS pump transitions and for the three different buffer gases. One key finding, in agreement with quasi-classical trajectory calculations, is that the spin-rotation changing transition rate in the A2Σ+ level of NO is almost zero for rotational quantum numbers ≥8. It was necessary to set this rate to near zero to obtain agreement with the TCPS spectra.

Chemistry↗

NASA NTRS1981

Structure of correlation tensors in homogeneous anisotropic turbulence

The theory of isotropic tensors, developed by Robertson (1940), Batchelor (1946), and Chandrasekhar (1951), is extended to cover the general case of turbulence with a pseudo-vector-preferred direction, without assuming mirror-reflection invariance. Attention is focused on two-point-correlation functions, and it is shown that the form of the decomposition into proper and pseudo-tensor contributions is restricted by the homogeneity requirement. The somewhat unexpected result that the vector- and pseudo-vector-preferred-direction cases yield different results is presented: A pseudo-vector-preferred direction allows the correlation matrix one more functional degree of freedom than does the proper vector case. An explicit form of the two-point-correlation tensor in the presence of a uniform mean magnetic field is presented which may be appropriate for use in analysis of magnetic fluctuations in plasma containment devices or the interplanetary medium. A procedure for determining the magnetic helicity from experimental data is presented.

Matthaeus, W. H.↗

NASA NTRS1991

Integrated mechanics for the passive damping of polymer-matrix composites and composite structures

Some recent developments on integrated damping mechanics for unidirectional composites, laminates, and composite structures are reviewed. Simplified damping micromechanics relate the damping of on-axis and off-axis composites to constituent properties, fiber volume ratio, fiber orientation, temperature, and moisture. Laminate and structural damping mechanics for thin composites are summarized. Discrete layer damping mechanics for thick laminates, including the effects of interlaminar shear damping, are developed and semianalytical predictions of modal damping in thick simply supported specialty composite plates are presented. Applications show the advantages of the unified mechanics, and illustrate the effect of fiber volume ratio, fiber orientation, structural geometry, and temperature on the damping. Additional damping properties for composite plates of various laminations, aspect ratios, fiber content, and temperature illustrate the merits and ranges of applicability of each theory (thin or thick laminates).

Saravanos, D. A.↗

NASA NTRS1981

An optical model for composite nuclear scattering

The optical model of composite particle scattering is considered and compared to the accuracies of other models. A nonrelativistic Schroedinger equation with two-body potentials is used for the scattering of a single particle by an energy-dependent local potential. The potential for the elastic channel is composed of matrix elements of a single scattering operator taken between the ground states of the projectile and the target; the coherent amplitude is considered as dominating the scattering in the forward direction. A multiple scattering series is analytically explored and formally summed by the solution of an equivalent Schroedinger equation. Cross sections of nuclear scattering are then determined for He-4 and C-12 nuclei at 3.6 GeV/nucleus and O-16 projectiles at 2.1 GeV/nucleus, and the optical model approximations are found to be consistently lower and more accurate than approximations made by use of Glauber's theory.

Wilson, J. W.↗

NASA NTRS1995

Galactic oscillator symmetry

Riemann ellipsoids model rotating galaxies when the galactic velocity field is a linear function of the Cartesian coordinates of the galactic masses. In nuclear physics, the kinetic energy in the linear velocity field approximation is known as the collective kinetic energy. But, the linear approximation neglects intrinsic degrees of freedom associated with nonlinear velocity fields. To remove this limitation, the theory of symplectic dynamical symmetry is developed for classical systems. A classical phase space for a self-gravitating symplectic system is a co-adjoint orbit of the noncompact group SP(3,R). The degenerate co-adjoint orbit is the 12 dimensional homogeneous space Sp(3,R)/U(3), where the maximal compact subgroup U(3) is the symmetry group of the harmonic oscillator. The Hamiltonian equations of motion on each orbit form a Lax system X = (X,F), where X and F are elements of the symplectic Lie algebra. The elements of the matrix X are the generators of the symplectic Lie algebra, viz., the one-body collective quadratic functions of the positions and momenta of the galactic masses. The matrix F is composed from the self-gravitating potential energy, the angular velocity, and the hydostatic pressure. Solutions to the hamiltonian dynamical system on Sp(3,R)/U(3) are given by symplectic isospectral deformations. The Casimirs of Sp(3,R), equal to the traces of powers of X, are conserved quantities.

Rosensteel, George↗

NASA NTRSDate not supplied

Unsteady Aerodynamic Modeling of Atmospheric Entry Vehicles in Subsonic and Incompressible flow: A Frequency Response Approach

The determination of time-varying lift force and pitch moment generated by a purely pitching Earth-entry capsule is investigated. Experiments were conducted in the 12-foot Low-Speed Tunnel at NASA Langley Research Center, testing a range of oscillation frequencies at a zero-mean angle of attack with a pitching amplitude of 10 degrees. In light of these measurements, a closed-form set of analytically derived equations for lift and moment was used to develop a semi-empirical formulation, incorporating empirically determined values from the experimental runs. The equations are grounded in potential flow theory, Theodorsen's classical theory of unsteady aerodynamics, and the Joukowski theorem of conformal mapping. The unsteady aerodynamics generated by the oscillating body are then modeled by constructing frequency response functions, with quasi-steady forces and moments serving as inputs and unsteady forces and moments as outputs. The experimentally determined gain and phase variations characterize the unsteady nature of the flow and the system's response and flow time-lag to input flow parameters for a blunt-body entry vehicle. The final semi-empirical model is validated with a set of parameters beyond the initial test matrix.

Entry Systems Modeling↗

NASA NTRSDate not supplied

Unsteady Aerodynamic Modeling of Atmospheric Entry Vehicles in Subsonic and Incompressible Flow: A Frequency Response Approach

Entry Systems Modeling↗

DOE OSTI2021

Isolating the Electrocatalytic Activity of a Confined NiFe Motif within Zirconium Phosphate

Unique classes of active-site motifs are needed for improved electrocatalysis. Herein, the activity of a new catalyst motif is engineered and isolated for the oxygen evolution reaction (OER) created by nickel–iron transition metal electrocatalysts confined within a layered zirconium phosphate matrix. It is found that with optimal intercalation, confined NiFe catalysts have an order of magnitude improved mass activity compared to more conventional surface-adsorbed systems in 0.1 m KOH. Interestingly, the confined environments within the layered structure also stabilize Fe-rich compositions (90%) with exceptional mass activity compared to known Fe-rich OER catalysts. Through controls and by grafting inert molecules to the outer surface, it is evidenced that the intercalated Ni/Fe species stay within the interlayer during catalysis and serve as the active site. After determining a possible structure (wycherproofite), density functional theory is shown to correlate with the observed experimental compositional trends. It is further demonstrated that the improved activity of this motif is correlated to the Fe and water content/composition within the confined space. This work highlights the catalytic enhancement possibilities available through zirconium phosphate and isolates the activity from the intercalated species versus surface/edge ones, thus opening new avenues to develop and understand catalysts within unique nanoscale chemical environments.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Transport phenomena in electrolyte solutions: Nonequilibrium thermodynamics and statistical mechanics

The theory of transport phenomena in multicomponent electrolyte solutions is presented here through the integration of continuum mechanics, electromagnetism, and nonequilibrium thermodynamics. The governing equations of irreversible thermodynamics, including balance laws, Maxwell's equations, internal entropy production, and linear laws relating the thermodynamic forces and fluxes, are derived. Green–Kubo relations for the transport coefficients connecting electrochemical potential gradients and diffusive fluxes are obtained in terms of the flux–flux time correlations. Furthermore, the relationship between the derived transport coefficients and those of the Stefan–Maxwell and infinitely dilute frameworks are presented, and the connection between the transport matrix and experimentally measurable quantities is described. To exemplify the application of the derived Green–Kubo relations in molecular simulations, the matrix of transport coefficients for lithium and chloride ions in dimethyl sulfoxide is computed using classical molecular dynamics and compared with experimental measurements.

42 ENGINEERING↗

DOE OSTI2026

Vibrational Spin-Orbit Coupling Contributions to Excited State Decay of Ligand-to-Ligand Charge Transfer States

Controlling excited state relaxation processes is important in a variety of photochemical and photophysical processes, including the generation of ground and excited state spin polarization for quantum information science applications. Here, we address how specific static distortions – based on vibrational spin-orbit active modes at C2v symmetry determined by group theory – in a series of low-symmetry ligand-to-ligand charge transfer complexes enable direct spin-orbit coupling contributions to T1 → S0 excited state decay. These results are used to address spin-vibronic coupling contributions to T1 → S0 decay in a high-symmetry (tBu2bpy)Pt(S,S) (tBu2bpy = 4,4’-di-tert-butyl-2,2’-bipyridine and S,S = benzene-1,2-dithiolate) ligand-to-ligand charge transfer complex with effective C2v symmetry, where T1 relaxation is both spin- and orbitally forbidden due to the direct spin-orbit coupling matrix element being zero by symmetry. Low-frequency vibrations that involve a pyridine-pyridine twisting motion within the bpy ligand generate large ∂/∂Qi values that will contribute significantly to T1 → S0 relaxation. The work advances ligand design strategies for the generation of tailored T1 → S0 relaxation rates, which can be utilized to optimize the generation of electron spin polarization in radical-elaborated ligand-to-ligand charge transfer complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1976

Two-electron bond-orbital model, 2

The two-electron bond-orbital model of tetrahedrally-coordinated solids is generalized and its application extended. All intrabond matrix elements entering the formalism are explicitly retained, including the direct overlap S between the anion and cation sp3 hybrid wavefunctions. Complete analytic results are obtained for the six two-electron eigenvalues and eigenstates of the anion-cation bond in terms of S, one-electron parameters V2 and V3, and two-electron correlation parameters V4, V5 and V6. Refined formulas for the dielectric constant and the nuclear exchange and pseudodipolar coefficients, as well as new expressions for the valence electron density, polarity of the bond and the cohesive energy, are then derived. The theory gives a good account of the experimentally observed trends in all properties considered and approximate quantitative agreement is achieved for the pseudodipolar coefficient.

Huang, C.↗

NASA NTRS2021

Theoretical Study of CH3Cl-N2 Line Shapes in the ν1 Band. Line Mixing Effects in QR Doublets and QQk Sub-Branches

Line coupling and line mixing effects of CH3Cl lines in the ν1 band perturbed by N2 have been investigated. We have taken into account the non-diagonality of the exp(−S2) operator within specified line spaces as well as the k-degeneracy of the transitions (due to the double degeneracy of the j,k levels with k ≠ 0). These transitions should be considered as doublets coupled by the line mixing process. A new problem appears in the calculation when the atom-atom potential model is introduced. In order to overcome this difficulty, a pragmatic approach is proposed. Comparisons of theoretically calculated matrix elements of W with measurements of QR(j,k) doublets as well as some QQk sub-branches, which are strongly affected by line mixing, have been made. The results show that the formalism improved in this way leads to rather accurate predictions.

Methyl chloride↗