In this thesis, we study the problem on how to find an equilibrium configuration of a knotted Kirchhoff elastic rod by evolving it within the same knot class. The evolution is given by the $L^2$ gradient flow of the total energy composed of the Kirchhoff elastic energy and the Möbius knot energy. We would prove the long-time existence of smooth solutions for the gradient flow and discuss the asymptotic limit.