Engineering Papers⌕ Search

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.

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

At least 1,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A matrix completion algorithm for efficient calculation of quantum and variational effects in chemical reactions

This work examines the viability of matrix completion methods as cost-effective alternatives to full nuclear Hessians for calculating quantum and variational effects in chemical reactions. The harmonic variety-based matrix completion (HVMC) algorithm, developed in a previous study (https://doi.org/10.1063/5.0018326), exploits the low-rank character of the polynomial expansion of potential energy to recover, using a small sample, vibrational frequencies (square roots of nuclear Hessian eigenvalues) constituting the reaction path. Furthermore, these frequencies are essential for calculating rate coefficients using variational transition state theory with multidimensional tunneling (VTST-MT). HVMC performance is examined for four SN2 reactions and five hydrogen transfer reactions, with each H-transfer reaction consisting of at least one vibrational mode strongly coupled to the reaction coordinate. HVMC is robust and captures zero-point energies, vibrational free energies, zero-curvature tunneling, and adiabatic ground state and free energy barriers as well as their positions on the reaction coordinate. For medium to large reactions involving H-transfer, with the exception of the most complex Ir catalysis system, less than 35% of total eigenvalue information is necessary for accurate recovery of key VTST-MT observables.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1978

Nonlinear flap-lag-axial equations of a rotating beam with arbitrary precone angle

In an attempt both to unify and extend the analytical basis of several aspects of the dynamic behavior of flexible rotating beams, the second-degree nonlinear equations of motion for the coupled flapwise bending, lagwise bending, and axial extension of an untwisted, torsionally rigid, nonuniform, rotating beam having an arbitrary angle of precone with the plane perpendicular to the axis of rotation are derived using Hamilton's principle. The derivation of the equations is based on the geometric nonlinear theory of elasticity and the resulting equations are consistent with the assumption that the strains are negligible compared to unity. No restrictions are imposed on the relative displacements or angular rotations of the cross sections of the beam other than those implied by the assumption of small strains. Illustrative numerical results, obtained by using an integrating matrix as the basis for the method of solution, are presented both for the purpose of validating the present method of solution and indicating the range of applicability of the equations of motion and the method of solution.

Kvaternik, R. G.↗
NASA NTRS1977

Derivation of the chemical-equilibrium rate coefficient using scattering theory

Scattering theory is applied to derive the equilibrium rate coefficient for a general homogeneous chemical reaction involving ideal gases. The reaction rate is expressed in terms of the product of a number of normalized momentum distribution functions, the product of the number of molecules with a given internal energy state, and the spin-averaged T-matrix elements. An expression for momentum distribution at equilibrium for an arbitrary molecule is presented, and the number of molecules with a given internal-energy state is represented by an expression which includes the partition function.

Mickens, R. E.↗
NASA NTRS1981

Nature and origin of basin-forming projectiles

The formation of the observed lunar multi-ring basins is discussed in the context of current theories of terrestrial planet formation, particularly those in which these planets formed by the accumulation of large planetesimals. The observed number, size, and timing of lunar basin-forming impacts is in the range expected for such theories. Tidal disruption during close encounters to earth and Venus can provide a single mechanism that explains a number of details concerning the number, size distribution, and stochastic nature of the timing of these impacts. A basin time scale is suggested in which Nectaris is associated with the 4.1 b.y. age of the Apollo 16 light matrix breccias. In accordance with the present consensus, Serenitatis is 3.86 b.y., Imbrium and Orientale 3.80-3.82 b.y. in age. Other nearside circular basins (e.g., Humorum and Crisium) are intermediate in age between 3.86 and 4.12 b.y. The large number of 3.8-3.9 b.y. ages is attributed primarily to the magnitude of the Imbrium and Serenitatis impacts, and sampling bias resulting from concentration of collection sites in proximity to these basins.

Wetherill, G. W.↗
NASA NTRS1983

Buckling of imperfect periodic lattice structures

A simplified buckling analysis is presented for a family of periodic lattice structures such as those proposed for large space structures. A transcendental 6 x 6 matrix of eigenvalues is shown to be sufficient for modeling buckling behavior because member stiffnesses are based on an exact solution of the beam-column equation. Exact stiffnesses are derived for a curved member, thus allowing modeling of imperfect lattice structures. Comparisons of predictions of the lattice model with those available from shell and beam theory underscore the inaccuracies introduced by treating the lattice structure as a continuum. Sample calculations are provided for an isogrid cylinder and a three element double-laced truss.

Anderson, M. S.↗
NASA NTRS1987

The capacity of the Hopfield associative memory

Techniques from coding theory are applied to study rigorously the capacity of the Hopfield associative memory. Such a memory stores n-tuple of + or - 1s. The components change depending on a hard-limited version of linear functions of all other components. With symmetric connections between components, a stable state is ultimately reached. By building up the connection matrix as a sum-of-outer products of m fundamental memories, it may be possible to recover a certain one of the m memories by using an initial n-tuple probe vector less than a Hamming distance n/2 away from the fundamental memory. If m fundamental memories are chosen at random, the maximum asymptotic value of m in order that most of the m original memories are exactly recoverable is n/(2 log n). With the added restriction that every one of the m fundamental memories be recoverable exactly, m can be no more than n/(4 log n) asymptotically as n approaches infinity. Extensions are also considered, in particular to capacity under quantization of the outer-product connection matrix. This quantized memory-capacity problem is closely related to the capacity of the quantized Gaussian channel.

Mceliece, Robert J.↗
NASA NTRS1995

Theoretical Infrared Spectra for Polycyclic Aromatic Hydrocarbon Neutrals, Cations and Anions

Calculations are carried out using density functional theory (DFT) to determine the harmonic frequencies and intensities of the neutrals and cations of thirteen polycyclic aromatic hydrocarbons (PAHs) up to the size of ovalene. Calculations are also carried out for a few PAH anions. The DFT harmonic frequencies, when uniformly scaled by the factor of 0.958 to account primarily for anharmonicity, agree with the matrix isolation fundamentals to within an average error of about 10 per centimeter. Electron correlation is found to significantly reduce the intensities of many of the cation harmonics, bringing them into much better agreement with the available experimental data. While the theoretical infrared spectra agree well with the experimental data for the neutral systems and for many of the cations, there are notable discrepancies with the experimental matrix isolation data for some PAH cations that are difficult to explain in terms of limitations in the calculations. In agreement with previous theoretical work, the present calculations show that the relative intensities for the astronomical unidentified infrared (UIR) bands agree reasonably well with those for a distribution of polycyclic aromatic hydrocarbon (PAH) cations, but not with a distribution of PAH neutrals. We also observe that the infrared spectra of highly symmetrical cations such as coronene agree much better with astronomical observations than do those of, for example, the polyacenes such as tetracene and pentacene. The total integrated intensities for the neutral species are found to increase linearly with size, while the total integrated intensities are much larger for the cations and scale more nearly quadratically with size. We conclude that emission from moderate-sized highly symmetric PAH cations such as coronene and larger could account for the UIR bands.

Langhoff, Stephen R.↗
NASA NTRS2017

Micromechanics of Composite Materials Governed by Vector Constitutive Laws

The high-fidelity generalized method of cells micromechanics theory has been extended for the prediction of the effective property tensor and the corresponding local field distributions for composites whose constituents are governed by vector constitutive laws. As shown, the shear analogy, which can predict effective transverse properties, is not valid in the general three-dimensional case. Consequently, a general derivation is presented that is applicable to both continuously and discontinuously reinforced composites with arbitrary vector constitutive laws and periodic microstructures. Results are given for thermal and electric problems, effective properties and local field distributions, ordered and random microstructures, as well as complex geometries including woven composites. Comparisons of the theory's predictions are made to test data, numerical analysis, and classical expressions from the literature. Further, classical methods cannot provide the local field distributions in the composite, and it is demonstrated that, as the percolation threshold is approached, their predictions are increasingly unreliable. XXXX It has been observed that the bonding between the fibers and matrix in composite materials can be imperfect. In the context of thermal conductivity, such imperfect interfaces have been investigated in micromechanical models by Dunn and Taya (1993), Duan and Karihaloo (2007), Nan et al. (1997) and Hashin (2001). The present HFGMC micromechanical method, derived for perfectly bonded composite materials governed by vector constitutive laws, can be easily generalized to include the effects of weak bonding between the constituents. Such generalizations, in the context of the mechanical micromechanics problem, involve introduction of a traction-separation law at the fiber/matrix interface and have been presented by Aboudi (1987), Bednarcyk and Arnold (2002), Bednarcyk et al. (2004) and Aboudi et al. (2013) and will be addressed in the future.

micromechanics↗
NASA NTRSDate not supplied

Multi-Domain Routing in Delay Tolerant Networks

The goal of Delay Tolerant Networking (DTN) is to provide the missing ingredient for the ever-growing collection of communicating nodes in our solar system to become a Solar System Internet (SSI). Great strides have been made in modeling particular types of DTNs, such as schedule- or discovery-based. Now, analogously to the Internet, these smaller DTNs can be considered routing domains which must be stitched together to form the overall SSI. In this paper, we propose a framework for cross-domain routing in DTNs as well as methodologies for detecting these sub-domains. Example time-varying networks are given to demonstrate the techniques proposed. A basic component is the mathematical theory of sheaves, which unifies the underlying model of DTN routing algorithms, by giving rise to routing sheaves – these can be defined for the dynamic and scheduled networks as noted above, and can also be used to define the interfaces between these domains in order to route across them. An immediate application would be routing across discovery-based networks connected by scheduled networks. These DTN subdomains remain elusive, however, and need to become well-defined and properly sized for tractable computability. In particular, a balance must be determined between areas that are too large (i.e. large matrix computations) versus areas that are too small (i.e. “many” single-noded domains). Moreover, the connections between the domains should, at least locally, be chosen to optimize data flow and connectivity: we address this in three ways. First, tools from persistent homology are given to understand underlying structures, reminiscent of hierarchies in the Internet Protocol (IP) addressing. Second, we construct a notion of temporal graph curvature based on network geometry to analyze flows induced by dynamical processes on these networks. Finally, Schrodinger Bridges, a tool arising from statistical physics, are proposed as a method of constructing flows on time-evolving networks with desirable properties such as speed, robustness, and load sensitivity. We construct an approach to temporal hypergraphs to simultaneously model unicast, multicast, and broadcast, using the language of scheme theory, and then consider DTN network coding as a way to achieve network-level computation and organization. The paper concludes with a discussion and ideas for future work.

Alan Hylton↗
NASA NTRS1995

Processing and property evaluation of metal matrix superconducting materials

Metal - superconductor (YBCO) systems have been prepared and characterized by resistivity, ac susceptibility and dc SQUID magnetic moment measurements. The silver composites showed superconducting transition for all the composites processed and the superconducting transition temperature tends to depend upon the concentration of the silver in the composite. Aluminum composites showed an unusual resistivity results with two transitions around 90 K and 120 K. The superconducting property of silver composites can be explained qualitatively in terms of the proximity theory that has been suggested for the low temperature superconductors.

Rao, Appajosula S.↗
DOE OSTI2020

Molecular Structure and Electron Affinity of Metal-Solvent Complexes: Insights from Density Functional Theory Simulations

A molecular level understanding of the structure and energetics of the monovalent and divalent metal ion complexes is of great importance for development of next-generation batteries. Here, Density Functional Theory (DFT) simulations at the ωb97xD/6-31 + G(d,p) level of theory are performed to investigate the interaction of metal ions (Li + , Na + , K + , Mg 2+ , Ca 2+ , Zn 2+ ) with 26 organic solvent molecules. The reduction energetics (electron affinity and reduction potential) and structural responses of the solvent molecules and the molecular complexes are discussed. The DFT calculations are carried out to investigate the structure, energetics and electron affinities of chelated complexes of water (H 2 O), tetrahydrofuran (THF) and di-methoxy ethane (DME) solvent molecules. Additionally, ab initio dynamic simulations (AIMD) at 298 K using atom centered density matrix propagation (ADMP) formalism are performed to understand the spontaneous structure formation upon electron attachment of the metal ion-solvent complexes. The ADMP simulations indicate the decomposition of Mg + -(DME) 3 complex via cleavage of C–O bond of one of the three DME molecules indicating irreversible decomposition of DME in the presence of the Mg + radical. We believe that the data collected as part of this investigation serves as a library of fundamental knowledge towards a deeper understanding of the electrode-electrolyte interfacial reactions.

25 ENERGY STORAGE↗
DOE OSTI2021

DMFTwDFT: An open-source code combining Dynamical Mean Field Theory with various density functional theory packages

We report Dynamical Mean Field Theory (DMFT) is a successful method to compute the electronic structure of strongly correlated materials, especially when it is combined with density functional theory (DFT). Here, we present an open-source computational package (and a library) combining DMFT with various DFT codes interfaced through the Wannier90 package. The correlated subspace is expanded as a linear combination of Wannier functions introduced in the DMFT approach as local orbitals. In particular, we provide a library mode for computing the DMFT density matrix. This library can be linked and then internally called from any DFT package, assuming that a set of localized orbitals can be generated in the correlated subspace. The existence of this library allows developers of other DFT codes to interface with our package and achieve the charge-self-consistency within DFT+DMFT loops. To test and check our implementation, we computed the density of states and the band structure of well-known solid-state correlated materials, namely LaNiO, SrVO, and NiO. The obtained results are compared to those obtained from other DFT+DMFT implementations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A highly stable Cu(OH) 2 -Poly(vinyl alcohol) nanocomposite membrane for dramatically enhanced direct borohydride fuel cell performance

Nano-additive aggregation, limited performance, and unclear modification mechanisms are the main obstacles in developing nanocomposite anion exchange membranes (AEMs). In this work, for the first time, an effective and highly stable Poly(vinyl alcohol) (PVA)-based AEMs with dispersive Cu(OH)$_2$ nanoclusters (Cu-AEMs) are prepared by a simple and eco-friendly three-step method: ‘CuCl$_2$ doping-casting-KOH immersing’. The doped Cu$^{2+}$ ions chemically combine with OH- ions to form anionic conductive Cu(OH)2 nanoclusters intermediated by attaching resins. The PVA skeletons wrap around Cu(OH)$_2$ while the hydroxyl groups expose to bulk water, forming Cu(OH)$_2$-PVA complex, which avoids nano-additive aggregation, increases anionic channels, and strengthens additive-matrix connection. A direct borohydride fuel cell using Cu-AEM with 0.56 wt % CuCl$_2$ possess the highest power density of 403.3 mW cm$^{-2}$ at 60 °C and a life span of over 200 h. The high-performance and durability come from the unique structure of Cu(OH)$_2$ nanoclusters-PVA complex. ‘Vehicle’ theory is considered to be the dominant mechanism for enhancing such nanocomposite Cu-AEMs. This work demonstrates a new concept for preparing the stable AEMs toward high-performance fuel cells. The synthetic chemistry involved can be broadly extended for fabricating versatile AEMs.

25 ENERGY STORAGE↗
DOE OSTI2023

Necessary and sufficient conditions for quadratic stabilizability of switched systems on non-uniform time domains

In this paper, we consider the quadratic stabilizability via state feedback for a particular class of switched systems that evolve on a non-uniform time domain by introducing time scales theory. The system considered switches between a continuous-time subsystem with variable lengths and a discrete-time subsystem with variable discrete step sizes. Necessary and sufficient conditions are derived to guarantee the quadratic stability of this class of switched systems via a switching state feedback law based on the existence of a common positive definite matrix satisfying the quadratic stabilizability condition by considering that the two subsystems are unstable. By state feedback, we mean that the switching among subsystems depends on the system states. Current results for this kind of state switching feedback control are derived only for switched systems evolving on a continuous time domain or a discrete time domain with fixed step’s size. These results are not applicable for the particular class of switched systems where there is a mixing between the continuous and discrete dynamics. This motivates the derivation of a new and more general state feedback control law for switched systems in this work. Here, a numerical example illustrating the results is presented.

42 ENGINEERING↗
NASA NTRS1974

Fluid-structure finite-element vibrational analysis

A fluid finite element has been developed for a quasi-compressible fluid. Both kinetic and potential energy are expressed as functions of nodal displacements. Thus, the formulation is similar to that used for structural elements, with the only differences being that the fluid can possess gravitational potential, and the constitutive equations for fluid contain no shear coefficients. Using this approach, structural and fluid elements can be used interchangeably in existing efficient sparse-matrix structural computer programs such as SPAR. The theoretical development of the element formulations and the relationships of the local and global coordinates are shown. Solutions of fluid slosh, liquid compressibility, and coupled fluid-shell oscillation problems which were completed using a temporary digital computer program are shown. The frequency correlation of the solutions with classical theory is excellent.

Feng, G. C.↗
NASA NTRS1988

An efficient Mindlin finite strip plate element based on assumed strain distribution

A simple two node, linear, finite strip plate bending element based on Mindlin-Reissner plate theory for the analysis of very thin to thick bridges, plates, and axisymmetric shells is presented. The new transverse shear strains are assumed for constant distribution in the two node linear strip. The important aspect is the choice of the points that relate the nodal displacements and rotations through the locking transverse shear strains. The element stiffness matrix is explicitly formulated for efficient computation and ease in computer implementation. Numerical results showing the efficiency and predictive capability of the element for analyzing plates with different supports, loading conditions, and a wide range of thicknesses are given. The results show no sign of the shear locking phenomenon.

Chulya, Abhisak↗
NASA NTRS1992

A parameter estimation algorithm for spatial sine testing - Theory and evaluation

This paper presents the theory and an evaluation of a spatial sine testing parameter estimation algorithm that uses directly the measured forced mode of vibration and the measured force vector. The parameter estimation algorithm uses an ARMA model and a recursive QR algorithm is applied for data reduction. In this first evaluation, the algorithm has been applied to a frequency response matrix (which is a particular set of forced mode of vibration) using a sliding frequency window. The objective of the sliding frequency window is to execute the analysis simultaneously with the data acquisition. Since the pole values and the modal density are obtained from this analysis during the acquisition, the analysis information can be used to help determine the forcing vectors during the experimental data acquisition.

Rost, R. W.↗
DOE OSTI2021

Fast and robust all-electron density functional theory calculations in solids using orthogonalized enriched finite elements

Here, we present a computationally efficient approach to perform systematically convergent real-space all-electron Kohn-Sham density functional theory calculations for solids using an enriched finite element (FE) basis. The enriched FE basis is constructed by augmenting the classical FE basis with atom-centered numerical basis functions, comprising of atomic solutions to the Kohn-Sham problem. Notably, to improve the conditioning, we orthogonalize the enrichment functions with respect to the classical FE basis, without sacrificing the locality of the resultant basis. In addition to improved conditioning, this orthogonalization procedure also renders the overlap matrix block diagonal, greatly simplifying its inversion. Subsequently, we use a Chebyshev polynomial based filtering technique to efficiently compute the occupied eigenspace in each self-consistent field iteration. We demonstrate the accuracy and efficiency of the proposed approach on periodic unit cells and supercells. The benchmark studies show a staggering 130× speedup of the orthogonalized enriched FE basis over the classical FE basis. We also present a comparison of the orthogonalized enriched FE basis with the linearized augmented plane-wave + local orbitals basis, both in terms of accuracy and efficiency. Notably, we demonstrate that the orthogonalized enriched FE basis can handle large system sizes of ~10 000 electrons. Finally, we observe good parallel scalability of our implementation with 92% efficiency at 22× speedup for a system with 620 electrons

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗