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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

SSME Bearing and Seal Tester Data Compilation, Analysis and Reporting; and Refinement of the Cryogenic Bearing Analysis Mathematical Model

SRS and NASA/MSFC have developed software with unique capabilities to couple bearing kinematic modeling with high fidelity thermal modeling. The core thermomechanical modeling software was developed by SRS and others in the late 1980's and early 1990's under various different contractual efforts. SRS originally developed software that enabled SHABERTH (Shaft Bearing Thermal Model) and SINDA (Systems Improved Numerical Differencing Analyzer) to exchange data and autonomously allowing bearing component temperature effects to propagate into the steady state bearing mechanical model. A separate contract was issued in 1990 to create a personal computer version of the software. At that time SRS performed major improvements to the code. Both SHABERTH and SINDA were independently ported to the PC and compiled. SRS them integrated the two programs into a single program that was named SINSHA. This was a major code improvement.

Moore, James↗
NASA NTRS2000

SSME Bearing and Seal Tester Data Compilation, Analysis, and Reporting; and Refinement of the Cryogenic Bearing Analysis Mathematical Model

SRS and NASA/MSFC have developed software with unique capabilities to couple bearing kinematic modeling with high fidelity thermal modeling. The core thermomechanical modeling software was developed by SRS and others in the late 1980's and early 1990's under various different contractual efforts. SRS originally developed software that enabled SHABERTH (Shaft Bearing Thermal Model) and SINDA (Systems Improved Numerical Differencing Analyzer) to exchange data and autonomously allowing bearing component temperature effects to propagate into the steady state bearing mechanical model. A separate contract was issued in 1990 to create a personal computer version of the software. At that time SRS performed major improvements to the code. Both SHABERTH and SINDA were independently ported to the PC and compiled. SRS them integrated the two programs into a single program that was named SINSHA. This was a major code improvement.

Moore, James↗
NASA NTRS2018

Coupled Thermomechanical Micromechanics Modeling of the Influence of Thermally Grown Oxide Layer in an Environmental Barrier Coating System

Environmental Barrier Coatings (EBCs) have emerged as a promising means of protecting critical components for high temperature applications (e.g., aircraft engines). EBCs are often used to protect an underlying material (substrate) from extreme thermal/chemical environments. However, systems that utilize EBCs are susceptible to a number of failure modes including oxidation/delamination. Environmental Barrier Coatings (EBCs) have emerged as a promising means of protecting silicon based ceramic matrix composite (CMC) components for high temperature applications (e.g., aircraft engines). EBCs are often used to protect an underlying material (substrate) such as silicon carbide from extreme thermal/chemical environments. In a typical CMC/EBC system, an EBC may or may not be adhered to an underlying substrate with a bond coat (e.g., silicon). Irrespective, systems that utilize EBCs are susceptible to a number of failure modes including oxidation/delamination, recession, chemical attack and dissolution, thermo-mechanical degradation, erosion, and foreign object damage. Current work at NASA Glenn Research Center is aimed at addressing these failure modes in EBC systems and developing robust analysis tools to aid in the design process. The Higher-Order Theory for Functionally Graded Materials (HOTFGM), a precursor to the High-Fidelity Generalized Method of Cells micromechanics approach, was developed to investigate the coupled thermo-mechanical behavior of functionally graded composites and will be used herein to assess the development and growth of a low-stiffness thermally grown oxide (TGO) layer in EBC/CMC systems without a silicon bond coat. To accomplish this a sensitivity study is conducted to examine the influence of uniformly and non-uniformly grown oxide layer on the associated driving forces leading to mechanical failure (spallation) of EBC layer when subjected to isothermal loading, recession, chemical attack and dissolution, thermomechanical degradation, erosion, and foreign object damage. Current work at NASA Glenn Research Center is aimed at addressing these failure modes in EBC systems and developing robust analysis tools to aid in the design process. The Higher-Order Theory for Functionally Graded Materials (HOTFGM), a precursor to the High-Fidelity Generalized Method of Cells micromechanics approach, was developed to investigate the coupled thermo-mechanical behavior of functionally graded composites (Aboudi et al., 1999, Composites B). For example, HOTFGM was previously used (Arnold et al, 1995, NASA CP 10178, paper 34), to assess interlaminar stresses (including free edge effects) in a substrate with a thermal barrier coating (TBC). In this study, HOTFGM micromechanics analyses will be used to assess the development and growth of a low-stiffness thermally grown oxide (TGO) layer between a silicon carbide substrate and a ytterbium disilicate EBC. In order to realistically simulate TGO growth, an evolution law will be incorporated into HOTFGM. In addition, the effect of TGO roughness will be explored consistent with previous TBC work (Pindera et al., 2000. Material Science and Engineering, A284, pp. 158-175). This model represents a first step in developing a robust analysis tool that can ultimately be used to design durable EBC systems. Additional failure modes will be considered as part of a future work.

environmental barrier coating↗
DOE OSTI2020

An Improved Fracture Mechanics-Informed Multiscale Thermomechanical Damage Model for Ceramic Matrix Composites

This paper extends recent work done by the authors in modeling length scale-dependent damage behavior of ceramic matrix composites (CMCs) to include effects of local anisotropy introduced by matrix cracking. This model captures scale-dependent damage initiation and propagation behavior of the brittle matrix by employing internal state variable (ISV) theory within a multiscale modeling framework to obtain damaged matrix stress/strain constitutive relationships at each length scale. The damage ISV captures the effects of matrix cracking and growth by using fracture mechanics and the self-consistent scheme to determine the reduced stiffness of the cracked matrix. Matrix cracks, which activate when stress intensity factors near manufacturing induced cavities exceed the fracture toughness of the material, are assumed to be transversely isotropic in the plane of the crack, and matrix anisotropy occurs when the damaged stiffness tensor is rotated from the crack plane to the global axes. The crack progression and temporal evolution of the damage ISV are governed by fracture mechanics and crack growth kinetics. The model effectively captures first matrix cracking, which is the first significant deviation from linear elasticity. The nonlinear predictive capabilities of the material model are demonstrated for monolithic silicon carbide (SiC) and a 2D woven five-harness satin (5HS) carbon fiber SiC matrix (C/SiC) CMC.

Skinner, Travis↗
DOE OSTI2020

A CALPHAD-informed approach to modeling constituent redistribution in Zr-based metallic fuels using BISON

Here, a CALPHAD-informed (Computer Coupling of Phase Diagrams and Thermochemistry) constituent redistribution model was developed for Zr-based metallic fuels and incorporated into the BISON fuel performance code. Three uncertain model parameters associated with β and γ phase kinetics were calibrated using integral test data from U-Zr fuel elements irradiated in Experimental Breeder Reactor II. The calibrated constituent redistribution model was shown to predict the behavior of U-Zr fuels with excellent accuracy. Model predictions for U-Pu-Zr fuels were physically reasonable but less accurate. Reduction of uncertainties in the ternary phase transition temperatures and collection of kinetic data for the ζ phase are expected to improve the model’s ternary predictions. Finally, the new model was coupled to existing thermomechanics models in BISON to simulate irradiation of an entire U-Zr fuel element, demonstrating its ability to accurately predict the behavior of U-Zr fuels with realistic geometries and mesh resolutions at the engineering scale.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Reliability and Lifetime Prediction Model of Sintered Silver Under High-Temperature Cycling

Although excellent reliability has been reported for sintered silver as a die-attach material under both thermal and power cycling loads in power electronics applications, the promise of this material as a large-area attachment at temperatures beyond 200 degrees C needs to be investigated. This paper presents insights into the thermomechanical behavior and reliability of sintered silver under extreme thermal cycling conditions. In this study, we bonded sintered silver samples and subjected it to a thermal cycling profile of -40 °C to 200 °C with high ramp rates. We periodically monitored samples under thermal cycling to detect the presence of any failure mechanisms using a scanning acoustic microscope. We also included 95Pb5Sn solder in the study to obtain reference data. Results show the occurrence of cracks in sintered silver followed by a rapid rate of crack growth that exceeded the failure criterion in just 50 cycles. The predominant failure mechanism we observed was adhesive failure. As a large-area attachment, solder exhibited a higher reliability than sintered silver but failed within 100 cycles. Finally, we performed thermomechanical modeling to compute strain energy density values and correlated these with the experimentally observed crack growth rates to formulate a lifetime prediction model for sintered silver.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2022

Bonded Interfaces for High-Temperature Power Electronics

Wide-bandgap devices have pushed the operational limit of semiconductor devices in automotive power electronics packages to higher temperatures (>200 degrees C). While a higher efficiency can be achieved through the use of wide-bandgap devices, the entire package must be re-designed with components that can withstand the higher temperature limits. Sintered silver and transient liquid phase bonds are potential candidates as bonded materials for use in higher temperatures however, the underlying mechanics of deformation under thermal loads and the resulting failure mechanisms in these materials are not well understood. Furthermore, thermomechanical modeling plays a key role in understanding the degradation behavior of bonded interface materials and for design optimization of power electronics packages.

ADVANCED PROPULSION SYSTEMS,DIRECT ENERGY CONVERSI↗
DOE OSTI2022

Implement and Test 3D Mortar Contact in BISON

We leverage the extension of the generation of mortar segment meshes to three dimensions in MOOSE’s framework to extend thermomechanical modeling capabilities to problems with three dimensions. A modular approach to gap heat transfer physics using the mortar finite element method was created and documented, mechanical contact was extended to three dimensions—including frictional behavior, performance and ease of use were improved, and steps towards scalability of solid mechanics problems involving contact were taken. Many of these new developments are demonstrated in the simulation of 3D light-water reactor (LWR) problems, where the thermomechanical interface problem is solved using the mortar finite element method. Usage of the mortar framework has improved convergence in 2D problems and has enabled employing friction in 3D problems, of which we show results of a short, local stack of 3D pellets. Consequently, the benefits of mortar in terms of solution convergence and quality are extended to three dimensions. Section 2 discusses fundamental developments that enabled the simulation of practical mortar problems in three dimensions and other general improvements, including the reduction of the derivative container size, the modification of dual basis computations when edge dropping (lack of secondary element projection) takes place, the improvement of conditioning when employing the VCP in-edge dropping conditions, and code usability and quality improvements. These latter code enhancements include the migration of tests using “old” mortar contact constraints to using dual mortar with a semi-smooth Newton solution strategy and the reuse of lower dimensional domains for straightforwardly setting up a mortar thermomechanical LWR problem, i.e. the MOOSE action is employed for mechanical contact and the thermal LWR action is employed to capture the gas conductance, contact, and radiation components of gap heat transfer physics. Independently of the mortar LWR thermal action, we developed a modular approach to gap heat transfer that resides in MOOSE and can be leveraged, e.g., in metallic fuel problems. This approach, whose code design based on MOOSE’s user objects to model specific physics was proposed by the maintenance activity, is detailed in Section 3. Based on the dual mortar finite element method, the frictional contact constraints were extended to three dimensions. A block sheared in two directions in and out of contact with a rigid plane is employed in Section 4 to show the way the approach handles changes in frictional states (e.g. stick to slip) within a competitive number of Newton iterations. Equations and numerical results on the use of the VCP with Cartesian Lagrange multipliers, whose combination enables their direct condensation, are described in Section 5.3. Two-dimensional and three-dimensional BISON LWR simulations are discussed in Section 6. Particularly, a stack of five eccentric pellets with a surface defect is simulated and the effect of pellet-cladding friction is assessed. Finally, conclusions are outlined in Section 7.

42 ENGINEERING↗
DOE OSTI2025

Multiphysics Simulation of KRUSTY Warm Critical Experiments Using MOOSE Tools

A series of critical experiments were performed during the Kilopower Reactor Using Stirling TechnologY (KRUSTY) project, producing valuable and unique data for the validation of corresponding microreactor models. In this paper, the development of a high-fidelity multiphysics model of the KRUSTY microreactor within the MOOSE ecosystem is presented, with a focus on the microreactor configuration during the warm critical experiments. The model was established by coupling a deterministic neutronics model using the Griffin code, with a BISON thermomechanical model. The multiphysics model successfully simulated the KRUSTY 15 ¢ and 30 ¢ warm critical reactivity insertion experiments in satisfactory agreement with the experimental measurement data. In conclusion, the modeling results demonstrated the performance of MOOSE-based multiphysics tools in microreactor simulations, and spotlighted the useful applications of these tools in supporting microreactor design, qualification, and licensing activities.

KRUSTY↗
DOE OSTI2021

Coupled physics simulation of fracture in nuclear fuel pellets induced by resistive heating

While it is well known that the UO 2 fuel used in light-water reactors cracks early in its life in the reactor, it is very difficult to collect in-reactor data regarding fuel cracking initiation and propagation to validate fracture models in fuel performance codes. Most fracture data comes from the post-irradiation examination of fuel that has been subjected to long and complex power histories. Replicating the thermal gradients experienced in this fuel to study crack initiation and propagation in a laboratory environment is challenging. This is primarily because the fuel is heated volumetrically by fission, which is difficult to replicate. One way to approximate this volumetric heating out of the reactor is to use electrical resistance (or Joule heating). An electrical conductivity model has been developed and coupled with the existing thermomechanical modeling capability in the BISON nuclear fuel performance code. This will permit simulation of experiments that use resistive heating for the purpose of validating fracture models in BISON. The modeling approach developed here is applied to simulate both historical and current resistive heating experiments on ceramic nuclear fuel.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Crystal plasticity including a phase-field deformation twinning model for the high-rate deformation of cyclotetramethylene tetranitramine

A new finite strain thermomechanical model for the high-rate deformation of the β-polymorph of cyclotetramethylene tetranitramine (β-HMX) has been developed and applied to simulations of plate impact experiments. The crystal plasticity model is based on a model developed previously for RDX (Luscher et al. 2017), which is extended to incorporate deformation twinning. Twinning during normal plate impacts is simulated with a phase-field twin model. First, material parameters governing the kinetics of dislocation slip are calibrated on the subset of simulations which had negative Schmid factors for the twin system. Second, a parametric study of the twin material parameters was performed to find suitable values. Overall, the results of the simulations with the phase-field twinning model are reported for impacts on several crystal orientations. We find that the twin growth decreases with increasing distance from the impact surface because of dissipation of the shock front via dislocation-mediated plasticity, and that the simulated interface velocity with and without phase-field twinning do not show appreciable differences. These modeling results suggest that the significance of twinning in β-HMX cannot be determined with traditional loading configurations and diagnostics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Modeling The Effects of Loading Scenario and Thermal Expansion Coefficient on Potential Failure of Cryo-compressed Hydrogen Vessels

A multiscale thermomechanical model for a simplified Type-3 cryogenic compressed hydrogen (H2) storage vessel is developed in this paper. The model accounts for the temperature-dependent elastic-plastic behavior of the vessel carbon/epoxy composite overwrap and aluminum alloy liner. The homogenized thermo-elastic-plastic behavior for the individual laminas of the vessel layup is obtained by an incremental Eshelby-Mori-Tanka approach associated with a micromechanical failure criterion to predict lamina failure while a standard elastic-plastic constitutive model is used to describe the behavior of a typical aluminum alloy assumed for the liner. The vessel response to external loadings is achieved by a finite element method. Four loading scenarios representing four thermomechanical cycles applied to the vessel are analyzed to evaluate constituent and lamina stresses as well as the associate failure criterion during the cycle according to these scenarios. The model can provide helpful guidance to mitigate thermal stresses by an adequate selection of loading scenario, optimizing the layup and by tailoring thermomechanical properties of the resin matrix.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Thermal shock resistance of lightweight cements developed for geothermal conditions

Cements are a critical component in well construction, as they act to prevent well fluid and gas escape, prevent corrosion of the casing, and strengthen the wellbore to prevent deformation. Under the high temperature/pressure conditions common in geothermal systems, the injection of cold water for energy production is expected to induce cyclic damage to the borehole cement through the rapid temperature fluctuations. These “thermal shocks” are expected to cause casing shrinkage, annulus formation, and cement tensile stresses. To understand the effect of cold water injection on the wellbore environment, a set of rock-cement-steel samples were created to simulate the structure of a geothermal well. Lightweight thermally-insulating cement blends were tested under thermal shock conditions in this study. In each test, the samples were pressurized to an effective pressure of ~3.5 MPa and placed at high temperatures. Thermal shocks were performed by injecting cold water (~10-15 °C) through the samples at a constant rate while keeping the samples at high temperatures until the sample temperature stopped decreasing and deformation ceased. Eight thermal shock tests were conducted with each sample – two at 100 °C and six at 200 °C. Post-tests analysis was then conducted by cutting open each sample to examine the damage in each component of the simulated wellbore. Experimental results suggest that all samples experienced similar degrees of axial and lateral contraction during cold water injection, but for the most part this contraction is recoverable when injection halts. Post-test analysis revealed that fly ash cenosphere pre-treatment had the best effect on improving thermal shock resistance in the cement blends. Thermomechanical modeling of likely stress paths experienced by the cements during heating/cooling cycles shows that elasto-plastic cement constitutive behavior results in most plastic strain occurring during the initial heating steps, with mostly elastic strain occurring during the thermal shock cycles. In conclusion, this agrees with experimental evidence, suggesting that cement damage from shocking occurs via other mechanisms such as chemical alteration, corrosion, and fatigue.

36 MATERIALS SCIENCE↗
NASA NTRS2005

Comprehensive Micromechanics-Analysis Code - Version 4.0

Version 4.0 of the Micromechanics Analysis Code With Generalized Method of Cells (MAC/GMC) has been developed as an improved means of computational simulation of advanced composite materials. The previous version of MAC/GMC was described in "Comprehensive Micromechanics-Analysis Code" (LEW-16870), NASA Tech Briefs, Vol. 24, No. 6 (June 2000), page 38. To recapitulate: MAC/GMC is a computer program that predicts the elastic and inelastic thermomechanical responses of continuous and discontinuous composite materials with arbitrary internal microstructures and reinforcement shapes. The predictive capability of MAC/GMC rests on a model known as the generalized method of cells (GMC) - a continuum-based model of micromechanics that provides closed-form expressions for the macroscopic response of a composite material in terms of the properties, sizes, shapes, and responses of the individual constituents or phases that make up the material. Enhancements in version 4.0 include a capability for modeling thermomechanically and electromagnetically coupled ("smart") materials; a more-accurate (high-fidelity) version of the GMC; a capability to simulate discontinuous plies within a laminate; additional constitutive models of materials; expanded yield-surface-analysis capabilities; and expanded failure-analysis and life-prediction capabilities on both the microscopic and macroscopic scales.

Arnold, S. M.↗
NASA NTRS2013

Constraints on the Composition and Evolution of the Lunar Crust from Meteorite NWA 3163

The lunar meteorite NWA 3163 (paired with NWA 4881, 4483) is a ferroan, feldspathic granulitic breccia characterized by pigeonite, augite, olivine, maskelynite and accessory Tichromite, ilmenite and troilite. Bulk rock geochemical signatures indicate the lack of a KREEP- derived component (Eu/Eu* = 3.47), consistent with previously studied lunar granulites and anorthosites. Bulk rock chondrite-normalized signatures are however distinct from the anorthosites and granulites sampled by Apollo missions and are relatively REE-depleted. In-situ analyses of maskelynite reveal little variation in anorthite content (average An% is 96.9 +/- 1.6, 2 sigma). Olivine is relatively ferroan and exhibits very little variation in forsterite content with mean Fo% of 57.7 +/- 2.0 (2 sigma). The majority of pyroxene is low-Ca pigeonite (En57Fs33Wo10). Augite (En46Fs21Wo33) is less common, comprising approximately 10% of analyzed pyroxene. Two pyroxene thermometry on co-existing orthopyroxene and augite yield an equilibrium temperature of 1070C which is in reasonable agreement with temperatures of 1096C estimated from pigeonite compositions. Rb-Sr isotopic systematics of separated fractions yield an average measured Sr-87/Sr-87 of 0.699282+/-0.000007 (2 sigma). Sr model ages are calculated using a modern day Sr-87/Sr-86 Basaltic Achondrite Best Initial (BABI) value of 0.70475, from an initial BABI value Sr-87/Sr-86 of 0.69891 and a corresponding Rb-87/Sr-97 of 0.08716. The Sr model Thermomechanical analysis (TMA) age, which represents the time of separation of a melt from a source reservoir having chondritic evolution, is 4.56+/-0.1 Ga. A Sr model T(sub RD) age, which is a Rb depletion age and assumes no contribution from Rb in the sample in the calculation, yields 4.34+/-0.1 Ga (i.e. a minimum age). The Ar-Ar dating of paired meteorite NWA 4881 reveals an age of c. 2 Ga, likely representing the last thermal event this meteorite experienced. An older Ar-40/Ar-39 age of c. 3.5 Ga may record the thermal event which produced the granulitic texture. Additional chronological constraints will be provided by Sm-Nd systematics. Ferroan Anorthosites like NWA 3163 have been interpreted to represent direct lunar magma ocean (LMO) crystallization products. If this is the case, trace element concentrations in NWA 3163 primary mineral phases should be in equilibrium with residual LMO liquids present during crystallization of those phases. Results from petrogenetic modeling suggest that the NWA 3163 protolith did not form from crystallization of an initially LREE depleted LMO but rather require an initially chondritic LMO with early garnet crystallization. Furthermore, a two-stage crystallization model where plagioclase crystalized prior to pyroxene (93% vs. 99.5% of LMO crystallization) is implied.

McLeod, C. L.↗
DOE OSTI2020

Single Crystal Plasticity Model with Deformation Twinning for the High Rate Deformation of β-HMX

Deformation twinning is an active deformation mechanism during loading along certain crystallographic directions of β-cyclotetramethylene tetranitramine (β-HMX). In this work, a finite strain thermomechanical model is extended to include twinning as a deformation mechanism in addition to plastic slip. The stress is derived from the free energy expression including a term representing the equation of state. The crystal plasticity framework is used to divide the total strain into inelastic and elastic component, where the elastic part is used in the expression for the free energy. The twin systems are treated as pseudo slip systems and the shear rate on the twin systems is evaluated in terms of the twin resistance and the projection of the stress tensor. The model parameters were calibrated against a set of plate impact experiments performed on β-HMX by Dick et al. The remaining plate impact experiments are used to evaluate the predictive capability of the model. The quality of the model fits and predictions is discussed from a physical and modeling perspective. Particularly, the role of twin modeling on results is highlighted. Finally, the model is used to explore the relationship between the propensity for twinning and crystal orientation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2001

Fabrication and Characterization of SMA Hybrid Composites

Results from an effort to fabrication shape memory alloy hybrid composite (SMAHC) test specimens and characterize the material system are presented in this study. The SMAHC specimens are conventional composite structures with an embedded SMA constituent. The fabrication and characterization work was undertaken to better understand the mechanics of the material system, address fabrication issues cited in the literature, and provide specimens for experimental validation of a recently developed thermomechanical model for SMAHC structures. Processes and hardware developed for fabrication of the SMAHC specimens are described. Fabrication of a SMA14C laminate with quasi-isotropic lamination and ribbon-type Nitinol actuators embedded in the 0' layers is presented. Beam specimens are machined from the laminate and are the focus of recent work, but the processes and hardware are readily extensible to more practical structures. Results of thermomechanical property testing on the composite matrix and Nitinol ribbon are presented. Test results from the Nitinol include stress-strain behavior, modulus versus temperature. and constrained recovery stress versus temperature and thermal cycle. Complex thermomechanical behaviors of the Nitinol and composite matrix are demonstrated, which have significant implications for modeling of SMAHC structures.

Turner, Travis L.↗