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 721 records · Page 40

Landsat-D thematic mapper simulation using aircraft multispectral scanner data

A simulation of imagery from the upcoming Landsat-D Thematic Mapper was accomplished by using selected channels of aircraft 24-channel multispectral scanner data. The purpose was to simulate Thematic Mapper 30-meter resolution imagery, to compare its spectral quality with the original aircraft MSS data, and to determine changes in thematic classification accuracy for the simulated imagery. The original resolution of approximately 7.5 meters IFOV and simulated resolution of 15, 30, and 60 meters were used to indicate the trend of spectral quality and classification accuracy. The study was based in a 6.5 square kilometer area of urban Los Angeles having a diversity of land use. The original imagery was reduced in resolution by two related methods: pixel matrix averaging, and matrix smoothing with a unity box filter, followed by matrix averaging. Thematic land use classification using training sites and a Bayesian maximum-likelihood algorithm was performed at three levels of standard deviation - 1.0, 2.0, and 3.0 sigma. Plots of relative standard deviation showed that for larger training sites with a normal distribution of data, as the resolution decreased, the distribution range of density values decreased. Also, the classification accuracies for three levels of standard deviation increased as resolution decreased. However, the indication is that a point of diminishing returns had been reached, and 30 meters IFOV should be the best for multispectral classification of urban scenes.

Clark, J.↗

Parity-violating electric-dipole transitions in helium

The paper examines parity-violating electric-dipole transitions in He in order to gain insight into the reliability of approximate calculations which are carried out for transitions in many-electron atoms. The contributions of the nearest-lying states are computed with a variety of wave functions, including very simple product wave functions, Hartree-Fock functions and Hylleraas-type wave functions with up to 84 parameters. It is found that values of the matrix elements of the parity-violating interaction can differ considerably from the values obtained from the good wave functions, even when these simple wave functions give accurate values for the matrix elements in question

Hiller, J.↗

A globally well-posed finite element algorithm for aerodynamics applications

A finite element CFD algorithm is developed for Euler and Navier-Stokes aerodynamic applications. For the linear basis, the resultant approximation is at least second-order-accurate in time and space for synergistic use of three procedures: (1) a Taylor weak statement, which provides for derivation of companion conservation law systems with embedded dispersion-error control mechanisms; (2) a stiffly stable second-order-accurate implicit Rosenbrock-Runge-Kutta temporal algorithm; and (3) a matrix tensor product factorization that permits efficient numerical linear algebra handling of the terminal large-matrix statement. Thorough analyses are presented regarding well-posed boundary conditions for inviscid and viscous flow specifications. Numerical solutions are generated and compared for critical evaluation of quasi-one- and two-dimensional Euler and Navier-Stokes benchmark test problems.

Iannelli, G. S.↗

Krylov subspace methods - Theory, algorithms, and applications

Projection methods based on Krylov subspaces for solving various types of scientific problems are reviewed. The main idea of this class of methods when applied to a linear system Ax = b, is to generate in some manner an approximate solution to the original problem from the so-called Krylov subspace span. Thus, the original problem of size N is approximated by one of dimension m, typically much smaller than N. Krylov subspace methods have been very successful in solving linear systems and eigenvalue problems and are now becoming popular for solving nonlinear equations. The main ideas in Krylov subspace methods are shown and their use in solving linear systems, eigenvalue problems, parabolic partial differential equations, Liapunov matrix equations, and nonlinear system of equations are discussed.

Sad, Youcef↗

Implicit method for the computation of unsteady flows on unstructured grids

An implicit method for the computation of unsteady flows on unstructured grids is presented. Following a finite difference approximation for the time derivative, the resulting nonlinear system of equations is solved at each time step by using an agglomeration multigrid procedure. The method allows for arbitrarily large time steps and is efficient in terms of computational effort and storage. Inviscid and viscous unsteady flows are computed to validate the procedure. The issue of the mass matrix which arises with vertex-centered finite volume schemes is addressed. The present formulation allows the mass matrix to be inverted indirectly. A mesh point movement and reconnection procedure is described that allows the grids to evolve with the motion of bodies. As an example of flow over bodies in relative motion, flow over a multi-element airfoil system undergoing deployment is computed.

Venkatakrishnan, V.↗

A T-matrix theory of galactic heavy-ion fragmentation

The theory of galactic heavy ion fragmentation is furthered by incorporating a T matrix approach into the description of the three step process of abrasion, ablation, and final state interations. The connection between this T matrix and the interaction potential is derived. For resonant states, the substitution of complex energies for real energies in the transition rate is formerly justified for up to third order processes. The previously developed abrasion-ablation fragmentation theory is rederived from first principles and is shown to result from time ordering, classical probability, and zero width resonance approximations. Improvements in the accuracy of the total fragmentation cross sections require an alternative to the latter two approximations. A Lorentz invariant differential abrasion-ablation cross section is derived which explicitly includes the previously derived abrasion total cross sections. It is demonstrated that spectral and angular distributions can be obtained from the general Lorentz invariant form.

Norbury, J. W.↗

New Constraints on the Slate Islands Impact Structure, Ontario, Canada

The Slate Islands in northern Lake Superior represent the eroded remains of a complex impact crater, originally approximately 32 km in diameter. New field studies there reveal allogenic crater fill deposits along the eastern and northern portions of the islands indicating that this 500-800 Ma impact structure is not as heavily eroded as previously thought. Near the crater center, on the western side or Patterson Island, massive blocks of target rocks, enclosed within a matrix of fine-grained polymict breccia, record the extensive deformation associated with the central uplift. Shatter cones are a common structural feature on the islands and range from less than 3 cm to over 10 m in length. Although shatter cones are powerful tools for recognizing and analyzing eroded impact craters, their origin remains poorly constrained.

Sharpton, Virgil L.↗

Cloning and regulation of rat tissue inhibitor of metalloproteinases-2 in osteoblastic cells

Rat tissue inhibitor of metalloproteinases-2 (TIMP-2) was cloned from a UMR 106-01 rat osteoblastic osteosarcoma cDNA library. The 969-bp full-length clone demonstrates 98 and 86% sequence identity to human TIMP-2 at the amino acid and nucleic acid levels, respectively. Parathyroid hormone (PTH), at 10(-8) M, stimulates an approximately twofold increase in both the 4.2- and 1.0-kb transcripts over basal levels in UMR cells after 24 h of exposure. The PTH stimulation of TIMP-2 transcripts was not affected by the inhibitor of protein synthesis, cycloheximide (10(-5) M), suggesting a primary effect of the hormone. This is in contradistinction to regulation of interstitial collagenase (matrix metalloproteinase-1) by PTH in these same cells. Nuclear run-on assays demonstrate that PTH causes an increase in TIMP-2 transcription that parallels the increase in message levels. Parathyroid hormone, in its stimulation of TIMP-2 mRNA, appears to act through a signal transduction pathway involving protein kinase A (PKA) since the increase in TIMP-2 mRNA is reproduced by treatment with the cAMP analogue, 8-bromo-cAMP (5 x 10(-3) M). The protein kinase C and calcium pathways do not appear to be involved due to the lack of effect of phorbol 12-myristate 13-acetate (2.6 x 10(-6) M) and the calcium ionophore, ionomycin (10(-7) M), on TIMP-2 transcript abundance. In this respect, regulation of TIMP-2 and collagenase in osteoblastic cells by PTH are similar. However, we conclude that since stimulation of TIMP-2 transcription is a primary event, the PKA pathway must be responsible for a direct increase in transcription of this gene.

Non-NASA Center↗

LISA Parameter Estimation using Numerical Merger Waveforms

Coalescing supermassive black holes are expected to provide the strongest sources for gravitational radiation detected by LISA. Recent advances in numerical relativity provide a detailed description of the waveforms of such signals. We present a preliminary study of LISA's sensitivity to waveform parameters using a hybrid numerical/analytic waveform describing the coalescence of two equal-mass, nonspinning black holes. The Synthetic LISA software package is used to simulate the instrument response and the Fisher information matrix method is used to estimate errors in the waveform parameters. Initial results indicate that inclusion of the merger signal can significantly improve the precision of some parameter estimates. For example, the median parameter errors for an ensemble of systems with total redshifted mass of 10(exp 6) deg M solar mass at a redshift of z is approximately 1 were found to decrease by a factor of slightly more than two when the merger was included.

Thorpe, J. I.↗

Quantitative EPMA Compositional Mapping of NWA 2995: Characterization, and Petrologic Interpretation of Mafic Clasts

We present the first fully quantitative compositional maps of lunar meteorite NWA 2995 using electron microprobe stage mapping, and compare selected clast mineralogy and chemistry. NWA 2995 is a feldspathic fragmental breccia containing numerous highland fine grained lithologies, including anorthosite, norite, olivine basalt, subophitic basalt, gabbro, KREEP-like basalt, granulitic and glassy impact melts, coarse-grained mineral fragments, Fe-Ni metal, and glassy matrix [1]. Chips of NWA 2995, representing these diverse materials, were analyzed by INAA and fused-bead electron-probe microanalysis (EPMA); comparison of analytical data suggests grouping of lunar meteorites NWA 2995, 2996, 3190, 4503, 5151, and 5152. The mean composition of NWA 2995 corresponds to a 2:1 mixture of feldspathic and mare material, with approximately 5% KREEP component [2]. Clast mineral chemistry and petrologic interpretation of paired stone NWA 2996 has been reported by Mercer et al. [3], and Gross et al. [4]. This study combines advances in quantitative EPMA compositional mapping and data analysis, as applied to selected mafic clasts in a polished section of NWA 2995, to investigate the origin of mafic lithic components and to demonstrate a procedural framework for petrologic analysis.

Carpenter, P. K.↗

Kinematics of the CS method for the treatment of molecular collisions

The problem of the quantum treatment of nonreactive collisions of two simple molecules or of an atom and a molecule is considered mathematically. Kinematical structure theorems analogous to the Wigner-Eckart theorem for scalar operators are derived for the transition operator (T-operator) of a two-molecule system, where, in addition to the usual conservation laws, certain additional conservation laws are presumed satisfied in a collision. The additional conservation laws were examined in order to develop a testable approximation scheme (CS method) for the specification of the kinematics of a system of two rigid diatomic molecules. The rand p-helicity basis states for a two-rigid rotor system are defined and unitary transformations established between these and conventional bases. It is found that only questionable criteria exist for testing the on-the-energy-shell matrix of a T(energy) operator for the presence of properties of two of the additional conservation laws.

Hahne, G. E.↗

Thermoviscoplastic analysis of fibrous periodic composites using triangular subvolumes

The nonlinear viscoplastic behavior of fibrous periodic composites is analyzed by discretizing the unit cell into triangular subvolumes. A set of these subvolumes can be configured by the analyst to construct a representation for the unit cell of a periodic composite. In each step of the loading history, the total strain increment at any point is governed by an integral equation which applies to the entire composite. A Fourier series approximation allows the incremental stresses and strains to be determined within a unit cell of the periodic lattice. The nonlinearity arising from the viscoplastic behavior of the constituent materials comprising the composite is treated as fictitious body force in the governing integral equation. Specific numerical examples showing the stress distributions in the unit cell of a fibrous tungsten/copper metal matrix composite under viscoplastic loading conditions are given. The stress distribution resulting in the unit cell when the composite material is subjected to an overall transverse stress loading history perpendicular to the fibers is found to be highly heterogeneous, and typical homogenization techniques based on treating the stress and strain distributions within the constituent phases as homogeneous result in large errors under inelastic loading conditions.

Walker, Kevin P.↗

Computational Simulation of the High Strain Rate Tensile Response of Polymer Matrix Composites

A research program is underway to develop strain rate dependent deformation and failure models for the analysis of polymer matrix composites subject to high strain rate impact loads. Under these types of loading conditions, the material response can be highly strain rate dependent and nonlinear. State variable constitutive equations based on a viscoplasticity approach have been developed to model the deformation of the polymer matrix. The constitutive equations are then combined with a mechanics of materials based micromechanics model which utilizes fiber substructuring to predict the effective mechanical and thermal response of the composite. To verify the analytical model, tensile stress-strain curves are predicted for a representative composite over strain rates ranging from around 1 x 10(exp -5)/sec to approximately 400/sec. The analytical predictions compare favorably to experimentally obtained values both qualitatively and quantitatively. Effective elastic and thermal constants are predicted for another composite, and compared to finite element results.

Goldberg, Robert K.↗

A Patch Density Recommendation based on Convergence Studies for Vehicle Panel Vibration Response resulting from Excitation by a Diffuse Acoustic Field

Producing fluid structural interaction estimates of panel vibration from an applied pressure field excitation are quite dependent on the spatial correlation of the pressure field. There is a danger of either over estimating a low frequency response or under predicting broad band panel response in the more modally dense bands if the pressure field spatial correlation is not accounted for adequately. It is a useful practice to simulate the spatial correlation of the applied pressure field over a 2d surface using a matrix of small patch area regions on a finite element model (FEM). Use of a fitted function for the spatial correlation between patch centers can result in an error if the choice of patch density is not fine enough to represent the more continuous spatial correlation function throughout the intended frequency range of interest. Several patch density assumptions to approximate the fitted spatial correlation function are first evaluated using both qualitative and quantitative illustrations. The actual response of a typical vehicle panel system FEM is then examined in a convergence study where the patch density assumptions are varied over the same model. The convergence study results illustrate the impacts possible from a poor choice of patch density on the analytical response estimate. The fitted correlation function used in this study represents a diffuse acoustic field (DAF) excitation of the panel to produce vibration response.

Smith, Andrew↗

Errors induced by the neglect of polarization in radiance calculations for Rayleigh-scattering atmospheres

Although neglecting polarization and replacing the rigorous vector radiative transfer equation by its approximate scalar counterpart has no physical background, it is a widely used simplification when the incident light is unpolarized and only the intensity of the reflected light is to be computed. We employ accurate vector and scalar multiple-scattering calculations to perform a systematic study of the errors induced by the neglect of polarization in radiance calculations for a homogeneous, plane-parallel Rayleigh-scattering atmosphere (with and without depolarization) above a Lambertian surface. Specifically, we calculate percent errors in the reflected intensity for various directions of light incidence and reflection, optical thicknesses of the atmosphere, single-scattering albedos, depolarization factors, and surface albedos. The numerical data displayed can be used to decide whether or not the scalar approximation may be employed depending on the parameters of the problem. We show that the errors decrease with increasing depolarization factor and/or increasing surface albedo. For conservative or nearly conservative scattering and small surface albedos, the errors are maximum at optical thicknesses of about 1. The calculated errors may be too large for some practical applications, and, therefore, rigorous vector calculations should be employed whenever possible. However, if approximate scalar calculations are used, we recommend to avoid geometries involving phase angles equal or close to 0 deg and 90 deg, where the errors are especially significant. We propose a theoretical explanation of the large vector/scalar differences in the case of Rayleigh scattering. According to this explanation, the differences are caused by the particular structure of the Rayleigh scattering matrix and come from lower-order (except first-order) light scattering paths involving right scattering angles and right-angle rotations of the scattering plane.

Mishchenko, M. I.↗

Materials, Manufacturing, and Test Development of a Composite Fan Blade Leading Edge Subcomponent for Improved Impact Resistance

Application of polymer matrix composite materials for jet engine fan blades is becoming attractive as an alternative to metallic blades; particularly for large engines where significant weight savings are recognized on moving to a composite structure. However, the weight benefit of the composite is offset by a reduction of aerodynamic efficiency resulting from a necessary increase in blade thickness; relative to the titanium blades. Blade dimensions are largely driven by resistance to damage on bird strike. Further development of the composite material is necessary to allow composite blade designs to approximate the dimensions of a metallic fan blade. The reduction in thickness over the state of the art composite blades is expected to translate into structural weight reduction, improved aerodynamic efficiency, and therefore reduced fuel consumption. This paper presents test article design, subcomponent blade leading edge fabrication, test method development, and initial results from ballistic impact of a gelatin projectile on the leading edge of composite fan blades. The simplified test article geometry was developed to realistically simulate a blade leading edge while decreasing fabrication complexity. Impact data is presented on baseline composite blades and toughened blades; where a considerable improvement to impact resistance was recorded.

Materials↗

Materials, Manufacturing and Test Development of a Composite Fan Blade Leading Edge Subcomponent for Improved Impact Resistance

Application of polymer matrix composite materials for jet engine fan blades is becoming attractive as an alternative to metallic blades; particularly for large engines where significant weight savings are recognized on moving to a composite structure. However, the weight benefit of the composite of is offset by a reduction of aerodynamic efficiency resulting from a necessary increase in blade thickness; relative to the titanium blades. Blade dimensions are largely driven by resistance to damage on bird strike. Further development of the composite material is necessary to allow composite blade designs to approximate the dimensions of a metallic fan blade. The reduction in thickness over the state of the art composite blades is expected to translate into structural weight reduction, improved aerodynamic efficiency, and therefore reduced fuel consumption. This paper presents test article design, subcomponent blade leading edge fabrication, test method development, and initial results from ballistic impact of a gelatin projectile on the leading edge of composite fan blades. The simplified test article geometry was developed to realistically simulate a blade leading edge while decreasing fabrication complexity. Impact data is presented on baseline composite blades and toughened blades; where a considerable improvement to impact resistance was recorded.

Ballistic Impact↗

Asymptotic (h tending to infinity) absolute stability for BDFs applied to stiff differential equations

Methods based on backward differentiation formulas (BDFs) for solving stiff differential equations require iterating to approximate the solution of the corrector equation on each step. One hope for reducing the cost of this is to make do with iteration matrices that are known to have errors and to do no more iterations than are necessary to maintain the stability of the method. This paper, following work by Klopfenstein, examines the effect of errors in the iteration matrix on the stability of the method. Application of the results to an algorithm is discussed briefly.

Krogh, F. T.↗