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 163 records · Page 9

Computationally efficient modeling and simulation of large scale systems

A method of simulating operation of a VLSI interconnect structure having capacitive and inductive coupling between nodes thereof. A matrix X and a matrix Y containing different combinations of passive circuit element values for the interconnect structure are obtained where the element values for each matrix include inductance L and inverse capacitance P. An adjacency matrix A associated with the interconnect structure is obtained. Numerical integration is used to solve first and second equations, each including as a factor the product of the inverse matrix X.sup.1 and at least one other matrix, with first equation including X.sup.1Y, X.sup.1A, and X.sup.1P, and the second equation including X.sup.1A and X.sup.1P.

Jain, Jitesh↗

Computationally efficient modeling and simulation of large scale systems

A method of simulating operation of a VLSI interconnect structure having capacitive and inductive coupling between nodes thereof. A matrix X and a matrix Y containing different combinations of passive circuit element values for the interconnect structure are obtained where the element values for each matrix include inductance L and inverse capacitance P. An adjacency matrix A associated with the interconnect structure is obtained. Numerical integration is used to solve first and second equations, each including as a factor the product of the inverse matrix X.sup.-1 and at least one other matrix, with first equation including X.sup.-1Y, X.sup.-1A, and X.sup.-1P, and the second equation including X.sup.-1A and X.sup.-1P.

Jain, Jitesh↗

Fibronectin is a survival factor for differentiated osteoblasts

The skeletal extracellular matrix produced by osteoblasts contains the glycoprotein fibronectin, which regulates the adhesion, differentiation and function of various adherent cells. Interactions with fibronectin are required for osteoblast differentiation in vitro, since fibronectin antagonists added to cultures of immature fetal calvarial osteoblasts inhibit their progressive differentiation. To determine if fibronectin plays a unique role in fully differentiated osteoblasts, cultures that had already formed mineralized nodules in vitro were treated with fibronectin antagonists. Fibronectin antibodies caused >95% of the cells in the mature cultures to display characteristic features of apoptosis (nuclear condensation, apoptotic body formation, DNA laddering) within 24 hours. Cells appeared to acquire sensitivity to fibronectin antibody-induced apoptosis as a consequence of differentiation, since antibodies failed to kill immature cells and the first cells killed were those associated with mature nodules. Intact plasma fibronectin, as well as fragments corresponding to the amino-terminal, cell-binding, and carboxy-terminal domains of fibronectin, independently induced apoptosis of mature (day-13), but not immature (day-4), osteoblasts. Finally, transforming growth factor-beta1 partially protected cells from the apoptotic effects of fibronectin antagonists. Thus, in the course of maturation cultured osteoblasts switch from depending on fibronectin for differentiation to depending on fibronectin for survival. These data suggest that fibronectin, together with transforming growth factor-beta1, may affect bone formation, in part by regulating the survival of osteoblasts.

NASA Center ARC↗

Parallelization of the Physical-Space Statistical Analysis System (PSAS)

Atmospheric data assimilation is a method of combining observations with model forecasts to produce a more accurate description of the atmosphere than the observations or forecast alone can provide. Data assimilation plays an increasingly important role in the study of climate and atmospheric chemistry. The NASA Data Assimilation Office (DAO) has developed the Goddard Earth Observing System Data Assimilation System (GEOS DAS) to create assimilated datasets. The core computational components of the GEOS DAS include the GEOS General Circulation Model (GCM) and the Physical-space Statistical Analysis System (PSAS). The need for timely validation of scientific enhancements to the data assimilation system poses computational demands that are best met by distributed parallel software. PSAS is implemented in Fortran 90 using object-based design principles. The analysis portions of the code solve two equations. The first of these is the "innovation" equation, which is solved on the unstructured observation grid using a preconditioned conjugate gradient (CG) method. The "analysis" equation is a transformation from the observation grid back to a structured grid, and is solved by a direct matrix-vector multiplication. Use of a factored-operator formulation reduces the computational complexity of both the CG solver and the matrix-vector multiplication, rendering the matrix-vector multiplications as a successive product of operators on a vector. Sparsity is introduced to these operators by partitioning the observations using an icosahedral decomposition scheme. PSAS builds a large (approx. 128MB) run-time database of parameters used in the calculation of these operators. Implementing a message passing parallel computing paradigm into an existing yet developing computational system as complex as PSAS is nontrivial. One of the technical challenges is balancing the requirements for computational reproducibility with the need for high performance. The problem of computational reproducibility is well known in the parallel computing community. It is a requirement that the parallel code perform calculations in a fashion that will yield identical results on different configurations of processing elements on the same platform. In some cases this problem can be solved by sacrificing performance. Meeting this requirement and still achieving high performance is very difficult. Topics to be discussed include: current PSAS design and parallelization strategy; reproducibility issues; load balance vs. database memory demands, possible solutions to these problems.

Larson, J. W.↗

Micromechanics-based Modeling of Laminated SiC/SiC Ceramic Matrix Composites

The behavior and response of ceramic matrix composites (CMCs), in particular silicon carbide fiber reinforced silicon carbide matrix (SiC/SiC), is affected by many factors such as variation of fiber volume fraction, residual stresses resulting from processing of the composites at high temperature, random microstructures, and the presence of matrix flaws (e.g., voids, pores, cracks etc.) as well as general material nonlinearity and heterogeneity that occurs randomly in a composite. Residual stresses arising from the phase change of constituents are evaluated in this paper and it is shown that they do influence composite strength and need to be properly accounted for. Additionally, the microstructures (location of fiber centers, coating thickness etc.) of advanced CMCs are usually disordered (or random) and fiber diameter and strength typically have a distribution. They rarely resemble the ordered fiber packing (square, rectangular, or hexagonal) that is generally assumed in micromechanics-based models with periodic boundary conditions for computational expediency. These issues raise the question of how should one model such systems effectively? Can an ordered hexagonal packed repeating unit cell (RUC) accurately represent the random microstructure behavior? How many fibers need to be included to enable accurate representation? Clearly, the number of fibers within an RUC must be limited to insure a balance between accuracy and efficiency. NASA’s in-house micromechanics-based code MAC/GMC provides a framework to analyze such RUCs for the overall composite behavior and the FEAMAC computer code provides linkage of MAC/GMC to the commercial FEA code, ABAQUS. The appropriate level of discretization of the RUC as well as the analysis method employed, i.e., Generalized Method of Cells (GMC) or High Fidelity Generalized Method of Cells (HFGMC), is investigated in this paper in the context of a unidirectional as well as a cross-ply laminated CMC. Results including effective composite properties, proportional limit stress (an important design parameter) and fatigue are shown utilizing both GMC as well as HFGMC. Finally, a few multiscale analyses are performed on smooth bar test coupons as well as test coupons with features such as open-hole and double notches using FEAMAC. Best practices and guidance are provided to take these phenomena into account and keep a proper balance between fidelity (accuracy) and efficiency. Following these guidelines can account for important physics of the problem and provide significant advantages when performing large multiscale composite structural analyses.

Ceramic Matrix Composites↗

Primary igneous carbon in ureilites - Petrological implications

The ureilite meteorites are carbonaceous olivine-pyroxene achondrites. They typically contain up to 4 wt.% carbon (carbonaceous matrix) as graphite, diamond, and lonsdaleite. Shock degradation has effectively obliterated primary textures in the carbonaceous matrix of previously described ureilites, a factor that has hampered efforts to explain the origin of this material. In contrast, the Antarctic ureilite ALHA78019 displays perfectly preserved primary textures in the carbonaceous matrix characterized by euhedral graphite blades intergrown with Fe-Ni metal and sulfide (diamonds are absent). This petrographic feature suggests that most graphite in ureilites originated by crystallization from a C-rich metallic phase. Assuming that fO2 is controlled by C-CO-CO2 reactions, the compositions of silicates and metals in ureilites imply a two-stage redox history. The noble gases and rare earths of ureilites are discussed in light of this model.

Berkely, J. L.↗

Thermal expansion of composites: Methods and results

The factors controlling the dimensional stability of various components of large space structures were investigated. Cyclic, thermal and mechanical loading were identified as the primary controlling factors of the dimensional stability of cables. For organic matrix composites, such as graphite-epoxy, it was found that these factors include moisture desorption in the space environment, thermal expansion as the structure moves from the sunlight to shadow in its orbit, mechanical loading, and microyielding of the material caused by microcracking of the matrix material. The major focus was placed on the thermal expansion of composites and in particular the development and testing of a method for its measurement.

Bowles, D. E.↗

Using Temporal Fill Factor to Reduce Frame Reconstruction Rates

The newer active matrix display technologies such as TFT-LCD, DMD, PDP maintain their pixel values through the entire frame time, presenting a 100% temporal fill factor, in contrast to the duty cycle produced by the phosphor impulse response of the CRT. This sample-and-hold characteristic can be exploited to lower the displayed frame rate without affecting visual quality. The lower frame rate results in significantly lower transmission bandwidth, power, and cost.

Larimer, James↗

Optimal inputs for system identification

Identification criteria are presented for linear dynamic systems with and without process noise. With process noise, the state equations are replaced by the Kalman filter equations. If the identification performance index is expanded in a Taylor's series with respect to the parameters to be identified, then maximizing the weighting factor of the quadratic term with respect to the inputs will insure that an identification algorithm will converge more rapidly and to a more accurate result than with nonoptimal inputs. The expectation of this weighting factor is the Fisher information matrix, and its inverse is a lower bound for the covariance of the parameters. Direct and indirect methods of calculating the information matrix are presented for systems with and without process noise.

Reid, D. B.↗

Multirate sampled-data systems analysis via vector operators

The primary difficulties of both the time-domain switch decomposition method and the frequency-domain decomposition method are overcome by introducing certain matrix operators and performing spectral factorization of resulting matrices of polynomials in the z-transform variable. Topological operations of the switch-decomposition method are simplified. This new approach eliminates the need to solve a system of equations with rational polynomial coefficients such as arises in the frequency-decomposition method. The determination of a multirate sampled-data system's characteristic polynomial no longer requires the evaluation of a determinant of rational polynomial elements. New results on obtaining modified z-transforms from standard z-transforms at a faster rate and vice versa are presented.

Boykin, W. H.↗

An analysis of thermal response factors and how to reduce their computational time requirement

Te RESFAC2 version of the Thermal Response Factor Program (RESFAC) is the result of numerous modifications and additions to the original RESFAC. These modifications and additions have significantly reduced the program's computational time requirement. As a result of this work, the program is more efficient and its code is both readable and understandable. This report describes what a thermal response factor is; analyzes the original matrix algebra calculations and root finding techniques; presents a new root finding technique and streamlined matrix algebra; supplies ten validation cases and their results.

Wiese, M. R.↗

Random field estimation approach to multibody dynamics

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

Rodriguez, G.↗

Random field estimation approach to robot dynamics

The difference equations of Kalman filtering and smoothing recursively factor and invert the covariance of the output of a linear state-space system driven by a white-noise process. Here it is shown that similar recursive techniques factor and invert the inertia matrix of a multibody robot system. The random field models are based on the assumption that all of the inertial (D'Alembert) forces in the system are represented by a spatially distributed white-noise model. They are easier to describe than the models based on classical mechanics, which typically require extensive derivation and manipulation of equations of motion for complex mechanical systems. With the spatially random models, more primitive locally specified computations result in a global collective system behavior equivalent to that obtained with deterministic models. The primary goal of applying random field estimation is to provide a concise analytical foundation for solving robot control and motion planning problems.

Rodriguez, Guillermo↗

A methodology to predict damage initiation, damage growth and residual strength in titanium matrix composites

In this research, a methodology to predict damage initiation, damage growth, fatigue life, and residual strength in titanium matrix composites (TMC) is outlined. Emphasis was placed on micromechanics-based engineering approaches. Damage initiation was predicted using a local effective strain approach. A finite element analysis verified the prevailing assumptions made in the formulation of this model. Damage growth, namely, fiber-bridged matrix crack growth, was evaluated using a fiber bridging (FB) model which accounts for thermal residual stresses. This model combines continuum fracture mechanics and micromechanics analyses yielding stress-intensity factor solutions for fiber-bridged matrix cracks. It is assumed in the FB model that fibers in the wake of the matrix crack are idealized as a closure pressure, and an unknown constant frictional shear stress is assumed to act along the debond length of the bridging fibers. This frictional shear stress was used as a curve fitting parameter to the available experimental data. Fatigue life and post-fatigue residual strength were predicted based on the axial stress in the first intact 0 degree fiber calculated using the FB model and a three-dimensional finite element analysis.

Bakuckas, J. G., Jr.↗

Accounting for the photochemical variation in stratospheric NO2 in the SAGE III/ISS solar occultation retrieval

The Stratospheric Aerosol and Gas Experiment (SAGE) III has been operating on the International Space Station (ISS) since mid-2017. Nitrogen dioxide (NO2) number density profiles are routinely retrieved from SAGE III/ISS solar occultation measurements in the middle atmosphere. Although NO2 density varies throughout the day due to photochemistry, the standard SAGE NO2 retrieval algorithm neglects these variations along the instrument’s line of sight by assuming that the number density has a constant gradient within a given vertical layer of the atmosphere. This assumption will result in a retrieval bias for a species like NO2 that changes rapidly across the terminator. In this work we account for diurnal variations in retrievals of NO2 from the SAGE III/ISS measurements, and we determine the impact of this algorithm improvement on the resulting NO2 number densities. The first step in applying the diurnal correction is to use publicly available SAGE III/ISS products to convert the retrieved number density profiles to optical depth profiles. The retrieval is then re-performed with a new matrix that applies photochemical scale factors for each point along the line of sight according to the changing solar zenith angle. In general NO2 that is retrieved by accounting for these diurnal variations is more than 10% lower than the standard algorithm below 30 km. This effect is greatest in winter at high latitudes and generally greater for sunrise occultations than sunset. Comparisons with coincident profiles from the Optical Spectrograph and InfraRed Imager System (OSIRIS) show that NO2 from SAGE III/ISS is generally biased high; however the agreement improves by up to 20% in the mid-stratosphere when diurnal variations are accounted for in the retrieval. We conclude that diurnal variations along the SAGE III/ISS line of sight are an important term to consider for NO2 analyses at altitudes below 30 km.

Stratospheric NO2↗

Analysis of the Effects of Residual Strains and Defects on Skin/Stiffener Debonding using Decohesion Elements

Delamination is one of the predominant forms of failure in laminated composites especially when there is no reinforcement in the thickness direction. To develop composite structures that are more damage tolerant, it is necessary to understand how delamination develops and how it can affect the residual performance. A number of factors such as residual thermal strains, matrix curing shrinkage, and manufacturing defects affect how damage will grow in a composite structure. It is important to develop analysis methods that are computationally efficient that can account for all such factors. The objective of the current work is to apply a newly developed decohesion element to investigate the debond strength of skin/stiffener composite specimens. The process of initiation of delaminations and the propagation of delamination fronts is investigated. The numerical predictions are compared with published experimental results.

Davila, Carlos G.↗