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 307 records · Page 17

Analyzing Multirate-Sampled Systems

New method available for monitoring stability of system on basis of data sampled at multiple rates--in particular, at two rates, one of which is twice frequency of other. Method is called Multirate Matrix Frequency Response (MMFR) analysis. MMFR was used in Space Shuttle to evaluate effect of decreasing sample rate of error loop in ascent-phase digital autopilot. Should also be useful to designers of other control systems and to structural, civil, and mechanical engineers for structure and vibration analysis.

Hendrix, N.↗

A fracture mechanics analysis for delamination growth during impact on composite plates

A fracture mechanics analysis has been developed that describes the progress of delamination damage in composite plates struck by a hard spherical object. The analysis is based on large deflection plate mechanics for circular isotropic plates wherein multiple axisymmetric delaminations grow. Test data show that the analysis predicts the influence of plate thickness, support conditions, and matrix toughness on the onset and propagation of delamination.

Bostaph, G. M.↗

Diffraction Modeling of Finite Subband EFC Probing on Dark Hole Contrast with WFIRST-CGI Shaped Pupil Coronagraph

The current NASA WFIRST (Wide-Field InfraRed Survey Telescope) coronagraph instrument (CGI) design allocates two subband filters per each full science band in order to contain system complexity and cost. We present our detailed investigation results on the adequacy of such a limited number of finite subband filters in achieving broadband dark hole contrast with a shaped pupil coronagraph (SPC). The study is based on diffraction propagation modeling with realistic WFIRST optics, where each subband’s image plane electric field is estimated from pairwise pupil plane deformable mirror (DM) probing and image plane intensity averaging of the resulting fields of multiple (subband) wavelengths. Multiple subband choices and probing and control strategies are explored, including standard subband probing; mixed wavelength and/or weighted Jacobian matrix; subband probing with intensity subtraction; and extended subband probing with intensity subtraction. Overall, the investigation shows that the achievable contrast with a limited number of finite subbands is about 2~2.5x worse than the designed contrast for the current SPC. The result suggests that future shaped pupil designs should be optimized for slightly broader bandwidths than the intended science bandpasses if limited subbands are used for wavefront sensing via probing.

WFIRST↗

The Development of a Real-Time Optical Angle-of-Attack Measurement Capability at the NASA Ames Unitary Plan Wind Tunnel

The following details the implementation of a real-time optical angle-of-attack capability under development at the NASA Ames Unitary Plan Wind Tunnel (UPWT). Recent advancements in the integration of high bandwidth imaging systems at the Ames UPWT has laid the groundwork for the deployment of a calibrated, imaging based measurement capability. Some of the measurement goals of this data system include: an estimation of the orientation of a wind tunnel model to six degrees of freedom, measurements made at framerates high enough to time resolve model dynamics, triangulate corresponding points of interest from multiple camera views in three-dimensional space, acquire imagery with no impact to the productivity of the test matrix, provide results in real-time, and finally compute an estimation of the errors associated with the measurement. The following will outline the imaging hardware and data systems architecture, describe the data flow and image processing routines currently implemented, and explore example data and results from previous wind tunnel test entries.

Ground Testing↗

Implementation of a digital optical matrix-vector multiplier using a holographic look-up table and residue arithmetic

The design and implementation of a digital (numerical) optical matrix-vector multiplier are presented. The objective is to demonstrate the operation of an optical processor designed to minimize computation time in performing a practical computing application. This is done by using the large array of processing elements in a Hughes liquid crystal light valve, and relying on the residue arithmetic representation, a holographic optical memory, and position coded optical look-up tables. In the design, all operations are performed in effectively one light valve response time regardless of matrix size. The features of the design allowing fast computation include the residue arithmetic representation, the mapping approach to computation, and the holographic memory. In addition, other features of the work include a practical light valve configuration for efficient polarization control, a model for recording multiple exposures in silver halides with equal reconstruction efficiency, and using light from an optical fiber for a reference beam source in constructing the hologram. The design can be extended to implement larger matrix arrays without increasing computation time.

Habiby, Sarry F.↗

Infrared spectoscopy of Triton and Pluto ice analogs: The case for saturated hydrocarbons

The infrared transmission spectra and photochemical behavior of various organic compounds isolated in solid N2 ices, appropriate for applications to Triton ad Pluto, are presented. It is shown that excess absorption in the surface spectra of Triton and Pluto, i.e., absorption not explained by present models incorporating molecules already identified on these bodies (N2, CH4, CO, and CO2), that starts near 4450/cm (2.25 microns) and extends to lower frequencies, may be due to alkanes (C(n)H(2n+2)) and related molecules frozen in the nitrogen. Branched and linear alkanes may be responsible. Experiments in which the photochemstry of N2: CH4 and N2: CH4: CO ices was explored demonsrtrate that the surface ices of Triton and Pluto may contain a wide variety of additional species containing H, C, O, and N. Of these, the reactive molecule diazomethane, CH2N2, is particularly important since it may be largely responsible for the synthesis of larger alkanes from CH4 and other small alkanes. Diazomethane would also be expected to drive chemical reactions involving organics in the surface ices of Triton and Pluto toward saturation, i.e., to reduce multiple CC bonds. The positions and intrinsic strengths (A values) of many of the infrared absorption bands of N2 matrix-isolated molecules of relevance to Triton and Pluto have also been determined. These can be used to aid in their search and to place constraints on their abundances.

Bohn, Robert B.↗

Attitude Estimation Using Modified Rodrigues Parameters

In this paper, a Kalman filter formulation for attitude estimation is derived using the Modified Rodrigues Parameters. The extended Kalman filter uses a gyro-based model for attitude propagation. Two solutions are developed for the sensitivity matrix in the Kalman filter. One is based upon an additive error approach, and the other is based upon a multiplicative error approach. It is shown that the two solutions are in fact equivalent. The Kalman filter is then used to estimate the attitude of a simulated spacecraft. Results indicate that then new algorithm produces accurate attitude estimates by determining actual gyro biases.

Crassidis, John L.↗

Shape and Structure

The convenient availability and simplicity of the Lorenz-Mie theory has resulted in a widespread practice of treating nonspherical particles (especially those in random orientation) as if they were spheres to which Lorenz-Mie results are applicable. However, the assumption of sphericity is rarely made after first having studied the effects of nonsphericity and having concluded that they are negligible. On the contrary, overwhelming evidence suggests that scattering properties of nonspherical particles, including those in random orientation, can significantly differ from those of volume- or surface-equivalent spheres. Hence, the last few decades have demonstrated major research efforts aimed at a significantly better understanding of the effects of particle shape and morphology on electromagnetic scattering. The goal of this presentation is to provide a concise summary of these efforts. The recent availability of theoretical techniques for computing single and multiple scattering of light by realistic polydispersions of spherical and nonspherical particles and the strong dependence of the Stokes scattering matrix on particle size, shape, and refractive index make polarization and depolarization measurements a powerful particle characterization tool. This presentation will focus on recent applications of photopolarimetric and lidar depolarization measurements to remote sensing characterization of tropospheric aerosols, polar stratospheric clouds (PSCs) and contrails. The talk will include (1) a short theoretical overview of the effects of particle microphysics on particle single-scattering characteristics; (2) the use of multi-angle multi-spectral photopolarimetry to retrieve the optical thickness, size distribution, refractive index, and number concentration of tropospheric aerosols over the ocean surface; and (3) the application of the T-matrix method to constraining the PSC and contrail particle microphysics using multi-spectral measurements of lidar backscatter and depolarization

Mishchenko, Michael↗

Retrievals of Aerosol and Cloud Particle Microphysics Using Polarization and Depolarization Techniques

The recent availability of theoretical techniques for computing single and multiple scattering of light by realistic polydispersions of spherical and nonspherical particles and the strong dependence of the Stokes scattering matrix on particle size, shape, and refractive index make polarization and depolarization measurements a powerful particle characterization tool. In this presentation I will describe recent applications of photopolarimetric and lidar depolarization measurements to remote sensing characterization of tropospheric aerosols, polar stratospheric clouds (PSCs), and contrails. The talk will include (1) a short theoretical overview of the effects of particle microphysics on particle single-scattering characteristics; (2) the use of multi-angle multi-spectral photopolarimetry to retrieve the optical thickness, size distribution, refractive index, and number concentration of tropospheric aerosols over the ocean surface; and (3) the application of the T-matrix method to constraining the PSC and contrail particle microphysics using multi-spectral measurements of lidar backscatter and depolarization.

Mishchenko, Michael↗

Genetic Optimization of a Tensegrity Structure

Marshall Space Flight Center (MSFC) is charged with developing advanced technologies for space telescopes. The next generation of space optics will be very large and lightweight. Tensegrity structures are built of compressive members (bars), and tensile members (strings). For most materials, the tensile strength of a longitudinal member is larger than its buckling strength; therefore a large stiffness to mass ratio can be achieved by increasing the use of tensile members. Tensegrities are the epitome of lightweight structures, since they take advantage of the larger tensile strength of materials. The compressive members of tensegrity structures are disjoint allowing compact storage of the structure. The structure has the potential to eliminate the requirement for assembly by man in space; it can be deployed by adjustments in its cable tension. A tensegrity structure can be more reliably modeled since none of the individual members experience bending moments. (Members that experience deformation in more than one dimension are much harder to model.) Structures that can be more precisely modeled can be more precisely controlled. Furthermore, an astoundingly wide variety of natural systems, including carbon atoms, water molecules, proteins, viruses, cells, tissues and even human and other living creatures are tensegrity structures. Through the process of evolution, nature continually improves the design of living creatures for the environment they live in. Since tensegrities are nature's structure of choice, it is conceivable that they have other benefits we are unaware of. A. Keane and S. Brown designed a satellite boom truss system with an enhanced vibration performance. They started with a standard truss system, then used a genetic algorithm to alter the design, while optimizing the vibration performance. An improvement of over 20,000% in frequency-averaged energy levels was obtained using this approach. In this report an introduction to tensegrity structures is given, along with a description of how to generate the nodal coordinates and connectivity of a multiple stage cylindrical tensegrity structure. A description of how finite elements can be used to develop a stiffness and mass matrix so that the modes of vibration can be determined from the eigenvalue problem is shown. A brief description of a micro genetic algorithm is then presented.

Jaime R. Taylor↗

Retrieval of Atmospheric and Surface Parameters from AIRS/AMSU/HSB Data Under Cloudy Conditions

New state of the art methodology is described to analyze AIRS/AMSU/HSB data in the presence of multiple cloud formations. The methodology forms the basis for the AIRS Science Team algorithm which will be used to analyze AIRS/AMSU/HSB data on EOS Aqua. The cloud clearing methodology requires no knowledge of the spectral properties of the clouds. The basic retrieval methodology is general and extracts the maximum information from the radiances, consistent with the channel noise covariance matrix. The retrieval methodology minimizes the dependence of the solution on the first guess field and the first guess error characteristics. Results are shown for AIRS Science Team simulation studies with multiple cloud formations. These simulation studies imply that clear column radiances can be reconstructed under partial cloud cover with an accuracy comparable to single spot channel noise in the temperature and water vapor sounding regions, temperature soundings can be produced under partial cloud cover with RMS errors on the order of, or better than, 1deg K in 1 km thick layers from the surface to 700 mb, 1 km layers from 700 mb to 300 mb, 3 km layers from 300 mb to 30 mb, and 5 km layers from 30 mb to 1 mb, and moisture profiles can be obtained with an accuracy better than 20% absolute errors in 1 km layers from the surface to nearly 200 mb.

Susskind, Joel↗

Approaches to optimization of SS/TDMA time slot assignment

Reduction techniques for traffic matrices are explored in some detail. These matrices arise in satellite switched time-division multiple access (SS/TDMA) techniques whereby switching of uplink and downlink beams is required to facilitate interconnectivity of beam zones. A traffic matrix is given to represent that traffic to be transmitted from n uplink beams to n downlink beams within a TDMA frame typically of 1 ms duration. The frame is divided into segments of time and during each segment a portion of the traffic is represented by a switching mode. This time slot assignment is characterized by a mode matrix in which there is not more than a single non-zero entry on each line (row or column) of the matrix. Investigation is confined to decomposition of an n x n traffic matrix by mode matrices with a requirement that the decomposition be 100 percent efficient or, equivalently, that the line(s) in the original traffic matrix whose sum is maximal (called critical line(s)) remain maximal as mode matrices are subtracted throughout the decomposition process. A method of decomposition of an n x n traffic matrix by mode matrices results in a number of steps that is bounded by n(2) - 2n + 2. It is shown that this upper bound exists for an n x n matrix wherein all the lines are maximal (called a quasi doubly stochastic (QDS) matrix) or for an n x n matrix that is completely arbitrary. That is, the fact that no method can exist with a lower upper bound is shown for both QDS and arbitrary matrices, in an elementary and straightforward manner.

Wade, T. O.↗

How well can gravity be recovered using Topex and GPS data?

When Topex is launched in mid-1992 it will carry a high quality GPS receiver which will operate in concert with a worldwide network of precision GPS ground receivers. The data from these receivers can be used to recover new information about the earth's gravity field at longer wavelengths. Software and algorithms have been developed which will allow this gravity field information to be recovered with much greater efficiency than with traditional techniques. The basis for these algorithms is the gravity bin formulation and related filtering techniques that exploit the repeat orbit of Topex and the sparse matrix structure of the problem. This new software has been used to evaluate the expected improvement in the gravity field using multiple ten-day arcs of GPS data from Topex.

Bertiger, Willy↗

Reflection and emission properties on Venus - Alpha Regio

Procedures developed from Hagfors' (1964) scattering model have been used with success to extract various properties (reflectivity, emissivity, slope, etc.) from both Pioneer Venus and Magellan data. However, under certain circumstances, such as physically rough areas or areas of nonhomogeneous composition, Hagfors' model may lead to serious errors in the interpretation of radar and radiometry data. The Magellan-derived surface properties are examined for some areas of below average emissivity in Alpha Regio and in an area to the east of Alpha, Crater Stuart, in order to determine if Hagfors' model is applicable. Methods of classical multiple scattering theory are also applied to show that enhanced radar backscattering and depressed thermal emission can be expected with a matrix of silicic or mafic rocks embedded in low loss soils.

Tryka, Kimberly A.↗

Component modes damping assignment methodology for articulated, multiflexible body structures

To simulate the dynamical motion of articulated, multiflexible body structures, one can use multibody simulation packages such as DISCOS. To this end, one must supply appropriate reduced-order models for all of the flexible components involved. The component modes projection and assembly model reduction (COMPARE) methodology is one way to construct these reduced-order component models, which when reassembled capture important system input-to-output mapping of the full-order model at multiple system configurations of interest. In conjunction, we must also supply component damping matrices which when reassembled generate a system damping matrix that has certain desirable properties. The problem of determining the damping factors of components' modes to achieve a given system damping matrix is addressed here. To this end, we must establish from first principles a matrix-algebraic relation between the system's modal damping matrix and the components' modal damping matrices. An unconstrained/constrained optimization problem can then be formulated to determine the component modes' damping factors that best satisfy that matrix-algebraic relation. The effectiveness of the developed methodology, called ModeDamp, has been successfully demonstrated on a high-order, finite element model of the Galileo spacecraft.

Lee, Allan Y.↗

Automated Guided-Wave Scanning Developed to Characterize Materials and Detect Defects

The Nondestructive Evaluation (NDE) Group of the Optical Instrumentation Technology Branch at the NASA Glenn Research Center has developed a scanning system that uses guided waves to characterize materials and detect defects. The technique uses two ultrasonic transducers to interrogate the condition of a material. The sending transducer introduces an ultrasonic pulse at a point on the surface of the specimen, and the receiving transducer detects the signal after it has passed through the material. The aim of the method is to correlate certain parameters in both the time and frequency domains of the detected waveform to characteristics of the material between the two transducers. The scanning system is shown. The waveform parameters of interest include the attenuation due to internal damping, waveform shape parameters, and frequency shifts due to material changes. For the most part, guided waves are used to gauge the damage state and defect growth of materials subjected to various mechanical or environmental loads. The technique has been applied to polymer matrix composites, ceramic matrix composites, and metal matrix composites as well as metallic alloys. Historically, guided wave analysis has been a point-by-point, manual technique with waveforms collected at discrete locations and postprocessed. Data collection and analysis of this type limits the amount of detail that can be obtained. Also, the manual movement of the sensors is prone to user error and is time consuming. The development of an automated guided-wave scanning system has allowed the method to be applied to a wide variety of materials in a consistent, repeatable manner. Experimental studies have been conducted to determine the repeatability of the system as well as compare the results obtained using more traditional NDE methods. The following screen capture shows guided-wave scan results for a ceramic matrix composite plate, including images for each of nine calculated parameters. The system can display up to 18 different wave parameters. Multiple scans of the test specimen demonstrated excellent repeatability in the measurement of all the guided-wave parameters, far exceeding the traditional point-by-point technique. In addition, the scan was able to detect a subsurface defect that was confirmed using flash thermography This technology is being further refined to provide a more robust and efficient software environment. Future hardware upgrades will allow for multiple receiving transducers and the ability to scan more complex surfaces. This work supports composite materials development and testing under the Ultra-Efficient Engine Technology (UEET) Project, but it also will be applied to other material systems under development for a wide range of applications.

Martin, Richard E.↗

Radiative interaction of atmosphere and surface: write up with elements of code

In passive satellite remote sensing of the Earth, separation of the path radiance (atmosphere-only contribution) from the surface reflection remains a “significant challenge”. Recent literature names it among the gaps in radiative transfer (RT) topics that “require continued research in the near future”. The challenge comes from multiple reflections (bouncing) between the atmosphere and surface – radiative interaction. In this paper we use a known RT technique, the matrix-operator method (MOM), and a new modification of the monochromatic vector RT (vRT) code IPOL (Intensity and POLarization) to simulate the interaction of a plane-parallel atmosphere and a few widely used surface reflection models. Following the idea of the Green’s function method, IPOL no longer takes the surface model parameters on input. Instead, it provides the path radiance, and the atmospheric reflection and transmission matrices as output. Despite many RT codes use the MOM formalism, this output does not seem common. The surface reflection matrix is computed externally. Therefore, this paper extends the Green’s function atmospheric correction technique to the case of polarized light. Aiming clarity rather than performance, we explain in Python the structure of the surface matrices for the isotropic (Lambertian), directional unpolarized, and polarized ocean reflection models. We then combine these surface matrices and the precomputed IPOL output to get numerically accurate signal at the top of atmosphere (TOA) and test it vs. published benchmarks. Then, for each benchmark scenario we show how to get the surface from the TOA signal, i.e. perform the RT-based atmospheric correction.

radiative transfer↗

Multi-Scale Thermo-Mechanical Modeling of Porous 3D Woven TPS Materials

This work summarizes the process to compute and analyze the thermal conductivity and mechanical properties of 3D woven TPS materials such as 3MDCP (3D Mid-Density Carbon-Phenolic). The PuMA [1,2] software, developed at NASA Ames, was used to characterize different 3MDCP samples from their constituents' data, averaging their thermal conductivity and elastic properties in the three main directions to obtain their effective orthotropic thermo-mechanical properties. This was performed in multiple steps: firstly, TPS samples were digitally reconstructed using micro computed tomography (µCT) and their constituents were segmented; then, the porous matrix phase was analyzed at the micro-scale and these results were used, along with the fibers’ constituents information, to model the tows at the meso-scale; finally, results from the constituents were homogenized and used to model the thermo-mechanical behavior of the unit cell at the macro-scale. Additionally, to gain a better understanding of the tows’ morphology and distribution of the carbon-phenolic blended fibers, microscopy images of the tows’ cross-section were segmented using deep learning techniques and analyzed with PuMA. This workflow can be applied to any TPS material to computationally obtain its thermo-mechanical properties, enabling more informed TPS design and manufacturing choices.

Multi-Scale↗