Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix multiplication”

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 253 records · Page 14

An Efficient Spectral Method for Ordinary Differential Equations with Rational Function Coefficients

We present some relations that allow the efficient approximate inversion of linear differential operators with rational function coefficients. We employ expansions in terms of a large class of orthogonal polynomial families, including all the classical orthogonal polynomials. These families obey a simple three-term recurrence relation for differentiation, which implies that on an appropriately restricted domain the differentiation operator has a unique banded inverse. The inverse is an integration operator for the family, and it is simply the tridiagonal coefficient matrix for the recurrence. Since in these families convolution operators (i.e. matrix representations of multiplication by a function) are banded for polynomials, we are able to obtain a banded representation for linear differential operators with rational coefficients. This leads to a method of solution of initial or boundary value problems that, besides having an operation count that scales linearly with the order of truncation N, is computationally well conditioned. Among the applications considered is the use of rational maps for the resolution of sharp interior layers.

Coutsias, Evangelos A.↗

Meshless Local Petrov-Galerkin Euler-Bernoulli Beam Problems: A Radial Basis Function Approach

A radial basis function implementation of the meshless local Petrov-Galerkin (MLPG) method is presented to study Euler-Bernoulli beam problems. Radial basis functions, rather than generalized moving least squares (GMLS) interpolations, are used to develop the trial functions. This choice yields a computationally simpler method as fewer matrix inversions and multiplications are required than when GMLS interpolations are used. Test functions are chosen as simple weight functions as in the conventional MLPG method. Compactly and noncompactly supported radial basis functions are considered. The non-compactly supported cubic radial basis function is found to perform very well. Results obtained from the radial basis MLPG method are comparable to those obtained using the conventional MLPG method for mixed boundary value problems and problems with discontinuous loading conditions.

Raju, I. S.↗

A Meshless Method Using Radial Basis Functions for Beam Bending Problems

A meshless local Petrov-Galerkin (MLPG) method that uses radial basis functions (RBFs) as trial functions in the study of Euler-Bernoulli beam problems is presented. RBFs, rather than generalized moving least squares (GMLS) interpolations, are used to develop the trial functions. This choice yields a computationally simpler method as fewer matrix inversions and multiplications are required than when GMLS interpolations are used. Test functions are chosen as simple weight functions as they are in the conventional MLPG method. Compactly and noncompactly supported RBFs are considered. Noncompactly supported cubic RBFs are found to be preferable. Patch tests, mixed boundary value problems, and problems with complex loading conditions are considered. Results obtained from the radial basis MLPG method are either of comparable or better accuracy than those obtained when using the conventional MLPG method.

Raju, I. S.↗

Electromagnetic energy coupling mechanism with matrix architecture control

The present invention relates generally to reconfigurable, solid-state matrix arrays comprising multiple rows and columns of reconfigurable secondary mechanisms that are independently tuned. Specifically, the invention relates to reconfigurable devices comprising multiple, solid-state mechanisms characterized by at least one voltage-varied parameter disposed within a flexible, multi-laminate film, which are suitable for use as magnetic conductors, ground surfaces, antennas, varactors, ferrotunable substrates, or other active or passive electronic mechanisms.

Knowles, Gareth↗

Evaluating Secondary Inorganic Aerosols in Three Dimensions

The spatial distribution of aerosols and their chemical composition dictates whether aerosols have a cooling or a warming effect on the climate system. Hence, properly modeling the three-dimensional distribution of aerosols is a crucial step for coherent climate simulations. Since surface measurement networks only give 2-D data, and most satellites supply integrated column information, it is thus important to integrate aircraft measurements in climate model evaluations. In this study, the vertical distribution of secondary inorganic aerosol (i.e., sulfate, ammonium, and nitrate) is evaluated against a collection of 14 AMS flight campaigns and surface measurements from 2000 to 2010 in the USA and Europe. GISS ModelE2 is used with multiple aerosol microphysics (MATRIX, OMA) and thermodynamic (ISORROPIA II, EQSAM) configurations. Our results show that the MATRIX microphysical scheme improves the model performance for sulfate, but that there is a systematic underestimation of ammonium and nitrate over the USA and Europe in all model configurations. In terms of gaseous precursors, nitric acid concentrations are largely underestimated at the surface while overestimated in the higher levels of the model. Heterogeneous reactions on dust surfaces are an important sink for nitric acid, even high in the troposphere. At high altitudes, nitrate formation is calculated to be ammonia limited. The underestimation of ammonium and nitrate in polluted regions is most likely caused by a too simplified treatment of the NH3/NH4(+) partitioning which affects the HNO3/NO3(-) partitioning.

aerosols↗

A Bi-Modal Distribution of ALHA77307 Matrix Olivine: Evidence for Fine-Grained Mixing from Multiple Reservoirs in the CO Formation Zone

The CO 3.0 chondrite ALHA77307 is thought to be the least metamorphosed of all the CO chondrites [1]. As such, the fine-grained (<30 m) olivine found in its matrix is a valuable resource for investigating the CO formation environment since its compositions should be primary. In the CO matrix, we indeed find a wide range of major element compositions (Fa(0.5-71)). However, more importantly, we find that the olivines make up two compositionally distinct populations (Fa(0.5-5) and Fa(21-71)). Grains from both populations are found within an extremely close proximity and we see no obvious evidence of two distinct lithologies within our samples. Therefore, we conclude that the olivine grains found in the ALHA77307 matrix must have crystallized within two unique formation conditions and were later mixed at a very fine scale during the accretion epoch. Here, we propose a possible explanation based on Cr and Mn concentrations in the olivine.

Frank, D.↗

Development of a microwave 20 x 20 switch matrix for 30/20 GHz SS-TDMA application

The design and fabrication of a 3-8 GHz, 20 x 20 Satellite Switched-Time Division Multiple Access IF switch matrix applicable to a 30/20 GHz communications satellite are described. An assessment of switch architecture in 1980 concluded that the GaAs FET-based coupled crossbar switch matrix, incorporating high speed CMOS LSI logic for switch crosspoint addressing, would be the optimum technology available for communications satellite switching by 1982. This assessment was based on such factors as switching speed, bandwidth, off-state isolation, and reliability, over a 10-year mission life. A proof-of-concept model's construction and testing are presented.

Cory, B. J.↗

Response of a Rotating Propeller to Aerodynamic Excitation

The flexural vibration of a rotating propeller blade with clamped shank is analyzed with the object of presenting, in matrix form, equations for the elastic bending moments in forced vibration resulting from aerodynamic forces applied at a fixed multiple of rotational speed. Matrix equations are also derived which define the critical speeds end mode shapes for any excitation order and the relation between critical speed and blade angle. Reference is given to standard works on the numerical solution of matrix equations of the forms derived. The use of a segmented blade as an approximation to a continuous blade provides a simple means for obtaining the matrix solution from the integral equation of equilibrium, so that, in the numerical application of the method presented, the several matrix arrays of the basic physical characteristics of the propeller blade are of simple form, end their simplicity is preserved until, with the solution in sight, numerical manipulations well-known in matrix algebra yield the desired critical speeds and mode shapes frame which the vibration at any operating condition may be synthesized. A close correspondence between the familiar Stodola method and the matrix method is pointed out, indicating that any features of novelty are characteristic not of the analytical procedure but only of the abbreviation, condensation, and efficient organization of the numerical procedure made possible by the use of classical matrix theory.

Arnoldi, Walter E.↗

Parallel Embeddings: a Visualization Technique for Contrasting Learned Representations

We introduce ``Parallel Embeddings'', a new technique that generalizes the classical Parallel Coordinates visualization technique from a sequence of features within a single matrix to a sequence of multiple feature matrices. This is accomplished by separately embedding each matrix in the sequence into a 2-D frame, juxtaposing these 2-D representations, then connecting points between frames. We outline how this visualization can be used for several ``model comparison'' tasks in machine learning. We compare user performance with Parallel Embeddings to a baseline using the TensorFlow Embedding Projector for estimating model accuracy and understanding conceptual model differences. We found that users were more accurate with Parallel Embeddings for model comparison and that users learned how to use Parallel Embeddings more quickly than TensorFlow Embedding Projector. Furthermore, users' analytical process using Embedding Projector was positively affected by using Parallel Embeddings beforehand.

Model comparison, machine learing, machine learnin↗

Moments of parton distribution functions of any order from lattice QCD

We describe a procedure to determine moments of parton distribution functions of any order in lattice quantum chromodynamics (QCD). The procedure is based on the gradient flow for fermion and gauge fields. The flowed matrix elements of twist-2 operators renormalize multiplicatively, and the matching with the physical matrix elements can be obtained using continuum symmetries and the irreducible representations of Euclidean 4-dimensional rotations. We calculate the matching coefficients at one-loop in perturbation theory for moments of any order in the flavor nonsinglet case. We also give specific examples of operators that could be used in lattice QCD computations. It turns out that it is possible to choose operators with identical Lorentz indices and still have a multiplicative matching. One can thus use twist-2 operators exclusively with temporal indices, thus substantially improving the signal-to-noise ratio in the computation of the hadronic matrix elements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Glauber-theory analysis of nuclear reactions on a 12 C target with variational Monte Carlo wave functions

The application of Glauber theory has been playing an increasingly important role with the study of unstable or exotic nuclei. Its adaptation to medium and high-energy nucleus-nucleus collisions is severely limited because one has to evaluate the matrix elements of multiple-scattering operators. The extraction of physical observables has been done using ‘approximate’ Glauber theory whose validity is hard to evaluate. Here, we perform a full calculation of the matrix elements using Monte Carlo integration and analyze the elastic differential cross sections and the total reaction cross sections for p+¹²C, ⁴,⁶He+¹²C, and ¹²C+¹²C collisions. We use the variational Monte Carlo wave functions for ⁴,⁶He and ¹²C obtained by using realistic two- and three-nucleon potentials. We demonstrate the performance of the Glauber-theory calculations by comparing with available experimental data. We further discuss the accuracy of the conventional approximate methods in the light of the cumulant expansion for Glauber’s phase-shift function.

Horiuchi, W. [Osaka Metropolitan University (Japan↗

Elastic and Piezoelectric Properties of Boron Nitride Nanotube Composites: Finite Element Model - Part II

This paper is the second part of a two-part series where the first part presents a molecular dynamics model of a single Boron Nitride Nanotube (BNNT) and this paper scales up to multiple BNNTs in a polymer matrix. This paper presents finite element (FE) models to investigate the effective elastic and piezoelectric properties of (BNNT) nanocomposites. The nanocomposites studied in this paper are thin films of polymer matrix with aligned co-planar BNNTs. The FE modelling approach provides a computationally efficient way to gain an understanding of the material properties. We examine several FE models to identify the most suitable models and investigate the effective properties with respect to the BNNT volume fraction and the number of nanotube walls. The FE models are constructed to represent aligned and randomly distributed BNNTs in a matrix of resin using 2D and 3D hollow and 3D filled cylinders. The homogenisation approach is employed to determine the overall elastic and piezoelectric constants for a range of volume fractions. These models are compared with an analytical model based on Mori-Tanaka formulation suitable for finite length cylindrical inclusions. The model applies to primarily single-wall BNNTs but is also extended to multi-wall BNNTs, for which preliminary results will be presented. Results from the Part 1 of this series can help to establish a constitutive relationship for input into the finite element model to enable the modeling of multiple BNNTs in a polymer matrix.

Kim, H. Alicia↗

Properties and microstructure evolution of silicon nitride and zirconium nitride following Ni ion irradiation

We report that silicon nitride and zirconium nitride have been proposed as potential materials for multiple nuclear applications (inert matrix fuels, accident tolerant fuels, space nuclear power, fusion reactor design), yet knowledge on their behavior under irradiation remains limited. Ion irradiations were performed using 15 MeV Ni 5+ ions on Si 3 N 4 and ZrN samples, with midrange doses (around 3 µm) from 1 to 50 dpa and temperatures from 300 to 700°C. Volumetric lattice swelling was determined by grazing incidence X-ray diffraction, defect production and evolution were tracked using Transmission Electron Microscopy, and nanoindentation was performed to quantify the ceramics’ mechanical properties evolution. The results from these irradiation studies on nitride ceramics help fill the current gap present in the literature. Behavior consistent with past work on irradiated Si 3 N 4 was observed with respect to mechanical properties and defect formation up to 15 dpa and 500°C. Failure of the grain boundary sintering aid in Si 3 N 4 was observed above these conditions. Different behavior was observed in both nitrides at 50 dpa and 700°C, where lattice swelling increased past potential saturation values. Unreported cavity formation was witnessed in both materials under all irradiation conditions, with stable number density and slight size increase above 15 dpa. The mechanism for the cavity formation remains to be determined.

36 MATERIALS SCIENCE↗

Block Krylov Subspace Methods for Functions of Matrices II: Modified Block FOM

We analyze an expansion of the generalized block Krylov subspace framework of [Electron. Trans. Numer. Anal., 47 (2017), pp. 100--126]. This expansion allows the use of low-rank modifications of the matrix projected onto the block Krylov subspace and contains, as special cases, the block GMRES method and the new block Radau--Arnoldi method. Within this general setting, we present results that extend the interpolation property from the nonblock case to a matrix polynomial interpolation property for the block case, and we relate the eigenvalues of the projected matrix to the latent roots of these matrix polynomials. Some error bounds for these modified block FOM methods for solving linear systems are presented. We then show how cospatial residuals can be preserved in the case of families of shifted linear block systems. This result is used to derive computationally practical restarted algorithms for block Krylov approximations that compute the action of a matrix function on a set of several vectors simultaneously. Finally, we prove some error bounds and present numerical results showing that two modifications of FOM, the block harmonic and the block Radau--Arnoldi methods for matrix functions, can significantly improve the convergence behavior.

97 MATHEMATICS AND COMPUTING↗

Eigenvalues of the Laplacian of a graph

Let G be a finite undirected graph with no loops or multiple edges. The Laplacian matrix of G, Delta(G), is defined by Delta sub ii = degree of vertex i and Delta sub ij = -1 if there is an edge between vertex i and vertex j. The structure of the graph G is related to the eigenvalues of Delta(G); in particular, it is proved that all the eigenvalues of Delta(G) are nonnegative, less than or equal to the number of vertices, and less than or equal to twice the maximum vertex degree. Precise conditions for equality are given.

Anderson, W. N., Jr.↗

Improved impact-resistant boron-aluminum composites for use as turbine engine fan blades

Efforts to improve the impact resistance of B/Al are reviewed and analyzed. Thin sheet Charpy and Izod impact tests and standard full size Charpy impact tests were conducted on unidirectional and angleply composites containing 4, 5.6 and 8 mil boron in 1100, 2024, 5052 and 6061 Al matrices. Impact failure modes of B/Al are proposed in an attempt to describe the mechanisms involved and to provide insight for maximizing impact resistance. The impact strength of B/Al was significantly increased by proper selection of materials and processing. The use of more ductile matrices (1100 Al) and larger diameter (8 mil) boron fibers gave the highest impact strengths by allowing matrix shear deformation and multiple fiber breakage.

Mcdanels, D. L.↗

Improved impact-resistant boron/aluminum composites for use as turbine engine fan blades

Thin-sheet Charpy and Izod impact tests and standard full-size Charpy impact tests were conducted on unidirectional and angleply composites containing 4, 5.6 and 8 mil boron in 1100, 2024, 5052 and 6061 Al matrices. Impact failure modes of B/Al are proposed to describe the mechanisms involved and to provide insight for maximizing impact resistance. The impact strength of B/Al was significantly increased by proper selection of materials and processing. The use of more ductile matrices and larger diameter boron fibers gave the highest impact strengths by allowing matrix shear deformation and multiple fiber breakage. Pendulum impact test results of improved B/Al were higher than notched titanium and indicate sufficient foreign object damage protection to warrant consideration of B/Al for application to fan blades in aircraft gas turbine engines.

Mcdanels, D. L.↗

Effect of fiber diameter and matrix alloys on impact-resistant boron/aluminum composites

Efforts to improve the impact resistance of B/Al are reviewed and analyzed. Nonstandard thin-sheet charpy and Izod impact tests and standard full-size Charpy impact tests were conducted on composites containing unidirectional 0.10mm, 0.14mm, and 0.20mm diameter boron fibers in 1100, 2024, 5052, and 6061 Al matrices. Impact failure modes of B/Al are proposed in an attempt to describe the mechanisms involved and to provide insight for maximizing impact resistance. The impact strength of B/Al was significantly increased by proper selection of materials and processing. The use of a ductile matrix and large diameter boron fibers gave the highest impact strengths. This combination resulted in improved energy absorption through matrix shear deformation and multiple fiber breakage.

Mcdanels, D. L.↗