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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

An efficient approximate factorization implicit scheme for the equations of gasdynamics

An efficient implicit finite-difference algorithm for the gas dynamic equations utilizing matrix reduction techniques is presented. A significant reduction in arithmetic operations is achieved while maintaining the same favorable stability characteristics and generality found in the Beam and Warming approximate factorization algorithm. Steady-state solutions to the conservative Euler equations in generalized coordinates are obtained for transonic flows about a NACA 0012 airfoil. The theoretical extension of the matrix reduction technique to the full Navier-Stokes equations in Cartesian coordinates is presented in detail. Linear stability, using a Fourier stability analysis, is demonstrated and discussed for the one-dimensional Euler equations. It is shown that the method offers advantages over the conventional Beam and Warming scheme and can retrofit existing Beam and Warming codes with minimal effort.

Barth, T. J.↗
NASA NTRS1971

Fundamental Mechanisms of Tensile Fracture in Aluminum Sheet Unidirectionally Reinforced with Boron Filament

Results are presented from an experimental research effort to gain a more complete understanding of the physics of tensile fracture in unidirectionally reinforced B-Al composite sheet. By varying the degree of filament degradation resulting from fabrication, composite specimens were produced which failed in tension by the cumulative mode, the noncumulative mode, or by any desired combination of the two modes. Radiographic and acoustic emission techniques were combined to identify and physically describe a previously unrecognized fundamental fracture mechanism which was responsible for the noncumulative mode. The tensile strength of the composite was found to be severely limited by the noncumulative mechanism which involved the initiation and sustenance of a chain reaction of filament fractures at a relatively low stress level followed by ductile fracture of the matrix. The minimum average filament stress required for initiation of the fracture mechanism was shown to be approximately 170 ksi, and appeared to be independent of filament diameter, number of filament layers, and the identity of the matrix alloy.

Herring, H. W.↗
NASA NTRS1991

Modeling the spatiotemporal organization of velocity storage in the vestibuloocular reflex by optokinetic studies

1. A generalized three-dimensional state space model of visual vestibular interaction was developed. Matrix and dynamical system operators associated with inputs from the semicircular canals, otolith velocity estimator, and the visual system have been incorporated into the model, which focus on their relationship to the velocity storage integrator. 2. A relationship was postulated between the eigenvalues and the direction of the eigenvectors of the system matrix and the orientation of the spatial vertical. It was assumed that the system matrix for a tilted position was a composition of two linear transformations of the system matrix for the upright position. One transformation modifies the eigenvalues of the system matrix, whereas another rotates the eigenvectors of the system matrix. The pitch and roll eigenvectors rotate with the head, whereas the yaw axis eigenvector remains approximately spatially invariant. 3. Based on the three-dimensional model, a computational procedure was formulated to identify the eigenvalues and eigenvectors of the system matrix with the use of a modification of the marquardt algorithm. With the use of data obtained from a monkey, it was shown that the three-dimensional behavior of velocity storage cannot be predicted solely in terms of its time constants, i.e., the inverse of its eigenvalues. With the use of the same eigenvalues the data could either be fit or not fit, depending on the eigenvector directions. Therefore, it is necessary to specify eigenvector directions when characterizing velocity storage in three dimensions. 4. Parameters found with the use of the Marquardt algorithm were incorporated into the model. Diagonal matrices in a head coordinate frame were introduced for coupling the visual system to the integrator and to the direct optokinetic pathway. Simulations of optokinetic nystagmus (OKN) and optokinetic after-nystagmus (OKAN) were run. The model predicted the behavior of yaw and pitch OKN and OKAN when the animal is upright. It also predicted the cross-coupling in the side down position. The trajectories in velocity space were also accurately simulated. 5. One of the predictions of the model is that when the stimulus direction is along an eigenvector, the trajectory in velocity space is a straight line. Using the "spectral width" of the residuals from a straight line sequence during OKAN, we developed a methodology to estimate how close the OKAN decay was to an eigenvector trajectory. 6. Thus we have developed a model-based approach for studying and interpreting the response characteristics of velocity storage in three dimensions.(ABSTRACT TRUNCATED AT 400 WORDS).

Non-NASA Center↗
NASA NTRS1974

Solution of an eigenvalue problem for the Laplace operator on a spherical surface

Methods for obtaining approximate solutions for the fundamental eigenvalue of the Laplace-Beltrami operator (also referred to as the membrane eigenvalue problem for the vibration equation) on the unit spherical surface are developed. Two specific types of spherical surface domains are considered: (1) the interior of a spherical triangle, i.e., the region bounded by arcs of three great circles, and (2) the exterior of a great circle arc extending for less than pi radians on the sphere (a spherical surface with a slit). In both cases, zero boundary conditions are imposed. In order to solve the resulting second-order elliptic partial differential equations in two independent variables, a finite difference approximation is derived. The symmetric (generally five-point) finite difference equations that develop are written in matrix form and then solved by the iterative method of point successive overrelaxation. Upon convergence of this iterative method, the fundamental eigenvalue is approximated by iteration utilizing the power method as applied to the finite Rayleigh quotient.

Walden, H.↗
NASA NTRS1993

Accuracy of least-squares methods for the Navier-Stokes equations

Recently there has been substantial interest in least-squares finite element methods for velocity-vorticity-pressure formulations of the incompressible Navier-Stokes equations. The main cause for this interest is the fact that algorithms for the resulting discrete equations can be devised which require the solution of only symmetric, positive definite systems of algebraic equations. On the other hand, it is well-documented that methods using the vorticity as a primary variable often yield very poor approximations. Thus, here we study the accuracy of these methods through a series of computational experiments, and also comment on theoretical error estimates. It is found, despite the failure of standard methods for deriving error estimates, that computational evidence suggests that these methods are, at the least, nearly optimally accurate. Thus, in addition to the desirable matrix properties yielded by least-squares methods, one also obtains accurate approximations.

Bochev, Pavel B.↗
NASA NTRS1993

Accuracy of least-squares methods for the Navier-Stokes equations

Recently there has been substantial interest in least-squares finite element methods for velocity-vorticity-pressure formulations of the incompressible Navier-Stokes equations. The main cause for this interest is the fact that algorithms for the resulting discrete equations can be devised which require the solution of only symmetric, positive definite systems of algebraic equations. On the other hand, it is well-documented that methods using the vorticity as a primary variable often yield very poor approximations. Thus, here we study the accuracy of these methods through a series of computational experiments, and also comment on theoretical error estimates. It is found, despite the failure of standard methods for deriving error estimates, that computational evidence suggests that these methods are, at the least, nearly optimally accurate. Thus, in addition to the desirable matrix properties yielded by least-squares methods, one also obtains accurate approximations.

Bochev, Pavel B.↗
NASA NTRS1991

The asymptotic spectra of banded Toeplitz and quasi-Toeplitz matrices

Toeplitz matrices occur in many mathematical, as well as, scientific and engineering investigations. This paper considers the spectra of banded Toeplitz and quasi-Toeplitz matrices with emphasis on non-normal matrices of arbitrarily large order and relatively small bandwidth. These are the type of matrices that appear in the investigation of stability and convergence of difference approximations to partial differential equations. Quasi-Toeplitz matrices are the result of non-Dirichlet boundary conditions for the difference approximations. The eigenvalue problem for a banded Toeplitz or quasi-Toeplitz matrix of large order is, in general, analytically intractable and (for non-normal matrices) numerically unreliable. An asymptotic (matrix order approaches infinity) approach partitions the eigenvalue analysis of a quasi-Toeplitz matrix into two parts, namely the analysis for the boundary condition independent spectrum and the analysis for the boundary condition dependent spectrum. The boundary condition independent spectrum is the same as the pure Toeplitz matrix spectrum. Algorithms for computing both parts of the spectrum are presented. Examples are used to demonstrate the utility of the algorithms, to present some interesting spectra, and to point out some of the numerical difficulties encountered when conventional matrix eigenvalue routines are employed for non-normal matrices of large order. The analysis for the Toeplitz spectrum also leads to a diagonal similarity transformation that improves conventional numerical eigenvalue computations. Finally, the algorithm for the asymptotic spectrum is extended to the Toeplitz generalized eigenvalue problem which occurs, for example, in the stability of Pade type difference approximations to differential equations.

Beam, Richard M.↗
NASA NTRS1979

A new approach to the solution of large, full matrix equations: A two-dimensional potential flow feasibility study

An approach to the solution of matrix problems resulting from integral equations of mathematical physics is presented. Based on the inherent smoothness in such equations, the problem is reformulated using a set of orthogonal basis vectors, leading to an equivalent coefficient problem which can be of lower order without significantly impairing the accuracy of the solution. This approach was evaluated using a two-dimensional Neumann problem describing the inviscid, incompressible flow over an airfoil. Two different kinds of mode functions were investigated, namely eigenfunction series and Fourier series. The method using Fourier series was found preferable. It uses all of the coefficients from a Fast Fourier Transform algorithm in an approximate method which exploits the known structure of the transformed coefficient matrix and very promising results for the flow over a realistic airfoil are obtained. On the basis of the results presented here, an order of magnitude reduction in this computer time can be expected for such problems as compared with the time for a direct matrix solution.

James, R. M.↗
NASA NTRS1989

Atomic data for opacity calculations. X - Oscillator strengths and photoionisation cross sections for O III

Energy levels and electric dipole radiative transitions were determined for O III. The wavefunctions for the bound and continuum states were derived by solving coupled integrodifferential equations in the close-coupling approximation with the aid of two methods, the R-matrix and linear algebraic methods. Configuration interaction wavefunctions are presented for eight states of the O IV target with configurations of 2s(x)sp(y) (x + y = 3). Oscillator strengths are calculated for the transitions between, and photoionization cross sections from, bound states of O III with configurations of 2s(x)2p(y)nl (for n of not greater than 10 and l of not greater than 3).

Luo, D.↗
NASA NTRS1989

The prediction of pseudo-resonance positions in the Schwinger variational principle

The Schwinger Variational Principle is appled to s-wave electron-hydrogen atom scattering. It is shown computationally, that, consistent with a paper by B. Apagyi, P. Levay, and K. Ladanyi, there are pseudo-resonances at the static exchange level of approximation, but not at the static level. The T-matrix was employed, as well as, the K-matrix version of the Schwinger Principle, with a real Slater basis, and obtained the same results in both. The origin of the pseudo-resonances as resulting from singularities in the separable potential that is effectively employed in the Lippmann-Schwinger equation from which the Schwinger Variational Principle can be derived. The determination of the pseudo-resonance parameters from the separable potential is computationally inexpensive and may be used to predict the pseudo-resonance parameters for the scattering calculations so that they may be avoided.

Weatherford, Charles A.↗
NASA NTRS1990

Matrix-based, finite-difference algorithms for computational acoustics

A compact numerical algorithm is introduced for simulating multidimensional acoustic waves. The algorithm is expressed in terms of a set of matrix coefficients on a three-point spatial grid that approximates the acoustic wave equation with a discretization error of O(h exp 5). The method is based on tracking a local phase variable and its implementation suggests a convenient coordinate splitting along with natural intermediate boundary conditions. Results are presented for oblique plane waves and compared with other procedures. Preliminary computations of acoustic diffraction are also considered.

Davis, Sanford↗
NASA NTRS1989

A study of pseudoresonances in the application of the Schwinger variational principle to electron scattering from atomic hydrogen

The Schwinger variational principle is applied to s-wave electron-hydrogen atom scattering. Computationally, consistent with a recent paper by Apagyi et al. (1988), there are pseudoresonances at the static exchange level of approximation, but not at the static level. The T-matrix as well as the K-matrix version of the Schwinger principle is used with a real Slater basis, and the same results are obtained in both. The origin of the pseudoresonances from singularities in the separable potential that is effectively employed in the Lippman-Schwinger equation from which the Schwinger variational principle can be derived. The determination of the pseudoresonance parameters from the separable potential is computationally inexpensive and may be used to predict the pseudoresonance parameters for the scattering calculations so as to avoid them.

Weatherford, Charles A.↗
NASA NTRS1990

Influence of interfacial shear strength on the mechanical properties of SiC fiber reinforced reaction-bonded silicon nitride matrix composites

The influence of fiber/matrix interface microstructure and interfacial shear strength on the mechanical properties of a fiber-reinforced ceramic composite was evaluated. The composite consisted of approximately 30 vol percent uniaxially aligned 142 microns diameter SiC fibers (Textron SCS-6) in a reaction-bonded Si3N4 matrix (SiC/RBSN). The interface microstructure was varied by controlling the composite fabrication conditions and by heat treating the composite in an oxidizing environment. Interfacial shear strength was determined by the matrix crack spacing method. The results of microstructural examination indicate that the carbon-rich coating provided with the as-produced SiC fibers was stable in composites fabricated at 1200 C in a nitrogen or in a nitrogen plus 4 percent hydrogen mixture for 40 hr. However this coating degraded in composites fabricated at 1350 C in N2 + 4 percent H2 for 40 and 72 hr and also in composites heat treated in an oxidizing environment at 600 C for 100 hr after fabrication at 1200 C in a nitrogen. It was determined that degradation occurred by carbon removal which in turn had a strong influence on interfacial shear strength and other mechanical properties. Specifically, as the carbon coating was removed, the composite interfacial shear strength, primary elastic modulus, first matrix cracking stress, and ultimate tensile strength decreased, but the first matrix cracking strain remained nearly the same.

Bhatt, Ramakrishna T.↗
NASA NTRS2000

Effects of Interface Coating and Nitride Enhancing Additive on Properties of Hi-Nicalon SiC Fiber Reinforced Reaction-Bonded Silicon Nitride Composites

Strong and tough Hi-Nicalon SiC fiber reinforced reaction-bonded silicon nitride matrix composites (SiC/ RBSN) have been fabricated by the fiber lay-up approach. Commercially available uncoated and PBN, PBN/Si-rich PBN, and BN/SiC coated SiC Hi-Nicalon fiber tows were used as reinforcement. The composites contained approximately 24 vol % of aligned 14 micron diameter SiC fibers in a porous RBSN matrix. Both one- and two-dimensional composites were characterized. The effects of interface coating composition, and the nitridation enhancing additive, NiO, on the room temperature physical, tensile, and interfacial shear strength properties of SiC/RBSN matrix composites were evaluated. Results indicate that for all three coated fibers, the thickness of the coatings decreased from the outer periphery to the interior of the tows, and that from 10 to 30 percent of the fibers were not covered with the interface coating. In the uncoated regions, chemical reaction between the NiO additive and the SiC fiber occurs causing degradation of tensile properties of the composites. Among the three interface coating combinations investigated, the BN/SiC coated Hi-Nicalon SiC fiber reinforced RBSN matrix composite showed the least amount of uncoated regions and reasonably uniform interface coating thickness. The matrix cracking stress in SiC/RBSN composites was predicted using a fracture mechanics based crack bridging model.

Bhatt, Ramakrishana T.↗
NASA NTRS1977

4.2-4.3 AE anorthositic soil fragments - Equilibrated or unequilibrated polycomponent systems

The petrographies and phase chemistries of five lunar crystalline anorthositic fragments with Ar-40/Ar-39 plateau 'ages' ranging from 4.1 to 4.3 AE are reported. An investigation is conducted to determine whether the individual fragments are monomict or polymict breccias. For polymict breccias it is attempted to establish whether clasts and matrix have equilibrated chemically. On the basis of the results of the investigation it is suggested that sample 73263,1,1 was a breccia (possibly a clast-laden melt) containing old clasts, incorporated in an ejecta blanket that was annealed at about 4.1 AE in a transient high-temperature event. The sample 73263,1,6 studied is an annealed anorthositic breccia containing approximately 20% clasts dominated by single-crystal plagioclase fragments. The obtained data are consistent with a two-component system containing old plagioclase clasts in a young (about 4.0 AE) annealed, matrix. The annealing event at approximately 4.0 AE was insufficient to equilibrate the argon. It is found that four of the breccias having ages of about 4.2 AE or more are polymict. The most striking feature of all of the samples studied is the presence of a significant amount of 4.2-4.3 AE argon.

Delano, J. W.↗
NASA NTRS2016

K-Shell Photoabsorption and Photoionisation of Trace Elements I. Isoelectronic Sequences With Electron Number 3< or = N < or = 11

Context. With the recent launching of the Hitomi X-ray space observatory, K lines and edges of chemical elements with low cosmic abundances, namely F, Na, P, Cl, K, Sc, Ti, V, Cr, Mn, Co, Cu and Zn, can be resolved and used to determine important properties of supernova remnants, galaxy clusters and accreting black holes and neutron stars.Aims. The second stage of the present ongoing project involves the computation of the accurate photoabsorption and photoionisation cross sections required to interpret the X-ray spectra of such trace elements.Methods. Depending on target complexity and computer tractability, ground-state cross sections are computed either with the close-coupling Breit-Pauli R-matrix method or with the autostructure atomic structure code in the isolated-resonance approximation. The intermediate-coupling scheme is used whenever possible. In order to determine a realistic K-edge behaviour for each species, both radiative and Auger dampings are taken into account, the latter being included in the R-matrix formalism by means of an optical potential.Results. Photoabsorption and total and partial photoionisation cross sections are reported for isoelectronic sequences with electron numbers 3< or = N< or = 11. The Na sequence (N=11) is used to estimate the contributions from configurations with a 2s hole (i.e. [2s]) and those containing 3d orbitals, which will be crucial when considering sequences with N 11.Conclusions. It is found that the [2s/u] configurations must be included in the target representations of species with N> 11 as they contribute significantly to the monotonic background of the cross section between the L and K edges. Configurations with 3d orbitals are important in rendering an accurate L edge, but they can be practically neglected in the K-edge region.

atomic data / X-rays: general↗
NASA NTRS1995

Cytoskeletal filament assembly and the control of cell spreading and function by extracellular matrix

This study was undertaken to analyze how cell binding to extracellular matrix produces changes in cell shape. We focused on the initial process of cell spreading that follows cell attachment to matrix and, thus, cell 'shape' changes are defined here in terms of alterations in projected cell areas, as determined by computerized image analysis. Cell spreading kinetics and changes in microtubule and actin microfilament mass were simultaneously quantitated in hepatocytes plated on different extracellular matrix substrata. The initial rate of cell spreading was highly dependent on the matrix coating density and decreased from 740 microns 2/h to 50 microns 2/h as the coating density was lowered from 1000 to 1 ng/cm2. At approximately 4 to 6 hours after plating, this initial rapid spreading rate slowed and became independent of the matrix density regardless of whether laminin, fibronectin, type I collagen or type IV collagen was used for cell attachment. Analysis of F-actin mass revealed that cell adhesion to extracellular matrix resulted in a 20-fold increase in polymerized actin within 30 minutes after plating, before any significant change in cell shape was observed. This was followed by a phase of actin microfilament disassembly which correlated with the most rapid phase of cell extension and ended at about 6 hours; F-actin mass remained relatively constant during the slow matrix-independent spreading phase. Microtubule mass increased more slowly in spreading cells, peaking at 4 hours, the time at which the transition between rapid and slow spreading rates was observed. However, inhibition of this early rise in microtubule mass using either nocodazole or cycloheximide did not prevent this transition. Use of cytochalasin D revealed that microfilament integrity was absolutely required for hepatocyte spreading whereas interference with microtubule assembly (using nocodazole or taxol) or protein synthesis (using cycloheximide) only partially suppressed cell extension. In contrast, cell spreading could be completely inhibited by combining suboptimal doses of cytochalasin D and nocodazole, suggesting that intact microtubules can stabilize cell form when the microfilament lattice is partially compromised. The physiological relevance of the cytoskeleton and cell shape in hepatocyte physiology was highlighted by the finding that a short exposure (6 hour) of cells to nocodazole resulted in production of smaller cells 42 hours later that exhibited enhanced production of a liver-specific product (albumin). These data demonstrate that spreading and flattening of the entire cell body is not driven directly by net polymerization of either microfilaments or microtubules.(ABSTRACT TRUNCATED AT 400 WORDS).

NASA Discipline Cell Biology↗
NASA NTRS1991

Application of Krylov exponential propagation to fluid dynamics equations

This paper presents an application of matrix exponentiation via Krylov subspace projection, to the solution of fluid dynamics problems. The main idea is to approximate the operation exp(A)v by means of a projection-like process onto a Krylov subspace. This results in a computation of an exponential matrix vector product similar to the one above but of a much smaller size. Time integration schemes can then be devised to exploit this basic computational kernel. The motivation of this approach is to provide time-integration schemes that are essentially of an explicit nature but which have good stability properties.

Saad, Y.↗