Simulation of a highly birefringent photonic crystal fiber in terahertz frequency region