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 685 records · Page 38

Autonomous satellite navigation methods using the Global Positioning Satellite System

This investigation considers the problem of autonomous satellite navigation using the NAVSTAR Global Positioning System (GPS). The major topics covered include the design, implementation, and validation of onboard navigation filter algorithms by means of computer simulations. The primary errors that the navigation filter design must minimize are computational effects and modeling inaccuracies due to limited capability of the onboard computer. The minimization of the effect of these errors is attained by applying the sequential extended Kalman filter using a factored covariance implementation with Q-matrix or dynamical model compensations. Peformance evaluation of the navigation filter design is carried out using both the CDC Cyber 170/750 computer and the PDP-11/60 computer. The results are obtained assuming the Phase I GPS constellation, consisting of six satellites, and a Landsat-D type spacecraft as the model for the user satellite orbit.

Murata, M.↗

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

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↗

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

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

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

Topological terms with qubit regularization and relativistic quantum circuits

Qubit regularization provides a rich framework to explore quantum field theories. The freedom to choose how the important symmetries of the theory are embedded in the qubit regularization scheme allows us to construct new lattice models with rich phase diagrams. Some of the phases can contain topological terms which lead to critical phases. In this work we introduce and study the SU(3)-F qubit regularization scheme to embed the SO(3) spin-symmetry. We argue that qubit models in this regularization scheme contain several phases including a critical phase which describes the k = 1 Wess-Zumino-Witten (WZW) conformal field theory (CFT) at long distances, and two massive phases one of which is trvially gapped and the other which breaks the lattice translation symmetry. We construct a simple space-time Euclidean lattice model with a single coupling U and study it using the Monte Carlo method. We show the model has a critical phase at small U and a trivially massive phase at large U with a first order transition separating the two. Another feature of our model is that it is symmetric under space-time rotations, which means the temporal and spatial lattice spacing are connected to each other. The unitary time evolution operator obtained by a Wick rotation of the transfer matrix of our model can help us compute the physics of the k = 1 WZW CFT in real time without the need for tuning the temporal lattice spacing to zero. We use this idea to introduce the concept of a relativistic quantum circuit on a discrete space-time lattice.

Bhattacharya, Tanmoy↗

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↗

Crystal Plasticity Model of Reactor Pressure Vessel Embrittlement in GRIZZLY

The integrity of reactor pressure vessels (RPVs) is of utmost importance to ensure safe operation of nuclear reactors under extended lifetime. Microstructure-scale models at various length and time scales, coupled concurrently or through homogenization methods, can play a crucial role in understanding and quantifying irradiation-induced defect production, growth and their influence on mechanical behavior of RPV steels. A multi-scale approach, involving atomistic, meso- and engineering-scale models, is currently being pursued within the GRIZZLY project to understand and quantify irradiation-induced embrittlement of RPV steels. Within this framework, a dislocation-density based crystal plasticity model has been developed in GRIZZLY that captures the effect of irradiation-induced defects on the flow stress behavior and is presented in this report. The present formulation accounts for the interaction between self-interstitial loops and matrix dislocations. The model predictions have been validated with experiments and dislocation dynamics simulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

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↗

On the intrinsic diversity of type II-plateau supernovae

Hydrogen-rich Type II-Plateau supernovae (SNe) exhibit correlations between the plateau luminosity L{sub pl}, the nickel mass M{sub Ni}, the explosion energy E{sub exp}, and the ejecta mass M{sub ej}. Using our global, self-consistent, multi-band model of nearby well-observed SNe, we find that the covariances of these quantities are strong and that the confidence ellipsoids are oriented in the direction of the correlations, which reduces their significance. By proper treatment of the covariance matrix of the model, we discover a significant intrinsic width to the correlations between L{sub pl}, E{sub exp} and M{sub Ni}, where the uncertainties due to the distance and the extinction dominate. For fixed E{sub exp}, the spread in M{sub Ni} is about 0.25 dex, which we attribute to the differences in the progenitor internal structure. We argue that the effects of incomplete γ-ray trapping are not important in our sample. Similarly, the physics of the Type II-Plateau SN light curves leads to inherently degenerate estimates of E{sub exp} and M{sub ej}, which makes their observed correlation weak. Ignoring the covariances of SN parameters or the intrinsic width of the correlations causes significant biases in the slopes of the fitted relations. Our results imply that Type II-Plateau SN explosions are not described by a single physical parameter or a simple one-dimensional trajectory through the parameter space, but instead reflect the diversity of the core and surface properties of their progenitors. We discuss the implications for the physics of the explosion mechanism and possible future observational constraints.

79 ASTRONOMY AND ASTROPHYSICS↗

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

Bounds on the number of bound states in the transfer matrix spectrum for some weakly correlated lattice models

We consider the interaction of particles in weakly correlated lattice quantum field theories. In the imaginary time functional integral formulation of these theories there is a relative coordinate lattice Schroedinger operator H which approximately describes the interaction of these particles. Scalar and vector spin, QCD and Gross-Neveu models are included in these theories. In the weakly correlated regime H=H{sub o}+W where H{sub o}=-{gamma}{Delta}{sub l}, 0 < {gamma} Much-Less-Than 1 and {Delta}{sub l} is the d-dimensional lattice Laplacian: {gamma}={beta}, the inverse temperature for spin systems and {gamma}={kappa}{sup 3} where {kappa} is the hopping parameter for QCD. W is a self-adjoint potential operator which may have non-local contributions but obeys the bound Double-Vertical-Line W(x, y) Double-Vertical-Line Less-Than-Or-Slanted-Equal-To cexp ( -a( Double-Vertical-Line x Double-Vertical-Line + Double-Vertical-Line y Double-Vertical-Line )), a large: exp-a={beta}/{beta}{sub o}{sup (1/2)}({kappa}/{kappa}{sub o}) for spin (QCD) models. H{sub o}, W, and H act in l{sub 2}(Z{sup d}), d Greater-Than-Or-Slanted-Equal-To 1. The spectrum of H below zero is known to be discrete and we obtain bounds on the number of states below zero. This number depends on the short range properties of W, i.e., the long range tail does not increase the number of states.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Open quantum systems and random matrix theory

A simple model for open quantum systems is analyzed with RMT. The system is coupled to the continuum in a minimal way. In this paper we see the effect of opening the system on the level statistics, in particular the level spacing, width distribution and Δ{sub 3}(L) statistic are examined as a function of the strength of this coupling. The usual super-radiant state is observed, and it is seen that as it is formed, the level spacing and Δ{sub 3}(L) statistic exhibit the signatures of missed levels.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Large covariance matrices: accurate models without mocks

Abstract Covariance matrix estimation is a persistent challenge for cosmology. We focus on a class of model covariance matrices that can be generated with high accuracy and precision, using a tiny fraction of the computational resources that would be required to achieve comparably precise covariance matrices using mock catalogues. In previous work, the free parameters in these models were determined using sample covariance matrices computed using a large number of mocks, but we demonstrate that those parameters can be estimated consistently and with good precision by applying jackknife methods to a single survey volume. This enables model covariance matrices that are calibrated from data alone, with no reference to mocks.

79 ASTRONOMY AND ASTROPHYSICS↗