Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix models”

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 1,495 records · Page 83

Heating distribution comparison between asymmetric and symmetric blunt cones

An experiment was performed to compare the heating distribution between symmetric and asymmetric large-angle blunt cones, with cone angles of 100, 120, and 140 deg. These hot-wall data were obtained from models made from typical thermal protection insulation for proposed aeroassisted orbital transfer vehicles. Experimental data are compared with predictions using a boundary-layer integral matrix procedure with kinetics to determine how well the heating distribution over an asymmetric cone could be approximated using axisymmetric solutions for a cone and spherical segment. In addition, a relationship between the stagnation-point heat-transfer rate and the bow-shock standoff distance for these cones is discussed. The heat-distribution data from the symmetric and asymmetric cones were very similar. Numerical results compared well with the measured wall temperatures at the stagnation point but slightly underpredicted them over the conical portion of the models.

Stewart, D. A.↗

The effect of matrix properties and fiber properties on impact failure mechanics

The low-velocity impact problem in graphite/epoxy composite sheets must be solved before large amounts of that material can be used in commercial aircraft. Many of the low-velocity impacts that affect aircraft parts occur during normal ground operations and maintenance. Service equipment and tools have masses above 1 kg, and at velocities of less than 3 m/s can impact structural parts with energies higher than composites can endure without degradation of stiffness or strength. Simple solutions were developed for large-mass, low-velocity impacts which can be modeled as quasi-static events. Static test data and impact data show that the fiber properties control the impact energy which can be absorbed before penetration. Matrix shear strength and peel resistance control the extent of delamination. Comparison of results from tough matrix and brittle matrix composites show that although tough matrices reduce the extent of delamination, they lead to more fiber damage in the contact area.

Elber, W.↗

Modeling of Melt-Infiltrated SiC/SiC Composite Properties

The elastic properties of a two-dimensional five-harness melt-infiltrated silicon carbide fiber reinforced silicon carbide matrix (MI SiC/SiC) ceramic matrix composite (CMC) were predicted using several methods. Methods used in this analysis are multiscale laminate analysis, micromechanics-based woven composite analysis, a hybrid woven composite analysis, and two- and three-dimensional finite element analyses. The elastic properties predicted are in good agreement with each other as well as with the available measured data. However, the various methods differ from each other in three key areas: (1) the fidelity provided, (2) the efforts required for input data preparation, and (3) the computational resources required. Results also indicate that efficient methods are also able to provide a reasonable estimate of local stress fields.

Mital, Subodh K.↗

Dynamic Mode Decomposition of Unsteady Pressure-Sensitive Paint Measurements for the NASA Unitary Plan Wind Tunnel Tests

This paper describes the Dynamic Mode Decomposition (DMD) of the pressures on the scale model of the Space Launch System (SLS) Block 1 cargo vehicle with the Unsteady Pressure-Sensitive Paint (uPSP) measurements, which were collected in the Ascent Transient Aerodynamics Tests with the Unitary Plan Wind Tunnel 11-by-11-foot Transonic Wind Tunnel in September 2019 at NASA Ames Research Center. The work described in this paper is a part of NASA’s development of a new state-of-the-art uPSP capability in production wind tunnels. The conventional DMD algorithm is based on the Singular Value Decomposition (SVD) of the data matrix. For the matrix of the uPSP measurements of the SLS ATAT, the number of rows is equal to the number of nodes in the grid of the scale model, and the number of columns is equal to the number of frames in the videos taken with 4 Phantom high-speed cameras. In this paper, it is verified that, for the time series with zero mean value, the DMD is equivalent to the decomposition with the Discrete Fourier Transform (DFT). Considering the uPSP is mainly used in the assessment of the unsteady, aerodynamic phenomena, the DMD of the uPSP measurements can be implemented in two steps: (1) subtract the mean value from the uPSP measurement on each of the grid nodes; (2) apply the Fast Fourier Transform (FFT) on the resulting zero-mean time series. The DMD of the uPSP measurements with FFT has two advantages: (1) the computational complexity of FFT is O(N*logN), where N is the length of the time series; (2) compared to the SVD-based DMD algorithm, the DMD with FFT can be easily implemented in parallel processing. A sample matrix of uPSP measurements, at the size of 341 grid nodes and 128 frames, is generated. Figures 1 and 2 show the eigenvalues and the ratios of the eigenvectors, respectively, of the sample matrix, without and with the mean value removed on each of the grid nodes, computed with the SVD-based DMD and the FFT. The figures demonstrate the equivalence of the SVD-based DMD and the decomposition with DFT/FFT for the time series with zero mean value. The results of DMD of the uPSP measurements of the SLS ATAT in September 2019 are presented in the paper. The DMD modes at different frequencies are shown, the aerodynamic phenomena (e.g. shockwave and vortex shedding) are demonstrated and the correlation of the DMD modes with the test configuration parameter (e.g., the Mach Number) is discussed. Figure 3 shows a software tool to visualize the DMD modes. The code to implement the algorithm described in this paper was written in C, with libraries of FFTW for FFT and MPI/OpenMP for parallel processing, and executed on the NASA Pleiades supercomputer. Funding for this research was provided by the NASA Aerosciences Evaluation and Test Capabilities Project.

Pressure-Sensitive Paint↗

A scattering model for defoliated vegetation

A scattering model for defoliated vegetation is conceived as a layer of dielectric, finite-length cylinders with specified size and orientation distributions above an irregular ground surface. The scattering phase matrix of a single cylinder is computed, then the radiative transfer technique is applied to link volume scattering from vegetation to surface scattering from the soil surface. Polarized and depolarized scattering are computed and the effects of the cylinder size and orientation distributions are illustrated. It is found that size and orientation distributions have significant effects on the backscattered signal. The model is compared with scattering from defoliated trees and agricultural crops.

Karam, M. A.↗

Integrated High Payoff Rocket Propulsion Technology (IHPRPT) SiC Recession Model

SiC stability and recession rates were modeled in hydrogen/oxygen combustion environments for the Integrated High Payoff Rocket Propulsion Technology (IHPRPT) program. The IHPRPT program is a government and industry program to improve U.S. rocket propulsion systems. Within this program SiC-based ceramic matrix composites are being considered for transpiration cooled injector faceplates or rocket engine thrust chamber liners. Material testing under conditions representative of these environments was conducted at the NASA Glenn Research Center, Cell 22. For the study described herein, SiC degradation was modeled under these Cell 22 test conditions for comparison to actual test results: molar mixture ratio, MR (O2:H2) = 6, material temperatures to 1700 C, combustion gas pressures between 0.34 and 2.10 atm, and gas velocities between 8,000 and 12,000 fps. Recession was calculated assuming rates were controlled by volatility of thermally grown silica limited by gas boundary layer transport. Assumptions for use of this model were explored, including the presence of silica on the SiC surface, laminar gas boundary layer limited volatility, and accuracy of thermochemical data for volatile Si-O-H species. Recession rates were calculated as a function of temperature. It was found that at 1700 C, the highest temperature considered, the calculated recession rates were negligible, about 200 m/h, relative to the expected lifetime of the material. Results compared favorably to testing observations. Other mechanisms contributing to SiC recession are briefly described including consumption of underlying carbon and pitting. A simple expression for liquid flow on the material surface was developed from a one-dimensional treatment of the Navier-Stokes Equation. This relationship is useful to determine under which conditions glassy coatings or thermally grown silica would flow on the material surface, removing protective layers by shear forces. The velocity of liquid flow was found to depend on the gas velocity, the viscosity of gas and liquid, as well as the thickness of the gas boundary layer and the liquid layer. Calculated flow rates of a borosilicate glass coating compared well to flow rates observed for this coating tested on a SiC panel in Cell 22.

Opila, E. J.↗

Optical Modeling Activities for the James Webb Space Telescope (JWST) Project: Determining Image Motion and Wavefront Error Over an Extended Field of View with a Segmented Optical System - II

This is part two of a series on the optical modeling activities for JWST. Starting with the linear optical model discussed in part one, we develop centroid and wavefront error sensitivities for the special case of a segmented optical system such as JWST, where the primary mirror consists of 18 individual segments. Our approach extends standard sensitivity matrix methods used for systems consisting of monolithic optics, where the image motion is approximated by averaging ray coordinates at the image and residual wavefront error is determined with global tip/tilt removed. We develop an exact formulation using the linear optical model, and extend it to cover multiple field points for performance prediction at each instrument aboard JWST. This optical model is then driven by thermal and dynamic structural perturbations in an integrated modeling environment. Results are presented.

Howard, Joseph M.↗

An improved error assessment for the GEM-T1 gravitational model

Several tests were designed to determine the correct error variances for the Goddard Earth Model (GEM)-T1 gravitational solution which was derived exclusively from satellite tracking data. The basic method employs both wholly independent and dependent subset data solutions and produces a full field coefficient estimate of the model uncertainties. The GEM-T1 errors were further analyzed using a method based upon eigenvalue-eigenvector analysis which calibrates the entire covariance matrix. Dependent satellite and independent altimetric and surface gravity data sets, as well as independent satellite deep resonance information, confirm essentially the same error assessment. These calibrations (utilizing each of the major data subsets within the solution) yield very stable calibration factors which vary by approximately 10 percent over the range of tests employed. Measurements of gravity anomalies obtained from altimetry were also used directly as observations to show that GEM-T1 is calibrated. The mathematical representation of the covariance error in the presence of unmodeled systematic error effects in the data is analyzed and an optimum weighting technique is developed for these conditions. This technique yields an internal self-calibration of the error model, a process which GEM-T1 is shown to approximate.

Lerch, F. J.↗

Simplified micromechanical equations for thermal residual stress analysis of coated fiber composites

The fabrication of metal matrix composites poses unique problems to the materials engineer. The large thermal expansion coefficient (CTE) mismatch between the fiber and matrix leads to high tensile residual stresses at the fiber/matrix (F/M) interface which could lead to premature matrix cracking during cooldown. Fiber coatings could be used to reduce thermal residual stresses. A simple closed form analysis, based on a three phase composite cylinder model, was developed to calculate thermal residual stresses in a fiber/interphase/matrix system. Parametric studies showed that the tensile thermal residual stresses at the F/M interface were very sensitive to the CTE and thickness of the interphase layer. The modulus of the layer had only a moderate effect on tensile residual stresses. For a silicon carbide titanium aluminide composite, the tangential stresses were 20 to 30 pct. larger than the axial stresses, over a wide range of interphase layer properties, indicating a tendency to form radial matrix cracks during cooldown. Guidelines for the selection of appropriate material properties of the fiber coating were also derived in order to minimize thermal residual stresses in the matrix during fabrication.

Naik, R. A.↗

Non-Boltzmann Modeling for Air Shock-Layer Radiation at Lunar-Return Conditions

This paper investigates the non-Boltzmann modeling of the radiating atomic and molecular electronic states present in lunar-return shock-layers. The Master Equation is derived for a general atom or molecule while accounting for a variety of excitation and de-excitation mechanisms. A new set of electronic-impact excitation rates is compiled for N, O, and N2+, which are the main radiating species for most lunar-return shock-layers. Based on these new rates, a novel approach of curve-fitting the non-Boltzmann populations of the radiating atomic and molecular states is developed. This new approach provides a simple and accurate method for calculating the atomic and molecular non-Boltzmann populations while avoiding the matrix inversion procedure required for the detailed solution of the Master Equation. The radiative flux values predicted by the present detailed non-Boltzmann model and the approximate curve-fitting approach are shown to agree within 5% for the Fire 1634 s case.

Johnston, Christopher O.↗

Multiscale Modeling of Thermal Conductivity of Polymer/Carbon Nanocomposites

Molecular dynamics simulation was used to estimate the interfacial thermal (Kapitza) resistance between nanoparticles and amorphous and crystalline polymer matrices. Bulk thermal conductivities of the nanocomposites were then estimated using an established effective medium approach. To study functionalization, oligomeric ethylene-vinyl alcohol copolymers were chemically bonded to a single wall carbon nanotube. The results, in a poly(ethylene-vinyl acetate) matrix, are similar to those obtained previously for grafted linear hydrocarbon chains. To study the effect of noncovalent functionalization, two types of polyethylene matrices. -- aligned (extended-chain crystalline) vs. amorphous (random coils) were modeled. Both matrices produced the same interfacial thermal resistance values. Finally, functionalization of edges and faces of plate-like graphite nanoparticles was found to be only modestly effective in reducing the interfacial thermal resistance and improving the composite thermal conductivity

Clancy, Thomas C.↗

Optimization of Residual Stresses in MMC's Using Compensating/Compliant Interfacial Layers. Part 2: OPTCOMP User's Guide

A user's guide for the computer program OPTCOMP is presented in this report. This program provides a capability to optimize the fabrication or service-induced residual stresses in uni-directional metal matrix composites subjected to combined thermo-mechanical axisymmetric loading using compensating or compliant layers at the fiber/matrix interface. The user specifies the architecture and the initial material parameters of the interfacial region, which can be either elastic or elastoplastic, and defines the design variables, together with the objective function, the associated constraints and the loading history through a user-friendly data input interface. The optimization procedure is based on an efficient solution methodology for the elastoplastic response of an arbitrarily layered multiple concentric cylinder model that is coupled to the commercial optimization package DOT. The solution methodology for the arbitrarily layered cylinder is based on the local-global stiffness matrix formulation and Mendelson's iterative technique of successive elastic solutions developed for elastoplastic boundary-value problems. The optimization algorithm employed in DOT is based on the method of feasible directions.

Pindera, Marek-Jerzy↗

Adaptive filtering with correlated state noise

An adaptive filter which uses a minimum variance criteria to estimate state noise covariance is presented. It is not necessary to assume white state noise in order to implement the filter. Simulation results are given which demonstrate that the filter tracks a satellite in the presence of modeling errors better than a conventional minimum variance filter with state noise. It is also shown that the propagated convariance matrix of the filter is an accurate indicator of the filter's performance.

Argentiero, P.↗

The containment set approach to upsets in digital systems

An approach to the fault analysis of digital systems is proposed which combines faults, hardware, and software into one overall model. The model is useful for the consideration of intermittent/transient faults, and supports a new method, based on the concept of a containment set, for realizing transient fault tolerance without massive redundancy. The model also serves as the basis for a new approach to system fault tolerance evaluation and validation using a transition matrix defined in terms of the containment set. The application of upset theory to the design of transient fault-tolerant microprocessor controllers is discussed.

Glaser, R. E.↗

Organic compounds in meteorites and their origins

The current investigation represents an extensively updated version of a review conducted by Anders et al. (1973). The investigation takes into account the literature through mid-1980. It is pointed out that Type 1 carbonaceous chondrites (C1) contain 6% of their cosmic complement of carbon, mainly in the form of organic matter. Most authors now agree that this material represents primitive prebiotic matter. The principal questions remaining are what abiotic processes formed the organic matter, and to what extent these processes took place in locales other than the solar nebula, such as interstellar clouds or meteorite parent bodes. The problem is approached in three stages. It is attempted to reconstruct the physical conditions during condensation from the clues contained in the inorganic matrix of the meteorite. The condensation behavior of carbon under these conditions is determined on the basis of thermodynamic calculations. Model experiments on the condensation of carbon are performed, and the synthesized compounds are compared with those actually found in meteorites.

Hayatsu, R.↗

Mineralogy and petrology of the Abee enstatite chondrite breccia and its dark inclusions

A model is proposed for the petrogenesis of the Abee E4 enstatite chondrite breccia, which consists of clasts, dark inclusions and matrix, and whose dark inclusions are an unusual kind of enstatite chondritic material. When the maximum metamorphic temperature of the breccia parent material was greater than 840 C, euhedral enstatite crystals in metallic Fe, Ni, and sulfide-rich areas grew into pliable metal and sulfide. Breccia parent material was impact-excavated, admixed with dark inclusions, and rapidly cooled. During this cooling, the clast and matrix material acquired thermal remanent magnetization. A subsequent ambient magnetic field imparted a uniform net magnetic orientation to the matrix and caused the magnetic orientation of the clasts to be less random. The Abee breccia was later consolidated by shock or by shallow burial and long period, low temperature metamorphism.

Rubin, A. E.↗

The NASTRAN user's manual

All information directly associated with problem solving using the NASTRAN program is presented. This structural analysis program uses the finite element approach to structural modeling wherein the distributed finite properties of a structure are represented by a finite element of structural elements which are interconnected at a finite number of grid points, to which loads are applied and for which displacements are calculated. Procedures are described for defining and loading a structural model. Functional references for every card used for structural modeling, the NASTRAN data deck and control cards, problem solution sequences (rigid formats), using the plotting capability, writing a direct matrix abstraction program, and diagnostic messages are explained. A dictionary of mnemonics, acronyms, phrases, and other commonly used NASTRAN terms is included.

Source record↗