SEARCH · Engineering Papers

Results for “matrix multiplication”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 523 records · Page 29

All sources NASA NTRS DOE OSTI

NASA NTRS1976

Development of Advanced Fuel Cell System (Phase 4)

A multiple-task research and development program was performed to improve the weight, life, and performance characteristics of hydrogen-oxygen alkaline fuel cells for advanced power systems. During Phase 4, the lowest stabilized degradation rate observed in all the testing completed during four phases of the program, 1 microvolt/hour, was demonstrated. This test continues after 5,000 hours of operation. The cell incorporates a PPf anode, a 90Au/10Pt cathode, a hybrid frame, and a Fybex matrix. These elements were developed under this program to extend cell life. The result demonstrated that the 80Au/20Pt cathode is as stable as a 90Au/10Pt cathode of twice the precious metal loading, was confirmed in full-scale cells. A hybrid frame two-cell plaque with dedicated flow fields and manifolds for all fluids was demonstrated to prevent the cell-to cell electrolyte transfer that limited the endurance of multicell plaques. At the conclusion of Phase 4, more than 90,900 hours of testing had been completed and twelve different cell designs had been evaluated. A technology base has been established which is ready for evaluation at the powerplant level.

Meyer, A. P.↗

NASA NTRS1985

Orbital transfer vehicle concept definition and system analysis study, 1985. Volume 5: Work breakdown structure and dictionary

To establish consistency and visibility within the Orbital Transfer Vehicle (OTV) program, a preliminary work breakdown structure (WBS) and dictionary were developed. The dictionary contains definitions of terms to be used in conjunction with the WBS so that a clear understanding of the content of the hardware, function, and cost elements may be established. The OTV WBS matrix is a two-dimensional structure which shows the interrelationship of these dimensions: the hardware elements dimension and the phase and function dimension. The dimension of time cannot be shown graphically, but must be considered. Each cost entry varies with time so that it is necessary to know these cost values by year for budget planning and approval as well as for establishing cost streams for discounting purposes in the economic analysis. While a multiple dimensional approach may at first appear complex, it actually provides benefits which outweigh any concern. This structural interrelationship provides the capability to view and analyze the OTV costs from a number of different financial and management aspects. Cost may be summed by hardware groupings, phases, or functions. The WBS may be used in a number of dimensional or single listing format applications.

Nelson, James H.↗

NASA NTRS2005

Recent Advances in Radar Polarimetry and Polarimetric SAR Interferometry

The development of Radar Polarimetry and Radar Interferometry is advancing rapidly, and these novel radar technologies are revamping Synthetic Aperture Radar Imaging decisively. In this exposition the successive advancements are sketched; beginning with the fundamental formulations and high-lighting the salient points of these diverse remote sensing techniques. Whereas with radar polarimetry the textural fine-structure, target-orientation and shape, symmetries and material constituents can be recovered with considerable improvements above that of standard amplitude-only Polarization Radar ; with radar interferometry the spatial (in depth) structure can be explored. In Polarimetric-Interferometric Synthetic Aperture Radar (POL-IN-SAR) Imaging it is possible to recover such co-registered textural plus spatial properties simultaneously. This includes the extraction of Digital Elevation Maps (DEM) from either fully Polarimetric (scattering matrix) or Interferometric (dual antenna) SAR image data takes with the additional benefit of obtaining co-registered three-dimensional POL-IN-DEM information. Extra-Wide-Band POL-IN-SAR Imaging - when applied to Repeat-Pass Image Overlay Interferometry - provides differential background validation and measurement, stress assessment, and environmental stress-change monitoring capabilities with hitherto unattained accuracy, which are essential tools for improved global biomass estimation. More recently, by applying multiple parallel repeat-pass EWB-POL-D(RP)-IN-SAR imaging along stacked (altitudinal) or displaced (horizontal) flight-lines will result in Tomographic (Multi- Interferometric) Polarimetric SAR Stereo-Imaging , including foliage and ground penetrating capabilities. It is shown that the accelerated advancement of these modern EWB-POL-D(RP)-IN-SAR imaging techniques is of direct relevance and of paramount priority to wide-area dynamic homeland security surveillance and local-to-global environmental ground-truth measurement and validation, stress assessment, and stress-change monitoring of the terrestrial and planetary covers. In addition, various closely related topics of (i) acquiring additional and protecting existing spectral windows of the Natural Electromagnetic Spectrum (NES) pertinent to Remote Sensing; (ii) mitigating against common "Radio Frequency Interference (RFI)" and intentional Directive Jamming of Airborne & Space borne POL-IN-SAR Imaging Platforms are appraised.

Boerner, Wolfgang-Martin↗

NASA NTRS2007

Improved C/SiC Ceramic Composites Made Using PIP

Improved carbon-fiber-reinforced SiC ceramic-matrix composite (C/SiC CMC) materials, suitable for fabrication of thick-section structural components, are producible by use of a combination of raw materials and processing conditions different from such combinations used in the prior art. In comparison with prior C/SiC CMC materials, these materials have more nearly uniform density, less porosity, and greater strength. The majority of raw-material/processing-condition combinations used in the prior art involve the use of chemical vapor infiltration (CVI) for densifying the matrix. In contrast, in synthesizing a material of the present type, one uses a combination of infiltration with, and pyrolysis of, a preceramic polymer [polymer infiltration followed by pyrolysis (PIP)]. PIP processing is performed in repeated, tailored cycles of infiltration followed by pyrolysis. Densification by PIP processing takes less time and costs less than does densification by CVI. When one of these improved materials was tested by exposure to a high-temperature, inert-gas environment that caused prior C/SiC CMCs to lose strength, this material did not lose strength. (Information on the temperature and exposure time was not available at the time of writing this article.) A material of the present improved type consists, more specifically, of (1) carbon fibers coated with an engineered fiber/matrix interface material and (2) a ceramic matrix, containing SiC, derived from a pre-ceramic polymer with ceramic powder additions. The enhancements of properties of these materials relative to those of prior C/SiC CMC materials are attributable largely to engineering of the fiber/ matrix interfacial material and the densification process. The synthesis of a material of this type includes processing at an elevated temperature to a low level of open porosity. The approach followed in this processing allows one to fabricate not only simple plates but also more complexly shaped parts. The carbon fiber reinforcement in a material of this type can be in any of several alternative forms, including tow, fabric, or complex preforms containing fibers oriented in multiple directions.

Easler, Timothy↗

NASA NTRS2011

Numerical Implementation of a Multiple-ISV Thermodynamically-Based Work Potential Theory for Modeling Progressive Damage and Failure in Fiber-Reinforced Laminates

A thermodynamically-based work potential theory for modeling progressive damage and failure in fiber-reinforced laminates is presented. The current, multiple-internal state variable (ISV) formulation, enhanced Schapery theory (EST), utilizes separate ISVs for modeling the effects of damage and failure. Damage is considered to be the effect of any structural changes in a material that manifest as pre-peak non-linearity in the stress versus strain response. Conversely, failure is taken to be the effect of the evolution of any mechanisms that results in post-peak strain softening. It is assumed that matrix microdamage is the dominant damage mechanism in continuous fiber-reinforced polymer matrix laminates, and its evolution is controlled with a single ISV. Three additional ISVs are introduced to account for failure due to mode I transverse cracking, mode II transverse cracking, and mode I axial failure. Typically, failure evolution (i.e., post-peak strain softening) results in pathologically mesh dependent solutions within a finite element method (FEM) setting. Therefore, consistent character element lengths are introduced into the formulation of the evolution of the three failure ISVs. Using the stationarity of the total work potential with respect to each ISV, a set of thermodynamically consistent evolution equations for the ISVs is derived. The theory is implemented into commercial FEM software. Objectivity of total energy dissipated during the failure process, with regards to refinements in the FEM mesh, is demonstrated. The model is also verified against experimental results from two laminated, T800/3900-2 panels containing a central notch and different fiber-orientation stacking sequences. Global load versus displacement, global load versus local strain gage data, and macroscopic failure paths obtained from the models are compared to the experiments.

Pineda, Evan J.↗

NASA NTRS2012

A Thermodynamically-Based Mesh Objective Work Potential Theory for Predicting Intralaminar Progressive Damage and Failure in Fiber-Reinforced Laminates

Pineda, Evan J.↗

NASA NTRS2014

In-Flight Aeroelastic Stability of the Thermal Protection System on the NASA HIAD, Part I: Linear Theory

Conical shell theory and piston theory aerodynamics are used to study the aeroelastic stability of the thermal protection system (TPS) on the NASA Hypersonic Inflatable Aerodynamic Decelerator (HIAD). Structural models of the TPS consist of single or multiple orthotropic conical shell systems resting on several circumferential linear elastic supports. The shells in each model may have pinned (simply-supported) or elastically-supported edges. The Lagrangian is formulated in terms of the generalized coordinates for all displacements and the Rayleigh-Ritz method is used to derive the equations of motion. The natural modes of vibration and aeroelastic stability boundaries are found by calculating the eigenvalues and eigenvectors of a large coefficient matrix. When the in-flight configuration of the TPS is approximated as a single shell without elastic supports, asymmetric flutter in many circumferential waves is observed. When the elastic supports are included, the shell flutters symmetrically in zero circumferential waves. Structural damping is found to be important in this case. Aeroelastic models that consider the individual TPS layers as separate shells tend to flutter asymmetrically at high dynamic pressures relative to the single shell models. Several parameter studies also examine the effects of tension, orthotropicity, and elastic support stiffness.

Goldman, Benjamin D.↗

NASA NTRS1984

Characterization of erosion of metallic materials under cavitation attack in a mineral oil

Cavitation erosion and erosion rates of eight metallic materials representing three crystal structures were studied using a 20-kHz ultrasonic magnetostrictive oscillator in viscous mineral oil. The erosion rates of the metals with an fcc matrix were 10 to 100 times higher than that of an hcp-matrix titanium alloy. The erosion rates of iron and molybdenum, with bcc matrices, were higher than that of the titanium alloy but lower than those of the fcc metals. Scanning electron microscopy indicates that the cavitation pits are initially formed at the grain boundaries and precipitates and that the pits that formed at the triple points grew faster than the others. Transcrystalline craters formed by cavitation attack over the surface of grains and roughened the surfaces by multiple slip and twinning. Surface roughness measurements show that the pits that formed over the grain boundaries deepended faster than other pits. Computer analysis revealed that a geometric expression describes the nondimensional erosion curves during the time period 0.5 t(0) t 2.5 t(0), where t(0) is the incubation period. The fcc metals had very short incubation periods; the titanium alloy had the longest incubation period.

Rao, B. C. S.↗

NASA NTRS1992

Failure detection and fault management techniques for flush airdata sensing systems

A high-angle-of-attack flush airdata sensing system was installed and flight tested on the F-18 High Alpha Research Vehicle at NASA-Dryden. This system uses a matrix of pressure orifices arranged in concentric circles on the nose of the vehicle to determine angles of attack, angles of sideslip, dynamic pressure, and static pressure as well as other airdata parameters. Results presented use an arrangement of 11 symmetrically distributed ports on the aircraft nose. Experience with this sensing system data indicates that the primary concern for real-time implementation is the detection and management of overall system and individual pressure sensor failures. The multiple port sensing system is more tolerant to small disturbances in the measured pressure data than conventional probe-based intrusive airdata systems. However, under adverse circumstances, large undetected failures in individual pressure ports can result in algorithm divergence and catastrophic failure of the entire system. How system and individual port failures may be detected using chi sq. analysis is shown. Once identified, the effects of failures are eliminated using weighted least squares.

Whitmore, Stephen A.↗

NASA NTRS1997

Complexity of Kronecker Operations on Sparse Matrices with Applications to the Solution of Markov Models

We present a systematic discussion of algorithms to multiply a vector by a matrix expressed as the Kronecker product of sparse matrices, extending previous work in a unified notational framework. Then, we use our results to define new algorithms for the solution of large structured Markov models. In addition to a comprehensive overview of existing approaches, we give new results with respect to: (1) managing certain types of state-dependent behavior without incurring extra cost; (2) supporting both Jacobi-style and Gauss-Seidel-style methods by appropriate multiplication algorithms; (3) speeding up algorithms that consider probability vectors of size equal to the "actual" state space instead of the "potential" state space.

Buchholz, Peter↗

NASA NTRS1991

Unified aeroacoustics analysis for high speed turboprop aerodynamics and noise. Volume 4: Computer user's manual for UAAP turboprop aeroacoustic code

The Unified AeroAcoustic Program (UAAP) code calculates the airloads on a single rotation prop-fan, or propeller, and couples these airloads with an acoustic radiation theory, to provide estimates of near-field or far-field noise levels. The steady airloads can also be used to calculate the nonuniform velocity components in the propeller wake. The airloads are calculated using a three dimensional compressible panel method which considers the effects of thin, cambered, multiple blades which may be highly swept. These airloads may be either steady or unsteady. The acoustic model uses the blade thickness distribution and the steady or unsteady aerodynamic loads to calculate the acoustic radiation. The users manual for the UAAP code is divided into five sections: general code description; input description; output description; system description; and error codes. The user must have access to IMSL10 libraries (MATH and SFUN) for numerous calls made for Bessel functions and matrix inversion. For plotted output users must modify the dummy calls to plotting routines included in the code to system-specific calls appropriate to the user's installation.

Menthe, R. W.↗

NASA NTRS2010

Probabilistic Round Trip Contamination Analysis of a Mars Sample Acquisition and Handling Process Using Markovian Decompositions

A method for evaluating the probability of a Viable Earth Microorganism (VEM) contaminating a sample during the sample acquisition and handling (SAH) process of a potential future Mars Sample Return mission is developed. A scenario where multiple core samples would be acquired using a rotary percussive coring tool, deployed from an arm on a MER class rover is analyzed. The analysis is conducted in a structured way by decomposing sample acquisition and handling process into a series of discrete time steps, and breaking the physical system into a set of relevant components. At each discrete time step, two key functions are defined: The probability of a VEM being released from each component, and the transport matrix, which represents the probability of VEM transport from one component to another. By defining the expected the number of VEMs on each component at the start of the sampling process, these decompositions allow the expected number of VEMs on each component at each sampling step to be represented as a Markov chain. This formalism provides a rigorous mathematical framework in which to analyze the probability of a VEM entering the sample chain, as well as making the analysis tractable by breaking the process down into small analyzable steps.

Hudson, Nicolas↗

NASA NTRS1993

Scene segmentation of natural images using texture measures and back-propagation

Knowledge of the three-dimensional world is essential for many guidance and navigation applications. A sequence of images from an electro-optical sensor can be processed using optical flow algorithms to provide a sparse set of ranges as a function of azimuth and elevation. A natural way to enhance the range map is by interpolation. However, this should be undertaken with care since interpolation assumes continuity of range. The range is continuous in certain parts of the image and can jump at object boundaries. In such situations, the ability to detect homogeneous object regions by scene segmentation can be used to determine regions in the range map that can be enhanced by interpolation. The use of scalar features derived from the spatial gray-level dependence matrix for texture segmentation is explored. Thresholding of histograms of scalar texture features is done for several images to select scalar features which result in a meaningful segmentation of the images. Next, the selected scalar features are used with a neural net to automate the segmentation procedure. Back-propagation is used to train the feed forward neural network. The generalization of the network approach to subsequent images in the sequence is examined. It is shown that the use of multiple scalar features as input to the neural network result in a superior segmentation when compared with a single scalar feature. It is also shown that the scalar features, which are not useful individually, result in a good segmentation when used together. The methodology is applied to both indoor and outdoor images.

Sridhar, Banavar↗

NASA NTRS1995

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗

NASA NTRS1990

Unsteady transonic Navier-Stokes computations for an oscillating wing using single and multiple zones

Modern jet transports and maneuvering tactical fighters operating in the transonic regime often give rise to time-dependent fluid physics that interact with flexible structural components, e.g., vortical flow, shocks, and separation. Efficient computational fluid dynamic (CFD) methods are required to study such computationally intensive problems. In this work a numerical method is presented to address this problem. Time-dependent, compressible, Navier-Stokes equations are used to simulate unsteady transonic flow about a three-dimensional rigid wing undergoing a forced periodic motion in angle of attack. An efficient, implicit, diagonal algorithm is utilized because of its low operation count per time step compared to other methods that solve systems of block matrix equations. The formal time accuracy is addressed theoretically and demonstrated numerically by comparison of computational results with experimental data. A zonal grid approach, capable of treating complex geometries, is presented and its time accuracy is demonstrated by comparing a zonal computation with a single grid computation and experimental data.

Chaderjian, Neal M.↗

NASA NTRS1985

Characterization of erosion of metallic materials under cavitation attack in a mineral oil

Cavitation erosion and erosion rates of eight metallic materials representing three crystal structures were studied. The erosion experiments were conducted with a 20-kHz ultrasonic magnetostrictive oscillator in a viscous mineral oil. The erosion rates of the metals with an fcc matrix were 10 to 100 times higher than that of an hop-matrix titanium alloy. The erosion rates of iron and molybdenum, with bcc matrices, were higher than that of the titanium alloy but lower than those of those of the fcc materials. Studies with scanning electron microscopy indicated that the cavitation pits were initially formed at the grain boundaries and precipitates and that the pits formed at the junction of grain boundaries grew faster than the others. Transcrystalline craters formed by cavitation attack over the surface of grains and roughened the surfaces by multiple slip and twinning. Surface roughness measurements showed that the pits that formed over the grain boundaries deepened faster than pits. Computer analysis revealed that a geometric expression describes the nondimensional erosion curves during the time period 0.5 t (sub 0) t 2.5 t (sub 0), where t (sub 0) is the incubation period. The fcc metals had very short incubation periods; the titanium alloy had the longest incubation period.

Rao, B. C. S.↗

NASA NTRS2006

Exploring the Connection Between Sampling Problems in Bayesian Inference and Statistical Mechanics

The Bayesian and statistical mechanical communities often share the same objective in their work - estimating and integrating probability distribution functions (pdfs) describing stochastic systems, models or processes. Frequently, these pdfs are complex functions of random variables exhibiting multiple, well separated local minima. Conventional strategies for sampling such pdfs are inefficient, sometimes leading to an apparent non-ergodic behavior. Several recently developed techniques for handling this problem have been successfully applied in statistical mechanics. In the multicanonical and Wang-Landau Monte Carlo (MC) methods, the correct pdfs are recovered from uniform sampling of the parameter space by iteratively establishing proper weighting factors connecting these distributions. Trivial generalizations allow for sampling from any chosen pdf. The closely related transition matrix method relies on estimating transition probabilities between different states. All these methods proved to generate estimates of pdfs with high statistical accuracy. In another MC technique, parallel tempering, several random walks, each corresponding to a different value of a parameter (e.g. "temperature"), are generated and occasionally exchanged using the Metropolis criterion. This method can be considered as a statistically correct version of simulated annealing. An alternative approach is to represent the set of independent variables as a Hamiltonian system. Considerab!e progress has been made in understanding how to ensure that the system obeys the equipartition theorem or, equivalently, that coupling between the variables is correctly described. Then a host of techniques developed for dynamical systems can be used. Among them, probably the most powerful is the Adaptive Biasing Force method, in which thermodynamic integration and biased sampling are combined to yield very efficient estimates of pdfs. The third class of methods deals with transitions between states described by rate constants. These problems are isomorphic with chemical kinetics problems. Recently, several efficient techniques for this purpose have been developed based on the approach originally proposed by Gillespie. Although the utility of the techniques mentioned above for Bayesian problems has not been determined, further research along these lines is warranted

Pohorille, Andrew↗

NASA NTRS2016

Tribological Studies of Dynamic Thermal Seal Materials

Thermal seals are required on high-speed vehicles in many dynamic applications such as variable inlets in propulsion systems and control surfaces. These seals, often referred to as dynamic thermal seals, must not only mitigate inboard heat transfer, but must also exhibit sufficient durability when scrubbed against mating surfaces. For high-temperature high-speed vehicle applications, the mating surfaces are often made from thermal protection system (TPS) materials, which are typically rougher and more abrasive than TPS materials used at lower temperatures. The high-temperature TPS materials used can include non-ablative (e.g., lightweight porous oxides, ceramic matrix composites) andor ablative systems (e.g., phenolic systems). Due to the increased need for durable high-temperature dynamic seals, researchers working with the NASA Glenn Research Center embarked on an effort to (a) characterize the tribological performance of state-of-the-art thermal seal materials against a variety of TPS materials and (b) develop approaches for improved wear resistance. Tests were conducted using a recently upgraded high-temperature tribometer to assess wear resistance for a variety of tribopairs under multiple conditions. This data will begin to frame the challenges of using these materials and eventually permit an improved ability to design and implement these critical TPS components.

thermal protection↗