  We combine the total variation flow suitable for crystal modeling and image
analysis with the dynamic boundary conditions. We analyze the behavior of
facets at the parts of the boundary where these conditions are imposed. We
devote particular attention to the radially symmetric data. We observe that the
boundary layer detachment actually can happen at concave parts of the boundary
