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.
Tolerance levels of interference from oxygen, iron, aluminum, silicon, and magnesium in neutron activation of simulated moon matrix
During the earliest period of solar system formation, water played an important role in the evolution of primitive dust, both after accretion of planetesimals and possible before accretion within the protoplanetary disk. Many chondrites show evidence of variable degrees of aqueous alteration, the CM2 chondrites being among the most studied [1]. This group of chondrites is characterized by mineral assemblages of both primary and secondary alteration phases. Hence, these meteorites retain a particularly important record of the reactions that occurred between primary high temperature nebular phases and water. Studies of these chondrites can provide information on the conditions and environments of aqueous alteration and the mobility of elements during alteration. This latter question is at the core of a debate concerning the location of aqueous alteration, i.e. whether alteration occurred predominantly within a closed system after accretion (parent body alteration) or whether some degree of alteration occurred within the solar nebula or on ephemeral protoplanetary bodies prior to accretion. At the core of the parent body alteration model is the hypothesis that elemental exchange between different components, principally chondrules and matrix, must have occurred. chondrules and matrix, must have occurred. In this study, we focus on the behavior of the minor element, phosphorus. This study was stimulated by observations of the behavior of P during the earliest stages of alteration in glassy mesostasis in type II chondrules in CR chondrites and extends the preliminary observations of on Y791198 to other CM chondrites.
The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased use in the aerospace and automotive communities. The aerospace community has identified several key capabilities which are currently lacking in the available material models in commercial transient dynamic finite element codes. To attempt to improve the predictive capability of composite impact simulations, a next generation material model is being developed for incorporation within the commercial transient dynamic finite element code LS-DYNA. The material model, which incorporates plasticity, damage and failure, utilizes experimentally based tabulated input to define the evolution of plasticity and damage and the initiation of failure as opposed to specifying discrete input parameters such as modulus and strength. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. For the damage model, a strain equivalent formulation is used to allow for the uncoupling of the deformation and damage analyses. In the damage model, a semi-coupled approach is employed where the overall damage in a particular coordinate direction is assumed to be a multiplicative combination of the damage in that direction resulting from the applied loads in various coordinate directions. For the failure model, a tabulated approach is utilized in which a stress or strain based invariant is defined as a function of the location of the current stress state in stress space to define the initiation of failure. Failure surfaces can be defined with any arbitrary shape, unlike traditional failure models where the mathematical functions used to define the failure surface impose a specific shape on the failure surface. In the current paper, the complete development of the failure model is described and the generation of a tabulated failure surface for a representative composite material is discussed.
Enhanced finite elements are elements with an embedded analytical solution which can capture detailed local fields, enabling more efficient, mesh-independent finite element analysis. In the present study, an enhanced finite element is applied to generate a general framework capable of modeling an array of joint types. The joint field equations are derived using the principle of minimum potential energy, and the resulting solutions for the displacement fields are used to generate shape functions and a stiffness matrix for a single joint finite element. This single finite element thus captures the detailed stress and strain fields within the bonded joint, but it can function within a broader structural finite element model. The costs associated with a fine mesh of the joint can thus be avoided while still obtaining a detailed solution for the joint. Additionally, the capability to model non-linear adhesive constitutive behavior has been included within the method, and progressive failure of the adhesive can be modeled by using a strain-based failure criteria and re-sizing the joint as the adhesive fails. Results of the model compare favorably with experimental and finite element results.
The cores and reflectors in high-temperature gas-cooled reactors (HTGRs) are made of graphite materials, with the graphite acting as a moderator, a fuel host matrix, or the foundation for various structural components. This study aims to survey the models in the literature for graphite materials being used as host matrices in pebble/fuel compacts and to implement those surveyed models into Bison to conduct an early assessment of graphite's thermo-mechanical response under various reactor conditions. In this study, thermal (e.g., thermal conductivity, and specific heat capacity) and mechanical (e.g., elastic properties, thermal expansion, irradiation-induced dimensional changes, and irradiation-induced creep) material models for various graphite grades (e.g., H-451, IG-110, G-348, 2020, A3-3, and A3-27) are incorporated into Bison. Two benchmark problems are then exercised utilizing these new graphite-related capabilities: (1) modeling an Advanced Gas Reactor (AGR)-2 fuel compact, and (2) modeling the debonding of a particle-matrix interface.
This report reviews the application of telescoping mechanics to composites using recursive laminate theory. The elemental scale is the fiber-matrix slice, the behavior of which propagates to laminate. The results from using applications for typical, hybrid, and smart composites and composite-enhanced reinforced concrete structures illustrate the versatility and generality of telescoping scale mechanics. Comparisons with approximate, single-cell, and two- and three-dimensional finite-element methods demonstrate the accuracy and computational effectiveness of telescoping scale mechanics for predicting complex composite behavior.
First, the Cholesky factorization is extended to cover uniformly partitioned banded positive definite matrices of rank n which may be real symmetric or Hermitian. Then, two stratagems are given for the use of the algorithm in concurrent machines where the number of processing elements is less than required to factor the matrix in as few serial steps as possible, and where uniformly high efficiency is expected from all processing elements. Expressions are given for the efficiency factor e appearing in the speed-up expression q = eN, and these are specialized for the N node hypercube machine as a function of partition size s, the number N of processing elements of the hypercube machine, and the cost mu of interelement transmission relative to computation. It is shown that the efficiency factor e is inversely proportional to mu/s, and that e is almost independent of N when N is large and mu/s = 0. The task is completed in n/s serial steps with no limit on n. The half bandwidth b of the matrix is 2 Ns.
Fiber/matrix fracture and fiber-matrix interface debonding in a metal matrix composite (MMC) are computationally simulated. These simulations are part of a research activity to develop computational methods for microfracture, microfracture propagation and fracture toughness of the metal matrix composites. The three-dimensional finite element model used in the simulation consists of a group of nine unidirectional fibers in three by three unit cell array of SiC/Ti15 metal matrix composite with a fiber volume ration of 0.35. This computational procedure is used to predict the fracture process and establish the hierarchy of fracture modes based on strain energy release rate. It is also used to predict stress redistribution to surrounding matrix-fibers due to initial and progressive fracture of fiber/matrix and due to debonding of fiber-matrix interface. Microfracture results for various loading cases such as longitudinal, transverse, shear and bending are presented and discussed. Step-by-step procedures are outlined to evaluate composite microfracture for a given composite system.
Fiber/matrix fracture and fiber-matrix interface debonding in a metal matrix composite (MMC) are computationally simulated. These simulations are part of a research activity to develop computational methods for microfracture, microfracture propagation and fracture toughness of the metal matrix composites. The three-dimensional finite element model used in the simulation consists of a group of nine unidirectional fibers in three by three unit cell array of SiC/Ti15 metal matrix composite with a fiber volume ration of 0.35. This computational procedure is used to predict the fracture process and establish the hierarchy of fracture modes based on strain energy release rate. It is also used to predict stress redistribution to surrounding matrix-fibers due to initial and progressive fracture of fiber/matrix and due to debonding of fiber-matrix interface. Microfracture results for various loading cases such as longitudinal, transverse, shear and bending are presented and discussed. Step-by-step procedures are outlined to evaluate composite microfracture for a given composite system.
This paper presents techniques based on a three-dimensional finite-element analysis for the analysis of continuous fiber reinforced metal matrix composite. Examples are shown for specific metal matrix composites such as boron/aluminum and silicon-carbide/aluminum. Specimen stress-strain behavior and stress at first fiber failure were predicted for boron/aluminum laminates containing circular holes and crack-like slits. The predictions compared very well with test data for laminates containing 0 deg fibers and reasonably well for (+/-45)2s laminates. Mesh configuration was shown to have an effect on the calculation of stresses local to the notch. The presence of thin interface layers of matrix material had a significant influence on the slit tip stress state, causing sharper stress gradients near the notch. Interface layers reduced the slit-tip fibers stress in a (+/-45)2s silicon-carbide/aluminum laminate but increased them in a (0/90)2s laminate.
A fluid finite element has been developed for a quasi-compressible fluid. Both kinetic and potential energy are expressed as functions of nodal displacements. Thus, the formulation is similar to that used for structural elements, with the only differences being that the fluid can possess gravitational potential, and the constitutive equations for fluid contain no shear coefficients. Using this approach, structural and fluid elements can be used interchangeably in existing efficient sparse-matrix structural computer programs such as SPAR. The theoretical development of the element formulations and the relationships of the local and global coordinates are shown. Solutions of fluid slosh, liquid compressibility, and coupled fluid-shell oscillation problems which were completed using a temporary digital computer program are shown. The frequency correlation of the solutions with classical theory is excellent.
Last year we presented techniques for the detection of fossils during robotic missions to Mars using both structural and chemical signatures[Storrie-Lombardi and Hoover, 2004]. Analyses included lossless compression of photographic images to estimate the relative complexity of a putative fossil compared to the rock matrix [Corsetti and Storrie-Lombardi, 2003] and elemental abundance distributions to provide mineralogical classification of the rock matrix [Storrie-Lombardi and Fisk, 2004]. We presented a classification strategy employing two exploratory classification algorithms (Principal Component Analysis and Hierarchical Cluster Analysis) and non-linear stochastic neural network to produce a Bayesian estimate of classification accuracy. We now present an extension of our previous experiments exploring putative fossil forms morphologically resembling cyanobacteria discovered in the Orgueil meteorite. Elemental abundances (C6, N7, O8, Na11, Mg12, Ai13, Si14, P15, S16, Cl17, K19, Ca20, Fe26) obtained for both extant cyanobacteria and fossil trilobites produce signatures readily distinguishing them from meteorite targets. When compared to elemental abundance signatures for extant cyanobacteria Orgueil structures exhibit decreased abundances for C6, N7, Na11, All3, P15, Cl17, K19, Ca20 and increases in Mg12, S16, Fe26. Diatoms and silicified portions of cyanobacterial sheaths exhibiting high levels of silicon and correspondingly low levels of carbon cluster more closely with terrestrial fossils than with extant cyanobacteria. Compression indices verify that variations in random and redundant textural patterns between perceived forms and the background matrix contribute significantly to morphological visual identification. The results provide a quantitative probabilistic methodology for discriminating putatitive fossils from the surrounding rock matrix and &om extant organisms using both structural and chemical information. The techniques described appear applicable to the geobiological analysis of meteoritic samples or in situ exploration of the Mars regolith. Keywords: cyanobacteria, microfossils, Mars, elemental abundances, complexity analysis, multifactor analysis, principal component analysis, hierarchical cluster analysis, artificial neural networks, paleo-biosignatures
The bridging of matrix cracks by fibers is an important toughening mechanism in fiber reinforced brittle matrix composites. This paper presents the results of a nonlinear finite element analysis of the Mode-I propagation of a bridged matrix crack in a finite size specimen. The composite is modeled as an orthotropic continuum and the bridging due to the fibers is modeled as a distribution of tractions which resist crack opening. A critical stress intensity factor criterion is employed for matrix crack propagation while a critical crack opening condition is used for fiber failure. The structural response of the specimen (load-deflection curves) as well as the stress intensity factor of the propagating crack are calculated for various constituent properties and specimen configurations for both tensile and bending loading. By controlling the length of the bridged crack results are obtained which highlight the transition from stable to unstable behavior of the propagating crack.
The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites is becoming critical as these materials are gaining increased usage in the aerospace and automotive industries. While there are several composite material models currently available within LS-DYNA (Registered), there are several features that have been identified that could improve the predictive capability of a composite model. To address these needs, a combined plasticity and damage model suitable for use with both solid and shell elements is being developed and is being implemented into LS-DYNA as MAT_213. A key feature of the improved material model is the use of tabulated stress-strain data in a variety of coordinate directions to fully define the stress-strain response of the material. To date, the model development efforts have focused on creating the plasticity portion of the model. The Tsai-Wu composite failure model has been generalized and extended to a strain-hardening based orthotropic material model with a non-associative flow rule. The coefficients of the yield function, and the stresses to be used in both the yield function and the flow rule, are computed based on the input stress-strain curves using the effective plastic strain as the tracking variable. The coefficients in the flow rule are computed based on the obtained stress-strain data. The developed material model is suitable for implementation within LS-DYNA for use in analyzing the nonlinear response of polymer composites.
On September 24, 2023, NASA’s OSIRIS-REx mission returned a capsule to Earth carrying material from asteroid Bennu. This event was the first time a U.S. mission delivered pristine samples from an asteroid and is the largest asteroid sample return to date. As these samples represent some of the oldest, most primitive, and pristine materials available to us, and which originate from a known and well-studied asteroid, they allow us a rare opportunity to gain a better understanding of the formation and evolution of our solar system. Key to understanding the material returned from asteroid Bennu is establishing its bulk chemical composition. Previous studies have shown that each chondrite group has a distinct elemental composition. For the carbonaceous chondrites specifically, each group exhibits a distinct pattern of moderately and highly volatile elemental depletions, relative to CI chondrites. CI chondrites are considered the most primitive chondrite group and broadly represent the solar photosphere composition. The two most striking features of these depletion patterns is that the moderately volatile element depletions increase with decreasing 50% condensation temperature and then plateau out at abundances that roughly correlate with matrix abundance. Due to these distinctive patterns, bulk elemental composition has become an important classification tool for establishing the different chondrite groups and the connections between them. As such, determining the bulk chemical composition of the Bennu aggregates will help to test two mission hypothesis: “Bennu’s bulk elemental composition reflects that of its main parent asteroid and is similar to the composition of the Sun, with depletions in moderately to highly volatile elements” and “Bennu’s dominant lithologies are comparable in bulk mineralogy, petrology, and composition to the most aqueously altered carbonaceous chondrites” . Furthermore as the carbonaceous and non-carbonaceous chondrites are thought to have likely formed in the inner and outer protoplanetary disk, respectively (e.g., [8]), determining the bulk chemical compositions of the Bennu aggregates will also help test the major mission hypothesis that “Bennu's parent body formed beyond the snow line by accretion of material in the pro-toplanetary disk” . Regarding chondrite formation, observed elemental patterns have been successfully used to model how the mixing of volatile-rich and volatile-poor chondritic components can produce the observed carbonaceous chondrites groups. They have also been significant in investigating how volatilization processes influenced chondrite formation. As such, establishing bulk chemical compositions of the pristine Bennu samples is vital for understanding the asteroid, and solar system formation. To begin this processes, we analyzed the bulk major and trace elemental compositions of Bennu aggregates.
Here, this work describes the development of matrix-free GPU-accelerated solvers for high-order finite element problems in H(div). The solvers are applicable to grad-div and Darcy problems in saddle-point formulation, and have applications in radiation diffusion and porous media flow problems, among others. Using the interpolation–histopolation basis, efficient matrix-free preconditioners can be constructed for the (1, 1)-block and Schur complement of the block system. With these approximations, block-preconditioned MINRES converges in a number of iterations that is independent of the mesh size and polynomial degree. The approximate Schur complement takes the form of an M-matrix graph Laplacian and therefore can be well-preconditioned by highly scalable algebraic multigrid methods. High-performance GPU-accelerated algorithms for all components of the solution algorithm are developed, discussed, and benchmarked. Numerical results are presented on a number of challenging test cases, including the “crooked pipe” grad-div problem, the SPE10 reservoir modeling benchmark problem, and a nonlinear radiation diffusion test case.
It is beneficial to use a distributed parameter model for large space structures because the approach minimizes the number of model parameters. Holzer's transfer matrix method provides a useful means to simplify and standardize the procedure for solving the system of partial differential equations. Any large space structures can be broken down into sub-structures with simple elastic and dynamical properties. For each single element, such as beam, tether, or rigid body, we can derive the corresponding transfer matrix. Combining these elements' matrices enables the solution of the global system equations. The characteristics equation can then be formed by satisfying the appropriate boundary conditions. Then natural frequencies and mode shapes can be determined by searching the roots of the characteristic equation at frequencies within the range of interest. This paper applies this methodology, and the maximum likelihood estimation method, to refine the modal characteristics of the NASA Mini-Mast Truss by successively matching the theoretical response to the test data of the truss. The method is being applied to more complex configurations.
The factors that affect the preferred positions of cations in ionic solid solutions were investigated utilizing vibrational spectroscopy. Solid solutions of the sulfate and chromate ions codoped with La(+3) and Ca(+2) in a KBr host lattice were examined as a function of the polyvalent cation concentration. The cation-anion pairing process was found to be random for Ca(+2), whereas the formation of La(+3)-SO4(-2) ion pairs with a C2 sub v bonding geometry is highly preferential to any type of La(+3)-CrO4(-2) ion pair formation. The relative populations of ion pair site configurations are discussed in terms of an energy-entropy competition model which can be applied to the partition of trace elements during magmatic processes.