This paper presents a three-dimensional particle-in-cell (PIC) simulation of a Ka-band relativistic Cherenkov source with a slow wave structure (SWS) consisting of metal photonic band gap (PBG) structures. In the simulation, a perfect match layer boundary is employed to absorb passing band modes supported by the PBG lattice with an artificial metal boundary. The simulated axial field distributions in the cross section and surface of the SWS demonstrate that the device operates in the vicinity of the π point of a TM01-1ike mode. The Fourier transformation spectra of the axial fields as functions of time and space show that only a single frequency appears at 36.27 GHz, which is in good agreement with that of the intersection of the dispersion curve with the slow space charge wave generated on the beam. The simulation results demonstrate that the SWS has good mode selectivity.