Power flow models of Interline Power Flow Controllers (IPFC) for large-scale power systems are studied, in details. Mathematical models of the IPFC, using the d-q axis decompositions of control parameters are derived. In this framework, for each IPFC, only two control parameters are added to the unknown vector in the iteration formula and the quadratic convergence characteristic is preserved. Simulations results from several practical large-scale power systems embedded with multiple Convertible Static Compensators (CSCs) demonstrate the effectiveness of the proposed models. Comparisons with existing models are made to elucidate the performance of the convergence.