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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Propulsion Electrification Architecture Selection Process and Cost of Carbon Abatement Analysis for Heavy-Duty Off-Road Material Handler

The heavy-duty off-road industry continues to expand efforts to reduce fuel consumption and CO 2 e (carbon dioxide equivalent) emissions. Many manufacturers are pursuing electrification to decrease fuel consumption and emissions. Future policies will likely require electrification for CO2e savings, as seen in light-duty on-road vehicles. Electrified architectures vary widely in the heavy-duty off-road space, with parallel hybrids in some applications and series hybrids in others. The diverse applications for different types of equipment mean different electrified configurations are required. Companies must also determine the value in pursuing electrified architectures; this work analyzes a range of electrified architectures, from micro hybrids to parallel hybrids to series hybrids to a BEV, looking at the total cost, total CO 2 e, and cost per CO 2 e (cost of carbon abatement, or cost of carbon reduction) using data for the year 2021. This study is focused on a heavy-duty off-road material handler, the Pettibone Cary-Lift 204i. This machine’s specialty application, including events like unloading large oil pipes from a railcar, requires a unique electrified architecture that suits its specific needs. However, the results from this study may be extrapolated to similar machinery to inform fuel savings options across the heavy-duty off-road industry. In this study, a unique electrified architecture is determined for the Cary-Lift. This architecture is informed by multiple rounds of a Pugh matrix decision analysis to select a shortened list of desirable electrified architectures. The shortened list is modeled and simulated to determine CO 2 e, cost, and cost per CO 2 e. A final architecture is determined as a plug-in series hybrid that reduces fuel consumption by 65%, targeting the large fuel and CO 2 e savings that are likely to be required for the future of the heavy-duty off-road industry.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1987

A simple method to model truss-beams as equivalent continua

A new method is presented for obtaining stiffness constants for structures modeled as three-dimensional linear elastic beams. The three-dimensional lattice is comprised of multiple two-dimensional lattices. Matrix manipulation is used to determine the contribution of each two-dimensional substructure to the complete three dimensional lattice.

Barton, O., Jr.↗
NASA NTRS1991

Burst-Compression And -Expansion For TDMA Communication

Burst-compression and -expansion technique enables interconnection of users transmitting and receiving data at rates asynchronous with respect to clocks within ground terminals of satellite-switched, time-division-multiple-access (TDMA) communication network. Matrix switch aboard satellite routes bursts of data from source users received on uplink antennas to downlink antennas illuminating ground areas containing destination users. TDMA ground terminal compresses streams of data from source users into rapid bursts for transmission and reexpands bursts of received data into slower streams of data for delivery to destination users. Greater flexibility in interconnecting widely dispersed users achieved by use of hopping beams.

Budinger, James M.↗
NASA NTRS2010

InGaP Heterojunction Barrier Solar Cells

A new solar-cell structure utilizes a single, ultra-wide well of either gallium arsenide (GaAs) or indium-gallium-phosphide (InGaP) in the depletion region of a wide bandgap matrix, instead of the usual multiple quantum well layers. These InGaP barrier layers are effective at reducing diode dark current, and photogenerated carrier escape is maximized by the proper design of the electric field and barrier profile. With the new material, open-circuit voltage enhancements of 40 and 100 mV (versus PIN control systems) are possible without any degradation in short-circuit current. Basic tenets of quantum-well and quantum- dot solar cells are utilized, but instead of using multiple thin layers, a single wide well works better. InGaP is used as a barrier material, which increases open current, while simultaneously lowering dark current, reducing both hole diffusion from the base, and space charge recombination within the depletion region. Both the built-in field and the barrier profile are tailored to enhance thermionic emissions, which maximizes the photocurrent at forward bias, with a demonstrated voltage increase. An InGaP heterojunction barrier solar cell consists of a single, ultra-wide GaAs, aluminum-gallium-arsenide (AlGaAs), or lower-energy-gap InGaP absorber well placed within the depletion region of an otherwise wide bandgap PIN diode. Photogenerated electron collection is unencumbered in this structure. InGaAs wells can be added to the thick GaAs absorber layer to capture lower-energy photons.

Welser, Roger E.↗
NASA NTRS2018

Results of the Aerosol Sampling Experiment on the International Space Station

The indoor air quality of the International Space Station (ISS), in terms of particulate matter, has been assessed through a sampling experiment that collected airborne debris for return to Earth. The Aerosol Sampling Experiment took place on ISS in December 2016 and January 2017, with sample return in February 2017. For redundancy, two different aerosol samplers collected particles in two size ranges spanning from 10 nm to hundreds of micrometers. Multiple analysis techniques have been used to characterize the particles, such as reflective light microscopy, manual scanning electron microscopy (SEM) analysis (secondary and backscattered imaging) with energy dispersive spectroscopy (EDS) x-ray analysis, computer controlled scanning electron microscopy (CCSEM) and Raman spectroscopy. Sampling took place in five locations on ISS, with the greatest deposition of particles in Node 3, followed by Node 1, the US Lab and Node 2. The PMM had very few airborne particles, as expected, owing to the lack of human presence. A large proportion of the sampled material consisted of fibers from clothing. Many particles had complex morphologies or contained multiple metals in a carbonaceous matrix. Metal particles were analyzed for each sampling location in terms of both weight and number percent, with 27 heavy metals detected. Fiberglass and TiO2 particles were present in every location. Quantitative analysis results are presented, and where possible, particle emission sources are identified with the intention to fuel discussions on mitigation.

Meyer, Marit E.↗
NASA NTRS2022

Full-Wave Simulations of Scattering by Corn Fields at L-Band

In this paper, the Numerical Maxwell Model of 3D (NMM3D) full-wave simulation is performed over a corn field using a hybrid method to study the vegetation effect on the microwave. The commercial software of FEKO is used to extract T-matrix of single corn in the first step. Then the calculated T-matrix is combined with Wave Multiple Scattering Theory (W-MST) in the second step to consider the multiple scattering among different plants. The hybrid method is validated with HFSS by solving scattering from 2 corns. A corn field of 25 corn is simulated using the hybrid method and the transmission is calculated and compared with those obtained from the classical radiative transfer model.

Yueh, Simon↗
NASA NTRS2003

Continuous-Time Bilinear System Identification

The objective of this paper is to describe a new method for identification of a continuous-time multi-input and multi-output bilinear system. The approach is to make judicious use of the linear-model properties of the bilinear system when subjected to a constant input. Two steps are required in the identification process. The first step is to use a set of pulse responses resulting from a constant input of one sample period to identify the state matrix, the output matrix, and the direct transmission matrix. The second step is to use another set of pulse responses with the same constant input over multiple sample periods to identify the input matrix and the coefficient matrices associated with the coupling terms between the state and the inputs. Numerical examples are given to illustrate the concept and the computational algorithm for the identification method.

Juang, Jer-Nan↗
NASA NTRS1994

Linear and nonlinear dynamic analysis of redundant load path bearingless rotor systems

The goal of this research is to develop the transfer matrix method to treat nonlinear autonomous boundary value problems with multiple branches. The application is the complete nonlinear aeroelastic analysis of multiple-branched rotor blades. Once the development is complete, it can be incorporated into the existing transfer matrix analyses. There are several difficulties to be overcome in reaching this objective. The conventional transfer matrix method is limited in that it is applicable only to linear branch chain-like structures, but consideration of multiple branch modeling is important for bearingless rotors. Also, hingeless and bearingless rotor blade dynamic characteristics (particularly their aeroelasticity problems) are inherently nonlinear. The nonlinear equations of motion and the multiple-branched boundary value problem are treated together using a direct transfer matrix method. First, the formulation is applied to a nonlinear single-branch blade to validate the nonlinear portion of the formulation. The nonlinear system of equations is iteratively solved using a form of Newton-Raphson iteration scheme developed for differential equations of continuous systems. The formulation is then applied to determine the nonlinear steady state trim and aeroelastic stability of a rotor blade in hover with two branches at the root. A comprehensive computer program is developed and is used to obtain numerical results for the (1) free vibration, (2) nonlinearly deformed steady state, (3) free vibration about the nonlinearly deformed steady state, and (4) aeroelastic stability tasks. The numerical results obtained by the present method agree with results from other methods.

Murthy, V. R.↗
DOE OSTI2020

Source identification by non-negative matrix factorization combined with semi-supervised clustering

Machine-learning methods and apparatus are provided to solve blind source separation problems with an unknown number of sources and having a signal propagation model with features such as wave-like propagation, medium-dependent velocity, attenuation, diffusion, and/or advection, between sources and sensors. In exemplary embodiments, multiple trials of non-negative matrix factorization are performed for a fixed number of sources, with selection criteria applied to determine successful trials. A semi-supervised clustering procedure is applied to trial results, and the clustering results are evaluated for robustness using measures for reconstruction quality and cluster separation. The number of sources is determined by comparing these measures for different trial numbers of sources. Source locations and parameters of the signal propagation model can also be determined. Disclosed methods are applicable to a wide range of spatial problems including chemical dispersal, pressure transients, and electromagnetic signals, and also to non-spatial problems such as cancer mutation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Source identification by non-negative matrix factorization combined with semi-supervised clustering

Machine-learning methods and apparatus are provided to solve blind source separation problems with an unknown number of sources and having a signal propagation model with features such as wave-like propagation, medium-dependent velocity, attenuation, diffusion, and/or advection, between sources and sensors. In exemplary embodiments, multiple trials of non-negative matrix factorization are performed for a fixed number of sources, with selection criteria applied to determine successful trials. A semi-supervised clustering procedure is applied to trial results, and the clustering results are evaluated for robustness using measures for reconstruction quality and cluster separation. The number of sources is determined by comparing these measures for different trial numbers of sources. Source locations and parameters of the signal propagation model can also be determined. Disclosed methods are applicable to a wide range of spatial problems including chemical dispersal, pressure transients, and electromagnetic signals, and also to non-spatial problems such as cancer mutation.

Alexandrov, Boian S.↗
NASA NTRS2001

On-Orbit Calibration of Redundant Spacecraft Gyros by Optimal Reduction to Three Axes

The Aqua spacecraft will carry four single-axis gyros configured with three orthogonal axes and one skew axis. This redundancy presents a challenge for batch methods of on-orbit gyro calibration that use a spacecraft rotation model deterministically related to gyro data, in that sensor data can respond to at most three angular velocity components. When the number of gyros, N, is greater than 3, the 3xN matrix, G, that reduces the N gyro measurements to three body-frame angular-velocity components cannot be fully determined by such methods; there are many such matrices that produce essentially the same angular velocity history. In such a case, spacecraft operators require information about the Nx3 gyro linear response matrix, R, that relates gyro outputs to the body-frame angular velocities causing them. This matrix provides sufficient information to determine multiple reduced-dimension G-matrices for use in case of failure or degradation of one or more gyros, as well as to determine an optimal 3xN G for the fully-functional configuration. A viable proposal is to apply a 3xN pre-filter matrix, F, to the N gyro outputs before carrying out a conventional gyro calibration procedure. The angular-velocity history emerging from conventional calibration may then be used as input data, together with the same gyro data that generated it, to fit the alignment, scale-factor, and bias parameters of each gyro axis in turn. A difficulty of such a proposal is the arbitrariness in the choice of F. Due to gyro noise, different pre-filter matrices produce different calibrations. This paper presents a method of choosing F that is based on optimizing gyro consistency in the limit of infinite weight on gyro data, as compared to sensor data. The choice of F is independent of a priori alignment and is based on the gyro data alone. The method is applicable to any N of three or more, but reduces to conventional batch-estimation methodologies when N = 3. Results of computational comparison among calibration simulations using various choices of F will be presented for the Aqua gyro configuration with N = 4.

Radomski, M. S.↗
DOE OSTI2026

An Adaptive Multiparameter Penalty Selection Method for Multiconstraint and Multiblock ADMM

This work presents a new method for online selection of multiple penalty parameters for the alternating direction method of multipliers (ADMM) algorithm applied to optimization problems with multiple constraints or functions with block matrix components. ADMM is widely used for solving constrained optimization problems in a variety of fields, including signal and image processing. Implementations of ADMM often utilize a single hyperparameter, referred to as the penalty parameter, which needs to be tuned to control the rate of convergence. However, in problems with multiple constraints, ADMM may demonstrate slow convergence regardless of penalty parameter selection due to scale differences between constraints. Accounting for scale differences between constraints to improve convergence in these cases requires introducing a penalty parameter for each constraint. The proposed method is able to adaptively account for differences in scale between constraints, providing robustness with respect to problem transformations and initial selection of penalty parameters. It is also simple to understand and implement. Our numerical experiments demonstrate that the proposed method performs favorably compared to a variety of existing penalty parameter selection methods.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

Characterization of unnotched SCS-6/Ti-15-3 metal matrix composites at 650 C

Ti-15-3 reinforced with SCS-6 silicon carbide fibers, in five different layups, was tested at 650 C to determine monotonic and fatigue strengths, basic mechanical properties, and damage initiation and progression. The elevated temperature results were compared to those obtained at room temperature. Analytical predictions were made of the monotonic stress-strain response as well as cyclic stress-strain hysteresis. The fiber reinforcement was found to significantly increase the static and fatigue strengths of the laminates over that of the matrix material at elevated temperature while the increase was insignificant at room temperature. Initial damage, in either the fibers or the matrix, was partitioned as a function of the life and applied strain range in the constituents. High strains and short lives resulted in multiple fiber failure with no signs of matrix fatigue cracking. Low strains and long lives resulted in extensive matrix cracking and no fiber breaks away from the fracture surface. At 650 C the matrix was too weak to cause fiber-matrix interface failure prior to matrix yielding. Laminate fatigue lives were hypothesized to be a function of the 0 deg fiber stress. More scatter was found in the 0 deg fiber stress vs. high temperature fatigue life data than for the room temperature data. An initial unloading modulus that was greater than the initial loading modulus was observed in the elevated temperature fatigue tests.

Pollock, W. D.↗
NASA NTRS1992

Characterization of unnotched SCS-6/Ti-15-3 metal matrix composites at 650 C

Ti-15-3 reinforced with SCS-6 silicon carbide fibers, in five different layups, was tested at 650 C to determine monotonic and fatigue strengths, basic mechanical properties, and damage initiation and progression. The elevated temperature results were compared to those obtained at room temperature. Analytical predictions were made of the monotonic stress-strain response as well as cyclic stress-strain hysteresis. The fiber reinforcement was found to significantly increase the static and fatigue strengths of the laminates over that of the matrix material at elevated temperature while the increase was insignificant at room temperature. Initial damage, in either the fibers or the matrix, was partitioned as a function of the life and applied strain range in the constituents. High strains and short lives resulted in multiple fiber failure with no signs of matrix fatigue cracking. Low strains and long lives resulted in extensive matrix cracking and no fiber breaks away from the fracture surface. At 650 C the matrix was too weak to cause fiber-matrix interface failure prior to matrix yielding. Laminate fatigue lives were hypothesized to be a function of the 0 deg fiber stress. More scatter was found in the 0 deg fiber stress vs. high temperature fatigue life data than for the room temperature data. An initial unloading modulus that was greater than the initial loading modulus was observed in the elevated temperature fatigue tests.

Pollock, William D.↗
DOE OSTI2022

Introduction of the Adding and Doubling Method for Solving Bateman Equations for Nuclear Fuel Depletion

This paper introduces and evaluates the Adding and Doubling Method (ADM) for solving the Bateman equations for depletion systems with varying numbers of nuclides and compares it to the Chebyshev Rational Approximation Method (CRAM), both implemented in the reactor physics analysis application Griffin. ADM, when applied to the Crank-Nicolson Finite Difference method, can produce results comparable in accuracy and precision to CRAM with comparable run times for systems with 35 or 297 nuclides. For systems with more than 300 nuclides, the matrix-matrix operations required by ADM are significantly more costly than the matrix-vector operations required by CRAM, making CRAM the more efficient method for systems with large numbers of nuclides. ADM is an accurate method that maintains other advantages over CRAM in that it does not depend on pre-generated coefficients or require complex number operations. ADM also manages to outperform CRAM by a factor of more than 250 in terms of run time for depletion systems that require multiple Bateman solves while the depletion matrix and time step size remain constant over all depletion intervals.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Interface instabilities in hafnium hydride entrained iron metal matrix composites

The chemical interactions in Fe–HfH 2 metal matrix composites (MMCs) are studied across multiple length scales to elucidate the decomposition of the parent phases and corresponding reaction zone physics during direct current sintering. Fe–HfH 2 composites were synthesized with increasing as-mixed hydride contents of Fe–25% HfH 2 , Fe–40% HfH 2 , Fe–55% HfH 2 , and Fe–70% HfH 2 (all in vol. %) to demonstrate the ability to achieve sintered MMCs with target hydride contents. Samples were probed across multiple length scales through a multi-modal workflow employing x-ray diffraction, scanning electron microscopy and segmentation analysis, and synchrotron techniques including hard x-ray fluorescence mapping and nanoprobe x-ray absorption near-edge structure measurements. Under the selected sintering temperature and pressure conditions, hydrogen evolution is seen to evolve through parallel paths: thermal decomposition from during the transformation of HfH 2 to HfH x<2 and through subsequent reaction with the Fe matrix leading to intermetallic phase formation. Specifically, HfFe and HfFe 2 intermetallic formation accelerates the release of hydrogen with a subsequent HfO 2 phase forming at grain boundaries. For this MMC, the consumption or loss of hydrogen can be considerable in compacts with initial hydride loading of 25%–40% HfH 2 approaching 83% hydrogen loss for the lower volume fraction composites. Increasing the volume fraction of HfH 2 to 70% enhanced the retained hydrogen content to 53% and attributed to the reduced interfacial area intrinsic to the increased HfH 2 loading in this MMC.

36 MATERIALS SCIENCE↗
NASA NTRS1993

A numerical model for predicting crack path and modes of damage in unidirectional metal matrix composites

A finite element-based numerical technique has been developed to simulate damage growth in unidirectional composites. This technique incorporates elastic-plastic analysis, micromechanics analysis, failure criteria, and a node splitting and node force relaxation algorithm to create crack surfaces. Any combination of fiber and matrix properties can be used. One of the salient features of this technique is that damage growth can be simulated without pre-specifying a crack path. In addition, multiple damage mechanisms in the forms of matrix cracking, fiber breakage, fiber-matrix debonding and plastic deformation are capable of occurring simultaneously. The prevailing failure mechanism and the damage (crack) growth direction are dictated by the instantaneous near-tip stress and strain fields. Once the failure mechanism and crack direction are determined, the crack is advanced via the node splitting and node force relaxation algorithm. Simulations of the damage growth process in center-slit boron/aluminum and silicon carbide/titanium unidirectional specimens were performed. The simulation results agreed quite well with the experimental observations.

Bakuckas, J. G.↗
NASA NTRS1999

Five-Hole Flow Angle Probe Calibration for the NASA Glenn Icing Research Tunnel

A spring 1997 test section calibration program is scheduled for the NASA Glenn Research Center Icing Research Tunnel following the installation of new water injecting spray bars. A set of new five-hole flow angle pressure probes was fabricated to properly calibrate the test section for total pressure, static pressure, and flow angle. The probes have nine pressure ports: five total pressure ports on a hemispherical head and four static pressure ports located 14.7 diameters downstream of the head. The probes were calibrated in the NASA Glenn 3.5-in.-diameter free-jet calibration facility. After completing calibration data acquisition for two probes, two data prediction models were evaluated. Prediction errors from a linear discrete model proved to be no worse than those from a full third-order multiple regression model. The linear discrete model only required calibration data acquisition according to an abridged test matrix, thus saving considerable time and financial resources over the multiple regression model that required calibration data acquisition according to a more extensive test matrix. Uncertainties in calibration coefficients and predicted values of flow angle, total pressure, static pressure. Mach number. and velocity were examined. These uncertainties consider the instrumentation that will be available in the Icing Research Tunnel for future test section calibration testing.

Gonsalez, Jose C.↗