The present study attempted to numerically predict both the flow-induced and thermally-induced residual stresses and birefringence in injection or injection/compression molded center-gated disks. A numerical analysis system has been developed to simulate the entire process based on a physical modeling including a nonlinear viscoelastic fluid model, stress-optical law, a linear viscoelastic solid model, free volume theory for density relaxation phenomena and a photoviscoelasticity and so on. Part I presents physical modeling and typical numerical analysis results of residual stresses and birefringence in the injection molded center-gated disk. Typical distribution of thermal residual stresses indicates a tensile stress in the core and a compressive stress near the surface. However, depending on the processing condition and material properties, the residual stress sometimes becomes tensile on the surface, especially when fast cooling takes place near the mold surface, preventing the shrinkage from occurring. The birefringence distribution shows a double-hump profile across the thickness with nonzero value at the center: the nonzero birefringence is found to be thermally induced, the outer peak due to the shear flow and subsequent stress relaxation during the filling stage and the inner peak due to the additional shear flow and stress relaxation during the packing stage. The combination of the flow-induced and thermally-induced birefringence makes the shape of predicted birefringence distribution quite similar to the experimental one.