An efficient mesh model repair method is proposed in this paper, which can repair the holes with C~(k-1) boundary continuity. It is mainly composed of the following steps: Firstly, the hole boundary is triangulated to generate the initial surface patch, of which the edges are distributed along the boundary averagely similar to an curtain covering at a window. Secondly, the initial surface patch is refined, and the interior edges are relaxed while refining to get a Delaunay-like refinement surface patch, which approximates the density of the surrounding mesh. Finally, the refinement surface patch is reshaped based on the minimization of the surface energy solved by the k order discrete Euler-Lagrange equation in order to realize the C~(k-1) continuity repair. The proposed method is proved to be efficient and robust by experiments of different kinds of boundary conditions.