Engineering PapersSearch

SEARCH · Engineering Papers

Results for “materials modeling”

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 919 records · Page 51

Multiscale Modeling of Ablation and Pyrolysis in PICA-Like materials

During atmospheric entry of planetary probes, the thermal protection system (TIPS) of the probe is exposed to high temperatures under low pressures. In these conditions, carbonous fibrous TIPS materials may undergo oxidation leading to mass loss and wall recession called ablation. This work aims to improve the understanding of material/environment interactions through a study of the coupling between oxygen transport in the Knudsen regime, heterogeneous oxidation of carbon, and surface recession. A 3D Random Walk Monte Carlo simulation tool is used for this study. The fibrous architecture of a model material, consisting of high porosity random array of carbon fibers, is numerically represented on a 3D Cartesian grid. Mass transport in the Knudsen regime from the boundary layer to the surface, and inside this porous material is simulated by random walk. A reaction probability is used to simulate the heterogeneous oxidation reaction. The surface recession of the fibers is followed by front tracking using a simplified marching cube approach. The output data of the simulations are ablation velocity and dynamic evolution of the material porosity. A parametric study is carried out to analyze the material behavior as a function of Knudsen number for the porous media (length of the mean free path compared to the mean pore diameter) and the intrinsic reactivity of the carbon fibers. The model is applied to Stardust mission reentry conditions and explains the unexpected behavior of the TIPS material that underwent mass loss in volume.

Lachaud, Jean

Optimizing acoustical treatment

A mathematical linear programming model is presented for optimizing acoustical treatment and interior decoration of concert and other public halls. This method provides the designer with a range of acoustically correct solutions at increased economical efficiency. The mathematical model uses geometrical data about the room, recommended reverberation time values, the architect's choice of given sound absorbing structures and finishing materials. The model permits inclusion of aesthetical considerations about conditioning, proportioning, or, on the contrary, reciprocal exclusion of any classes of material and/or sound absorbing structure.

Beuran, N.

Anisotropic constitutive modeling for nickel-base single crystal superalloys

An anisotropic constitutive model was developed based on crystallographic slip theory for nickel base single crystal superalloys. The constitutive equations developed utilizes drag stress and back stress state variables to model the local inelastic flow. Specially designed experiments were conducted to evaluate the existence of back stress in single crystal superalloy Rene N4 at 982 C. The results suggest that: (1) the back stress is orientation dependent; and (2) the back stress state variable is required for the current model to predict material anelastic recovery behavior. The model was evaluated for its predictive capability on single crystal material behavior including orientation dependent stress-strain response, tension/compression asymmetry, strain rate sensitivity, anelastic recovery behavior, cyclic hardening and softening, stress relaxation, creep and associated crystal lattice rotation. Limitation and future development needs are discussed.

Sheh, Michael Y.

Surrogate modeling of Monte Carlo radiation transport with convolutional neural networks for shielding optimization

Here, we present a machine learning (ML)-based surrogate model using convolutional neural networks (CNN) designed to emulate the attenuation of neutron fields as they pass through various shielding materials. This model can compute the outgoing neutron flux almost instantaneously and achieves reasonable accuracy compared to traditional Monte Carlo (MC)-based codes, which are computationally intensive. This emulator alleviates the complexity of neutron radiation transport through shielding materials by reducing the dimensionality and enables shielding optimization for a known radiation environment. This optimization process, which would have taken an unrealistic timeline due to several complex radiation transport simulations, can now be achieved in minutes, thus increasing computational capabilities in radiation shielding assessment. We demonstrate the applications of this emulator in computing effective dose rates and optimizing shielding solutions for a heavy-ion accelerator facility, such as the Facility for Rare Isotope Beams, where secondary neutrons produced via beam interactions dominate the radiation environment.

accelerator shielding

Computational modelling of polymers

Polymeric materials and polymer/graphite composites show a very diverse range of material properties, many of which make them attractive candidates for a variety of high performance engineering applications. Their properties are ultimately determined largely by their chemical structure, and the conditions under which they are processed. It is the aim of computational chemistry to be able to simulate candidate polymers on a computer, and determine what their likely material properties will be. A number of commercially available software packages purport to predict the material properties of samples, given the chemical structures of their constituent molecules. One such system, Cerius, has been in use at LaRC. It is comprised of a number of modules, each of which performs a different kind of calculation on a molecule in the programs workspace. Particularly, interest is in evaluating the suitability of this program to aid in the study of microcrystalline polymeric materials. One of the first model systems examined was benzophenone. The results of this investigation are discussed.

Celarier, Edward A.

Spectrophotometry and organic matter on Iapetus. 1: Composition models

Iapetus shows a greater hemispheric albedo asymmetry than any other body in the solar system. Hapke scattering theory and optical constants measured in the laboratory are used to identify possible compositions for the dark material on the leading hemisphere of Iapetus. The materials considered are poly-HCN, kerogen, Murchison organic residue, Titan tholin, ice tholin, and water ice. Three-component mixtures of these materials are modeled in intraparticle mixture of 25% poly-HCN, 10% Murchison residue, and 65% water ice is found to best fit the spectrum, albedo, and phase behavior of the dark material. The Murchison residue and/or water ice can be replaced by kerogen and ice tholin, respectively, and still produce very good fits. Areal and particle mixtures of poly-HCN, Titan tholin, and either ice tholin or Murchison residue are also possible models. Poly-HCN is a necessary component in almost all good models. The presence of poly-HCN can be further tested by high-resolution observations near 4.5 micrometers.

Wilson, Peter D.

Tethers as Debris: Simulating Impacts of Kevlar Tethers on Shuttle Tiles

In a previous paper I examined the effects of impacts of polymer tethers on aluminum plates using the SPHC hydrodynamic code. In this paper I apply tether models to a new target - models of Space Shuttle tiles developed during the STS 107 accident investigation. In this three-dimensional simulation, a short tether fragment strikes a single tile supported on an aluminum backing plate. A tile of the LI-900 material is modeled. Penetration and damage to the tile and the backwall are characterized for three normal impact velocities. The tether is modeled as a bundle of eight 1-mm strands, with the bundle having dimensions 2-mm x 4-mm x 20-cm. The bulk material properties used are those of Kevlar(TradeMark) 49, for which a Mie-Gruneisen multiphase equation of state (eos) is used. In addition, the strength model is applied in a linear sense, such that tensile loads along the strand length are supported, but there is no strength in the lateral directions. Tile models include the various layers making up the tile structure. The outermost layer is a relatively dense borosilicate glass, known as RCG, 0.5-mm thick. The RCG layer is present on the top and four sides of the tile. Below this coating is the bulk of the tile, 1.8- in thick, made of LI-900, a product consisting of rigidized fiberous silica with a density of 9 lWft3. Below the main insulating layer is a bottom layer of the same material that has been treated to increase its density by approximately 69% to improve its strength. This densified layer is bonded to a Strain Isolation Pad (SIP), modeled as a refractory felt fabric. The SIP is bonded to an aluminum 2024 wall 0.1-in thick. The tile and backwall materials use a Me-Gruneisen multiphase eos, with the exception of the SIP felt, which uses a fabric equation of state. Fabrics must be crushed to the full bulk material density before bulk material properties and a Mie-Gruneisen eos are applied. Tether fragment impact speeds of 3,7, and 10 km/s are simulated, with impact velocities normal to the tile face. Damage results are presented in tabular format.

Evans, Steven W.

Micromechanical simulation of damage progression in carbon phenolic composites

Carbon/phenolic composites are used extensively as ablative insulating materials in the nozzle region of solid rocket motors. The current solid rocket motor (RSRM) on the space shuttle is fabricated from woven rayon cloth which is carbonized and then impregnated with the phenolic resin. These plies are layed up in the desired configuration and cured to form the finished part. During firing, the surface of the carbon/phenolic insulation is exposed to 5000 F gases from the rocket exhaust. The resin pyrolizes and the material chars to a depth which progresses with time. The rate of charring and erosion are generally predictable, and the insulation depth is designed to allow adequate safety margins over the firing time of the motor. However, anomalies in the properties and response of the carbon/phenolic materials can lead to severe material damage which may decrease safety margins to unacceptable levels. Three macro damage modes which were observed in fired nozzles are: ply lift, 'wedge out', and pocketing erosion. Ply lift occurs in materials with plies oriented nearly parallel to the surface. The damage occurs in a region below the charred material where material temperatures are relatively low - about 500 F. Wedge out occurs at the intersection of nozzle components whose plies are oriented at about 45 deg. The corner of the block of material breaks off along a ply interface. Pocketing erosion occurs in material with plies oriented normal to the surface. Thermal expansion is restrained in two directions resulting in large tensile strains and material failure normal to the surface. When a large section of material is removed as a result of damage, the insulation thickness is reduced which may lead to failure of the nozzle due to excessive heating of critical components. If these damage events cannot be prevented with certainty, the designer must increase the thickness of the insulator thus adding to both weight and cost. One of the difficulties in developing a full understanding of these macro damage mechanisms is that the loading environment and the material response to that environment are extremely complex. These types of damage are usually only observed in actual motor firings. Therefore, it is difficult and expensive to evaluate the reliability of new materials. Standard material tests which measure mechanical and thermal properties of test specimens can only provide a partial picture of how the material will respond in the service environment. The development of the ANALOG test procedure which can combine high heating rates and mechanical loads on a specimen will improve the understanding of the interactive effects of the various loads on the system. But a mechanistic model of material response which can account for the heterogeneity of the material, the progression of various micromechanical damage mechanisms, and the interaction of mechanical and thermal stresses on the material is required to accurately correlate material tests with response to service environments. A model based on fundamental damage mechanisms which is calibrated and verified under a variety of loading conditions will provide a general tool for predicting the response of rocket nozzles. The development of a micromechanical simulation technique was initiated and demonstrated to be effective for studying across-ply tensile failure of carbon/phenolic composites.

Slattery, Kerry T.

Model for thermal stress resistance of truly elastic materials containing more than one crack

The model was developed upon the physical properties of surface energy and intrinsic modulus of elasticity of a material containing a number of equal sized microcracks which are independent of one another. The effect of these cracks upon the strain energy per unit volume of material necessary to continue simultaneous crack growth as well as the measured physical properties was established, and the thermal stress resistance is developed in terms of this energy. The model is expressed in its final form in terms of the measured physical properties of fracture strength, effective modulus of elasticity, and coefficient of thermal expansion. The model was applied to existent thermal stress data of ceramic materials for which these physical properties had been measured. On the basis of these data it was concluded that the thermal stress resistance of a material may be improved by increasing the fracture strength.

Lineback, L. D.

Physics-based Entry, Descent and Landing Risk Model

A physics-based risk model was developed to assess the risk associated with thermal protection system failures during the entry, descent and landing phase of a manned spacecraft mission. In the model, entry trajectories were computed using a three-degree-of-freedom trajectory tool, the aerothermodynamic heating environment was computed using an engineering-level computational tool and the thermal response of the TPS material was modeled using a one-dimensional thermal response tool. The model was capable of modeling the effect of micrometeoroid and orbital debris impact damage on the TPS thermal response. A Monte Carlo analysis was used to determine the effects of uncertainties in the vehicle state at Entry Interface, aerothermodynamic heating and material properties on the performance of the TPS design. The failure criterion was set as a temperature limit at the bondline between the TPS and the underlying structure. Both direct computation and response surface approaches were used to compute the risk. The model was applied to a generic manned space capsule design. The effect of material property uncertainty and MMOD damage on risk of failure were analyzed. A comparison of the direct computation and response surface approach was undertaken.

descent and landing

Solar White Thermal Coating for Cryogenic Propulsion Systems

NASA is currently conducting research into the potential of storing cryogenic fluid in low Earth orbit (LEO). Having cryogenic propellant readily available for high-performance propulsion systems can be very beneficial for deep space missions in the near future. One of the key challenges to storing cryogenic fluid in LEO is minimizing boil-off. To address the challenge, NASA is evaluating new concepts in thermal insulation. One recent experimental study evaluated the feasibility of using Yttrium Oxide (Y2O3), formed into tiles or spray coating that can potentially be used as a thermal coating for cryogenic propellant storage applications in deep space. Due to its temperature and wavelength dependent optical properties, this “solar white” material can reflect a vast majority of the Sun’s radiative energy while having a very high infrared emissivity for rejecting heat to deep space. As a part of the material development and proof of concept testing, multiple tests have been run at KSC and GRC to demonstrate the performance of the material. At GRC, the tests were run using the Deep Space Solar Simulator (DS3) which contains a thermal vacuum chamber in which the solar white sample was exposed to a deep space environment (< 10 K, optically dense walls) while under full illumination via solar lamp. In order to improve the use of the test results and apply them to spacecraft, there is a need to be able to model the material properties within NASA’s standard thermal modeling tools. As such, it was set out to verify a thermal model of one of the experiments using Thermal Desktop with Sinda. The thermal model of the DS3 test setup includes the solar white sample along with a solar lamp simulator capable of outputting heat at varying wavelengths. The model was developed in order to validate the test results and also to help predict results that will be obtained in future tests. This paper will review the modeling methods and thermal analysis results for various test cases that were run within the DS3 facility.

Cryogenic

A Risk-Based Approach for Aerothermal/TPS Analysis and Testing

The current status of aerothermal and thermal protection system modeling for civilian entry missions is reviewed. For most such missions, the accuracy of our simulations is limited not by the tools and processes currently employed, but rather by reducible deficiencies in the underlying physical models. Improving the accuracy of and reducing the uncertainties in these models will enable a greater understanding of the system level impacts of a particular thermal protection system and of the system operation and risk over the operational life of the system. A strategic plan will be laid out by which key modeling deficiencies can be identified via mission-specific gap analysis. Once these gaps have been identified, the driving component uncertainties are determined via sensitivity analyses. A Monte-Carlo based methodology is presented for physics-based probabilistic uncertainty analysis of aerothermodynamics and thermal protection system material response modeling. These data are then used to advocate for and plan focused testing aimed at reducing key uncertainties. The results of these tests are used to validate or modify existing physical models. Concurrently, a testing methodology is outlined for thermal protection materials. The proposed approach is based on using the results of uncertainty/sensitivity analyses discussed above to tailor ground testing so as to best identify and quantify system performance and risk drivers. A key component of this testing is understanding the relationship between the test and flight environments. No existing ground test facility can simultaneously replicate all aspects of the flight environment, and therefore good models for traceability to flight are critical to ensure a low risk, high reliability thermal protection system design. Finally, the role of flight testing in the overall thermal protection system development strategy is discussed.

Wright, Michael J.

Accumulation of sedimentary debris in the south polar region of Mars and implications for climate history

A detailed modeling of the observed crater populations in stratigraphic unit mappings of the Martian south pole region is presently used to establish absolute constraints on the age of materials' emplacement. Modeling results indicate that cratered terrains poleward of 65 deg south latitude have subsequently experienced a steady-state net accumulation of materials at about 0.1 km/Ga; equatorial cratered terrains have by contrast retained a virtually pristine form.

Plaut, J. J.

Finite element analysis of a composite material interface

A finite element model of a composite material interface is developed to study the influence of the interface on the thermal strain in the composite. A plane stress model is used with an axisymmetric model as a check. The interface thickness, thermal coefficient, modulus, Poisson's ratio and the percent of mineral in the composite are variables in the study. The results confirmed the usability of the finite element model in studying the polymer-mineral interface.

Murray, K. H.

Effect of various approximations on predicted progressive failure in plain weave composites

Three-dimensional finite element analysis was used to simulate progressive failure of a plain weave composite subjected to in-plane extension. The loading was parallel to one of the tow directions. The effects of various characteristics of the finite element model on predicted behavior were examined. The predicted behavior was found to be sensitive to quadrature order, mesh refinement, and the material degradation model. Also the sensitivity of the predictions to the tow waviness was studied. The predicted strength decreased considerably with increased waviness. More numerical studies and comparisons with experimental data are needed to establish reliable guidelines for accurate progressive failure prediction.

Whitcomb, John

Analysis of Pyrolysis Products from Ablative Thermal Protection Systems

NASA’s state-of-the-art ablative materials are composed of a three-dimensional network of carbon fibers impregnated with polymers that dissipate thermal energy through pyrolysis. A fundamental understanding of the decomposition mechanisms and pyrolysis product distributions of various classes of polymers is instrumental in the design of new ablative materials. Furthermore, innovative experiments are essential to the continuous modernization of material response models by providing high-fidelity data. Thus, an apparatus has been designed to measure pyrolysis products from polymers and composite materials by implementing in-situ mass spectrometric techniques. Initial results from experiments performed on siloxane resins and a common phenolic resin will be discussed. Both classes of polymers exhibit heating-rate-dependent decomposition mechanisms. At the onset of heating, phenolic polymers decompose through competitive reactions to form gaseous products and a carbonaceous char. Gas phase products of phenolic resins are typically composed of molecular hydrogen, water, and aromatic hydrocarbons. Pyrolysis products from siloxane polymers include molecular hydrogen, small molecules, and cyclic oligomers from the polymer backbone.

Pyrolysis

Analysis of Pyrolysis Products from Ablative TPS

NASA’s state-of-the-art ablative materials are composed of a three-dimensional network of carbon fibers impregnated with polymers that dissipate thermal energy through pyrolysis. A fundamental understanding of the decomposition mechanisms and pyrolysis product distributions of various classes of polymers is instrumental in the design of new ablative materials. Furthermore, innovative experiments are essential to the continuous modernization of material response models by providing high-fidelity data. Thus, an apparatus has been designed to measure pyrolysis products from polymers and composite materials by implementing in-situ mass spectrometric techniques. Initial results from experiments performed on siloxane resins and a common phenolic resin will be discussed. Both classes of polymers exhibit heating-rate-dependent decomposition mechanisms. At the onset of heating, phenolic polymers decompose through competitive reactions to form gaseous products and a carbonaceous char. Gas phase products of phenolic resins are typically composed of molecular hydrogen, water, and aromatic hydrocarbons. Pyrolysis products from siloxane polymers include molecular hydrogen, small molecules, and cyclic oligomers from the polymer.

Pyrolysis

Active learning of a crystal plasticity flow rule from discrete dislocation dynamics simulations

Continuum-scale material deformation models, such as crystal plasticity (CP), can significantly enhance their predictive accuracy by incorporating input from lower-scale (i.e. mesoscale) models. The procedure to generate and extract the relevant information is however typically complex and ad hoc, involving decision and intervention by domain experts, leading to long development times. In this study, we develop a principled approach for calibration of continuum-scale models using lower scale information by representing a CP flow rule as a Gaussian process model. This representation allows for efficient parameter space exploration, guided by the uncertainty embedded in the model through a process known as Bayesian optimization (BO). We demonstrate a semi-autonomous BO loop which instantiates discrete dislocation dynamics simulations whose initial conditions are automatically chosen to optimize the uncertainty of a model CP flow rule. Our self-guided computational pipeline efficiently generated a dataset and corresponding model whose error, uncertainty, and physical feature sensitivities were validated with comparison to an independent dataset four times larger, demonstrating a valuable and efficient active learning implementation readily transferable to similar material systems.

36 MATERIALS SCIENCE