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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Testing of DLR C/C-SiC and C/C for HIFiRE 8 Scramjet Combustor

Ceramic Matrix Composites (CMCs) have been proposed for use as lightweight hot structures in scramjet combustors. Previous studies have calculated significant weight savings by utilizing CMCs (active and passive) versus actively cooled metallic scramjet structures. Both a carbon/carbon (C/C) and a carbon/carbon-silicon carbide (C/C-SiC) material fabricated by DLR (Stuttgart, Germany) are being considered for use in a passively cooled combustor design for Hypersonic International Flight Research Experimentation (HIFiRE) 8, a joint Australia / Air Force Research Laboratory hypersonic flight program, expected to fly at Mach 7 for approximately 30 sec, at a dynamic pressure of 55 kilopascals. Flat panels of the DLR C/C and C/C-SiC materials were installed downstream of a hydrogen-fueled, dual-mode scramjet combustor and tested for several minutes at conditions simulating flight at Mach 5 and Mach 6. Gaseous hydrogen fuel was used to fuel the scramjet combustor. The test panels were instrumented with embedded Type K and Type S thermocouples. Zirconia felt insulation was used during some of the tests to reduce heat loss from the back surface and thus increase the heated surface temperature of the C/C-SiC panel approximately 177 C (350 F). The final C/C-SiC panel was tested for three cycles totaling over 135 sec at Mach 6 enthalpy. Slightly more erosion was observed on the C/C panel than the C/C-SiC panels, but both material systems demonstrated acceptable recession performance for the HIFiRE 8 flight.

Glass, David E.↗
NASA NTRSDate not supplied

Expanding the Capability of A Legacy Combustion Flametube to Test High Temperature Engine Materials in Relevant Environments

New materials and component designs are needed to advance gas turbine engine technology and provide the performance and efficiency needs for future applications. In order to advance these materials, testing in combustion environments is a critical step prior to engine testing. In this work, we detail the design and the fabrication of a materials test sector in a flametube combustor facility. The facility simulates a combustion environment similar to that experienced by components in gas turbine engines. The flow regime is characterized by a combination of high-temperature, high-velocity, high-heat flux, and high-velocity that components experience in gas turbine engines. Exposure of components in this facility allows for the study of combined environmental effects and the impact on both coating and substrate durability. The test facility may operate across a wide range of pressures from 275-400 psig (1,896-2,758 kPa) and an air flow rate of 5 lb/s (2.27 kg/s). While combustion gas temperature is expected in excess of 3,000°F (1,649°C), 900°F (482°C) cooling air may be supplied to the backside of components or test articles. The flametube combustor was previously used to evaluate fuel injectors and combustion products, and the new test configuration will also allow for materials exposure to complex, engine-like conditions. The interior of the Test Section was additively manufactured from GRCop-84 and cryogenically fit and brazed to a stainless steel 304 housing. The use of a copper liner minimizes welds and with active cooling, is expected to provide better durability over traditional hardware using stainless steel or Inconel with a ceramic liner. The Test Section has two opposing removable windows approximately 230 mm x 80 mm that can accommodate articles up to 85 mm tall. This modular design allows for custom platforms to hold coupons, panels, or airfoil shapes to be tested with minimal re-engineering or fabrication. The bolted joint and sealing remains consistent, so any new testing only needs to work within the existing design footprint. This paper will provide an overview of the facility capabilities, design considerations, as well as thermal and structural analysis of the hardware. Future testing of ceramic matrix composite (CMC) airfoils and advanced environmental barrier coatings (EBCs) will also be discussed.

Combustion↗
NASA NTRS1984

Fine-grained millimeter-sized objects in type 3 ordinary chondrites and their relation to chondrules and matrix

Fine-grained, opaque, silicate-rich matrix material, occurring as rims on chondrules and clasts, discrete clasts, and between chondrule and mineral fragments is a primitive component of chondrites and may resemble the material from which chondrules formed. To elucidate the relationship between chondrules and matrix and to help identify chondrule precursor material, the nimeralogy and bulk compositions of 30 mm-sized, fine-grained objects in Tieschitz (h3.6), Sharps (H3.4), Allan Hills A77299 (H3.7), Yamato 74191 (L3.6), Semarkona (LL3.0), and St. Mary's Co. (LL3.3) were studied. All are composed of submicron to micron sized material and may contain larger mineral fragments. Three texturally and mineralogically distinct types of objects were identified; all are minor components of the chondrites studied ( 1% except Sharps where melt-breccia-textured objects are approximately minus 1%). Bulk compositions of network and melt-breccia-textured objects define a rather narrow trend on a plot of Si/Al vs. Mg/Al, emphasizing a possible genetic relationship. Data show that Si/Mg is higher for network and melt-breccia-textured objects than for porphyritic and barred olivine chondrules. Melting in vaccuo of network-textured objects would cause a loss of Si relative to Mg (3), suggesting that material compositionally similar to fine grained objects could have been precursors to some chondrules.

Recca, S. I.↗
NASA NTRS1993

LDEF materials overview

The flight and retrieval of the National Aeronautics and Space Administration's Long Duration Exposure Facility (LDEF) provided an opportunity for the study of the low-Earth orbit (LEO) environment and long-duration space environmental effects (SEE) on materials that is unparalleled in the history of the U.S. Space Program. The 5-year, 9-month flight of LDEF greatly enhanced the potential value of all materials on LDEF to the international SEE community, compared to that of the original 1-year flight plan. The remarkable flight attitude stability of LDEF enables specific analyses of individual and combined effects of LEO environmental parameters on identical materials on the same space vehicle. NASA recognized this potential by forming the LDEF Space Environmental Effects on Materials Special Investigation Group (MSIG) to address the greatly expanded materials and LEO space environment analysis opportunities available in the LDEF structure, experiment trays, and corollary measurements so that the combined value of all LDEF materials data to current and future space missions will be addressed and documented. An overview of the interim LDEF materials findings of the principal investigators and the Materials Special Investigation Group is provided. These revelations are based on observations of LEO environmental effects on materials made in space during LDEF retrieval and during LDEF tray deintegration at the Kennedy Space Center, and on findings of approximately 1.5 years of laboratory analyses of LDEF materials by the LDEF materials scientists. These findings were extensively reviewed and discussed at the MSIG-sponsored LDEF Materials Workshop '91. The results are presented in a format that categorizes the revelations as 'clear findings' or 'obscure preliminary findings' (and progress toward their resolution), plus resultant needs for new space materials developments and ground simulation testing/analytical modeling, in seven categories: materials/environmental parameters and data bases; LDEF contamination; thermal control coatings and protective treatments; polymers and films; polymer-matrix composites; metals, ceramics, and optical materials; and systems-related materials. The utilization of LDEF materials data for future low-Earth orbit missions is also discussed, concentrating on Space Station Freedom. In general, the LDEF data is remarkably consistent; LDEF will provide a 'benchmark' for materials design data bases for satellites in low-Earth orbit. Some materials were identified to be encouragingly resistant to LEO SEE for 5.8-years; other 'space qualified' materials displayed significant environmental degradation. General contamination levels on LDEF were low, but molecular contamination was widespread; LDEF offers an unprecedented opportunity to provide a unified perspective of unmanned LEO spacecraft contamination mechanisms. New material development requirements for long-term LEO missions were identified and current ground simulation testing methods/data for new, durable materials concepts can be validated with LDEF results. LDEF findings are already being integrated into the design of Space Station Freedom.

Stein, Bland A.↗
NASA NTRS1993

Nested Krylov methods and preserving the orthogonality

Recently the GMRESR inner-outer iteraction scheme for the solution of linear systems of equations was proposed by Van der Vorst and Vuik. Similar methods have been proposed by Axelsson and Vassilevski and Saad (FGMRES). The outer iteration is GCR, which minimizes the residual over a given set of direction vectors. The inner iteration is GMRES, which at each step computes a new direction vector by approximately solving the residual equation. However, the optimality of the approximation over the space of outer search directions is ignored in the inner GMRES iteration. This leads to suboptimal corrections to the solution in the outer iteration, as components of the outer iteration directions may reenter in the inner iteration process. Therefore we propose to preserve the orthogonality relations of GCR in the inner GMRES iteration. This gives optimal corrections; however, it involves working with a singular, non-symmetric operator. We will discuss some important properties, and we will show by experiments that, in terms of matrix vector products, this modification (almost) always leads to better convergence. However, because we do more orthogonalizations, it does not always give an improved performance in CPU-time. Furthermore, we will discuss efficient implementations as well as the truncation possibilities of the outer GCR process. The experimental results indicate that for such methods it is advantageous to preserve the orthogonality in the inner iteration. Of course we can also use iteration schemes other than GMRES as the inner method; methods with short recurrences like GICGSTAB are of interest.

Desturler, Eric↗
NASA NTRS1996

Elastic Moduli and Damping of Vibrational Modes of Aluminum/Silicon Carbide Composite Beams

Elastic and shear moduli were determined for two aluminum matrix composites containing 20 and 40 volume percent discontinuous silicon carbide, respectively, using transverse, longitudinal, and torsional vibrational modes of specimens prepared as thin beams. These moduli are consistent with those determined from stress-strain measurements. The damping factors for these modes were also determined. Thermal properties are used to show that part of the damping of transverse modes is caused by the transverse thermal currents discussed by C. Zener (thermo-elastic damping); this damping is frequency-dependent with a maximum damping factor of approximately 0.002. The remaining damping is frequency-independent, and has roughly similar values in transverse, longitudinal, and torsional modes: approximately 0.0001.

Leidecker, Henning↗
NASA NTRS2010

Development of Implicit Methods in CFD NASA Ames Research Center 1970's - 1980's

The focus here is on the early development (mid 1970's-1980's) at NASA Ames Research Center of implicit methods in Computational Fluid Dynamics (CFD). A class of implicit finite difference schemes of the Beam and Warming approximate factorization type will be addressed. The emphasis will be on the Euler equations. A review of material pertinent to the solution of the Euler equations within the framework of implicit methods will be presented. The eigensystem of the equations will be used extensively in developing a framework for various methods applied to the Euler equations. The development and analysis of various aspects of this class of schemes will be given along with the motivations behind many of the choices. Various acceleration and efficiency modifications such as matrix reduction, diagonalization and flux split schemes will be presented.

Pulliam, Thomas H.↗
NASA NTRS2016

Development of a Subcell Based Modeling Approach for Modeling the Architecturally Dependent Impact Response of Triaxially Braided Polymer Matrix Composites

Understanding the high velocity impact response of polymer matrix composites with complex architectures is critical to many aerospace applications, including engine fan blade containment systems where the structure must be able to completely contain fan blades in the event of a blade-out. Despite the benefits offered by these materials, the complex nature of textile composites presents a significant challenge for the prediction of deformation and damage under both quasi-static and impact loading conditions. The relatively large mesoscale repeating unit cell (in comparison to the size of structural components) causes the material to behave like a structure rather than a homogeneous material. Impact experiments conducted at NASA Glenn Research Center have shown the damage patterns to be a function of the underlying material architecture. Traditional computational techniques that involve modeling these materials using smeared homogeneous, orthotropic material properties at the macroscale result in simulated damage patterns that are a function of the structural geometry, but not the material architecture. In order to preserve heterogeneity at the highest length scale in a robust yet computationally efficient manner, and capture the architecturally dependent damage patterns, a previously-developed subcell modeling approach where the braided composite unit cell is approximated as a series of four adjacent laminated composites is utilized. This work discusses the implementation of the subcell methodology into the commercial transient dynamic finite element code LS-DYNA (Livermore Software Technology Corp.). Verification and validation studies are also presented, including simulation of the tensile response of straight-sided and notched quasi-static coupons composed of a T700/PR520 triaxially braided [0deg/60deg/-60deg] composite. Based on the results of the verification and validation studies, advantages and limitations of the methodology as well as plans for future work are discussed.

finite element method↗
NASA NTRS1991

Sensor-actuator placement for flexible structures with actuator dynamics

A novel approach for placement of sensors and actuators in control of flexible space structures is developed. Using an approximation of the control forces and output measurements by spatially continuous functions, the approach follows a nonlinear programming technique to determine optimal locations for sensors and actuators. Two different criteria are considered for the placement of sensors and actuators. The first criterion optimizes the location of the sensors and actuators in order to move the transmission zeros of the system farther to the left of the imaginary axis. The second criterion, however, places the sensors and actuators to optimize a function of the singular values of the Hankel matrix, which includes both measures of controllability and observability. Moreover, the effect of actuator dynamics in the placement of sensors and actuators is investigated.

Maghami, P. G.↗
NASA NTRS1972

Study of modal coupling procedures for the shuttle: A matrix method for damping synthesis

The damping method was applied successfully to real structures as well as analytical models. It depends on the ability to determine an appropriate modal damping matrix for each substructure. In the past, modal damping matrices were assumed diagonal for lack of being able to determine the coupling terms which are significant in the general case of nonproportional damping. This problem was overcome by formulating the damped equations of motion as a linear perturbation of the undamped equations for light structural damping. Damped modes are defined as complex vectors derived from the complex frequency response vectors of each substructure and are obtained directly from sinusoidal vibration tests. The damped modes are used to compute first order approximations to the modal damping matrices. The perturbation approach avoids ever having to solve a complex eigenvalue problem.

Hasselman, T. K.↗
NASA NTRS1982

Pole placement and order reduction in two-time-scale control systems through Riccati iteration

A transformation of variables taken from singular perturbations may be applied to two-time-scale linear systems in state space form to reduce the system to block-diagonal form with slow and fast modes decoupled. The transformation is easily computed by applying the new Riccati iteration. The iteration yields a solution to the nonsymmetric algebraic Riccati equation obtained by partitioning the original system matrix A. The numerical procedure is initiated with the trivial iterate L(0) = 0, and is globally convergent to the desired unique time scale decoupling solution. After transformation, the decoupled system may be used in controller design to achieve exact closed loop pole placement in the slow subsystem without altering the poles of the fast subsystem. The decoupled form may also be used to reduce system order by wetting a small parameter to zero. Provided the fast subsystem is stable, the order reduction can be expected to yield a good approximation to the original system. These methods are demonstrated using the 16th order linear model of a turbofan engine.

Anderson, L. R.↗
NASA NTRS1995

Design of a unidirectional composite momentum wheel rim

A preliminary study comparing twelve unidirectional-fiber composite systems to five metal materials conventionally used in momentum wheels is presented. Six different fibers are considered in the study: E-Glass, S-Glass, Boron, AS, T300, and Kevlar. Because of the possibility of high momentum requirements, and thus high stresses, only two matrix materials are considered: a high-modulus (HM) and a intermediate-modulus-high-strength (IMHS) matrix. Each of the six fibers are coupled with each of the two matrix materials. In an effort to optimize the composite system, each composite is considered while varying the fiber volume ratio from 0.0 to 0.7 in increments of 0.1. For fiber volume ratios above 0.2, all twelve unidirectional-fiber composite systems meet the study's requirements with higher factors of safety and less mass than the five conventional isotropic (metal) materials. For example, at a fiber volume ratio of 0.6, the Kevlar/IMHS composite system has a safety factor 4.5 times greater than that of a steel (maraging) system and an approximately 10 percent reduction in weight.

Shogrin, Bradley↗
NASA NTRS1995

CVD Fiber Coatings for Al2O3/NiAl Composites

While sapphire-fiber-reinforced nickel aluminide (Al2O3/NiAl) composites are an attractive candidate for high-temperature structures, the significant difference in the coefficient of thermal expansion between the NiAl matrix and the sapphire fiber creates substantial residual stresses in the composite. This study seeks to produce two fiber-coating systems with the potential to reduce the residual stresses in the sapphire/NiAl composite system. Chemical vapor deposition (CVD) was used to produce both the compensating and compliant-fiber coatings for use in sapphire/NiAl composites. A special reactor was designed and built to produce the FGM and to handle the toxic nickel precursors. This process was successfully used to produce 500-foot lengths of fiber with coating thicknesses of approximately 3 microns, 5 microns, and 10 microns.

Boss, Daniel E.↗
NASA NTRS2000

Far Ultraviolet Imaging from the Image Spacecraft: Wideband FUV Imaging - 2

The Far Ultraviolet Wideband Imaging Camera (WIC) complements the magnetospheric images taken by the IMAGE satellite instruments with simultaneous global maps of the terrestrial aurora. Thus, a primary requirement of WIC is to image the total intensity of the aurora in wavelength regions most representative of the aurora] source and least contaminated by dayglow, have sufficient field of view to cover the entire polar region from spacecraft apogee and have resolution that is Sufficient to resolve auroras on a scale of 1 to 2 latitude degrees, The instrument is sensitive in the spectral region from 140- 190 nm. The WIC is mounted on the rotating, IMAGE spacecraft viewing radially outward and has a field of view of 17 deg in the direction parallel to the spacecraft spin axis. Its field of view is 30 deg in the direction perpendicular to the spin axis, although only a 17 deg x 17 deg image of the Earth is recorded. The optics was an all-reflective, inverted Cassegrain Burch camera using concentric optics with a small convex primary and a large concave secondary mirror. The mirrors were coated by a special multi-layer coating, which has low reflectivity in the visible and near UV region, The detector consists of a MCP-Intensified CCD. The MCP is curved to accommodate the focal surface of the concentric optics. Tile phosphor of the image intensifier is deposited on a concave fiberoptic window, which is then Coupled to the CCD with a fiberoptic taper. The camera head operates in a fast frame transfer mode with the CCD being read approximately 30 full frames (512 by 256 pixel) per second with an exposure time of 0.033 s. The image motion (file to the satellite spin is minimal during such a short exposure. Each image is electronically distortion corrected using the look up table scheme. An offset is added to each memory address that is proportional to the image shift due to satellite rotation, and the charge signal is digitally summed in memory. On orbit, approximately 300 frames will be added to produce one WIC image in memory. The advantage of the electronic motion compensation and distortion correction is that it is extremely flexible, permitting several kinds of corrections including motions parallel and perpendicular to the predicted axis of rotation. File instrument was calibrated by applying ultraviolet light through a vacuum monochromator and measuring the absolute responsivity of the instrument. To obtain the data for the distortion look up table the camera was turned through various angles and the input angles corresponding to a pixel matrix were recorded. It was found that the spectral response peaked at 150 nm and fell off in either direction. The equivalent aperture of the camera, including mirror reflectivities and effective photocathode quantum efficiency, is about 0.04 sq cm. Thus, a 100 Rayleigh LBH aurora is expected to produce 23 equivalent counts per pixel per 10 s exposure at the peak of instrument response.

Mende, S. B.↗
NASA NTRS1992

Validation of three-dimensional incompressible spatial direct numerical simulation code: A comparison with linear stability and parabolic stability equation theories for boundary-layer transition on a flat plate

Spatially evolving instabilities in a boundary layer on a flat plate are computed by direct numerical simulation (DNS) of the incompressible Navier-Stokes equations. In a truncated physical domain, a nonstaggered mesh is used for the grid. A Chebyshev-collocation method is used normal to the wall; finite difference and compact difference methods are used in the streamwise direction; and a Fourier series is used in the spanwise direction. For time stepping, implicit Crank-Nicolson and explicit Runge-Kutta schemes are used to the time-splitting method. The influence-matrix technique is used to solve the pressure equation. At the outflow boundary, the buffer-domain technique is used to prevent convective wave reflection or upstream propagation of information from the boundary. Results of the DNS are compared with those from both linear stability theory (LST) and parabolized stability equation (PSE) theory. Computed disturbance amplitudes and phases are in very good agreement with those of LST (for small inflow disturbance amplitudes). A measure of the sensitivity of the inflow condition is demonstrated with both LST and PSE theory used to approximate inflows. Although the DNS numerics are very different than those of PSE theory, the results are in good agreement. A small discrepancy in the results that does occur is likely a result of the variation in PSE boundary condition treatment in the far field. Finally, a small-amplitude wave triad is forced at the inflow, and simulation results are compared with those of LST. Again, very good agreement is found between DNS and LST results for the 3-D simulations, the implication being that the disturbance amplitudes are sufficiently small that nonlinear interactions are negligible.

Joslin, Ronald D.↗
NASA NTRS2015

An Empirical State Error Covariance Matrix Orbit Determination Example

State estimation techniques serve effectively to provide mean state estimates. However, the state error covariance matrices provided as part of these techniques suffer from some degree of lack of confidence in their ability to adequately describe the uncertainty in the estimated states. A specific problem with the traditional form of state error covariance matrices is that they represent only a mapping of the assumed observation error characteristics into the state space. Any errors that arise from other sources (environment modeling, precision, etc.) are not directly represented in a traditional, theoretical state error covariance matrix. First, consider that an actual observation contains only measurement error and that an estimated observation contains all other errors, known and unknown. Then it follows that a measurement residual (the difference between expected and observed measurements) contains all errors for that measurement. Therefore, a direct and appropriate inclusion of the actual measurement residuals in the state error covariance matrix of the estimate will result in an empirical state error covariance matrix. This empirical state error covariance matrix will fully include all of the errors in the state estimate. The empirical error covariance matrix is determined from a literal reinterpretation of the equations involved in the weighted least squares estimation algorithm. It is a formally correct, empirical state error covariance matrix obtained through use of the average form of the weighted measurement residual variance performance index rather than the usual total weighted residual form. Based on its formulation, this matrix will contain the total uncertainty in the state estimate, regardless as to the source of the uncertainty and whether the source is anticipated or not. It is expected that the empirical error covariance matrix will give a better, statistical representation of the state error in poorly modeled systems or when sensor performance is suspect. In its most straight forward form, the technique only requires supplemental calculations to be added to existing batch estimation algorithms. In the current problem being studied a truth model making use of gravity with spherical, J2 and J4 terms plus a standard exponential type atmosphere with simple diurnal and random walk components is used. The ability of the empirical state error covariance matrix to account for errors is investigated under four scenarios during orbit estimation. These scenarios are: exact modeling under known measurement errors, exact modeling under corrupted measurement errors, inexact modeling under known measurement errors, and inexact modeling under corrupted measurement errors. For this problem a simple analog of a distributed space surveillance network is used. The sensors in this network make only range measurements and with simple normally distributed measurement errors. The sensors are assumed to have full horizon to horizon viewing at any azimuth. For definiteness, an orbit at the approximate altitude and inclination of the International Space Station is used for the study. The comparison analyses of the data involve only total vectors. No investigation of specific orbital elements is undertaken. The total vector analyses will look at the chisquare values of the error in the difference between the estimated state and the true modeled state using both the empirical and theoretical error covariance matrices for each of scenario.

Frisbee, Joseph H., Jr.↗
NASA NTRS1980

Electromagnetic scattering from two identical pseudospheres

In an attempt to gain some insight into the problem of how close particles have to be on the average before phase effects become important, the cross section and phase matrices for two identical particles as a function of the distance between their centers for both broadside and end-on illumination have been calculated. The point dipole approximation is employed to build two pseudospherical particles each consisting of 32 dipoles with an effective size parameter of 0.9283 and refractive index of 1.54-0.0i. In regard to cross sections, it is found that in certain cases several percent deviation from the converged value, which is twice the single particle cross section, could be found up to distances of separation of 25 single particle diameters. The maxima and minima in the cross section curves can be explained in terms of simple interference effects. Also the element P11 of the phase matrix showed a great deal of structure, and the position of maxima and minima could be explained by phase shift analysis.

Kattawar, G. W.↗
NASA NTRS1994

Correlated noise in the COBE DMR sky maps

The Cosmic Background Explorer Satellite Differential Radiometer (COBE DMR) sky maps contain low-level correlated noise. We obtain estimates of the amplitude and pattern of the correlated noise from three techniques: angular averages of the covariance matrix, Monte Carlo simulations of two-point correlation functions and direct analysis of the DMR maps. The results from the three methods are mutually consistent. The noise covariance matrix of a DMR sky maps is diagonal to an accuracy of better than 1%. For a given sky pixel, the dominant noise covariance occure with the ring of pixels at an angular separation of 60 deg due to the 60 deg separation of the DMR horns. The mean covariance at 60 deg is 0.45%((sup +0.18)(sub -0.14)) of the mean variance. Additionally, the variance in a given pixel is 0.7% greater than would be expected from a single beam experiment with the same noise properties. Autocorrelation functions suffer from a approximately 1.5 sigma positive bias at 60 deg while cross-correlations have no bias. Published COBE DMR results are not significantly affected by correlated noise.

Lineweaver, C. H.↗