Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Random matrix theory”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Radiative Transfer Theory Verified by Controlled Laboratory Experiments

We report the results of high-accuracy controlled laboratory measurements of the Stokes reflection matrix for suspensions of submicrometer-sized latex particles in water and compare them with the results of a numerically exact computer solution of the vector radiative transfer equation (VRTE). The quantitative performance of the VRTE is monitored by increasing the volume packing density of the latex particles from 2 to 10. Our results indicate that the VRTE can be applied safely to random particulate media with packing densities up to 2. VRTE results for packing densities of the order of 5 should be taken with caution, whereas the polarized bidirectional reflectivity of suspensions with larger packing densities cannot be accurately predicted. We demonstrate that a simple modification of the phase matrix entering the VRTE based on the so-called static structure factor can be a promising remedy that deserves further examination.

Radiative transfer↗
NASA NTRS2017

Micromechanics of Composite Materials Governed by Vector Constitutive Laws

The high-fidelity generalized method of cells micromechanics theory has been extended for the prediction of the effective property tensor and the corresponding local field distributions for composites whose constituents are governed by vector constitutive laws. As shown, the shear analogy, which can predict effective transverse properties, is not valid in the general three-dimensional case. Consequently, a general derivation is presented that is applicable to both continuously and discontinuously reinforced composites with arbitrary vector constitutive laws and periodic microstructures. Results are given for thermal and electric problems, effective properties and local field distributions, ordered and random microstructures, as well as complex geometries including woven composites. Comparisons of the theory's predictions are made to test data, numerical analysis, and classical expressions from the literature. Further, classical methods cannot provide the local field distributions in the composite, and it is demonstrated that, as the percolation threshold is approached, their predictions are increasingly unreliable. XXXX It has been observed that the bonding between the fibers and matrix in composite materials can be imperfect. In the context of thermal conductivity, such imperfect interfaces have been investigated in micromechanical models by Dunn and Taya (1993), Duan and Karihaloo (2007), Nan et al. (1997) and Hashin (2001). The present HFGMC micromechanical method, derived for perfectly bonded composite materials governed by vector constitutive laws, can be easily generalized to include the effects of weak bonding between the constituents. Such generalizations, in the context of the mechanical micromechanics problem, involve introduction of a traction-separation law at the fiber/matrix interface and have been presented by Aboudi (1987), Bednarcyk and Arnold (2002), Bednarcyk et al. (2004) and Aboudi et al. (2013) and will be addressed in the future.

micromechanics↗
NASA NTRS1987

Spatially random models, estimation theory, and robot arm dynamics

Spatially random models provide an alternative to the more traditional deterministic models used to describe robot arm dynamics. These alternative models can be used to establish a relationship between the methodologies of estimation theory and robot dynamics. A new class of algorithms for many of the fundamental robotics problems of inverse and forward dynamics, inverse kinematics, etc. can be developed that use computations typical in estimation theory. The algorithms make extensive use of the difference equations of Kalman filtering and Bryson-Frazier smoothing to conduct spatial recursions. The spatially random models are very easy to describe and are based on the assumption that all of the inertial (D'Alembert) forces in the system are represented by a spatially distributed white-noise model. The models can also be used to generate numerically the composite multibody system inertia matrix. This is done without resorting to the more common methods of deterministic modeling involving Lagrangian dynamics, Newton-Euler equations, etc. These methods make substantial use of human knowledge in derivation and minipulation of equations of motion for complex mechanical systems.

Rodriguez, G.↗
DOE OSTI2025

Dual-unitary shadow tomography

We introduce a classical shadow tomography scheme based on dual-unitary brick-wall circuits termed "dual-unitary shadow tomography" (DUST). For this we study operator spreading and Pauli weight dynamics in one-dimensional qubit systems, evolved by random two-local dual-unitary gates arranged in a brick-wall structure, ending with a final measurement layer. We do this by deriving general constraints on the Pauli weight transfer matrix and specializing to the case of dual-unitarity. We first show that dual-unitaries must have a minimal amount of entropy production. Remarkably, we find that operator spreading in these circuits have a rich structure resembling that of relativistic quantum field theories, with massless chiral excitations that can decay or fuse into each other, which we call left- or right-movers. We develop a mean-field description of the Pauli weight in terms of $\rho(x,t)$, which represents the probability of having nontrivial support at site $x$ and depth $t$ starting from a fixed weight distribution. We develop an equation of state for $\rho(x,t)$, and simulate it numerically using Monte Carlo simulations. Lastly, we demonstrate that the fast-thermalizing properties of dual-unitary circuits make them better at predicting large operators than shallow brick-wall Clifford circuits. Our results are robust to finite-size effects due to the chirality of dual-unitary brick-wall circuits.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2013

Coarsening in Solid-liquid Mixtures: Overview of Experiments on Shuttle and ISS

The microgravity environment on the Shuttle and the International Space Station (ISS) provides the ideal condition to perform experiments on Coarsening in Solid-Liquid Mixtures (CSLM) as deleterious effects such as particle sedimentation and buoyancy-induced convection are suppressed. For an ideal system such as Lead-Tin in which all the thermophysical properties are known, the initial condition in microgravity of randomly dispersed particles with local clustering of solid Tin in eutectic liquid Lead-Tin matrix, permitted kinetic studies of competitive particle growth for a range of volume fractions. Verification that the quenching phase of the experiment had negligible effect of the spatial distribution of particles is shown through the computational solution of the dynamical equations of motion, thus insuring quench-free effects from the coarsened microstructure measurements. The low volume fraction experiments conducted on the Shuttle showed agreement with transient Ostwald ripening theory, and the steady-state requirement of LSW theory was not achieved. More recent experiments conducted on ISS with higher volume fractions have achieved steady-state condition and show that the kinetics follows the classical diffusion limited particle coarsening prediction and the measured 3D particle size distribution becomes broader as predicted from theory.

Ostwald ripening↗
NASA NTRS1993

Probabilistic sizing of laminates with uncertainties

A reliability based design methodology for laminate sizing and configuration for a special case of composite structures is described. The methodology combines probabilistic composite mechanics with probabilistic structural analysis. The uncertainties of constituent materials (fiber and matrix) to predict macroscopic behavior are simulated using probabilistic theory. Uncertainties in the degradation of composite material properties are included in this design methodology. A multi-factor interaction equation is used to evaluate load and environment dependent degradation of the composite material properties at the micromechanics level. The methodology is integrated into a computer code IPACS (Integrated Probabilistic Assessment of Composite Structures). Versatility of this design approach is demonstrated by performing a multi-level probabilistic analysis to size the laminates for design structural reliability of random type structures. The results show that laminate configurations can be selected to improve the structural reliability from three failures in 1000, to no failures in one million. Results also show that the laminates with the highest reliability are the least sensitive to the loading conditions.

Shah, A. R.↗
NASA NTRS1993

Interactive Reliability Model for Whisker-toughened Ceramics

Wider use of ceramic matrix composites (CMC) will require the development of advanced structural analysis technologies. The use of an interactive model to predict the time-independent reliability of a component subjected to multiaxial loads is discussed. The deterministic, three-parameter Willam-Warnke failure criterion serves as the theoretical basis for the reliability model. The strength parameters defining the model are assumed to be random variables, thereby transforming the deterministic failure criterion into a probabilistic criterion. The ability of the model to account for multiaxial stress states with the same unified theory is an improvement over existing models. The new model was coupled with a public-domain finite element program through an integrated design program. This allows a design engineer to predict the probability of failure of a component. A simple structural problem is analyzed using the new model, and the results are compared to existing models.

Palko, Joseph L.↗
NASA NTRS2013

Comprehensive T-Matrix Reference Database: A 2012 - 2013 Update

The T-matrix method is one of the most versatile, efficient, and accurate theoretical techniques widely used for numerically exact computer calculations of electromagnetic scattering by single and composite particles, discrete random media, and particles imbedded in complex environments. This paper presents the fifth update to the comprehensive database of peer-reviewed T-matrix publications initiated by us in 2004 and includes relevant publications that have appeared since 2012. It also lists several earlier publications not incorporated in the original database, including Peter Waterman's reports from the 1960s illustrating the history of the T-matrix approach and demonstrating that John Fikioris and Peter Waterman were the true pioneers of the multi-sphere method otherwise known as the generalized Lorenz - Mie theory.

Complex scattering↗
NASA NTRS1976

Optimal control of a low wing-loading STOL aircraft

Linear optimal quadratic control theory is applied to a low wing-loading STOL aircraft for ride quality and flight path following. Design criteria include minimum rms response to wind turbulence and desired transient response characteristics. Design techniques include proper choosing of design versus evaluation models, choosing appropriate performance index responses, and use of classical evaluation techniques. Results are obtained through a combination of frequency response shaping and gust observation. Effects of control rate and authority saturation are examined with a new rapid calculation of random input describing functions. Parameter sensitivity is also evaluated using a Liapunov type matrix equation.

Cunningham, T. B.↗
NASA NTRS1996

Kinetics of Diffusional Droplet Growth in a Liquid/Liquid Two-Phase System

We address the problem of diffusional interactions in a finite sized cluster of spherical particles for volume fractions, V(sub v) in the range 0-0.01. We determined the quasi-static monopole diffusion solution for n particles distributed at random in a continuous matrix. A global mass conservation condition is employed, obviating the need for any external boundary condition. The numerical results provide the instantaneous (snapshot) growth or shrinkage rate of each particle, precluding the need for extensive time-dependent computations. The close connection between these snapshot results and the coarsegrained kinetic constants are discussed. A square-root dependence of the deviations of the rate constants from their zero volume fraction value is found for the higher V(sub v) investigated. This behavior is consistent with predictions from diffusion Debye-Huckel screening theory. By contrast, a cube-root dependence, reported in earlier numerical studies, is found for the lower V(sub v) investigated. The roll-over region of the volume fraction where the two asymptotics merge depends on the number of particles, n, alone. A theoretical estimate for the roll-over point predicts that the corresponding V(sub v) varies as n(sup -2), in good agreement with the numerical results.

Glicksman, M. E.↗
NASA NTRS2017

Modeling of Electromagnetic Scattering by Discrete and Discretely Heterogeneous Random Media by Using Numerically Exact Solutions of the Maxwell Equations

In this paper, we discuss some aspects of numerical modeling of electromagnetic scattering by discrete random medium by using numerically exact solutions of the macroscopic Maxwell equations. Typical examples of such media are clouds of interstellar dust, clouds of interplanetary dust in the Solar system, dusty atmospheres of comets, particulate planetary rings, clouds in planetary atmospheres, aerosol particles with numerous inclusions and so on. Our study is based on the results of extensive computations of different characteristics of electromagnetic scattering obtained by using the superposition T-matrix method which represents a direct computer solver of the macroscopic Maxwell equations for an arbitrary multisphere configuration. As a result, in particular, we clarify the range of applicability of the low-density theories of radiative transfer and coherent backscattering as well as of widely used effective-medium approximations.

Dlugach, Janna M.↗
NASA NTRS2015

A Numerical Solution Routine for Investigating Oxidation-Induced Strength Degradation Mechanisms in SiC/SiC Composites

The stress rupture strength of silicon carbide fiber-reinforced silicon carbide (SiCSiC) composites with a boron nitride (BN) fiber coating decreases with time within the intermediate temperature range of 700-950 C. Various theories have been proposed to explain the cause of the time dependent stress rupture strength. Some previous authors have suggested that the observed composite strength behavior is due to the inherent time dependent strength of the fibers, which is caused by the slow growth of flaws within the fibers. Flaw growth is supposedly enabled by oxidation of free carbon at the grain boundaries. The objective of this paper is to investigate the relative significance of the various theories for the time-dependent strength of SiCSiC composites. This is achieved through the development of a numerically-based progressive failure analysis routine and through the application of the routine to simulate the composite stress rupture tests. The progressive failure routine is a time marching routine with an iterative loop between a probability of fiber survival equation and a force equilibrium equation within each time step. Failure of the composite is assumed to initiate near a matrix crack and the progression of fiber failures occurs by global load sharing. The probability of survival equation is derived from consideration of the strength of ceramic fibers with randomly occurring and slow growing flaws as well as the mechanical interaction between the fibers and matrix near a matrix crack. The force equilibrium equation follows from the global load sharing presumption. The results of progressive failure analyses of the composite tests suggest that the relationship between time and stress-rupture strength is attributed almost entirely to the slow flaw growth within the fibers. Although other mechanisms may be present, they appear to have only a minor influence on the observed time dependent behavior.

degradation↗
NASA NTRS1987

The capacity of the Hopfield associative memory

Techniques from coding theory are applied to study rigorously the capacity of the Hopfield associative memory. Such a memory stores n-tuple of + or - 1s. The components change depending on a hard-limited version of linear functions of all other components. With symmetric connections between components, a stable state is ultimately reached. By building up the connection matrix as a sum-of-outer products of m fundamental memories, it may be possible to recover a certain one of the m memories by using an initial n-tuple probe vector less than a Hamming distance n/2 away from the fundamental memory. If m fundamental memories are chosen at random, the maximum asymptotic value of m in order that most of the m original memories are exactly recoverable is n/(2 log n). With the added restriction that every one of the m fundamental memories be recoverable exactly, m can be no more than n/(4 log n) asymptotically as n approaches infinity. Extensions are also considered, in particular to capacity under quantization of the outer-product connection matrix. This quantized memory-capacity problem is closely related to the capacity of the quantized Gaussian channel.

Mceliece, Robert J.↗
NASA NTRS2016

Time-Dependent Stress Rupture Strength Degradation of Hi-Nicalon Fiber-Reinforced Silicon Carbide Composites at Intermediate Temperatures

The stress rupture strength of silicon carbide fiber-reinforced silicon carbide composites with a boron nitride fiber coating decreases with time within the intermediate temperature range of 700 to 950 degree Celsius. Various theories have been proposed to explain the cause of the time-dependent stress rupture strength. The objective of this paper is to investigate the relative significance of the various theories for the time-dependent strength of silicon carbide fiber-reinforced silicon carbide composites. This is achieved through the development of a numerically based progressive failure analysis routine and through the application of the routine to simulate the composite stress rupture tests. The progressive failure routine is a time-marching routine with an iterative loop between a probability of fiber survival equation and a force equilibrium equation within each time step. Failure of the composite is assumed to initiate near a matrix crack and the progression of fiber failures occurs by global load sharing. The probability of survival equation is derived from consideration of the strength of ceramic fibers with randomly occurring and slow growing flaws as well as the mechanical interaction between the fibers and matrix near a matrix crack. The force equilibrium equation follows from the global load sharing presumption. The results of progressive failure analyses of the composite tests suggest that the relationship between time and stress-rupture strength is attributed almost entirely to the slow flaw growth within the fibers. Although other mechanisms may be present, they appear to have only a minor influence on the observed time-dependent behavior.

Progressive Failure Analysis↗
NASA NTRS1997

Visibility of wavelet quantization noise

The discrete wavelet transform (DWT) decomposes an image into bands that vary in spatial frequency and orientation. It is widely used for image compression. Measures of the visibility of DWT quantization errors are required to achieve optimal compression. Uniform quantization of a single band of coefficients results in an artifact that we call DWT uniform quantization noise; it is the sum of a lattice of random amplitude basis functions of the corresponding DWT synthesis filter. We measured visual detection thresholds for samples of DWT uniform quantization noise in Y, Cb, and Cr color channels. The spatial frequency of a wavelet is r 2-lambda, where r is display visual resolution in pixels/degree, and lambda is the wavelet level. Thresholds increase rapidly with wavelet spatial frequency. Thresholds also increase from Y to Cr to Cb, and with orientation from lowpass to horizontal/vertical to diagonal. We construct a mathematical model for DWT noise detection thresholds that is a function of level, orientation, and display visual resolution. This allows calculation of a "perceptually lossless" quantization matrix for which all errors are in theory below the visual threshold. The model may also be used as the basis for adaptive quantization schemes.

NASA Center ARC↗
NASA NTRS2016

Scattering of Gaussian Beams by Disordered Particulate Media

A frequently observed characteristic of electromagnetic scattering by a disordered particulate medium is the absence of pronounced speckles in angular patterns of the scattered light. It is known that such diffuse speckle-free scattering patterns can be caused by averaging over randomly changing particle positions and/or over a finite spectral range. To get further insight into the possible physical causes of the absence of speckles, we use the numerically exact superposition T-matrix solver of the Maxwell equations and analyze the scattering of plane-wave and Gaussian beams by representative multi-sphere groups. We show that phase and amplitude variations across an incident Gaussian beam do not serve to extinguish the pronounced speckle pattern typical of plane-wave illumination of a fixed multi-particle group. Averaging over random particle positions and/or over a finite spectral range is still required to generate the classical diffuse speckle-free regime.

radiative transfer theory↗
NASA NTRS1997

Modulation of statolith mass and grouping in white clover (Trifolium repens) growth in 1-g, microgravity and on the clinostat

Current models of gravity perception in higher plants focus on the buoyant weight of starch-filled amyloplasts as the initial gravity signal susceptor (statolith). However, no tests have yet determined if statolith mass is regulated to increase or decrease gravity stimulus to the plant. To this end, the root caps of white clover (Trifolium repens) grown in three gravity environments with three different levels of gravity stimulation have been examined: (i) 1-g control with normal static gravistimulation, (ii) on a slow clinostat with constant gravistimulation, and (iii) in the stimulus-free microgravity aboard the Space Shuttle. Seedlings were germinated and grown in the BioServe Fluid Processing Apparatus and root cap structure was examined at both light and electron microscopic levels, including three-dimensional cell reconstruction from serial sections. Quantitative analysis of the electron micrographs demonstrated that the starch content of amyloplasts varied with seedling age but not gravity condition. It was also discovered that, unlike in starch storage amyloplasts, all of the starch granules of statolith amyloplasts were encompassed by a fine filamentous, ribosome-excluding matrix. From light micrographic 3-D cell reconstructions, the absolute volume, number, and positional relationships between amyloplasts showed (i) that individual amyloplast volume increased in microgravity but remained constant in seedlings grown for up to three days on the clinostat, (ii) the number of amyloplasts per cell remained unchanged in microgravity but decreased on the clinostat, and (iii) the three-dimensional positions of amyloplasts were not random. Instead amyloplasts in microgravity were grouped near the cell centers while those from the clinostat appeared more dispersed. Taken together, these observations suggest that changing gravity stimulation can elicit feedback control over statolith mass by changing the size, number, and grouping of amyloplasts. These results support the starch-statolith theory of graviperception in higher plants and add to current models with a new feedback control loop as a mechanism for modulation of statolith responsiveness to inertial acceleration.

Flight Experiment↗
NASA NTRS1996

Perceptually Lossless Wavelet Compression

The Discrete Wavelet Transform (DWT) decomposes an image into bands that vary in spatial frequency and orientation. It is widely used for image compression. Measures of the visibility of DWT quantization errors are required to achieve optimal compression. Uniform quantization of a single band of coefficients results in an artifact that is the sum of a lattice of random amplitude basis functions of the corresponding DWT synthesis filter, which we call DWT uniform quantization noise. We measured visual detection thresholds for samples of DWT uniform quantization noise in Y, Cb, and Cr color channels. The spatial frequency of a wavelet is r 2(exp -1), where r is display visual resolution in pixels/degree, and L is the wavelet level. Amplitude thresholds increase rapidly with spatial frequency. Thresholds also increase from Y to Cr to Cb, and with orientation from low-pass to horizontal/vertical to diagonal. We propose a mathematical model for DWT noise detection thresholds that is a function of level, orientation, and display visual resolution. This allows calculation of a 'perceptually lossless' quantization matrix for which all errors are in theory below the visual threshold. The model may also be used as the basis for adaptive quantization schemes.

Watson, Andrew B.↗