A computational analysis of a ZrO2–SiO2 scale for a ZrB2–ZrC–Zr ultrahigh temperature ceramic composite system