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 55 records · Page 3

An analytical technique for approximating unsteady aerodynamics in the time domain

An analytical technique is presented for approximating unsteady aerodynamic forces in the time domain. The order of elements of a matrix Pade approximation was postulated, and the resulting polynomial coefficients were determined through a combination of least squares estimates for the numerator coefficients and a constrained gradient search for the denominator coefficients which insures stable approximating functions. The number of differential equations required to represent the aerodynamic forces to a given accuracy tends to be smaller than that employed in certain existing techniques where the denominator coefficients are chosen a priori. Results are shown for an aeroelastic, cantilevered, semispan wing which indicate a good fit to the aerodynamic forces for oscillatory motion can be achieved with a matrix Pade approximation having fourth order numerator and second order denominator polynomials.

Dunn, H. J.↗

Revisiting the Hail Radar Reflectivity–Kinetic Energy Flux Relation by Combining T-Matrix and Discrete Dipole Approximation Calculations to Size Distribution Observations

The retrieval of hail kinetic energy with weather radars or its simulation in numerical models is challenging because of the shape complexity and variable density of hailstones. We combine 3D scans of individual hailstones with measurements of the particle size distributions (PSD) and T-matrix calculations to understand how hail reflectivity Z changes when approximating hailstones as spheroids, as compared to the realistic shapes obtained by 3D scanning technology. Additionally, recent terminal velocity relations are used to compare Z to the hail kinetic energy flux E ˙ . We parameterize the hail backscattering cross sections at L, S, C, and X bands as a function of size between 0.5 and 5.0 cm, matching the range of the observed PSDs. The scattering calculations use the T-matrix method for size parameters below 1.0 and the discrete dipole approximation (DDA) method otherwise. The DDA calculations are done for 48 digital models of realistic hailstones of sizes between 1 and 5 cm. The DDA cross sections are calculated for multiple orientations and averaged assuming a fully random orientation distribution to provide a single value per hailstone. The T-matrix reflectivity assuming solid ice spheres presents negligible differences to DDA results for size parameters below 1.0. Therefore, T matrix was used to fill in the gaps left by the DDA calculations. The results are mapped to the same size bins of the observed PSDs, allowing the calculation of the radar reflectivity. This is then correlated to E ˙ , allowing a potential improvement of past retrieval methods of E ˙ from Z in multiple wavelengths.

54 ENVIRONMENTAL SCIENCES↗

Randomized Algorithms for Symmetric Nonnegative Matrix Factorization

Symmetric Nonnegative Matrix Factorization (SymNMF) is a technique in data analysis and machine learning that approximates a matrix with a product of a nonnegative, low-rank matrix and it transpose. To design faster and more scalable algorithms for SymNMF we develop two randomized algorithms for its computation. The first method uses randomized matrix sketching to compute an initial low-rank approximation to the input matrix and proceeds to uses this as a low-rank input to rapidly compute a SymNMF. The second methods uses randomized leverage score sampling to approximately solve constrained least squares problems. Many successful methods for SymNMF rely on (approximately) solving sequences of constrained least squares problems. Here, we prove theoretically that leverage score sampling can approximately solve constrained least squares problems to e-accuracy. Finally we demonstrate both methods work in practice by applying them to graph clustering tasks on large real world data sets. These experiments show that our methods approximately maintain solution quality and achieve significant speed ups for both large dense and large sparse problems.

97 MATHEMATICS AND COMPUTING↗

A computational procedure for automated flutter analysis

A direct solution procedure for computing the flutter Mach number and the flutter frequency is applied to the aeroelastic analysis of propfans using a finite element structural model and an unsteady aerodynamic model based on a three-dimensional subsonic compressible lifting surface theory. An approximation to the Jacobian matrix that improves the efficiency of the iterative process is presented. The Jacobian matrix is indirectly approximated from approximate derivatives of the flutter matrix. Examples are used to illustrate the convergence properties. The direct solution procedure facilitates the automated flutter analysis in addition to contributing to the efficient use of computer time as well as the analyst's time.

Murthy, Durbha V.↗

A computational procedure for automated flutter analysis

A direct solution procedure for computing the flutter Mach number and the flutter frequency is applied to the aeroelastic analysis of propfans using an unsteady aerodynamic model based on a three-dimensional subsonic compressible lifting surface theory. An approximation to the Jacobian matrix that improves the efficiency of the iterative process is presented. The Jacobian matrix is indirectly approximated from approximate derivatives of the flutter matrix. Examples are used to illustrate the convergence properties. The direct solution procedure facilitates the automated flutter analysis in addition to contributing to the efficient use of computer time as well as the analyst's time.

Murthy, Durbha V.↗

A computational procedure for automated flutter analysis

A direct solution procedure for computing the flutter Mach number and the flutter frequency is applied to the aeroelastic analysis of propfans using a finite element structural model and an unsteady aerodynamic model based on a three-dimensional subsonic compressible lifting surface theory. An approximation to the Jacobian matrix that improves the efficiency of the iterative process is presented. The Jacobian matrix is indirectly approximated from approximate derivatives of the flutter matrix, which are updated only in the direction of the last move. Examples are used to illustrate the convergence properties. The direct solution procedure facilitates the automated flutter analysis in addition to contributing to the efficient use of computer time as well as the analyst's time.

Murthy, Durbha V.↗

Alteration in CM carbonaceous chondrites inferred from modal and chemical variations in matrix

The modal abundance of matrix in CM chondrites appears to vary from 57-85 vol%. The concentrations of volatile elements that should occur in matrix remain approximately constant despite differences in the proportions of matrix, suggesting that the differing matrix contents may not be real primary variations but are optical effects due to aqueous alteration processes that make other petrologic components unrecognizable. Apparent matrix content can be used as a qualitative measure of the degree of alteration experienced by each CM chondrite. Fe/Si ratios in matrices decrease progressively with increasing alteration due to the formation of new phyllosilicate phases with higher Mg/Fe ratios and optically recognizable opaque minerals that are not counted as matrix. The aqueous alteration process in CM chondrites appears to have been largely isochemical if the bulk meteorites are considered as the reacting systems, although depletion patterns and isotopic anomalies indicate open-system behavior for a few highly mobile components.

Mcsween, H. Y., Jr.↗

Progress on a generalized coordinates tensor product finite element 3DPNS algorithm for subsonic

A generalized coordinates form of the penalty finite element algorithm for the 3-dimensional parabolic Navier-Stokes equations for turbulent subsonic flows was derived. This algorithm formulation requires only three distinct hypermatrices and is applicable using any boundary fitted coordinate transformation procedure. The tensor matrix product approximation to the Jacobian of the Newton linear algebra matrix statement was also derived. Tne Newton algorithm was restructured to replace large sparse matrix solution procedures with grid sweeping using alpha-block tridiagonal matrices, where alpha equals the number of dependent variables. Numerical experiments were conducted and the resultant data gives guidance on potentially preferred tensor product constructions for the penalty finite element 3DPNS algorithm.

Baker, A. J.↗

Integrated Circuit For Simulation Of Neural Network

Ballast resistors deposited on top of circuit structure. Cascadable, programmable binary connection matrix fabricated in VLSI form as basic building block for assembly of like units into content-addressable electronic memory matrices operating somewhat like networks of neurons. Connections formed during storage of data, and data recalled from memory by prompting matrix with approximate or partly erroneous signals. Redundancy in pattern of connections causes matrix to respond with correct stored data.

Thakoor, Anilkumar P.↗

Design for active and passive flutter suppression and gust alleviation

Analytical design techniques for active and passive control of aeroelastic systems are based on a rational approximation of the unsteady aerodynamic loads in the entire Laplace domain, which yields matrix equations of motion with constant coefficients. Some existing schemes are reviewed, the matrix Pade approximant is modified, and a technique which yields a minimal number of augmented states for a desired accuracy is presented. The state-space aeroelastic model is used to design an active control system for simultaneous flutter suppression and gust alleviation. The design target is for a continuous controller which transfers some measurements taken on the vehicle to a control command applied to a control surface. Structural modifications are formulated in a way which enables the treatment of passive flutter suppression system with the same procedures by which active control systems are designed.

Karpel, M.↗

Applications of Laplace transform methods to airfoil motion and stability calculations

This paper reviews the development of generalized unsteady aerodynamic theory and presents a derivation of the generalized Possio integral equation. Numerical calculations resolve questions concerning subsonic indicial lift functions and demonstrate the generation of Kutta waves at high values of reduced frequency, subsonic Mach number, or both. The use of rational function approximations of unsteady aerodynamic loads in aeroelastic stability calculations is reviewed, and a reformulation of the matrix Pade approximation technique is given. Numerical examples of flutter boundary calculations for a wing which is to be flight tested are given. Finally, a simplified aerodynamic model of transonic flow is used to study the stability of an airfoil exposed to supersonic and subsonic flow regions.

Edwards, J. W.↗

Synthesis of Zeolites LiX and LiAgX

In this work, we produced approximately 50 grams of zeolite LiX and 50 grams of zeolite LiAgX from commercially available zeolite X. The production of the two specialized zeolite phases was accomplished via a series of ion exchange steps followed by washing, drying, and calcining to produce the final materials. The ion exchange was confirmed by several methods including inductively-coupled plasma emission spectroscopy of exchange solution aliquots, gravimetric analysis, energy dispersive x-ray spectroscopy, and X-ray diffraction. We synthesized approximately 115 grams of zeolite LiX in a large-scale production step after completing a proof-of-concept experiment demonstrating the methodology. The production of zeolite LiX from zeolite X yielded an ion exchange (Li + →Na + ) of approximately 82%. The production of zeolite LiAgX from the previously synthesized zeolite LiX was targeted to exchange approximately 20% of the Li + ions with Ag + ions. We first performed a proof-of-concept experiment demonstrating the reaction pathway and selectivity of the Ag + ion, and then performed a large-scale synthesis of zeolite LiAgX using approximately 50 grams of the zeolite LiX. Our final exchange percentage of silver into the zeolite LiX matrix was approximately 20.8%.

36 MATERIALS SCIENCE↗

Program Determines Minimum-State Approximations To Unsteady Aerodynamic Forces

MIST implements minimum-state method for determining rational approximations of matrix of aerodynamic force coefficients. Accepts complex tabular data representing generalized unsteady aerodynamic forces over set of reduced frequencies. Determines approximations to tabular data in Laplace domain by use of rational functions. Provides capability to select coefficients of denominator in rational approximations, to constrain approximation selectably without increasing size of problem, and determines and emphasizes criticial frequency ranges in determining approximations. Written in ANSI FORTRAN 77.

Karpel, Mordechay↗

Proposed framework for thermomechanical life modeling of metal matrix composites

The framework of a mechanics of materials model is proposed for thermomechanical fatigue (TMF) life prediction of unidirectional, continuous-fiber metal matrix composites (MMC's). Axially loaded MMC test samples are analyzed as structural components whose fatigue lives are governed by local stress-strain conditions resulting from combined interactions of the matrix, interfacial layer, and fiber constituents. The metallic matrix is identified as the vehicle for tracking fatigue crack initiation and propagation. The proposed framework has three major elements. First, TMF flow and failure characteristics of in situ matrix material are approximated from tests of unreinforced matrix material, and matrix TMF life prediction equations are numerically calibrated. The macrocrack initiation fatigue life of the matrix material is divided into microcrack initiation and microcrack propagation phases. Second, the influencing factors created by the presence of fibers and interfaces are analyzed, characterized, and documented in equation form. Some of the influences act on the microcrack initiation portion of the matrix fatigue life, others on the microcrack propagation life, while some affect both. Influencing factors include coefficient of thermal expansion mismatch strains, residual (mean) stresses, multiaxial stress states, off-axis fibers, internal stress concentrations, multiple initiation sites, nonuniform fiber spacing, fiber debonding, interfacial layers and cracking, fractured fibers, fiber deflections of crack fronts, fiber bridging of matrix cracks, and internal oxidation along internal interfaces. Equations exist for some, but not all, of the currently identified influencing factors. The third element is the inclusion of overriding influences such as maximum tensile strain limits of brittle fibers that could cause local fractures and ensuing catastrophic failure of surrounding matrix material. Some experimental data exist for assessing the plausibility of the proposed framework.

Halford, Gary R.↗

Towards excitations and dynamical quantities in correlated lattices with density matrix embedding theory

Density matrix embedding theory (DMET) provides a framework to describe ground-state expectation values in strongly correlated systems, but its extension to dynamical quantities is still an open problem. We show one route to obtaining excitations and dynamical spectral functions by using the techniques of DMET to approximate the matrix elements that arise in a single-mode inspired excitation ansatz. We demonstrate this approach in the one-dimensional Hubbard model, comparing the neutral excitations, single-particle density of states, charge, and spin dynamical structure factors to benchmarks from the Bethe ansatz and density matrix renormalization group. Finally, our work highlights the potential of these ideas in building computationally efficient approaches for dynamical quantities.

1-dimensional systems↗