Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MATRICES”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Controlling matric potential in microfluidics to examine microbial dynamics in unsaturated porous media

The use of microfluidics for the study of soil microbial ecology is an emerging field. Most microfluidic studies of biological systems, however, have been performed under fully saturated conditions that are not representative of natural soil. Therefore, while microfluidics offer many unique capabilities that other methodologies cannot, they are not currently suited to address the effects of matric potential, an important variable defining the microbial moisture niche. Here, a methodology is presented that allows the user to control the aqueous conditions within microfluidic networks by manipulating matric potential using a hanging water column. The method relies on hydrophilic surface treatment of the microfluidic device using polyvinyl alcohol (PVA) and incorporating a bed of small pores at the network boundaries, which serve as a porous ceramic plate analogue (PPA). The method was validated on a simple capillary bundle and then on a more complex pore network. A water retention curve, exhibiting hysteresis, was generated for the pore network over a narrow matric potential range of 0 to – 5 kPa. Both the drainage and wetting curves were reproducible, as were the spatial configuration and the number of fragmented moisture niches in the pore network, particularly on the drainage curve. In contrast, the wetting curve exhibited greater variability in spatial configuration due to the “ink bottle effect,” where capillarity was interrupted by wider pore bodies. Ultimately, the methodology provides realistic pore-scale moisture conditions that can be easily manipulated and maintained, enabling new opportunities to explore soil biophysics and microbial biogeography in unsaturated porous media. As a brief example, images showing the localization of fluorescently tagged Pantoea sp. YR343 at −4.3 kPa are presented, highlighting bacterial distributions in water films and air-water interfaces.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Universal Polarization Transformations: Spatial Programming of Polarization Scattering Matrices Using a Deep Learning‐Designed Diffractive Polarization Transformer

Abstract Controlled synthesis of optical fields having nonuniform polarization distributions presents a challenging task. Here, a universal polarization transformer is demonstrated that can synthesize a large set of arbitrarily‐selected, complex‐valued polarization scattering matrices between the polarization states at different positions within its input and output field‐of‐views (FOVs). This framework comprises 2D arrays of linear polarizers positioned between isotropic diffractive layers, each containing tens of thousands of diffractive features with optimizable transmission coefficients. After its deep learning‐based training, this diffractive polarization transformer can successfully implement N i N o = 10 000 different spatially‐encoded polarization scattering matrices with negligible error, where N i and N o represent the number of pixels in the input and output FOVs, respectively. This universal polarization transformation framework is experimentally validated in the terahertz spectrum by fabricating wire‐grid polarizers and integrating them with 3D‐printed diffractive layers to form a physical polarization transformer. Through this set‐up, an all‐optical polarization permutation operation of spatially‐varying polarization fields is demonstrated, and distinct spatially‐encoded polarization scattering matrices are simultaneously implemented between the input and output FOVs of a compact diffractive processor. This framework opens up new avenues for developing novel devices for universal polarization control and may find applications in, e.g., remote sensing, medical imaging, security, material inspection, and machine vision.

Optical neural networks↗
DOE OSTI2024

Reduced-order model to approximate response matrices for filter stack spectrometers

We present a reduced-order model to calculate response matrices rapidly for filter stack spectrometers (FSSs). The reduced-order model allows response matrices to be built modularly from a set of pre-computed photon and electron transport and scattering calculations through various filter and detector materials. While these modular response matrices are not appropriate for high-fidelity analysis of experimental data, they encode sufficient physics to be used as a forward model in design optimization studies of FSSs, particularly for machine learning approaches that require sampling and testing a large number of FSS designs.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Perturbation theory approach to predict the covariance matrices of the galaxy power spectrum and bispectrum in redshift space

In this paper, we predict the covariance matrices of both the power spectrum and the bispectrum, including full non-Gaussian contributions, redshift space distortions, linear bias effects, and shot-noise corrections, using perturbation theory (PT). To quantify the redshift-space distortion effect, we focus mainly on the monopole and quadrupole components of both the power and bispectra. We, for the first time, compute the 5- and 6-point spectra to predict the cross-covariance between the power and bispectra, and the autocovariance of the bispectrum in redshift space. Further, we test the validity of our calculations by comparing them with the covariance matrices measured from the MultiDark-Patchy mock catalogues that are designed to reproduce the galaxy clustering measured from the Baryon Oscillation Spectroscopic Survey Data Release 12. We argue that the simple, leading-order PT works because the shot-noise corrections for the Patchy mocks are more dominant than other higher order terms we ignore. In the meantime, we confirm some discrepancies in the comparison, especially of the cross-covariance. We discuss potential sources of such discrepancies. We also show that our PT model reproduces well the cumulative signal-to-noise ratio of the power spectrum and the bispectrum as a function of maximum wavenumber, implying that our PT model captures successfully essential contributions to the covariance matrices.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Validation of semi-analytical, semi-empirical covariance matrices for two-point correlation function for early DESI data

ABSTRACT We present an extended validation of semi-analytical, semi-empirical covariance matrices for the two-point correlation function (2PCF) on simulated catalogs representative of luminous red galaxies (LRGs) data collected during the initial 2 months of operations of the Stage-IV ground-based Dark Energy Spectroscopic Instrument (DESI). We run the pipeline on multiple effective Zel’dovich (EZ) mock galaxy catalogs with the corresponding cuts applied and compare the results with the mock sample covariance to assess the accuracy and its fluctuations. We propose an extension of the previously developed formalism for catalogs processed with standard reconstruction algorithms. We consider methods for comparing covariance matrices in detail, highlighting their interpretation and statistical properties caused by sample variance, in particular, non-trivial expectation values of certain metrics even when the external covariance estimate is perfect. With improved mocks and validation techniques, we confirm a good agreement between our predictions and sample covariance. This allows one to generate covariance matrices for comparable data sets without the need to create numerous mock galaxy catalogs with matching clustering, only requiring 2PCF measurements from the data itself. The code used in this paper is publicly available at https://github.com/oliverphilcox/RascalC.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Skew-Symmetric adjacency matrices for clustering directed graphs

Graph clustering methods often critically rely on the symmetry of graph matrices. Developing analogous methods for digraphs often proves more challenging, because digraph matrices are typically asymmetric and not orthogonally diagonalizable. However, researchers have recently proposed several complex-valued Hermitian digraph matrices. In particular, one such representation has been utilized as an input to an algorithm for finding imbalanced cuts. In this work, we establish an algebraic relationship between this matrix and an associated real-valued matrix. We show using this real-valued matrix for imbalanced cut-finding algorithms is not only sufficient but advantageous. Our algorithm uses less memory and asymptotically less computation while provably preserving solution quality. We also show our method can be easily implemented using standing computational building blocks, possesses better numerical properties, and loans itself to a natural interpretation via an objective function relaxation argument. We empirically demonstrate these advantages on real world data sets and show how our algorithm can uncover meaningful cluster structure.

Hayashi, Koby↗
DOE OSTI2022

Controlled precision QUBO-based algorithm to compute eigenvectors of symmetric matrices

We describe an algorithm to compute the extremal eigenvalues and corresponding eigenvectors of a symmetric matrix which is based on solving a sequence of Quadratic Binary Optimization problems. This algorithm is robust across many different classes of symmetric matrices; It can compute the eigenvector/eigenvalue pair to essentially any arbitrary precision, and with minor modifications, can also solve the generalized eigenvalue problem. Performance is analyzed on small random matrices and selected larger matrices from practical applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Implementation of a feature selection algorithm in FARM to identify important state variables and time-invariant matrices

The FARM (Feasible Actuator Range Modifier) software module is a component of the RAVEN-based FORCE framework for analysis of Integrated Energy Systems (IES). FARM aids the HERON software module in the evaluation of the optimal dispatch for the different IES components. Set-point trajectories are required to meet limits on both production variables (i.e., the variables to be optimized such as the electrical power, the hydrogen production rate, etc.) and process variables tied to the service life of equipment (e.g., steam flowrate, vessel pressure, turbine firing temperature, etc.). To evaluate the feasibility of HERON generated set-points and to do so in an acceptable time, FARM employs reduced order models to represent the dynamic behavior of the systems to be dispatched. These surrogate models take the form of a linear dynamic system with sets of Linear Parameter Varying (LPV) matrices that are mapped to the system operating space. These matrices are derived from the trajectories of system state variables and system output variables during transients. The accuracy of LPV matrices depends on the selection of state variables. In previous reports, state variables were selected by adopting a complicated workflow requiring multiple software licenses and an advanced level of user expertise. In this report, a new workflow that automates the state variable selection process is presented. It significantly reduces the frequency of user interventions and does not require multiple software licenses. Each module in the new workflow is described in detail, and the input / output examples in each step of the workflow are provided. It was demonstrated that this workflow can greatly reduce the complexity of the state variable selection process, and that the updated FARM-Gamma and FARM-Delta validators can benefit from this workflow when solving the power dispatch problem of a representative IES test case. Finally, some code improvements that can further enhance the efficiency are suggested.

42 ENGINEERING↗
NASA NTRS1974

Stiffness and mass matrices for shells of revolution (SAMMSOR II)

Utilizing element properties, structural stiffness and mass matrices are generated for as many as twenty harmonics and stored on magnetic tape. Matrices generated constitute input data to be used by other stiffness of revolution programs. Variety of boundary and loading conditions can be employed without having to create new mass and stiffness matrices for each case.

Tillerson, J. R.↗
NASA NTRS1976

An algorithm for the weighting matrices in the sampled-data optimal linear regulator problem

The sampled-data optimal linear regulator problem provides a means whereby a control designer can use an understanding of continuous optimal regulator design to produce a digital state variable feedback control law which satisfies continuous system performance specifications. A basic difficulty in applying the sampled-data regulator theory is the requirement that certain digital performance index weighting matrices, expressed as complicated functions of system matrices, be computed. Infinite series representations are presented for the weighting matrices of the time-invariant version of the optimal linear sampled-data regulator problem. Error bounds are given for estimating the effect of truncating the series expressions after a finite number of terms, and a method is described for their computer implementation. A numerical example is given to illustrate the results.

Armstrong, E. S.↗
NASA NTRS1977

On the representation matrices of the spin permutation group

A method is presented for the determination of the representation matrices of the spin permutation group (symmetric group), a detailed knowledge of these matrices being required in the study of the electronic structure of atoms and molecules. The method is characterized by the use of two different coupling schemes. Unlike the Yamanouchi spin algebraic scheme, the method is not recursive. The matrices for the fundamental transpositions can be written down directly in one of the two bases. The method results in a computationally significant reduction in the number of matrix elements that have to be stored when compared with, say, the standard Young tableaux group theoretical approach.

Wilson, S.↗
NASA NTRS1980

Graphite composites with advanced resin matrices

The effect of processing variables on the flammability and mechanical properties for state-of-the-art and advanced resin matrices for graphite composites were studied. Resin matrices which were evaluated included state-of-the-art epoxy, phenolic-novolac, phenolic-xylok, two types of bismaleimides, benzyl, polyethersulfone, and poly(p-phenylene sulfone). Comparable flammability and thermochemical data on graphite-reinforced laminates prepared with these resin matrices are presented, and the relationship of some of these properties to the anaerobic char yield of the resins is described.

Kourtides, D. A.↗
NASA NTRS1981

Two-dimensional integrating matrices on rectangular grids

The use of integrating matrices in solving differential equations associated with rotating beam configurations is examined. In vibration problems, by expressing the equations of motion of the beam in matrix notation, utilizing the integrating matrix as an operator, and applying the boundary conditions, the spatial dependence is removed from the governing partial differential equations and the resulting ordinary differential equations can be cast into standard eigenvalue form. Integrating matrices are derived based on two dimensional rectangular grids with arbitrary grid spacings allowed in one direction. The derivation of higher dimensional integrating matrices is the initial step in the generalization of the integrating matrix methodology to vibration and stability problems involving plates and shells.

Lakin, W. D.↗
NASA NTRS1981

Transfer-matrices for series-type microwave antenna circuits

Transfer matrices are developed which permit analysis and computer evaluation of certain series type microwave antenna circuits associated with an L-Band microwave radiometer (LBMR) under investigation at Goddard Space Flight Center. This radiometer is one of several diverse instrument designs to be used for the determination of soil moisture, sea state, salinity, and temperature data. Four port matrix notation is used throughout for the evaluation of LBMR circuits with mismatched couplers and lossy transmission lines. Matrix parameters in examples are predicted on an impedance analysis and an assumption of an array aperture distribution. The notation presented is easily adapted to longer and more varied chains of matrices, and to matrices of larger dimension.

Schmidt, R. F.↗
NASA NTRS1984

Processing and flammability parameters of bismaleimide and some other thermally stable resin matrices for composites

The effect of processing variables on the flammability and mechanical properties of state-of-the-art and advanced resin matrices for graphite composites were studied. Resin matrices which were evaluated included state-of-the-art epoxy, phenolic-novolac, phenolic-Xylok, two types of bismaleimides, benzyl, polyethersulfone, and polyphenylsulfone. Comparable flammability and thermochemical data on graphite-reinforced laminates prepared with these resin matrices are presented, and the relationship of some of these properties to the anaerobic char yield of the resins is described.

Kourtides, D. A.↗
NASA NTRS1985

Experimental measurements of hydrodynamic radial forces and stiffness matrices for a centrifugal pump-impeller

Measurements of the steady-state hydrodynamic forces on a centrifugal pump impeller are presented as a function of position within two geometrically different volutes. These correspond to the forces experienced by the impeller at zero whirl frequency. The hydrodynamic force matrices derived from these measurements exhibit both diagonal and off-diagonal terms of substantial magnitude. These terms are of the form which would tend to excite a whirl motion in a rotordynamic analysis of the pump; this may be the cause of 'rough running' reported in many pumps. Static pressure measurements in the impeller discharge flow show that the hydrodynamic force on the impeller contains a substantial component due to the nonisotropy of the net momentum flux leaving the impeller. A similar breakdown of the contributions to the stiffness matrices reveals that the major component of these matrices results from the nonisotropy of the momentum flux.

Chamieh, D. S.↗
NASA NTRS1986

Preconditioning matrices for Chebyshev derivative operators

The problem of preconditioning the matrices arising from pseudo-spectral Chebyshev approximations of first order operators is considered in both one and two dimensions. In one dimension a preconditioner represented by a full matrix which leads to preconditioned eigenvalues that are real, positive, and lie between 1 and pi/2, is already available. Since there are cases in which it is not computationally convenient to work with such a preconditioner, a large number of preconditioners were studied which were more sparse (in particular three and four diagonal matrices). The eigenvalues of such preconditioned matrices are compared. The results were applied to the problem of finding the steady state solution to an equation of the type u sub t = u sub x + f, where the Chebyshev collocation is used for the spatial variable and time discretization is performed by the Richardson method. In two dimensions different preconditioners are proposed for the matrix which arises from the pseudo-spectral discretization of the steady state problem. Results are given for the CPU time and the number of iterations using a Richardson iteration method for the unpreconditioned and preconditioned cases.

Rothman, Ernest E.↗
NASA NTRS1987

Approximations to eigenvalues of modified general matrices

The reanalysis of non-self-adjoint dynamic models is computationally very expensive in design optimization applications. This paper describes several approximations that can be applied to eigenvalues of non-hermitian matrices to reduce that computational cost. Approximations based on eigenvalue derivatives, generalized Rayleigh quotient and the trace theorem are presented and their accuracy and computational cost are estimated. The accuracy and cost estimates are verified by applying the approximations to random matrices and matrices arising in flutter analysis of compressor blades. Recommendations are made for selection of the best approximation when the derivatives are available and when they are not. In particular, it is concluded that the quadratic approximation for eigenvalues should never be used as higher order approximations are always more accurate as well as more efficient.

Murthy, Durbha V.↗