Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix factorization”

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 433 records · Page 24

Signatures of the attractive interaction in spin spectra of one-dimensional cuprate chains

Identifying the minimal model for cuprates is crucial for explaining the high- T c pairing mechanism. Recent photoemission experiments have suggested a significant near-neighbor attractive interaction V in cuprate chains, favoring pairing instability. To determine its strength, we systematically investigate the dynamical spin structure factors S ( q , ω ) using the density matrix renormalization group. Our analysis quantitatively reveals a notable softening in the two-spinon continuum, particularly evident in the intense spectrum at large momentum. This softening is primarily driven by the renormalization of the superexchange interaction, as determined by a comparison with the slave-boson theory. We also demonstrate the feasibility of detecting this spectral shift in thin-film samples using resonant inelastic x-ray scattering. Therefore, this provides a distinctive fingerprint for the attractive interaction, motivating future experiments to unveil essential ingredients in cuprates. Published by the American Physical Society 2024

Shen, Zecheng↗

The TPR domain of PgaA is a multifunctional scaffold that binds PNAG and modulates PgaB-dependent polymer processing

The synthesis of exopolysaccharides as biofilm matrix components by pathogens is a crucial factor for chronic infections and antibiotic resistance. Many periplasmic proteins involved in polymer processing and secretion in Gram-negative synthase dependent exopolysaccharide biosynthetic systems have been individually characterized. The operons responsible for the production of PNAG, alginate, cellulose and the Pel polysaccharide each contain a gene that encodes an outer membrane associated tetratricopeptide repeat (TPR) domain containing protein. While the TPR domain has been shown to bind other periplasmic proteins, the functional consequences of these interactions for the polymer remain poorly understood. Herein, we show that the C-terminal TPR region of PgaA interacts with the de-N-acetylase domain of PgaB, and increases its deacetylase activity. Additionally, we found that when the two proteins form a complex, the glycoside hydrolase activity of PgaB is also increased. To better understand structure-function relationships we determined the crystal structure of a stable TPR module, which has a conserved groove formed by three repeat motifs. Tryptophan quenching, mass spectrometry analysis and molecular dynamics simulation studies suggest that the crystallized TPR module can bind PNAG/dPNAG via its electronegative groove on the concave surface, and potentially guide the polymer through the periplasm towards the porin for export. Our results suggest a scaffolding role for the TPR domain that combines PNAG/dPNAG translocation with the modulation of its chemical structure by PgaB.

59 BASIC BIOLOGICAL SCIENCES↗

Analysis of Covariance Intersection For Triangulation

In this document I will discuss different implementations of Covariance Intersection (CI) for object triangulation as well as the robustness of CI. For CI methods we will compare the performance of different methods and discuss edge cases which must be considered. For robustness we will focus on the impact of removing sensors on the final fused estimate. Here we look at factors which influence the final fused covariance matrix.

97 MATHEMATICS AND COMPUTING↗

Binning is Sinning: Redemption for Hubble Diagram Using Photometrically Classified Type Ia Supernovae

Bayesian Estimation Applied to Multiple Species (BEAMS) is implemented in the BEAMS with Bias Corrections (BBC) framework to produce a redshift-binned Hubble diagram (HD) for Type Ia supernovae (SNe Ia). BBC corrects for selection effects and non–SN Ia contamination, and systematic uncertainties are described by a covariance matrix with dimension matching the number of BBC redshift bins. For spectroscopically confirmed SN Ia samples, a recent "Binning is Sinning" article showed that an unbinned HD and covariance matrix reduces the systematic uncertainty by a factor of ~1.5 compared to the binned approach. Here we extend their analysis to obtain an unbinned HD for a photometrically identified sample processed with BBC. To test this new method, we simulate and analyze 50 samples corresponding to the Dark Energy Survey (DES) with a low-redshift anchor; the simulation includes SNe Ia, and contaminants from core-collapse SNe and peculiar SNe Ia. The analysis includes systematic uncertainties for calibration and measures the dark energy equation of state parameter (w). Compared to a redshift-binned HD, the unbinned HD with nearly 2000 events results in a smaller systematic uncertainty, in qualitative agreement with BHS21, and averaging results among the 50 samples we find no evidence for a w-bias. To reduce computation time for fitting an unbinned HD with large samples, we propose an HD-rebinning method that defines the HD in bins of redshift, color, and stretch; the rebinned HD results in similar uncertainty as the unbinned case, and shows no evidence for a w-bias.

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis and fifteen-year projection of the market for LANDSAT data

The potential market for LANDSAT products through the 1990's was determined. Results are presented in a matrix format. Improved resolution is a major factor in the marketability of LANDSAT data, the 10 meter resolution (projected for 1995) having a significant impact on the federal, private, and international users, and on the agricultural, minerals, and national defense applications. Data delivery time and competition from the French remote sensing system are considered.

Source record↗

Fracture toughness of fibrous composite materials

Laminates with various proportions of 0 deg, 45 deg, and 90 deg plies were fabricated from T300/5208 and T300/BP-907 graphite/epoxy prepreg tape material. The fracture toughness of each laminate orientation or lay-up was determined by testing center-cracked specimens, and it was also predicted with the general fracture-toughness parameter. The predictions were good except when crack-tip splitting was large, at which time the toughness and strengths tended to be underpredicted. By using predictions, a parametric study was also made of factors that influence fracture toughness. Fiber and matrix properties as well as lay-up were investigated. Without crack-tip splitting, fracture toughness increases in proportion to fiber strength and fiber volume fraction, increases linearly with E(22)/E(11), is largest when the modulus for non-0 deg fibers is greater than that of 0 deg fibers, and is smallest for 0(m)/90(p)(s) lay-ups. (The E(11) and E(22) are Young's moduli of the lamina parallel to and normal to the direction of the fibers, respectively). For a given proportion of 0 deg plies, the most notch-sensitive lay-ups are 0(m)/90(p)(s) and the least sensitive are 0(m)/45(n)(s) and alpha(s). Notch sensitivity increases with the proportion of 0 deg plies and decreases with alpha. Strong, tough matrix materials, which inhibit crack-tip splitting, generally lead to minimum fracture toughness.

Poe, C. C., Jr.↗

Composite interlaminar fracture - Effect of matrix fracture energy

A major result of the use of the double cantilever beam test to provide a measure of the interlaminar fracture energy in composites subjected to simple Mode I-type loading is the observation that increasing the toughness of the matrix resin by as much as a factor of 20 produced a major (approximately four-fold) increase in interlaminar fracture energy. The data presently analyzed show that, with brittle polymers, resin toughness is fully transferred to the composite, while in the case of tougher polymers the resin toughness is only partially transferred due to the fibers' restriction of the crack tip deformation zone in the polymer, together with their changing of local stress field details. Factors which tend to increase interlaminar toughness are fiber nesting and bridging, as well as fiber breakage and pull-out during crack growth.

Hunston, D. L.↗

Parallel solution of closely coupled systems

An odd-even permutation and a nested dissection technique were used to circumvent the strong seriality of a system of closely coupled equations. The effect of transforming the n x n Hermitian definite positive matrix coefficient on the topology of Cholesky factors is discussed. A series of directed graphs is constructed in order to show the computational steps required for the odd-even permutation. Numerical expressions for the speed-up and efficiency of parallel N-processing techniques and sequential processing by a single computer are derived. Similar expressions are derived for the case of insufficient processing capacity. The application of the odd-even permutation to the ensemble class of computer architectures is demonstrated.

Utku, S.↗

Assessment of space environment induced microdamage in toughened composite materials

The effects of simulated space environments on the microdamage in a series of commercially available toughened matrix composite systems was determined. Low-earth orbit (LEO) exposures were simulated by thermal cycling; geosynchronous orbit (GEO) exposures were simulated by electron irradiation plus thermal cycling. Material response was characterized by assessing the induced microcracking and its influence on chemical and mechanical property changes. All materials, including several advanced, tough thermoplastics microcracked when exposed to the simulated LEO environment except a 177 C cured single phase toughened epoxy composite. The GEO simulated environment produced microdamage in all materials. The results suggest that increased matrix toughness may not be the overriding factor leading to improved durability in the space environment.

Sykes, George F.↗

Thomas precession, Wigner rotations and gauge transformations

The exact Lorentz kinematics of the Thomas precession is discussed in terms of Wigner's O(3)-like little group which describes rotations in the Lorentz frame in which the particle is at rest. A Lorentz-covariant form for the Thomas factor is derived. It is shown that this factor is a Lorentz-boosted rotation matrix, which becomes a gauge transformation in the infinite-momentum or zero-mass limit.

Han, D.↗

Modeling, model simplification and stability robustness with aeroelastic vehicles

The linearization and simplification of a nonlinear, literal model for flexible aircraft is highlighted. Areas of model fidelity that are critical if the model is to be used for control system synthesis are developed, and several simplification techniques that can deliver the necessary model fidelity are discussed. These techniques include both numerical and analytical approaches. An analytical approach, based on first-order sensitivity theory is shown to lead not only to excellent numerical results, but also to closed-form analytical expressions for key system dynamic properties such as the pole/zero factors of the vehicle transfer-function matrix. The analytical results are expressed in terms of vehicle mass properties, vibrational characteristics, and rigid-body and aeroelastic stability derivatives, thus leading to the underlying causes for critical dynamic characteristics.

Schmidt, David K.↗

Flexible aircraft dynamic modeling for dynamic analysis and control synthesis

The linearization and simplification of a nonlinear, literal model for flexible aircraft is highlighted. Areas of model fidelity that are critical if the model is to be used for control system synthesis are developed and several simplification techniques that can deliver the necessary model fidelity are discussed. These techniques include both numerical and analytical approaches. An analytical approach, based on first-order sensitivity theory is shown to lead not only to excellent numerical results, but also to closed-form analytical expressions for key system dynamic properties such as the pole/zero factors of the vehicle transfer-function matrix. The analytical results are expressed in terms of vehicle mass properties, vibrational characteristics, and rigid-body and aeroelastic stability derivatives, thus leading to the underlying causes for critical dynamic characteristics.

Schmidt, David K.↗

Efficient computer algebra algorithms for polynomial matrices in control design

The theory of polynomial matrices plays a key role in the design and analysis of multi-input multi-output control and communications systems using frequency domain methods. Examples include coprime factorizations of transfer functions, cannonical realizations from matrix fraction descriptions, and the transfer function design of feedback compensators. Typically, such problems abstract in a natural way to the need to solve systems of Diophantine equations or systems of linear equations over polynomials. These and other problems involving polynomial matrices can in turn be reduced to polynomial matrix triangularization procedures, a result which is not surprising given the importance of matrix triangularization techniques in numerical linear algebra. Matrices with entries from a field and Gaussian elimination play a fundamental role in understanding the triangularization process. In the case of polynomial matrices, matrices with entries from a ring for which Gaussian elimination is not defined and triangularization is accomplished by what is quite properly called Euclidean elimination. Unfortunately, the numerical stability and sensitivity issues which accompany floating point approaches to Euclidean elimination are not very well understood. New algorithms are presented which circumvent entirely such numerical issues through the use of exact, symbolic methods in computer algebra. The use of such error-free algorithms guarantees that the results are accurate to within the precision of the model data--the best that can be hoped for. Care must be taken in the design of such algorithms due to the phenomenon of intermediate expressions swell.

Baras, J. S.↗

Numerical and literal aeroelastic-vehicle-model reduction for feedback control synthesis

The simplification of a high-order, literal model for large flexible aircraft is discussed. Areas of model fidelity that are critical if the model is to be used for control law synthesis are presented. Several simplification techniques, some new and some widely available, that can deliver the necessary model fidelity are presented and applied to a model from the literature. The techniques include both numerical and analytical approaches. An analytical approach, based on first-order sensitivity theory, is shown to lead not only to excellent numerical results, but also to closed-form analytical expressions for key system dynamic properties such as the pole/zero factors of the vehicle transfer-function matrix. The analytical results are expressed in terms of vehicle vibrational characteristics and rigid-body and aeroelastic stability derivatives, thus providing insight into the underlying causes for critical dynamic characteristics.

Newman, Brett↗

Parallel CFD Algorithms for Aerodynamical Flow Solvers on Unstructured Meshes

The Advisory Group for Aerospace Research and Development (AGARD) has requested my participation in the lecture series entitled Parallel Computing in Computational Fluid Dynamics to be held at the von Karman Institute in Brussels, Belgium on May 15-19, 1995. In addition, a request has been made from the US Coordinator for AGARD at the Pentagon for NASA Ames to hold a repetition of the lecture series on October 16-20, 1995. I have been asked to be a local coordinator for the Ames event. All AGARD lecture series events have attendance limited to NATO allied countries. A brief of the lecture series is provided in the attached enclosure. Specifically, I have been asked to give two lectures of approximately 75 minutes each on the subject of parallel solution techniques for the fluid flow equations on unstructured meshes. The title of my lectures is "Parallel CFD Algorithms for Aerodynamical Flow Solvers on Unstructured Meshes" (Parts I-II). The contents of these lectures will be largely review in nature and will draw upon previously published work in this area. Topics of my lectures will include: (1) Mesh partitioning algorithms. Recursive techniques based on coordinate bisection, Cuthill-McKee level structures, and spectral bisection. (2) Newton's method for large scale CFD problems. Size and complexity estimates for Newton's method, modifications for insuring global convergence. (3) Techniques for constructing the Jacobian matrix. Analytic and numerical techniques for Jacobian matrix-vector products, constructing the transposed matrix, extensions to optimization and homotopy theories. (4) Iterative solution algorithms. Practical experience with GIVIRES and BICG-STAB matrix solvers. (5) Parallel matrix preconditioning. Incomplete Lower-Upper (ILU) factorization, domain-decomposed ILU, approximate Schur complement strategies.

Barth, Timothy J.↗

Maximally Informative Statistics for Localization and Mapping

This paper presents an algorithm for localization and mapping for a mobile robot using monocular vision and odometry as its means of sensing. The approach uses the Variable State Dimension filtering (VSDF) framework to combine aspects of Extended Kalman filtering and nonlinear batch optimization. This paper describes two primary improvements to the VSDF. The first is to use an interpolation scheme based on Gaussian quadrature to linearize measurements rather than relying on analytic Jacobians. The second is to replace the inverse covariance matrix in the VSDF with its Cholesky factor to improve the computational complexity. Results of applying the filter to the problem of localization and mapping with omnidirectional vision are presented.

Deans, Matthew C.↗

Mapping Successional Stages in a Wet Tropical Forest Using Landsat ETM+ and Forest Inventory Data

In this study, we test whether an existing classification technique based on the integration of Landsat ETM+ and forest inventory data enables detailed characterization of successional stages in a wet tropical forest site. The specific objectives were: (1) to map forest age classes across the La Selva Biological Station in Costa Rica; and (2) to quantify uncertainties in the proposed approach in relation to field data and existing vegetation maps. Although significant relationships between vegetation height entropy (a surrogate for forest age) and ETM+ data were detected, the classification scheme tested in this study was not suitable for characterizing spatial variation in age at La Selva, as evidenced by the error matrix and the low Kappa coefficient (12.9%). Factors affecting the performance of the classification at this particular study site include the smooth transition in vegetation structure between intermediate and advanced successional stages, and the low sensitivity of NDVI to variations in vertical structure at high biomass levels.

forest structure↗

Scaled ILU Smoothers for Navier-Stokes Pressure Projection

Incomplete LU (ILU) smoothers are effective in the algebraic multigrid (AMG) V-cycle for reducing high-frequency components of the error. However, the requisite direct triangular solves are comparatively slow on GPUs. Previous work has demonstrated the advantages of Jacobi iteration as an alternative to direct solution of these systems. Depending on the threshold and fill-level parameters chosen, the factors can be highly nonnormal and Jacobi is unlikely to converge in a low number of iterations. We demonstrate that row scaling can reduce the departure from normality, allowing us to replace the inherently sequential solve with a rapidly converging Richardson iteration. There are several advantages beyond the lower compute time. Scaling is performed locally for a diagonal block of the global matrix because it is applied directly to the factor. Further, an ILUT Schur complement smoother maintains a constant GMRES iteration count as the number of MPI ranks increases, and thus parallel strong-scaling is improved. Our algorithms have been incorporated into hypre, and we demonstrate improved time to solution for linear systems arising in the Nalu-Wind and PeleLM pressure solvers. For large problem sizes, GMRES+AMG executes at least five times faster when using iterative triangular solves compared with direct solves on massively parallel GPUs.

algebraic multigrid↗