Reconfigurable Control for the Formation Flying of Multiple Spacecraft
Several results on the reconfigurable control architecture for the formation flying of multiple spacecraft are presented.
formation↗
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.
Distributed Sensing of Carbon-Epoxy Composites and Filament Wound Pressure Vessels Using Fiber-Bragg Gratings
Multiple Fiber Bragg-gratings are embedded in carbon-epoxy laminates as well as in composite wound pressure vessel. Structural properties of such composites are investigated. The measurements include stress-strain relation in laminates and Poisson's ratio in several specimens with varying orientation of the optical fiber Bragg-sensor with respect to the carbon fiber in an epoxy matrix. Additionally, fiber Bragg gratings are bonded on the surface of these laminates and cylinders fabricated out of carbon-epoxy composites and multiple points are monitored and compared for strain measurements at several locations.
Grant, J.↗
Multiple-Flat-Panel System Displays Multidimensional Data
The NASA Ames hyperwall is a display system designed to facilitate the visualization of sets of multivariate and multidimensional data like those generated in complex engineering and scientific computations. The hyperwall includes a 77 matrix of computer-driven flat-panel video display units, each presenting an image of 1,280 1,024 pixels. The term hyperwall reflects the fact that this system is a more capable successor to prior computer-driven multiple-flat-panel display systems known by names that include the generic term powerwall and the trade names PowerWall and Powerwall. Each of the 49 flat-panel displays is driven by a rack-mounted, dual-central-processing- unit, workstation-class personal computer equipped with a hig-hperformance graphical-display circuit card and with a hard-disk drive having a storage capacity of 100 GB. Each such computer is a slave node in a master/ slave computing/data-communication system (see Figure 1). The computer that acts as the master node is similar to the slave-node computers, except that it runs the master portion of the system software and is equipped with a keyboard and mouse for control by a human operator. The system utilizes commercially available master/slave software along with custom software that enables the human controller to interact simultaneously with any number of selected slave nodes. In a powerwall, a single rendering task is spread across multiple processors and then the multiple outputs are tiled into one seamless super-display. It must be noted that the hyperwall concept subsumes the powerwall concept in that a single scene could be rendered as a mosaic image on the hyperwall. However, the hyperwall offers a wider set of capabilities to serve a different purpose: The hyperwall concept is one of (1) simultaneously displaying multiple different but related images, and (2) providing means for composing and controlling such sets of images. In place of elaborate software or hardware crossbar switches, the hyperwall concept substitutes reliance on the human visual system for integration, synthesis, and discrimination of patterns in complex and high-dimensional data spaces represented by the multiple displayed images. The variety of multidimensional data sets that can be displayed on the hyperwall is practically unlimited. For example, Figure 2 shows a hyperwall display of surface pressures and streamlines from a computational simulation of airflow about an aerospacecraft at various Mach numbers and angles of attack. In this display, Mach numbers increase from left to right and angles of attack increase from bottom to top. That is, all images in the same column represent simulations at the same Mach number, while all images in the same row represent simulations at the same angle of attack. The same viewing transformations and the same mapping from surface pressure to colors were used in generating all the images.
Gundo, Daniel↗
Micromechanical Modeling of the Thermal Expansion of Graphite/copper Composites with Nonuniform Microstructure
Two micromechanical models were developed to investigate the thermal expansion of graphite/copper (Gr/Cu) composites. The models incorporate the effects of temperature-dependent material properties, matrix inelasticity, initial residual stresses due to processing history, and nonuniform fiber distribution. The first model is based on the multiple concentric cylinder geometry, with each cylinder treated as a two-phase composite with a characteristic fiber volume fractions. By altering the fiber volume fraction of the individual cylinders, unidirectional composites with radially nonuniform fiber distributions can be investigated using this model. The second model is based on the inelastic lamination theory. By varying the fiber content in the individual laminae, composites with nonuniform fiber distribution in the thickness direction can be investigated. In both models, the properties of the individual regions (cylinders or laminae) are calculated using the method of cells micromechanical model. Classical incremental plasticity theory is used to model the inelastic response of the copper matrix at the microlevel. The models were used to characterize the effects of nonuniform fiber distribution on the thermal expansion of Gr/Cu. These effects were compared to the effects of matrix plasticity, choice of stress-free temperature, and slight fiber misalignment. It was found that the radially nonuniform fiber distribution has little effect on the thermal expansion of Gr/Cu but could become significant for composites with large fiber-matrix transverse CTE and Young's modulus mismatch. The effect of nonuniform fiber distribution in the through-thickness direction of a laminate was more significant, but only approached that of the stress-free temperature for the most extreme cases that include large amounts of bending. Subsequent comparison with experimental thermal expansion data indicated the need for more accurate characterization of the graphite fiber thermomechanical properties. Correlation with cyclic data revealed the presence of a mechanism not considered in the developed models. The predicted response did, however, exhibit ratcheting behavior that has been observed experimentally in Gr/Cu. Finally, simulation of the actual fiber distribution of particular specimens had little effect on the predicted thermal expansion.
Bednarcyk, Brett A.↗
North Alabama Lightning Mapping Array (LMA): VHF Source Retrieval Algorithm and Error Analyses
Two approaches are used to characterize how accurately the North Alabama Lightning Mapping Array (LMA) is able to locate lightning VHF sources in space and in time. The first method uses a Monte Carlo computer simulation to estimate source retrieval errors. The simulation applies a VHF source retrieval algorithm that was recently developed at the NASA Marshall Space Flight Center (MSFC) and that is similar, but not identical to, the standard New Mexico Tech retrieval algorithm. The second method uses a purely theoretical technique (i.e., chi-squared Curvature Matrix Theory) to estimate retrieval errors. Both methods assume that the LMA system has an overall rms timing error of 50 ns, but all other possible errors (e.g., multiple sources per retrieval attempt) are neglected. The detailed spatial distributions of retrieval errors are provided. Given that the two methods are completely independent of one another, it is shown that they provide remarkably similar results. However, for many source locations, the Curvature Matrix Theory produces larger altitude error estimates than the (more realistic) Monte Carlo simulation.
Koshak, W. J.↗
A high-accuracy optical linear algebra processor for finite element applications
Optical linear processors are computationally efficient computers for solving matrix-matrix and matrix-vector oriented problems. Optical system errors limit their dynamic range to 30-40 dB, which limits their accuray to 9-12 bits. Large problems, such as the finite element problem in structural mechanics (with tens or hundreds of thousands of variables) which can exploit the speed of optical processors, require the 32 bit accuracy obtainable from digital machines. To obtain this required 32 bit accuracy with an optical processor, the data can be digitally encoded, thereby reducing the dynamic range requirements of the optical system (i.e., decreasing the effect of optical errors on the data) while providing increased accuracy. This report describes a new digitally encoded optical linear algebra processor architecture for solving finite element and banded matrix-vector problems. A linear static plate bending case study is described which quantities the processor requirements. Multiplication by digital convolution is explained, and the digitally encoded optical processor architecture is advanced.
Casasent, D.↗
Implementation of a block Lanczos algorithm for Eigenproblem solution of gyroscopic systems
The details of implementation of a general numerical procedure developed for the accurate and economical computation of natural frequencies and associated modes of any elastic structure rotating along an arbitrary axis are described. A block version of the Lanczos algorithm is derived for the solution that fully exploits associated matrix sparsity and employs only real numbers in all relevant computations. It is also capable of determining multiple roots and proves to be most efficient when compared to other, similar, exisiting techniques.
Gupta, Kajal K.↗
Implementation of a block Lanczos algorithm for eigenproblem solution of gyroscopic systems
This paper describes the details of implementation of a general numerical procedure developed for the accurate and economical computation of natural frequencies and associated modes of any elastic structure rotating along an arbitrary axis. A block version of the Lanczos algorithm is derived for the solution that fully exploits associated matrix sparsity and employs only real numbers in all relevant computations. It is also capable of determining multiple roots and proves to be most efficient when compared to other, similar, existing techniques.
Gupta, K. K.↗
Determining water cloud particle sizes from lidar depolarization measurements and time dependent multiple scattering coefficients
We present a technique of extracting water cloud particle size information from lidar measurements in conjunction with double scattering calculations. In this presentation, we describe the technique and give examples using data taken with the Air Force Phillips Laboratory's (Geophysics Directorate) low altitude Nd:YAG, elastic backscatter lidar. In a related presentation we describe the double scattering lidar model which we developed for this work. The technique uses simultaneous measurements of two orthogonal linear polarization components of lidar returns from water clouds or other media composed of spherical particles. Any depolarization of the incident lidar radiation backscattered by such a media can only be due to multiple scattering. The amount of depolarization is dependent on the extinction coefficient and the single scatter phase matrix, both of which are functions of position in the medium. The phase matrix is dependent on the index of refraction of the particles and the particle size distribution. Our technique is a modification of a procedure presented in Reference 1. There, particle sizes of water clouds are determined from double scattering calculations together with measurements of radiation scattered from volumes outside the lidar receiver's field of view (which can only be multiply scattered radiation). The methodology of our technique is similar but our 'probe' of the scattering phase function (and thus the particle size distribution) is different.
Garner, R. C.↗
Measuring attitude with a gradiometer
This paper explores using a gravity gradiometer to measure the attitude of a satellite, given that the gravity field is accurately known. Since gradiometers actually measure a combination of the gradient and attitude rate and acceleration terms, the answer is far from obvious. The paper demonstrates that it can be done and at microradian accuracy. The technique employed is dynamic estimation, based on the momentum biased Euler equations. The satellite is assumed nominally planet pointed, and subject to control, gravity gradient, and partly radom drag torques. The attitude estimator is unusual. While the standard method of feeding back measurement residuals is used, the feedback gain matrix isn't derived from Kalman theory. instead, it's chosen to minimize a measure of the terminal covariance of the error in the estimate. This depends on the gain matrix and the power spectra of all the process and measurement noises. An integration is required over multiple solutions of Lyapunov equations.
Sonnabend, David↗
Model-Based Fatigue Prognosis of Fiber-Reinforced Laminates Exhibiting Concurrent Damage Mechanisms
Prognostics of large composite structures is a topic of increasing interest in the field of structural health monitoring for aerospace, civil, and mechanical systems. Along with recent advancements in real-time structural health data acquisition and processing for damage detection and characterization, model-based stochastic methods for life prediction are showing promising results in the literature. Among various model-based approaches, particle-filtering algorithms are particularly capable in coping with uncertainties associated with the process. These include uncertainties about information on the damage extent and the inherent uncertainties of the damage propagation process. Some efforts have shown successful applications of particle filtering-based frameworks for predicting the matrix crack evolution and structural stiffness degradation caused by repetitive fatigue loads. Effects of other damage modes such as delamination, however, are not incorporated in these works. It is well established that delamination and matrix cracks not only co-exist in most laminate structures during the fatigue degradation process but also affect each other's progression. Furthermore, delamination significantly alters the stress-state in the laminates and accelerates the material degradation leading to catastrophic failure. Therefore, the work presented herein proposes a particle filtering-based framework for predicting a structure's remaining useful life with consideration of multiple co-existing damage-mechanisms. The framework uses an energy-based model from the composite modeling literature. The multiple damage-mode model has been shown to suitably estimate the energy release rate of cross-ply laminates as affected by matrix cracks and delamination modes. The model is also able to estimate the reduction in stiffness of the damaged laminate. This information is then used in the algorithms for life prediction capabilities. First, a brief summary of the energy-based damage model is provided. Then, the paper describes how the model is embedded within the prognostic framework and how the prognostics performance is assessed using observations from run-to-failure experiments
Prognostics↗
Stochastic Nonlinear Response of Woven CMCs
It is well known that failure of a material is a locally driven event. In the case of ceramic matrix composites (CMCs), significant variations in the microstructure of the composite exist and their significance on both deformation and life response need to be assessed. Examples of these variations include changes in the fiber tow shape, tow shifting/nesting and voids within and between tows. In the present work, the influence of scale specific architectural features of woven ceramic composite are examined stochastically at both the macroscale (woven repeating unit cell (RUC)) and structural scale (idealized using multiple RUCs). The recently developed MultiScale Generalized Method of Cells methodology is used to determine the overall deformation response, proportional elastic limit (first matrix cracking), and failure under tensile loading conditions and associated probability distribution functions. Prior results showed that the most critical architectural parameter to account for is weave void shape and content with other parameters being less in severity. Current results show that statistically only the post-elastic limit region (secondary hardening modulus and ultimate tensile strength) is impacted by local uncertainties both at the macro and structural level.
Kuang, C. Liu↗
A probabilistic model for the evaluation of fault-tolerant multiprocessor systems using concurrent error detection
A probabilistic model to evaluate fault-tolerant multiprocessor systems has been developed. The matrix-based model and the analysis algorihtms based on it are described. Various probabilities associated with an algorithm-based fault tolerance system are discussed and the fault coverage of a given check is derived analytically and illustrated with examples. The probability matrices that are formed by introducing the spatial probabilities into the matrix model are considered. Based on these matrices, a technique is developed to determine the combined coverage of multiple numbers of checks. Examples of the analysis of systems using the model are given.
Nair, V. S. S.↗
Software Computes Tape-Casting Parameters
Tcast2 is a FORTRAN computer program that accelerates the setup of a process in which a slurry containing metal particles and a polymeric binder is cast, to a thickness regulated by a doctor blade, onto fibers wound on a rotating drum to make a green precursor of a metal-matrix/fiber composite tape. Before Tcast2, setup parameters were determined by trial and error in time-consuming multiple iterations of the process. In Tcast2, the fiber architecture in the final composite is expressed in terms of the lateral distance between fibers and the thickness-wise distance between fibers in adjacent plies. The lateral distance is controlled via the manner of winding. The interply spacing is controlled via the characteristics of the slurry and the doctor-blade height. When a new combination of fibers and slurry is first cast and dried to a green tape, the shrinkage from the wet to the green condition and a few other key parameters of the green tape are measured. These parameters are provided as input to Tcast2, which uses them to compute the doctor-blade height and fiber spacings needed to obtain the desired fiber architecture and fiber volume fraction in the final composite.
deGroh, Henry C., III↗
Formation Flying Control of Multiple Spacecraft: Graph Theoretic Properties and Switching Schemes
We present several results on the leader-following (LF) paradigm in the formation flying of multiple spacecraft in free space.
Formation↗
Results on the two population feature selection problem using probability of correct classification as a criterion
Variational equations are presented for maximizing the probability of correct classification as a function of a 1xn feature selection matrix B for the two-population problem. For the special case of equal covariance matrices the optimal B is unique up to scalar multiples and rank one sufficient. For equal population means, the best 1xn B is an eigenvector corresponding either to the largest or smallest eigenvalue of sigma sub 2 to the minus 1 power sigma sub 1 where sigma sub 1 and sigma sub 2 are the nxn covariance matrices of the two populations. The transformed probability of correct classification depends only on the eigenvalue. Finally, a procedure is proposed for constructing an optimal or nearly optimal kxn matrix of rank k without solving the k-dimensional variational equation.
Peters, B. C.↗
Degree and plane of polarization of multiple scattered light. 1: Homogeneous cloud layers
The degree of polarization and the direction of the plane of polarization are calculated by a Monte Carlo method for homogeneous layers. Two solar zenith angles and a range of optical thicknesses up to 10 are considered. The results are compared with calculations for single scattered photons. For a given pair of incident and scattered directions, there are only two possible values for the direction of the plane of polarization differing by 90 deg for single scattering from spherical aerosols. The choice between these two values depends only on the sign of the element M(-) in the first row and second column of the scattering matrix in the I, Q, U, V representation. In most cases there is little change in the direction of the plane of polarization when multiple scattering is taken into account, so that this quantity can usually be predicted from a very simple trigonometric relationship to good accuracy. Measurements of the direction of the plane of polarization at appropriately chosen angles provides information about the size distribution of the scattering centers.
Kattawar, G. W.↗
Lectures on the scattering of light
The exact (Mie) theory for the scattering of a plane wave by a dielectric sphere is presented. Since this infinite series solution is computationally impractical for large spheres, another formulation is given in terms of an integral equation valid for a bounded, but otherwise general array of scatterers. This equation is applied to the scattering by a single sphere, and several methods are suggested for approximating the scattering cross section in closed form. A tensor scattering matrix is introduced, in terms of which some general scattering theorems are derived. The application of the formalism to multiple scattering is briefly considered.
Saxon, D. S.↗