In this paper an accurate analysis of two-dimensional (2D) Photonic Crystal (PhC) based multimode resonant cavities is carried out. The analysis is performed with a robust and accurate Finite Volume Time Domain (FVTD) technique. The analysis proves the ability of the FVTD method to extract different resonant modes from a multimode PhC resonant cavity with the use of appropriate source profiles. A detailed explanation on how the source is engineered and used to excite different modes is given. Furthermore, parameters such as resonant frequency and quality factor for each resonant mode are accurately calculated.