Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermomechanical 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 235 records · Page 13

Micro-scale Artificial Weave Generation Capabilities for Thermal Protection System Material Modeling

Thermal Protection System (TPS) modeling requires accurate representation and prediction of the thermomechanical behavior of ablative materials. State-of-the-art TPS materials such as Phenolic Impregnated Carbon Ablator (PICA) have a proven flight record and demonstrate exceptional capabilities for handling extreme aerothermal heating conditions. The constant push for lightweight materials that are flexible in their design and performance, and hence allow for a wide range of mission profiles, has led NASA over the past years to develop its Heatshield for Extreme Entry Environment Technology (HEEET). HEEET is based primarily on a dual layer woven carbon fiber architecture and the technology has successfully been tested in arc-jet facilities. These recent developments have sparked interest in the accurate micro-scale modeling of composite weave architectures, to predict the structural response of macro-scale heatshields upon atmospheric entry. This effort can be extended to incorporate in-depth failure mechanics analyses as a result of local thermal gradients or high-velocity particle impact.

Visser, Sander J.↗

Software Architecture and Hierarchy of the Nasa Multiscale Analysis Tool

The NASA Multiscale Analysis Tool (NASMAT) serves as a state-of-the-art, “plug and play,” software package which utilizes multiscale recursive micromechanics as a platform for massively multiscale modeling for hierarchical materials and structures subjected to thermomechanical loads on high performance computing systems. This paper is intended to give an overview of the design of NASMAT and how the design supports modularity, upgradability and maintainability, interoperability, and utility. First, the software architecture and hierarchy will be explored. Details on each of the 11 NASMAT procedures and the arrangement of NASMAT data will be presented. Finally, application program interfaces (APIs) that were developed to facilitate the communication of NASMAT with other programs will be described.

multiscale modeling↗

Toward improved durability in advanced combustors and turbines: Progress in the prediction of thermomechanical loads

NASA is sponsoring the Turbine Engine Hot Section Technology (HOST) Project to address the need for improved durability in advanced combustors and turbines. Analytical and experimental activities aimed at more accurate prediction of the aerothermal environment, the thermomechanical loads, the material behavior and structural responses to such loading, and life predictions for high temperature cyclic operation have been underway for several years and are showing promising results. Progress is reported in the development of advanced instrumentation and in the improvement of combustor aerothermal and turbine heat transfer models that will lead to more accurate prediction of thermomechanical loads.

Sokolowski, Daniel E.↗

A mechanistic damage model for solid oxide fuel cell ceramic materials - Part I: Constitutive modeling

A multiscale mechanistic damage model for solid oxide fuel cell (SOFC) ceramics that combines micromechanics of stiffness reduction due to material porosity change and microcracking with a continuum damage mechanics (CDM) description for the evolution of microcracks up to fracture is developed in this work. The model also accounts for volumetric swelling of the anode due to redox cycling. Porosity change and swelling during redox combined with normal operating thermomechanical loads could increase the stress, strain and damage distributions in SOFC stacks leading them to failure. The proposed model involves three governing parameters: (i) the porosity in terms of the pore or void volume fraction, (ii) material swelling magnitude prescribed at a given loading step, and (iii) the damage variable that describes microcracking caused by thermomechanical loads and swelling. Pores and microcracks in the ceramics are modeled as randomly distributed ellipsoidal inclusions with negligible stiffness by an Eshelby-Mori-Tanaka formulation combined with an inclusion orientation distribution method. Microcracking damage evolution is described by a CDM formulation. Volumetric swelling is treated in a similar way to thermal expansion in the constitutive relations. After validation and numerical checks, the damage model is used to analyze a simplified SOFC positive-electrode/electrolyte/negative-electrode (PEN) structure subjected a redox cycle in addition to thermomechanical loads experienced during normal operation of a SOFC stack.

30 DIRECT ENERGY CONVERSION↗

Development and application of a uranium mononitride (UN) potential: Thermomechanical properties and Xe diffusion

We report that atomic-scale modeling of thermophysical and defect properties of uranium mononitride (UN) plays an important role in establishing a better understanding and improved models of UN fuel performance. Having an accurate interatomic potential is crucial for generating reliable data at finite temperatures using molecular dynamic simulations. We report a new interatomic potential for UN, based on a combination of many-body and pairwise interactions, a simple form that we later show could be easily adapted to include Xe-U and Xe-N interactions, i.e., generating a U-N-Xe interatomic potential. The potential was fitted to experimental thermal expansion and single crystal elastic constants, as well as Frenkel, Schottky, anti-Schottky, and antisite pair reaction energies from density functional theory (DFT) calculations. Using the potential, we successfully reproduced experimental lattice parameters, thermal expansion, single crystal elastic constants, and temperature dependent heat capacity. The potential also performs reasonably well in reproducing the energy of the aforementioned stoichiometric defect reactions and defect migration barriers calculated using DFT. However, the potential underestimates the energy difference between the tetrahedral and dumbbell uranium interstitials, and a more complex potential form might be needed to overcome this issue. The potential was also used to predict UN single crystal elastic constants and elastic properties at different temperatures, showing that UN becomes softer and more compressible with increasing temperature. We also compare our potential against literature data from previous empirical potentials, demonstrating similar or better behavior depending on the property of interest. To enable the simulation of Xe in UN a Buckingham potential has been fitted to DFT-derived Xe incorporation energies. The potential was then used to determine the activation energy for Xe diffusion due to various Xe-containing defects, with {XeU : VU} exhibiting the lowest activation energy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling Framework to Analyze Performance and Structural Reliability of Solid Oxide Electrolysis Cells

Solid oxide electrolysis cells (SOEC) have been receiving significant attention recently because of their high energy efficiency and fast hydrogen production. In this study a multi-physics model to simulate the SOEC performance and structural reliability of a state-of-the-art planar SOEC design was developed. The electrochemical reactions, fluid dynamics, species transport, electron transfer, and heat transfer were modeled in the commercial computational fluid dynamics (CFD) software STAR-CCM+. The thermomechanical analysis and the associated structural reliability evaluations were conducted using the commercial finite element analysis software ANSYS. The electrochemistry model was validated by using the experimentally obtained current-voltage (I-V) characteristics of the electrode-supported SOECs. The reliability analysis using a risk-of-rupture approach showed low failure probabilities under standard operating conditions considered in this study. For cells operated at voltages well above a thermoneutral voltage, the reliability evaluations indicated a potential risk of cell failure, but the damage was concentrated locally in specific areas of the cell which typically do not lead to total loss of cell function. The presented approach provides insights for evaluating representative cell and stack performances and structural reliability without intensive testing and for developing optimally performing and structurally reliable SOECs for efficient hydrogen generation.

25 ENERGY STORAGE↗

High Temperature Thick Film Sensor Development Based on Refractory Oxide Semiconductors

Solid-state, high-temperature sensing devices are required for accurate temperature, strain/stress, and failure monitoring for a range of advanced manufacturing, transportation, and military applications. High temperature conditions limit sensing strategies, where typically traditional metal and semiconductor materials are unstable, and the sensing options are limited to optical spectroscopy methods. In addition, strategies to embed the sensors directly within active components and protective thermal refractory are of interest that permit near-environment sensing which is not possible with the traditional electronic materials and components. Rare-earth chromite semiconductors exhibit some properties of interest for high temperature sensing technologies, such as: high microstructure and chemical stability, relatively high electronic conductivity at high temperatures, and relative thermal expansion coefficient matching to refractory thermal protective insulation (such as high-zirconia and -alumina ceramics). In this work, n-type and p-type doped lanthanum chromites compositions were computationally modelled and prepared, and the electrical, thermoelectric, and thermomechanical properties were measured for the obtained material systems. A few compositions were then processed into thick film health and temperature sensors and tested at high temperatures. To establish relative stability under harsh environment conditions, post- mortem microstructure and chemical composition was characterized after extended and cycle testing at high temperatures.

20 FOSSIL-FUELED POWER PLANTS↗

Computational Simulation of Damage Propagation in Three-Dimensional Woven Composites

Three dimensional (3D) woven composites have demonstrated multi-directional properties and improved transverse strength, impact resistance, and shear characteristics. The objective of this research is to develop a new model for predicting the elastic constants, hygrothermal effects, thermomechanical response, and stress limits of 3D woven composites; and to develop a computational tool to facilitate the evaluation of 3D woven composite structures with regard to damage tolerance and durability. Fiber orientations of weave and braid patterns are defined with reference to composite structural coordinates. Orthotropic ply properties and stress limits computed via micromechanics are transformed to composite structural coordinates and integrated to obtain the 3D properties. The various stages of degradation, from damage initiation to collapse of structures, in the 3D woven structures are simulated for the first time. Three dimensional woven composite specimens with various woven patterns under different loading conditions, such as tension, compression, bending, and shear are simulated in the validation process of this research. Damage initiation, growth, accumulation, and propagation to fracture are included in these simulations.

Huang, Dade↗

Estimation Filter for Alignment of the Spitzer Space Telescope

A document presents a summary of an onboard estimation algorithm now being used to calibrate the alignment of the Spitzer Space Telescope (formerly known as the Space Infrared Telescope Facility). The algorithm, denoted the S2P calibration filter, recursively generates estimates of the alignment angles between a telescope reference frame and a star-tracker reference frame. At several discrete times during the day, the filter accepts, as input, attitude estimates from the star tracker and observations taken by the Pointing Control Reference Sensor (a sensor in the field of view of the telescope). The output of the filter is a calibrated quaternion that represents the best current mean-square estimate of the alignment angles between the telescope and the star tracker. The S2P calibration filter incorporates a Kalman filter that tracks six states - two for each of three orthogonal coordinate axes. Although, in principle, one state per axis is sufficient, the use of two states per axis makes it possible to model both short- and long-term behaviors. Specifically, the filter properly models transient learning, characteristic times and bounds of thermomechanical drift, and long-term steady-state statistics, whether calibration measurements are taken frequently or infrequently. These properties ensure that the S2P filter performance is optimal over a broad range of flight conditions, and can be confidently run autonomously over several years of in-flight operation without human intervention.

Bayard, David↗

Material Characterization and Geometric Segmentation of a Composite Structure Using Microfocus X-Ray Computed Tomography Image-Based Finite Element Modeling

This study utilizes microfocus x-ray computed tomography (CT) slice sets to model and characterize the damage locations and sizes in thermal protection system materials that underwent impact testing. ScanIP/FE software is used to visualize and process the slice sets, followed by mesh generation on the segmented volumetric rendering. Then, the local stress fields around several of the damaged regions are calculated for realistic mission profiles that subject the sample to extreme temperature and other severe environmental conditions. The resulting stress fields are used to quantify damage severity and make an assessment as to whether damage that did not penetrate to the base material can still result in catastrophic failure of the structure. It is expected that this study will demonstrate that finite element modeling based on an accurate three-dimensional rendered model from a series of CT slices is an essential tool to quantify the internal macroscopic defects and damage of a complex system made out of thermal protection material. Results obtained showing details of segmented images; three-dimensional volume-rendered models, finite element meshes generated, and the resulting thermomechanical stress state due to impact loading for the material are presented and discussed. Further, this study is conducted to exhibit certain high-caliber capabilities that the nondestructive evaluation (NDE) group at NASA Glenn Research Center can offer to assist in assessing the structural durability of such highly specialized materials so improvements in their performance and capacities to handle harsh operating conditions can be made.

Abdul-Aziz, Ali↗

Constitutive modeling for isotropic materials

The third and fourth years of a 4-year research program, part of the NASA HOST Program, are described. The program goals were: (1) to develop and validate unified constitutive models for isotropic materials, and (2) to demonstrate their usefulness for structural analysis of hot section components of gas turbine engines. The unified models selected for development and evaluation were those of Bodner-Partom and of Walker. The unified approach for elastic-viscoplastic constitutive equations is a viable method for representing and predicting material response characteristics in the range where strain rate and temperature dependent inelastic deformations are experienced. This conclusion is reached by extensive comparison of model calculations against the experimental results of a test program of two high temperature Ni-base alloys, B1900+Hf and Mar-M247, over a wide temperature range for a variety of deformation and thermal histories including uniaxial, multiaxial, and thermomechanical loading paths. The applicability of the Bodner-Partom and the Walker models for structural applications has been demonstrated by implementing these models into the MARC finite element code and by performing a number of analyses including thermomechanical histories on components of hot sections of gas turbine engines and benchmark notch tensile specimens. The results of the 4-year program have been published in four annual reports. The results of the base program are summarized in this report. The tasks covered include: (1) development of material test procedures, (2) thermal history effects, and (3) verification of the constitutive model for an alternative material.

Chan, K. S.↗

Application of powder densification models to the consolidation processing of composites

Unidirectional fiber reinforced metal matrix composite tapes (containing a single layer of parallel fibers) can now be produced by plasma deposition. These tapes can be stacked and subjected to a thermomechanical treatment that results in a fully dense near net shape component. The mechanisms by which this consolidation step occurs are explored, and models to predict the effect of different thermomechanical conditions (during consolidation) upon the kinetics of densification are developed. The approach is based upon a methodology developed by Ashby and others for the simpler problem of HIP of spherical powders. The complex problem is devided into six, much simpler, subproblems, and then their predicted contributions are added to densification. The initial problem decomposition is to treat the two extreme geometries encountered (contact deformation occurring between foils and shrinkage of isolated, internal pores). Deformation of these two geometries is modelled for plastic, power law creep and diffusional flow. The results are reported in the form of a densification map.

Wadley, H. N. G.↗

Thermomechanical characterization of Hastelloy-X under uniaxial cyclic loading

In most high-temperature engineering applications, components are subjected to complex combinations of thermal and mechanical loading during service. A number of viscoplastic constitutive models were proposed which potentially can provide mathematical descriptions of material response under such conditions. Implementation of these models into large finite element codes such as MARC has already resulted in much improved inelastic analysis capability for hot-section aircraft engine components. However, a number of questions remain regarding the validity of methods adopted in characterizing these constitutive models for particular high-temperature materials. One area of concern is that the majority of experimental data available for this purpose are determined under isothermal conditions. This is in contrast to service conditions which, as noted above, almost always involve some form of thermal cycling. The obvious question arises as to whether a constitutive model characterized using an isothermal data base can adequately predict material response under thermomechanical conditions. An experimental program was initiated within the HOST program to address this particular concern. The results of the most recent isothermal and thermomechanical experiments are described.

Ellis, J. R.↗

Model form and sensitivity analysis of CALPHAD-based nucleation models in b-stabilized Ti alloys

Accurate prediction of α-phase nucleation and growth in β-stabilized titanium alloys is crucial for designing heat treatments to optimize mechanical properties in additively manufactured lightweight components. Ideally, predictions of nucleation and growth would incorporate both top-down observations of past experimental heat treatments and bottom-up modeling of phase transformations; however, the appropriate method of combining these information sources is not self-evident. Combining top-down and bottom-up information requires a unified form of model that can connect between spatiotemporal scales, as well as sets of fitting parameters that can be identified by each data source. The selection of which parameters to fit to which data source can be made based on expert opinion, or by performing a sensitivity analysis. In solid-solid nucleation, direct observation of the nucleation and growth process is challenging. Most data on the heat treatment-controlled phase transformations are not in-situ. To predict the process and outcome of the nucleation, growth and coarsening of precipitates, theoretical models of the nucleation pathway are used to bridge the gap. Many sources of uncertainty affect the modeling of this nucleation process. It can be influenced by small variations in the thermomechanical processing history, chemical composition, and initial microstructure. If molecular dynamics (MD) simulations are used to determine thermodynamic quantities and inform CALPHAD modeling, additional uncertainty can be introduced and accounted for using Bayesian methods. Top-down uncertainties require additional steps to quantify. The influence of nucleation model form on the sensitivity of predictions to input parameters and physical conditions is the focus of this study. Classical nucleation theory (CNT) allows modeling to formulate the nucleation as homogeneous or, more commonly, heterogeneous. Non-classical nucleation models are also increasingly explored as a means of reconciling top-down and bottom-up data. In this study, the sensitivity of the intragranular nucleation of α in a β-annealed, slow-cooled aging (BASCA) heat treatment of β-stabilized Ti5553 alloy is explored using CNT and both heterogeneous and homogeneous assumptions. The Kampmann-Wagner Numerical model of precipitate nucleation and growth is employed. Using open-source tools (pyCalphad and thermodynamic modeling of TiMo as a surrogate system, a sensitivity analysis is performed to measure variations in key parameters, including chemical driving force, interfacial energy, and diffusivity, as they relate to predictions of precipitate number density. The inclusion of top-down and bottom-up data in selection of nucleation model form is discussed.

Rodriguez Negron, A. M.↗

Damage Modeling of Solid Oxide Fuel Cells Accounting for Redox Effects

This work applies a multiscale mechanistic damage model developed for brittle ceramics and implemented in commercial finite element (FE) packages via user subroutines to study progressive damage in solid oxide fuel cells (SOFC) subjected to thermomechanical loading under normal operating and shutdown conditions including redox effects. The damage model captures the micromechanics of stiffness reduction due to material porosity change and microcracking and integrates the as-obtained stiffness reduction law into a continuum damage mechanics (CDM) formulation for the evolution of microcracks up to fracture. The volumetric “swelling” that occurs during redox is treated in constitutive modeling similarly to thermal expansion, but swelling strains are irreversible. Furthermore, this damage model was first validated through predictions of strength and stress-strain response for the SOFC electrode materials. Next, it has been applied to predict the potential for degradation in a generic planar SOFC stack with large active area cells. Multicell stack models were simulated in both co-flow and counter-flow configurations. In addition, a constant temperature redox cycle was also simulated to capture overall cell electrode damage due to volumetric swelling of the nickel (Ni)-based anode in the anode-supported cells.

25 ENERGY STORAGE↗

Material Characterization and Modeling of Room Temperature Vulcanizing Silicone

Room Temperature Vulcanizing silicone (RTV) is a high-temperature adhesive that has successfully been used as a gap-filler between Thermal Protection System (TPS) tiles for heatshields on numerous missions. It is also used to bond instrumentation plugs such as temperature and pressure sensors into the heatshields. While RTV has been traditionally assumed to be a non-porous and non-ablating material, numerous experiments have shown that RTV pyrolyzes and becomes highly porous as it is heated. Heating RTV has also shown swelling, or intumescence, which can pose unique problems that lead to roughness induced boundary-layer transition, surface oxide formation and contamination of heat shield sensors. Therefore, it is crucial to understand and model the intumescence phenomenon of RTV. As data for RTV material properties is limited, the first step in modeling RTV is to collect material properties such as pyrolysis mass-loss, microstructure change, virgin and char porosity, etc. which was performed in our initial study. Additionally, thermomechanical properties such as Young’s modulus and Poisson ratio are required for modeling the intumescence of RTV, which were taken from literature and the coefficient of thermal expansion was collected using in-situ heating and Micro Computed Tomography (µ-CT) in previous studies. Finally, numerous other properties such as pyrolysis gas properties, virgin and char thermal conductivity and specific heat were compiled from previous experiments and literature into a material database that can be used for simulations. In Porous Material Analysis Toolbox based on OpenFOAM (PATO) [4], structural mechanics coupled with material response was used for simulating the intumescence of RTV as it is heated. However, since the permeability of the material is very low, the pyrolysis gas creates an internal pressure build-up as the material is being heated, significantly contributing to the deformation of the material. To correctly characterize this phenomenon, additional physics models were implemented into PATO's stress analysis solver, and results were compared with RTV dilatometry test data as a preliminary verification case. Future work will include experiments of RTV at the Plasmatron X facility and the in-situ heating cell with µ-CT, and improvement of simulation tools to more accurately model RTV intumescence.

TPS↗

Material Properties and Modeling of Room Temperature Vulcanizing Silicone

Room Temperature Vulcanizing silicone (RTV) is a high-temperature adhesive that has successfully been used as a gap-filler between Thermal Protection System (TPS) tiles for heatshields on numerous missions. It is also used to bond instrumentation plugs such as temperature and pressure sensors into the heatshields. While RTV has been traditionally assumed to be a non-porous and non-ablating material, numerous experiments have shown that RTV pyrolyzes and becomes highly porous as it is heated. Heating RTV has also shown swelling, or intumescence, which can pose unique problems that lead to roughness induced boundary-layer transition, surface oxide formation and contamination of heat shield sensors. Therefore, it is crucial to understand and model the intumescence phenomenon of RTV. As data for RTV material properties is limited, the first step in modeling RTV is to collect material properties such as pyrolysis mass-loss, microstructure change, virgin and char porosity, etc. which was performed in our initial study. Additionally, thermomechanical properties such as Young’s modulus and Poisson ratio are required for modeling the intumescence of RTV, which were taken from literature and the coefficient of thermal expansion was collected using in-situ heating and Micro Computed Tomography (µ-CT) in previous studies. Finally, numerous other properties such as pyrolysis gas properties, virgin and char thermal conductivity and specific heat were compiled from previous experiments and literature into a material database that can be used for simulations. In Porous Material Analysis Toolbox based on OpenFOAM (PATO) [4], structural mechanics coupled with material response was used for simulating the intumescence of RTV as it is heated. However, since the permeability of the material is very low, the pyrolysis gas creates an internal pressure build-up as the material is being heated, significantly contributing to the deformation of the material. To correctly characterize this phenomenon, additional physics models were implemented into PATO's stress analysis solver, and results were compared with RTV dilatometry test data as a preliminary verification case. Future work will include experiments of RTV at the Plasmatron X facility and the in-situ heating cell with µ-CT, and improvement of simulation tools to more accurately model RTV intumescence.

PATO↗

Modeling and life prediction methodology for Titanium Matrix Composites subjected to mission profiles

Titanium matrix composites (TMC) are being evaluated as structural materials for elevated temperature applications in future generation hypersonic vehicles. In such applications, TMC components are subjected to complex thermomechanical loading profiles at various elevated temperatures. Therefore, thermomechanical fatigue (TMF) testing, using a simulated mission profile, is essential for evaluation and development of life prediction methodologies. The objective of the research presented in this paper was to evaluate the TMF response of the (0/90)2s SCS-6/Timetal-21S subjected to a generic hypersonic flight profile and its portions with a temperature ranging from -130 C to 816 C. It was found that the composite modulus, prior to rapid degradation, had consistent values for all the profiles tested. A micromechanics based analysis was used to predict the stress-strain response of the laminate and of the constituents in each ply during thermomechanical loading conditions by using only constituent properties as input. The fiber was modeled as elastic with transverse orthotropic and temperature dependent properties. The matrix was modeled using a thermoviscoplastic constitutive relation. In the analysis, the composite modulus degradation was assumed to result from matrix cracking and was modeled by reducing the matrix modulus. Fatigue lives of the composite subjected to the complex generic hypersonic flight profile were well correlated using the predicted stress in 0 degree fibers.

Mirdamadi, M.↗