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,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Hierarchical microstructures: a potential route to enhanced stability in structural materials for advanced nuclear reactors

The drive to increase efficiency in nuclear energy systems is leading to the need for materials that operate at higher temperatures and stress levels for extended periods, while maintaining stable microstructures to ensure their performance is not compromised. A novel route to producing materials that can perform well in such environments is the creation of hierarchical microstructures. A hierarchical microstructure is a microstructure in which features are present at multiple length scales simultaneously. In this work, a hierarchical microstructure is fabricated in a nickel-base superalloy, featuring nanometer-size gamma precipitates inside larger gamma-prime particles, which are in turn embedded in the gamma matrix phase. The hierarchical features of the microstructure lead to enhanced stability of the gamma-prime precipitates during annealing; the particle size does not follow the expected t^(1/3) growth law predicted by the classic LSW theory. Phase-field simulations are used to understand the unexpected stability of the gamma-prime precipitates.

36 - MATERIALS SCIENCE↗
DOE OSTI2021

Discrete Green’s functions and spectral graph theory for computationally efficient thermal modeling

Here, this work concerns solutions of the heat equation with the spectral graph method, for which the temperature is defined at discrete points in the domain and the spatial relationship among the points is described by a graph. The heat equation on the graph is solved using matrix techniques involving the eigenvectors and eigenvalues of the Laplacian matrix. The spectral graph approach precludes the computationally intensive meshing and numerous time-integration steps of the finite element method. In the present work, the spectral graph method is extended to include heat loss at the boundaries with a generalized boundary condition, and physics-based edge weights are introduced which simplify the calibration process. From this approach a discrete Green’s function is defined which allows for solutions under a variety of heating conditions including: space-varying initial conditions; time-and-space varying internal heating; and, time-and-space-varying heating at boundaries of type 1 (Dirichlet), type 2 (Neumann) and type 3 (Robin). Results are provided for benchmark heat transfer problems in one spatial dimension and in three spatial dimensions, and verification is provided by comparison with exact analytical solutions and finite difference solutions. The spectral graph method converges within 0.4% error of the analytical solution. The practical utility of the approach is demonstrated by thermal simulation of a multilayer additive manufacturing process. The spectral graph results are compared to experimentally-obtained temperature data for two metal parts, with error less than 5% of the experimental measurements, with computation time less than one minute on a desktop computer.

36 MATERIALS SCIENCE↗
NASA NTRS1974

On the theory of large amplitude Alfven waves

Large amplitude Alfvenic disturbances of arbitrary spatial shape and polarization are described by MHD equations, without resort to the usual assumption of planarity. However, because of their nonplanar nature, the direction of propagation of these disturbances cannot, in general, be determined by looking for minima in a variance matrix constructed from observed field fluctuations. When such minima exist, one is observing that subset of interplanetary Alfven waves that is essentially planar.

Goldstein, M. L.↗
NASA NTRS1976

A counter-example in linear feature selection theory

The paper shows that it is possible to construct two k x n matrices, both of which maximize divergence in the transformed space of the linear feature selection problem in multiclass pattern recognition, and which are not row equivalent. Thus, even under extremely strong conditions, it is not possible to assume that all matrix solutions which maximize transformed divergence are row equivalent.

Brown, D. R.↗
DOE OSTI2020

Microstructural development in DED stainless steels: applying welding models to elucidate the impact of processing and alloy composition

Austenitic stainless steel microstructures produced by directed energy deposition (DED)are analogous to those developed during welding, particularly high energy density welding. To better understand microstructural development during DED, theories of microstructural evolution,which have been established to contextualize weld microstructures, are applied in this study to microstructural development in DED austenitic stainless steels. Phenomenological welding models that describe the development of oxide inclusions, compositional microsegregation, ferrite,matrix austenite grains, and dislocation substructures are utilized to clarify microstructural evolution during deposition of austenitic stainless steels. Two different alloys, 304L and 316L, arecompared to demonstrate the broad applicability of this framework for understanding microstmctural development during the DED process. Despite differences in grain morphology and solidification mode for these two alloys (which can be attributed to compositional differences),similar tensile properties are achieved. It is the fine-scale compositional segregation and dislocation structures that ultimately determine the strength of these materials. The evolution of microsegregation and dislocation structures is shown to be dependent on the rapid solidification and thermomechanical history of the DED processing method and not the composition of the starting material.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Computational Aspects of Single-Molecule Kinetics for Coupled Catalytic Cycles: A Spectral Analysis

Catalysis from single active sites is analyzed using methods developed from single molecule kinetics. Using a stochastic Markov state description, the observable properties of general catalytic networks of reactions are expressed using an eigenvalue decomposition of the transition matrix for the Markov process. By the use of a sensitivity analysis, the necessary eigenvalues and eigenvectors are related to the energies of controlling barriers and wells located along the reaction routes. A generalization of the energetic span theory allows the eigenvalues to be computed from several activation energies corresponding to distinct barrier-well pairings. The formalism is demonstrated accurately for model problems for a physically realistic mechanism for an alkene hydrogenation reaction on a single atom catalyst.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Upscaling Reactive Transport and Clogging in Shale Microcracks by Deep Learning

Fracture networks in shales exhibit multiscale features. A rock system may contain a few main fractures and thousands of microcracks, whose length and aperture are orders of magnitude smaller than the former. It is computationally prohibitive to resolve all the fractures explicitly for such multiscale fracture networks. One traditional approach is to model the small-scale features (e.g., microcracks in shales) as an effective medium. Although this fracture-matrix conceptualization significantly reduces the problem complexity, there are classes of physical processes that cannot be accurately upscaled by effective medium approximations, e.g., microcrack clogging during mineral reactions. Here, we employ deep learning in place of effective medium theory to upscale physical processes in small-scale features. Specifically, we consider reactive transport in a fracture-microcrack network where microcracks can be clogged by precipitation. A deep learning multiscale algorithm is developed, in which the microcracks are upscaled as a wall boundary condition of the main fractures. The wall boundary condition is constructed by recurrent neural networks, which take concentration histories as input and predict the solute transport from main fractures to microcracks. The deep learning multiscale algorithm is firstly employed in specific scenarios, then a general model is developed which can work under various conditions. The new approach is validated against fully resolved simulations and an analytical solution, providing a reliable and efficient solution for problems that cannot be upscaled by effective medium models.

58 GEOSCIENCES↗
DOE OSTI2025

Spectral Clustering-Based Partitioning of Large-Scale Power Electronics-Based Power Systems for Small-Signal Stability Analysis

The nodal admittance matrix (NAM)-based approach is well-suited for small-signal stability analysis of large-scale power electronics-based power systems (PEPSs), as it preserves the system structure through its admittance matrix. Previous studies have explored partitioning such systems into subareas and interconnections to reduce computational burden; however, they lacked a formal algorithmic procedure for determining feasible partitions. While several grid partitioning methods, such as those based on graph theory or machine learning, exist in the literature, they cannot be directly applied to NAM-based analysis due to differing objectives and constraints. Here, this paper addresses this gap by presenting a systematic, step-by-step procedure for applying a spectral partitioning algorithm that yields a division of the system into subareas suitable for NAM-based analysis. The computational complexity of the proposed method is also derived to demonstrate its efficiency and justify the practicality of the resulting subarea decomposition. The performance of the partitioning method is evaluated by applying the spectral clustering-derived subareas and interconnections to the NAM-based partitioning approach on a 140-bus system. Computational times for the full-system and partitioned NAM analyses are compared using MATLAB. Additionally, PSCAD simulations of the complete system and partitioned subareas are carried out to verify the effectiveness of the proposed method.

Nupur [Univ. of Tennessee, Knoxville, TN (United S↗
DOE OSTI2024

All about the neutron from lattice QCD

I describe how simulations of lattice QCD using the path integral formulation provide the two basic quantum mechanical properties of QCD, its ground state in which correlation functions are calculated, and Fock state wavefunctions between which matrix elements of operators are calculated. Both constructs are stochastic, so unfortunately one gets no intuitive picture or even qualitative understanding of what they look like, nevertheless they contain and display all the subtleties of the quantum field theory. Today, these simulations provide many quantities that are impacting phenomenology and experiments. I illustrate the methods and the steps in the analysis using, as examples, three observables: the isovector charges of the nucleon, the contribution of the quark’s intrinsic spin to the nucleon spin, and the pion-nucleon sigma term.

Gupta, Rajan (ORCID:0000000317843058)↗
NASA NTRS1962

GENERAL INSTABILITY OF ORTHOGONALLY STIFFENED CYLINDRICAL SHELLS

Earlier research at the National Aeronautics Research Institute (N.L.R.), Amsterdam, which forms the basis of recent work is reviewed. This early work refers to 2 schemes: the orthotropic shell and, in view of buckling modes where the half wave length is of the order of the ring distance, the shell with continuously distributed stringers and discrete rings. Linear theory is considered to be adequate for these structures, where the imperfections are small in comparison to the height of the ring sections. Recent developments account for pressure difference in addition to axial compression, for the correct stiffness matrix of skin panels in the post-buckling stage and for stringer bending due to hoop stresses in the skin, which are of importance as has been shown by the investigation of the post-buckling behaviour. Numerical data for the stiffness matrix of skin panels have been established. Numerical evaluation of the stability equation has not been performed as yet.

SHELL STABILITY↗
NASA NTRS1972

Quantum-mechanical transport equation for atomic systems.

A quantum-mechanical transport equation (QMTE) is derived which should be applicable to a wide range of problems involving the interaction of radiation with atoms or molecules which are also subject to collisions with perturber atoms. The equation follows the time evolution of the macroscopic atomic density matrix elements of atoms located at classical position R and moving with classical velocity v. It is quantum mechanical in the sense that all collision kernels or rates which appear have been obtained from a quantum-mechanical theory and, as such, properly take into account the energy-level variations and velocity changes of the active (emitting or absorbing) atom produced in collisions with perturber atoms. The present formulation is better suited to problems involving high-intensity external fields, such as those encountered in laser physics.

Berman, P. R.↗
NASA NTRS1987

Optimal control concepts in design sensitivity analysis

A close link is established between open loop optimal control theory and optimal design by noting certain similarities in the gradient calculations. The resulting benefits include a unified approach, together with physical insights in design sensitivity analysis, and an efficient approach for simultaneous optimal control and design. Both matrix displacement and matrix force methods are considered, and results are presented for dynamic systems, structures, and elasticity problems.

Belegundu, Ashok D.↗
NASA NTRS1987

Crosslinking-property relationships in PMR polyimide composites. I

The thermooxidatively-induced crosslinking/ physical and mechanical property relationships of graphite fiber-reinforced PMR polyimide-matrix composites were studied during isothermal exposure of the composite specimens at 288 C in air for periods of up to 5000 hr. The crosslinking densities due to this treatment were estimated on the basis of the kinetic theory of rubber elasticity and shifts in the glass transition temperature T(g). Several linear relationships are noted between crosslink density and physical and mechanical properties: T(g), initial weight loss, and elevated temperature interlaminar shear strength increase with crosslink density, while initial moisture absorption decreases. After achieving the highest crosslink density, several of the composite properties begin to decrease from their maximum values.

Pater, R. H.↗
NASA NTRS1989

Efficient computer algebra algorithms for polynomial matrices in control design

The theory of polynomial matrices plays a key role in the design and analysis of multi-input multi-output control and communications systems using frequency domain methods. Examples include coprime factorizations of transfer functions, cannonical realizations from matrix fraction descriptions, and the transfer function design of feedback compensators. Typically, such problems abstract in a natural way to the need to solve systems of Diophantine equations or systems of linear equations over polynomials. These and other problems involving polynomial matrices can in turn be reduced to polynomial matrix triangularization procedures, a result which is not surprising given the importance of matrix triangularization techniques in numerical linear algebra. Matrices with entries from a field and Gaussian elimination play a fundamental role in understanding the triangularization process. In the case of polynomial matrices, matrices with entries from a ring for which Gaussian elimination is not defined and triangularization is accomplished by what is quite properly called Euclidean elimination. Unfortunately, the numerical stability and sensitivity issues which accompany floating point approaches to Euclidean elimination are not very well understood. New algorithms are presented which circumvent entirely such numerical issues through the use of exact, symbolic methods in computer algebra. The use of such error-free algorithms guarantees that the results are accurate to within the precision of the model data--the best that can be hoped for. Care must be taken in the design of such algorithms due to the phenomenon of intermediate expressions swell.

Baras, J. S.↗
NASA NTRS2006

HFGMC Enhancement of MAC/GMC

Additional information about a mathematical model denoted the high-fidelity generalized method of cells (HFGMC) and implementation of the HFGMC within version 4.0 of the MAC/GMC software has become available. MAC/GMC (Micromechanics Analysis Code With Generalized Method of Cells) was a topic of several prior NASA Tech Briefs articles, version 4.0 having been described in "Comprehensive Micromechanics-Analysis Code - Version 4.0" (LEW-17495-1), NASA Tech Briefs, Vol. 29, No. 9 (September 2005), page 54. MAC/GMC predicts elastic and inelastic thermomechanical responses of composite materials. MAC/GMC utilizes the generalized method of cells (GMC) - a model of micromechanics that predicts macroscopic responses of a composite material as functions of the properties, sizes, shapes, and responses of its constituents (e.g., matrix and fibers). The accuracy of the GMC is limited by neglect of coupling between normal and shear stresses. The HFGMC was developed by combining elements of the GMC and a related model, denoted the higher-order theory for functionally graded materials (HOTFGM), that can account for this coupling. Hence, the HFGMC enables simulation of stress and strain with greater accuracy. Some alterations of the MAC/GMC data structure were necessitated by the greater computational complexity of the HFGMC.

Arnold, Steven M.↗
NASA NTRSDate not supplied

Progressive Damage Response of 3D Woven Composites via the Multiscale Recursive Micromechanics Solution with Tailored Fidelity

Progressive failure simulations have been performed for orthogonal 3D woven composites consisting of RTM6 resin matrix and AS4 carbon fibers. The Multiscale Recursive Micromechanics approach has been used, which, while being computationally efficient, captures the primary effects of the microstructure at each considered length scale. This approach also enables use of any micromechanics theory at any length scale, and herein, the fidelity of the chosen theories across the scales has been tailored to strike a balance with computational efficiency. The Mori-Tanaka method is employed at the lowest length scale, the Generalized Method of Cells is used at intermediate scales, and the High-Fidelity Generalized Method of Cells is used at the highest woven composite repeating unit cell scale. Furthermore, two different damage models, also with different levels of fidelity and efficiency, have been used for the resin material at the lowest length scale. Results for the mechanical behavior in response to loading in various directions are compared for the two damage models and with available test data.

3D woven↗
NASA NTRS2001

Higher-Order Theory for Functionally Graded Materials

Functionally graded materials (FGM's) are a new generation of engineered materials wherein the microstructural details are spatially varied through nonuniform distribution of the reinforcement phase(s). Engineers accomplish this by using reinforcements with different properties, sizes, and shapes, as well as by interchanging the roles of the reinforcement and matrix phases in a continuous manner (ref. 1). The result is a microstructure that produces continuously or discretely changing thermal and mechanical properties at the macroscopic or continuum scale. This new concept of engineering the material's microstructure marks the beginning of a revolution both in the materials science and mechanics of materials areas since it allows one, for the first time, to fully integrate the material and structural considerations into the final design of structural components. Functionally graded materials are ideal candidates for applications involving severe thermal gradients, ranging from thermal structures in advanced aircraft and aerospace engines to computer circuit boards. Owing to the many variables that control the design of functionally graded microstructures, full exploitation of the FGM's potential requires the development of appropriate modeling strategies for their response to combined thermomechanical loads. Previously, most computational strategies for the response of FGM's did not explicitly couple the material's heterogeneous microstructure with the structural global analysis. Rather, local effective or macroscopic properties at a given point within the FGM were first obtained through homogenization based on a chosen micromechanics scheme and then subsequently used in a global thermomechanical analysis.

Aboudi, J.↗
DOE OSTI2021

Hydrodynamics of nonintegrable systems from a relaxation-time approximation

We develop a general kinetic theory framework to describe the hydrodynamics of strongly interacting, nonequilibrium quantum systems in which integrability is weakly broken, leaving a few residual conserved quantities. This framework is based on a generalized relaxation-time approximation; it gives a simple, but surprisingly accurate, prescription for computing nonequilibrium transport even in strongly interacting systems. We validate the predictions of this approximation against matrix product operator calculations on chaotic quantum spin chains, finding surprisingly good agreement. Lastly, we show that despite its simplicity, our framework can capture phenomena distinctive to strongly interacting systems, such as widely separated charge and energy diffusion constants.

1-dimensional spin chains↗