Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “materials modelling”

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 703 records · Page 39

Full-Scale MSL Heatshield Material Response Using DSMC and CFD to Compute the Aerothermal Environments

During Mars atmospheric entry, the Mars Science Laboratory (MSL) was protected by a 4.5 meters diameter ablative heatshield assembled in 113 tiles [1]. The heatshield was made of NASA's flagship ablative material, the Phenolic Impregnated Carbon Ablator (PICA) [2]. Prior work [3] compared the traditional one-dimensional and three-dimensional material response models at different locations in the heatshield. It was observed that the flow was basically one-dimensional in the nose and flank regions, but three-dimensional flow effects were observed in the outer flank. Additionally, the effects of tiled versus monolithic heatshield models were also investigated. It was observed that the 3D tiled and 3D monolithic configurations yielded relative differences for in-depth material temperature up to 18% and 28%, respectively, when compared to the a 1D model.

Aerothermal↗

Fatigue life prediction modeling for turbine hot section materials

This paper presents a summary of the life prediction methods developed under the NASA Lewis Research Center's Hot Section Technology (HOST) program. A major objective of the fatigue and fracture efforts under the HOST program was to significantly improve the analytic life prediction tools used by the aeronautical gas turbine engine industry. This has been achieved in the areas of high-temperature thermal and mechanical fatigue of bare and coated high-temperature superalloys. Such technical improvements will eventually reduce life cycle costs. The cyclic crack initiation and propagation resistance of nominally isotropic polycrystalline alloys and highly anisotropic single crystal alloys has been addressed. A sizeable data base has been generated for three alloys [cast PWA 1455 (B-1900 + Hf), wrought Inconel 718, and cast single crystal PWA 1480] in bare and coated conditions. Two coatings systems, diffusion aluminide (PWA 273) and plasma sprayed MCrAIY overlay (PWA 286) were employed. Life prediction modeling efforts were devoted to creep-fatigue interaction, oxidation, coatings interactions, multiaxiality of stress-strain states, mean stress effects, cumulative damage, and thermomechanical fatigue. The fatigue crack initiation life models developed to date include the Cyclic Damage Accumulation (CDA) Model of Pratt & Whitney and the Total Strain Version of Strainrange Partitioning (TS-SRP) of NASA Lewis for nominally isotropic materials, and the Tensile Hysteretic Energy Model of Pratt & Whitney for anisotropic superalloys. The fatigue model being developed by the General Electric Company is based upon the concepts of Path-lndependent Integrals (PII) for describing cyclic crack growth under complex nonlinear response at the crack tip due to thermomechanical loading conditions. A micromechanistic oxidation crack extension model has been derived by researchers at Syracuse University. The models are described and discussed in the paper. Only limited verification has been achieved to-date as several of the technical programs are still in progress and the verification tasks are scheduled, quite naturally, near the conclusion of the program. To-date, efforts have concentrated on development of independent models for cyclic constitutive behavior, cyclic crack initiation, and cyclic crack propagation. The transition between crack initiation and crack propagation has not been thoroughly researched as yet, and the integration of these models into a unified life prediction method has not been addressed.

Metal fatigue↗

Particulate erosion mechanisms

Particulate damage and erosion of ductile metals are today plaguing design and field engineers in diverse fields of engineering and technology. It was found that too many models and theories were proposed leading to much speculation from debris analysis and failure mechanism postulations. Most theories of solid particle erosion are based on material removal models which do not fully represent the actual physical processes of material removal. The various mechanisms proposed thus far are: melting, low-cycle fatigue, extrusion, delamination, shear localization, adhesive material transfer, etc. The experimental data on different materials highlighting the observed failure modes of the deformation and cutting wear processes using optical and scanning electron microscopy are presented. The most important mechanisms proved from the experimental observations of the specimens exposed to both spherical and angular particles are addressed, and the validity of the earlier theories discussed. Both the initial stages of damage and advanced stages of erosion were studied to gain a fundamental understanding of the process.

Veerabhadrarao, P.↗

A thermalized ion explosion model for high energy sputtering and track registration

A velocity spectrum of neutral sputtered particles as well as a low resolution mass spectrum of sputtered molecular ions was measured for 4.74 MeV F-19(+2) incident of UF4. The velocity spectrum is dramatically different from spectra taken with low energy (keV) bombarding ions, and is shown to be consistent with a hot plasma of atoms in thermal equilibrium inside the target. A thermalized ion explosion model is proposed for high energy sputtering which is expected to describe track formation in dielectric materials. The model is shown to be consistent with the observed total sputtering yield and the dependence of the yield on the primary ionization rate of the incident ion.

Seiberling, L. E.↗

Multiphysics Nuclear Thermal Rocket Thrust Chamber Analysis

The objective of this effort is t o develop an efficient and accurate thermo-fluid computational methodology to predict environments for hypothetical thrust chamber design and analysis. The current task scope is to perform multidimensional, multiphysics analysis of thrust performance and heat transfer analysis for a hypothetical solid-core, nuclear thermal engine including thrust chamber and nozzle. The multiphysics aspects of the model include: real fluid dynamics, chemical reactivity, turbulent flow, and conjugate heat transfer. The model will be designed to identify thermal, fluid, and hydrogen environments in all flow paths and materials. This model would then be used to perform non- nuclear reproduction of the flow element failures demonstrated in the Rover/NERVA testing, investigate performance of specific configurations and assess potential issues and enhancements. A two-pronged approach will be employed in this effort: a detailed analysis of a multi-channel, flow-element, and global modeling of the entire thrust chamber assembly with a porosity modeling technique. It is expected that the detailed analysis of a single flow element would provide detailed fluid, thermal, and hydrogen environments for stress analysis, while the global thrust chamber assembly analysis would promote understanding of the effects of hydrogen dissociation and heat transfer on thrust performance. These modeling activities will be validated as much as possible by testing performed by other related efforts.

Wang, Ten-See↗

Tangle-Free Mesh Motion for Ablation Simulations

Problems involving mesh motion-which should not be mistakenly associated with moving mesh methods, a class of adaptive mesh redistribution techniques-are of critical importance in numerical simulations of the thermal response of melting and ablative materials. Ablation is the process by which material vaporizes or otherwise erodes due to strong heating. Accurate modeling of such materials is of the utmost importance in design of passive thermal protection systems ("heatshields") for spacecraft, the layer of the vehicle that ensures survival of crew and craft during re-entry. In an explicit mesh motion approach, a complete thermal solve is first performed. Afterwards, the thermal response is used to determine surface recession rates. These values are then used to generate boundary conditions for an a posteriori correction designed to update the location of the mesh nodes. Most often, linear elastic or biharmonic equations are used to model this material response, traditionally in a finite element framework so that complex geometries can be simulated. A simple scheme for moving the boundary nodes involves receding along the surface normals. However, for all but the simplest problem geometries, evolution in time following such a scheme will eventually bring the mesh to intersect and "tangle" with itself, inducing failure. This presentation demonstrates a comprehensive and sophisticated scheme that analyzes the local geometry of each node with help from user-provided clues to eliminate the tangle and enable simulations on a wide-class of difficult problem geometries. The method developed is demonstrated for linear elastic equations but is general enough that it may be adapted to other modeling equations. The presentation will explicate the inner workings of the tangle-free mesh motion algorithm for both two and three-dimensional meshes. It will show abstract examples of the method's success, including a verification problem that demonstrates its accuracy and correctness. The focus of the presentation will be on the algorithm; specifics on how the techniques may be used in spacecraft design will be not discussed.

Droba, Justin↗

The ultraviolet spectrum of an oxygen-rich supernova remnant in the Small Magellanic Cloud

Ultraviolet and optical spectra data are presented for the oxygen-rich supernova remnant 1E 0102-7219 in the SMC. These UV data are the first to seriously constrain the UV emission from an oxygen-rich remnant. Emission lines of O I, forbidden O II, semiforbidden O III, semiforbidden O IV, C IV, forbidden Ne IV, and Mg II have been unambiguously detected in the UV. Shock models for material with abundances appropriate to massive stellar ejecta are calculated, and a number of differences between these models and the observations are found. The UV lines are observed to be weaker relative to optical lines than expected from models. Models of photoionization by the remnant's X-ray emission are calculated, and these provide a better qualitative match to the observed spectra. Approximate abundances in the ejecta are derived and the implications these have for the precursor star are discussed.

Blair, William P.↗

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.↗

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.↗