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 163 records · Page 9

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

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↗

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↗

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

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

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

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

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

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

MATRIX-VBS (v1.0): Implementing an Evolving Organic Aerosol Volatility in an Aerosol Microphysics Model

The gas-particle partitioning and chemical aging of semi-volatile organic aerosol are presented in a newly developed box model scheme, where its effect on the growth, composition, and mixing state of particles is examined. The volatility-basis set (VBS) framework is implemented into the aerosol microphysical scheme MATRIX (Multiconfiguration Aerosol TRacker of mIXing state), which resolves mass and number aerosol concentrations and in multiple mixing-state classes. The new scheme, MATRIX-VBS, has the potential to significantly advance the representation of organic aerosols in Earth system models by improving upon the conventional representation as non-volatile particulate organic matter, often also with an assumed fixed size distribution. We present results from idealized cases representing Beijing, Mexico City, a Finnish forest, and a southeastern US forest, and investigate the evolution of mass concentrations and volatility distributions for organic species across the gas and particle phases, as well as assessing their mixing state among aerosol populations. Emitted semi-volatile primary organic aerosols evaporate almost completely in the intermediate-volatility range, while they remain in the particle phase in the low-volatility range. Their volatility distribution at any point in time depends on the applied emission factors, oxidation by OH radicals, and temperature. We also compare against parallel simulations with the original scheme, which represented only the particulate and non-volatile component of the organic aerosol, examining how differently the condensed-phase organic matter is distributed across the mixing states in the model. The results demonstrate the importance of representing organic aerosol as a semi-volatile aerosol, and explicitly calculating the partitioning of organic species between the gas and particulate phases.

volatility-basis set↗

AR-40 AR-39 Age of an Impact-Melt Lithology in DHOFAR 961

The South Pole-Aitken (SPA) basin is the stratigraphically oldest identifiable lunar basin and is therefore one of the most important targets for absolute age-dating to help understand whether ancient lunar bombardment history smoothly declined or was punctuated by a cataclysm. The SPA basin also has another convenient property, a geochemically distinct interior, unobscured by extensive mare basalt fill. A case has been made for the possible origin of the Dhofar 961 lunar meteorite in the South Pole-Aitken (SPA) basin, based on comparing its composition with Lunar Prospector gamma-ray data for the interior of the SPA basin. Dhofar 961 contains several different impact-melt (IM) lithologies. Jolliff et al. described two classes of mafic impact-melt lithologies, one dominated by olivine (Lithology A) and the other by plagioclase (An 95-96.5) (Lithology B). Broad-beam analyses of these lithologies yielded (is) approximately 14.0 wt% FeO, 11.7 wt% MgO, and 15.4 wt% Al2O3. Lithologies A and B differ by approximately 2.5% Al2O3, 1.5% FeO and 1.5% MgO, consistent with the occurrence of olivine phenocrysts in A and plagioclase clasts in B. Both lithologies are considerably more mafic than the Apollo mafic impact-melt breccias, corresponding to olivine gabbronorite. Joy et al. used U-Pb dating to investigate phosphate fragments in the Dhofar 961 matrix and impact-melt clasts. Matrix phosphates have 4.34 to 4 Ga ages, consistent with ancient KREEP-driven magmatic episodes and Pre-Nectarian ((is) greater than 3.92 Ga). Phosphates found within Dhofar 961 crystalline impact melt breccia clasts range from 4.26 to 3.89 Ga, potentially recording events throughout the basin forming epoch of lunar history. The youngest reset ages in the Dhofar 961 sample represent an upper limit for the time of formation of the meteorite. Joy et al suggested this age represents the final impact that mixed and consolidated several generations of precursor rocks into the Dhofar meteorite group, although they note that further age dating of all the stones is required to test this hypothesis. We received a split of Dhofar 961 from R. Zeigler consisting of a large clast of IM Lithology B, with some light-colored, friable matrix clinging to the external margins of the impact-melt clast. This lithology was not present in the samples investigated by Joy et al. and thus does not have corresponding U-Pb ages on it. We created multiple subsplits of both the IM and matrix lithologies, each weighing several tens of micrograms. We conducted Ar-40 Ar-39 dating of this candidate SPA material by high-resolution step heating and comparing it with the regolith that surrounds it.

Lithology↗

Experimental and theoretical investigation of stress wave attenuation in fiber reinforced composites.

The propagation of an initialrcidrical pressure pulse through a linear elastic fiber reinforced composite medium is analysed, both experimentally and analytically. In the experiment, tests were performed on plates with single and multiple circular inclusions embedded in a matrix of lower characteristic impedance. Sharp compression pulses were generated at an edge of the plate. Strain gages were mounted on various positions of the plate to determine the attenuation of the transient stress in the fiber reinforced composite. The qualitative analytical treatment is based on the methods of propagating stress discontinuities. Computer programs were written to numerically determine the changes in the shape of the leading wave front and the stresses immediately behind it. Experimental results for the attenuation of stress wave on steel-aluminum and steel-brass fiber-matrix composites compared very well with the computed analytical results when the applied pressure is generated by small explosive charges. The results did not compare well when the applied pressure is generated by projectile impact.

Yang, J. C. S.↗

Free-Vibration Analysis of Structures

Improved numerical procedure more than twice as fast as previous methods. Unified numerical algorithm efficiently solves free-vibration problems of stationary or spinning structures with or without viscous or structural damping. Algorithm used to solve static problems involving multiple loads and to solve quadratic matrix eigenvalue problems associated with finite-dynamic-element structural discretization.

Gupta, K. K.↗

MatLab Script and Functional Programming

MatLab Script and Functional Programming: MatLab is one of the most widely used very high level programming languages for scientific and engineering computations. It is very user-friendly and needs practically no formal programming knowledge. Presented here are MatLab programming aspects and not just the MatLab commands for scientists and engineers who do not have formal programming training and also have no significant time to spare for learning programming to solve their real world problems. Specifically provided are programs for visualization. The MatLab seminar covers the functional and script programming aspect of MatLab language. Specific expectations are: a) Recognize MatLab commands, script and function. b) Create, and run a MatLab function. c) Read, recognize, and describe MatLab syntax. d) Recognize decisions, loops and matrix operators. e) Evaluate scope among multiple files, and multiple functions within a file. f) Declare, define and use scalar variables, vectors and matrices.

Shaykhian, Gholam Ali↗

Evaluation of Higher-order Quadrature Schemes in Improving Computational Efficiency for Orientation-averaged Single-Scattering Properties of Nonspherical Ice Particles

We evaluate several high-order quadrature schemes for accuracy and efficacy in obtaining orientation-averaged single-scattering properties (SSPs). We use the highly efficient MIDAS to perform electromagnetic scattering calculations to evaluate the gain in efficiency from these schemes. MIDAS is shown to be superior to DDSCAT, a popular discrete dipole approximation (DDA) method. This study is motivated by the fact that quality physical precipitation retrievals rely on using accurate orientation-averaged SSPs derived from realistic hydrometeors as input to radiative transfer Models (RTMs). The DDA has been a popular choice for single-scattering calculations, due to its versatility with respect to target geometry. However, being iterative-solver-based (ISB), the most used DDA codes, e.g. DDSCAT and ADDA, must solve the scattering problem for each orientation of the target separately. As the size parameter and geometric anisotropy of the hydrometeor increase, the number of orientations needed to obtain accurate orientation-averages can increase drastically and so does the computation cost incurred by the ISB-DDA methods. MIDAS is a Direct-Solver-Based (DSB) code, its decomposition of the original large matrix with a high rank into multiple more manageable smaller matrices of lower ranks makes it much more computationally efficient while maintaining excellent accuracy. In addition, direct solvers consider all requested orientations at once, giving MIDAS further advantage over popular ISB-DDA methods. MIDAS, when combined with high-order quadrature for orientation averaging, can be greater than three orders of magnitude more efficient in obtaining RTM-ready SSPs of complex-shaped hydrometeors than existing ISB-DDA methods, with the native quadrature schemes they offer.

Ines Fenni↗

Effect of Thin-Film Adhesives on Mode II Interlaminar Fracture Toughness in Carbon Fiber Composites with Shape Memory Alloy Inserts

A single sheet of nickel-titanium (NiTi) shape memory alloy (SMA) was introduced within an IM7/8552 polymer matrix composite (PMC) panel in conjunction with multiple thin film adhesives to promote the interfacial bond strength between the SMA and PMC. End notched flexure (ENF) testing was performed in accordance to ASTM D7905 method for evaluation of mode II interlaminar fracture toughness (GIIC) of unidirectional fiber-reinforced polymer matrix composites. Acoustic emissions (AE) were monitored during testing with two acoustic sensors attached to the specimens. The composite panels were subjected to C-scan before testing, and examined using optical and scanning electron microscopy techniques after part failure. GIIC values for the control composite samples were found to be higher than those of samples with embedded SMA sheets. The presence of adhesives bonded to SMA sheets further diminished the GIIC values. AE values revealed poor bonding of the panels, with little to no signals during testing.

End Notched Flexure↗

The Effect of Thin Film Adhesives on Mode II Interlaminar Fracture Toughness in Carbon Fiber Composites with Shape Memory Alloy Inserts

A single sheet of nickel-titanium (NiTi) shape memory alloy (SMA) was introduced within an IM7/8552 polymer matrix composite (PMC) panel in conjunction with multiple thin film adhesives to promote the interfacial bond strength between the SMA and PMC. End notched flexure (ENF) testing was performed in accordance to ASTM D7905 method for evaluation of mode II interlaminar fracture toughness (GIIC) of unidirectional fiber-reinforced polymer matrix composites. Acoustic emissions (AE) were monitored during testing with two acoustic sensors attached to the specimens. The composite panels examined using scanning electron microscopy techniques after part failure. GIIC values for the control composite samples were found to be higher than those of samples with embedded SMA sheets. The presence of adhesives bonded to SMA sheets further diminished the GIIC values. AE values revealed poor bonding of the panels, with little to no signals during testing.

End Notched Flexure↗