A new type of cluster algorithm is developed to simulate the two-dimensional spin one half quantum Heisenberg antiferromagnet (QHAF) on a hexagonal lattice. The numerical results show that the method can be applied to study efficiently the QHAF systems at lowertemperatures, with large Trotter number, on the general bipartite lattices.