In this paper, we introduce a high-order accurate constrained transport type finite volume method to solve ideal magnetohydrodynamic equations on two-dimensional triangular meshes. A new divergence-free WENO-based reconstruction method is developed to maintain exactly divergence-free evolution of the numerical magnetic field. In this formulation, the normal component of the magnetic field at each face of a triangle is reconstructed uniquely and with the desired order of accuracy. Additionally, a new weighted flux interpolation approach is also developed to compute the z-component of the electric field at vertices of grid cells. We also present numerical examples to demonstrate the accuracy and robustness of the proposed scheme.
