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 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Experiment K-6-06. Morphometric and EM analyses of tibial epiphyseal plates from Cosmos 1887 rats

Light and electron microscopy studies were carried out on decalcified tibial epiphyseal plates of rats flown aboard Cosmos 1887 (12.5d flight plus 53.5h recovery). Analysis of variance showed that the proliferative zone of flight animals was significantly higher than that of synchronous controls, while the hypertrophic/calcification zone was significantly reduced. Flight animals had more cells than synchronous controls in the proliferative zone, and less in the hypertrophic/calcification region. The total number of cells, however, was significantly higher in flight animals. No differences were found for perimeter or shape factor of growth plates, but area was significantly lower in flight animals in comparison to synchronous controls. Collagen fibrils in flight animals were shorter and wider than in synchronous controls. The time required for a cell to cycle through the growth plate is 2 to 3 days, so most of the cells and matrix present were formed after the animals had returned to 1 g, and probably represent stages of recovery from microgravity exposure, which in itself is an interesting question.

Duke, P. J.↗
NASA NTRS1990

Computation of the factorized error covariance of the difference between correlated estimators

A state estimation problem where some of the measurements may be common to two or more data sets is considered. Two approaches for computing the error covariance of the difference between filtered estimates (for each data set) are discussed. The first algorithm is based on postprocessing of the Kalman gain profiles of two correlated estimators. It uses UD factors of the covariance of the relative error. The second algorithm uses a square root information filter applied to relative error analysis. In the absence of process noise, the square root information filter is computationally more efficient and more flexible than the Kalman gain (covariance update) method. Both the algorithms (covariance and information matrix based) are applied to a Venus orbiter simulation, and their performances are compared.

Wolff, Peter J.↗
NASA NTRS1991

Controllable surface-plasmon resonance in engineered nanometer epitaxial silicide particles embedded in silicon

Epitaxial CoSi2 particles in a single-crystal silicon matrix are grown by molecular-beam epitaxy using a technique that allows nanometer control over particle size in three dimensions. These composite layers exhibit resonant absorption predicted by effective-medium theory. Selection of the height and diameter of disklike particles through a choice of growth conditions allows tailoring of the depolarization factor and hence of the surface-plasmon resonance energy. Resonant absorption from 0.49 to 1.04 eV (2.5 to 1.2 micron) is demonstrated and shown to agree well with values predicted by the Garnett (1904, 1906) theory using the bulk dielectric constants for CoSi2 and Si.

Fathauer, R. W.↗
NASA NTRS1995

A Hybrid Parallel Preconditioning Algorithm For CFD

A new hybrid preconditioning algorithm will be presented which combines the favorable attributes of incomplete lower-upper (ILU) factorization with the favorable attributes of the approximate inverse method recently advocated by numerous researchers. The quality of the preconditioner is adjustable and can be increased at the cost of additional computation while at the same time the storage required is roughly constant and approximately equal to the storage required for the original matrix. In addition, the preconditioning algorithm suggests an efficient and natural parallel implementation with reduced communication. Sample calculations will be presented for the numerical solution of multi-dimensional advection-diffusion equations. The matrix solver has also been embedded into a Newton algorithm for solving the nonlinear Euler and Navier-Stokes equations governing compressible flow. The full paper will show numerous examples in CFD to demonstrate the efficiency and robustness of the method.

Barth,Timothy J.↗
DOE OSTI2020

Matrix-free subcell residual distribution for Bernstein finite elements: Monolithic limiting

This paper is focused on the aspects of limiting in residual distribution (RD) schemes for high-order finite element approximations to advection problems. Both continuous and discontinuous Galerkin methods are considered in this work. Discrete maximum principles are enforced using algebraic manipulations of element contributions to the global nonlinear system. The required modifications can be carried out without calculating the element matrices and assembling their global counterparts. The components of element vectors associated with the standard Galerkin discretization are manipulated directly using localized subcell weights to achieve optimal accuracy. Low-order nonlinear RD schemes of this kind were originally developed to calculate local extremum diminishing predictors for flux-corrected transport (FCT) algorithms. In the present paper, we incorporate limiters directly into the residual distribution procedure, which makes it applicable to stationary problems and leads to well-posed nonlinear discrete problems. To circumvent the second-order accuracy barrier, the correction factors of monolithic limiting approaches and FCT schemes are adjusted using smoothness sensors based on second derivatives. Finally, the convergence behavior of presented methods is illustrated by numerical studies for two-dimensional test problems.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1975

The inclusion problem with a crack crossing the boundary

A solution is given to the problems of a crack in an elastic inclusion with one or both ends approaching and terminating at the interface, of two collinear cracks (one in the inclusion and one in the matrix), and of a crack crossing the interface. The problems are formulated in terms of a system of singular integral equations. In the second case, the dominant parts of the kernels become generalized Cauchy kernels giving rise to stress singularities of powers other than one over the square root. For this case stress intensity factors are defined and some detailed results are presented for various crack-inclusion geometries and material combinations.

Erdogan, F.↗
DOE OSTI2025

Convergence Analysis of the Alternating Anderson–Picard Method for Nonlinear Fixed-Point Problems

Anderson acceleration (AA) has been widely used to solve nonlinear fixed-point problems due to its rapid convergence. This work focuses on a variant of AA in which multiple Picard iterations are performed between each AA step, referred to as the Alternating Anderson–Picard (AAP) method. Furthermore, despite introducing more “slow” Picard iterations, this method has been shown to be efficient and even more robust in both linear and nonlinear cases. However, there is a lack of theoretical analysis for AAP in the nonlinear case. In this paper, we address this gap by establishing the equivalence between AAP and a multisecant-GMRES method that uses GMRES to solve a multisecant linear system at each iteration. From this perspective, we show that AAP “converges” to the Newton-GMRES method. Specifically, as the residual approaches zero, the multisecant matrix, the approximate Jacobian inverse, the search direction, and the optimization gain of AAP converge to their counterparts in the Newton-GMRES method. These connections provide insights for analyzing the asymptotic convergence properties of AAP. Consequently, we show that AAP is locally 𝑞-linear convergent and provide an upper bound for the convergence factor of AAP. To validate the theoretical results, numerical examples are provided.

Anderson acceleration↗
NASA NTRS1990

High-performance equation solvers and their impact on finite element analysis

The role of equation solvers in modern structural analysis software is described. Direct and iterative equation solvers which exploit vectorization on modern high-performance computer systems are described and compared. The direct solvers are two Cholesky factorization methods. The first method utilizes a novel variable-band data storage format to achieve very high computation rates and the second method uses a sparse data storage format designed to reduce the number of operations. The iterative solvers are preconditioned conjugate gradient methods. Two different preconditioners are included; the first uses a diagonal matrix storage scheme to achieve high computation rates and the second requires a sparse data storage scheme and converges to the solution in fewer iterations that the first. The impact of using all of the equation solvers in a common structural analysis software system is demonstrated by solving several representative structural analysis problems.

Poole, Eugene L.↗
NASA NTRS1992

High-performance equation solvers and their impact on finite element analysis

The role of equation solvers in modern structural analysis software is described. Direct and iterative equation solvers which exploit vectorization on modern high-performance computer systems are described and compared. The direct solvers are two Cholesky factorization methods. The first method utilizes a novel variable-band data storage format to achieve very high computation rates and the second method uses a sparse data storage format designed to reduce the number od operations. The iterative solvers are preconditioned conjugate gradient methods. Two different preconditioners are included; the first uses a diagonal matrix storage scheme to achieve high computation rates and the second requires a sparse data storage scheme and converges to the solution in fewer iterations that the first. The impact of using all of the equation solvers in a common structural analysis software system is demonstrated by solving several representative structural analysis problems.

Poole, Eugene L.↗
NASA NTRS2004

CMC Property Variability and Life Prediction Methods for Turbine Engine Component Application

The ever increasing need for lower density and higher temperature-capable materials for aircraft engines has led to the development of Ceramic Matrix Composites (CMCs). Today's aircraft engines operate with >3000"F gas temperatures at the entrance to the turbine section, but unless heavily cooled, metallic components cannot operate above approx.2000 F. CMCs attempt to push component capability to nearly 2700 F with much less cooling, which can help improve engine efficiency and performance in terms of better fuel efficiency, higher thrust, and reduced emissions. The NASA Glenn Research Center has been researching the benefits of the SiC/SiC CMC for engine applications. A CMC is made up of a matrix material, fibers, and an interphase, which is a protective coating over the fibers. There are several methods or architectures in which the orientation of the fibers can be manipulated to achieve a particular material property objective as well as a particular component geometric shape and size. The required shape manipulation can be a limiting factor in the design and performance of the component if there is a lack of bending capability of the fiber as making the fiber more flexible typically sacrifices strength and other fiber properties. Various analysis codes are available (pcGINA, CEMCAN) that can predict the effective Young's Moduli, thermal conductivities, coefficients of thermal expansion (CTE), and various other properties of a CMC. There are also various analysis codes (NASAlife) that can be used to predict the life of CMCs under expected engine service conditions. The objective of this summer study is to utilize and optimize these codes for examining the tradeoffs between CMC properties and the complex fiber architectures that will be needed for several different component designs. For example, for the pcGINA code, there are six variations of architecture available. Depending on which architecture is analyzed, the user is able to specify the fiber tow size, tow spacing, weave parameter, and angle of orientation of fibers. By holding the volume fraction of the fibers constant, variations in tow spacing can be explored for different architectures. The CMC material properties are usually calculated assuming the component is manufactured perfectly. However, this is typically not the case so that a quantification of the material property variability is needed to account for processing and/or manufacturing imperfections. The overall inputs and outputs are presented using a regression software to rapidly investigate the tradeoffs associated with fiber architecture, material properties, and ultimately cost. This information is then propagated through lifing models and Larson-Miller data to assess timehemperature-dependent CMC strength. In addition, a first order cost estimation will be quantified from a current qualitative perspective. This cost estimation includes the manufacturing challenges, such as tooling, as well as the component cost for a particular application. Ultimately, a cost to performance ratio should be established that compares the effectiveness of CMCs to their current rival, nickel superalloys.

Cheplak, Matthew L.↗
NASA NTRS2008

Reliability and Confidence Interval Analysis of a CMC Turbine Stator Vane

High temperature ceramic matrix composites (CMC) are being explored as viable candidate materials for hot section gas turbine components. These advanced composites can potentially lead to reduced weight, enable higher operating temperatures requiring less cooling and thus leading to increased engine efficiencies. However, these materials are brittle and show degradation with time at high operating temperatures due to creep as well as cyclic mechanical and thermal loads. In addition, these materials are heterogeneous in their make-up and various factors affect their properties in a specific design environment. Most of these advanced composites involve two- and three-dimensional fiber architectures and require a complex multi-step high temperature processing. Since there are uncertainties associated with each of these in addition to the variability in the constituent material properties, the observed behavior of composite materials exhibits scatter. Traditional material failure analyses employing a deterministic approach, where failure is assumed to occur when some allowable stress level or equivalent stress is exceeded, are not adequate for brittle material component design. Such phenomenological failure theories are reasonably successful when applied to ductile materials such as metals. Analysis of failure in structural components is governed by the observed scatter in strength, stiffness and loading conditions. In such situations, statistical design approaches must be used. Accounting for these phenomena requires a change in philosophy on the design engineer s part that leads to a reduced focus on the use of safety factors in favor of reliability analyses. The reliability approach demands that the design engineer must tolerate a finite risk of unacceptable performance. This risk of unacceptable performance is identified as a component's probability of failure (or alternatively, component reliability). The primary concern of the engineer is minimizing this risk in an economical manner. The methods to accurately determine the service life of an engine component with associated variability have become increasingly difficult. This results, in part, from the complex missions which are now routinely considered during the design process. These missions include large variations of multi-axial stresses and temperatures experienced by critical engine parts. There is a need for a convenient design tool that can accommodate various loading conditions induced by engine operating environments, and material data with their associated uncertainties to estimate the minimum predicted life of a structural component. A probabilistic composite micromechanics technique in combination with woven composite micromechanics, structural analysis and Fast Probability Integration (FPI) techniques has been used to evaluate the maximum stress and its probabilistic distribution in a CMC turbine stator vane. Furthermore, input variables causing scatter are identified and ranked based upon their sensitivity magnitude. Since the measured data for the ceramic matrix composite properties is very limited, obtaining a probabilistic distribution with their corresponding parameters is difficult. In case of limited data, confidence bounds are essential to quantify the uncertainty associated with the distribution. Usually 90 and 95% confidence intervals are computed for material properties. Failure properties are then computed with the confidence bounds. Best estimates and the confidence bounds on the best estimate of the cumulative probability function for R-S (strength - stress) are plotted. The methodologies and the results from these analyses will be discussed in the presentation.

Murthy, Pappu L. N.↗
DOE OSTI2021

Cellular mRNA triggers structural transformation of Ebola virus matrix protein VP40 to its essential regulatory form

The Ebola virus matrix protein VP40 forms distinct structures linked to distinct functions in the virus life cycle. Dimeric VP40 is a structural protein associated with virus assembly, while octameric, ring-shaped VP40 is associated with transcriptional control. In this study, we show that suitable nucleic acid is sufficient to trigger a dynamic transformation of VP40 dimer into the octameric ring. Deep sequencing reveals a binding preference of the VP40 ring for the 3' untranslated region of cellular mRNA and a guanine- and adenine-rich binding motif. Complementary analyses of the nucleic-acid-induced VP40 ring by native mass spectrometry, electron microscopy, and X-ray crystal structures at 1.8 and 1.4 Å resolution reveal the stoichiometry of RNA binding, as well as an interface involving a key guanine nucleotide. The host factor-induced structural transformation of protein structure in response to specific RNA triggers in the Ebola virus life cycle presents unique opportunities for therapeutic inhibition.

3′ UTR↗
DOE OSTI2024

Defect-engineered monolayer MoS 2 with enhanced memristive and synaptic functionality for neuromorphic computing

Two-dimensional transition metal dichalcogenides (TMDs)-based memristors are promising candidates for realizing artificial synapses in next-generation computing. However, practical implementation faces several challenges, such as high non-linearity and asymmetry in synaptic weight updates, limited dynamic range, and cycle-to-cycle variability. Here, utilizing optimal-power argon plasma treatment, we significantly enhance the performance matrix of memristors fabricated from monolayer MoS 2 . Our approach not only improves linearity and symmetry in synaptic weight updates but also increases the number of available synaptic weight updates and enhances Spike-Time Dependent Plasticity. Notably, it broadens the switching ratio by two orders, minimizes cycle-to-cycle variability, reduces non-linear factors, and achieves an energy consumption of ~30 fJ per synaptic event. Implementation of these enhancements is demonstrated through Artificial Neural Network simulations, yielding a learning accuracy of ~97% on the MNIST hand-written digits dataset. Our findings underscore the significance of defect engineering as a powerful tool in advancing the synaptic functionality of memristors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Accelerating $\mathrm{LHC}$ event generation with simplified pilot runs and fast $\mathrm{PDF}$s

Poor computing efficiency of precision event generators for LHC physics has become a bottleneck for Monte-Carlo event simulation campaigns. We provide solutions to this problem by focusing on two major components of general-purpose event generators: The PDF evaluator and the matrix-element generator. For a typical production setup in the ATLAS experiment, we show that the two can consume about 80% of the total runtime. Using NLO simulations of pp→ℓ + ℓ - +jets and pp→$t\bar{t}$+jets as an example, we demonstrate that the computing footprint of LHAPDF and SHERPA can be reduced by factors of order 10, while maintaining the formal accuracy of the event sample. The improved codes are made publicly available.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1972

A fast numerical method for determining the optimum SNR of an array subject to a Q factor constraint.

The problem of finding the optimum current excitations for a point-detector array which extracts a signal from background noise using the criterion of maximizing the SNR subject to a constraint on the supergain ratio (or Q factor) has recently been solved. The formulation reduced the problem to one of solving for the roots of a polynomial whose coefficients are, in general, complex. We demonstrate that by making a certain transformation, it becomes possible to employ a state variable technique which changes the problem into one of finding the eigenvalues of a real matrix which is considerably faster to do.

Winkler, L. P.↗
NASA NTRS1993

A transfer matrix approach to vibration localization in mistuned blade assemblies

A study of mode localization in mistuned bladed disks is performed using transfer matrices. The transfer matrix approach yields the free response of a general, mono-coupled, perfectly cyclic assembly in closed form. A mistuned structure is represented by random transfer matrices, and the expansion of these matrices in terms of the small mistuning parameter leads to the definition of a measure of sensitivity to mistuning. An approximation of the localization factor, the spatially averaged rate of exponential attenuation per blade-disk sector, is obtained through perturbation techniques in the limits of high and low sensitivity. The methodology is applied to a common model of a bladed disk and the results verified by Monte Carlo simulations. The easily calculated sensitivity measure may prove to be a valuable design tool due to its system-independent quantification of mistuning effects such as mode localization.

Ottarson, Gisli↗
NASA NTRS1999

Angiotensin-converting enzyme and matrix metalloproteinase inhibition with developing heart failure: comparative effects on left ventricular function and geometry

The progression of congestive heart failure (CHF) is left ventricular (LV) myocardial remodeling. The matrix metalloproteinases (MMPs) contribute to tissue remodeling and therefore MMP inhibition may serve as a useful therapeutic target in CHF. Angiotensin converting enzyme (ACE) inhibition favorably affects LV myocardial remodeling in CHF. This study examined the effects of specific MMP inhibition, ACE inhibition, and combined treatment on LV systolic and diastolic function in a model of CHF. Pigs were randomly assigned to five groups: 1) rapid atrial pacing (240 beats/min) for 3 weeks (n = 8); 2) ACE inhibition (fosinopril, 2.5 mg/kg b.i.d. orally) and rapid pacing (n = 8); 3) MMP inhibition (PD166793 2 mg/kg/day p.o.) and rapid pacing (n = 8); 4) combined ACE and MMP inhibition (2.5 mg/kg b.i.d. and 2 mg/kg/day, respectively) and rapid pacing (n = 8); and 5) controls (n = 9). LV peak wall stress increased by 2-fold with rapid pacing and was reduced in all treatment groups. LV fractional shortening fell by nearly 2-fold with rapid pacing and increased in all treatment groups. The circumferential fiber shortening-systolic stress relation was reduced with rapid pacing and increased in the ACE inhibition and combination groups. LV myocardial stiffness constant was unchanged in the rapid pacing group, increased nearly 2-fold in the MMP inhibition group, and was normalized in the ACE inhibition and combination treatment groups. Increased MMP activation contributes to the LV dilation and increased wall stress with pacing CHF and a contributory downstream mechanism of ACE inhibition is an effect on MMP activity.

Non-NASA Center↗
NASA NTRS1996

Immunohistochemical evidence of rapid extracellular matrix remodeling after iron-particle irradiation of mouse mammary gland

High-LET radiation has unique physical and biological properties compared to sparsely ionizing radiation. Recent studies demonstrate that sparsely ionizing radiation rapidly alters the pattern of extracellular matrix expression in several tissues, but little is known about the effect of heavy-ion radiation. This study investigates densely ionizing radiation-induced changes in extracellular matrix localization in the mammary glands of adult female BALB/c mice after whole-body irradiation with 0.8 Gy 600 MeV iron particles. The basement membrane and interstitial extracellular matrix proteins of the mammary gland stroma were mapped with respect to time postirradiation using immunofluorescence. Collagen III was induced in the adipose stroma within 1 day, continued to increase through day 9 and was resolved by day 14. Immunoreactive tenascin was induced in the epithelium by day 1, was evident at the epithelial-stromal interface by day 5-9 and persisted as a condensed layer beneath the basement membrane through day 14. These findings parallel similar changes induced by gamma irradiation but demonstrate different onset and chronicity. In contrast, the integrity of epithelial basement membrane, which was unaffected by sparsely ionizing radiation, was disrupted by iron-particle irradiation. Laminin immunoreactivity was mildly irregular at 1 h postirradiation and showed discontinuities and thickening from days 1 to 9. Continuity was restored by day 14. Thus high-LET radiation, like sparsely ionizing radiation, induces rapid-remodeling of the stromal extracellular matrix but also appears to alter the integrity of the epithelial basement membrane, which is an important regulator of epithelial cell proliferation and differentiation.

NASA Discipline Radiation Health↗