Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix theories”

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 541 records · Page 30

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.↗

Improving inference with matrix elements and machine learning

Particle physics processes bring together a high-energy amplitude described by quantum field theory and nonperturbative effects and detector interactions described by complex computer simulations. We review some recently developed multivariate inference techniques that leverage this structure and combine matrix-element information with machine learning. Automated by the MadMiner package, the new techniques have been applied to multiple problems in particle physics, allowing for stronger limits than traditional analysis methods and showing their potential to improve the sensitivity of the LHC legacy measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

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↗

Field Theories for type-II fractons

We derive an effective field theory for a type-II fracton starting from the Haah code on the lattice. The effective topological theory is not given exclusively in terms of an action; it must be supplemented with a condition that selects physical states. Without the constraint, the action only describes a type-I fracton. The constraint emerges from a condition that cube operators multiply to the identity, and it cannot be consistently implemented in the continuum theory at the operator level, but only in a weaker form, in terms of matrix elements of physical states. Informed by these studies and starting from the opposite end, i.e., the continuum, we discuss a Chern-Simons-like theory that does not need a constraint or projector, and yet has no mobile excitations. Whether this continuum theory admits a lattice counterpart remains unanswered.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

Eigenvalue Placement and Stabilization by Constrained Optimization

A pole placement algorithm is proposed which uses constrained nonlinear optimization techniques on a finite dimensional model of a linear n degree of freedom system. Low order feedback control is assumed where r poles may be assigned; r being the rank of the sensor coefficient matrix. It is shown that by combining feedback control theory methods with optimization techniques, one can ensure the stability characteristics of a system, and can alter its transient response.

Decaro, S. M.↗

A slewing control experiment for flexible structures

A hardware set-up has been developed to study slewing control for flexible structures including a steel beam and a solar panel. The linear optimal terminal control law is used to design active controllers which are implemented in an analog computer. The objective of this experiment is to demonstrate and verify the dynamics and optimal terminal control laws as applied to flexible structures for large angle maneuver. Actuation is provided by an electric motor while sensing is given by strain gages and angle potentiometer. Experimental measurements are compared with analytical predictions in terms of modal parameters of the system stability matrix and sufficient agreement is achieved to validate the theory.

Juang, J.-N.↗

Composite blade structural analyzer (COBSTRAN) user's manual

The installation and use of a computer code, COBSTRAN (COmposite Blade STRuctrual ANalyzer), developed for the design and analysis of composite turbofan and turboprop blades and also for composite wind turbine blades was described. This code combines composite mechanics and laminate theory with an internal data base of fiber and matrix properties. Inputs to the code are constituent fiber and matrix material properties, factors reflecting the fabrication process, composite geometry and blade geometry. COBSTRAN performs the micromechanics, macromechanics and laminate analyses of these fiber composites. COBSTRAN generates a NASTRAN model with equivalent anisotropic homogeneous material properties. Stress output from NASTRAN is used to calculate individual ply stresses, strains, interply stresses, thru-the-thickness stresses and failure margins. Curved panel structures may be modeled providing the curvature of a cross-section is defined by a single value function. COBSTRAN is written in FORTRAN 77.

Aiello, Robert A.↗

Composite Blade Structural Analyzer (COBSTRAN) theoretical/programmer's manual

This manual describes the organization and flow of data and analysis in the computer code, COBSTRAN (COmposite Blade STRuctural ANalyzer). This code combines composite mechanics and laminate theory with an internal data base of fiber and matrix properties and was developed for the design and analysis of composite turbofan and turboprop blades and composite wind turbine blades. Inputs to the code are constituent fiber and matrix material properties, factors reflecting the fabrication process, composite geometry and blade geometry. COBSTRAN performs the micromechanics and laminate analyses of these fiber composites and generates a NASTRAN finite element model of the blade. This manual describes the equations formulated and solved in the code and the function of each of the seventy-two subroutines. COBSTRAN is written in FORTRAN 77.

Aiello, Robert A.↗

Numerical simulation of spatially-evolving instability

A computational study of the spatial stability of plane Poiseuille flow is presented. The numerical scheme employs a time-splitting method to integrate the full Navier-Stokes equations using spectral collocation/finite-difference discretization on a non-staggered mesh. The eigenvalue decomposition procedure is applied for the solution of the Poisson equations using the capacitance matrix technique. The buffer domain method is incorporated for the outflow boundary conditions. The input perturbation velocities are obtained by solving the Orr-Sommerfeld equation for the nonlinear eigenvalue problem employing the companion matrix method. Computational results are compared with the linear theory for two-dimensional disturbances.

Danabasoglu, G.↗

On the thermodynamic framework of generalized coupled thermoelastic-viscoplastic-damage modeling

A complete potential based framework using internal state variables is put forth for the derivation of reversible and irreversible constitutive equations. In this framework, the existence of the total (integrated) form of either the (Helmholtz) free energy or the (Gibbs) complementary free energy are assumed a priori. Two options for describing the flow and evolutionary equations are described, wherein option one (the fully coupled form) is shown to be over restrictive while the second option (the decoupled form) provides significant flexibility. As a consequence of the decoupled form, a new operator, i.e., the Compliance operator, is defined which provides a link between the assumed Gibb's and complementary dissipation potential and ensures a number of desirable numerical features, for example the symmetry of the resulting consistent tangent stiffness matrix. An important conclusion reached, is that although many theories in the literature do not conform to the general potential framework outlined, it is still possible in some cases, by slight modifications of the used forms, to restore the complete potential structure.

Arnold, S. M.↗

Unsymmetric Lanczos model reduction and linear state function observer for flexible structures

This report summarizes part of the research work accomplished during the second year of a two-year grant. The research, entitled 'Application of Lanczos Vectors to Control Design of Flexible Structures' concerns various ways to use Lanczos vectors and Krylov vectors to obtain reduced-order mathematical models for use in the dynamic response analyses and in control design studies. This report presents a one-sided, unsymmetric block Lanczos algorithm for model reduction of structural dynamics systems with unsymmetric damping matrix, and a control design procedure based on the theory of linear state function observers to design low-order controllers for flexible structures.

Su, Tzu-Jeng↗

Probabilistic sizing of laminates with uncertainties

A reliability based design methodology for laminate sizing and configuration for a special case of composite structures is described. The methodology combines probabilistic composite mechanics with probabilistic structural analysis. The uncertainties of constituent materials (fiber and matrix) to predict macroscopic behavior are simulated using probabilistic theory. Uncertainties in the degradation of composite material properties are included in this design methodology. A multi-factor interaction equation is used to evaluate load and environment dependent degradation of the composite material properties at the micromechanics level. The methodology is integrated into a computer code IPACS (Integrated Probabilistic Assessment of Composite Structures). Versatility of this design approach is demonstrated by performing a multi-level probabilistic analysis to size the laminates for design structural reliability of random type structures. The results show that laminate configurations can be selected to improve the structural reliability from three failures in 1000, to no failures in one million. Results also show that the laminates with the highest reliability are the least sensitive to the loading conditions.

Shah, A. R.↗