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 325 records · Page 18

Telescoping Mechanics: A New Paradigm for Composite Behavior Simulation

This report reviews the application of telescoping mechanics to composites using recursive laminate theory. The elemental scale is the fiber-matrix slice, the behavior of which propagates to laminate. The results from using applications for typical, hybrid, and smart composites and composite-enhanced reinforced concrete structures illustrate the versatility and generality of telescoping scale mechanics. Comparisons with approximate, single-cell, and two- and three-dimensional finite-element methods demonstrate the accuracy and computational effectiveness of telescoping scale mechanics for predicting complex composite behavior.

Chamis, C. C.↗

Micromechanical modeling of damage growth in titanium based metal-matrix composites

The thermomechanical behavior of continuous-fiber reinforced titanium based metal-matrix composites (MMC) is studied using the finite element method. A thermoviscoplastic unified state variable constitutive theory is employed to capture inelastic and strain-rate sensitive behavior in the Timetal-21s matrix. The SCS-6 fibers are modeled as thermoplastic. The effects of residual stresses generated during the consolidation process on the tensile response of the composites are investigated. Unidirectional and cross-ply geometries are considered. Differences between the tensile responses in composites with perfectly bonded and completely debonded fiber/matrix interfaces are discussed. Model simulations for the completely debonded-interface condition are shown to correlate well with experimental results.

Sherwood, James A.↗

Research in Theoretical Particle Physics: Neutrino Masses and The Origin of CP-Violation

The main goal of this project was to investigate the theories of neutrino masses at the low scale and the origin of CP-violation in physics beyond the Standard Model. The origin of neutrino masses is one of the most pressing issues in particle physics. In this proposal, we investigated the theories for neutrino masses based on local total lepton number. In these theories anomaly cancellation predicts the existence of extra fermions with lepton number, and one of them is a good dark matter (DM) candidate. The cosmological constraints on the DM relic density implies that the lepton number symmetry breaking scale must be below the multi-TeV scale. Therefore, one can hope to test the origin of neutrino masses in current or future experiments. We investigated in great detail the anomaly cancellation in these gauge theories, study the different mechanisms for neutrino masses, study the predictions for direct and indirect dark matter experiments. Since these theories predict new sources of CP violation, we investigated the predictions for the electric dipole moments (EDM). The Higgs decays and the different signatures at the Large Hadron Collider were investigated in great detail. Finally, we studied the possible baryogenesis mechanisms in these theories in agreement with the EDM, DM and collider constraints. The origin of CP violation in the Standard Model is unknown. In this proposal, we investigated the simplest mechanisms for spontaneous CP-violation to explain the CP-violation in the CKM matrix and the value of the QCD vacuum angle. We discussed the Nelson-Barr mechanism in gauge theories predicting vector-like quarks from anomaly cancellation such as theories for local baryon number. We discussed the Bento-Branco-Parada mechanism in gauge theories to explain the CP-Violation in the CKM matrix. We will investigate models with two Higgs doublets in the context of the minimal theory for quark-lepton unification. We will show the non-decoupling effects in the Higgs sector, the main constraints coming from flavour violating processes when the Yukawa couplings are related by the gauge symmetry in theories for quark-lepton unification. We will investigate the predictions for electric dipole moments in these theories. The possibility to have successful baryogenesis was investigated. Finally, we investigated the relation between CP-violation in the quark and leptonic sectors. The future results from these studies could help us to understand two main issues in physics beyond the Standard Model: The Origin of Neutrino Masses and CP-violation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Weighted graph based ordering techniques for preconditioned conjugate gradient methods

We describe the basis of a matrix ordering heuristic for improving the incomplete factorization used in preconditioned conjugate gradient techniques applied to anisotropic PDE's. Several new matrix ordering techniques, derived from well-known algorithms in combinatorial graph theory, which attempt to implement this heuristic, are described. These ordering techniques are tested against a number of matrices arising from linear anisotropic PDE's, and compared with other matrix ordering techniques. A variation of RCM is shown to generally improve the quality of incomplete factorization preconditioners.

Clift, Simon S.↗

Discretized partial differential equations - Examples of control systems defined on modules

The purpose of this paper is to show how the important problems of linear system theory can be solved concisely for a particular class of linear systems, namely block circulant systems, by exploiting the algebraic structure. This type of system arises in lumped approximations to linear partial differential equations. The computation of the transition matrix, the variation of constants formula, observability, controllability, pole allocation, realization theory, stability and quadratic optimal control are discussed. In principle, all questions which are solved here could also be solved by standard methods; the present paper clearly exposes the structure of the solution, and thus permits various savings in computational effort.

Brockett, R. W.↗

A fiber matrix model for fluid flow and streaming potentials in the canaliculi of an osteon

A theoretical model is developed to predict the fluid shear stress and streaming potential at the surface of osteocytic processes in the lacunar-canalicular porosity of an osteon when the osteon is subject to mechanical loads that are parallel or perpendicular to its axis. The theory developed in Weinbaum et al. (31) for the flow through a proteoglycan matrix in a canaliculus is employed in a poroelastic model for the osteon. Our formulation is a generalization of that of Petrov et al. (17). Our model predicts that, in order to satisfy the measured frequency dependence of the phase and magnitude of the SGP in macroscopic bone samples, the fiber spacing in the fluid annulus must lie in the narrow range 6-7 nm typical of the spacing of GAG sidechains along a protein monomer. The model predictions for the local SGP profiles in the osteon agree with the experimental observations of Starkebaum et al. (24). The theory predicts that the pore pressure relaxation time, tau d, for a 150-300 microns diameter osteon with the foregoing matrix structure is approximately 0.03-0.13 sec, and that the amplitude of the mean fluid shear stress on the membrane of the osteocytic process at the mean areal radius of the osteon has a maximum at 28 Hz if tau d = 0.06 sec. This maximum, which is independent of the magnitude of the loading, could be important in vivo since the recent experiments of Turner et al. (28) and McLeod et al. (15) have a peak in the strain frequency spectrum between 20 and 30 Hz that also appears to be independent of the type (magnitude) of loading. Numerical predictions for the amplitude of the average fluid shear stress on the osteocytic membrane at the mean areal radius of the osteon show that the fluid shear stress associated with the low amplitude 20-30 Hz spectral strain component is at least as large as the average fluid shear stress associated with the high amplitude 1 Hz stride component, although the latter loading is an order of magnitude larger, and has a magnitude that lies within the middle of the range, 6-30 dynes/cm2, where fluid shear stresses in tissue culture studies with osteoblast monolayers have elicited an intracellular Ca++ response (31). The implications of these results for intracellular electrical communication are discussed.

Non-NASA Center↗

Laminate behavior for SiC fiber-reinforced reaction-bonded silicon nitride matrix composites

The room temperature mechanical properties of SiC fiber reinforced reaction-bonded silicon nitride matrix composite laminates (SiC/RBSN) have been measured. The laminates contained approx 30 volume fraction of aligned 142-micron diameter SiC fiber in a porous RBSN matrix. Three types of laminate studied were unidirectional: (1) (0) sub 8, (2) (10) sub 8, and (3) (45) sub 8, and (90) sub 8; cross plied laminates (0 sub 2/90 sub 2); and angle plied laminates: (+45 sub 2/-45 sub 2). Each laminate contained eight fiber plies. Results of the unidirectionally reinforced composites tested at various angles to the reinforcement direction indicate large anisotropy in in-plane properties. In addition, strength properties of these composites along the fiber direction were independent of specimen gage length and were unaffected by notches normal to the fiber direction. Splitting parallel to the fiber at the notch tip appears to be the dominant crack blunting mechanism responsible for notch insensitive behavior of these composites. In-plane properties of the composites can be improved by 2-D laminate construction. Mechanical property results for (0 sub 2/90 sub 2)sub s and (+45/-45 sub 2) sub s laminates showed that their matrix failure strains were similar to that for (0) sub 8 laminates, but their primary elastic moduli, matrix cracking strengths, and ultimate composite strengths were lower. The elastic properties of unidirectional, cross-ply, and angle-ply composites can be predicted from modified constitutive equations and laminate theory. Further improvements in laminate properties may be achieved by reducing the matrix porosity and by optimizing the bond strength between the SiC fiber and RBSN matrix.

Rhatt, R. T.↗

Laminate behavior for SiC fiber-reinforced reaction-bonded silicon nitride matrix composites

The room temperature mechanical properties of SiC fiber reinforced reaction-bonded silicon nitride matrix composite laminates (SiC/RBSN) have been measured. The laminates contained approx 30 volume fraction of aligned 142-micron diameter SiC fiber in a porous RBSN matrix. Three types of laminate studied were unidirectional: (1) (0) sub 8, (2) (10) sub 8, and (3) (45) sub 8, and (90) sub 8; cross plied laminates (0 sub 2/90 sub 2); and angle plied laminates: (+45 sub 2/-45 sub 2). Each laminate contained eight fiber plies. Results of the unidirectionally reinforced composites tested at various angles to the reinforcement direction indicate large anisotropy in in-plane properties. In addition, strength properties of these composites along the fiber direction were independent of specimen gage length and were unaffected by notches normal to the fiber direction. Splitting parallel to the fiber at the notch tip appears to be the dominant crack blunting mechanism responsible for notch insensitive behavior of these composites. In-plane properties of the composites can be improved by 2-D laminate construction. Mechanical property results for (0 sub 2/90 sub 2) sub s and (+45/-45 sub 2) sub s laminates showed that their matrix failure strains were similar to that for (0) sub 8 laminates, but their primary elastic moduli, matrix cracking strengths, and ultimate composite strengths were lower. The elastic properties of unidirectional, cross-ply, and angle-ply composites can be predicted from modified constitutive equations and laminate theory. Further improvements in laminate properties may be achieved by reducing the matrix porosity and by optimizing the bond strength between the SiC fiber and RBSN matrix.

Bhatt, Ramakrishna T.↗

A computational procedure to analyze metal matrix laminates with nonlinear lamination residual strains

An approximate computational procedure is described for the analysis of angleplied laminates with residual nonlinear strains. The procedure consists of a combination of linear composite mechanics and incremental linear laminate theory. The procedure accounts for initial nonlinear strains, unloading, and in-situ matrix orthotropic nonlinear behavior. The results obtained in applying the procedure to boron/aluminum angleplied laminates show that this is a convenient means to accurately predict the initial tangent properties of angleplied laminates in which the matrix has been strained nonlinearly by the lamination residual stresses. The procedure predicted initial tangent properties results which were in good agreement with measured data obtained from boron/aluminum angleplied laminates.

Chamis, C. C.↗

A computational procedure to analyze metal matrix laminates with nonlinear lamination residual strains

An approximate computational procedure is described for the analysis of angleplied laminates with residual nonlinear strains. The procedure consists of a combination of linear composite mechanics and incremental linear laminate theory. The procedure accounts for initial nonlinear strains, unloading, and in-situ matrix orthrotropic nonlinear behavior. The results obtained in applying the procedure to boron/aluminum angleplied laminates show that this is a convenient means to accurately predict the initial tangent properties of angleplied laminates in which the matrix has been strained nonlinearly by the lamination residual stresses. The procedure predicted initial tangent properties results which were in good agreement with measured data obtained from boron/aluminum angleplied laminates.

Chamis, C. C.↗

A Study of Feature Extraction Using Divergence Analysis of Texture Features

An empirical study of texture analysis for feature extraction and classification of high spatial resolution remotely sensed imagery (10 meters) is presented in terms of specific land cover types. The principal method examined is the use of spatial gray tone dependence (SGTD). The SGTD method reduces the gray levels within a moving window into a two-dimensional spatial gray tone dependence matrix which can be interpreted as a probability matrix of gray tone pairs. Haralick et al (1973) used a number of information theory measures to extract texture features from these matrices, including angular second moment (inertia), correlation, entropy, homogeneity, and energy. The derivation of the SGTD matrix is a function of: (1) the number of gray tones in an image; (2) the angle along which the frequency of SGTD is calculated; (3) the size of the moving window; and (4) the distance between gray tone pairs. The first three parameters were varied and tested on a 10 meter resolution panchromatic image of Maryville, Tennessee using the five SGTD measures. A transformed divergence measure was used to determine the statistical separability between four land cover categories forest, new residential, old residential, and industrial for each variation in texture parameters.

Hallada, W. A.↗

Transient Creep of a Composite Lower Crust: Constitutive Theory - 1

A composite model is proposed to describe the time-dependent response of the Earth's lower crust. The motivation for such it model is twofold: First, new observations of widespread postseismic deformation indicate that the deep continental crust responds viscoelastically, having both long-and short-term decay times. Second, by any number of observationally based rationales, the lower crust is compositionally and structurally heterogeneous over many length scales. For heterogeneities that have much smaller characteristic lengths than the minimum deformation wavelength of interest, the aggregate rheology can be described by composite media theory. For wavelengths of the order of the thickness of the lower crust (approx. = 25-40 km) and larger, composite theory may be applied to heterogeneities that are smaller than about several hundred meters, or equivalent to the vertical extent of a thick lower crustal mylonitic shear zone. The composite media theory developed here is constructed using both Eshelhy-Mori-Tanaka theory for aligned generalized spheroidal inclusions and a generalized self-consistent method. The inclusions and matrix are considered to be Maxwellian viscoelastic: a rheology that is consistent with past homogeneous models of postseismic stress relaxation. The composite theory presented here introduces a transient response to a suddenly imposed stress field which does not appear in homogeneous Maxwell models. Analytic expressions for the amplitude and duration of the transient and for the effective long-and short-term viscosities of the composite are given which describe the sensitivity to inclusion concentration (phi), to shape, and to ratio of inclusion-to-matrix viscosity (R).

Ivins, Erik R.↗

Effective field theories for collective excitations of atomic nuclei

Collective modes emerge as the relevant degrees of freedom that govern low-energy excitations of atomic nuclei. These modes—rotations, pairing rotations, and vibrations—are separated in energy from non-collective excitations, making it possible to describe them in the framework of effective field theory. Rotations and pairing rotations are the remnants of Nambu–Goldstone modes from the emergent breaking of rotational symmetry and phase symmetries in finite deformed and finite superfluid nuclei, respectively. The symmetry breaking severely constrains the structure of low-energy Lagrangians and thereby clarifies what is essential and simplifies the description. The approach via effective field theories exposes the essence of nuclear collective excitations and is defined with a breakdown scale in mind. This permits one to make systematic improvements and to estimate and quantify uncertainties. Effective field theories of collective excitations have been used to compute spectra, transition rates, and other matrix elements of interest. In particular, predictions of the nuclear matrix element for neutrinoless double beta decay then come with quantified uncertainties. This review summarizes these results and also compares the approach via effective field theories to well-known models and ab initio computations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measurements and interpretations of W ± Z production cross-sections in pp collisions at $\sqrt{s}=13$ TeV with the ATLAS detector

Measurements of integrated and differential cross-sections for W ± Z production in proton-proton collisions are presented. The data collected by the ATLAS detector at the Large Hadron Collider from 2015 to 2018 at a centre-of-mass energy of $\sqrt{s}=13$ TeV are used, corresponding to an integrated luminosity of 140 fb −1 . The W ± Z candidate events are reconstructed using leptonic decay modes of the gauge bosons into electrons or muons. The integrated cross-section per lepton flavour for the production of W ± Z is measured in the detector fiducial region with a relative precision of 4%. The measured value is compared with the Standard Model prediction at a precision of up to next-to-next-to-leading-order in QCD and next-to-leading-order in electroweak. Cross-sections for W + Z and W − Z production and their ratio are presented. The W ± Z production is also measured differentially as functions of various kinematic variables, including new observables sensitive to CP-violation effects. All measurements are compared with state-of-the-art Standard Model predictions from fixed-order calculations or Monte Carlo generators based on next-to-leading-order matrix elements interfaced with parton showers. An effective field theory interpretation of the measurements is performed, considering both CP-conserving and CP-violating dimension-6 operators modifying the W ± Z production. In the absence of observed deviations from the Standard Model, limits on CP-conserving Wilson coefficients are extracted using the transverse mass of the W ± Z system. For CP-violating coefficients a machine learning approach is used to construct an observable with enhanced sensitivity to CP-violation effects.

hadron-hadron scattering↗

Random Phase Approximation Correlation Energy Using Real-Space Density Functional Perturbation Theory

We present a real-space method for computing the random phase approximation (RPA) correlation energy within Kohn–Sham density functional theory, leveraging the low-rank nature of the frequency-dependent density response operator. In particular, we employ a cubic-scaling formalism based on density functional perturbation theory that circumvents the calculation of the response function matrix, instead relying on the ability to compute its product with a vector through the solution of the associated Sternheimer linear systems. We develop a large-scale parallel implementation of this formalism using the subspace iteration method in conjunction with the spectral quadrature method while employing the Kronecker product-based method for the application of the Coulomb operator and the conjugate orthogonal conjugate gradient method for the solution of the linear systems. We demonstrate convergence with respect to key parameters and verify the method’s accuracy by comparing with plane-wave results. We show that the framework achieves good strong scaling to many thousands of processors, reducing the time to solution for a lithium hydride system with 128 electrons to around 150 s on 4608 processors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Applications of flow models to the generation of correlated lattice QCD ensembles

Machine-learned normalizing flows can be used in the context of lattice quantum field theory to generate statistically correlated ensembles of lattice gauge fields at different action parameters. This work demonstrates how these correlations can be exploited for variance reduction in the computation of observables. Three different proof-of-concept applications are demonstrated using a novel residual flow architecture: continuum limits of gauge theories, the mass dependence of QCD observables, and hadronic matrix elements based on the Feynman–Hellmann approach. In all three cases, it is shown that statistical uncertainties are significantly reduced when machine-learned flows are incorporated as compared with the same calculations performed with uncorrelated ensembles or direct reweighting. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Probabilistic Error Bounds for Low-Rank Tensor Decompositions Used in Large-Scale Data Analysis Applications (LDRD Final Report)

This report documents a research project on analyzing low-rank tensor models for data analysis that took place at Sandia National Laboratories from October 2023–September 2025. The focus of this work was to extend theoretical frameworks from statistics and probability theory for use with models for scalar, vector, and matrix data to models with tensor, or general multi-dimensional array, data. Through this work, we have provided a new set of tools for bounding errors on low-rank tensor models of both complete and sampled data. The remainder of this report is organized as follows. In Section 1, we describe the proposed work at the start of the project. Section 2 describes the research advances made as part of the project. Other research contributions in the form of conference presentations and software development is provided in Section 3. Workforce development at Sandia and Florida Atlantic University (via a subcontract on this project) is provided in Section 4.

97 MATHEMATICS AND COMPUTING↗