A novel mathematical model has been developed to understand the mechanism of blister initiation and propagation. The model employs a two-part theoretical approach encompassing the debondment of a coating film from the substrate, coupled with the design components incorporating diffusion and fracture mechanics, where the latter is derived from equi-biaxial tensile loading. Integrating the two components, a comprehensive mathematical design for the propagation of blister boundaries based on specific toughness functions and mode adjustment parameters has been developed. This approach provided a reliable and efficient prediction method for blister growth rate and mechanisms. The model provided a foundation for holistic design based on diffusion and mechanic components to enable better understanding of the debondment of thin elastic films bonded to a metallic substrate.