Modelings of the interface distribution and flow-induced residual stresses and birefringence in the sequential co-injection molding (CIM) of a center-gated disk were carried out using a numerical scheme based on a hybrid finite element/finite difference/control volume method. A nonlinear viscoelastic constitutive equation and stress-optical rule were used to model the frozen-in flow stresses in disks. The compressibility of melts is included in modeling of the packing and cooling stages and not in the filling stage. The thermally induced residual birefringence was calculated using the linear viscoelastic and photoviscoelastic constitutive equations combined with the first-order rate equation for volume relaxation and the master curves for the relaxation modulus and strain-optical coefficient functions of each polymer. The influence of the processing variables including melt and mold temperatures and volume of skin melt on the birefringence and interface distribution was analyzed for multilayered PS-PC-PS, PS-PMMA-PS, and PMMA-PC-PMMA molded disks obtained by CIM. The interface distribution and residual birefringence in the molded disks were measured. The measured interface distributions and the gapwise birefringence distributions in CIM disks were found to be in a fair agreement with the predicted interface distributions and the total residual birefringence obtained by the summation of the predicted frozen-in flow and thermal birefringence. POLYM. ENG. SCI., 55:88-106, 2015. (c) 2014 Society of Plastics Engineers