Finite element simulations based on an interface cohesive zone model (CZM) have been developed to mimic the interfacial cracking behavior between the α-Al2O3 thermally grown oxide (TGO) and the aluminum rich Pt-Al metallic bond coat (BC) during cooling from high temperature to ambient temperature. A two dimensional half-periodic sinusoidal geometry corresponding to interface undulation is modelled. The effects of TGO thickness and interface asperity on the stress distribution and the cracking behavior are examined by parametric studies. The simulation results show that cracking behavior due to residual stress and interface asperity during cooling process leads to stress redistribution around the rough interface. The TGO thickness has strong influence on the maximum tensile stress of TGO and the interfacial crack development. For the sinusoidal asperities, there exist a critical amplitude above which interfacial cracking is energetically favored. For any specific TGO thickness, crack initiation is dominated by the amplitude while crack propagation is restricted to the combine actions of the wavelength and the amplitude of the sinusoidal asperity.