SEARCH · Engineering Papers
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.
A viewgraph presentation on the correlation of Polymer Matrix Composites (PMC) and Metal Matrix Composites (MMC) bonded joints using three dimensional finite element analyses with materials tests is shown.
A parallel-processing computer program finds a few eigenvalues in a sparse Hermitian matrix that contains as many as 100 million diagonal elements. This program finds the eigenvalues faster, using less memory, than do other, comparable eigensolver programs. This program implements a Lanczos algorithm in the American National Standards Institute/ International Organization for Standardization (ANSI/ISO) C computing language, using the Message Passing Interface (MPI) standard to complement an eigensolver in PARPACK. [PARPACK (Parallel Arnoldi Package) is an extension, to parallel-processing computer architectures, of ARPACK (Arnoldi Package), which is a collection of Fortran 77 subroutines that solve large-scale eigenvalue problems.] The eigensolver runs on Beowulf clusters of computers at the Jet Propulsion Laboratory (JPL).
Producing fluid structural interaction estimates of panel vibration from an applied pressure field excitation are quite dependent on the spatial correlation of the pressure field. There is a danger of either over estimating a low frequency response or under predicting broad band panel response in the more modally dense bands if the pressure field spatial correlation is not accounted for adequately. It is a useful practice to simulate the spatial correlation of the applied pressure field over a 2d surface using a matrix of small patch area regions on a finite element model (FEM). Use of a fitted function for the spatial correlation between patch centers can result in an error if the choice of patch density is not fine enough to represent the more continuous spatial correlation function throughout the intended frequency range of interest. Several patch density assumptions to approximate the fitted spatial correlation function are first evaluated using both qualitative and quantitative illustrations. The actual response of a typical vehicle panel system FEM is then examined in a convergence study where the patch density assumptions are varied over the same model. The convergence study results illustrate the impacts possible from a poor choice of patch density on the analytical response estimate. The fitted correlation function used in this study represents a diffuse acoustic field (DAF) excitation of the panel to produce vibration response.
Multifunctional Boron Nitride nanotube-Boron Nitride (BN-BN) nanocomposites for energy transducers, thermal conductors, anti-penetrator/wear resistance coatings, and radiation hardened materials for harsh environments. An all boron-nitride structured BN-BN composite is synthesized. A boron nitride containing precursor is synthesized, then mixed with boron nitride nanotubes (BNNTs) to produce a composite solution which is used to make green bodies of different forms including, for example, fibers, mats, films, and plates. The green bodies are pyrolized to facilitate transformation into BN-BN composite ceramics. The pyrolysis temperature, pressure, atmosphere and time are controlled to produce a desired BN crystalline structure. The wholly BN structured materials exhibit excellent thermal stability, high thermal conductivity, piezoelectricity as well as enhanced toughness, hardness, and radiation shielding properties. By substituting with other elements into the original structure of the nanotubes and/or matrix, new nanocomposites (i.e., BCN, BCSiN ceramics) which possess excellent hardness, tailored photonic bandgap and photoluminescence, result.
Nuclear microreactors promise to open new markets for nuclear industry. This is due to their potential competitive cost in non-traditional market segments, e.g., mines, forward military basis, extraterrestrial surfaces and remote areas, and inherently safe characteristics that make them deployable where other power sources are not available or difficult to exploit. Transition metal dihydrides have been considered among the most promising candidates for moderating nuclear microreactors. In particular, yttrium hydride (YH\textsubscript{x}) has been selected for high temperature applications due to its high thermal stability and relatively high hydrogen retention at temperatures exceeding $870^\circ C$. One of the main issues associated with the use of hydrides is that, when exposed to temperature, stress, or concentration gradients, the hydrogen contained in the metallic matrix tends to redistribute and leak from the moderating elements, potentially leading to reactivity losses and power swings. The purpose of this paper is to gain a better physical understanding of the neutronic feedback associated with the hydrogen redistribution in YH\textsubscript{x}, by using Griffin and Bison. This feedback is inherently multiphysics, since the hydrogen distribution is strongly dependent on the moderator temperature spatial distribution. In particular, we want to understand the sign (+/-) of the neutronic feedback, its order of magnitude, and its underlying physical causes.
Yttrium hydride is the main candidate for moderation of high temperature nuclear microreactors . This is due to its high thermal stability and relatively high hydrogen retention at temperatures exceeding 870 C. One of the main issues associated with the use of yttrium hydride is that, when exposed to temperature, stress, or concentration gradients, the hydrogen contained in the metallic matrix tends to redistribute and leak from the moderating elements, potentially leading to reactivity losses and power swings. This paper is focused on presenting selected results from the asymptotic hydrogen redistribution analysis performed on the Simplified Microreactor Benchmark Assessment (SiMBA) problem, a microreactor core conceptual design developed at INL. The analysis is based on full core analysis, relies uniquely on Bison to perform hydrogen redistribution calculations, and discusses the physical causes of the reactivity feedback. It was found that the feedback is negative, but it is one order of magnitude lower than the one found for the empire reactor unit-cell. This is due to the lower axial temperature gradient together with the effect of the reflector.
The abundances of the moderately volatile elements (those which condense or evaporate in the temperature range 650-1350 K) in chondritic meteorites deviate from solar (or 'cosmic', as defined by CI meteorite composition) in a manner that strongly suggests that they were determined primarily by volatility, independent of chemical affinity. Wasson and collaborators have long argued that the observed correlation of relative abundance with condensation temperature reflects a systematically selective process which favored the accretion of refractory material over volatile material from a nebula cooling from a hot initial state. This view is supported by the detailed examination of chondrules and matrix in carbonaceous meteorites and the contrast between meteorite trace element abundance patterns with those produced in heating experiments. However, the idea that the moderately volatile abundances reflect global condensation has not previously been quantitatively tested or used as a constraint on nebula evolution. We have therefore constructed models of the solar nebula designed specifically for addressing meteoritic data. Calculations so far indicate that the abundance patterns of the moderately volatile elements in chondritic meteorites can be produced naturally in an evolving nebula that cools with the reduction in opacity associated with the accumulation of meteoritic bodies. Additional information is contained in the original extended abstract.
The intent of modal analysis is to develop a reliable model of a structure by working with the analytical and experimental modal properties of frequency, damping and mode shape. In addition to identifying these modal properties, it would be desirable to determine spatially which parts of the structure are modelled poorly or well. It is shown how the pattern of discrepancies in the analytical and experimental test values for the pole and the driving point zero frequencies of a structure can be linked to discrepancies in the mass or stiffness of the structural elements. The success of the procedure depends on the numerical conditioning of a modal reference matrix. Strategies to insure adequate numerical conditioning require a formulation which avoids geometric and energy storage symmetries of the structure, and ignores structural elements which contribute negligibly small potential or kinetic energy to the excited modes. Physical insight into the numerical conditioning problem is provided by a numerical example and by localization of a mass discrepancy in a real structure based on lab tests.
Stiffness matrix for structural stability analyses on beam columns, frames, and arches
Technologies pertaining to the calibration of ultra-wide-band active-electronically-scanned arrays to compensate for transmit and receive mutual coupling between array elements and/or non-ideal isolated channels are described herein. A plurality of near-field measurements are taken by way of probe or probes that may be moved among a plurality of positions aligned with respective elements in the array. For each position of the probe, each of the elements of the array is stimulated to transmit or receive a calibration signal to or from the probe, respectively. Frequency-domain transfer functions are computed from the received signals by the probe or the element for each of the array elements in each of the positions of the probe. The inverse of the matrix of transfer functions comprise frequency-domain transmit and receive correction factors that are used to modify desired array inputs/outputs such that the modified signals correct for mutual coupling between elements in the array.
Interfacial reactions, strength degradation, and fracture morphology have been characterized in B/Al(1100) and B/Al(6061) composite systems at a volume fraction of 0.55. Reaction was promoted by isothermal exposure at 350 C or 500 C for times up to 860,000 s. Tetragonal AlB12 was identified as the reaction product in the 6061 matrix and hexagonal AlB2 in the 1100 matrix. Interface development is analyzed based on a moving boundary model modified by the addition of ternary elements. While strength decreases in both systems, the time dependence of elevated-temperature exposure is different in the two matrix materials. In the 45 deg fiber orientation, failure occurs primarily by matrix shear, whereas in the 90 deg fiber orientation, fracture morphology is characterized by a mixture of fiber splitting, ductile matrix fracture, and failure at fiber-matrix interfaces.
A finite element analysis was used to quantitatively predict the effect of matrix microcracks in the 90 deg plies of graphite/epoxy composites on the coefficient of thermal expansion in the 0 deg direction, alpha (y) (perpendicular to the cracks). Results were generated for (0m/90n)s, (0/+ or 45/90)s and (0/90/ + or - 45)s graphite/epoxy laminate configurations. Analytical predictions were compared with experimental results for the two quasi-isotropic laminate configurations. Both analytical and experimental results showed that microcracks reduced the effective stiffness of the 90 deg plies, thus causing the laminates, thermal response to be more like that of a (0) laminate. The change in alpha(y) was a function of lamina material properties, microcrack density, fiber orientation, and stacking sequence. A combination of classical lamination theory and finite element analysis was used to predict the effect of microcracks in both the 90 deg and 0 deg plies. Analytical results showed that the addition of microcracks in the 0 deg plies do affect alpha(y), but to a lesser extent than those in the 90 deg plies.
The effect of the stress-free edge on the growth of local delaminations initiating from a matrix crack in (O sub 2/90 sub 4) sub s and (+/- 45.90 sub 4) sub s glass epoxy laminates is investigated using 3D finite element analysis. The presence of high interlaminar normal stresses at the intersection (corner) of the matrix crack with the stress-free edge, suggests that a mode I delamination may initiate at the corners. The strain energy release rates (G) were calculated by modeling a uniform through-width delamination and two inclined delaminations at 10.6 deg and 45 deg to the matrix crack. All components of G have high values near the free edges. The mode I component of G is high at small delamination length and becomes zero for a delamination length of one-ply thickness. The total G values near the free edge agreed well with previously derived closed form solution. The quasi-3D solutions agreed well with the 3D interior solutions.
This report documents a series of investigations to develop an approach for structural sizing of various morphing wing concepts. For the purposes of this report, a morphing wing is one whose planform can make significant shape changes in flight - increasing wing area by 50% or more from the lowest possible area, changing sweep 30 or more, and / or increasing aspect ratio by as much as 200% from the lowest possible value. These significant changes in geometry mean that the underlying load-bearing structure changes geometry. While most finite element analysis packages provide some sort of structural optimization capability, these codes are not amenable to making significant changes in the stiffness matrix to reflect the large morphing wing planform changes. The investigations presented here use a finite element code capable of aeroelastic analysis in three different optimization approaches -a "simultaneous analysis" approach, a "sequential" approach, and an "aggregate" approach.
An indentation test technique for compressively loading the ends of individual fibers to produce debonding has been applied to metal, glass, and glass-ceramic matrix composites; bond strength values at debond initiation are calculated using a finite-element model. Results are correlated with composite longitudinal and interlaminar shear behavior for carbon and Nicalon fiber-reinforced glasses and glass-ceramics including the effects of matrix modifications, processing conditions, and high-temperature oxidation embrittlement. The data indicate that significant bonding to improve off-axis and shear properties can be tolerated before the longitudinal behavior becomes brittle. Residual stress and other mechanical bonding effects are important, but improved analyses and multiaxial interfacial failure criteria are needed to adequately interpret bond strength data in terms of composite performance.
The general non-linear three-dimensional equation for acoustic potential is derived by using a perturbation technique. The linearized axisymmetric equation is then solved by using a finite element algorithm based on the Galerkin formulation for a harmonic time dependence. The solution is carried out in complex number notation for the acoustic velocity potential. Linear, isoparametric, quadrilateral elements with non-uniform distribution across the duct section are implemented. The resultant global matrix is stored in banded form and solved by using a modified Gauss elimination technique. Sound pressure levels and acoustic velocities are calculated from post element solutions. Different duct geometries are analyzed and compared with experimental results.
A research program is in progress to develop strain rate dependent deformation and failure models for the analysis of polymer matrix composites subject to high strain rate impact loads. Strain rate dependent inelastic constitutive equations have been developed to model the polymer matrix, and have been incorporated into a micromechanics approach to analyze polymer matrix composites. The Hashin failure criterion has been implemented within the micromechanics results to predict ply failure strengths. The deformation model has been implemented within LS-DYNA, a commercially available transient dynamic finite element code. The deformation response and ply failure stresses for the representative polymer matrix composite AS4/PEEK have been predicted for a variety of fiber orientations and strain rates. The predicted results compare favorably to experimentally obtained values.
Elemental abundances (C6, N7, O8, Na11, Mg12, Al3, P15, S16, Cl17, K19, Ca20, Ti22, Mn25, Fe26, and Ni28) were obtained for a set of terrestrial fossils and the rock matrix surrounding them. Principal Component Analysis extracted five factors accounting for the 92.5% of the data variance, i.e. information content, of the elemental abundance data. Hierarchical Cluster Analysis provided unsupervised sample classification distinguishing fossil from matrix samples on the basis of either raw abundances or PCA input that agreed strongly with visual classification. A stochastic, non-linear Artificial Neural Network produced a Bayesian probability of correct sample classification. The results provide a quantitative probabilistic methodology for discriminating terrestrial fossils from the surrounding rock matrix using chemical information. To demonstrate the applicability of these techniques to the assessment of meteoritic samples or in situ extraterrestrial exploration, we present preliminary data on samples of the Orgueil meteorite. In both systems an elemental signature produces target classification decisions remarkably consistent with morphological classification by a human expert using only structural (visual) information. We discuss the possibility of implementing a complexity analysis metric capable of automating certain image analysis and pattern recognition abilities of the human eye using low magnification optical microscopy images and discuss the extension of this technique across multiple scales.