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 1,765 records · Page 98

Fundamental Study on Sorption Behavior of Europium to Mica Flakes - 19116

It is well known that biotite exists with flake form in actual underground environment, dominating sorption of cationic radionuclides in a geological disposal system consisting of crystalline rock. The behavior of sorption of radionuclides will differ depending on the difference in shape, while most of the previous studies have used powder biotite. This study focused on the influence of flake structure of phyllosilicate minerals such as biotite on the diffusion and sorption processes of radionuclides, examining the sorption behaviors of Eu on biotite, phlogopite and muscovite in flake form. These flakes (5 mm x 6 mm) were prepared for the sorption experiments and the initial concentrations of the Eu ions were set to 1 mM. The sorption experiments were conducted by mixing 0.4 g of biotite flakes and 4 mL at pH 5 for 7 days. As a result, the sorption amount of Eu increased little by little with internal diffusion in all minerals, respectively. The diffusion process of Eu into the layer structure could be also observed by secondary ion-microprobe mass spectrometry (SIMS) analysis. The apparent diffusion coefficient of Eu into each flake was estimated by fitting a 2-D diffusion model to the sorption data. The estimated diffusion values were clearly smaller than that of matrix diffusion in granite. Hence it was suggested that the so-called Kd-model, which is a local equilibrium model, is not always applicable to the estimation of the retardation effect (i.e., migration rate) of radionuclides in rock matrix. In other words, the migration rate might depend also on the diffusion process of radionuclide in the inside of phyllosilicate minerals in addition of the macro diffusion process in rock matrix. Besides, sorption amount and diffusion coefficient of biotite, phlogopite and muscovite were largely different. Not only the tetrahedral structure adjacent to the potassium but also the octahedral structure might affect the restraint of the interlayer ion (K ions) and the sorption behavior. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Rotor dynamic state and parameter identification from simulated forward flight transients, part 1

State and parameter identifications from simulated forward flight blade flapping measurements are presented. The transients were excited by progressing cyclic pitch stirring or by hub stirring with constant stirring acceleration. Rotor dynamic inflow models of varying degree of sophistication were used from a one parameter inflow model (equivalent Lock number) to an eight parameter inflow model. The maximum likelihood method with assumed fixed measurement error covariance matrix was applied. The rotor system equations for both fixed hub and tilting hub are given. The identified models were verified by comparing true responses with predicted responses. An optimum utilization of the simulated measurement data can be defined. From the numerical results it can be anticipated that brief periods of either accelerated cyclic pitch stirring or of hub stirring are sufficient to extract with adequate accuracy up to 8 rotor dynamic inflow parameters plus the blade Lock number from the transients.

Hohenemser, K. H.↗

Radiative transfer in realistic planetary atmospheres

Some 40 publications that appeared in scientific journals from 1973 to 1981 as well as 45 scientific reports issued during the grant period are listed by title. Topics cover the development of a matrix operator theory of radiative transfer which made possible the exact model calculations of the radiance as a function of height in planetary atmospheres; calculation of the Mie phase matrix for various types of particles as well as for radiance and polarization in planetary atmospheres; analysis of high dispersion spectroscopic observations of Venus; calculation of curves of growth for Venus; the development of a theory for calculating radiative transfer in spherical shell atmospheres; investigations of zonal winds on Venus; and examination of Rayleigh scattering.

Plass, G. N.↗

Supernova-Explosion Mechanism Involving Neutrinos

A major part of the energy released upon the gravitational collapse of massive-star cores is carried away by neutrinos. Neutrinos play a crucial role in collapsing supernovae (SNe). At the present time, mathematical models of core-collapse SNe are based on multidimensional gas dynamics and thermonuclear reactions, whereas the neutrino transport is frequently treated in simplified way. An accurate analysis of neutrinos in a spherically symmetric gravitational collapse is performed on the basis of Boltzmann kinetic equations including all weak-interaction reactions with exact quantum-mechanical matrix elements. The role of multidimensional effects is studied bymeans of multidimensional gas dynamics allowing for the neutrino transport via diffusion treated by employing flux limiters. The possibility of largescale convection, which is of interest both from the point of view of explaining a type II supernova (SN) and from the point of view of implementing an experiment aimed at detecting possible energetic (≳10 MeV) neutrinos from an SN, is discussed. Thermonuclear burning leads to the explosion of a type I SN. A hot central region and the subsequent large-scale convection may also play an important role in the SN mechanism. If neutrinos and convection play a key role for a type II SN, then, in order to explain gamma radiation from product radioactive elements, convection is of importance in the case of SNe belonging to both types. In addition, convection may be important for bright type I SNe. Original methods are presented for multidimensional gas dynamics involving thermonuclear burning and for multitemperature gas dynamics involving radiative transfer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Supersymmetry in the Jaynes-Cummings model

A review is presented of the Darboux method and its relation to the supersymmetric quantum mechanics, together with the embedding of a n-dimensional scalar Hamiltonian into a supersymmetric matrix. It is also shown that the Jaynes-Cummings model, with or without rotating wave approximation, admit a supersymmetric quantum mechanics description.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Vector-Ordering Filter Procedure for Data Reduction

The vector-ordering filter (VOF) technique involves a procedure for sampling a large population of data vectors to select a subset of data vectors that fully characterize the state space of the large population. The VOF technique enables a large reduction of the volume of data that must be handled in the automated monitoring system and method discussed in the two immediately preceding articles. In so doing, the VOF technique enables the development of data-driven mathematical models of a monitored asset from sets of data that would otherwise exceed the memory capacities of conventional engineering computers. Data-driven mathematical models have been shown to offer high fidelity for purposes of control and monitoring of assets. In practice, a collection of asset-operating observations is acquired with the intention that the collection contain observations characteristic of the full dynamic range of operation of the asset. Often, such a collection contains an extremely large number of observations, many of which are redundant. The VOF technique fills the need for a means to extract, from the original collection of observational data, a reduced data matrix that excludes redundant data while maintaining the full statistical character and dynamic range of the original data. The reduced data matrix can then be used as the input data for development of a mathematical model of the monitored asset, or as training data for a neural-network substitute for an explicit mathematical model of the asset. Alternatively, the reduced data matrix can, itself, be used directly as a mathematical model of the monitored asset, as is commonly done in multivariate state-estimation techniques. The original data are collected from the asset over a range of operating states and are put in matrix form. Each column vector in the original data matrix represents the signal values acquired at a particular operational state of the asset. Thus, the number of columns of the original data matrix equals the number of observed states and the number of rows in this matrix equals the number of signals acquired at each observation. In the VOF technique, one extracts the reduced data matrix from the original data matrix through the selection of a representative subset of the column (state) vectors.

Bickford, Randall L.↗

Studies in Quantum Field Theory and Astroparticle Physics

OVERVIEW OF THE RESEARCH PLAN: Task 1: Theoretical Research in High Energy Physics Prof. Alford investigated the phases and properties of matter at the highest densities achieved in nature. His main focus was on quark matter and its manifestations and signatures in neutron stars. His work is directly relevant to both the Cosmic and the Energy Frontiers. Prof. Bender developed PT-symmetric quantum theory, which is a complex generalization of conventional Hermitian quantum theory. During the grant period, he constructed and analyzed new PT-symmetric field theoretic models that exhibit symmetry breaking. He also worked on other fundamental problems in theoretical physics. His work falls under the Intensity Frontier classification. Prof. Bernard continued to work on lattice gauge theory and associated issues in chiral perturbation theory. His research focused primarily on the calculation of weak matrix elements relevant to tests of the standard model and the search for new physics at the Intensity Frontier. This work included both direct numerical simulations in QCD and the development of theoretical tools to analyze the simulations and facilitate the comparison with Intensity-Frontier experiments. Prof. Ferrer continued working on dark matter and other topics in theoretical particle cosmology. His research focused on weakly interacting massive particles relevant to both the Cosmic and the Energy Frontiers, and on dark photons relevant to the Intensity Frontier. Prof. Ogilvie continued his work on lattice gauge theory and the phase structure of gauge theories. His research on the properties of confined, conformal and Higgs phases is relevant to the Energy Frontier, as well as addressing fundamental issues of broad theoretical interest. Task 2: Experimental Research at the Cosmic Frontier in High Energy Physics. Over the three year grant period, Profs. Buckley and Krawczynski continued research under the Cosmic Frontier. Prof. Buckley's group worked on two Cosmic Frontier projects (VERITAS and LZ), and conducted research in dark matter science and other areas of particle astrophysics. Buckley’s group worked on the conceptual design of the next generation ground-based gamma-ray experiment (CTA) and on the development of camera electronics for a 2-mirror SCT prototype telescope. Buckley's group also continued work on GaN photodetector detector development as part of the Instrumentation Frontier (Detector R&D) including work with the LAPPD group on the development robust amorphous photocathodes for very large area MCP-based photodetectors and ongoing development of deep UV photosensors for liquid Ar and Xenon Direct Dark matter detection. Prof. Krawczynski's group worked on two Cosmic Frontier topics: astrophysical tests of Lorentz Invariance and CPT symmetries at the Planck energy scale with the VERITAS gamma-ray experiment, and the search for dark matter particles with indirect detection experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Theoretical Development of an Orthotropic Elasto-Plastic Generalized Composite Material Model

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites is becoming critical as these materials are gaining increased usage in the aerospace and automotive industries. While there are several composite material models currently available within LS-DYNA (Registered), there are several features that have been identified that could improve the predictive capability of a composite model. To address these needs, a combined plasticity and damage model suitable for use with both solid and shell elements is being developed and is being implemented into LS-DYNA as MAT_213. A key feature of the improved material model is the use of tabulated stress-strain data in a variety of coordinate directions to fully define the stress-strain response of the material. To date, the model development efforts have focused on creating the plasticity portion of the model. The Tsai-Wu composite failure model has been generalized and extended to a strain-hardening based orthotropic material model with a non-associative flow rule. The coefficients of the yield function, and the stresses to be used in both the yield function and the flow rule, are computed based on the input stress-strain curves using the effective plastic strain as the tracking variable. The coefficients in the flow rule are computed based on the obtained stress-strain data. The developed material model is suitable for implementation within LS-DYNA for use in analyzing the nonlinear response of polymer composites.

Plasticity↗

Relationship between the structure and electrical characteristics of diamond-like carbon films

To elucidate the relationship between the structure and the electrical characteristics of diamond-like carbon (DLC) films, DLC films were synthesized in a well-controlled glow discharge with the aid of photoelectrons in an argon/methane atmosphere. The dielectric constant and breakdown strength of the films exhibited opposite behaviors, depending on the total pressure during the synthesis. The product of these two values decreased monotonically as the pressure increased. The Raman spectra were analyzed with a Voigt-type formula. Based on the results, the authors propose the “sp{sup 2} cluster model” for the DLC structure. This model consists of conductive clusters of sp{sup 2} carbons surrounded by a dielectric matrix sea of sp{sup 2} carbon, sp{sup 3} carbon, and hydrogen, and indicates that the dielectric constant of the whole DLC film is determined by the balance between the dielectric constant of the matrix and the total size of the clusters, while the breakdown strength is determined by the reciprocal of the cluster size. The model suggests that a high-κ DLC film can be synthesized at a middle pressure and consists of well-grown sp{sup 2} clusters and a dense matrix. A low-κ DLC film can be synthesized both at low and high pressures. The sp{sup 2} cluster model explains that a low-κ DLC film synthesized at low pressure consists of a dense matrix and a low density of sp{sup 2} clusters, and exhibits a high breakdown strength. On the other hand, a low-κ film synthesized at high pressure consists of a coarse matrix and a high density of clusters and exhibits a low breakdown strength.

36 MATERIALS SCIENCE↗

Distributed-Memory Sparse Deep Neural Network Inference Using Global Arrays

Partitioned Global Address Space (PGAS) models exhibit tremendous promise in developing efficient and productive distributed-memory parallel applications. They have been used extensively in scientific computations due to conveniently offering a ``shared-memory''-like model and convenient interfaces that separate communication with synchronization. Traditionally, PGAS communication models have been applied to dense/contiguously distributed data, but most modern applications depict varied levels of sparsity. Existing PGAS models require certain adaptations to support distributed sparse computations, since associated computations often require matrix arithmetic, in addition to data movement. The Global Arrays toolkit from Pacific Northwest National Laboratory (PNNL) is one of the earliest PGAS models to combine one-sided data communication and distributed matrix operations and is still used in the popular NWChem quantum chemistry suite. Recently, we have expanded the Global Arrays toolkit to support common sparse operations, like sparse matrix-dense matrix multiplies (SpMM), sparse matrix-sparse matrix multiplication (SpGEMM) and Sampled Dense-Dense Matrix Multiplication (SDDMM). As it turns out, these operations are the bedrock of sparse Deep Learning (DL); sparse deep neural networks and Graph Neural Networks (GNNs) have gained increasing attention recently in achieving speedups on training and inference with reduced memory footprints. Unlike scientific applications in High Performance Computing (HPC), modern (distributed-memory capable) DL toolkits often rely on non-standardized and closed-source vendor software optimizations, creating challenges in software-hardware co-design at scale. Our goal is to support a variety of distributed-memory sparse matrix operations and helper functions in the newly created Sparse Global Arrays (SGA), such that it is possible to build portable and productive Machine Learning scenarios for algorithm/software and hardware codesign purposes. Contemporary data-parallel schemes for training/inference are undergoing a major overhaul since model replication limits scalability and causes resource inefficiencies. As such, we have adopted tensor parallelism in decomposing the model and inputs, to mitigate memory issues. Current implementation is built on top of MPI and uses CPUs to maximize the portability across the platforms.

Distributed computing, machine learning↗

The Influence of SLD Drop Size Distributions on Ice Accretion in the NASA Icing Research Tunnel

An ice shape database has been created to document ice accretions on a 21-inch chord NACA0012 model and a 72-inch chord NACA 23012 airfoil model resulting from an exposure to a Supercooled Large Drop (SLD) icing cloud with a bimodal drop size distribution. The ice shapes created were documented with photographs, laser scanned surface measurements over a section of the model span, and measurement of the ice mass over the same section of each accretion. The icing conditions used in the test matrix were based upon previously used conditions on the same models but with an alternate approach to evaluation of drop distribution effects. Ice shapes resulting from the bimodal distribution as well as from equivalent monomodal drop size distributions were obtained and compared. Results indicate that the ice shapes resulting from the monomodal and bimodal drop size distributions had similar shapes, but the bimodal distributions had greater mass and volume measurements and icing limits that extended further back on the chord of the model.

Potapczuk, Mark G.↗

Implementation of thermal residual stresses in the analysis of fiber bridged matrix crack growth in titanium matrix composites

In this research, thermal residual stresses were incorporated in an analysis of fiber-bridged matrix cracks in unidirectional and cross-ply titanium matrix composites (TMC) containing center holes or center notches. Two TMC were investigated, namely, SCS-6/Timelal-21S laminates. Experimentally, matrix crack initiation and growth were monitored during tension-tension fatigue tests conducted at room temperature and at an elevated temperature of 200 C. Analytically, thermal residual stresses were included in a fiber bridging (FB) model. The local R-ratio and stress-intensity factor in the matrix due to thermal and mechanical loadings were calculated and used to evaluate the matrix crack growth behavior in the two materials studied. The frictional shear stress term, tau, assumed in this model was used as a curve-fitting parameter to matrix crack growth data. The scatter band in the values of tau used to fit the matrix crack growth data was significantly reduced when thermal residual stresses were included in the fiber bridging analysis. For a given material system, lay-up and temperature, a single value of tau was sufficient to analyze the crack growth data. It was revealed in this study that thermal residual stresses are an important factor overlooked in the original FB models.

Bakuckas, John G., Jr.↗

Different manifestations of S-matrix poles

Making use of the analytical properties of the S-matrix and a theorem of Mittag-Leffler, model independent non-relativistic expressions for cross sections in single channel elastic scattering, scattering phase shifts and survival probabilities of resonances are derived. Provided certain conditions are satisfied by the poles, the residues can also be estimated analytically. Considerations of the low energy behaviour of the S-matrix and cross sections reveal additional conditions on the residues of the poles appearing in the Mittag-Leffler expansions. The exact expressions for the resonant cross section and phase shift are shown to reduce to the commonly used Breit–Wigner formula plus corrections. The latter is shown to approach the exact result with the example of a meson and a baryon resonance. Finally, a comparison of the exact expressions with some realistic examples is presented. The calculations of survival probabilities in particular reveal the reason behind the non-observability of non-exponential decay.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Incorporation of Plasticity and Damage Into an Orthotropic Three-Dimensional Model with Tabulated Input Suitable for Use in Composite Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased usage in the aerospace and automotive industries. While there are several composite material models currently available within commercial transient dynamic finite element codes, several features have been identified as being lacking in the currently available material models that could substantially enhance the predictive capability of the impact simulations. A specific desired feature pertains to the incorporation of both plasticity and damage within the material model. Another desired feature relates to using experimentally based tabulated stress-strain input to define the evolution of plasticity and damage as opposed to specifying discrete input properties (such as modulus and strength) and employing analytical functions to track the response of the material. To begin to address these needs, a combined plasticity and damage model suitable for use with both solid and shell elements is being developed for implementation within the commercial code LS-DYNA. The plasticity model is based on extending the Tsai-Wu composite failure model into a strain-hardening based orthotropic plasticity model with a non-associative flow rule. The evolution of the yield surface is determined based on tabulated stress-strain curves in the various normal and shear directions and is tracked using the effective plastic strain. The effective plastic strain is computed by using the non-associative flow rule in combination with appropriate numerical methods. To compute the evolution of damage, a strain equivalent semi-coupled formulation is used, in which a load in one direction results in a stiffness reduction in multiple coordinate directions. A specific laminated composite is examined to demonstrate the process of characterizing and analyzing the response of a composite using the developed model.

Impact↗

Modeled Microgravity Disrupts Collagen I/Integrin Signaling During Osteoblastic Differentiation of Human Mesenchymal Stem Cells

Spaceflight leads to reduced bone mineral density in weight bearing bones that is primarily attributed to a reduction in bone formation. We have previously demonstrated severely reduced osteoblastogenesis of human mesenchymal stem cells (hMSC) following seven days culture in modeled microgravity. One potential mechanism for reduced osteoblastic differentiation is disruption of type I collagen-integrin interactions and reduced integrin signaling. Integrins are heterodimeric transmembrane receptors that bind extracellular matrix proteins and produce signals essential for proper cellular function, survival, and differentiation. Therefore, we investigated the effects of modeled microgravity on integrin expression and function in hMSC. We demonstrate that seven days of culture in modeled microgravity leads to reduced expression of the extracellular matrix protein, type I collagen (Col I). Conversely, modeled microgravity consistently increases Col I-specific alpha2 and beta1 integrin protein expression. Despite this increase in integrin sub-unit expression, autophosphorylation of adhesion-dependent kinases, focal adhesion kinase (FAK) and proline-rich tyrosine kinase 2 (PYK2), is significantly reduced. Activation of Akt is unaffected by the reduction in FAK activation. However, reduced downstream signaling via the Ras-MAPK pathway is evidenced by a reduction in Ras and ERK activation. Taken together, our findings indicate that modeled microgravity decreases integrin/MAPK signaling, which likely contributes to the observed reduction in osteoblastogenesis.

Meyers, Valerie E.↗

Attention-based quantum tomography

Abstract With rapid progress across platforms for quantum systems, the problem of many-body quantum state reconstruction for noisy quantum states becomes an important challenge. There has been a growing interest in approaching the problem of quantum state reconstruction using generative neural network models. Here we propose the ‘attention-based quantum tomography’ (AQT), a quantum state reconstruction using an attention mechanism-based generative network that learns the mixed state density matrix of a noisy quantum state. AQT is based on the model proposed in ‘Attention is all you need’ by Vaswani et al (2017 NIPS ) that is designed to learn long-range correlations in natural language sentences and thereby outperform previous natural language processing (NLP) models. We demonstrate not only that AQT outperforms earlier neural-network-based quantum state reconstruction on identical tasks but that AQT can accurately reconstruct the density matrix associated with a noisy quantum state experimentally realized in an IBMQ quantum computer. We speculate the success of the AQT stems from its ability to model quantum entanglement across the entire quantum system much as the attention model for NLP captures the correlations among words in a sentence.

97 MATHEMATICS AND COMPUTING↗

Prediction of payload vibration environments by mechanical admittance test techniques

A series of experiments was conducted with simple beam and mass launch vehicle and payload models in order to determine the validity of mechanical admittance/impedance techniques applied to development of improved payload vibration tests. Admittances and impedances were measured from tests of the individual components to form matrices which were combined analytically to allow prediction of responses for the complete system. Results were computed for a transmission matrix approach and an admittance matrix approach. Both a rigid body and a flexible payload model were considered. The results clearly demonstrate that the transmission matrix method is too sensitive to measurement error to be practical for this application, while the pure admittance matrix method produces quite satisfactory results. The effects of various errors on the final results are demonstrated.

Kana, D. D.↗

Universal Kardar-Parisi-Zhang Dynamics in Integrable Quantum Systems

Although the Bethe ansatz solution of the spin-1/2 Heisenberg model dates back nearly a century, the anomalous nature of its high-temperature transport dynamics has only recently been uncovered. Indeed, numerical and experimental observations have demonstrated that spin transport in this paradigmatic model falls into the Kardar-Parisi-Zhang (KPZ) universality class. This has inspired the significantly stronger conjecture that KPZ dynamics, in fact, occur in all integrable spin chains with non-Abelian symmetry. Here, in this letter, we provide extensive numerical evidence affirming this conjecture. Moreover, we observe that KPZ transport is even more generic, arising in both supersymmetric and periodically driven models. Motivated by recent advances in the realization of SU(N)-symmetric spin models in alkaline-earth-based optical lattice experiments, we propose and analyze a protocol to directly investigate the KPZ scaling function in such systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗