Thermal resistance across the interface between touching surfaces is critical for many industrial applications. We developed a network model to predict the macroscopic thermal resistance of mechanically contacting surfaces. Contacting interfaces are fractally rough, with small islands of locally intimate contact separated by regions with a wider gas filled boundary gap. Heat flow across the interface is therefore heterogeneous and thus the contact model is based on a network of thermal resistors representing boundary resistance at local contacts and the access resistance for lateral transport to contacts. Molecular dynamics simulations have been performed to characterize boundary resistance of Silicon Alumina interfaces for testing the sensitivity of thermal resistance to contact opening. Boltzmann transport simulations of access resistance in Si are conducted in the ballistic transport regime.