This thesis presents the consistency, stability and convergence analysis of an operator splitting method, namely the iterative operator splitting method, using various approaches for parabolic partial differential equations. The idea of the method is based first on splitting the complex problems into simpler equations. Then, each sub-problem is combined with iterative schemes and efficiently solved with suitable integrators. The analyses are based on the type of the operators of the problems. When the operators are bounded, the consistency is proved in two ways: first from derived explicit local error bounds and the second using the Taylor series expansion after combining iterative schemes with midpoint rule. As for the unbounded operators, since the Taylor series expansion is no longer valid, the consistency is derived using C0 semigroup theory. The stability is presented by constructing stability functions for each iterative schemes when the operators are bounded. For the unbounded, two stability analyses are offered: first one uses the continuous Fourier transform and the second uses semigroup theory. Lax- Richtmyer equivalence theorem and Lady Windermere’s fan argument which combine the stability and consistency are proposed for the convergence. In the computational part, the method is applied to three linear parabolic PDEs and to Korteweg-de Vries equation. These three equations are capillary formation model in tumor angiogenesis, solute transport problem and heat equation. Finally, numerical results are presented to illustrate the high accuracy and efficiency of the method relative to other classical methods. These numerical results align with the obtained theoretical results.