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 991 records · Page 55

All sources NASA NTRS DOE OSTI

NASA NTRS1996

A Note on Substructuring Preconditioning for Nonconforming Finite Element Approximations of Second Order Elliptic Problems

In this paper an algebraic substructuring preconditioner is considered for nonconforming finite element approximations of second order elliptic problems in 3D domains with a piecewise constant diffusion coefficient. Using a substructuring idea and a block Gauss elimination, part of the unknowns is eliminated and the Schur complement obtained is preconditioned by a spectrally equivalent very sparse matrix. In the case of quasiuniform tetrahedral mesh an appropriate algebraic multigrid solver can be used to solve the problem with this matrix. Explicit estimates of condition numbers and implementation algorithms are established for the constructed preconditioner. It is shown that the condition number of the preconditioned matrix does not depend on either the mesh step size or the jump of the coefficient. Finally, numerical experiments are presented to illustrate the theory being developed.

Maliassov, Serguei↗

NASA NTRS2015

3D Progressive Damage Modeling for Laminated Composite Based on Crack Band Theory and Continuum Damage Mechanics

A simple continuum damage mechanics (CDM) based 3D progressive damage analysis (PDA) tool for laminated composites was developed and implemented as a user defined material subroutine to link with a commercially available explicit finite element code. This PDA tool uses linear lamina properties from standard tests, predicts damage initiation with an easy-to-implement Hashin-Rotem failure criteria, and in the damage evolution phase, evaluates the degradation of material properties based on the crack band theory and traction-separation cohesive laws. It follows Matzenmiller et al.'s formulation to incorporate the degrading material properties into the damaged stiffness matrix. Since nonlinear shear and matrix stress-strain relations are not implemented, correction factors are used for slowing the reduction of the damaged shear stiffness terms to reflect the effect of these nonlinearities on the laminate strength predictions. This CDM based PDA tool is implemented as a user defined material (VUMAT) to link with the Abaqus/Explicit code. Strength predictions obtained, using this VUMAT, are correlated with test data for a set of notched specimens under tension and compression loads.

Wang, John T.↗

NASA NTRS1981

Creep and creep rupture of laminated graphite/epoxy composites

An incremental numerical procedure based on lamination theory is developed to predict creep and creep rupture of general laminates. Existing unidirectional creep compliance and delayed failure data is used to develop analytical models for lamina response. The compliance model is based on a procedure proposed by Findley which incorporates the power law for creep into a nonlinear constitutive relationship. The matrix octahedral shear stress is assumed to control the stress interaction effect. A modified superposition principle is used to account for the varying stress level effect on the creep strain. The lamina failure model is based on a modification of the Tsai-Hill theory which includes the time dependent creep rupture strength. A linear cumulative damage law is used to monitor the remaining lifetime in each ply.

Dillard, D. A.↗

NASA NTRS2015

Microstructural Influence on Deformation and Fatigue Life of Composites Using the Generalized Method of Cells

A fully coupled deformation and damage approach to modeling the response of composite materials and composite laminates is presented. It is based on the semi-‐analytical generalized method of cells (GMC) micromechanics model as well as its higher fidelity counterpart, HFGMC, both of which provide closed-form constitutive equations for composite materials as well as the micro scale stress and strain fields in the composite phases. The provided constitutive equations allow GMC and HFGMC to function within a higher scale structural analysis (e.g., finite element analysis or lamination theory) to represent a composite material point, while the availability of the micro fields allow the incorporation of lower scale sub‐models to represent local phenomena in the fiber and matrix. Further, GMC's formulation performs averaging when applying certain governing equations such that some degree of microscale field accuracy is surrendered in favor of extreme computational efficiency, rendering the method quite attractive as the centerpiece in a integrated computational material engineering (ICME) structural analysis; whereas HFGMC retains this microscale field accuracy, but at the price of significantly slower computational speed. Herein, the sensitivity of deformation and the fatigue life of graphite/epoxy PMC composites, with both ordered and disordered microstructures, has been investigated using this coupled deformation and damage micromechanics based approach. The local effects of fiber breakage and fatigue damage are included as sub‐models that operate on the microscale for the individual composite phases. For analysis of laminates, classical lamination theory is employed as the global or structural scale model, while GMC/HFGMC is embedded to operate on the microscale to simulate the behavior of the composite material within each laminate layer. A key outcome of this study is the statistical influence of microstructure and micromechanics idealization (GMC or HFGMC) on the overall accuracy of unidirectional and laminated composite deformation and fatigue response.

Microscale Models↗

DOE OSTI2024

Spectral scheme for atomic structure calculations in density functional theory

In this study, we present a spectral scheme for atomic structure calculations in pseudopotential Kohn-Sham density functional theory. In particular, after applying an exponential transformation of the radial coordinates, we employ global polynomial interpolation on a Chebyshev grid, with derivative operators approximated using the Chebyshev differentiation matrix, and integrations using Clenshaw-Curtis quadrature. We demonstrate the accuracy and efficiency of the scheme through spin-polarized and unpolarized calculations for representative atoms, while considering local, semilocal, and hybrid exchange-correlation functionals. In particular, we find that $\mathcal{O}$(200) grid points are sufficient to achieve an accuracy of 1 microhartree in the eigenvalues for optimized norm conserving Vanderbilt pseudopotentials spanning the periodic table from atomic number Ζ = 1 to 83.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1997

ICAN Computer Code Adapted for Building Materials

The NASA Lewis Research Center has been involved in developing composite micromechanics and macromechanics theories over the last three decades. These activities have resulted in several composite mechanics theories and structural analysis codes whose applications range from material behavior design and analysis to structural component response. One of these computer codes, the Integrated Composite Analyzer (ICAN), is designed primarily to address issues related to designing polymer matrix composites and predicting their properties - including hygral, thermal, and mechanical load effects. Recently, under a cost-sharing cooperative agreement with a Fortune 500 corporation, Master Builders Inc., ICAN was adapted to analyze building materials. The high costs and technical difficulties involved with the fabrication of continuous-fiber-reinforced composites sometimes limit their use. Particulate-reinforced composites can be thought of as a viable alternative. They are as easily processed to near-net shape as monolithic materials, yet have the improved stiffness, strength, and fracture toughness that is characteristic of continuous-fiber-reinforced composites. For example, particlereinforced metal-matrix composites show great potential for a variety of automotive applications, such as disk brake rotors, connecting rods, cylinder liners, and other hightemperature applications. Building materials, such as concrete, can be thought of as one of the oldest materials in this category of multiphase, particle-reinforced materials. The adaptation of ICAN to analyze particle-reinforced composite materials involved the development of new micromechanics-based theories. A derivative of the ICAN code, ICAN/PART, was developed and delivered to Master Builders Inc. as a part of the cooperative activity.

Murthy, Pappu L. N.↗

DOE OSTI2021

Magnetic states of the quasi-one-dimensional iron chalcogenide Ba 2 FeS 3

Quasi-one-dimensional iron-based ladders and chains, with the 3d iron electronic density n=6, are attracting considerable attention. Recently, a new iron chain system Ba 2 FeS 3 , also with n=6, was prepared under high-pressure and high-temperature conditions. Here the magnetic and electronic phase diagrams are theoretically studied for this quasi-one-dimensional compound. Based on first-principles calculations, a strongly anisotropic one-dimensional electronic band behavior near the Fermi level was observed. In addition, a three-orbital electronic Hubbard model for this chain was constructed. Introducing the Hubbard and Hund couplings and studying the model via the density matrix renormalization group (DMRG) method, we studied the ground-state phase diagram. A robust staggered ↑-↓-↑-↓ AFM region was unveiled in the chain direction, consistent with our density functional theory (DFT) calculations. Furthermore, at intermediate Hubbard U coupling strengths, this system was found to display an orbital selective Mott phase (OSMP) with one localized orbital and two itinerant metallic orbitals. At very large U/W (W=bandwidth), the system displays Mott insulator characteristics, with two orbitals half-filled and one doubly occupied. Our results for high pressure Ba 2 FeS 3 provide guidance to experimentalists and theorists working on this one-dimensional iron chalcogenide chain material.

1-dimensional systems↗

NASA NTRS1987

A simplified computer solution for the flexibility matrix of contacting teeth for spiral bevel gears

A computer code, FLEXM, was developed to calculate the flexibility matrices of contacting teeth for spiral bevel gears using a simplified analysis based on the elementary beam theory for the deformation of gear and shaft. The simplified theory requires a computer time at least one order of magnitude less than that needed for the complete finite element method analysis reported earlier by H. Chao, and it is much easier to apply for different gear and shaft geometries. Results were obtained for a set of spiral bevel gears. The teeth deflections due to torsion, bending moment, shearing strain and axial force were found to be in the order 10(-5), 10(-6), 10(-7), and 10(-8) respectively. Thus, the torsional deformation was the most predominant factor. In the analysis of dynamic load, response frequencies were found to be larger when the mass or moment of inertia was smaller or the stiffness was larger. The change in damping coefficient had little influence on the resonance frequency, but has a marked influence on the dynamic load at the resonant frequencies.

Hsu, C. Y.↗

DOE OSTI2021

Holographic duality for Ising CFT with boundary

We extend the holographic duality between 3d pure gravity and 2d Ising CFT proposed in ref. [1] to CFTs with boundaries. Besides the usual asymptotic boundary, the dual bulk spacetime now has a real cutoff, on which live branes with finite tension, giving Neumann boundary condition on the metric tensor. The strongly coupled bulk theory requires that we dress the well-known semiclassical AdS/BCFT answer with boundary gravitons, turning the partition function into the form of Virasoro characters. Using this duality, we relate the brane tensions to the modular S-matrix elements of the dual BCFT and derive the transformation between gravitational solutions with different brane tensions under modular S action.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

T$ \overline{T} $ deformation in SCFTs and integrable supersymmetric theories

We calculate the $ \mathcal{S} $ -multiplets for two-dimensional Euclidean $ \mathcal{N} $ = (0 , 2) and $ \mathcal{N} $ = (2 , 2) superconformal field theories under the T $ \overline{T} $ deformation at leading order of perturbation theory in the deformation coupling. Then, from these $ \mathcal{N} $ = (0 , 2) deformed multiplets, we calculate two- and three-point correlators. We show the $ \mathcal{N} $ = (0 , 2) chiral ring’s elements do not flow under the T $ \overline{T} $ deformation. Specializing to integrable supersymmetric seed theories, such as $ \mathcal{N} $ = (2 , 2) Landau-Ginzburg models, we use the thermodynamic Bethe ansatz to study the S-matrices and ground state energies. From both an S-matrix perspective and Melzer’s folding prescription, we show that the deformed ground state energy obeys the inviscid Burgers’ equation. Finally, we show that several indices independent of D -term perturbations including the Witten index, Cecotti-Fendley-Intriligator-Vafa index and elliptic genus do not flow under the T $ \overline{T} $ deformation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Magnon interactions in the quantum paramagnetic phase of CoNb 2 O 6

In this work, we study effects of magnon interactions in the excitation spectrum of CoNb 2 O 6 in the quantum paramagnetic phase in transverse field, where the 1/S spin-wave theory exhibits unphysical divergences at the critical field. We propose a self-consistent Hartree-Fock approach that eliminates such unphysical singularities while preserving the integrity of the singular threshold phenomena of magnon decay and spectrum renormalization that are present in both theory and experiment. With the microscopic parameters adopted from previous studies, this method yields a close quantitative agreement with the available experimental data for CoNb 2 O 6 in the relevant regime. Furthermore, insights into the general structure of the spin-anisotropic model of CoNb 2 O 6 and related zigzag chain materials are also provided and a discussion of the effects of additional longitudinal field on the spectrum is given.

1-dimensional spin chains↗

DOE OSTI2021

Symbol alphabets from plabic graphs III: n = 9

Symbol alphabets of n-particle amplitudes in N = 4 super-Yang-Mills theory are known to contain certain cluster variables of G(4, n) as well as certain algebraic functions of cluster variables. In this paper we solve the C Z = 0 matrix equations associated to several cells of the totally non-negative Grassmannian, combining methods of arXiv:2012.15812 for rational letters and arXiv:2007.00646 for algebraic letters. We identify sets of parameterizations of the top cell of G + (5, 9) for which the solutions produce all of (and only) the cluster variable letters of the 2-loop nine-particle NMHV amplitude, and identify plabic graphs from which all of its algebraic letters originate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Spin-phonon dispersion in magnetic materials

Microscopic coupling between the electron spin and the lattice vibration is responsible for an array of exotic properties from morphic effects in simple non magnets to magnetodielectric coupling in multiferroic spinels and hematites. Traditionally, a single spin–phonon coupling constant is used to characterize how effectively the lattice can affect the spin, but it is hardly enough to capture novel electromagnetic behaviors to the full extent. Here, we introduce a concept of spin–phonon dispersion to project the spin moment change along the phonon crystal momentum direction, so the entire spin change can be mapped out. Different from the phonon dispersion, the spin–phonon dispersion has both positive and negative frequency branches even in the equilibrium ground state, which correspond to the spin enhancement and spin reduction, respectively. Our study of bcc Fe and hcp Co reveals that the spin force matrix, that is, the second-order spatial derivative of spin moment, is similar to the vibrational force matrix, but its diagonal elements are smaller than the off-diagonal ones. This leads to the distinctive spin–phonon dispersion. Furthermore, the concept of spin–phonon dispersion expands the traditional Elliott-Yafet theory in nonmagnetic materials to the entire Brillouin zone in magnetic materials, thus opening the door to excited states in systems such as CoF 2 and NiO, where a strong spin-lattice coupling is detected in the THz regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2000

Nonlinearity and Strain-Rate Dependence in the Deformation Response of Polymer Matrix Composites Modeled

There has been no accurate procedure for modeling the high-speed impact of composite materials, but such an analytical capability will be required in designing reliable lightweight engine-containment systems. The majority of the models in use assume a linear elastic material response that does not vary with strain rate. However, for containment systems, polymer matrix composites incorporating ductile polymers are likely to be used. For such a material, the deformation response is likely to be nonlinear and to vary with strain rate. An analytical model has been developed at the NASA Glenn Research Center at Lewis Field that incorporates both of these features. A set of constitutive equations that was originally developed to analyze the viscoplastic deformation of metals (Ramaswamy-Stouffer equations) was modified to simulate the nonlinear, rate-dependent deformation of polymers. Specifically, the effects of hydrostatic stresses on the inelastic response, which can be significant in polymers, were accounted for by a modification of the definition of the effective stress. The constitutive equations were then incorporated into a composite micromechanics model based on the mechanics of materials theory. This theory predicts the deformation response of a composite material from the properties and behavior of the individual constituents. In this manner, the nonlinear, rate-dependent deformation response of a polymer matrix composite can be predicted.

Goldberg, Robert K.↗

NASA NTRS1974

Acoustic theory of axisymmetric multisectioned ducts

Equations are developed for the acoustic field in a duct system which is made up of a number of connected circular and annular ducts. These equations are suitable for finding the acoustic field inside of and radiated from an aircraft turbofan engine. Acoustic modes are used as generalized coordinates in order to develop a set of matrix equations for the acoustic field. Equations for these modes are given for circular and annular ducts with uniform flow. Modal source equations are derived for point acoustic sources. General equations for mode transmission and reflection are developed and detailed equations are derived for ducts with multiple sections of acoustic treatment and for ducts with circumferential splitter rings. The general theory is applied to the special case of a uniform area circular duct with multisection liners and it is shown that the mode reflection effects are proportional to differences of the acoustic admittances of adjacent liners. A numerical example is given which shows that multisection liners may provide greater noise suppression than uniform liners.

Zorumski, W. E.↗

DOE OSTI2024

Flattening of the EFT-hedron: supersymmetric positivity bounds and the search for string theory

We examine universal positivity constraints on 2 → 2 scattering in 4d planar N = 4 supersymmetric Yang-Mills theory with higher-derivative corrections. We present numerical evidence that the convex region of allowed Wilson coefficients (the “EFT-hedron”) flattens completely along about one-third of its dimensions when an increasing number of constraints on the spectral density from crossing-symmetry are included. Our analysis relies on the formulation of the positivity constraints as a linear optimization problem, which we implement using two numerical solvers, SDPB and CPLEX. Motivated by the flattening, we propose a novel partially resummed low-energy expansion of the 2 → 2 amplitude. As part of the analysis, we provide additional evidence in favor of the conjecture [1] that the Veneziano amplitude is the only amplitude compatible with both S-matrix bootstrap constraints and string monodromy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

S-matrix bootstrap in 3+1 dimensions: regularization and dual convex problem

The S-matrix bootstrap maps out the space of S-matrices allowed by analyticity, crossing, unitarity, and other constraints. For the 2 → 2 scattering matrix S 2→2 such space is an infinite dimensional convex space whose boundary can be determined by maximizing linear functionals. On the boundary interesting theories can be found, many times at vertices of the space. Here we consider 3 + 1 dimensional theories and focus on the equivalent dual convex minimization problem that provides strict upper bounds for the regularized primal problem and has interesting practical and physical advantages over the primal problem. Its variables are dual partial waves k ℓ (s) that are free variables, namely they do not have to obey any crossing, unitarity or other constraints. Nevertheless they are directly related to the partial waves f ℓ (s), for which all crossing, unitarity and symmetry properties result from the minimization. Numerically, it requires only a few dual partial waves, much as one wants to possibly match experimental results. We consider the case of scalar fields which is related to pion physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Determining the Ensemble N -Representability of Reduced Density Matrices

The N-representability problem for reduced density matrices remains a fundamental challenge in electronic structure theory. Following our previous work that employs a unitary-evolution algorithm based on an adaptive derivative-assembled pseudo-Trotter variational quantum algorithm to probe pure-state N-representability of reduced density matrices [J. Chem. Theory Comput. 2024, 20, 9968], in this work we propose a practical framework for determining the ensemble N-representability of a p-body matrix. This is accomplished using a purification strategy that embeds an ensemble state into a pure state defined on an extended Hilbert space, such that the reduced density matrices of the purified state reproduce those of the original ensemble. By iteratively applying variational unitaries to an initial purified state, the proposed algorithm minimizes the Hilbert-Schmidt distance between its p-body reduced density matrix and a specified target p-body matrix, which serves as a measure of the N-representability of the target. This methodology facilitates both error correction of defective ensemble reduced density matrices and quantum-state reconstruction on a quantum computer, offering a route for density-matrix refinement. We validate the algorithm with numerical simulations on systems of two, three, and four electrons in both simple models as well as molecular systems at finite temperature, demonstrating its robustness.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗