Varicella-zoster virus(VZV) is known to be an unstable virus. Infectious cells infect directly adjacent cells and VZV has the strongly cell-associated nature. Kinetics analysis of VZV proliferation shows different experimental values due to the environmental sensitivity of VZV. So the prediction of exact yield of VZV is difficult. In order to solve this problem, a theoretical model is required to estimate more precise yield of VZV. The yield of VZV is related to multiplicity of infection(MOI) and time of culture in the proliferation of VZV. A model was developed by using the finite element method considering the cell-to-cell infection characteristic of VZV. The spreading time of plaque forming cells to the neighboring cells was determined as 8 hours by several experiments. It was assumed that plaque forming cells lost infectivity within 48 hours after it was infected. From the simulation of two-dimensional model used in this work, it was confirmed that the maximum yield of VZV could be obtained at the MOI of 1:10, which corresponded with experimental results. The VZV propagation model can be extended to the future 3-dimensional model which is closer to real VZV propagation.