Viscosity-based numerical model for fault-zone development in drape folding