  In this paper, we construct and analyze a multiscale (finite element) method
for parabolic problems with heterogeneous dynamic boundary conditions. As
origin, we consider a reformulation of the system in order to decouple the
discretization of bulk and surface dynamics. This allows us to combine
multiscale methods on the boundary with standard Lagrangian schemes in the
interior. We prove convergence and quantify explicit rates for low-regularity
solutions, independent of the oscillatory behavior of the heterogeneities. As a
result, coarse discretization parameters, which do not resolve the fine scales,
can be considered. The theoretical findings are justified by a number of
numerical experiments including dynamic boundary conditions with random
diffusion coefficients.
