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,351 records · Page 75

Fuzzy Modeling and Parallel Distributed Compensation for Aircraft Flight Control from Simulated Flight Data

A method is described that combines fuzzy system identification techniques with Parallel Distributed Compensation (PDC) to develop nonlinear control methods for aircraft using minimal a priori knowledge, as part of NASA’s Learn-to-Fly initiative. A fuzzy model was generated with simulated flight data, and consisted of a weighted average of multiple linear time invariant state-space cells having parameters estimated using the equation-error approach and a least-squares estimator. A compensator was designed for each subsystem using Linear Matrix Inequalities (LMI) to guarantee closed-loop stability and performance requirements. This approach is demonstrated using simulated flight data to automatically develop a fuzzy model and design control laws for a simplified longitudinal approximation of the F-16 nonlinear flight dynamics simulation. Results include a comparison of flight data with the estimated fuzzy models and simulations that illustrate the feasibility and utility of the combined fuzzy modeling and control approach.

Weinstein, Rose↗

Rarefied solids

One important limit to creating low density materials is the objects' own weight. As a solid or colloidal matrix becomes more rarefied, gravity acts destructively to compress its suporting skeleton. We describe experimental results and propose a model which matches the low gravity behavior of rarefied or fractal solids. On parabolic airplane flights, we sought to demonstrate a key component of producing higher surface area fractals. Flight paths were selected to give a range of gravity levels: 0.01 g/g(sub 0) (low), 0.16 g(sub 0) (Lunar), 0.33 g/g(sub 0) (Martian), 1 g/g(sub 0) (Earth) and 1.8 g/g(sub 0) (high) (where g(sub 0) = 980 cm/sq s). Results using the model material of hydrophobic silica indicated that stable agglomeration of such tenuous objects can increase markedly in reduced gravity. Optical characterization revealed that fractal dimension changed directly with varying gravity. As measured by fractal dimension, effective surface area and roughness increased by 40% in low gravity. This finding supports the conclusion that relieving internal weight stresses on delicate aggregates can enhance their overall size (by two orders of magnitude) and internal surface area. We conclude that gravitational restructuring limits the overall size and void content of low-density solids. These sparse colloidal regimes may present new and technologically attractive physics, ranging from improved insulators, liquid-like tension in a 'solid' matrix, and characteristically low conductivities for sound and (8 to 14 micrometers wavelength) infrared radiation.

Noever, D. A.↗

A Technique Using Principal Component Analysis to Compare Seasonal Cycles of Earth Radiation from CERES and Model Computations

A method for quantitatively comparing the seasonal cycles of two global data sets is presented. The seasonal cycles of absorbed solar radiation (ASR) and outgoing longwave radiation (OLR) have been computed from an eight-year data set from the Clouds and Earth's Radiant Energy System (CERES) scanning radiometers and from a model data set produced by the NASA Goddard Space Flight Center Global Modeling and Assimilation Office. To compare the seasonal cycles from these two data sets, principal component (PC) analysis is used, where the PCs express the time variations and the corresponding empirical orthogonal functions (EOFs) describe the geographic variations. Ocean has a long thermal response time compared to land, so land and ocean are separated for the analysis. The root-mean square values for the seasonal cycles of ASR and OLR are extremely close for the two data sets. The first three PCs are quite close, showing that the time responses and magnitudes over the globe are very similar. The agreement between the two sets of PCs is quantified by computing the matrix of inner products of the two sets. For ASR over land, the first PCs of CERES and the model agree to better than 99.9%. The EOF maps are similar for most of the globe, but differ in a few places, and the agreement of the EOF maps is likewise quantified. Maps of differences between the annual cycles show regions of agreement and disagreement.

Smith, G. Louis↗

In-Situ Scanning Electron Microscope Experiments for Microscale Mechanical Testing and Validated Modeling of Fiber Reinforced Thermoplastics

A novel, in-situ, scanning electron microscope (SEM) mechanical testing capability for materials at the microscale which provides experimental validation to a machine learning (ML) toolset for full-field validation of physics-based micromechanics models is being developed by researchers at NASA Glenn Research Center. These are enabling technologies for the integration of multiscale digital twins for materials into system level models which will result in the improved performance, material discovery, reduced production cost and time, rapid characterization, and prognostic structural health monitoring (SHM) for materials and structures for extreme environments in support of NASA space exploration missions. In order to bridge the material structure-to-system gap for digital twins, physics-based models must be experimentally validated at multiple length scales. Seminal microscale experiments, conducted at the Air Force Research Laboratory (AFRL), were limited to transverse compression of single-layer, unidirectional thermoset polymer matrix composite (PMC) micropillar specimens [1]. The early phases of the current project followed those initial results and setup to reproduce the compression testing of PMC material on the custom-built piezoelectric actuated micromechanical testing rig built by MicroTesting Solutions LLC. In this work, samples of thermoplastic PMC material were first machined into 3 mm cubes, and then further machining and final milling was done using a Focused Ion Beam (FIB). The initial experiment was done on a pillar roughly 20 µm x 20 µm x 40 µm tall. Additional pillars were milled with final sizes ranging from 20 µm x 20 µm x 40 µm tall to 40 µm x 40 µm x 65 µm tall. A speckle pattern for in-situ full-field measurements using Digital Image Correlation (DIC) was applied with platinum, which was coated on the surface, and then the FIB was used to mill away some of the coating to produce an irregular pattern of Pt on the pillar surface. The samples were loaded into the custom testing rig and placed into the SEM and loaded under compression until failure. Images were collected in the SEM during testing. Post-processing of the images was conducted using DIC to obtain full-field displacement and strain measurements elucidating the role of the matrix as well as fiber-fiber interaction at the microscale within the composite subjected to compression loading well into the non-linear regime of the material. Moreover, the evolution of fiber-matrix debonding and matrix cracking is observed in-situ at the microscale. This data, along with images segmented with a newly developed ML toolset [2], was used to create and validate physics-based micromechanics models. An image of the failed micropillar is shown in Figure 1. The techniques developed in the initial compression experiment was tailored to the validation needs of the models and expanded to include different sized samples as well as possibly tension and fatigue.

Laura Wilson↗

Fracture behavior of unidirectional boron/aluminum composite laminates

An experimental investigation of the fracture behavior of unidirectional boron/aluminum composite laminates was conducted in order to verify the results of mathematical models. These models predict the fiber stresses and displacements and the amount of damage growth in a center-notched lamina as a function of the applied remote stress and the matrix and fiber material properties. A brittle lacquer coating was used to detect the yielding in the matrix while X-ray techniques were used to determine the number of broken fibers in the laminate. The notched strengths and the amounts of damage found in the specimens agreed well with those predicted by the models. It was shown that for thin laminates the amount of damage and the fiber displacements do not depend strongly on the number of plies for a given notch width.

Jones, F. W.↗

Pyridoxine deficiency affects biomechanical properties of chick tibial bone

The mechanical integrity of bone is dependent on the bone matrix, which is believed to account for the plastic deformation of the tissue, and the mineral, which is believed to account for the elastic deformation. The validity of this model is shown in this study based on analysis of the bones of vitamin B6-deficient and vitamin B6-replete chick bones. In this model, when B6-deficient and control animals are compared, vitamin B6 deficiency has no effect on the mineral content or composition of cortical bone as measured by ash weight (63 +/- 6 vs. 58 +/- 3); mineral to matrix ratio of the FTIR spectra (4.2 +/- 0.6 vs. 4.5 +/- 0.2), line-broadening analyses of the X-ray diffraction 002 peak (beta 002 = 0.50 +/- 0.1 vs. 0.49 +/- 0.01), or other features of the infrared spectra. In contrast, collagen was significantly more extractable from vitamin B6-deficient chick bones (20 +/- 2% of total hydroxyproline extracted vs. 10 +/- 3% p < or = 0.001). The B6-deficient bones also contained an increased amount of the reducible cross-links DHLNL, dehydro-dihydroxylysinonorleucine, (1.03 +/- 0.07 vs. 0.84 +/- 0.13 p < or = 0.001); and a nonsignificant increase in HLNL, dehydro-hydroxylysinonorleucine, (0.51 +/- 0.03 vs. 0.43 +/- 0.03, p < or = 0.10). There were no significant changes in bone length, bone diameter, or area moment of inertia. In four-point bending, no significant changes in elastic modulus, stiffness, offset yield deflection, or fracture deflection were detected. However, fracture load in the B6-deficient animals was decreased from 203 +/- 35 MPa to 151 +/- 23 MPa, p < or = 0.01, and offset yield load was decreased from 165 +/- 9 MPa to 125 +/- 14 MPa, p < or = 0.05. Since earlier histomorphometric studies had demonstrated that the B6-deficient bones were osteopenic, these data suggest that although proper cortical bone mineralization occurred, the alterations of the collagen resulted in changes to bone mechanical performance.

Non-NASA Center↗

Multi-Scale Thermo-Mechanical Modeling of Porous 3D Woven TPS Materials

This work summarizes the process to compute and analyze the thermal conductivity and mechanical properties of 3D woven TPS materials such as 3MDCP (3D Mid-Density Carbon-Phenolic). The PuMA [1,2] software, developed at NASA Ames, was used to characterize different 3MDCP samples from their constituents' data, averaging their thermal conductivity and elastic properties in the three main directions to obtain their effective orthotropic thermo-mechanical properties. This was performed in multiple steps: firstly, TPS samples were digitally reconstructed using micro computed tomography (µCT) and their constituents were segmented; then, the porous matrix phase was analyzed at the micro-scale and these results were used, along with the fibers’ constituents information, to model the tows at the meso-scale; finally, results from the constituents were homogenized and used to model the thermo-mechanical behavior of the unit cell at the macro-scale. Additionally, to gain a better understanding of the tows’ morphology and distribution of the carbon-phenolic blended fibers, microscopy images of the tows’ cross-section were segmented using deep learning techniques and analyzed with PuMA. This workflow can be applied to any TPS material to computationally obtain its thermo-mechanical properties, enabling more informed TPS design and manufacturing choices.

Multi-Scale↗

Experimental Waterflow Determination of the Dynamic Hydraulic Transfer Function for the J-2X Oxidizer Turbopump: Results and Interpretation - Part Two

Experimental results describing the hydraulic dynamic pump transfer matrix (Yp) for a cavitating J-2X oxidizer turbopump inducer+impeller tested in subscale waterflow are presented. The transfer function is required for integrated vehicle pogo stability analysis as well as optimization of local inducer pumping stability. Dynamic transfer functions across widely varying pump hydrodynamic inlet conditions are extracted from measured data in conjunction with 1D-model based corrections. Derived Dynamic transfer functions are initially interpreted relative to traditional Pogo pump equations. Water-to-liquid oxygen scaling of measured cavitation characteristics are discussed. Comparison of key dynamic transfer matrix terms derived from waterflow testing are made with those implemented in preliminary Ares Upper Stage Pogo stability modeling. Alternate cavitating pump hydraulic dynamic equations are suggested which better reflect frequency dependencies of measured transfer matrices.

Zoladz, Tom↗

Study of Anti-Vortex Baffle Effect in Suppressing Swirling Flow in LOX Tank

Experimental results describing the hydraulic dynamic pump transfer matrix (Yp) for a cavitating J-2X oxidizer turbopump inducer+impeller tested in subscale waterflow are presented. The transfer function is required for integrated vehicle pogo stability analysis as well as optimization of local inducer pumping stability. Dynamic transfer functions across widely varying pump hydrodynamic inlet conditions are extracted from measured data in conjunction with 1D-model based corrections. Derived Dynamic transfer functions are initially interpreted relative to traditional Pogo pump equations. Water-to-liquid oxygen scaling of measured cavitation characteristics are discussed. Comparison of key dynamic transfer matrix terms derived from waterflow testing are made with those implemented in preliminary Ares Upper Stage Pogo stability modeling. Alternate cavitating pump hydraulic dynamic equations are suggested which better reflect frequency dependencies of measured transfer matrices.

Yang, H. Q.↗

Algorithms for Performance, Dependability, and Performability Evaluation using Stochastic Activity Networks

Modeling tools and technologies are important for aerospace development. At the University of Illinois, we have worked on advancing the state of the art in modeling by Markov reward models in two important areas: reducing the memory necessary to numerically solve systems represented as stochastic activity networks and other stochastic Petri net extensions while still obtaining solutions in a reasonable amount of time, and finding numerically stable and memory-efficient methods to solve for the reward accumulated during a finite mission time. A long standing problem when modeling with high level formalisms such as stochastic activity networks is the so-called state space explosion, where the number of states increases exponentially with size of the high level model. Thus, the corresponding Markov model becomes prohibitively large and solution is constrained by the the size of primary memory. To reduce the memory necessary to numerically solve complex systems, we propose new methods that can tolerate such large state spaces that do not require any special structure in the model (as many other techniques do). First, we develop methods that generate row and columns of the state transition-rate-matrix on-the-fly, eliminating the need to explicitly store the matrix at all. Next, we introduce a new iterative solution method, called modified adaptive Gauss-Seidel, that exhibits locality in its use of data from the state transition-rate-matrix, permitting us to cache portions of the matrix and hence reduce the solution time. Finally, we develop a new memory and computationally efficient technique for Gauss-Seidel based solvers that avoids the need for generating rows of A in order to solve Ax = b. This is a significant performance improvement for on-the-fly methods as well as other recent solution techniques based on Kronecker operators. Taken together, these new results show that one can solve very large models without any special structure.

Deavours, Daniel D.↗

Static aeroelastic analysis of a three-dimensional generic wing

A continuation of research on the static aeroelastic analysis of a generic wing configuration is presented. Results of the study of the asymmetric oblique wing model developed by Rockwell International, in conjunction with the NASA Oblique Wing Research Aircraft Program, are reported. The capability to perform static aeroelastic analyses of an oblique wing at arbitrary skew positions is demonstrated by applying the MSC/NASTRAN static analysis scheme modified by the aerodynamic influence coefficient matrix created by the NASA Ames aerodynamic panel codes. The oblique wing is studied at two skew angles, and in particular, the capability to calculate 3-D thickness effects on the aerodynamic properties of the wing is investigated. The ability to model asymmetric wings in both subsonic and supersonic Mach numbers is shown. The aerodynamic influence coefficient matrix computed by the external programs is inserted in MSC/NASTRAN static aeroelasticity analysis run stream to compute the aeroelastic deformation and internal forces. Various aerodynamic coefficients of the oblique wing were computed for two Mach numbers, 0.7 and 1.4, and the angle of attach -5 through 15 deg.

Green, John A.↗

Topic Modeling Tool for PeTaL (Periodic Table of Life)

A topic modeling tool is constructed for the purpose of providing insights from biology to the engineer within the framework of PeTaL (Periodic Table of Life). The machine learning text mining tools–latent Dirichlet allocation (LDA) and nonnegative matrix factorization (NMF) with Kullback-Leibler (KL) divergence—are used to provide topic clusters to the user. Topic clusters are the underlying themes of a paper. For the text modeling problem, NMF-KL is the equivalent of probabilistic latent semantic analysis. Both LDA and NMF-KL are top-performing modeling tools. These tools are used to identify biological specimens relevant to the user. Various organisms solve a particular survival problem in nature differently. The topic clusters allow people without domain expertise to find these cross-topic themes in the body of documents and then branch out and examine papers whose target organisms solve the engineer’s problem. Abstracts from the Journal of Experimental Biology were used as input for the clustering tool in addition to a curated set of articles for validation. The tool is able to accept alternate input sources.

Machine learning↗

Radiative interaction of atmosphere and surface: write up with elements of code

In passive satellite remote sensing of the Earth, separation of the path radiance (atmosphere-only contribution) from the surface reflection remains a “significant challenge”. Recent literature names it among the gaps in radiative transfer (RT) topics that “require continued research in the near future”. The challenge comes from multiple reflections (bouncing) between the atmosphere and surface – radiative interaction. In this paper we use a known RT technique, the matrix-operator method (MOM), and a new modification of the monochromatic vector RT (vRT) code IPOL (Intensity and POLarization) to simulate the interaction of a plane-parallel atmosphere and a few widely used surface reflection models. Following the idea of the Green’s function method, IPOL no longer takes the surface model parameters on input. Instead, it provides the path radiance, and the atmospheric reflection and transmission matrices as output. Despite many RT codes use the MOM formalism, this output does not seem common. The surface reflection matrix is computed externally. Therefore, this paper extends the Green’s function atmospheric correction technique to the case of polarized light. Aiming clarity rather than performance, we explain in Python the structure of the surface matrices for the isotropic (Lambertian), directional unpolarized, and polarized ocean reflection models. We then combine these surface matrices and the precomputed IPOL output to get numerically accurate signal at the top of atmosphere (TOA) and test it vs. published benchmarks. Then, for each benchmark scenario we show how to get the surface from the TOA signal, i.e. perform the RT-based atmospheric correction.

radiative transfer↗

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 mismatch between the fiber and matrix leads to high tensile residual stresses at the fiber/matrix interface which could lead to premature matrix cracking during cooldown. Fiber coating 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/interface/matrix system. Guidelines, in the form of simple equations, for the selection of appropriate material properties of the fiber coating, were also derived to minimize thermal residual stresses in the matrix during fabrication.

Naik, Rajiv A.↗

Supersonic Retro-Propulsion Experimental Design for Computational Fluid Dynamics Model Validation

The development of supersonic retro-propulsion, an enabling technology for heavy payload exploration missions to Mars, is the primary focus for the present paper. A new experimental model, intended to provide computational fluid dynamics model validation data, was recently designed for the Langley Research Center Unitary Plan Wind Tunnel Test Section 2. Pre-test computations were instrumental for sizing and refining the model, over the Mach number range of 2.4 to 4.6, such that tunnel blockage and internal flow separation issues would be minimized. A 5-in diameter 70-deg sphere-cone forebody, which accommodates up to four 4:1 area ratio nozzles, followed by a 10-in long cylindrical aftbody was developed for this study based on the computational results. The model was designed to allow for a large number of surface pressure measurements on the forebody and aftbody. Supplemental data included high-speed Schlieren video and internal pressures and temperatures. The run matrix was developed to allow for the quantification of various sources of experimental uncertainty, such as random errors due to run-to-run variations and bias errors due to flow field or model misalignments. Some preliminary results and observations from the test are presented, although detailed analyses of the data and uncertainties are still on going.

Berry, Scott A.↗

The effect of weak interface on transverse properties of a ceramic matrix composite

Experimental studies conducted at NASA Lewis on SiC reaction-bonded Si3N4 composite system showed that transverse stiffness and strength were much lower than those predicted from existing analytical models based on good interfacial bonding. It was believed that weakened interfaces were responsible for the decrease in tranverse properties. To support this claim, a two-dimensional FEM analysis was performed for a transverse representative volume element. Specifically, the effect of fiber/matrix displacement compatibility at the interface was studied under both tensile and compressive transverse loadings. Interface debonding was represented using active gap elements connecting the fiber and matrix. The analyses show that the transverse tensile strength and stiffness are best predicted when a debonded interface is assumed for the composite. In fact, the measured properties can be predicted by simply replacing the fibers by voids. Thus, it is found that little or no interfacial bonding exists in the composite, and that an elastic analysis can predict the transverse stiffness and strength.

Shimansky, R. A.↗

Kinematic modeling of the Galaxy. I - The Yale Bright Star Catalogue

A code has been developed to investigate the kinematic properties of stellar populations of the Galaxy by modeling the distributions of proper motion and line-of-sight velocity. The coordinate frames of reference are discussed in detail and the general matrix transformations valid beyond the solar neighborhood are given. As a first application, the kinematic properties of the magnitude-selected sample of 3993 stars in the Bright Star Catalogue with V less than 6.0 and (B-V) greater than 0 are selected. The observed kinematic distributions are compared with a conventional model that has an exponential disk with scale heights, components of velocity dispersion ellipsoid, and asymmetric drift taken from the literature and defined as functions of spectral type. The comparisons suggest that a significant fraction of the Bright Star Catalogue belongs to a relatively young population that is not properly represented in the model.

Ratnatunga, Kavan U.↗

The Development of 2700-3000 F Environmental Barrier Coatings for SiC/SiC Ceramic Matrix Composites: Challenges and Opportunities

Environmental barrier coatings (EBCs) and SiCSiC ceramic matrix composites (CMCs) systems will play a crucial role in future turbine engines for hot-section component applications because of their ability to significantly increase engine operating temperatures, reduce engine weight and cooling requirements. The development of prime-reliant environmental barrier coatings is a key to enable the applications of the envisioned 2700-3000F EBC - CMC systems to help achieve next generation engine performance and durability goals. This paper will primarily address the performance requirements and design considerations of environmental barrier coatings for turbine engine applications. The emphasis is placed on current NASA candidate environmental barrier coating systems for SiCSiC CMCs, their performance benefits and design limitations in long-term operation and combustion environments. The efforts have been also directed to developing prime-reliant, self-healing 2700F EBC bond coat; and high stability, lower thermal conductivity, and durable EBC top coats. Major technical barriers in developing environmental barrier coating systems, the coating integrations with next generation CMCs having the improved environmental stability, cyclic durability, erosion-impact resistance, and long-term system performance will be described. The research and development opportunities for turbine engine environmental barrier coating systems by utilizing improved compositions, state-of-the-art processing methods, and simulated environment testing and durability modeling will be discussed.

CMAS Resistance↗