  We extend the doubly degenerate Cahn-Hilliard (DDCH) models for isotropic
surface diffusion, which yield more accurate approximations than classical
degenerate Cahn-Hilliard (DCH) models, to the anisotropic case. We consider
both weak and strong anisotropies and demonstrate the capabilities of the
approach for these cases numerically. The proposed model provides a variational
and energy dissipative approach for anisotropic surface diffusion, enabling
large scale simulations with material-specific parameters.
