Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Parallel Matrix Multiplication”

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.

62 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

MATRIX-VBS (v1.0): Implementing an Evolving Organic Aerosol Volatility in an Aerosol Microphysics Model

The gas-particle partitioning and chemical aging of semi-volatile organic aerosol are presented in a newly developed box model scheme, where its effect on the growth, composition, and mixing state of particles is examined. The volatility-basis set (VBS) framework is implemented into the aerosol microphysical scheme MATRIX (Multiconfiguration Aerosol TRacker of mIXing state), which resolves mass and number aerosol concentrations and in multiple mixing-state classes. The new scheme, MATRIX-VBS, has the potential to significantly advance the representation of organic aerosols in Earth system models by improving upon the conventional representation as non-volatile particulate organic matter, often also with an assumed fixed size distribution. We present results from idealized cases representing Beijing, Mexico City, a Finnish forest, and a southeastern US forest, and investigate the evolution of mass concentrations and volatility distributions for organic species across the gas and particle phases, as well as assessing their mixing state among aerosol populations. Emitted semi-volatile primary organic aerosols evaporate almost completely in the intermediate-volatility range, while they remain in the particle phase in the low-volatility range. Their volatility distribution at any point in time depends on the applied emission factors, oxidation by OH radicals, and temperature. We also compare against parallel simulations with the original scheme, which represented only the particulate and non-volatile component of the organic aerosol, examining how differently the condensed-phase organic matter is distributed across the mixing states in the model. The results demonstrate the importance of representing organic aerosol as a semi-volatile aerosol, and explicitly calculating the partitioning of organic species between the gas and particulate phases.

volatility-basis set↗
NASA NTRS2000

Ordering Unstructured Meshes for Sparse Matrix Computations on Leading Parallel Systems

The ability of computers to solve hitherto intractable problems and simulate complex processes using mathematical models makes them an indispensable part of modern science and engineering. Computer simulations of large-scale realistic applications usually require solving a set of non-linear partial differential equations (PDES) over a finite region. For example, one thrust area in the DOE Grand Challenge projects is to design future accelerators such as the SpaHation Neutron Source (SNS). Our colleagues at SLAC need to model complex RFQ cavities with large aspect ratios. Unstructured grids are currently used to resolve the small features in a large computational domain; dynamic mesh adaptation will be added in the future for additional efficiency. The PDEs for electromagnetics are discretized by the FEM method, which leads to a generalized eigenvalue problem Kx = AMx, where K and M are the stiffness and mass matrices, and are very sparse. In a typical cavity model, the number of degrees of freedom is about one million. For such large eigenproblems, direct solution techniques quickly reach the memory limits. Instead, the most widely-used methods are Krylov subspace methods, such as Lanczos or Jacobi-Davidson. In all the Krylov-based algorithms, sparse matrix-vector multiplication (SPMV) must be performed repeatedly. Therefore, the efficiency of SPMV usually determines the eigensolver speed. SPMV is also one of the most heavily used kernels in large-scale numerical simulations.

Oliker, Leonid↗
NASA NTRS2004

Error Modeling of Multi-baseline Optical Truss: Application to SIM Metrology Truss Field Dependent Error - Part II

The current design of the Space Interferometry Mission (SIM) employs a 19 laser-metrology-beam system (also called L19 external metrology truss) to monitor changes of distances between the fiducials of the flight system's multiple baselines. The function of the external metrology truss is to aid in the determination of the time-variations of the interferometer baseline. The largest contributor to truss error occurs in SIM wide-angle observations when the articulation of the siderostat mirrors (in order to gather starlight from different sky coordinates) brings to light systematic errors due to offsets at levels of instrument components (which include comer cube retro-reflectors, etc.). This error is labeled external metrology wide-angle field-dependent error. Physics-based model of field-dependent error at single metrology gauge level is developed and linearly propagated to errors in interferometer delay. In this manner delay error sensitivity to various error parameters or their combination can be studied using eigenvalue/eigenvector analysis. Also validation of physics-based field-dependent model on SIM testbed lends support to the present approach. As a first example, dihedral error model is developed for the comer cubes (CC) attached to the siderostat mirrors. Then the delay errors due to this effect can be characterized using the eigenvectors of composite CC dihedral error. The essence of the linear error model is contained in an error-mapping matrix. A corresponding Zernike component matrix approach is developed in parallel, first for convenience of describing the RMS of errors across the field-of-regard (FOR), and second for convenience of combining with additional models. Average and worst case residual errors are computed when various orders of field-dependent terms are removed from the delay error. Results of the residual errors are important in arriving at external metrology system component requirements. Double CCs with ideally co-incident vertices reside with the siderostat. The non-common vertex error (NCVE) is treated as a second example. Finally combination of models, and various other errors are discussed.

comer cube retro-reflector↗
NASA NTRS1994

Errors induced by the neglect of polarization in radiance calculations for Rayleigh-scattering atmospheres

Although neglecting polarization and replacing the rigorous vector radiative transfer equation by its approximate scalar counterpart has no physical background, it is a widely used simplification when the incident light is unpolarized and only the intensity of the reflected light is to be computed. We employ accurate vector and scalar multiple-scattering calculations to perform a systematic study of the errors induced by the neglect of polarization in radiance calculations for a homogeneous, plane-parallel Rayleigh-scattering atmosphere (with and without depolarization) above a Lambertian surface. Specifically, we calculate percent errors in the reflected intensity for various directions of light incidence and reflection, optical thicknesses of the atmosphere, single-scattering albedos, depolarization factors, and surface albedos. The numerical data displayed can be used to decide whether or not the scalar approximation may be employed depending on the parameters of the problem. We show that the errors decrease with increasing depolarization factor and/or increasing surface albedo. For conservative or nearly conservative scattering and small surface albedos, the errors are maximum at optical thicknesses of about 1. The calculated errors may be too large for some practical applications, and, therefore, rigorous vector calculations should be employed whenever possible. However, if approximate scalar calculations are used, we recommend to avoid geometries involving phase angles equal or close to 0 deg and 90 deg, where the errors are especially significant. We propose a theoretical explanation of the large vector/scalar differences in the case of Rayleigh scattering. According to this explanation, the differences are caused by the particular structure of the Rayleigh scattering matrix and come from lower-order (except first-order) light scattering paths involving right scattering angles and right-angle rotations of the scattering plane.

Mishchenko, M. I.↗
NASA NTRS2008

Update on Development of SiC Multi-Chip Power Modules

Progress has been made in a continuing effort to develop multi-chip power modules (SiC MCPMs). This effort at an earlier stage was reported in 'SiC Multi-Chip Power Modules as Power-System Building Blocks' (LEW-18008-1), NASA Tech Briefs, Vol. 31, No. 2 (February 2007), page 28. The following recapitulation of information from the cited prior article is prerequisite to a meaningful summary of the progress made since then: 1) SiC MCPMs are, more specifically, electronic power-supply modules containing multiple silicon carbide power integrated-circuit chips and silicon-on-insulator (SOI) control integrated-circuit chips. SiC MCPMs are being developed as building blocks of advanced expandable, reconfigurable, fault-tolerant power-supply systems. Exploiting the ability of SiC semiconductor devices to operate at temperatures, breakdown voltages, and current densities significantly greater than those of conventional Si devices, the designs of SiC MCPMs and of systems comprising multiple SiC MCPMs are expected to afford a greater degree of miniaturization through stacking of modules with reduced requirements for heat sinking; 2) The stacked SiC MCPMs in a given system can be electrically connected in series, parallel, or a series/parallel combination to increase the overall power-handling capability of the system. In addition to power connections, the modules have communication connections. The SOI controllers in the modules communicate with each other as nodes of a decentralized control network, in which no single controller exerts overall command of the system. Control functions effected via the network include synchronization of switching of power devices and rapid reconfiguration of power connections to enable the power system to continue to supply power to a load in the event of failure of one of the modules; and, 3) In addition to serving as building blocks of reliable power-supply systems, SiC MCPMs could be augmented with external control circuitry to make them perform additional power-handling functions as needed for specific applications. Because identical SiC MCPM building blocks could be utilized in such a variety of ways, the cost and difficulty of designing new, highly reliable power systems would be reduced considerably. This concludes the information from the cited prior article. The main activity since the previously reported stage of development was the design, fabrication, and testing a 120- VDC-to-28-VDC modular power-converter system composed of eight SiC MCPMs in a 4 (parallel)-by-2 (series) matrix configuration, with normally-off controllable power switches. The SiC MCPM power modules include closed-loop control subsystems and are capable of operating at high power density or high temperature. The system was tested under various configurations, load conditions, load-transient conditions, and failure-recovery conditions. Planned future work includes refinement of the demonstrated modular system concept and development of a new converter hardware topology that would enable sharing of currents without the need for communication among modules. Toward these ends, it is also planned to develop a new converter control algorithm that would provide for improved sharing of current and power under all conditions, and to implement advanced packaging concepts that would enable operation at higher power density.

Lostetter, Alexander↗
NASA NTRS1999

Integration of a Decentralized Linear-Quadratic-Gaussian Control into GSFC's Universal 3-D Autonomous Formation Flying Algorithm

A decentralized control is investigated for applicability to the autonomous formation flying control algorithm developed by GSFC for the New Millenium Program Earth Observer-1 (EO-1) mission. This decentralized framework has the following characteristics: The approach is non-hierarchical, and coordination by a central supervisor is not required; Detected failures degrade the system performance gracefully; Each node in the decentralized network processes only its own measurement data, in parallel with the other nodes; Although the total computational burden over the entire network is greater than it would be for a single, centralized controller, fewer computations are required locally at each node; Requirements for data transmission between nodes are limited to only the dimension of the control vector, at the cost of maintaining a local additional data vector. The data vector compresses all past measurement history from all the nodes into a single vector of the dimension of the state; and The approach is optimal with respect to standard cost functions. The current approach is valid for linear time-invariant systems only. Similar to the GSFC formation flying algorithm, the extension to linear LQG time-varying systems requires that each node propagate its filter covariance forward (navigation) and controller Riccati matrix backward (guidance) at each time step. Extension of the GSFC algorithm to non-linear systems can also be accomplished via linearization about a reference trajectory in the standard fashion, or linearization about the current state estimate as with the extended Kalman filter. To investigate the feasibility of the decentralized integration with the GSFC algorithm, an existing centralized LQG design for a single spacecraft orbit control problem is adapted to the decentralized framework while using the GSFC algorithm's state transition matrices and framework. The existing GSFC design uses both reference trajectories of each spacecraft in formation and by appropriate choice of coordinates and simplified measurement modeling is formulated as a linear time-invariant system. Results for improvements to the GSFC algorithm and a multiple satellite formation will be addressed. The goal of this investigation is to progressively relax the assumptions that result in linear time-invariance, ultimately to the point of linearization of the non-linear dynamics about the current state estimate as in the extended Kalman filter. An assessment will then be made about the feasibility of the decentralized approach to the realistic formation flying application of the EO-1/Landsat 7 formation flying experiment.

Folta, David C.↗
NASA NTRS1993

Relationships between pigment composition variation and reflectance for plant species from a coastal savannah in California

Advances in imaging spectroscopy have indicated that remotely sensed reflectance measurements of the plant canopy may be used to identify and qualify some classes of canopy biochemicals; however, the manner in which differences in biochemical compositions translate into differences is not well understood. Most frequently, multiple linear regression routines have been used to correlate narrow band reflectance values with measured biochemical concentrations. Although some success has been achieved with such methods for given data sets, the bands selected by multiple regression are not consistent between data sets, nor is it always clear what physical or biological basis underlies the correlation. To examine the relationship between biochemical concentration and leaf reflectance signal we chose to focus on the visible spectrum where the primary biochemical absorbances are due to photosynthetic pigments. Pigments provide a range of absorbance features, occur over a range of concentrations in natural samples, and are ecophysiologically important. Concentrations of chlorophyll, for example, have been strongly correlated to foliar nitrogen levels within a species and to photosynthetic capacity across many species. In addition pigments effectively absorb most of the photosynthetically active radiation between 400-700 nm, a spectral region for which silicon detectors have good signal/noise characteristics. Our strategy has been to sample a variety of naturally occurring species to measure leaf reflectance and pigment compositions. We hope to extend our understanding of pigment reflectance effects to interpret small overlapping absorbances of other biochemicals in the infrared region. For this reason, selected samples were also tested to determine total nitrogen, crude protein, cellulose, and lignin levels. Leaf reflectance spectra measured with AVIRIS bandwidths and wavelengths were compared between species and within species and for differences between seasons, for changes in the the shape of the spectra. We attempt to statistically correlate these shape changes with differences in pigment compositions. In parallel with our comparisons of pigment composition and leaf reflectance, we have modified the PROSPECT leaf reflectance model to test the contributions of pigments or pigment group concentrations. PROSPECT considers a leaf as a multi-layer dielectric plane with an uneven surface. Jacquemoud adapted the basic analysis of Allen for surface effects, a leaf thickness factor, and the absorption of water and chlorophyll (actually all pigments) and the plant matrix. Our modifications to PROSPECT in the forward direction include breaking out the pigment concentration parameter into separate components for chlorophyll a and b and a number of xanthophylls and carotenes, and introducing a shift and convolution function to model the spread and shift from their in vitro measurements to their in vivo state. Further, we have considered how the matrix elements (i.e., all biochemicals and structural effects not modeled explicity) vary with species.

Ustin, Susan L.↗
NASA NTRS1997

Breccia Formation at a Complex Impact Crater: Slate Islands, Lake Superior, Ontario, Canada

The Slate Islands impact structure is the eroded remnant of a approximately 30-32 km-diameter complex impact structure located in northern Lake Superior, Ontario, Canada. Target rocks are Archean supracrustal and igneous rocks and Proterozoic metavolcanics, metasediments, and diabase. A wide variety of breccias occurs on the islands, many of which contain fragments exhibiting shock metamorphic features. Aphanitic, narrow and inclusion-poor pseudotachylite veins, commonly with more or less parallel boundaries and apophyses branching off them, represent the earliest breccias formed during the compression stage of the impact process. Coarse-grained, polymictic elastic matrix breccias form small to very large, inclusion-rich dikes and irregularly shaped bodies that may contain altered glass fragments. These breccias have sharp contacts with their host rocks and include a wide range of fragment types some of which were transported over minimum distances of approximately 2 km away from the center of the structure. They cut across pseudotachylite veins and contain inclusions of them. Field and petrographic evidence indicate that these polymictic breccias formed predominantly during the excavation and central uplift stages of the impact process. Monomictic breccias, characterized by angular fragments and transitional contacts with their host rocks, occur in parautochthonous target rocks, mainly on the outlying islands of the Slate Islands archipelago. A few contain fragmented and disrupted, coarse-grained, polymictic clastic matrix breccia dikes. This is an indication that at least some of these monomictic breccias formed late in the impact process and that they are probably related to a late crater modification stage. A small number of relatively large occurrences of glass-poor, suevitic breccias occur at the flanks of the central uplift and along the inner flank of the outer ring of the Slate Islands complex crater. A coarse, glass-free, allogenic breccia, containing shatter-coned fragments derived from Proterozoic target rocks (upper target strata), observed at two locations may be analogous to the 'Bunt Breccia' of the Ries crater in Germany. At one of these locations this breccia lies close to a crater suevite deposit. At the other, it overlies parautochthonous, monomictic breccia. The State Islands impact breccias are superbly exposed, much better than breccias in most other terrestrial impact structures. Observations, including those indicative of multiple and and sequential processes, provide insight on how impact breccias form and how they relate to the various phases of the impact process. Eventually they will lead to an improved understanding of planetary impact processes.

Dressler, B. O.↗