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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Flap-Lag-Torsion Stability in Forward Flight

An aeroelastic stability of three-degree flap-lag-torsion blade in forward flight is examined. Quasisteady aerodynamics with a dynamic inflow model is used. The nonlinear time dependent periodic blade response is calculated using an iterative procedure based on Floquet theory. The periodic perturbation equations are solved for stability using Floquet transition matrix theory as well as constant coefficient approximation in the fixed reference frame. Results are presented for both stiff-inplane and soft-inplane blade configurations. The effects of several parameters on blade stability are examined, including structural coupling, pitch-flap and pitch-lag coupling, torsion stiffness, steady inflow distribution, dynamic inflow, blade response solution and constant coefficient approximation.

Panda, B.↗
NASA NTRS1991

Critical speeds and forced response solutions for active magnetic bearing turbomachinery, part 1

The prediction of critical speeds and forced response of active magnetic bearing turbomachinery is of great interest due to the increased use of this new and promising technology. Calculating the system undamped critical speeds and forced response is important to all those who are involved in the design of the active magnetic bearing system. An extended Jeffcott model which was used as an approximate solution to a more accurate transfer matrix procedure is presented. Theory behind a two-degree-of freedom extended Jeffcoat model is presented. Results of the natural frequency calculation are shown followed by the results of the forced response calculation. The system response was predicted for two types of forcing. A constant magnitude excitation with a wide frequency variation was applied at the bearings as one forcing function. The normal unbalance force at the midspan was the second source of excitation. The results of this extended Jeffcott solution gives useful design guidance for the influence of the first and third modes of a symmetric rotor system.

Keesee, J.↗
NASA NTRS1989

Optimal Chebyshev polynomials on ellipses in the complex plane

The design of iterative schemes for sparse matrix computations often leads to constrained polynomial approximation problems on sets in the complex plane. For the case of ellipses, we introduce a new class of complex polynomials which are in general very good approximations to the best polynomials and even optimal in most cases.

Fischer, Bernd↗
NASA NTRS1990

Bearingless rotor aeromechanical stability measurements and correlations using nonlinear aerodynamics

The aeromechanical stability of a 1/8th Froude scaled bearingless rotor model was investigated experimentally in a wind tunnel. Both shaft-fixed and shaft-free conditions were examined to study the aeroelastic stability of a bearingless rotor without the incorporation of auxiliary dampers. This wind tunnel investigation generated a set of stability data for four different advance ratios, and a wide range of collective pitch settings. Theoretical analysis was performed using the newly developed University of Maryland Advanced Rotorcraft Code (UMARC). For analysis, the blade is modeled as an elastic beam undergoing flap bending, lag bending, elastic twist, and axial deformation. Blade response is calculated using a finite element method in time. Nonlinear aerodynamic effects are included by using a semiempirical stall modeling. The linearized periodic rotor perturbation equations in the nonrotating frame are solved for stability roots using Floquet transition matrix theory, as well as constant coefficient approximation. The predicted results are compared with the experimental data.

Wang, James M.↗
NASA NTRS1992

Electron impact polarization of atomic spectral lines. I - A general theoretical scheme

A suitable theoretical scheme able to describe, in a wide variety of astrophysical situations, the phenomenon of atomic line polarization by electron impact is developed. Starting from the general principles of quantum mechanics and assuming the Born approximation, the rate equations for the density matrix elements of a multilevel atomic system, interacting with a nonrelativistic electron beam having any kind of angular distribution, are derived in full generality. The resulting theory generalizes the previous ones by accounting for the collisional rates and the cross sections concerning both inelastic and superelastic collisions (in any geometrical situation), and, moreover, by taking into account the coherences among Zeeman sublevels split by a magnetic field. As an example of particular relevance, the general formulas derived in the first sections of the paper are subsequently particularized to the case of the electric dipole interaction.

Fineschi, Silvano↗
NASA NTRS2003

Many-Body Effects in Quantum-Well Intersubband Transitions

Intersubband polarization couples to collective excitations of the interacting electron gas confined in a semiconductor quantum well (Qw) structure. Such excitations include correlated pair excitations (repellons) and intersubband plasmons (ISPs). The oscillator strength of intersubband transitions (ISBTs) strongly varies with QW parameters and electron density because of this coupling. We have developed a set of kinetic equations, termed the intersubband semiconductor Bloch equations (ISBEs), from density matrix theory with the Hartree-Fock approximation, that enables a consistent description of these many-body effects. Using the ISBEs for a two-conduction-subband model, various many-body effects in intersubband transitions are studied in this work. We find interesting spectral changes of intersubband absorption coefficient due to interplay of the Fermi-edge singularity, subband renormalization, intersubband plasmon oscillation, and nonparabolicity of bandstructure. Our results uncover a new perspective for ISBTs and indicate the necessity of proper many-body theoretical treatment in order for modeling and prediction of ISBT line shape.

Li, Jian-Zhong↗
NASA NTRS2011

On Time Delay Margin Estimation for Adaptive Control and Optimal Control Modification

This paper presents methods for estimating time delay margin for adaptive control of input delay systems with almost linear structured uncertainty. The bounded linear stability analysis method seeks to represent an adaptive law by a locally bounded linear approximation within a small time window. The time delay margin of this input delay system represents a local stability measure and is computed analytically by three methods: Pade approximation, Lyapunov-Krasovskii method, and the matrix measure method. These methods are applied to the standard model-reference adaptive control, s-modification adaptive law, and optimal control modification adaptive law. The windowing analysis results in non-unique estimates of the time delay margin since it is dependent on the length of a time window and parameters which vary from one time window to the next. The optimal control modification adaptive law overcomes this limitation in that, as the adaptive gain tends to infinity and if the matched uncertainty is linear, then the closed-loop input delay system tends to a LTI system. A lower bound of the time delay margin of this system can then be estimated uniquely without the need for the windowing analysis. Simulation results demonstrates the feasibility of the bounded linear stability method for time delay margin estimation.

Nguyen, Nhan T.↗
NASA NTRS2014

Generating Inviscid and Viscous Fluid Flow Simulations over a Surface Using a Quasi-simultaneous Technique

A fluid-flow simulation over a computer-generated surface is generated using a quasi-simultaneous technique. The simulation includes a fluid-flow mesh of inviscid and boundary-layer fluid cells. An initial fluid property for an inviscid fluid cell is determined using an inviscid fluid simulation that does not simulate fluid viscous effects. An initial boundary-layer fluid property a boundary-layer fluid cell is determined using the initial fluid property and a viscous fluid simulation that simulates fluid viscous effects. An updated boundary-layer fluid property is determined for the boundary-layer fluid cell using the initial fluid property, initial boundary-layer fluid property, and an interaction law. The interaction law approximates the inviscid fluid simulation using a matrix of aerodynamic influence coefficients computed using a two-dimensional surface panel technique and a fluid-property vector. An updated fluid property is determined for the inviscid fluid cell using the updated boundary-layer fluid property.

Sturdza, Peter↗
NASA NTRS2023

Stand Alone Electrical Equipment PRL 227a (Formerly PRL 227 LOC3) Confirmatory Sampling Report Addendum

HydroGeoLogic, Inc. (HGL) prepared this Confirmatory Sampling Report Addendum (CSRA) for the National Aeronautics and Space Administration (NASA), Kennedy Space Center (KSC), Florida. PRL 227a Location of Concern (LOC) 3 is located on 10th Street, west of the former Vertical Processing Facility (M7-1469), SWMU 077. The Load Break Switches (LBS) were installed in 1986 and are suspected to have formerly contained dielectric fluid that contained polychlorinated biphenyls (PCBs). The LBS were replaced in 1993 with new LBSs that do not contain PCBs. The LBS are located on a concrete pad surrounded by grass that is mowed and maintained. A total of 330 soil samples were collected between June 2020 and December 2022 from 88 locations at depths ranging from 0-0.5 to 19-20 ft bls. In general, soil samples collected from depths less than 5 ft were collected via hand auger and samples collected beyond 5 ft bls were collected using a direct-push technology drill rig. Soil samples collected from greater than 5 ft bls to 20 ft bls were collected to confirm PCBs were not present at depth in the aquifer matrix. Delineation samples were collected using an approximately 20 ft horizontal spacing. PCB impacted soil has been completely delineated to the State of Florida residential soil cleanup target level (SCTL) at PRL 227a horizontally and vertically. PCB concentrations have been identified above the residential and industrial SCTL as well as concentrations greater than or equal to 50 ppm requiring coordination with EPA and disposal at a facility authorized under TSCA. Based upon the CSRA sampling, an interim measure is recommended to remove PCB impacted soil above the residential SCTL. The total volume of PCB impacted soil to be removed is 591.6 cubic yards. This includes approximately 9.3 cubic yards of soil with PCB concentrations greater than or equal to 50 ppm. An interim measure work plan is included as an attachment to this report.

electrical equipment↗
DOE OSTI2024

Symmetry-Based Structured Matrices for Efficient Approximately Equivariant Networks

There has been much recent interest in designing symmetry-aware neural networks (NNs) exhibiting relaxed equivariance. Such NNs aim to interpolate between being exactly equivariant and being fully flexible, affording consistent performance benefits. In a separate line of work, certain structured parameter matrices -- those with displacement structure, characterized by low displacement rank (LDR) -- have been used to design small-footprint NNs. Displacement structure enables fast function and gradient evaluation, but permits accurate approximations via compression primarily to classical convolutional neural networks (CNNs). In this work, we propose a general framework -- based on a novel construction of symmetry-based structured matrices -- to build approximately equivariant NNs with significantly reduced parameter counts. Our framework integrates the two aforementioned lines of work via the use of so-called Group Matrices (GMs), a forgotten precursor to the modern notion of regular representations of finite groups. GMs allow the design of structured matrices -- resembling LDR matrices -- which generalize the linear operations of a classical CNN from cyclic groups to general finite groups and their homogeneous spaces. We show that GMs can be employed to extend all the elementary operations of CNNs to general discrete groups. Further, the theory of structured matrices based on GMs provides a generalization of LDR theory focussed on matrices with cyclic structure, providing a tool for implementing approximate equivariance for discrete groups. We test GM-based architectures on a variety of tasks in the presence of relaxed symmetry. We report that our framework consistently performs competitively compared to approximately equivariant NNs, and other structured matrix-based compression frameworks, sometimes with a one or two orders of magnitude lower parameter count.

Samudre, Ashwin↗
DOE OSTI2026

First Detection of the Baryon Acoustic Oscillation (BAO) Feature in the 3-Point Correlation Function of DESI DR1 Luminous Red Galaxies

We present the first detection of the 3-Point Correlation Function (3PCF) Baryon Acoustic Oscillation (BAO) signal from the DESI Data Release 1 (DR1) sample of Luminous Red Galaxies (LRGs), which contains over 2.1 million galaxies. Our analysis is based on a tree-level redshift-space bispectrum template, which is then transformed to position space using the Fast Fourier Transform on Logarithmic scales (FFTLog) algorithm. We detect the BAO feature with a significance of approximately $8.1σ$ using the EZmock covariance matrix and $8.5σ$ using the analytical covariance matrix, for the full LRG redshift range ($0.4

Kamalinejad, Farshad [Florida U.] (ORCID:000000017↗
DOE OSTI2023

FastACE

SAND2024-01893O The Fast Adaptive Cosine Estimator (FastACE) algorithm modifies the well-known Adaptive Cosine Estimator for target detection in hyperspectral imagery. Specifically, FastACE modifies the computation of the background precision matrix (C_b^-1) under a Vecchia approximation. That is, each spectral band, conditioned on a local neighborhood of bands around that band, is independent of the other spectral bands. The FastACE algorithm leverages a parameterizable, auto-regressive neighborhood for the conditional independence assumption. The underlying math used for computing detection scores is equivalent to ACE, but it is implemented within FastACE. The software implements a target detection algorithm and associated utilities for target detection in hyperspectral imagery. Provided with a hyperspectral image and corresponding target signature, it produces relevant background statistics and the corresponding detection scores for the target signature in the image. The software is designed to integrate into other end-user applications or processing. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

VanderLaan, John↗
DOE OSTI2020

A computationally-efficient method for flamelet calculations

A new open-source code for the simulation of the diffusion flamelet equations is proposed. Emphasis is placed on using an approximate Jacobian to reduce the computational cost of the matrix operations. Performance of the proposed solvers is tested by performing flamelet calculations with kinetic mechanisms of varying sizes. For the unity Lewis number equations, the present iterative Newton solver using an approximate Jacobian greatly outperforms direct Newton solvers using exact Jacobians. The computation cost scales linearly with the number of species, leading to a reduction in solution times by two orders of magnitude for mechanisms containing thousands of species. The applicability of the Jacobian approximations to the solution of the non-unity Lewis number flamelet equations is assessed. The approximations are generally inadequate to solve the full non-unity Lewis number equations but can be used in some applications depending on the balance of terms in the flamelet equations. As an example, the flamelet solver is applied to the study of sooting tendencies in laminar co-flow diffusion flames where modified non-unity Lewis number flamelet equations, previously shown to accurately reproduce experimentally-measured Yield Sooting Indices (YSI), are solved. Here, the accelerated flamelet solver is well suited for sensitivity analysis and uncertainty quantification with large detailed kinetic mechanisms, tasks for which the computational cost was previously prohibitive.

42 ENGINEERING↗
NASA NTRS1996

Bit Error Probability for Maximum Likelihood Decoding of Linear Block Codes

In this paper, the bit error probability P(sub b) for maximum likelihood decoding of binary linear codes is investigated. The contribution of each information bit to P(sub b) is considered. For randomly generated codes, it is shown that the conventional approximation at high SNR P(sub b) is approximately equal to (d(sub H)/N)P(sub s), where P(sub s) represents the block error probability, holds for systematic encoding only. Also systematic encoding provides the minimum P(sub b) when the inverse mapping corresponding to the generator matrix of the code is used to retrieve the information sequence. The bit error performances corresponding to other generator matrix forms are also evaluated. Although derived for codes with a generator matrix randomly generated, these results are shown to provide good approximations for codes used in practice. Finally, for decoding methods which require a generator matrix with a particular structure such as trellis decoding or algebraic-based soft decision decoding, equivalent schemes that reduce the bit error probability are discussed.

Lin, Shu↗
NASA NTRS2003

Fabrication And Evaluation Of Sic/Sic Tubes With Various Fiber Architectures

SiC/SiC composites are excellent material candidates for high temperature applications where the performance requirements are high strength, high creep-rupture resistance, high environmental durability, and high thermal conductivity. In the past, the NASA UEET program has demonstrated fabrication of high-performance SiC/SiC flat panels reinforced by Sylramic-iBN SiC fibers. Currently NASA UEET is scaling up this SiC/SiC system by fabrication of more complex shaped components using the same fiber type. This paper reports the effects of various fiber architectures on the processing, mechanical, and durability behavior of small-diameter 0.5" ID SiC/SiC tubes, which are potential sub-elements for leading edges and cooling channels in turbine vanes and blades. Nine different fiber architectures were utilized for construction of seamless tube preforms, from simple 2D jelly-rolling to complex braiding, pin-weaving, filament-winding and 3D orthogonal weaving with approximately 5% fibers in the thru-thickness direction. Using the BN interphase and Sic matrix processing steps established for the flat panels, SiC/SiC tubes were fabricated with wall thicknesses of approximately 60 mils and total fiber fractions of approximately 35%. The "D" split ring tests for hoop tensile properties, micro-structural examinations for relationship between fiber architecture formation and matrix infiltration, and the low-pressure burner rig tests for the high temperature durability under thru-thickness thermal gradient were conducted. The better matrix infiltration and higher hoop strength were achieved using the tri-axial braided and the three-float pin woven SiC/SiC tubes. In general, it needs not only higher hoop direction fibers but also axial direction fibers for the higher hoop strength and the better infiltration, respectively. These results are analyzed to offer general guidelines for selecting fiber pre-form architectures and SiC/SiC processes that maximize tube hoop strength, thru-thickness thermal conductivity, and burner-rig durability under a high thermal gradient.

Yun, H. M.↗
NASA NTRS2014

Phosphates at the Surface of Mars: Primary Deposits and Alteration Products

Phosphorus is an essential element in terrestrial organisms and thus characterizing the occurrences of phosphate phases at the martian surface is crucial in the assessment of habitability. The Alpha Particle X-Ray Spectrometers onboard Spirit, Opportunity and Curiosity discovered a variety of primary and secondary phosphate phases allowing direct comparisons across the three landing sites. The Spirit rover at Gusev Crater encountered the "Wishstone/Watchtower" class of P-rich (up to 5.2 wt% P2O5) rocks interpreted to be alkaline volcanic rocks with a physical admixture of approximately 10 to 20% merrillite [Usui et al 2008]. These rocks are characterized by elevated Ti and Y and anomalously low Cr and Ni, which could largely reflect the nature of the protoliths: Evolved magmatic rocks. Many of these chemical signatures are also found in pyroclastic deposits at nearby "Home Plate" and in phosphate precipitates derived from fluid interactions with these rocks ("Paso Robles" soils). The Opportunity rover at Meridiani Planum recently analyzed approximately 4 cm clast in a fine-grained matrix, one of numerous rocks of similar appearance at the rim of Endeavour Crater. This clast, "Sarcobatus," has minor enrichments in Ca and P relative to the matrix, and like the P-rich rocks at Gusev, Sarcobatus also shows elevated Al and Ti. On the same segment of the Endeavour rim, subsurface samples were found with exceptional levels of Mn (approximately 3.5 wt% MnO). These secondary and likely aqueous deposits contain strong evidence for associated Mg-sulfate and Ca-phosphate phases. Finally, the Curiosity traverse at Gale crater encountered P-rich rocks compositionally comparable to Wishstone at Gusev, including elevated Y. Phosphorous-rich rocks with similar chemical characteristics are prevalent on Mars, and the trace and minor element signatures provide constraints on whether these are primary deposits, secondary products of physical weathering or secondary products of chemical weathering.

Yen, Albert S.↗
DOE OSTI2024

The Effect of the Prior and the Experimental Design on the Inference of the Precision Matrix in Gaussian Chain Graph Models

Here, we investigate whether (and how) experimental design could aid in the estimation of the precision matrix in a Gaussian chain graph model, especially the interplay between the design, the effect of the experiment and prior knowledge about the effect. Estimation of the precision matrix is a fundamental task to infer biological graphical structures like microbial networks. We compare the marginal posterior precision of the precision matrix under four priors: flat, conjugate Normal-Wishart, Normal-MGIG and a general independent. Under the flat and conjugate priors, the Laplace-approximated posterior precision is not a function of the design matrix rendering useless any efforts to find an optimal experimental design to infer the precision matrix. In contrast, the Normal-MGIG and general independent priors do allow for the search of optimal experimental designs, yet there is a sharp upper bound on the information that can be extracted from a given experiment. We confirm our theoretical findings via a simulation study comparing (i) the KL divergence between prior and posterior and (ii) the Stein’s loss difference of MAPs between random and no experiment. Our findings provide practical advice for domain scientists conducting experiments to better infer the precision matrix as a representation of a biological network.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1993

Visual optimization of DCT quantization matrices for individual images

Many image compression standards (JPEG, MPEG, H.261) are based on the Discrete Cosine Transform (DCT). However, these standards do not specify the actual DCT quantization matrix. We have previously provided mathematical formulae to compute a perceptually lossless quantization matrix. Here I show how to compute a matrix that is optimized for a particular image. The method treats each DCT coefficient as an approximation to the local response of a visual 'channel'. For a given quantization matrix, the DCT quantization errors are adjusted by contrast sensitivity, light adaptation, and contrast masking, and are pooled non-linearly over the blocks of the image. This yields an 8x8 'perceptual error matrix'. A second non-linear pooling over the perceptual error matrix yields total perceptual error. With this model we may estimate the quantization matrix for a particular image that yields minimum bit rate for a given total perceptual error, or minimum perceptual error for a given bit rate. Custom matrices for a number of images show clear improvement over image-independent matrices. Custom matrices are compatible with the JPEG standard, which requires transmission of the quantization matrix.

Watson, Andrew B.↗