Search NASA⌕ Search

SEARCH · Search NASA

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.

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

At least 217 records · Page 12

Thermal expansion behavior of graphite/glass and graphite/magnesium

The thermal expansion behavior of n (+/- 8)s graphite fiber reinforced magnesium laminate and four graphite reinforced glass-matrix laminates (a unidirectional laminate, a quasi-isotropic laminate, a symmetric low angle-ply laminate, and a random chopped-fiber mat laminate) was determined, and was found, in all cases, to not be significantly affected by thermal cycling. Specimens were cycled up to 100 times between -200 F and 100 F, and the thermal expansion coefficients determined for each material as a function of temperature were found to be low. Some dimensional changes as a function of thermal cycling, and some thermal-strain hysteresis, were observed.

Tompkins, Stephen S.↗

Optics of Water Microdroplets with Soot Inclusions: Exact Versus Approximate Results

We use the recently generalized version of the multi-sphere superposition T-matrix method (STMM) to compute the scattering and absorption properties of microscopic water droplets contaminated by black carbon. The soot material is assumed to be randomly distributed throughout the droplet interior in the form of numerous small spherical inclusions. Our numerically-exact STMM results are compared with approximate ones obtained using the Maxwell-Garnett effective-medium approximation (MGA) and the Monte Carlo ray-tracing approximation (MCRTA). We show that the popular MGA can be used to calculate the droplet optical cross sections, single-scattering albedo, and asymmetry parameter provided that the soot inclusions are quasi-uniformly distributed throughout the droplet interior, but can fail in computations of the elements of the scattering matrix depending on the volume fraction of soot inclusions. The integral radiative characteristics computed with the MCRTA can deviate more significantly from their exact STMM counterparts, while accurate MCRTA computations of the phase function require droplet size parameters substantially exceeding 60.

cloud droplets↗

Comprehensive Thematic T-Matrix Reference Database: a 2017-2019 Update

Since its inception in the mid-1960s, the T-matrix method has been one of the most versatile and efficient numerically exact computer solvers of the time-harmonic macroscopic Maxwell equations. It has been widely used for the computation of electromagnetic scattering by single and composite particles, discrete random media, periodic structures (including metamaterials), and particles in the vicinity of plane or rough interfaces separating media with different refractive indices. This compilation is the ninth update to the comprehensive thematic database of peer-reviewed T-matrix publications initiated in 2004 and lists relevant publications that have appeared since 2017. It also includes a few earlier publications that have so far been overlooked.

Michael I Mishchenko↗

Radiative transfer theory for active remote sensing of a layer of nonspherical particles

The radiative transfer theory is applied to calculate the scattering by a layer of randomly positioned and oriented nonspherical particles. The scattering amplitude functions of each individual particle are calculated with Waterman's T matrix method, which utilizes vector spherical wave functions for expansion of incident, scattered, and surface fields. The orientation of the particles is described by a probability density function of the Eulerian angles of rotation. A rotation matrix is used to relate the T matrix of the principal frame to that of the natural frame of the particle. The extinction matrix and phase matrix of the radiative transfer equations are expressed in terms of the T matrix elements. The extinction matrix for nonspherical particles is generally nondiagonal. There are only two attenuation rates in a specified direction of propagation. The radiative transfer equations are solved by an iterative method to first order in albedo. Numerical results are illustrated as functions of incidence angle and frequency with applications to active remote sensing.

Tsang, L.↗

Volatile/mobile trace elements in Bholghati howardite

Ag, Au, Bi, Cd, Co, Cs, Ga, In, Rb, Sb, Se, Te, Tl, U, and Zn by RNAA are determined in two kinds of material from the Bholghati howardite: a carbonaceous chondrite clast (BH-4) that originated by primary nebular condensation/accretion processes; whole-rock samples and clasts resulting from parent-body processes, especially igneous differentiation. Trace element contents in Bholghati whole-rock (matrix) samples and/or white eucritic clasts in it are high for an achondrite, and variable. These characteristics are indicative of random condensation of volcanic emanations admixed with preexisting eucrite, diogenite, and howardite layers in the HED parent body. Other HED samples show similar patterns, and the conclusion that these reflect extensive volcanism in an asteroidal-sized object is supported by other Consortium noble gas and cosmochronologic measurements.

Wang, Ming-Sheng↗

Active and Passive Radiative Transfer Modeling of the Olympic Mountains Experiment

Sensor forward models are an important tool for interpreting remote sensing observations of geophysical phenomena. By implementing a three-dimensional framework, we can simulate and analyze observations from various sensors on disparate platforms. To demonstrate our model framework, we simulate observations from the Olympic Mountains Experiment (OLYMPEX). The use of cloud model simulations allows us to understand sensor response to cloud ice, falling snow, and other processes and features, and the application of model tools to observations allows us to quantify precipitation.MIIST 3D Forward ModelThe Multi-Instrument Inverse Solver Testbed(MIIST) uses the Atmospheric Radiative TransferSimulator (ARTS) for solving the vector radiativetransfer (RT) equation in up to three spatialdimensions within a spherical geometry• Gas absorptiono Line-by-line calculationso Fast transmittance tables• Hydrometeor scattering solverso Discrete ordinateo RT4 (Evans, 1D)o Radar Single Scattering (1D or 3D)o Monte Carlo (3D)Scattering TablesHigh-fidelity hydrometeor scatteringtables are necessary for accurateand consistent forward modeling ofmulti-frequency observations• Requires full Stokes matriceso And absorption vector• Randomly oriented particleso Discrete Dipole Approximationo Characteristic Basis Function Method(coming soon)• Horizontally-oriented plateso Invariant Imbedding T-matrix MethodCloud Resolving SimulationsCloud resolving simulations (e.g.,NU-WRF) supply output consistentwith ARTS needs• Atmospheric Informationo Temperatureo Pressure / heighto Water vapor• Hydrometeor Profileso ARTS architecture ripe for explicit binmicrophysics• Examples use Morrison 2M schemeThe Olympic Mountains Experiment (OLYMPEX)Validation for GPM of mid-latitudefrontal systems approaching nearcoastalmountains from the ocean• Large collection of ground-based andairborne sensorso Radarso Radiometerso In situ• Contemporaneous with RADEXo Two sets of radar at same frequenciesRadiometer Simulation (3 km NUWRF, 20151203, 15:00)2018.12.14 7Simulate 166 GHz polarizationdifference• Corresponds to the presence of aligned icecrystals• Look at trends for both simulations andobservations• Simulations can tolerate lower resolutiono Larger domainSimulations from Observations: OLYMPEXSimulate sensor response usinggeophysical retrievals as input• Single frequency radar retrievals• Multiple scattering enhancementapparent at W band• Spatially dependent phenomenonModeling Application: 1D Retrievals03 December 2015• DC-8 and ER-2 flightso Focus on APR-3 (DC-8)• Citationo Stacked microphysics legso Qualitative comparisonso Range of frozen habitso Presence of supercooledliquid cloudsResults• Retrievals match probeso Good qualitative match• Bands of increasedreflectivity correspond tolarge Dm and highaggregate fraction• Significant amounts ofsupercooled liquid water

Adams, Ian S.↗

Microscale Constitutive Model Sensitivity on Multiscale Modeling of Fiber Reinforced Composites

Fiber reinforced composites are desirable in applications where low weight and high strength are needed, but are susceptible to microscale variability during manufacturing, making failure predictions difficult. The impact microscale variability has on macroscale mechanical response is difficult to predict due to the computational efficiency needed to simulate many, large, high fidelity, microscale models. In this study, a multiscale approach was taken to model 3-point bend, 4-point bend, and tensile experiments of a unidirectional composite from only having microstructure scans of these samples and constituent properties from literature. These scans were sampled with different sized windows, and statistically equivalent microstructures were generated, then simulated for stiffness, strength, and fracture toughness using an efficient micromechanical. Mesoscale models were created where element sizes equaled microstructure size, and properties were assigned through sampling of microscale simulation results. First, this study showed the effect of using Weibull scaling on constituent matrix strength on macroscale response. Then, a comparison was made between different element sizes and experiments. Finally, model dimensions were fixed, and the effect of randomly distributed local properties alone was examined. Results showed that the scatter of strength and stiffness in the experiments could be predicted well using images of the microscale fiber morphologies and that using stochastic properties produced a 3% coefficient of variation of strength for all experiments.

statistical microstructure↗

Comprehensive Thematic T-Matrix Reference Database: A 2015-2017 Update

The T-matrix method pioneered by Peter C. Waterman is one of the most versatile and efficient numerically exact computer solvers of the time-harmonic macroscopic Maxwell equations. It is widely used for the computation of electromagnetic scattering by single and composite particles, discrete random media, periodic structures (including metamaterials), and particles in the vicinity of plane or rough interfaces separating media with different refractive indices. This paper is the eighth update to the comprehensive thematic database of peer-reviewed T-matrix publications initiated in 2004 and lists relevant publications that have appeared since 2015. It also references a small number of earlier publications overlooked previously.

T-matrix method↗

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↗

Interpretation of satellite-measured bidirectional reflectance from Cirrus cloudy atmospheres

The interpretation of the observed bidirectional reflectance from cirrus cloudy atmospheres is presented. A theoretical model was developed for the computation of the transfer of solar radiation in an anisotropic medium with particular applications to oriented ice crystals in cirrus clouds. In this model, the adding principle for radiative transfer was used with modifications to account for the anisotropy of scattering particles and the associated scattering phase matrix. The single-scattering properties, including the phase function, single-scattering albedo, and extinction cross sections, were used for randomly and horizontally oriented hexagonal ice crystals in radiative transfer computations. The radiative transfer model developed for the cirrus clouds was modified to account for the scattering contributions from the atmosphere and the surface. In order to test the relevance and significance of the ice crystal model for the interpretation of observed bidirectional reflectance from satellites, visible radiances collected on the half hour by the GOES series were selected. The data are calibrated and corrected with the proper filter functions, then navigated to match selected landmark data. A number of clear and cloudy cases during the cirrus IFO of the Fire experiment were chosen for theoretical analyses. The cloud particle shape and size distributions that were taken during satellite overpasses are used in radiative transfer calculations. The sensitivities of the shape, orientation, and size distribution of ice crystals on the reflected intensities at the top of the atmosphere are investigated. Finally, the relative importance of these cloud microphysical properties in the interpretation of satellite bidirectional reflectance are assessed and presented.

Takano, Y.↗

A statistical model of carbon/carbon composite failure

A failure model which considers the stochastic nature of the damage accumulation process is essential to assess reliability and to accurately scale the results from standard test specimens to composite structures. A superior filamentary composite for high temperature applications is composed of carbon fibers in a carbon matrix. Carbon-carbon composites are the strongest known material at very high temperatures. Since there appears to be a significant randomness in C-C material strength which cannot be controlled or detected with current technology, a better model of the material failure based upon statistical principles should be used. Simple applications of the model based upon the limited data provide encouraging results that indicate that better design of test specimens would provide a substantially higher prediction for the design strength of C-C composites. An A-basis strength for the C-C tensile rings from a first stage D-5 billets was estimated. A statistical failure model was developed for these rings which indicates that this strength may be very conservative for larger C-C parts. The analysis may be improved by use of a heterogeneous/noncontinuum finite element approach on the minimechanical level.

Slattery, Kerry T.↗

On Kalman filter solution of space-time interpolation

The approximate Kalman filtering algorithm presented in [1] for image sequence processing can introduce unacceptable negative eigenvalues in the information matrix and can have degraded performance in some applications. The improved algorithm presented in this note guarantees a positive definite information matrix, leading to more stable filter performance.

data↗

Comparison of the Compressive Strengths for Stitched and Toughened Composite Systems

The compression strength of a stitched and a toughened matrix graphite/epoxy composite was determined and compared to a baseline unstitched untoughened composite. Two different layups with a variety of test lengths were tested under both ambient and hot/wet conditions. No significant difference in strength was seen for the different materials when the gage lengths of the specimens were long enough to lead to a buckling failure. For shorter specimens, a 30% reduction in strength from the baseline was seen due to stitching for both a 48-ply quasi-isotropic and a (0/45/0/-45/90/-45/0/45/0)s laminate. Analysis of the results suggested that the decrease in strength was due to increased fiber misalignment due to the stitches. An observed increasing strength with decreasing gage length, which was seen for all materials, was explained with a size effect model. The model assumed a random distribution of flaws (misaligned fibers). The toughened material showed a small increase in strength over the baseline material for both laminates presumably due to the compensating effects of a more compliant matrix and straighter fibers in the toughened material. The hot/wet strength of the stitched and baseline material fell 30% below their ambient strengths for shorter, non-buckling specimen, while the strength of the toughened matrix material only fell 20%. Video images of the failing specimen were recorded and showed local failures prior to global collapse of the specimen. These images support the theory of a random distribution of flaws controlling composite failure. Failed specimen appearance however, seems to be a misleading indication of the cause of failure.

Reeder, James R.↗

Advanced measurement techniques in quantum Monte Carlo: The permutation matrix representation approach

In a typical finite temperature quantum Monte Carlo (QMC) simulation, estimators for simple static observables such as specific heat and magnetization are known. With a great deal of system-specific manual labor, one can sometimes also derive more complicated non-local or even dynamic observable estimators. In contrast, we show that arbitrary static observables can be estimated within the permutation matrix representation (PMR) flavor for any Hamiltonian. We then generalize these results to general imaginary-time correlation functions and non-trivial integrated susceptibilities thereof. Finally, we demonstrate the practical versatility of our method by estimating various non-local, random observables for the transverse-field Ising model on a square lattice and a toy random model.

Permutation matrix representation↗

Comparison of the compressive strengths for stitched and toughened composite systems

The compression strength of a stitched and a toughened matrix graphite/epoxy composite was determined and compared to a baseline unstitched untoughened composite. Two different layups with a variety of test lengths were tested under both ambient and hot/wet conditions. No significant difference in strength was seen for the different materials when the gage lengths of the specimens were long enough to lead to a buckling failure. For shorter specimens, a 30 percent reduction in strength from the baseline was seen due to stitching for both a 48-ply quasi-isotropic and a (0/45/0/-45/90/-45/0/45/0)s laminate. Analysis of the results suggested that the decrease in strength was due to increased fiber misalignment due to the stitches. An observed increasing strength with decreasing gage length, which was seen for all materials, was explained with a size effect model. The model assumed a random distribution of flaws (misaligned fibers). The toughened materials showed a small increase in strength over the baseline material for both laminates presumably due to the compensating effects of a more compliant matrix and straighter fibers in the toughened material. The hot/wet strength of the stitched and baseline material fell 30 percent below their ambient strengths for shorter, nonbuckling specimen, while the strength of the toughened matrix material only fell 20 percent. Video images of the failing specimen were recorded and showed local failures prior to global collapse of the specimen. These images support the theory of a random distribution of flaws controlling composite failure. Failed specimen appearance, however, seems to be a misleading indication of the cause of failure.

Reeder, James R.↗

Improving the chi-squared approximation for bivariate normal tolerance regions

Let X be a two-dimensional random variable distributed according to N2(mu,Sigma) and let bar-X and S be the respective sample mean and covariance matrix calculated from N observations of X. Given a containment probability beta and a level of confidence gamma, we seek a number c, depending only on N, beta, and gamma such that the ellipsoid R = (x: (x - bar-X)'S(exp -1) (x - bar-X) less than or = c) is a tolerance region of content beta and level gamma; i.e., R has probability gamma of containing at least 100 beta percent of the distribution of X. Various approximations for c exist in the literature, but one of the simplest to compute -- a multiple of the ratio of certain chi-squared percentage points -- is badly biased for small N. For the bivariate normal case, most of the bias can be removed by simple adjustment using a factor A which depends on beta and gamma. This paper provides values of A for various beta and gamma so that the simple approximation for c can be made viable for any reasonable sample size. The methodology provides an illustrative example of how a combination of Monte-Carlo simulation and simple regression modelling can be used to improve an existing approximation.

Feiveson, Alan H.↗

Deconsolidation and Leach Burn Leach of Seven As Irradiated AGR 5/6/7 TRISO Fuel Compacts from Capsules 2, 3, 4, and 5

Seven as-irradiated Advanced Gas Reactor (AGR) 5/6/7 compacts underwent destructive post-irradiation examination via deconsolidation-leach-burn leach at Idaho National Laboratory (INL). The selection of the compacts extended the upper and lower limits of time-average volume-average (TAVA) temperature for compacts that had gone through deconsolidation-leach-burn leach so far. The measured inventories of fission products and actinides in the compact matrix and outer pyrolytic carbon were reported. Results indicated unexpectedly higher rates of fuel kernel leaching compared to compacts from AGR-1 and AGR-2. These failure rates were attributed to damage during post-irradiation sample handling, rather than irradiation itself. AGR-5/6/7 compacts have little or no matrix coverage for some particles at the top and bottom ends of cylindrical fuel compacts, making them more fragile. Sixty particles were randomly sampled from each compact, and the gamma results were reported. Three SiC shells from were identified, one of which showed signs of chemical attack in the high-irradiation temperature compact.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Exploring the Connection Between Sampling Problems in Bayesian Inference and Statistical Mechanics

The Bayesian and statistical mechanical communities often share the same objective in their work - estimating and integrating probability distribution functions (pdfs) describing stochastic systems, models or processes. Frequently, these pdfs are complex functions of random variables exhibiting multiple, well separated local minima. Conventional strategies for sampling such pdfs are inefficient, sometimes leading to an apparent non-ergodic behavior. Several recently developed techniques for handling this problem have been successfully applied in statistical mechanics. In the multicanonical and Wang-Landau Monte Carlo (MC) methods, the correct pdfs are recovered from uniform sampling of the parameter space by iteratively establishing proper weighting factors connecting these distributions. Trivial generalizations allow for sampling from any chosen pdf. The closely related transition matrix method relies on estimating transition probabilities between different states. All these methods proved to generate estimates of pdfs with high statistical accuracy. In another MC technique, parallel tempering, several random walks, each corresponding to a different value of a parameter (e.g. "temperature"), are generated and occasionally exchanged using the Metropolis criterion. This method can be considered as a statistically correct version of simulated annealing. An alternative approach is to represent the set of independent variables as a Hamiltonian system. Considerab!e progress has been made in understanding how to ensure that the system obeys the equipartition theorem or, equivalently, that coupling between the variables is correctly described. Then a host of techniques developed for dynamical systems can be used. Among them, probably the most powerful is the Adaptive Biasing Force method, in which thermodynamic integration and biased sampling are combined to yield very efficient estimates of pdfs. The third class of methods deals with transitions between states described by rate constants. These problems are isomorphic with chemical kinetics problems. Recently, several efficient techniques for this purpose have been developed based on the approach originally proposed by Gillespie. Although the utility of the techniques mentioned above for Bayesian problems has not been determined, further research along these lines is warranted

Pohorille, Andrew↗