Search NASA⌕ Search

SEARCH · Search NASA

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

ICAN - Integrated composites analyzer

The ICAN computer program performs all the essential aspects of mechanics/analysis/design of multilayered fiber composites. Modular, open-ended and user friendly, the program can handle a variety of composite systems having one type of fiber and one matrix as constituents as well as intraply and interply hybrid composite systems. It can also simulate isotropic layers by considering a primary composite system with negligible fiber volume content. This feature is specifically useful in modeling thin interply matrix layers. Hygrothermal conditions and various combinations of in-plane and bending loads can also be considered. Usage of this code is illustrated with a sample input and the generated output. Some key features of output are stress concentration factors around a circular hole, locations of probable delamination, a summary of the laminate failure stress analysis, free edge stresses, microstresses and ply stress/strain influence coefficients. These features make ICAN a powerful, cost-effective tool to analyze/design fiber composite structures and components.

Murthy, P. L. N.↗
DOE OSTI2024

First measurement in a magnetic confinement fusion experiment of the H 3 + H 3 → He 5 + n intermediate two-body resonant reaction

We report on the first experimental measurements made at a magnetic confinement fusion device of the tritium(T)-tritium(T) reaction T + T → He 4 + 2 n indicating the presence of the intermediate two-body resonant reaction T + T → He 5 + n . During the second deuterium-tritium campaign (DTE2) at the Joint European Torus, measurements of fusion plasmas with high tritium concentrations, n T / ( n T + n D ) ≈ 0.99 , heated with tritium neutral beam injection, were performed using the neutron time-of-flight (TOF) spectrometer TOFOR. We detect a peak in the neutron emission TOF spectrum consistent with the two-body resonant reaction. The TT neutron emission energy spectrum is modeled using an R -matrix framework where the distributions of the most likely model parameters given our experimental TOF data are determined utilizing a Markov chain Monte Carlo approach. We compare our best estimate of the T + T neutron emission energy spectrum with results obtained at inertial confinement fusion experiments at the OMEGA facility and find a spectral shape that is consistent with the energy dependency in the neutron spectrum observed at OMEGA. Published by the American Physical Society 2024

Physics↗
NASA NTRS1994

Mathematical modeling of a class of multibody flexible spacecraft structures

A mathematical model for a general multibody flexible spacecraft is obtained. The generic spacecraft considered consists of a flexible central body to which a number of flexible multibody structures are attached. The coordinate systems used in the derivation allow effective decoupling of the translational motion of the entire spacecraft from its rotational motion about its center of mass. The derivation assumes that the deformations in the bodies are only due to elastic motions. The dynamic model derived is a closed-form vector-matrix differential equation. The model developed can be used for analysis and simulation of many realistic spacecraft configurations.

Kelkar, Atul, G.↗
NASA NTRS2005

Simulation of the ATIC-2 Silicon Matrix for Protons and Helium GCR Primaries at 0.3, 10, and 25 TeV/Nucleon

The energy deposition distribution for protons and helium galactic cosmic ray primaries at 0.3, 10, and 25 TeV/nucleon in the ATIC-2 silicon matrix detector are simulated with GEANT4. The GEANT3 geometrical model of ATIC developed by the University of Maryland was combined with a GEANT4 application developed for the Deep Space Test Bed (DSTB) detector package. The new code included relatively minor modifications to completely describe the ATIC materials and a more detailed model of the Silicon Matrix detector. For this analysis all particles were started as a unidirectional beam at a single point near the center of the Silicon Matrix front surface. The point was selected such that each primary passed through at least two of the overlapping silicon pixels.

Watts, J.↗
NASA NTRS2000

Implementation of Laminate Theory Into Strain Rate Dependent Micromechanics Analysis of Polymer Matrix Composites

A research program is in progress to develop strain rate dependent deformation and failure models for the analysis of polymer matrix composites subject to impact loads. Previously, strain rate dependent inelastic constitutive equations developed to model the polymer matrix were implemented into a mechanics of materials based micromechanics method. In the current work, the computation of the effective inelastic strain in the micromechanics model was modified to fully incorporate the Poisson effect. The micromechanics equations were also combined with classical laminate theory to enable the analysis of symmetric multilayered laminates subject to in-plane loading. A quasi-incremental trapezoidal integration method was implemented to integrate the constitutive equations within the laminate theory. Verification studies were conducted using an AS4/PEEK composite using a variety of laminate configurations and strain rates. The predicted results compared well with experimentally obtained values.

Goldberg, Robert K.↗
DOE OSTI2025

Architecture-Aware Models of AI Engines for High-Performance Matrix Matrix Multiplication

The AI Engine (AIE) architecture, available in systems from mobile SoCs to server-class FPGAs, aims to efficiently execute AI/ML tasks through a two-dimensional array of compute tiles. Previous work on AIEs has explored different approaches to mapping computation across spatial arrays, but the compute kernel running on each tile has not been the focus. Additionally, the AIE-ML architecture introduces memory tiles and omits programmable logic, requiring new approaches to staging and moving data throughout the array. In this work we update analytical models developed for CPUs to produce the design of high performance kernels while introducing new model considerations such as memory structure, throughput, and latency as required by the AIE hardware. We evaluate our models by developing AIE-ML kernels for matrix multiplication in low-precision data types showing performance up to 95% of compute peak for the kernel when data resides in local memory and above 90% of compute peak when data resides in main memory.

Binder, Elliott D. [Carnegie Mellon University, Pi↗
NASA NTRS2015

MATRIX-VBS Condensing Organic Aerosols in an Aerosol Microphysics Model

The condensation of organic aerosols is represented in a newly developed box-model scheme, where its effect on the growth and composition of particles are examined. We implemented the volatility-basis set (VBS) framework into the aerosol mixing state resolving microphysical scheme Multiconfiguration Aerosol TRacker of mIXing state (MATRIX). This new scheme is unique and advances the representation of organic aerosols in models in that, contrary to the traditional treatment of organic aerosols as non-volatile in most climate models and in the original version of MATRIX, this new scheme treats them as semi-volatile. Such treatment is important because low-volatility organics contribute significantly to the growth of particles. The new scheme includes several classes of semi-volatile organic compounds from the VBS framework that can partition among aerosol populations in MATRIX, thus representing the growth of particles via condensation of low volatility organic vapors. Results from test cases representing Mexico City and a Finish forrest condistions show good representation of the time evolutions of concentration for VBS species in the gas phase and in the condensed particulate phase. Emitted semi-volatile primary organic aerosols evaporate almost completely in the high volatile range, and they condense more efficiently in the low volatility range.

Aerosols↗
NASA NTRS1972

Determination of stability derivatives from flight data using a Newton-Raphson minimization technique

A modified Newton-Raphson or quasilinearization minimization technique for determining stability derivatives from flight data was developed and compared with simple-equations, analog-matching, least-squares, and Shinbrot methods of analysis. For the data analyzed, the solutions computed by using the estimates obtained from the Newton-Raphson technique fit the data and determined coefficients adequately. A further modification to include a priori information was found to be useful. A model statistically similar to the flight data was analyzed using the same methods (excluding analog matching), and the Newton-Raphson technique was found to yield superior estimates. An approximate Cramer-Rao bound was compared with the error covariance matrix of the model and was found to provide information about the reliability of the individual estimates obtained. The technique was successfully applied to data obtained from a light airplane, a large supersonic airplane, and a lifting body vehicle. It was shown that the reliability of the estimates of a given coefficient obtained from these vehicles depends upon the data analyzed.

Iliff, K. W.↗
NASA NTRS2001

Implementation of Fiber Substructuring Into Strain Rate Dependent Micromechanics Analysis of Polymer Matrix Composites

A research program is in progress to develop strain rate dependent deformation and failure models for the analysis of polymer matrix composites subject to impact loads. Previously, strain rate dependent inelastic constitutive equations developed to model the polymer matrix were incorporated into a mechanics of materials based micromechanics method. In the current work, the micromechanics method is revised such that the composite unit cell is divided into a number of slices. Micromechanics equations are then developed for each slice, with laminate theory applied to determine the elastic properties, effective stresses and effective inelastic strains for the unit cell. Verification studies are conducted using two representative polymer matrix composites with a nonlinear, strain rate dependent deformation response. The computed results compare well to experimentally obtained values.

Goldberg, Robert K.↗
DOE OSTI2025

Uncertainty estimation of bifurcated solutions in the Rayleigh–Bénard problem for advanced nuclear reactors applications

Multiphysics models of nuclear reactors frequently comprise nonlinear systems of equations. The nonlinear nature of these models could lead to solution bifurcations, where a small change in a certain parameter, e.g., the thermophysical properties of the coolant, can lead to a sudden change in the system’s behavior. At the point in parameter space where this happens, called a critical point, the Jacobian matrix of the model’s nonlinear operator becomes singular potentially permitting multiple solutions to coexist. In this paper, we perform uncertainty estimation (UE) in a parameter range that includes bifurcated solutions within the context of Rayleigh–Bénard problem. We perform this analysis assuming uncertain temperature difference, and tilt angle for the iterative solution algorithm with a unit Prandtl number (Pr = 1). Also, we perform this analysis under uncertain thermophysical properties for both FLiBe molten salt and liquid sodium as working fluid. We deploy two approaches to compute statistical moments for the resulting distributions of selected flow-field variables. The first approach is the blind computation of the mean and the standard deviation without any consideration of solution bifurcation, while the second approach utilizes k-means clustering to cluster each branch’s solutions together and compute separate statistical moments for each branch. The statistical distributions are obtained by perturbing the selected parameters about nominal values that correspond to a solution on one of the valid branches, and that solution is used as initial guess for the iterative solution algorithm. We found that perturbation of any parameter when its nominal value is close to its critical point always leads to branch jumping, i.e., the iterations converge to a solution on a branch different from the branch of the initial guess. This produces a statistical ensemble comprised of fundamentally different solutions leading to wrong mean values and uncertainty estimates, whereas clustering provides an efficient way to deal with this type of computation. This work is important for developing Gen IV nuclear systems because many of these systems rely on natural convection for cooling especially in accident conditions.

97 - MATHEMATICS AND COMPUTING↗
NASA NTRS2003

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

The main objective of this NASA Faculty Awards for Research (FAR) project is to conduct ultrasonic assessment of impact-induced damage and microcracking in fiber-metal laminated (FML) composites at various temperatures. It is believed that the proposed study of impact damage assessment on FML 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. 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 2024-T3 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.

Gyekanyesi, John↗
NASA NTRS2006

An Analytical Thermal Model for Autonomous Soaring Research

A viewgraph presentation describing an analytical thermal model used to enable research on autonomous soaring for a small UAV aircraft is given. The topics include: 1) Purpose; 2) Approach; 3) SURFRAD Data; 4) Convective Layer Thickness; 5) Surface Heat Budget; 6) Surface Virtual Potential Temperature Flux; 7) Convective Scaling Velocity; 8) Other Calculations; 9) Yearly trends; 10) Scale Factors; 11) Scale Factor Test Matrix; 12) Statistical Model; 13) Updraft Strength Calculation; 14) Updraft Diameter; 15) Updraft Shape; 16) Smoothed Updraft Shape; 17) Updraft Spacing; 18) Environment Sink; 19) Updraft Lifespan; 20) Autonomous Soaring Research; 21) Planned Flight Test; and 22) Mixing Ratio.

Allen, Michael↗
NASA NTRS1982

Evaluation of trends in wheat yield models

Trend terms in models for wheat yield in the U.S. Great Plains for the years 1932 to 1976 are evaluated. The subset of meteorological variables yielding the largest adjusted R(2) is selected using the method of leaps and bounds. Latent root regression is used to eliminate multicollinearities, and generalized ridge regression is used to introduce bias to provide stability in the data matrix. The regression model used provides for two trends in each of two models: a dependent model in which the trend line is piece-wise continuous, and an independent model in which the trend line is discontinuous at the year of the slope change. It was found that the trend lines best describing the wheat yields consisted of combinations of increasing, decreasing, and constant trend: four combinations for the dependent model and seven for the independent model.

Ferguson, M. C.↗
NASA NTRS2013

Piloted Simulation Assessment of the Impact of Flexible Structures on Handling Qualities of Generic Supersonic Aircraft

The NASA Langley Research Center Cockpit Motion Facility (CMF) was used to conduct a piloted simulation assessment of the impact of flexible structures on flying qualities. The CMF was used because of its relatively high bandwidth, six degree-of-freedom motion capability. Previous studies assessed and attempted to mitigate the effects of multiple dynamic aeroservoelastic modes (DASE). Those results indicated problems existed, but the specific cause and effect was difficult to ascertain. The goal of this study was to identify specific DASE frequencies, damping ratios, and gains that cause degradation in handling qualities. A generic aircraft simulation was developed and designed to have Cooper-Harper Level 1 handling qualities when flown without DASE models. A test matrix of thirty-six DASE modes was implemented. The modes had frequencies ranging from 1 to 3.5 Hz and were applied to each axis independently. Each mode consisted of a single axis, frequency, damping, and gain, and was evaluated individually by six subject pilots with test pilot backgrounds. Analysis completed to date suggests that a number of the DASE models evaluated degrade the handling qualities of this class of aircraft to an uncontrollable condition.

Stringer, Mary T.↗
NASA NTRS1991

Analysis of metal-matrix composite structures. I - Micromechanics constitutive theory. II - Laminate analyses

The micromechanical constitutive theory is used to examine the nonlinear behavior of continuous-fiber-reinforced metal-matrix composite structures. Effective lamina constitutive relations based on the Abouli micromechanics theory are presented. The inelastic matrix behavior is modeled by the unified viscoplasticity theory of Bodner and Partom. The laminate constitutive relations are incorporated into a first-order deformation plate theory. The resulting boundary value problem is solved by utilizing the finite element method. Attention is also given to computational aspects of the numerical solution, including the temporal integration of the inelastic strains and the spatial integration of bending moments. Numerical results the nonlinear response of metal matrix composites subjected to extensional and bending loads are presented.

Arenburg, R. T.↗
DOE OSTI2025

Predicting Flow in Fracture Networks With Quantum Algorithms

Uncertainty quantification plays a crucial role in the modeling of subsurface flow. For instance, uncertainties in the properties of geologic fracture networks significantly impact flow, requiring numerous simulations to accurately estimate quantities of interest. However, each simulation is computationally expensive because it requires solving a large linear system to capture features that involve both small and large fractures. An example is in percolation, where the interaction of many small fractures (which cumulatively can have a large surface area) with the rock matrix must be modeled precisely. Quantum computing is an emerging tool with the potential to address this issue. Quantum algorithms offer a significant speedup in solving linear systems, achieving efficiencies that are challenging to match with classical approaches. These classical approaches include direct solvers, such as LU decomposition, and iterative methods, notably preconditioned conjugate gradient, commonly used in subsurface modeling to solve large sparse systems. However, applying quantum algorithms to geologic fracture flow requires careful attention to algorithmic and problem-specific constraints to fully realize this quantum advantage. In this work we describe a quantum algorithm for generalized Monte Carlo applications with a quadratic speedup over the classical approaches which can be combined with the quantum speedup, currently under investigation, for solving quantum linear systems for subsurface flow. We show that for quantum algorithms the computational cost of estimating a quantity of interest for a statistical ensemble of networks is roughly the same as that of a single realization, essentially implying that one can get uncertainty quantification for free.

58 GEOSCIENCES↗
NASA NTRS2004

Spectroscopic Determination of the Physical Conditions in Hot Optically Thin Sources

Astrophysical Plasma Emission Database - APED consists of atomic data, primarily theoretical, needed to calculate X-ray through optical emission spectra of hot thermal plasmas. These data are supplemented by experimental values, such as laboratory plasma wave- lengths, and are validated by experiment only rarely. Thus the comparison of predicted spectra to astrophysical spectra serves as invaluable feedback on the quality of the models. Following up on unfavorable comparisons, we have communicated atomic data needs to atomic theorists and experimental plasma physicists so that improvements can be made to the database. The database is now reasonably complete from approx. 1 to 30 A, and from - 90 to 180 A, with significant gaps remaining in the range 30 to 90 A. Astrophysical Plasma Emission Code - APEC uses the APED data as inputs and computes level populations and line emissivities using a rate matrix solver. Tabulated models for a fine grid of temperatures and densities can be applied to observational data. The current public version of APEC output uses an input ionization balance model, assuming CIE; however, we have a working version to calculate the ionization state of the gas. We will begin testing non-equilibrium ionization (NEI) models over the next year - we are now finalizing the atomic database needed to perform the ionization state calculations. At that time, the code will be made public. We have developed optimized error codes, additional testing protocols and substantial documentation, in preparation for this public release. Secondarily, we have begun to investigate the use of APED and APEC for X-ray photoionized plasma, beginning with opacity modeling (e.g. warm absorbers). Application to bright AGN Chandra grating spectra suggests that the accurate wavelengths in APED are of great importance to fitting the data. We are exploring additional subroutines to APEC to extend its usefulness for X-ray photoionized plasma, e.g. to include collisional or photo-excitation of metastable levels for absorption.

Oliversen, Ronald J.↗
NASA NTRS1999

Optimal Frequency-Domain System Realization with Weighting

Several approaches are presented to identify an experimental system model directly from frequency response data. The formulation uses a matrix-fraction description as the model structure. Frequency weighting such as exponential weighting is introduced to solve a weighted least-squares problem to obtain the coefficient matrices for the matrix-fraction description. A multi-variable state-space model can then be formed using the coefficient matrices of the matrix-fraction description. Three different approaches are introduced to fine-tune the model using nonlinear programming methods to minimize the desired cost function. The first method uses an eigenvalue assignment technique to reassign a subset of system poles to improve the identified model. The second method deals with the model in the real Schur or modal form, reassigns a subset of system poles, and adjusts the columns (rows) of the input (output) influence matrix using a nonlinear optimizer. The third method also optimizes a subset of poles, but the input and output influence matrices are refined at every optimization step through least-squares procedures.

Juang, Jer-Nan↗