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 343 records · Page 19

Asymptotic Solutions for Optical Properties of Large Particles with Strong Absorption

The transverse wave condition is not applicable to the refracted electromagnetic wave within the context of geometrical optics when absorption is involved. Either the transverse magnetic (TM) or the transverse electric (TE) wave condition can be assumed for the wave to locally satisfy the electromagnetic boundary condition in a ray-tracing calculation. The assumed wave mode affects both the reflection and the refraction coefficients. As a result, nonunique solutions for these coefficients are inevitable. In this study the appropriate solutions for the Fresnel reflection-refraction coefficients are identified in light-scattering calculations based on the ray-tracing technique. In particular, a 3x2 refraction or transmission matrix is derived to account for the inhomogeneity of the refracted wave in an absorbing medium. An asymptotic solution that completely includes the effect of medium absorption on Fresnel coefficients is obtained for the scattering properties of a general polyhedral particle. Numerical results are presented for hexagonal plates and columns with both preferred and random orientations.

Yang, Ping↗

The effect of adiabatic focusing upon charged-particle propagation in random magnetic fields

The charged particles considered are scattered by random fields while they propagate along the diverging lines of force of a spatially inhomogeneous guiding field. Their longitudinal transport is described in terms of the eigenfunctions of a Sturm-Liouville operator which incorporates the effect of adiabatic focusing along with that of scattering. The relaxation times and characteristic velocities which appear in this matrix formulation of the transport problem are graphed and tabulated. Explicit formulas which describe the particle-density profile that results from a localized impulsive injection are derived for two different regimes. In the first regime, where focusing is relatively weak, a diffusive mode of propagation is dominant, but coherent modes are also present, and they become prominent as the intensity of focusing increases. In the second regime, where focusing is strong and where diffusion does not occur, the propagation is purely coherent. The existence of this supercoherent mode of particle transport opens up many possibilities for the interpretation of astrophysical phenomena.

Earl, J. A.↗

Surrogate models for linear response

Linear response theory is a well-established method in physics and chemistry for exploring excitations of many-body systems. In particular, the quasiparticle random-phase approximation (QRPA) provides a powerful microscopic framework by building excitations on top of the mean-field vacuum; however, its high computational cost limits model calibration and uncertainty quantification studies. Here, we present two complementary QRPA surrogate models and apply them to study response functions of finite nuclei. One is a reduced-order model that exploits the underlying QRPA structure, while the other utilizes the recently developed parametric matrix model algorithm to construct a map between the system’s Hamiltonian and observables. Our benchmark applications, the calculation of the electric dipole polarizability of 180 Yb and the 𝛽-decay half-life of 80 Ni, show that both emulators can achieve 0.1%–1% accuracy while offering a 6–7 orders of magnitude speedup compared to state-of-the-art QRPA solvers. These results demonstrate that the developed QRPA emulators are well positioned to enable Bayesian calibration and large-scale studies of computationally expensive physics models describing the properties of many-body systems.

Beta decay↗

Absorption and scattering of light by nonspherical particles

Using the example of the polarization of scattered light, it is shown that the scattering matrices for identical, randomly ordered particles and for spherical particles are unequal. The spherical assumptions of Mie theory are therefore inconsistent with the random shapes and sizes of atmospheric particulates. The implications for corrections made to extinction measurements of forward scattering light are discussed. Several analytical methods are examined as potential bases for developing more accurate models, including Rayleigh theory, Fraunhoffer Diffraction theory, anomalous diffraction theory, Rayleigh-Gans theory, the separation of variables technique, the Purcell-Pennypacker method, the T-matrix method, and finite difference calculations.

Bohren, C. F.↗

Dermal Aged and Fetal Fibroblasts Realign in Response to Mechanical Strain

Integrins specifically recognize and bind extracellular matrix components, providing physical anchor points and functional setpoints. Focal adhesion complexes, containing integrin and cytoskeletal proteins, are potential mechanoreceptors, poised to distribute applied forces through the cytoskeleton. Pursuing the hypothesis that cells both perceive and respond to external force, we applied a stretch/relaxation regimen to normal human fetal and aged dermal fibroblast monolayers cultured on flexible membranes. The frequency and magnitude of the applied force is precisely controlled by the Flexercell Unit(Trademark). A protocol of stretch (20% elongation of the monolayer) at a frequency of 6 cycles/min caused a progressive change from a randomly distributed pattern of cells to a symmetric, radial distribution with cells aligned parallel to the applied force. We have coined the term 'orienteering' as the process of active alignment of cells in response to applied force. Cytochalasin D was added in graded doses to investigate the role of the actin cytoskeleton in force perception and transmission. A clear dose response was found; at high concentrations orienteering was abolished; and the drug's impact was reversible. The two cell strains used were similar in their alignment behavior and in their responses to cytochalasin D. Orienteering was influenced by cell density, and the cell strains studied differed in this respect. Fetal cells, unlike their aged counterparts, failed to orient at high cell density. In both cell strains, mid-density cultures aligned rapidly and sparse cultures lagged. These results indicate that both cell-cell adhesion and cytoskeleton integrity are critical in mediating the orienteering response. Differences between these two cell strains may relate to their expression of extracellular matrix molecules (fibronectin, collagen type 1) integrins and their relative binding affinities.

Sawyer, Christine↗

Computation of convective flow with gravity modulation in rectangular cavities

In this work, a computational study is presented for the investigation of gravity modulation (g-jitter) effects in thermally driven cavity flows at terrestrial and microgravity environments. The two-dimensional, time-dependent Navier-Stokes equations are numerically integrated by a time-split method using direct matrix solvers. Computations at terrestrial gravity are utilized to assess the effects of adiabatic side-wall boundary conditions as well as the full nonlinearity of the governing equations on the sinusoidally forced Benard problem studied by Gresho and Sani. The low-g calculations focus on the establishment of critical frequency ranges and consider the effects of modulation direction and randomness. The applicability of linear analysis in the excitable frequency range at low g is also discussed.

Biringen, S.↗

A reliable algorithm for optimal control synthesis

In recent years, powerful design tools for linear time-invariant multivariable control systems have been developed based on direct parameter optimization. In this report, an algorithm for reliable optimal control synthesis using parameter optimization is presented. Specifically, a robust numerical algorithm is developed for the evaluation of the H(sup 2)-like cost functional and its gradients with respect to the controller design parameters. The method is specifically designed to handle defective degenerate systems and is based on the well-known Pade series approximation of the matrix exponential. Numerical test problems in control synthesis for simple mechanical systems and for a flexible structure with densely packed modes illustrate positively the reliability of this method when compared to a method based on diagonalization. Several types of cost functions have been considered: a cost function for robust control consisting of a linear combination of quadratic objectives for deterministic and random disturbances, and one representing an upper bound on the quadratic objective for worst case initial conditions. Finally, a framework for multivariable control synthesis has been developed combining the concept of closed-loop transfer recovery with numerical parameter optimization. The procedure enables designers to synthesize not only observer-based controllers but also controllers of arbitrary order and structure. Numerical design solutions rely heavily on the robust algorithm due to the high order of the synthesis model and the presence of near-overlapping modes. The design approach is successfully applied to the design of a high-bandwidth control system for a rotorcraft.

Vansteenwyk, Brett↗

Estimation of 3D Woven Design Sensitivities Using a Rapid Multiscale Analysis Technique

Highly-refined finite element models of three-dimension (3D) woven composite systems currently require excessive computational demands that limit their use in sensitivity analysis, uncertainty quantification, and optimization. An alternative analysis methodology was developed using the NASA Multiscale Analysis Tool (NASMAT) where multiscale models of a 3D woven composite (including inter-tow matrix voids and constituent failure) can be completed on a single central processing unit (CPU) on the order of ~30s. To develop inputs and validation data for the NASMAT model, coupon and acid-digesting testing and x-ray computed tomography were performed. The NASMAT inputs were parameterized using a set of 25 input variables and distributions. These inputs were randomly sampled to generate a total of 100,000 NASMAT analyses that could be used to understand the influence of different material and geometric properties on the warp and weft-direction stiffness and strength. These analyses (including pre/post-processing) were performed in less than eight hours on a 120 CPU cluster. The computational efficiency of the NASMAT model enabled a sensitivity analysis to be performed, and dominant input variables were able to be identified. Key results were consistent with theoretical and experimental observations for the specific 3D woven system studied in this work.

NASMAT↗

Estimation of 3D Woven Design Sensitivities Using a Rapid Multiscale Analysis Technique

Highly-refined finite element models of three-dimension (3D) woven composite systems currently require excessive computational demands that limit their use in sensitivity analysis, uncertainty quantification, and optimization. An alternative analysis methodology was developed using the NASA Multiscale Analysis Tool (NASMAT) where multiscale models of a 3D woven composite (including inter-tow matrix voids and constituent failure) can be completed on a single central processing unit(CPU)on the order of ~30 s. To develop inputs and validation data for the NASMAT model, coupon and acid-digesting testing and x-ray computed tomography were performed. The NASMAT inputs were parameterized using a set of 25 input variables and distributions. These inputs were randomly sampled to generate a total of 100,000 NASMAT analyses that could be used to understand the influence of different material and geometric properties on the warp and weft-direction stiffness and strength. These analyses (including pre/post-processing) were performed in less than eight hours on a 120 CPU cluster. The computational efficiency of the NASMAT model enabled a sensitivity analysis to be performed, and dominant input variables were able to be identified. Key results were consistent with theoretical and experimental observations for the specific 3D woven system studied in this work.

NASMAT↗

Energy partitioning of gaseous ions in an electric field.

The partitioning of ion energy among thermal energy, drift energy, and random-field energy is studied by solution of the Boltzmann equation. An expansion in powers of the square of the electric field strength is obtained by Kihara's method. Numerical calculations for several ion-neutral force laws show that Wannier's constant mean-free-time model gives a reasonable first approximation. The formal extension to multicomponent mixtures is also given. The matrix elements obtained are tabulated, and can be used to study the field dependence of other moments of the ion-distribution function.

Hahn, H.-S.↗

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.↗

Changes in gravity inhibit lymphocyte locomotion through type I collagen

Immunity relies on the circulation of lymphocytes through many different tissues including blood vessels, lymphatic channels, and lymphoid organs. The ability of lymphocytes to traverse the interstitium in both nonlymphoid and lymphoid tissues can be determined in vitro by assaying their capacity to locomote through Type I collagen. In an attempt to characterize potential causes of microgravity-induced immunosuppression, we investigated the effects of simulated microgravity on human lymphocyte function in vitro using a specialized rotating-wall vessel culture system developed at the Johnson Space Center. This very low shear culture system randomizes gravitational vectors and provides an in vitro approximation of microgravity. In the randomized gravity of the rotating-wall vessel culture system, peripheral blood lymphocytes did not locomote through Type I collagen, whereas static cultures supported normal movement. Although cells remained viable during the entire culture period, peripheral blood lymphocytes transferred to unit gravity (static culture) after 6 h in the rotating-wall vessel culture system were slow to recover and locomote into collagen matrix. After 72 h in the rotating-wall vessel culture system and an additional 72 h in static culture, peripheral blood lymphocytes did not recover their ability to locomote. Loss of locomotory activity in rotating-wall vessel cultures appears to be related to changes in the activation state of the lymphocytes and the expression of adhesion molecules. Culture in the rotating-wall vessel system blunted the ability of peripheral blood lymphocytes to respond to polyclonal activation with phytohemagglutinin. Locomotory response remained intact when peripheral blood lymphocytes were activated by anti-CD3 antibody and interleukin-2 prior to introduction into the rotating-wall vessel culture system. Thus, in addition to the systemic stress factors that may affect immunity, isolated lymphocytes respond to gravitational changes by ceasing locomotion through model interstitium. These in vitro investigations suggest that microgravity induces non-stress-related changes in cell function that may be critical to immunity. Preliminary analysis of locomotion in true microgravity revealed a substantial inhibition of cellular movement in Type I collagen. Thus, the rotating-wall vessel culture system provides a model for analyzing the microgravity-induced inhibition of lymphocyte locomotion and the investigation of the mechanisms related to lymphocyte movement.

Weightlessness↗

Adjoints and Low-rank Covariance Representation

Quantitative measures of the uncertainty of Earth System estimates can be as important as the estimates themselves. Second moments of estimation errors are described by the covariance matrix, whose direct calculation is impractical when the number of degrees of freedom of the system state is large. Ensemble and reduced-state approaches to prediction and data assimilation replace full estimation error covariance matrices by low-rank approximations. The appropriateness of such approximations depends on the spectrum of the full error covariance matrix, whose calculation is also often impractical. Here we examine the situation where the error covariance is a linear transformation of a forcing error covariance. We use operator norms and adjoints to relate the appropriateness of low-rank representations to the conditioning of this transformation. The analysis is used to investigate low-rank representations of the steady-state response to random forcing of an idealized discrete-time dynamical system.

Tippett, Michael K.↗

Toward memory-efficient melt pool monitoring: a classification framework using event-based imaging and sparse sensing technique

Vision sensors like CMOS and CCD cameras are often used for in-process monitoring of melt pools in laser-based additive and welding processes, but they require transferring large amounts of data and computational processing resources. Event-based neuromorphic imagery, on the other hand, detects only the change in pixel intensity, thus potentially reducing the data amount and latency. With an event imager, this study develops a framework for melt pool condition classification, including image construction, time scale selection, optimal pixel selection, and sparse classification, to achieve a highly memory-efficient scheme. These are based on sparse sensing techniques with singular value decomposition (SVD) and QR pivoting, the two fundamental matrix transformations for linear dimensionality reduction. The framework is then validated by classifying a controlled experiment by exciting various mode shapes of liquid gallium pools of varying depths (3, 6, and 8 mm). At 200 pixels, the classifier can reach overall accuracy of 75%, while at 2000 pixels (0.013% of the total possible pixels), the accuracy is nearly 90% (89.86%). At the same number of pixels, random selection can only achieve 46% and 67%, respectively. The memory savings of the sparsely sampled event data compared to a conventional imager is about 500 times. In addition to performance, implementation and limitations of the framework are also discussed.

42 ENGINEERING↗

Effective propagation constants in dense random media under effective medium approximation

The effective medium approximation is applied to study the effective propagation constants in a dense random medium. The dyadic Green's function is introduced to establish the effective medium approximation formalism for electromagnetic waves. The multiple scattering equations and the Lippmann-Schwinger equations for the transition operator are obtained in the configuration average form. The dispersion equations of multiple scattering is derived by using a standard method in quantum mechanics. To obtain an expression for the effective propagation constants the matrix elements of the configuration average dyadic transiton operator are calculated in momentum representation. Numerical illustations are carried out to demonstrate the difference in the effective propagation constants between the use of this approximation and the well-known quasi-crystalline approximation. A comparison is made with measured loss tangent in dry snow.

Zhu, P. Y.↗

Nominal 30-M Cropland Extent Map of Continental Africa by Integrating Pixel-Based and Object-Based Algorithms Using Sentinel-2 and Landsat-8 Data on Google Earth Engine

A satellite-derived cropland extent map at high spatial resolution (30-m or better) is a must for food and water security analysis. Precise and accurate global cropland extent maps, indicating cropland and non-cropland areas, is a starting point to develop high-level products such as crop watering methods (irrigated or rainfed), cropping intensities (e.g., single, double, or continuous cropping), crop types, cropland fallows, as well as assessment of cropland productivity (productivity per unit of land), and crop water productivity (productivity per unit of water). Uncertainties associated with the cropland extent map have cascading effects on all higher-level cropland products. However, precise and accurate cropland extent maps at high spatial resolution over large areas (e.g., continents or the globe) are challenging to produce due to the small-holder dominant agricultural systems like those found in most of Africa and Asia. Cloud-based Geospatial computing platforms and multi-date, multi-sensor satellite image inventories on Google Earth Engine offer opportunities for mapping croplands with precision and accuracy over large areas that satisfy the requirements of broad range of applications. Such maps are expected to provide highly significant improvements compared to existing products, which tend to be coarser in resolution, and often fail to capture fragmented small-holder farms especially in regions with high dynamic change within and across years. To overcome these limitations, in this research we present an approach for cropland extent mapping at high spatial resolution (30-m or better) using the 10-day, 10 to 20-m, Sentinel-2 data in combination with 16-day, 30-m, Landsat-8 data on Google Earth Engine (GEE). First, nominal 30-m resolution satellite imagery composites were created from 36,924 scenes of Sentinel-2 and Landsat-8 images for the entire African continent in 2015-2016. These composites were generated using a median-mosaic of five bands (blue, green, red, near-infrared, NDVI) during each of the two periods (period 1: January-June 2016 and period 2: July-December 2015) plus a 30-m slope layer derived from the Shuttle Radar Topographic Mission (SRTM) elevation dataset. Second, we selected Cropland/Non-cropland training samples (sample size 9791) from various sources in GEE to create pixel-based classifications. As supervised classification algorithm, Random Forest (RF) was used as the primary classifier because of its efficiency, and when over-fitting issues of RF happened due to the noise of input training data, Support Vector Machine (SVM) was applied to compensate for such defects in specific areas. Third, the Recursive Hierarchical Segmentation (RHSeg) algorithm was employed to generate an object-oriented segmentation layer based on spectral and spatial properties from the same input data. This layer was merged with the pixel-based classification to improve segmentation accuracy. Accuracies of the merged 30-m crop extent product were computed using an error matrix approach in which 1754 independent validation samples were used. In addition, a comparison was performed with other available cropland maps as well as with LULC maps to show spatial similarity. Finally, the cropland area results derived from the map were compared with UN FAO statistics. The independent accuracy assessment showed a weighted overall accuracy of 94, with a producers accuracy of 85.9 (or omission error of 14.1), and users accuracy of 68.5 (commission error of 31.5) for the cropland class. The total net cropland area (TNCA) of Africa was estimated as 313 Mha for the nominal year 2015.

Cropland mapping; cropland areas; 30-m; Landsat-8;↗

Phenylquinoxaline copolymers

The preparation of a series of homopolymers and block and random type phenylquinoxaline copolymers is described together with their important characteristics and preliminary measurements of composite properties. Phenylquinoxaline homopolymers consisting of low, intermediate, and high molecular weight versions were prepared from the reaction of 3,3'-diaminobenzidine (DAB) and 3,3',4,4'-tetraaminobenzophenone (TAB), each with p-bis(phenylglyoxalyl)benzene as solutions (18% solids content) in a 1:1 mixture of m-cresol and xylene. Three random and three block copolymers were prepared from the reaction of various amounts of DAB and TAB with p-bis(phenylglyoxalyl)benzene. The copolymers were prepared with DAB to TAB distribution ratios of 3:1, 1:1, and 1:3. The work was performed primarily to determine if phenylquinoxaline copolymers exhibit any advantage over the homopolymers for use as functional or structural resins, especially as a matrix material in combination with a high-modulus graphite fiber reinforcement.

Hergenrother, P. M.↗

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↗