Search NASA⌕ Search

Engineering topics

Wen, You-Hai

Publications and source records attributed to Wen, You-Hai.

Phase-field modeling of thermally-grown oxide and damage evolution in environmental barrier coatings

Silicon carbide-based ceramic matrix composites protected by environmental barrier coatings (EBCs) present a promising materials solution for next-generation gas turbines. Developing more robust and efficient EBCs is therefore of significant technological importance. During the service in high-temperature oxidative environments, there is a thermally grown oxide (TGO) layer, spontaneously formed in the EBC system. TGO is recognized as a critical factor for the degradation and failure of EBCs, yet the detailed mechanisms of TGO growth and its effect on EBC failure remain unclear. In this study we develop a comprehensive chemo-mechano-phase-field model to simulate growth of the TGO in EBCs, factoring in creep and deformation, and especially the cracking behaviors. The volume expansion due to TGO growth and the resulting large inelastic deformation are addressed by using our recently developed, so-called incremental realization of inelastic deformation (IRID) algorithm, in combination with an adapted Hu-Chen spectral solver for elasticity. Simulations of TGO growth are performed considering different growth modes of TGOs determined mainly by the ratio of oxidant permeability in the topcoat to that in the TGO itself. Large-scale three-dimensional (3D) simulations are performed to model the formation of interconnecting vertical/channel cracks (often called ‘mud cracks’). The simulated crack morphology are in excellent agreement with the experimental observations from the literature. The simulations also provide insights into the cracking of EBCs and its dependence on the structure and constituent properties of the coating system. Furthermore, these results demonstrate the developed damage model can be a useful tool for design of more durable EBCs.

Cracking↗

Phase-field modeling of crack growth under coupled creep-fatigue

Crack growth under coupled creep-fatigue strongly influences service life of metallic components operating at high temperatures. However, there lack numerical models that enable direct simulations of crack growth and plasticity development under coupled creep-fatigue. Here, in this study, a unified phase-field model is developed for simulating crack growth under fatigue, creep, and coupled creep-fatigue. This model is able to reproduce the Paris laws for both brittle and ductile materials in pure fatigue mode, which shows that the degradation rate due to cyclic stress in ductile materials is generally larger than that in brittle materials. Fatigue simulations with different yield strengths demonstrate that plasticity facilitates crack growth in typical ductile materials. High-throughput simulations are performed under coupled creep-fatigue conditions, with varying hold-stress levels and hold-time. By analyzing crack growth rate under different hold-time, a creep-fatigue interaction term is obtained. Simulation results demonstrate that the creep-fatigue interaction is caused by smaller stress gradient and enhanced degradation rate ahead of the crack tips under coupled creep-fatigue, compared to the situation of pure fatigue. This work reveals the origin of different Paris-law exponents in brittle and ductile materials and demonstrates how the interplay between creep and fatigue affects crack growth.

36 MATERIALS SCIENCE↗

Phase-field modeling of alloy oxidation at high temperatures

Oxide growth is a complex process involving transport of reactive species, heterogeneous reactions, and microstructure evolution. Predicting oxidation kinetics and especially the oxide morphological change has been a long-standing challenge. Here we develop a phase-field model for predicting the oxide growth kinetics of a multicomponent alloy during high temperature oxidation, focusing on internal oxidation (non-protective) and its transition to external oxidation (protective). In this work, the predicted kinetics and oxide morphology are analyzed and compared to the classical Wagner’s theory and an existing analytical model by Zhao and Gleeson. Some assumptions used in the analytical models and the limitation are discussed. In addition, it is demonstrated that the morphology and distribution of the initial oxide nuclei play an important role in the later stage oxide connectivity and thus the transition to external oxidation.

36 MATERIALS SCIENCE↗

Phase-field modeling of crack growth and mitigation in solid oxide cells

Fracture and crack growth is one of the main degradation mechanisms in solid oxide cells (SOCs). However, the modeling of crack growth in SOCs is challenging due to their complex microstructures and possible plasticity development within the Ni particles in Ni-based SOC electrodes. In this study, a phase-field fracture model is developed, which incorporates the SOC microstructures and phase-dependent material properties, including yield strength, fracture toughness in the bulk and at the interphase boundaries. The model is employed to study crack initiation and growth under thermal and redox cycling on the hydrogen electrode side of SOCs. The simulation results demonstrate that under thermal cycling, work-zone cracking dominates in electrolyte-supported SOCs with cracks initiated at the triple-phase boundaries, while only minor mechanical degradation occurs in hydrogen-electrode-supported SOCs after hundreds of thermal cycles. Under redox cycling, through-cracking of yttria-stabilized zirconia (YSZ) in the hydrogen electrode and electrolyte layers dominates. Finally, the simulation results suggest several crack-mitigation strategies, including decreasing the porosity in the hydrogen electrode support layer and synchronizing thermal strain to balance oxidation strain.

08 HYDROGEN↗

Modeling Ni redistribution in the hydrogen electrode of solid oxide cells through Ni(OH) 2 diffusion and Ni-YSZ wettability change

Ni redistribution in the hydrogen electrodes of solid oxide cells is an important degradation mechanism. Its driving force is still under debate. This work focuses on the interplay between Ni(OH) 2 diffusion and Ni-YSZ wettability change. With a Ni(OH) 2 diffusion model developed in our previous work, we further employ three models for the Ni-YSZ contact angle to quantify the Ni-YSZ wettability change. The microstructure evolutions in a reconstructed Ni-YSZ electrode are simulated under selected experimental operating conditions. It is shown that the phenomenological model can capture the driving force of Ni spreading/detachment on YSZ surface and the Ni migration reported in experiments can be qualitatively reproduced through the competition between Ni(OH) 2 diffusion and Ni-YSZ wettability change. It is also found that the initial microstructure of the active layer and the microstructure characteristics in the support layer can strongly affect the distribution of steam partial pressure and overpotential, hence the Ni redistribution, in the active layer, which may explain the inconsistencies in experiments. Further, our results show that there are two missing pieces in current theory of Ni redistribution: a mechanism of fast Ni diffusion under humid atmosphere and the physical reason behind the Ni spreading/detachment.

25 ENERGY STORAGE↗

A Phase-Field Study on Internal to External Oxidation Transition in High-Temperature Structural Alloys

Structural alloys applied at high temperatures rely on an external dense layer of oxide scale for protection. As some reactive alloy components are selectively oxidized internally, understanding how these dispersed metal oxide particles reach the surface to form a protective scale—the so-called internal to external oxidation transition—is crucial for designing these alloys. While the literature is replete with experimental studies on oxidation of alloys, there is a lack of computational studies in this realm due to the complex nature of coupled reaction and diffusion processes in multicomponent multi-phase alloy systems. In this work, we apply a recently developed phase-field model to simulate the oxidation processes under different compositions and nucleation scenarios to gain insights into how a continuous oxide scale can be established. The results show that while alloy composition is critical for internal to external oxidation transition, the oxide nuclei size, shape and distribution also have significant impact on the transition kinetics.

36 MATERIALS SCIENCE↗

Phase-field framework with constraints and its applications to ductile fracture in polycrystals and fatigue

Abstract Modeling of ductile fracture in polycrystalline structures is challenging, since it requires integrated modeling of cracks, crystal plasticity, and grains. Here we extend the typical phase-field framework to the situations with constraints on the order parameters, and formulate two types of phase-field models on ductile fracture. The Type-I model incorporates three sets of order parameters, which describe the distributions of cracks, plastic strain, and grains, respectively. Crystal plasticity is employed within grain interiors accommodated by J 2 plasticity at grain boundaries. The applications of the Type-I model to single crystals and bicrystals demonstrate the influences of grain orientations and grain boundaries on crack growth. In the Type-II model, J 2 plasticity is assumed for the whole system and grain structures are neglected. Taking advantage of the efficiency of the fast Fourier transform, our Type-II model is employed to study low cycle fatigue. Crack closure and striation-like patterning of plastic strain are observed in the simulations. Crack growth rate is analyzed as a function of the J -integral, and the simulated fatigue life as a function of plastic strain agrees with the Coffin–Manson relation without a priori assumption.

36 MATERIALS SCIENCE↗

Phase-field model of precipitation processes with coherency loss

A phase-field model is proposed to simulate coherency loss coupled with microstructure evolution. A special field variable is employed to describe the degree of coherency loss of each particle and its evolution is governed by a Ginzburg-Landau type kinetic equation. For the sake of computational efficiency, a flood-fill algorithm is introduced that can drastically reduce the required number of field variables, which allows the model to efficiently simulate a large number of particles sufficient for characterizing their statistical features during Ostwald ripening. The model can incorporate size dependence of coherency loss, metastability of coherent particles, and effectively incorporate the underlying mechanisms of coherency loss by introducing a so-called differential energy criterion. The model is applied to simulate coarsening of Al 3 Sc precipitates in aluminum alloy and comprehensively compared with experiments. Our results clearly show how the particle size distribution is changed during coherency loss and affects the coarsening rate.

36 MATERIALS SCIENCE↗