Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “random matrices”

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 73 records · Page 4

Design of optimal partial state feedback controllers for linear systems in stochastic environments

The problem of obtaining an optimal control law, which is constrained to be a feedback of the available measurements, is considered for both continuous and discrete time linear systems subjected to additive white process noise and measurement noise. Necessary conditions are obtained for minimizing a quadratic performance function for both finite and infinite duration cases. The feedback gain matrices are constrained to be constant for the infinite duration cases. For all the cases considered, algorithms are derived for generating sequences of feedback gain matrices which successively improve the performance function. Computational aspects are discussed via application to two continuous time processes, including a helicopter/slung load system subjected to measurement noise and random wind gust input.

Joshi, S. M.↗

A probabilistic model of a porous heat exchanger

This paper presents a probabilistic one-dimensional finite element model for heat transfer processes in porous heat exchangers. The Galerkin approach is used to develop the finite element matrices. Some of the submatrices are asymmetric due to the presence of the flow term. The Neumann expansion is used to write the temperature distribution as a series of random variables, and the expectation operator is applied to obtain the mean and deviation statistics. To demonstrate the feasibility of the formulation, a one-dimensional model of heat transfer phenomenon in superfluid flow through a porous media is considered. Results of this formulation agree well with the Monte-Carlo simulations and the analytical solutions. Although the numerical experiments are confined to parametric random variables, a formulation is presented to account for the random spatial variations.

Agrawal, O. P.↗

Sparse distributed memory and related models

Described here is sparse distributed memory (SDM) as a neural-net associative memory. It is characterized by two weight matrices and by a large internal dimension - the number of hidden units is much larger than the number of input or output units. The first matrix, A, is fixed and possibly random, and the second matrix, C, is modifiable. The SDM is compared and contrasted to (1) computer memory, (2) correlation-matrix memory, (3) feet-forward artificial neural network, (4) cortex of the cerebellum, (5) Marr and Albus models of the cerebellum, and (6) Albus' cerebellar model arithmetic computer (CMAC). Several variations of the basic SDM design are discussed: the selected-coordinate and hyperplane designs of Jaeckel, the pseudorandom associative neural memory of Hassoun, and SDM with real-valued input variables by Prager and Fallside. SDM research conducted mainly at the Research Institute for Advanced Computer Science (RIACS) in 1986-1991 is highlighted.

Kanerva, Pentti↗

Polarimetric radar data decomposition and interpretation

Significant efforts have been made to decompose polarimetric radar data into several simple scattering components. The components which are selected because of their physical significance can be used to classify SAR (Synthetic Aperture Radar) image data. If particular components can be related to forest parameters, inversion procedures may be developed to estimate these parameters from the scattering components. Several methods have been used to decompose an averaged Stoke's matrix or covariance matrix into three components representing odd (surface), even (double-bounce) and diffuse (volume) scatterings. With these decomposition techniques, phenomena, such as canopy-ground interactions, randomness of orientation, and size of scatters can be examined from SAR data. In this study we applied the method recently reported by van Zyl (1992) to decompose averaged backscattering covariance matrices extracted from JPL SAR images over forest stands in Maine, USA. These stands are mostly mixed stands of coniferous and deciduous trees. Biomass data have been derived from field measurements of DBH and tree density using allometric equations. The interpretation of the decompositions and relationships with measured stand biomass are presented in this paper.

Sun, Guoqing↗

Dual Extended Kalman Filter for the Identification of Time-Varying Human Manual Control Behavior

A Dual Extended Kalman Filter was implemented for the identification of time-varying human manual control behavior. Two filters that run concurrently were used, a state filter that estimates the equalization dynamics, and a parameter filter that estimates the neuromuscular parameters and time delay. Time-varying parameters were modeled as a random walk. The filter successfully estimated time-varying human control behavior in both simulated and experimental data. Simple guidelines are proposed for the tuning of the process and measurement covariance matrices and the initial parameter estimates. The tuning was performed on simulation data, and when applied on experimental data, only an increase in measurement process noise power was required in order for the filter to converge and estimate all parameters. A sensitivity analysis to initial parameter estimates showed that the filter is more sensitive to poor initial choices of neuromuscular parameters than equalization parameters, and bad choices for initial parameters can result in divergence, slow convergence, or parameter estimates that do not have a real physical interpretation. The promising results when applied to experimental data, together with its simple tuning and low dimension of the state-space, make the use of the Dual Extended Kalman Filter a viable option for identifying time-varying human control parameters in manual tracking tasks, which could be used in real-time human state monitoring and adaptive human-vehicle haptic interfaces.

manual control↗

Statistical analysis of effective singular values in matrix rank determination

A major problem in using SVD (singular-value decomposition) as a tool in determining the effective rank of a perturbed matrix is that of distinguishing between significantly small and significantly large singular values to the end, conference regions are derived for the perturbed singular values of matrices with noisy observation data. The analysis is based on the theories of perturbations of singular values and statistical significance test. Threshold bounds for perturbation due to finite-precision and i.i.d. random models are evaluated. In random models, the threshold bounds depend on the dimension of the matrix, the noisy variance, and predefined statistical level of significance. Results applied to the problem of determining the effective order of a linear autoregressive system from the approximate rank of a sample autocorrelation matrix are considered. Various numerical examples illustrating the usefulness of these bounds and comparisons to other previously known approaches are given.

Konstantinides, Konstantinos↗

Architectures and random properties of symplectic quantum circuits

Parametrized and random unitary (or orthogonal) n-qubit circuits play a central role in quantum information. As such, one could naturally assume that circuits implementing symplectic transformations would attract similar attention. However, this is not the case, as $\mathbb{SP}(d/2)$—the group of d × d unitary symplectic matrices—has thus far been overlooked. In this work, we aim at starting to fill this gap. We begin by presenting a universal set of generators $\mathcal{G}$ for the symplectic algebra $\mathfrak{sp}(d/2)$, consisting of one- and two-qubit Pauli operators acting on neighboring sites in a one-dimensional lattice. Here, we uncover two critical differences between such set, and equivalent ones for unitary and orthogonal circuits. Namely, we find that the operators in $\mathcal{G}$ cannot generate arbitrary local symplectic unitaries and that they are not translationally invariant. We then review the Schur–Weyl duality between the symplectic group and the Brauer algebra, and use tools from Weingarten calculus to prove that Pauli measurements at the output of Haar random symplectic circuits can converge to Gaussian processes. As a by-product, such analysis provides us with concentration bounds for Pauli measurements in circuits that form t-designs over $\mathbb{SP}(d/2)$. To finish, we present tensor-network tools to analyze shallow random symplectic circuits, and we use these to numerically show that computational-basis measurements anti-concentrate at logarithmic depth.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Get Non-Real: Randomized Sketching for High-Dimensional Non-Real Valued Data (Final Report)

In our final report for DE-C0022186, we describe the work we did on this grant towards the goals we proposed. Our first goal was characterizing fundamental limits for sketching of discrete high-dimensional matrices with low-dimensional structures. Our second main goal was designing algorithms for data reconstruction from sketches. We focus on approaches that are either specifically designed for non-real-valued data (binary, finite field) or that will translate more readily to that setting.

97 MATHEMATICS AND COMPUTING↗

Direct Demonstration of the Concept of Unrestricted Effective-Medium Approximation

The modified unrestricted effective-medium refractive index is defined as one that yields accurate values of a representative set of far-field scattering characteristics (including the scattering matrix) for an object made of randomly heterogeneous materials. We validate the concept of the modified unrestricted effective-medium refractive index by comparing numerically exact superposition T-matrix results for a spherical host randomly filled with a large number of identical small inclusions and Lorenz-Mie results for a homogeneous spherical counterpart. A remarkable quantitative agreement between the superposition T-matrix and Lorenz-Mie scattering matrices over the entire range of scattering angles demonstrates unequivocally that the modified unrestricted effective-medium refractive index is a sound (albeit still phenomenological) concept provided that the size parameter of the inclusions is sufficiently small and their number is sufficiently large. Furthermore, it appears that in cases when the concept of the modified unrestricted effective-medium refractive index works, its actual value is close to that predicted by the Maxwell-Garnett mixing rule.

Scattering↗

Estimation of full modal damping matrices from complex test modes

This paper describes the refinements of a previously published method for estimating a full modal damping matrix from complex test modes. It also documents application of the refined method to a structure where complex test modes were derived by the ERA method from multi-input random vibration test data. A numerical example based on simulated test data is presented to demonstrate the validity of the method. The application using real data was not successful, presumably because of noise in the small phase angles of the measured complex modes. Alternative test and data reduction procedures are suggested as possible remedies to the problem. A careful analysis of measurement and data processing errors should be made to examine basic feasibility before implementing the alternative procedures. The ability to estimate a full modal damping matrix is considered important for the preflight estimation of on-orbit damping, and for the synthesis of structural damping from substructure tests.

Hasselman, T. K.↗

Brazing SiC/SiC Composites to Metals

Experiments have shown that active brazing alloys (ABAs) can be used to join SiC/SiC composite materials to metals, with bond strengths sufficient for some structural applications. The SiC/SiC composite coupons used in the experiments were made from polymerbased SiC fiber preforms that were chemical-vapor-infiltrated with SiC to form SiC matrices. Some of the metal coupons used in the experiments were made from 304 stainless steel; others were made from oxygen-free, high-conductivity copper. Three ABAs were chosen for the experiments: two were chosen randomly from among a number of ABAs that were on hand at the time; the third ABA was chosen because its titanium content (1.25 percent) is less than those of the other two ABAs (1.75 and 4.5 percent, respectively) and it was desired to evaluate the effect of reducing the titanium content, as described below. The characteristics of ABAs that are considered to be beneficial for the purpose of joining SiC/SiC to metal include wettability, reactivity, and adhesion to SiC-based ceramics. Prior to further development, it was verified that the three chosen ABAs have these characteristics. For each ABA, suitable vacuum brazing process conditions were established empirically by producing a series of (SiC/SiC)/ABA wetting samples. These samples were then sectioned and subjected to scanning electron microscopy (SEM) and energy-dispersive x-ray spectrometry (EDS) for analysis of their microstructures and compositions. Specimens for destructive mechanical tests were fabricated by brazing of lap joints between SiC/SiC coupons 1/8-in. (.3.2- mm) thick and, variously, stainless steel or copper tabs. The results of destructive mechanical tests and the SEM/EDS analysis were used to guide the development of a viable method of brazing the affected materials.

Steffier, Wayne S.↗

Coherent Eigenmodes in Homogeneous MHD Turbulence

The statistical mechanics of Fourier models of ideal, homogeneous, incompressible magnetohydrodynamic (MHD) turbulence is discussed, along with their relevance for dissipative magnetofluids. Although statistical theory predicts that Fourier coefficients of fluid velocity and magnetic field are zero-mean random variables, numerical simulations clearly show that certain coefficients have a non-zero mean value that can be very large compared to the associated standard deviation, i.e., we have coherent structure. We use eigenanalysis of the modal covariance matrices in the probability density function to explain this phenomena in terms of `broken ergodicity', which is defined to occur when dynamical behavior does not match ensemble predictions on very long time-scales. We provide examples from 2-D and 3-D magnetohydrodynamic simulations of homogeneous turbulence, and show new results from long-time simulations of MHD turbulence with and without a mean magnetic field

Shebalin, John V.↗

Remote sensing of Earth Terrain

The Mueller matrix and polarization covariance matrix for polarimetric radar systems were studied. The clutter is modelled by a layer of random permittivity, described by a three dimensional correlation function, with variance, and horizontal and vertical correlation lengths. This model is applied, using the wave theory with Born approximations carried to the second order, to find the backscattering elements of the polarimetric matrices. It is found that 8 out of 16 elements of the Mueller matrix are identically zero, corresponding to a covariance matrix with four zero elements. Theoretical predictions are matched with experimental data for vegetation fields. A general mixing formula is derived for discrete scatterers immersed in a host medium. The inclusion particles are assumed to be ellipsoidal. The electric field inside the scatterers is determined by quasistatic analysis, assuming the diameter of the inclusion particles to be much smaller than the wavelength. The results are applicable to general multiphase mixtures.

Kong, J. A.↗

Radiative Transfer in a Discrete Random Medium Adjacent to a Half-Space with a Rough Interface

For a macroscopically plane-parallel discrete random medium, the boundary conditions for the specific coherency dyadic at a rough interface are derived. The derivation is based on a modification of the Twersky approximation for a scattering system consisting of a group of particles and the rough surface, and reduces to the solution of the scattering problem for a rough surface illuminated by a plane electromagnetic wave propagating in a discrete random medium with non-scattering boundaries. In a matrix-form setting, the boundary conditions for the specific coherency dyadic imply the boundary conditions for specific intensity column vectors which in turn, yield the expressions for the reflection and transmission matrices. The derived expressions are shown to be identical to those obtained by applying a phenomenological approach based on a facet model to the solution of the scattering problem for a rough surface illuminated by a plane electromagnetic wave.

Adrian Doicu↗

Linear discriminant analysis with misallocation in training samples

Linear discriminant analysis for a two-class case is studied in the presence of misallocation in training samples. A general appraoch to modeling of mislocation is formulated, and the mean vectors and covariance matrices of the mixture distributions are derived. The asymptotic distribution of the discriminant boundary is obtained and the asymptotic first two moments of the two types of error rate given. Certain numerical results for the error rates are presented by considering the random and two non-random misallocation models. It is shown that when the allocation procedure for training samples is objectively formulated, the effect of misallocation on the error rates of the Bayes linear discriminant rule can almost be eliminated. If, however, this is not possible, the use of Fisher rule may be preferred over the Bayes rule.

Chhikara, R.↗

Theoretical models for polarimetric radar clutter

The Mueller matrix and polarization covariance matrix are described for polarimetric radar systems. The clutter is modeled by a layer of random permittivity, described by a three-dimensional correlation function, with variance, and horizontal and vertical correlation lengths. This model is applied, using the wave theory with Born approximations carried to the second order, to find the backscattering elements of the polarimetric matrices. It is found that 8 out of 16 elements of the Mueller matrix are identically zero, corresponding to a covariance matrix with four zero elements. Theoretical predictions are matched with experimental data for vegetation fields.

Borgeaud, M.↗

Structural dynamics analysis using an unsymmetric block Lanczos algorithm

A method for reducing the order of a dynamical model of a large structure with arbitrary damping is developed analytically and demonstrated. A Lanczos algorithm is described which can reduce square unsymmetric system matrices to block-tridiagonal form, and a procedure for defining the reduced-order model from the right and left Lanczos vectors is outlined. Results for sample problems involving the 8-DOF FEM model of a beam-rotor assembly subjected to random and stepped external forces are presented in extensive graphs and briefly characterized.

Craig, Roy R., Jr.↗

NASA's New Orbital Debris Engineering Model, ORDEM2010

This paper describes the functionality and use of ORDEM2010, which replaces ORDEM2000, as the NASA Orbital Debris Program Office (ODPO) debris engineering model. Like its predecessor, ORDEM2010 serves the ODPO mission of providing spacecraft designers/operators and debris observers with a publicly available model to calculate orbital debris flux by current-state-of-knowledge methods. The key advance in ORDEM2010 is the input file structure of the yearly debris populations from 1995-2035 of sizes 10 micron - 1 m. These files include debris from low-Earth orbits (LEO) through geosynchronous orbits (GEO). Stable orbital elements (i.e., those that do not randomize on a sub-year timescale) are included in the files as are debris size, debris number, material density, random error and population error. Material density is implemented from ground-test data into the NASA breakup model and assigned to debris fragments accordingly. The random and population errors are due to machine error and uncertainties in debris sizes. These high-fidelity population files call for a much higher-level model analysis than what was possible with the populations of ORDEM2000. Population analysis in the ORDEM2010 model consists of mapping matrices that convert the debris population elements to debris fluxes. One output mode results in a spacecraft encompassing 3-D igloo of debris flux, compartmentalized by debris size, velocity, pitch, and yaw with respect to spacecraft ram direction. The second output mode provides debris flux through an Earth-based telescope/radar beam from LEO through GEO. This paper compares the new ORDEM2010 with ORDEM2000 in terms of processes and results with examples of specific orbits.

Krisko, Paula H.↗