Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix approximation”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

CVD silicon carbide monofilament reinforced SrO-Al2O3-2SiO2 (SAS) glass-ceramic composites

Unidirectional CVD SiC fiber-reinforced SrO.Al2O3.2SiO2 (SAS) glass-ceramic matrix composites have been fabricated by hot pressing at various combinations of temperature, pressure and time. Both carbon-rich surface coated SCS-6 and uncoated SCS-0 fibers were used as reinforcements. Almost fully dense composites have been obtained. Monoclinic celsian, SrAl2Si2O8, was the only crystalline phase observed in the matrix from x-ray diffraction. During three point flexure testing of composites, a test span to thickness ratio of approximately 25 or greater was necessary to avoid sample delamination. Strong and tough SCS-6/SAS composites having a first matrix crack stress of approximately 300 MPa and an ultimate bend strength of approximately 825 MPa were fabricated. No chemical reaction between the SCS-6 fibers and the SAS matrix was observed after high temperature processing. The uncoated SCS-0 fiber-reinforced SAS composites showed only limited improvement in strength over SAS monolithic. The SCS-0/SAS composite having a fiber volume fraction of 0.24 and hot pressed at 1400 deg C exhibited a first matrix cracking stress of approximately 231 +/- 20 MPa and ultimate strength of 265 +/- 17 MPa. From fiber push-out tests, the fiber/matrix interfacial debonding strength (tau(sub debond)) and frictional sliding stress (tau(sub friction)) in the SCS-6/SAS system were evaluated to be approximately 6.7 +/- 2.3 MPa and 4.3 +/- 0.6 MPa, respectively, indicating a weak interface. However, for the SCS-0/SAS composite, much higher values of approximately 17.5 +/- 2.7 MPa for tau(sub debond) and 11.3 +/- 1.6 MPa for tau(sub friction) respectively, were observed; some of the fibers were so strongly bonded to the matrix that they could not be pushed out. Examination of fracture surfaces revealed limited short pull-out length of SCS-0 fibers. The applicability of various micromechanical models for predicting the values of first matrix cracking stress and ultimate strength of these composites were examined.

Bansal, Narottam P.↗
NASA NTRS1977

Failure detection and correction for turbofan engines

In this paper, a failure detection and correction strategy for turbofan engines is discussed. This strategy allows continuing control of the engines in the event of a sensor failure. An extended Kalman filter is used to provide the best estimate of the state of the engine based on currently available sensor outputs. Should a sensor failure occur the control is based on the best estimate rather than the sensor output. The extended Kalman filter consists of essentially two parts, a nonlinear model of the engine and up-date logic which causes the model to track the actual engine. Details on the model and up-date logic are presented. To allow implementation, approximations are made to the feedback gain matrix which result in a single feedback matrix which is suitable for use over the entire flight envelope. The effect of these approximations on stability and response is discussed. Results from a detailed nonlinear simulation indicate that good control can be maintained even under multiple failures.

Corley, R. C.↗
NASA NTRS1996

Effects of Fiber Content on Mechanical Properties of CVD SiC Fiber-Reinforced Strontium Aluminosilicate Glass-Ceramic Composites

Unidirectional CVD SiC(f)(SCS-6) fiber-reinforced strontium aluminosilicate (SAS) glass-ceramic matrix composites containing various volume fractions, approximately 16 to 40 volume %, of fibers were fabricated by hot pressing at 1400 C for 2 h under 27.6 MPa. Monoclinic celsian, SrAl2Si2O8, was the only crystalline phase formed, with complete absence of the undesired hexacelsian phase, in the matrix. Room temperature mechanical properties were measured in 3-point flexure. The matrix microcracking stress and the ultimate strength increased with increase in fiber volume fraction, reached maximum values for V(sub f) approximately equal to 0.35, and degraded at higher fiber loadings. This degradation in mechanical properties is related to the change in failure mode, from tensile at lower V(sub f) to interlaminar shear at higher fiber contents. The extent of fiber loading did not have noticeable effect on either fiber-matrix debonding stress, or frictional sliding stress at the interface. The applicability of micromechanical models in predicting the mechanical properties of the composites was also examined. The currently available theoretical models do not appear to be useful in predicting the values of the first matrix cracking stress, and the ultimate strength of the SCS-6/SAS composites.

Bansal, Narottam P.↗
NASA NTRS1977

Structured estimation - Sample size reduction for adaptive pattern classification

The Gaussian two-category classification problem with known category mean value vectors and identical but unknown category covariance matrices is considered. The weight vector depends on the unknown common covariance matrix, so the procedure is to estimate the covariance matrix in order to obtain an estimate of the optimum weight vector. The measure of performance for the adapted classifier is the output signal-to-interference noise ratio (SIR). A simple approximation for the expected SIR is gained by using the general sample covariance matrix estimator; this performance is both signal and true covariance matrix independent. An approximation is also found for the expected SIR obtained by using a Toeplitz form covariance matrix estimator; this performance is found to be dependent on both the signal and the true covariance matrix.

Morgera, S.↗
NASA NTRS2016

Hydrogen Isotopic Composition of Apatite in Northwest Africa 7034: A Record of the "Intermediate" H-Isotopic Reservoir in the Martian Crust?

Northwest Africa (NWA) 7034 and its pairings comprise a regolith breccia with a basaltic bulk composition [1] that yields a better match than any other martian meteorite to visible-infrared reflectance spectra of the martian surface measured from orbit [2]. The composition of the fine-grained matrix within NWA 7034 bears a striking resemblance to the major element composition estimated for the martian crust, with several exceptions. The NWA 7034 matrix is depleted in Fe, Ti, and Cr and enriched in Al, Na, and P [3]. The differences in Al and Fe are the most substantial, but the Fe content of NWA 7034 matrix falls within the range reported for the southern highlands crust [6]. It was previously suggested by [4] that NWA 7034 was sourced from the southern highlands based on the ancient 4.4 Ga ages recorded in NWA 7034/7533 zircons [4, 5]. In addition, the NWA 7034 matrix material is enriched in incompatible trace elements by a factor of 1.2-1.5 [7] relative to estimates of the bulk martian crust. The La/Yb ratio of the bulk martian crust is estimated to be approximately 3 [7], and the La/Yb of the NWA 7034 matrix materials ranges from approximately 3.9 to 4.4 [3, 8], indicating a higher degree of LREE enrichment in the NWA 7034 matrix materials. This elevated La/Yb ratio and enrichment in incompatible lithophile trace elements is consistent with NWA 7034 representing a more geochemically enriched crustal terrain than is represented by the bulk martian crust, which would be expected if NWA 7034 represents the bulk crust from the southern highlands. Given the similarities between NWA 7034 and the martian crust, NWA 7034 may represent an important sample for constraining the composition of the martian crust, particularly the ancient highlands. In the present study, we seek to constrain the H isotopic composition of the martian crust using Cl-rich apatite in NWA 7034. Usui et al., [9] recently proposed that a H isotopic reservoir exists within the martian crust that has a H-isotopic composition that is intermediate (δD of 1000-2000per mille) between an isotopically light mantle (Delta D is less than 275per mille [10]) and an isotopically heavy atmosphere (δD of 2500-6100per mille [11, 12]). Apatites in NWA 7034 occur in a number of lithologic domains, however apatites across all lithologic domains have been affected by a Pb-loss event at about 1.5 Ga before present [5], so they are unlikely to have retained a primary composition and are more likely to have equilibrated with fluids within the martian crust that may or may not have exchanged with the martian atmosphere. Equilibration of apatite with crustal fluids is further supported by the chlorine-rich compositions exhibited by apatites in NWA 7034 in comparison to apatites from other martian meteorites (Figure 1; [13]). Cl is more hydrophilic than F, which promotes formation of Cl-rich apatite compositions in fluid-rich systems [e.g., 14, 15-17].

McCubbin, F. M.↗
DOE OSTI2023

Systems and methods for customizing kernel machines with deep neural networks

A method including receiving an input data set. The input data set can include one of a feature domain set or a kernel matrix. The method also can include constructing dense embeddings using: (i) Nyström approximations on the input data set when the input data set comprises the kernel matrix, and (ii) clustered Nyström approximations on the input data set when the input data set comprises the feature domain set. The method additionally can include performing representation learning on each of the dense embeddings using a multi-layer fully-connected network for each of the dense embeddings to generate latent representations corresponding to each of the dense embeddings. The method further can include applying a fusion layer to the latent representations corresponding to the dense embeddings to generate a combined representation. The method additionally can include performing classification on the combined representation. Other embodiments of related systems and methods are also disclosed.

Song, Huan↗
NASA NTRS1988

Statistical analysis of effective singular values in matrix rank determination

A major problem in using SVD (singular-value decomposition) as a tool in determining the effective rank of a perturbed matrix is that of distinguishing between significantly small and significantly large singular values to the end, conference regions are derived for the perturbed singular values of matrices with noisy observation data. The analysis is based on the theories of perturbations of singular values and statistical significance test. Threshold bounds for perturbation due to finite-precision and i.i.d. random models are evaluated. In random models, the threshold bounds depend on the dimension of the matrix, the noisy variance, and predefined statistical level of significance. Results applied to the problem of determining the effective order of a linear autoregressive system from the approximate rank of a sample autocorrelation matrix are considered. Various numerical examples illustrating the usefulness of these bounds and comparisons to other previously known approaches are given.

Konstantinides, Konstantinos↗
NASA NTRS1974

A computational procedure to analyze metal matrix laminates with nonlinear lamination residual strains

An approximate computational procedure is described for the analysis of angleplied laminates with residual nonlinear strains. The procedure consists of a combination of linear composite mechanics and incremental linear laminate theory. The procedure accounts for initial nonlinear strains, unloading, and in-situ matrix orthotropic nonlinear behavior. The results obtained in applying the procedure to boron/aluminum angleplied laminates show that this is a convenient means to accurately predict the initial tangent properties of angleplied laminates in which the matrix has been strained nonlinearly by the lamination residual stresses. The procedure predicted initial tangent properties results which were in good agreement with measured data obtained from boron/aluminum angleplied laminates.

Chamis, C. C.↗
NASA NTRS1974

A computational procedure to analyze metal matrix laminates with nonlinear lamination residual strains

An approximate computational procedure is described for the analysis of angleplied laminates with residual nonlinear strains. The procedure consists of a combination of linear composite mechanics and incremental linear laminate theory. The procedure accounts for initial nonlinear strains, unloading, and in-situ matrix orthrotropic nonlinear behavior. The results obtained in applying the procedure to boron/aluminum angleplied laminates show that this is a convenient means to accurately predict the initial tangent properties of angleplied laminates in which the matrix has been strained nonlinearly by the lamination residual stresses. The procedure predicted initial tangent properties results which were in good agreement with measured data obtained from boron/aluminum angleplied laminates.

Chamis, C. C.↗
NASA NTRS2007

Spatial Pattern of Cell Damage in Tissue from Heavy Ions

A new Monte Carlo algorithm was developed that can model passage of heavy ions in a tissue, and their action on the cellular matrix for 2- or 3-dimensional cases. The build-up of secondaries such as projectile fragments, target fragments, other light fragments, and delta-rays was simulated. Cells were modeled as a cell culture monolayer in one example, where the data were taken directly from microscopy (2-d cell matrix). A simple model of tissue was given as abstract spheres with close approximation to real cell geometries (3-d cell matrix), as well as a realistic model of tissue was proposed based on microscopy images. Image segmentation was used to identify cells in an irradiated cell culture monolayer, or slices of tissue. The cells were then inserted into the model box pixel by pixel. In the case of cell monolayers (2-d), the image size may exceed the modeled box size. Such image was is moved with respect to the box in order to sample as many cells as possible. In the case of the simple tissue (3-d), the tissue box is modeled with periodic boundary conditions, which extrapolate the technique to macroscopic volumes of tissue. For real tissue, specific spatial patterns for cell apoptosis and necrosis are expected. The cell patterns were modeled based on action cross sections for apoptosis and necrosis estimated based on BNL data, and other experimental data.

Ponomarev, Artem L.↗
NASA NTRS2002

Charge Resolution of the Silicon Matrix of the ATIC Experiment

ATIC (Advanced Thin Ionization Calorimeter) is a balloon borne experiment designed to measure the cosmic ray composition for elements from hydrogen to iron and their energy spectra from approx.50 GeV to near 100 TeV. It consists of a Si-matrix detector to determine the charge of a CRT particle, a scintillator hodoscope for tracking, carbon interaction targets and a fully active BGO calorimeter. ATIC had its first flight from McMurdo, Antarctica from 28/12/2000 to 13/01/2001. The ATIC flight collected approximately 25 million events. The silicon matrix of the ATIC spectrometer is designed to resolve individual elements from proton to iron. To provide this resolution careful calibration of each pixel of the silicon matrix is required. Firstly, for each electronic channel of the matrix the pedestal value was subtracted taking into account its drift during the flight. The muon calibration made before the flight was used then to convert electric signals (in ADC channel number) to energy deposits in each pixel. However, the preflight muon calibration was not accurate enough for the purpose, because of lack of statistics in each pixel. To improve charge resolution the correction was done for the position of Helium peak in each pixel during the flight . The other way to set electric signals in electronics channels of the Si-matrix to one scale was correction for electric channel gains accurately measured in laboratory. In these measurements it was found that small different nonlinearities for different channels are present in the region of charge Z > 20. The correction for these non-linearities was not done yet. In linear approximation the method provides practically the same resolution as muon calibration plus He-peak correction. For searching a pixel with the signal of primary particle an indication from the cascade in the calorimeter was used. For this purpose a trajectory was reconstructed using weight centers of energy deposits in BGO layers. The point of intersection of this trajectory with Si-matrix and its RMS was determined. The pixel with maximal signal in 3sigma region was taken as sought. The signal in this pixel was corrected by trajectory zenith angle. The preliminary results on charge resolution of the Si-matrix in the range from protons to iron are presented.

Zatsepin, V. I.↗
DOE OSTI2024

Spectral scheme for atomic structure calculations in density functional theory

In this study, we present a spectral scheme for atomic structure calculations in pseudopotential Kohn-Sham density functional theory. In particular, after applying an exponential transformation of the radial coordinates, we employ global polynomial interpolation on a Chebyshev grid, with derivative operators approximated using the Chebyshev differentiation matrix, and integrations using Clenshaw-Curtis quadrature. We demonstrate the accuracy and efficiency of the scheme through spin-polarized and unpolarized calculations for representative atoms, while considering local, semilocal, and hybrid exchange-correlation functionals. In particular, we find that $\mathcal{O}$(200) grid points are sufficient to achieve an accuracy of 1 microhartree in the eigenvalues for optimized norm conserving Vanderbilt pseudopotentials spanning the periodic table from atomic number Ζ = 1 to 83.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Model-based and Model-free Designs for an Extended Continuous-time LQR with Exogenous Inputs

We present an extended linear quadratic regulator (LQR) design for continuous-time linear time-invariant (LTI) systems in the presence of exogenous inputs. We first propose a model-based solution with cost minimization guarantees for states and inputs using dynamic programming (DP). The control law consists of a combination of the optimal state feedback and an additional optimal term dependent on the exogenous inputs. The control gains for the two components are obtained by solving a set of matrix differential equations. We provide these solutions for both finite horizons and steady-state cases. In the second part of the paper, we formulate a reinforcement learning (RL) based algorithm which does not need any model information except the input matrix, and can compute an approximate steady-state LQR gain using measurements of the states, the control inputs, and the exogenous inputs. Both model-based and data-driven optimal control algorithms are tested with a numerical example under different exogenous inputs showcasing the effectiveness of the designs.

Mukherjee, Sayak↗
DOE OSTI2020

Covariance generation for the prompt neutron multiplicity of 239 Pu including the (n, γ f) process

Fission cross section of 239 Pu can be seen as a sum of the “immediate" fission and “two-step" (n,γf) reactions. In the Resolved Resonance Range of the reaction cross sections, the contribution of the (n,γf) process has an impact on the determination of the partial widths magnitude involved in the Reich-Moore approximation of the R-matrix theory. The present work aims to investigate this impact by using the CONRAD code and the partial width Γ γf for the (n,γf) reaction calculated by Lynn et al. A special attention will be paid to the covariance matrix obtained on ν p .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Expanded capabilities of the CarMa code in modeling resistive wall mode dynamics with 3-D conductors

In this work, an improved version of the CarMa code is presented, called CarMa-D, for the analysis of resistive wall modes (RWMs) in fusion devices, simultaneously considering the effects of volumetric three-dimensional conducting structures, plasma dynamics, toroidal rotation or drift-kinetic damping. Unlike static CarMa, CarMa-D does not rely on the simplifying assumptions such as neglecting the plasma mass, toroidal rotation and kinetic damping. The new coupling strategy is based on matrix-based Padé rational functions approximation of plasma a response. The arising mathematical model is formally equivalent to the original CarMa model, but with a higher number of degrees of freedom to model the dynamics of the plasma. CarMa-D overcomes the main limitations of the original CarMa, in particular: (i) the massless assumption for the plasma is removed, allowing modeling of global modes growing both on ideal kink time scales and in the typical RWM growth rate regime, with a suitable treatment of the model; (ii) the effects of toroidal plasma flow and drift kinetic damping can be included into the new model, providing a powerful tool to study macroscopic phenomena where both plasma dynamics and 3-D conducting structures play important roles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Fully self-consistent finite-temperature $GW$ in Gaussian Bloch orbitals for solids

In this work, we present algorithmic and implementation details for the fully self-consistent finite-temperature $GW$ method in Gaussian Bloch orbitals for solids. Our implementation is based on the finite-temperature Green's function formalism in which all equations are solved on the imaginary axis, without resorting to analytical continuation during the self-consistency. No quasiparticle approximation is employed and all matrix elements of the self-energy are explicitly evaluated. The method is tested by evaluating the band gaps of selected semiconductors and insulators. We show agreement with other, differently formulated, finite-temperature sc ⁢$GW$ implementations when finite-size corrections and basis-set errors are taken into account. By migrating computationally intensive calculations to graphics processing units, we obtain scalable results on large supercomputers with nearly optimal performance. Our work demonstrates the applicability of Gaussian orbital based sc⁢ $GW$ for ab initio correlated material simulations and provides a sound starting point for embedding methods built on top of $GW$.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

The XFaster Power Spectrum and Likelihood Estimator for the Analysis of Cosmic Microwave Background Maps

We present the XFaster analysis package, XFaster is a fast, iterative angular power spectrum estimator based on a diagonal approximation to the quadratic Fisher matrix estimator. XFaster uses Monte Carlo simulations to compute noise biases and filter transfer functions and is thus a hybrid of both Monte Carlo and quadratic estimator methods. In contrast to conventional pseudo-C ℓ based methods, the algorithm described here requires a minimal number of simulations, and does not require them to be precisely representative of the data to estimate accurate covariance matrices for the bandpowers. The formalism works with polarization-sensitive observations and also data sets with identical, partially overlapping, or independent survey regions. The method was first implemented for the analysis of BOOMERanG data (Netterfield et al. 2002; Jones et al. 2006), and also used as part of the Planck analysis (Rocha et al. 2011). Here, we describe the full, publicly available analysis package, written in Python, as developed for the analysis of data from the 2015 flight of the SPIDER instrument (SPIDER Collaboration 2021). The package includes extensions for self-consistently estimating null spectra and for estimating fits for Galactic foreground contributions. We show results from the extensive validation of XFaster using simulations, and its application to the SPIDER data set.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1980

Spectral methods for problems in complex geometries

The properties of spectral methods are surveyed and their extension to solve problems in complex geometries is developed. A new iteration procedure is introduced to solve efficiently the full matrix equations resulting from spectral approximations to nonconstant coefficient boundary-value problems in complex geometries. It is shown that the work required to solve these spectral equations exceeds that of solving the lowest-order finite-difference approximation to the same problem by only O(N log N).

Orszag, S. A.↗