The physics-informed neural networks (PINNs) have shown great potential in solving a variety of high-dimensional partial differential equations (PDEs), but the complexity of a realistic problem still restricts the practical application of the PINNs for solving most complicated PDEs. In this paper, we propose a parallel framework for PINNs that is capable of modeling two-phase flows with complicated interface evolution. The proposed framework divides the problem into several simplified subproblems and solves them through training several PINNs on corresponding subdomains simultaneously. To enhance the accuracy of the parallel training framework in two-phase flow, the overlapping domain decomposition method is adopted. The optimal subnetwork sizes and partitioned method are systematically discussed, and a series of cases including a bubble rising, droplet splashing, and the Rayleigh–Taylor instability are applied for quantitative validation. The maximum relative error of quantitative values in these cases is 0.1319. Our results show that the proposed framework not only can accelerate the training procedure of PINNs, but also can capture the spatiotemporal evolution of the interface between various phases. This framework overcomes the difficulties of training PINNs to solve a forward problem in two-phase flow, and it is expected to model more realistic dynamic systems in nature.