Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix models”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Ground-state phase diagram of the t-t ′ -J model

Significance Theoretical efforts to understand the superconducting cuprates have focused on a few simplified but challenging models: notably the Hubbard and the t - t ′ - J models. Using density matrix renormalization group (DMRG) simulations on width-8 cylinders, we study the ground-state phase diagram of the t - t ′- J model. The next nearest-neighbor hopping t ′ allows one to distinguish between hole ( t ′ < 0 ) and electron doping ( t ′ > 0 ). We find for the magnetic, charge, and single-particle properties that this model roughly captures the differences between the electron- and hole-doped cuprates. However, the superconducting properties differ. For hole doping ( t ′ < 0 ), there is an absence of superconductivity, while for electron doping ( t ′ > 0), we find d -wave superconductivity coexisting with antiferromagnetism and a secondary triplet p -wave pairing component.

Jiang, Shengtao↗
DOE OSTI2024

Accuracy optimized neural networks do not effectively model optic flow tuning in brain area MSTd

Accuracy-optimized convolutional neural networks (CNNs) have emerged as highly effective models at predicting neural responses in brain areas along the primate ventral stream, but it is largely unknown whether they effectively model neurons in the complementary primate dorsal stream. We explored how well CNNs model the optic flow tuning properties of neurons in dorsal area MSTd and we compared our results with the Non-Negative Matrix Factorization (NNMF) model, which successfully models many tuning properties of MSTd neurons. To better understand the role of computational properties in the NNMF model that give rise to optic flow tuning that resembles that of MSTd neurons, we created additional CNN model variants that implement key NNMF constraints – non-negative weights and sparse coding of optic flow. While the CNNs and NNMF models both accurately estimate the observer's self-motion from purely translational or rotational optic flow, NNMF and the CNNs with nonnegative weights yield substantially less accurate estimates than the other CNNs when tested on more complex optic flow that combines observer translation and rotation. Despite its poor accuracy, NNMF gives rise to tuning properties that align more closely with those observed in primate MSTd than any of the accuracy-optimized CNNs. This work offers a step toward a deeper understanding of the computational properties and constraints that describe the optic flow tuning of primate area MSTd.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Multi deep learning-based stochastic microstructure reconstruction and high-fidelity micromechanics simulation of time-dependent ceramic matrix composite response

A multi deep learning-based framework is developed for efficient, automated microstructure reconstruction and generation of stochastic representative volume elements (SRVEs) with periodic boundary conditions (PBCs) for accurate modeling of ceramic matrix composite (CMC) response. The methodology comprises a convolutional neural network coupled with regression layers to act as a vanilla regression network for semantic segmentation of the microstructure, allowing accurate characterization of the phases and their distributions at the microscale. Scanning electron microscope and confocal microscope are used to obtain C/SiNC and SiC/SiNC CMCs micrographs for vanilla regression testing. Microstructure variability in terms of fiber volume fraction and porosity are quantified through the output regression layer, ensuring accurate representation of material variability in SRVE construction. Generative adversarial network (GAN) and its variants are designed to produce high-fidelity SRVE, spanning CMCs microstructure variability space. A circular padding algorithm is developed to generate SRVEs with PBCs during training of GANs. The accuracy of the generated SRVEs is established through micromechanics simulations, where an efficient formulation of the high-fidelity generalized methods of cells (HFGMC) approach is used to compute the effective mechanical properties. Furthermore, an iterative algorithm is implemented in the HFGMC solver to simulate time-dependent deformation of SiC/SiNC subjected to creep loading conditions.

36 MATERIALS SCIENCE↗
NASA NTRS1972

Structural model optimization using statistical evaluation

The results of research in applying statistical methods to the problem of structural dynamic system identification are presented. The study is in three parts: a review of previous approaches by other researchers, a development of various linear estimators which might find application, and the design and development of a computer program which uses a Bayesian estimator. The method is tried on two models and is successful where the predicted stiffness matrix is a proper model, e.g., a bending beam is represented by a bending model. Difficulties are encountered when the model concept varies. There is also evidence that nonlinearity must be handled properly to speed the convergence.

Collins, J. D.↗
NASA NTRS1973

Decomposition-aggregation stability analysis

This report presents the development and description of the decomposition aggregation approach to stability investigations of high dimension mathematical models of dynamic systems. The high dimension vector differential equation describing a large dynamic system is decomposed into a number of lower dimension vector differential equations which represent interconnected subsystems. Then a method is described by which the stability properties of each subsystem are aggregated into a single vector Liapunov function, representing the aggregate system model, consisting of subsystem Liapunov functions as components. A linear vector differential inequality is then formed in terms of the vector Liapunov function. The matrix of the model, which reflects the stability properties of the subsystems and the nature of their interconnections, is analyzed to conclude over-all system stability characteristics. The technique is applied in detail to investigate the stability characteristics of a dynamic model of a hypothetical spinning Skylab.

Siljak, D. D.↗
NASA NTRS1988

Nonlinear modeling, identification and estimation of slew induced structural deformation

This proposal discusses a new nonlinear, nonparametric method for off-line modeling and on-line estimation of the deformation of a flexible structure undergoing rapid retargeting maneuvers. In these circumstances, the structural stiffness and damping coefficients depend on the angular acceleration omega(dot), the angular rate omega, and the square of the angular rate omega. In the single axis case, the excitation of the structure is represented by the vector u(exp T) = (omega(dot), omega(exp 2), 2(omega)), to which the structural dynamics responds as a 'bilinear' (i.e., parametrically excited) system. A similar technique for multiaxial rotations yields a bilinear model with respect to matrix valued excitations. Three methods of estimation and modeling are described in this proposal to achieve deformation state determination: (1) a method based on a feedback linearized procedure which gives an estimate of the state by means of observers installed in the deformable body; (2) off-line modeling of the deformation state of the structure by means of topological interpolators; and (3) an on-line structural state estimation method based on a combination of the two previous techniques.

Dwyer, Thomas A., III↗
NASA NTRS2001

Ultrasonic Assessment of Impact-Induced Damage and Microcracking in Polymer Matrix Composites

The main objective of this NASA FAR project is to conduct ultrasonic assessment of impact-induced damage and microcracking in polymer matrix composites at various temperatures. It is believed that the proposed study of impact damage assessment on polymer matrix composites will benefit several NASA missions and current interests, such as ballistic impact testing of composite fan containment and high strain rate deformation modeling of polymer matrix composites. Impact-induced damage mechanisms in GLARE and ARALL fiber-metal laminates subject to instrumented drop-weight impacts at various temperatures were studied. GLARE and ARALL are hybrid composites made of alternating layers of aluminum and glass (for GLARE) and aramid- (for ARALL) fiber-reinforced epoxy. Damage in pure aluminum panels impacted by foreign objects was mainly characterized by large plastic deformation surrounding a deep penetration dent. On the other hand, plastic deformation in fiber-metal laminates was often not as severe although the penetration dent was still produced. The more stiff fiber-reinforced epoxy layers provided better bending rigidity; thus, enhancing impact damage tolerance. Severe cracking, however, occurred due to the use of these more brittle fiber-reinforced epoxy layers. Fracture patterns, e.g., crack length and delamination size, were greatly affected by the lay-up configuration rather than by the number of layers, which implies that thickness effect was not significant for the panels tested in this study. Immersion ultrasound techniques were then used to assess damages generated by instrumented drop-weight impacts onto these fiber-metal laminate panels as well as 6061-T6 aluminum/cast acrylic sandwich plates adhered by epoxy. Depending on several parameters, such as impact velocity, mass, temperature, laminate configuration, sandwich construction, etc., various types of impact damage were observed, including plastic deformation, radiating cracks emanating from the impact site, ring cracks surrounding the impact site, partial and full delamination, and combinations of these damages.

Liaw, Benjamin↗
DOE OSTI2021

Method and apparatus for shape-based energy analysis of solids

A computer-readable medium stores instructions including storing a model of a physical structure and defining a mesh for the model. The mesh includes finite elements, each defined by a respective set of edges. The instructions include, for each finite element, identifying a governing differential equation and a set of complementary functions that exactly satisfy the governing differential equation. The instructions include determining an applied physical stimulus for the physical structure. The instructions include generating an energy optimization model that minimizes a difference between internal energy and external energy of the finite elements in response to the applied physical stimulus. The model includes a matrix of respective scalar multipliers for the complementary functions for each finite element. The instructions include transforming the matrix and calculating a physical parameter of interest. The instructions include, in response to the physical parameter not satisfying a design parameter, updating the model of the physical structure.

Spears, Robert E.↗
NASA NTRS2006

An Initial Non-Equilibrium Porous-Media Model for CFD Simulation of Stirling Regenerators

The objective of this paper is to define empirical parameters for an initial thermal non-equilibrium porous-media model for use in Computational Fluid Dynamics (CFD) codes for simulation of Stirling regenerators. The two codes currently used at Glenn Research Center for Stirling modeling are Fluent and CFD-ACE. The codes porous-media models are equilibrium models, which assume solid matrix and fluid are in thermal equilibrium. This is believed to be a poor assumption for Stirling regenerators; Stirling 1-D regenerator models, used in Stirling design, use non-equilibrium regenerator models and suggest regenerator matrix and gas average temperatures can differ by several degrees at a given axial location and time during the cycle. Experimentally based information was used to define: hydrodynamic dispersion, permeability, inertial coefficient, fluid effective thermal conductivity, and fluid-solid heat transfer coefficient. Solid effective thermal conductivity was also estimated. Determination of model parameters was based on planned use in a CFD model of Infinia's Stirling Technology Demonstration Converter (TDC), which uses a random-fiber regenerator matrix. Emphasis is on use of available data to define empirical parameters needed in a thermal non-equilibrium porous media model for Stirling regenerator simulation. Such a model has not yet been implemented by the authors or their associates.

Tew, Roy C.↗
NASA NTRS2004

Modeling the Nonlinear, Strain Rate Dependent Deformation of Woven Ceramic Matrix Composites With Hydrostatic Stress Effects Included

An analysis method based on a deformation (as opposed to damage) approach has been developed to model the strain rate dependent, nonlinear deformation of woven ceramic matrix composites with a plain weave fiber architecture. In the developed model, the differences in the tension and compression response have also been considered. State variable based viscoplastic equations originally developed for metals have been modified to analyze the ceramic matrix composites. To account for the tension/compression asymmetry in the material, the effective stress and effective inelastic strain definitions have been modified. The equations have also been modified to account for the fact that in an orthotropic composite the in-plane shear stiffness is independent of the stiffness in the normal directions. The developed equations have been implemented into a commercially available transient dynamic finite element code, LS-DYNA, through the use of user defined subroutines (UMATs). The tensile, compressive, and shear deformation of a representative plain weave woven ceramic matrix composite are computed and compared to experimental results. The computed values correlate well to the experimental data, demonstrating the ability of the model to accurately compute the deformation response of woven ceramic matrix composites.

Goldberg, Robert K.↗
DOE OSTI2026

Data from a four-day long microcosm experiment addressing the destabilization of artificial mineral-associated organic matter by model root exudates embedded in a soil matrix from the Rocky Mountain Biological Laboratory (Gothic, CO, USA), 2019

This dataset provides data collected during a four-day long laboratory soil microcosm experiment testing the efficacy of root exudate-driven mineral-associated organic matter destabilization. This dataset contains four data files in comma-separate values (*.csv). The files provide the metadata and the experimental results on microbial respiration, MAOM-derived respiration, and sequential mineral-extractions. This data was used to produce the figures in Bölscher et al., 2026. The results of the experiment can be found in the open access article Bölscher et al., 2026 (https://doi.org/10.1016/j.soilbio.2026.110276). Abstract: Mineral-associated organic matter (MAOM) is often considered stable, but root exudates can destabilize MAOM via various pathways. Theory and model system studies suggest that direct MAOM destabilization by strong ligands, like oxalic acid, or reducing agents, like catechol, is more effective than indirect, microbial-mediated MAOM destabilization, stimulated by less reactive compounds like glucose. Here, we demonstrate that the presence of a soil matrix alters the efficacy of exudate-driven MAOM destabilization pathways. Glucose and catechol destabilized significantly greater amounts of MAOM from ferrihydrite and aluminum hydroxide (Al (OH)3) embedded in a soil matrix than oxalic acid. Our findings indicate that indirect, microbial-mediated MAOM destabilization may play a larger role than direct MAOM destabilization in soil environments.

Destabilization↗
NASA NTRS2000

Ultrasonic Assessment of Impact-Induced Damage and Microcracking in Polymer Matrix Composites

The main objective of this NASA FAR project is to conduct ultrasonic assessment of impact-induced damage and microcracking in polymer matrix composites at various temperatures. It is believed that the proposed study of impact damage assessment on polymer matrix composites will benefit several NASA's missions and current interests, such as ballistic impact testing of composite fan containment and high strain rate deformation modeling of polymer matrix composites. Currently, impact-induced delamination and fracture in 6061-T6 aluminum/cast acrylic sandwich plates adhered by epoxy were generated in an instrumented drop-weight impact machine. Although only a small dent was produced on the aluminum side when a hemispherical penetrator tup was dropped onto it from a couple of inches, a large ring of delamination at the interface was observed. The delamination damage was often accompanied by severe shattering in the acrylic substratum. Damage patterns in the acrylic layer include radial and ring cracks and, together with delamination at the interface, may cause peeling-off of acrylic material from the sandwich plate. Theory of stress-wave propagation can be used to explain these damage patterns. The impact tests were conducted at various temperatures. The results also show clearly that temperature effect is very important in impact damage. For pure cast acrylic nil-ductile transition (NDT) occurs between 185-195 F Excessive impact energy was dissipated into fracture energy when tested at temperature below this range or through plastic deformation when tested at temperature above the NDT temperature. Results from this study will be used as baseline data for studying fiber-metal laminates, such as GLARE and ARALL for advanced aeronautical and astronautical applications.

Liaw, Benjamin↗
NASA NTRS2004

Bayesian Approach to the Joint Inversion of Gravity and Magnetic Data, with Application to the Ismenius Area of Mars

This viewgraph presentation reviews a Bayesian approach to the inversion of gravity and magnetic data with specific application to the Ismenius Area of Mars. Many inverse problems encountered in geophysics and planetary science are well known to be non-unique (i.e. inversion of gravity the density structure of a body). In hopes of reducing the non-uniqueness of solutions, there has been interest in the joint analysis of data. An example is the joint inversion of gravity and magnetic data, with the assumption that the same physical anomalies generate both the observed magnetic and gravitational anomalies. In this talk, we formulate the joint analysis of different types of data in a Bayesian framework and apply the formalism to the inference of the density and remanent magnetization structure for a local region in the Ismenius area of Mars. The Bayesian approach allows prior information or constraints in the solutions to be incorporated in the inversion, with the "best" solutions those whose forward predictions most closely match the data while remaining consistent with assumed constraints. The application of this framework to the inversion of gravity and magnetic data on Mars reveals two typical challenges - the forward predictions of the data have a linear dependence on some of the quantities of interest, and non-linear dependence on others (termed the "linear" and "non-linear" variables, respectively). For observations with Gaussian noise, a Bayesian approach to inversion for "linear" variables reduces to a linear filtering problem, with an explicitly computable "error" matrix. However, for models whose forward predictions have non-linear dependencies, inference is no longer given by such a simple linear problem, and moreover, the uncertainty in the solution is no longer completely specified by a computable "error matrix". It is therefore important to develop methods for sampling from the full Bayesian posterior to provide a complete and statistically consistent picture of model uncertainty, and what has been learned from observations. We will discuss advanced numerical techniques, including Monte Carlo Markov

data analysis↗
DOE OSTI2022

Constructing Hubbard models for the hydrogen chain using sliced-basis density matrix renormalization group

Sliced-basis DMRG(sb-DMRG) is used to simulate a chain of hydrogen atoms and to construct low-energy effective Hubbard-like models. The downfolding procedure first involves a change of basis to a set of atom-centered Wannier functions constructed from the natural orbitals of the exact DMRG one-particle density matrix. The Wannier function model is then reduced to a fewer-parameter Hubbard-like model, whose parameters are determined by minimizing the expectation value of the Wannier Hamiltonian in the ground state of the Hubbard Hamiltonian. This indirect variational procedure not only yields compact and simple models for the hydrogen chain, but also allows us to explore the importance of constraints in the effective Hamiltonian, such as the restricting the range of the single-particle hopping and two-particle interactions, and to assess the reliability of more conventional downfolding. The entanglement entropy for a model's ground state, cut in the middle, is an important property determining the ability of DMRG and tensor networks to simulate the model, and here we study its variation with the range of the interactions. Counterintuitively, we find that shorter ranged interactions often have larger entanglement.

36 MATERIALS SCIENCE↗
NASA NTRS2004

A Method for Incorporating Changing Structural Characteristics Due to Propellant Mass Usage in a Launch Vehicle Ascent Simulation

Launch vehicles consume large quantities of propellant quickly, causing the mass properties and structural dynamics of the vehicle to change dramatically. Currently, structural load assessments account for this change with a large collection of structural models representing various propellant fill levels. This creates a large database of models complicating the delivery of reduced models and requiring extensive work for model changes. Presented here is a method to account for these mass changes in a more efficient manner. The method allows for the subtraction of propellant mass as the propellant is used in the simulation. This subtraction is done in the modal domain of the vehicle generalized model. Additional computation required is primarily for constructing the used propellant mass matrix from an initial propellant model and further matrix multiplications and subtractions. An additional eigenvalue solution is required to uncouple the new equations of motion; however, this is a much simplier calculation starting from a system that is already substantially uncoupled. The method was successfully tested in a simulation of Saturn V loads. Results from the method are compared to results from separate structural models for several propellant levels, showing excellent agreement. Further development to encompass more complicated propellant models, including slosh dynamics, is possible.

McGhee, D. S.↗
NASA NTRS1985

Short-wavelength buckling and shear failures for compression-loaded composite laminates

The short-wavelength buckling (or the microbuckling) and the interlaminar and inplane shear failures of multi-directional composite laminates loaded in uniaxial compression are investigated. A laminate model is presented that idealizes each lamina. The fibers in the lamina are modeled as a plate, and the matrix in the lamina is modeled as an elastic foundation. The out-of-plane w displacement for each plate is expressed as a trigonometric series in the half-wavelength of the mode shape for laminate short-wavelength buckling. Nonlinear strain-displacement relations are used. The model is applied to symmetric laminates having linear material behavior. The laminates are loaded in uniform end shortening and are simply supported. A linear analysis is used to determine the laminate stress, strain, and mode shape when short-wavelength buckling occurs. The equations for the laminate compressive stress at short-wavelength buckling are dominated by matrix contributions.

Shuart, M. J.↗
NASA NTRS1993

Evolution of damage and plasticity in titanium-based, fiber-reinforced composites

The inelastic deformation mechanisms were evaluated for a model titanium-based, fiber-reinforced composite: a beta titanium alloy (Ti-15V-3Al-3Cr-3Sn) reinforced with SiC (SCS-6) fibers. The primary emphasis of this article is to illustrate the sequence in which damage and plasticity evolved for this system. The mechanical responses and the results of detailed microstructural evaluations for the 0(8), 90(8), and +/- 45(2s) line oriented laminates are provided. It is shown that the characteristics of the reaction zone around the fiber play a very important role in the way damage and plasticity evolve, particularly in the microyield regime of deformation, and must be included in any realistic constitutive model. Fiber-matrix debonding was a major damage mode for the off-axis systems. The tension test results are also compared with the predictions of a few constitutive models.

Majumdar, B. S.↗