Weld lines are one of the typical quality issues of injection molded parts. They are unavoidable even for moderate complex products. To explore the formation of weld lines in injection molding, a flow model based on viscoelastic constitutive equation is presented in this article to predict and evaluate of weld lines location and properties, in which Rolie-Poly constitutive equation is used to describe the viscoelastic properties of polymer melts and level set method is used to capture melt front. Furthermore, the formation and evolution processes of weld lines are investigated by co-injection molding technique in simulation, since the processes are determined by the flow front and rheological history of the melt streams. Numerical results show that weld lines can be predicted accurately, and the formation and evolution processes of weld line can be shown more clearly by co-injection molding technique. In addition, properties of polymer melt such as velocity, pressure, principal stress difference, and molecular orientation in the weld line regions are investigated and presented.