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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Algorithm Development for Predicting Biodiversity Based on Phytoplankton Absorption

Ocean color remote sensing has provided the scientific community with unprecedented global coverage of chlorophyll a, an indicator of phytoplankton biomass. Together, satellite-derived chlorophyll a and knowledge of Phytoplankton Functional Types (PFTs) will improve our limited understanding of marine ecosystem responses to physiochemical climate drivers involved in carbon cycle dynamics and linkages. Using cruise data from the Gulf of Maine and the Middle Atlantic Bight (N 269 pairs of HPLC and phytoplankton absorption samples), two modeling approaches were utilized to predict phytoplankton absorption and pigments. Algorithm I predicts the chlorophyll-specific absorption coefficient (a*(sub ph) (m(exp 2) mg chl a(exp - 1)) using inputs of temperature, light, and chlorophyll a. Modeled r(exp 2) values (400700 nm) ranged from 0.79 to 0.99 when compared to in situ observations with approximately 25 percent lower r(exp 2) values in the UV region. Algorithm II-a utilizes matrix inversion analysis to predict a(sub ph)(m(exp -1), 400700 nm) and r(exp 2) values ranged from 0.89 to 0.99. The prediction of phytoplankton pigments with Algorithm II-b produced r(exp 2) values that ranged from 0.40 to 0.93. When used in combination, Algorithm I, and Algorithm II-a are able to use satellite products of SST, PAR, and chlorophyll a (Algorithm I) to predict pigment concentrations and ratios to describe the phytoplankton community. The results of this study demonstrate that the spatial variation in modeled pigment ratios differ significantly from the 10-year SeaWiFS average chlorophyll a data set. Contiguous observations of chlorophyll a and phytoplankton biodiversity will elucidate ecosystem responses with unprecedented complexity.

phytoplankton↗
NASA NTRS1981

New approach to the solution of large, full matrix equations

A new approach to the solution of matrix equations resulting from integral equations is presented and applied to the solution of two-dimensional Neumann problems describing the inviscid, incompressible flow past an airfoil. The problem is reformulated in terms of a preselected set of mode functions giving an equivalent matrix equation to be solved for the mode-function expansion coefficients. Because of the inherent smoothness of the original problem, the coefficient problem can be solved approximately without significantly affecting the accuracy of the final solution. Very promising two-dimensional results are obtained and the extension of the method to three-dimensional problems is investigated. On the basis of these results it is shown that the computing time for the matrix solution for a large three-dimensional panel method calculation could be reduced by an order of magnitude compared with that required for a direct solution.

Clark, R. W.↗
NASA NTRS1993

A study of fiber volume fraction effects in notched unidirectional SCS-6/Ti-15V-3Cr-3Al-3Sn composite

Notched unidirectional SCS-6/Ti-15-3 composite of three different fiber volume fractions (vf = 0.15, 0.37, and 0.41) was investigated for various room temperature microstructural and material properties including: fatigue crack initiation, fatigue crack growth, and fracture toughness. While the matrix hardness is similar for all fiber volume fractions, the fiber/matrix interfacial shear strength and matrix residual stress increases with fiber volume fraction. The composite fatigue crack initiation stress is shown to be matrix controlled and occurs when the net maximum matrix stress approaches the endurance limit stress of the matrix. A model is presented which includes residual stresses and presents the composite initiation stress as a function of fiber volume fraction. This model predicts a maximum composite initiation stress at vf approximately 0.15 which agrees with the experimental data. The applied composite stress levels were increased as necessary for continued crack growth. The applied Delta(K) values at crack arrest increase with fiber volume fraction by an amount better approximated using an energy based formulation rather than when scaled linear with modulus. After crack arrest, the crack growth rate exponents for vf37 and vf41 were much lower and toughness much higher, when compared to the unreinforced matrix, because of the bridged region which parades with the propagating fatigue crack. However, the vf15 material exhibited a higher crack growth rate exponent and lower toughness than the unreinforced matrix because once the bridged fibers nearest the crack mouth broke, the stress redistribution broke all bridged fibers, leaving an unbridged crack. Degraded, unbridged behavior is modeled using the residual stress state in the matrix ahead of the crack tip. Plastic zone sizes were directly measured using a metallographic technique and allow prediction of an effective matrix stress intensity which agrees with the fiber pressure model if residual stresses are considered. The sophisticated macro/micro finite element models of the 0.15 and 0.37 fiber volume fractions presented show good agreement with experimental data and the fiber pressure model when an estimated effective fiber/matrix debond length is used.

Covey, Steven J.↗
NASA NTRS1982

A compliant, high failure strain, fibre-reinforced glass-matrix composite

A glass-matrix composite reinforced by discontinuous graphite fibers was produced by hot pressing glass-powder-impregnated two-dimensional arrays of in-plane randomly oriented graphite fibers held together by approximately 5-10% by weight of organic binder (generally polyester). The composite tensile behavior is characterized by a highly nonlinear stress-strain curve which differs markedly from that of either unreinforced glass or a similarly reinforced epoxy-matrix composite. By virtue of this nonlinearity, the composite is able to redistribute applied stresses to achieve a high load-carrying capacity. The fibrous microstructure and the low fiber-matrix bond provide a mechanism for achieving high fracture toughness and unusually high compliance. For a 96%-silica-matrix composite, the strength is retained to over 1000 C.

Prewo, K. M.↗
NASA NTRS1985

A fast efficient implicit scheme for the gasdynamic equations using a matrix reduction technique

An efficient implicit finite-difference algorithm for the gasdynamic equations utilizing matrix reduction techniques is presented. A significant reduction in arithmetic operations is achieved without loss of the stability characteristics 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 and used to show that the method offers computational advantages over the conventional Beam and Warming scheme. Existing Beam and Warming codes can be retrofit with minimal effort. 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.

Barth, T. J.↗
NASA NTRS1990

A New Approximate Fracture Mechanics Analysis Methodology for Composites with a Crack or Hole

A new approximate theory which links the inherent flaw concept with the theory of crack tip stress singularities at a bi-material interface was developed. Three assumptions were made: (1) the existence of inherent flaw (i.e., damage zone) at the tip of the crack, (2) a fracture of the filamentary composites initiates at a crack lying in the matrix material at the interface of the matrix/filament, and (3) the laminate fails whenever the principal load-carrying laminae fails. This third assumption implies that for a laminate consisting of 0 degree plies, cracks into matrix perpendicular to the 0 degree filaments are the triggering mechanism for the final failure. Based on this theory, a parameter bar K sub Q which is similar to the stress intensity factor for isotropic materials but with a different dimension was defined. Utilizing existing test data, it was found that bar K sub Q can be treated as a material constant. Based on this finding a fracture mechanics analysis methodology was developed. The analytical results are correlated well with test results. This new approximate theory can apply to both brittle and metal matrix composite laminates with crack or hole.

Tsai, H. C.↗
NASA NTRS1967

Matrix Analysis of Longitudinal and Torsional Vibrations in Nonuniform Multibranch Beams

Since longitudinal modes and frequencies provide basic data for dynamic analyses of arbitrary beam-like structures and since closed-form solutions for the modes are generally not feasible to obtain, an approximate method is developed for computing the natural frequencies and the corresponding mode shapes for a variable-section, unconstrained multibranch beam. A lumped mass analogy employing influence coefficients is used to represent the beam. The simultaneous equations of motion for the lumped mass system are derived in matrix form and algebraically manipulated to yield a classical eigenvalue equation solvable by standard procedures. The orthogonality relationship of the natural modes is derived and used to form the basis of an orthogonal sweeping process for determination of modes above the fundamental. Numerical examples including an application to a solid-fuel launch system are presented. Also, a detailed discussion is devoted to the theoretical verifications of the approximate modes and frequencies.

Wingate, Robert T.↗
NASA NTRS1974

Aeroelastic stability of periodic systems with application to rotor blade flutter

The dynamics of a helicopter blade in forward flight are described by a system of linear differential equations with periodic coefficients. The stability of this periodic aeroelastic system is determined, using multivariable Floquet-Liapunov theory. The transition matrix at the end of the period is evaluated by: (1) direct numerical integration, and (2) a new, approximate method, which consists in approximating a periodic function by a series of step functions. The numerical accuracy and efficiency of the methods is compared, and the second method is shown to be superior by far. Results illustrating the effect of the periodic coefficients and various blade parameters are presented.

Friedmann, P.↗
NASA NTRS1977

Problems in characterizing atmospheric aerosols by lidar alone

In this paper we review the scattering properties of real atmospheric aerosols as they relate to laser radar measurements. We address two major subjects: theoretical approximations for modeling lidar experiments and laboratory measurements for characterizing real aerosols. In the first category we review: (1) the Mueller algebra and Stokes vectors, (2) approximations for the single scattering properties of various aerosols, (3) the information content of the Mueller or scattering matrix and related limitations of lidar measurements, and (4) the effects of multiple scattering on lidar returns. In the second category we review: (1) polar nephelometers, (2) calibration procedures for nephelometers, and (3) characterization of the scattering medium. Finally, we summarize the outstanding problems in measuring the atmospheric aerosols by lidar alone.

Holland, A. C.↗
NASA NTRS1980

Steady, Oscillatory, and Unsteady Subsonic and Supersonic Aerodynamics, production version (SOUSSA-P 1.1). Volume 1: Theoretical manual

Recent developments of the Green's function method and the computer program SOUSSA (Steady, Oscillatory, and Unsteady Subsonic and Supersonic Aerodynamics) are reviewed and summarized. Applying the Green's function method to the fully unsteady (transient) potential equation yields an integro-differential-delay equation. With spatial discretization by the finite-element method, this equation is approximated by a set of differential-delay equations in time. Time solution by Laplace transform yields a matrix relating the velocity potential to the normal wash. Premultiplying and postmultiplying by the matrices relating generalized forces to the potential and the normal wash to the generalized coordinates one obtains the matrix of the generalized aerodynamic forces. The frequency and mode-shape dependence of this matrix makes the program SOUSSA useful for multiple frequency and repeated mode-shape evaluations.

Morino, L.↗
NASA NTRS1992

Preliminary Design Of Adhesively Bonded Composite Joints

Joints symmetrical and long enough to withstand cyclic loads and hygrothermal effects. Step-by-step procedure for preliminary design of adhesively bonded joints between laminated fiber-reinforced-matrix composite bars based on simplified methods for prediction of microstresses and local strengths, including interlaminar strengths. Involves approximate calculations of stresses in various components of joint. Includes calculations to determine such critical parameters as minimum length, maximum shear stress in adhesive, and peeloff stress. Deleterious effects of cyclic loading, increased temperature, and humidity also taken into account.

Chamis, C. C.↗
NASA NTRS1998

A Parallel Non-Overlapping Domain-Decomposition Algorithm for Compressible Fluid Flow Problems on Triangulated Domains

This paper considers an algebraic preconditioning algorithm for hyperbolic-elliptic fluid flow problems. The algorithm is based on a parallel non-overlapping Schur complement domain-decomposition technique for triangulated domains. In the Schur complement technique, the triangulation is first partitioned into a number of non-overlapping subdomains and interfaces. This suggests a reordering of triangulation vertices which separates subdomain and interface solution unknowns. The reordering induces a natural 2 x 2 block partitioning of the discretization matrix. Exact LU factorization of this block system yields a Schur complement matrix which couples subdomains and the interface together. The remaining sections of this paper present a family of approximate techniques for both constructing and applying the Schur complement as a domain-decomposition preconditioner. The approximate Schur complement serves as an algebraic coarse space operator, thus avoiding the known difficulties associated with the direct formation of a coarse space discretization. In developing Schur complement approximations, particular attention has been given to improving sequential and parallel efficiency of implementations without significantly degrading the quality of the preconditioner. A computer code based on these developments has been tested on the IBM SP2 using MPI message passing protocol. A number of 2-D calculations are presented for both scalar advection-diffusion equations as well as the Euler equations governing compressible fluid flow to demonstrate performance of the preconditioning algorithm.

Barth, Timothy J.↗
NASA NTRS1991

Laboratory studies of the newly discovered infrared band at 4705.2 cm-1 (2.1253 micrometers) in the spectrum of Io: the tentative identification of CO2

We discuss over 120 laboratory experiments pertaining to the identification of the new absorption band discovered by Trafton et al. (1991) at 4705.2 cm-1 (2.1253 micrometers) in the spectrum of Io. It is shown that this band is not due to overtones or combinations of the fundamental bands associated with the molecules (or their chemical complexes) already identified on Io, namely, SO2, H2S, and H2O. Thus, this band is due to a new, previously unidentified, component of Io. Experiments also demonstrate that the band is not due to molecular H2 frozen in SO2 frosts. Since the frequency of this band is very close to the first overtone of the nu 3 asymmetric stretching mode of CO2, we have investigated the spectral behavior of CO2 under a variety of conditions appropriate for Io. The profile of the Io band is not consistent with the rotational envelope expected for single, freely rotating, gaseous CO2 under Io-like conditions. It was found that pure, solid CO2 and CO2 intimately mixed in a matrix of solid SO2 and H2S produce bands with similar widths (5-10 cm-1), but that these bands consistently fall at frequencies about 10-20 cm-1 (approximately 0.007 micrometer) lower than the Io band. CO2 in SO2 : H2S ices also produces several additional bands that are not in the Io spectra. The spectral fit improves, however, as the CO2 concentration in SO2 increases, suggesting that CO2-CO2 interactions might be involved. A series of Ar : CO2 and Kr : CO2 matrix isolation experiments, as well as laboratory work done elsewhere, show that CO2 clustering shifts the band position to higher frequencies and provides a better fit to the Io band. Various laboratory experiments have shown that gaseous CO2 molecules have a propensity to cluster between 80 and 100 K, temperatures similar to those found on the colder regions of Io. We thus tentatively identify the newly discovered Io band at 4705.2 cm-1 (2.1253 micrometers) with CO2 multimers or "clusters" on Io. Whether these clusters are buried within an SO2 frost, reside on the surface, or are in a residual, steady-state "atmospheric aerosol" population over local coldtraps is not entirely clear, although we presently favor the latter possibility. The size of these clusters is not well defined, but evidence suggests groups of more than four molecules are required. The absorption strength of the 2 nu 3 CO2 cluster overtone determined in the laboratory, in conjunction with the observed strength of the Io band, suggests that the disk-integrated abundance of CO2 is less than 1% that of the SO2. Studies of the sublimation behavior of CO2 indicate that it probably resides predominantly in the cooler areas (< 100 K) of Io. The relative constancy of the Io feature over a variety of orbital phases suggests that the polar regions may contain much of the material. Some consequences of the physical properties of CO2 under conditions pertinent to Io are discussed. The presence of CO2 clusters on Io could be verified by the detection of any one of several other infrared bands associated with the CO2 molecule, of which the strongest are the nu 3 12CO2 asymmetric stretch fundamental near 2350 cm-1 (4.25 micrometers) and the nu 2 bending mode fundamental near 660 cm-1 (15.1 micrometers). Weaker bands that may also be detectable include the nu 3 13CO2 asymmetric stretch fundamental near 2280 cm-1 (4.39 micrometers), the 2 nu 2 + nu 3 combination/overtone band near 3600 cm-1 (2.78 micrometers), and the nu 1 + nu 3 combination band near 3705 cm-1 (2.70 micrometers).

NASA Center ARC↗
NASA NTRSDate not supplied

Matrix Theory for Data Association in PVS

Consider a collection of data generated by sensors from a set of aircraft. Data association is the process of connecting each sensor measurement with its corresponding aircraft. Furthermore once the data association has taken place, the state of the aircraft can be approximated using a Kalman filter. This talk aims to explore formal specification and verification of data association in the Prototype Verification System (PVS). Formal specification and verification of data association includes development of Kalman filters, Mahalanobis distance, and other topics of matrix analysis in PVS.

Linear Algebra↗
NASA NTRS1993

Improving the chi-squared approximation for bivariate normal tolerance regions

Let X be a two-dimensional random variable distributed according to N2(mu,Sigma) and let bar-X and S be the respective sample mean and covariance matrix calculated from N observations of X. Given a containment probability beta and a level of confidence gamma, we seek a number c, depending only on N, beta, and gamma such that the ellipsoid R = (x: (x - bar-X)'S(exp -1) (x - bar-X) less than or = c) is a tolerance region of content beta and level gamma; i.e., R has probability gamma of containing at least 100 beta percent of the distribution of X. Various approximations for c exist in the literature, but one of the simplest to compute -- a multiple of the ratio of certain chi-squared percentage points -- is badly biased for small N. For the bivariate normal case, most of the bias can be removed by simple adjustment using a factor A which depends on beta and gamma. This paper provides values of A for various beta and gamma so that the simple approximation for c can be made viable for any reasonable sample size. The methodology provides an illustrative example of how a combination of Monte-Carlo simulation and simple regression modelling can be used to improve an existing approximation.

Feiveson, Alan H.↗
NASA NTRS1976

The small motor plume/material impingement test program at MSFC utilizing nine different propellant formulations (April and May 1974), section 2

The relative effects of several candidate SSRM propellant formulations and their plume impingement effects on HRSI and RCC materials were evaluated. Nine solid propellant formulations were tested. The selected propellant matrix allowed an evaluation of propellants with and without metal additives, with and without burning rate catalyst, and low (approximately 1927 C) and high (approximately 2649 C) combustion temperatures. Motors were fired at a simulated SRB staging altitude of 3.96 km (130,000 ft) (nominal). The altitude pressure was predicted to drop approximately 0.6 km (20,000 ft) during a motor firing. All motors were loaded with 1.8 to 2.3 kg (4 to 5 lb) of propellant and burned for approximately 2 s.

Source record↗
NASA NTRS1990

Convergence of the SMI and the diagonally loaded SMI algorithms with weak interference

Approximations for the power levels at the output of an adaptive array that uses the diagonally loaded sample-matrix inversion (SMI) algorithm are derived. Diagonal loading is a technique where the diagonal of the covariance matrix is augmented with a positive or negative constant prior to inversion. This paper examines how the signal-to-interference-plus-noise ratio (SINR) and signal-to-interference ratio (SIR) at the array output vary with the number of samples taken when the input signals are CW. It is shown that positive loading produces more rapid convergence with a reduction in output SIR. Negative loading provides an improved SIR level, but positive loading produces more rapid convergence. Negative loading provides an improved SIR level, but the output power levels are erratic and slow to converge. Simulation results which verify the theoretical procedure are given.

Ganz, Matthew W.↗
NASA NTRS1997

Wave Number Selection for Incompressible Parallel Jet Flows Periodic in Space

The temporal instability of a spatially periodic parallel flow of an incompressible inviscid fluid for various jet velocity profiles is studied numerically using Floquet Analysis. The transition matrix at the end of a period is evaluated by direct numerical integration. For verification, a method based on approximating a continuous function by a series of step functions was used. Unstable solutions were found only over a limited range of wave numbers and have a band type structure. The results obtained are analogous to the behavior observed in systems exhibiting complexity at the edge of order and chaos.

Miles, Jeffrey Hilton↗