In this paper, we report our numerical study on a novel poly-rubbing process which employs more than one rubbing condition on the same boundary surface of liquid crystal for implementing a multi-domain. A compact surface anchoring model was implemented in a commercially available 3D finite element method (FEM) numerical solver, TechWiz LCD, for calculating the director distribution of molecules and optical characteristics of a pixel structure. Our surface anchoring model was tested by investigating the director distribution of the Vertical Alignment (VA) and Electrically Controlled Birefringence (ECB) mode cells with poly-rubbed surface. The size of each pixel was 88 x 264 mu m(2) while the cell gap is mu m. Our numerical simulation implied that LC mode with poly-rubbed surface can provide a viewing angle as good as that of the conventional PVA cell having a complicated electrode structure.