Search NASA⌕ Search

SEARCH · Search NASA

Results for “growth law”

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 73 records · Page 4

Impurity radiation seeding of neoclassical tearing mode growth

Abstract The physics of neoclassical tearing modes (NTMs) is of great concern to tokamak plasma stability and performance, especially in the burning plasma regime. Whereas in many situations the seed events can be clearly identified, such as sawteeth and edge localized modes, the potential seeding mechanism of NTMs due to the resistive tearing instability driven by the impurity radiation cooling still needs more study. Recent NIMROD simulations have demonstrated that local impurity radiation cooling can drive the seed island growth and trigger the subsequent onset of NTM instability. The seed island is mainly driven by the local helical perturbation of the diamagnetic current induced by the perturbed pressure gradient as a result of the impurity radiative cooling on the rational surface. A heuristic closure for the neoclassical viscosity is adopted, and the seed island is further driven by the perturbed bootstrap current induced from the neoclassical electron viscous stress in the extended Ohm’s law. The growth rate of the NTM in simulations is found proportional to the electron neoclassical viscosity, and a theoretical neoclassical driving term is adopted to account for the nonlinear neoclassical island growth in the simulations.

Zeng, Shiyong (ORCID:000000017098818X)↗

Algebraically growing waves in ducts with sheared mean flow

Analysis of the behavior of standing and traveling acoustic waves in a smooth duct with a fluid flow having a sheared mean velocity profile, when the waves grow algebraically as they travel along the duct axis. It is shown that standing waves growing algebraically with the axial distance cannot exist in a smooth duct when the duct wall have a finite resistance. The existence of traveling waves subject to the same law of growth is also dismissed under realistic flow conditions.

Nayfeh, A. H.↗

Coarsening in binary solid-liquid mixtures

A theory of Ostwald ripening has been developed for a solid-liquid mixture cosisting of a low volume fraction array of spherical solid particles in a liquid wherein the coarsening process proceeds via the transport of both heat and mass. It is found that the simultaneous transport of heat and mass during ripening does not alter the exponents of the temporal power laws governing the ripening process from their classical values but does alter the amplitudes of these power laws. The growth rate of the cube of the average particle radius, the rate constant, is found to depend both on the alloy solute concentration and the ratio of the thermal to solutal diffusivities. In most metallic systems, a large decrease in the rate constant can be expected with small additions of solute to a pure metal. Possible extensions of this theory to the analogous problem of ripening in isothermal ternary alloys are also discussed.

Voorhees, P. W.↗

Diffusive crack growth at a bimaterial interface

The diffusional microcrack growth behavior in a bimaterial system is investigated with an aim at estimating service life of advanced ceramic composites under creep-rupture conditions. The crack is assumed to grow via a coupled surface and grain-boundary diffusion under steady state conditions. The tensile stress distribution along the interface ahead of the moving crack tip is solved, and it is found that a new length parameter exists as a scaling factor for which the solution becomes identical to the single phase case when plotted on the nondimensional physical plane. In contrast to the elastic stress solution which shows singularity at the tip, together with oscillatory character away from the tip, the creep stresses have a peak value away from the tip due to a wedging effect and interfacial sliding eliminates stress oscillation resulting in a decoupling between mode I and mode II loading. This solution ties the far-field loading parameter to the crack tip conditions in terms of the unknown crack velocity to give a specific V-K(1) relationship. It is shown that an exponent of 12 in the conventional crack growth power law emerges at the higher applied stress range.

Chuang, Tze-Jer↗

Fatigue crack growth under variable amplitude loading

Fatigue crack growth tests were conducted on an Fe 510 E C-Mn steel and a submerged arc welded joint from the same material under constant, variable, and random loading amplitudes. Paris-Erdogan's crack growth rate law was tested for the evaluation of m and C using the stress intensity factor K, the J-integral, the effective stress intensity factor K(sub eff), and the root mean square stress intensity factor K(sub rms) fracture mechanics concepts. The effect of retardation and residual stresses resulting from welding was also considered. It was found that all concepts gave good life predictions in all cases.

Sidawi, Jihad A.↗

Durability of YSZ Coated Ti2AlC in 1300 °C Mach 0.3 Burner Rig Tests

A thermal barrier coating system survived burner rig testing at 1300 °C for 500 h. A 160 µm thick yttria stabilized zirconia (YSZ) coating was applied to a Ti2AlC MAX phase bar sample by plasma spray physical vapor deposition (PS-PVD) and tested face-on in an atmospheric Mach 0.3 jet fuel burner, using 5-h thermal cycles. No thermal barrier coating (TBC) spallation or recession was observed, only a 2.4 mg/cm2 mass gain. The modest weight gain precluded severe volatility losses under high velocity burner conditions. The coating surface exhibited colonies of (111)flourite fiber-textured columns separated by craze patterns, with no visible moisture attack. The metastable tetragonal t' YSZ phase was obtained initially, transitioning to equilibrium teq and cubic YSZ, but with little detrimental monoclinic. The thickness of the alumina TGO was ~21 to 23 μm under the heated YSZ face and ~13 to 15 μm on the uncoated, cooler backside. The backside exhibited removal of initial transient TiO2 nodules and partial etching of the underlying Al2O3 scale by volatile hydroxides formed in high temperature, high velocity water vapor. Aerodynamic forces produced some bending of the cantilevered sample via creep. The test indicated exceptional stability of YSZ coatings on Ti2AlC under turbine conditions, with thermal expansion matching playing a key role. The purpose of this study was to demonstrate long term durability of YSZ/MAX phase system in aggressive high temperature burner rig testing. MAX phases have been keenly studied because of their unique crystal structure and intriguing properties (Refs. 1 and 2). Having Mn+1(Al,Si)(C,N)n general composition, they are defined as ceramics, but possess unusual desirable attributes such as high conductivity, thermal shock resistance, easy machinability, and deformation tolerance. The mechanical properties derive from weak M-(Al,Si) bonding in the basal plane that leads to sliding and kinking in preference to catastrophic crack growth. Like most ceramics they are phase stable at high temperatures, generally up to 1500 °C. High temperature oxidation resistance is excellent for alumina-forming Ti3AlC2, Ti2AlC, and Cr2AlC, as reviewed by Tallman, et al. (Ref. 3). Compatibility with α-Al2O3 scales is further enhanced in cyclic exposures by a close matching of thermal expansion coefficients, (Ref. 4) i.e., (~9.3, 10.2, 11.3×10–6/K for Al2O3, Ti2AlC, and YSZ, to be discussed). Turbine environments generally contain 10 percent water vapor in the combustion gases, therefore moisture effects can be a concern for some materials (Ref. 5). Furnace tests of MAX phases in high temperature steam generally showed little effect on Al2O3 scale growth (Ref. 6). However, high velocity and high pressure gas can influence scale losses by the formation of volatile reaction products, such as TiO(OH)2 and Al(OH)3 (Refs. 7 to 10). This phenomenon had been discussed for 1100 to 1300 °C high pressure burner rig tests of Ti2AlC (Ref. 11). A single cubic growth rate parameter kcubic was measurably lower than comparable furnace TGA data, but it could be matched reasonably well if corrected for a slight volatility term. In general, a two-parameter cubic-linear growth-volatility law was believed to apply. Corresponding scale volatility loss rates, directly measured at 1300 °C on a pre-oxidized sample, were moderate (0.012 mg/cm2/h) and largely attributed to removal of the initial TiO2 transient scale. A related CH4 burner study of high purity Cr2AlC MAX phase demonstrated 1200 °C durability after 500 rapid (5 min. heat and 2 min. cool) thermal shock cycling (29 h hot time) (Ref. 12). Heating and cooling rates were ~1000 and 500 °C per minute, with a gas velocity of 5 m/s, producing a 75 °C/mm gradient. A 7 μm Al2O3 surface scale and a 13 μm Cr7C3 depletion zone formed with no signs of failure. No evidence of scale volatility was evident, although weight change was not provided, the velocity was moderate, and the total hot time was not extensive. The same high gradient BRT was used to produce 1400 °C surface temperatures for a YSZ/Cr2AlC/IN738 system in the first study of MAX phases used as bond coats for thermal barrier coatings (TBC) (Ref. 25). Here TBC failure was reported after 745 cycles, with only a 1.5 μm Al2O3 scale entrained within a porous, Cr7C3 bondcoat depletion phase. YSZ thermal barrier coatings have been considered to be a compatible complement to Al-MAX phases because of thermal expansion matching and extremely low volatility in water vapor. Initial studies showed superior oxidative stability up to 1300 °C, for long times (at least 500 h) for Ti2AlC substrates and less (268 h) for Cr2AlC, while withstanding large alumina TGO scale thickness (~35 to 40 μm) (Refs. 13 and 14). By comparison, typical superalloy systems can only survive 1150 °C maximum interface temperatures for extended periods, with a maximum sustained TGO below 10 μm (Ref. 15). High temperature SiC based systems are known to form slow-growing SiO2 scales. But these are subject to rate enhancement and volatile Si(OH)4 products in the presence of water vapor, as described comprehensively by Opila, et al. (Refs. 5, 16 to 19). Net weight losses are generally observed in high velocity, high pressure burner rig studies (e.g., 0.084 mg/cm2/h at 1300 °C) (Ref. 20). Furthermore, the loss rates have been shown from chemical physics to scale with v1/2 and pH2O 2 (Ref. 16). Low activity, moisture-resistant environmental barrier coatings (EBC), such as rare earth silicates, are needed to prevent substrate recession under turbine conditions (Refs. 21 to 23).

Smialek, James L.↗

Scaling laws for two-dimensional dendritic crystal growth in a narrow channel

Here, we investigate analytically and computationally the dynamics of two-dimensional needle crystal growth from the melt in a narrow channel. Our analytical theory predicts that, in the low supersaturation limit, the growth velocity $\textit{V}$ decreases in time $\textit{t}$ as a power law $V ~t^{–2/3}$, which we validate by phase-field and dendritic-needle-network simulations. Simulations further reveal that, above a critical channel width $Λ ≈ 5l_D$, where $l_D$ is the diffusion length, needle crystals grow with a constant $V < V_s$, where $V_s$ is the free-growth needle crystal velocity, and approaches $V_s$ in the limit $Λ \gg l_D$.

36 MATERIALS SCIENCE↗

Slow Crack Growth and Fracture Toughness of Sapphire for the International Space Station Fluids and Combustion Facility

The fracture toughness, inert flexural strength, and slow crack growth parameters of the r- and a-planes of sapphire grown by the Heat Exchange Method were measured to qualify sapphire for structural use in the International Space Station. The fracture toughness in dry nitrogen, K(sub Ipb), was 2.31 +/- 0.12 MPa(square root of)m and 2.47 +/- 0.15 MPa(squre root of)m for the a- and r-planes, respectively. Fracture toughness measured in water via the operational procedure in ASTM C1421 was significantly lower, K(sub Ivb) = 1.95+/- 0.03 MPa(square root of)m, 1.94 +/- 0.07 and 1.77 +/- 0.13 MPa(square root of)m for the a- , m- and r-planes, respectively. The mean inert flexural strength in dry nitrogen was 1085 +/- 127 MPa for the r-plane and 1255 +/- 547 MPa for the a-plane. The power law slow crack growth exponent for testing in water was n = 21 +/- 4 for the r-plane and n (greater than or equal to) 31 for the a-plane. The power law slow crack growth coefficient was A = 2.81 x 10(exp -14) m/s x (MPa(squre root of)m)/n for the r-plane and A (approx. equals)2.06 x 10(exp -15) m/s x (MPa(square root of)m)/n for the a-plane. The r- and a-planes of sapphire are relatively susceptible to stress corrosion induced slow crack growth in water. However, failure occurs by competing modes of slow crack growth at long failure times and twinning for short failure time and inert environments. Slow crack growth testing needs to be performed at low failure stress levels and long failure times so that twinning does not affect the results. Some difficulty was encountered in measuring the slow crack growth parameters for the a-plane due to a short finish (i.e., insufficient material removal for elimination of the damage generated in the early grinding stages). A consistent preparation method that increases the Weibull modulus of sapphire test specimens and components is needed. This would impart higher component reliability, even if higher Weibull modulus is gained at the sacrifice of absolute strength of the component. The current specification frequently used for the preparation of sapphire test specimens and components (e.g., a "60/40" scratch-dig finish) is inadequate to avoid a short finish.

Salem, Jonathan A.↗

The Transition from Complex Crater to Peak-Ring Basin on the Moon: New Observations from the Lunar Orbiter Laser Altimeter (LOLA) Instrument

Impact craters on planetary bodies transition with increasing size from simple, to complex, to peak-ring basins and finally to multi-ring basins. Important to understanding the relationship between complex craters with central peaks and multi-ring basins is the analysis of protobasins (exhibiting a rim crest and interior ring plus a central peak) and peak-ring basins (exhibiting a rim crest and an interior ring). New data have permitted improved portrayal and classification of these transitional features on the Moon. We used new 128 pixel/degree gridded topographic data from the Lunar Orbiter Laser Altimeter (LOLA) instrument onboard the Lunar Reconnaissance Orbiter, combined with image mosaics, to conduct a survey of craters >50 km in diameter on the Moon and to update the existing catalogs of lunar peak-ring basins and protobasins. Our updated catalog includes 17 peak-ring basins (rim-crest diameters range from 207 km to 582 km, geometric mean = 343 km) and 3 protobasins (137-170 km, geometric mean = 157 km). Several basins inferred to be multi-ring basins in prior studies (Apollo, Moscoviense, Grimaldi, Freundlich-Sharonov, Coulomb-Sarton, and Korolev) are now classified as peak-ring basins due to their similarities with lunar peak-ring basin morphologies and absence of definitive topographic ring structures greater than two in number. We also include in our catalog 23 craters exhibiting small ring-like clusters of peaks (50-205 km, geometric mean = 81 km); one (Humboldt) exhibits a rim-crest diameter and an interior morphology that may be uniquely transitional to the process of forming peak rings. Comparisons of the predictions of models for the formation of peak-ring basins with the characteristics of the new basin catalog for the Moon suggest that formation and modification of an interior melt cavity and nonlinear scaling of impact melt volume with crater diameter provide important controls on the development of peak rings. In particular, a power-law model of growth of an interior melt cavity with increasing crater diameter is consistent with power-law fits to the peak-ring basin data for the Moon and Mercury. We suggest that the relationship between the depth of melting and depth of the transient cavity offers a plausible control on the onset diameter and subsequent development of peak-ring basins and also multi-ring basins, which is consistent with both planetary gravitational acceleration and mean impact velocity being important in determining the onset of basin morphological forms on the terrestrial planets.

Baker, David M. H.↗

CARES/Life Ceramics Durability Evaluation Software Enhanced for Cyclic Fatigue

The CARES/Life computer program predicts the probability of a monolithic ceramic component's failure as a function of time in service. The program has many features and options for materials evaluation and component design. It couples commercial finite element programs--which resolve a component's temperature and stress distribution--to reliability evaluation and fracture mechanics routines for modeling strength-limiting defects. The capability, flexibility, and uniqueness of CARES/Life have attracted many users representing a broad range of interests and has resulted in numerous awards for technological achievements and technology transfer. Recent work with CARES/Life was directed at enhancing the program s capabilities with regards to cyclic fatigue. Only in the last few years have ceramics been recognized to be susceptible to enhanced degradation from cyclic loading. To account for cyclic loads, researchers at the NASA Lewis Research Center developed a crack growth model that combines the Power Law (time-dependent) and the Walker Law (cycle-dependent) crack growth models. This combined model has the characteristics of Power Law behavior (decreased damage) at high R ratios (minimum load/maximum load) and of Walker law behavior (increased damage) at low R ratios. In addition, a parameter estimation methodology for constant-amplitude, steady-state cyclic fatigue experiments was developed using nonlinear least squares and a modified Levenberg-Marquardt algorithm. This methodology is used to give best estimates of parameter values from cyclic fatigue specimen rupture data (usually tensile or flexure bar specimens) for a relatively small number of specimens. Methodology to account for runout data (unfailed specimens over the duration of the experiment) was also included.

Nemeth, Noel N.↗

The asymptotic state of rotating homogeneous turbulence at high Reynolds numbers

The long-time, asymptotic state of rotating homogeneous turbulence at high Reynolds numbers has been examined using large-eddy simulation of the incompressible Navier-Stokes equations. The simulations were carried out using 128 x 128 x 512 collocation points in a computational domain that is four times longer along the rotation axis than in the other directions. Subgrid-scale motions in the simulations were parameterized using a spectral eddy viscosity modified for system rotation. Simulation results show that in the asymptotic state the turbulence kinetic energy undergoes a power-law decay with an exponent which is independent of rotation rate, depending only on the low-wavenumber form of the initial energy spectrum. Integral lengthscale growth in the simulations is also characterized by power-law growth; the correlation length of transverse velocities exhibiting much more rapid growth than observed in non-rotating turbulence.

Squires, Kyle D.↗

A Novel Method for Characterizing Fatigue Delamination Growth Under Mode I Using the Double Cantilever Beam Specimen

A novel method is proposed to obtain Mode I delamination growth rate from a Double Cantilever Beam (DCB) specimen. In the proposed method, Unidirectional (UD) DCB specimens are tested in fatigue at different initial maximum energy release rates levels. The growth rate data obtained in the first increments of crack growth at each maximum energy release rate level are used to generate a Paris Law equation, which characterizes delamination growth rate without fiber-bridging, and can also be used to determine a delamination onset curve. The remaining delamination growth rate data from each test are used to determine a modified Paris law, which characterizes the delamination growth rate in a DCB specimen, explicitly accounting for fiber-bridging. The proposed expression captures well the scatter in experimental data obtained using the DCB specimens, suggesting its adequacy. The Paris Law characterizing delamination growth rate without fiber-bridging predicts higher delamination growth rates for the same maximum energy release rate applied, leading to a conservative estimate for delamination growth. This is particularly relevant, since in generic ply interfaces, fiber-bridging is less predominant than in UD DCB specimens. Failing to account for fiber-bridging in UD DCB specimens may underestimate the delamination growth rate, yielding non-conservative predictions.

Carvalho, Nelson↗

Power Law Versus Exponential Form of Slow Crack Growth of Advanced Structural Ceramics: Dynamic Fatigue

The life prediction analysis based on an exponential crack velocity formulation was examined using a variety of experimental data on glass and advanced structural ceramics in constant stress-rate ("dynamic fatigue") and preload testing at ambient and elevated temperatures. The data fit to the strength versus In (stress rate) relation was found to be very reasonable for most of the materials. It was also found that preloading technique was equally applicable for the case of slow crack growth (SCG) parameter n > 30. The major limitation in the exponential crack velocity formulation, however, was that an inert strength of a material must be known priori to evaluate the important SCG parameter n, a significant drawback as compared to the conventional power-law crack velocity formulation.

Choi, Sung R.↗

Dendritic Alloy Solidification Experiment (DASE)

A space experiment, and supporting ground-based research, is proposed to study the microstructural evolution in free dendritic growth from a supercooled melt of the transparent model alloy succinonitrile-acetone (SCN-ACE). The research is relevant to equiaxed solidification of metal alloy castings. The microgravity experiment will establish a benchmark for testing of equiaxed dendritic growth theories, scaling laws, and models in the presence of purely diffusive, coupled heat and solute transport, without the complicating influences of melt convection. The specific objectives are to: determine the selection of the dendrite tip operating state, i.e. the growth velocity and tip radius, for free dendritic growth of succinonitrile-acetone alloys; determine the growth morphology and sidebranching behavior for freely grown alloy dendrites; determine the effects of the thermal/solutal interactions in the growth of an assemblage of equiaxed alloy crystals; determine the effects of melt convection on the free growth of alloy dendrites; measure the surface tension anisotropy strength of succinon itrile -acetone alloys establish a theoretical and modeling framework for the experiments. Microgravity experiments on equiaxed dendritic growth of alloy dendrites have not been performed in the past. The proposed experiment builds on the Isothermal Dendritic Growth Experiment (IDGE) of Glicksman and coworkers, which focused on the steady growth of a single crystal from pure supercooled melts (succinonitrile and pivalic acid). It also extends the Equiaxed Dendritic Solidification Experiment (EDSE) of the present investigators, which is concerned with the interactions and transients arising in the growth of an assemblage of equiaxed crystals (succinonitrile). However, these experiments with pure substances are not able to address the issues related to coupled heat and solute transport in growth of alloy dendrites.

Beckermann, C.↗

Conceptualizing Biogeochemical Reactions With an Ohm's Law Analogy

Abstract In studying problems like plant‐soil‐microbe interactions in environmental biogeochemistry and ecology, one usually has to quantify and model how substrates control the growth of, and interaction among, biological organisms (and abiotic factors, e.g., adsorptive mineral soil surfaces). To address these substrate‐consumer relationships, many substrate kinetics and growth rules have been developed, including the famous Monod kinetics for single‐substrate‐based growth and Liebig's law of the minimum for multiple‐nutrient‐colimited growth. However, the mechanistic basis that leads to these various concepts and mathematical formulations and the implications of their parameters are often quite uncertain. Here, we show that an analogy based on Ohm's law in electric circuit theory is able to unify many of these different concepts and mathematical formulations. In this Ohm's law analogy, a resistor is defined by a combination of consumers’ and substrates’ kinetic traits. In particular, the resistance is equal to the mean first passage time that has been used to derive the Michaelis‐Menten kinetics under substrate replete conditions for a single substrate as well as the predation rate of individual organisms. We further show that this analogy leads to important insights on various biogeochemical problems, such as (a) multiple‐nutrient‐colimited biological growth, (b) denitrification, (c) fermentation under aerobic conditions, (d) metabolic temperature sensitivity, and (e) the legitimacy of Monod kinetics for describing bacterial growth. We expect that our approach will help both modelers and nonmodelers to better understand and formulate hypotheses when studying certain aspects of environmental biogeochemistry and ecology.

54 ENVIRONMENTAL SCIENCES↗

Dependency of Shear Strength on Test Rate in SiC/BSAS Ceramic Matrix Composite at Elevated Temperature

Both interlaminar and in-plane shear strengths of a unidirectional Hi-Nicalon(TM) fiber-reinforced barium strontium aluminosilicate (SiC/BSAS) composite were determined at 1100 C in air as a function of test rate using double notch shear test specimens. The composite exhibited a significant effect of test rate on shear strength, regardless of orientation which was either in interlaminar or in in-plane direction, resulting in an appreciable shear-strength degradation of about 50 percent as test rate decreased from 3.3 10(exp -1) mm/s to 3.3 10(exp -5) mm/s. The rate dependency of composite's shear strength was very similar to that of ultimate tensile strength at 1100 C observed in a similar composite (2-D SiC/BSAS) in which tensile strength decreased by about 60 percent when test rate varied from the highest (5 MPa/s) to the lowest (0.005 MPa/s). A phenomenological, power-law slow crack growth formulation was proposed and formulated to account for the rate dependency of shear strength of the composite.

Choi, Sung R.↗

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↗