SEARCH · Engineering Papers

Results for “Random Matrix”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS2001

Constraints on PSC Particle Microphysics Derived From Lidar Observations

Based on extensive T-matrix computations of light scattering by polydispersions of randomly oriented, rotationally symmetric nonspherical particles, we analyze existing lidar observations of polar stratospheric clouds (PSCs) and derive several constraints on PSC particle microphysical properties. We show that sharp-edged nonspherical particles (finite circular cylinders) exhibit less variability of lidar backscattering characteristics with particle size and aspect ratio than particles with smooth surfaces (spheroids). For PSC particles significantly smaller than the wavelength, the backscatter color index Alpha and the depolarization color index Beta are essentially shape-independent. Observations for type Ia PSCs can be reproduced by spheroids with aspect ratios larger than 1.2, oblate cylinders with diameter-to-length ratios greater than 1.6, and prolate cylinders with length-to-diameter ratios greater than 1.4. The effective equal-volume-sphere radius for type la PSCs is about 0.8 microns or larger. Type Ib PSCs are likely to be composed of spheres or nearly spherical particles with effective radii smaller than 0.8 microns. Observations for type II PSCs are consistent with large ice crystals (effective radius greater than 1 micron modeled as cylinders or prolate spheroids.

Liu, Li↗

NASA NTRS2015

Polarized Bidirectional Reflectance of Optically Thick Sparse Particulate Layers: an Efficient Numerically Exact Radiative-Transfer Solution

We describe a simple yet efficient numerical algorithm for computing polarized bidirectional reflectance of an optically thick (semi-infinite), macroscopically flat layer composed of statistically isotropic and mirror symmetric random particles. The spatial distribution of the particles is assumed to be sparse, random, and statistically uniform. The 44 Stokes reflection matrix is calculated by iterating the Ambartsumian's vector nonlinear integral equation. The result is a numerically exact solution of the vector radiative transfer equation and as such fully satisfies the energy conservation law and the fundamental reciprocity relation. Since this technique bypasses the computation of the internal radiation field, it is very fast and highly accurate. The FORTRAN implementation of the technique is publicly available on the World Wide Web at http://www.giss.nasa.gov/staff/ mmishchenko/brf. It can be combined with several existing computer programs providing the requisite single-scattering properties of spherical or morphologically complex particles and applied to a wide range of optical characterization problems. Benchmark results obtained with this program can be used for testing alternative solvers of the vector radiative transfer equation.

radiative transfer↗

NASA NTRS1990

Scattering from randomly oriented scatterers with strong permittivity fluctuations

Strong permittivity fluctuation theory is used to solve the problem of scattering from a medium composed of completely randomly oriented scatterers under a low frequency limit. Gaussian statistics are not assumed for the renormalized scattering sources. Numerical results on effective permittivity are illustrated for oblate and prolate spheroidal scatterers and compared with the results for spherical scatterers. The results are consistent with discrete scatterer theory. The effective permittivity of a random medium embedded with nonspherical scatterers shows a higher imaginary part than the spherical scatterer case with equal correlation volume. Under the distorted Born approximation, the polarimetric covariance matrix for the backscattered electric field is calculated for half-space randomly oriented scatterers. The nonspherical geometry of the scatterers shows significant effects on the cross-polarized backscattering returns, and the correlation coefficient between HH and VV returns. The polarimetric backscattering coefficients can provide useful information in distinguishing the geometry of scatterers.

Yueh, S. H.↗

NASA NTRS1988

Random field estimation approach to multibody dynamics

The use of recursive techniques similar to random field models to factor and invert the inertia matrix of a multibody system is discussed. An equivalence is established between the composite multibody system inertia matrix and the covariance of the output of a described linear system model. Conditional mean estimation and sequential estimation problems are solved along with problems of filtering and smoothing. Formulas are developed to compute the covariance of several relevant quantities. The foregoing results are used to obtain the inverse of the composite multibody system inertia in closed form.

Rodriguez, G.↗

NASA NTRS1992

Comparison of three controllers applied to helicopter vibration

A comparison was made of the applicability and suitability of the deterministic controller, the cautious controller, and the dual controller for the reduction of helicopter vibration by using higher harmonic blade pitch control. A randomly generated linear plant model was assumed and the performance index was defined to be a quadratic output metric of this linear plant. A computer code, designed to check out and evaluate these controllers, was implemented and used to accomplish this comparison. The effects of random measurement noise, the initial estimate of the plant matrix, and the plant matrix propagation rate were determined for each of the controllers. With few exceptions, the deterministic controller yielded the greatest vibration reduction (as characterized by the quadratic output metric) and operated with the greatest reliability. Theoretical limitations of these controllers were defined and appropriate candidate alternative methods, including one method particularly suitable to the cockpit, were identified.

Leyland, Jane A.↗

NASA NTRS1992

Redistribution of Ti and Al in deuterium charged TiAl

The redistribution of titanium and aluminum in a single-phase TiAl alloy that has been exposed to a high pressure of deuterium gas at high temperature is studied. The microstructure in the as-received, uncharged condition consisted of single-phase gamma TiAl grains and a random distribution of precipitates. Precipitates were distributed throughout the matrix and along the grain boundaries. The chemistry of the precipitates varied considerably; some were rich in Al, while other were mostly Ti with some Si and Al. The dislocation density in most grains was low, although in a few grains a high dislocation density was observed. FCC deuterides with a lattice parameter of 0.45 nm form in a Ti-52.1Al-2.1Ta (at. pct) alloy after exposure to 1.38 MPa of deuterium gas at 650 C for 213 hr. The structure and lattice parameter are consistent with the formation of Ti(l)D2. The deuterides that form in this alloy are enriched in Ti and deficient in Al and Ta compared to the deuteride-free matrix. Regions of the matrix contiguous with the deuterides have a correspondingly enhanced aluminum and tantalum concentration.

Legzdina, D.↗

NASA NTRS2005

Comparison of the Lithologies on the Surface of the Asteroid 4 Vesta Based on the Petrology of 91 and 92 Series Antarctic Achondrites

We studied meteorite sections PCA91179,9; EET92003,14; EET92025,5; EET92026,4; and EET92027,5 t determine lithological variations on the surface of the asteroid 4 Vesta. These are brecciated polymict eucrites. The clasts contained within these samples and the matrix which surrounds them in each sample provide a random sampling of the characteristic compositions found on the surface of asteroid 4 Vesta.

Source record↗

NASA NTRS2011

A Multiple Sphere T-Matrix Fortran Code for Use on Parallel Computer Clusters

A general-purpose Fortran-90 code for calculation of the electromagnetic scattering and absorption properties of multiple sphere clusters is described. The code can calculate the efficiency factors and scattering matrix elements of the cluster for either fixed or random orientation with respect to the incident beam and for plane wave or localized- approximation Gaussian incident fields. In addition, the code can calculate maps of the electric field both interior and exterior to the spheres.The code is written with message passing interface instructions to enable the use on distributed memory compute clusters, and for such platforms the code can make feasible the calculation of absorption, scattering, and general EM characteristics of systems containing several thousand spheres.

Mackowski, D. W.↗

NASA NTRS1992

Load redistribution considerations in the fracture of ceramic matrix composites

Using a macroscopic viewpoint, composite laminae are homogeneous orthotropic solids whose directional strengths are random variables. Incorporation of these random variable strengths into failure models, either interactive or noninteractive, allows for the evaluation of the lamina reliability under a given stress state. Using a noninteractive criterion for demonstration purposes, laminate reliabilities are calculated assuming previously established load sharing rules for the redistribution of load as the failure of laminae occur. The matrix cracking predicted by ACK theory is modeled to allow a loss of stiffness in the fiber direction. The subsequent failure in the fiber direction is controlled by a modified bundle theory. Results are compared with previous models which did not permit separate consideration of matrix cracking, as well as to results obtained from experimental data. The effects of variations from the ideal physical geometry which is normally used to depict the matrix cracking are also studied.

Thomas, David J.↗

NASA NTRS1988

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↗

NASA NTRS1973

On Digital Simulation of Multicorrelated Random Processes and Its Applications

Two methods are described to simulate, on a digital computer, a set of correlated, stationary, and Gaussian time series with zero mean from the given matrix of power spectral densities and cross spectral densities. The first method is based upon trigonometric series with random amplitudes and deterministic phase angles. The random amplitudes are generated by using a standard random number generator subroutine. An example is given which corresponds to three components of wind velocities at two different spatial locations for a total of six correlated time series. In the second method, the whole process is carried out using the Fast Fourier Transform approach. This method gives more accurate results and works about twenty times faster for a set of six correlated time series.

Sinha, A. K.↗

NASA NTRS1974

On the design of optimal input signals in system identification

The problem of designing optimal inputs in the identification of linear systems with unknown random parameters is considered using a Bayesian approach. The information matrix, which is positive definite for the class of systems analyzed, gives a measure of performance for the system inputs. The computation of the optimal closed-loop input mappings is shown to be a nontrivial exercise in adaptive control. Deterministic optimal inputs are shown to be easily computable. Numerical examples are given. A Kalman filter is used to estimate the parameters. A necessary condition for the Kalman filter not to diverge when applying linear feedback is also given.

Lopez-Toledo, A. A.↗

NASA NTRS1986

Vector solution for the mean electromagnetic fields in a layer of random particles

The mean electromagnetic fields are found in a layer of randomly oriented particles lying over a half space. A matrix-dyadic formulation of Maxwell's equations is employed in conjunction with the Foldy-Lax approximation to obtain equations for the mean fields. A two variable perturbation procedure, valid in the limit of small fractional volume, is then used to derive uncoupled equations for the slowly varying amplitudes of the mean wave. These equations are solved to obtain explicit expressions for the mean electromagnetic fields in the slab region in the general case of arbitrarily oriented particles and arbitrary polarization of the incident radiation. Numerical examples are given for the application to remote sensing of vegetation.

Lang, R. H.↗

NASA NTRS2018

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↗

NASA NTRS2005

Probabilisitc Geobiological Classification Using Elemental Abundance Distributions and Lossless Image Compression in Recent and Modern Organisms

Last year we presented techniques for the detection of fossils during robotic missions to Mars using both structural and chemical signatures[Storrie-Lombardi and Hoover, 2004]. Analyses included lossless compression of photographic images to estimate the relative complexity of a putative fossil compared to the rock matrix [Corsetti and Storrie-Lombardi, 2003] and elemental abundance distributions to provide mineralogical classification of the rock matrix [Storrie-Lombardi and Fisk, 2004]. We presented a classification strategy employing two exploratory classification algorithms (Principal Component Analysis and Hierarchical Cluster Analysis) and non-linear stochastic neural network to produce a Bayesian estimate of classification accuracy. We now present an extension of our previous experiments exploring putative fossil forms morphologically resembling cyanobacteria discovered in the Orgueil meteorite. Elemental abundances (C6, N7, O8, Na11, Mg12, Ai13, Si14, P15, S16, Cl17, K19, Ca20, Fe26) obtained for both extant cyanobacteria and fossil trilobites produce signatures readily distinguishing them from meteorite targets. When compared to elemental abundance signatures for extant cyanobacteria Orgueil structures exhibit decreased abundances for C6, N7, Na11, All3, P15, Cl17, K19, Ca20 and increases in Mg12, S16, Fe26. Diatoms and silicified portions of cyanobacterial sheaths exhibiting high levels of silicon and correspondingly low levels of carbon cluster more closely with terrestrial fossils than with extant cyanobacteria. Compression indices verify that variations in random and redundant textural patterns between perceived forms and the background matrix contribute significantly to morphological visual identification. The results provide a quantitative probabilistic methodology for discriminating putatitive fossils from the surrounding rock matrix and &om extant organisms using both structural and chemical information. The techniques described appear applicable to the geobiological analysis of meteoritic samples or in situ exploration of the Mars regolith. Keywords: cyanobacteria, microfossils, Mars, elemental abundances, complexity analysis, multifactor analysis, principal component analysis, hierarchical cluster analysis, artificial neural networks, paleo-biosignatures

Storrie-Lombardi, Michael C.↗

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 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 NTRS2011

Direct Simulation of Multiple Scattering by Discrete Random Media Illuminated by Gaussian Beams

The conventional orientation-averaging procedure developed in the framework of the superposition T-matrix approach is generalized to include the case of illumination by a Gaussian beam (GB). The resulting computer code is parallelized and used to perform extensive numerically exact calculations of electromagnetic scattering by volumes of discrete random medium consisting of monodisperse spherical particles. The size parameters of the scattering volumes are 40, 50, and 60, while their packing density is fixed at 5%. We demonstrate that all scattering patterns observed in the far-field zone of a random multisphere target and their evolution with decreasing width of the incident GB can be interpreted in terms of idealized theoretical concepts such as forward-scattering interference, coherent backscattering (CB), and diffuse multiple scattering. It is shown that the increasing violation of electromagnetic reciprocity with decreasing GB width suppresses and eventually eradicates all observable manifestations of CB. This result supplements the previous demonstration of the effects of broken reciprocity in the case of magneto-optically active particles subjected to an external magnetic field.

Mackowski, Daniel W.↗