3D fluid based numerical modelling is done for a VHF multi hollow cathode array plasma enhanced chemical vapor deposition system. In order to understand the fundamental characteristics of it, Ar plasma is analyzed with a condition of 40 MHz, 100 Vrf and 1 Torr. For hole array of 6 mm diameter and 20 mm inter-hole distance, plasma is well confined within the hole at an electrode gap of 10 mm. The peak plasma density was $5{ imes}10^{11}#/cm^3$ at the center of the hole. When the substrate was assumed at ground potential, electron temperature showed a peak at the vicinity of the grounded walls including the substrate and chamber walls. The reaction rate of metastable based two step ionization was 10 times higher than the direct electron impact ionization at this condition. For $H_2$, the spatial localization of discharge is harder to get than Ar due to various pathways of electron impact reactions other than ionization.