We tried to predict the multilayer parison shape at pinch-off stage in extrusion blow molding by nonisothermal and purely viscous non-Newtonian flow simulation using the finite element method (FEM). We assumed the parison deformation as a flow problem. The Carreau model was used as the constitutive equation and FEM was used for calculation method. Multi layer parison used in this simulation was composed of high-density polyethylene (HDPE) as inner and outer layers and low-density polyethylene (LOPE) of which viscosity is five times lower than HDPE as a middle layer. We discussed multilayer parison shape in pinch-off region. The results obtained are as follows; the parison shape of each layer was clearly visible in the pinch-off during the mold closing. In addition, the distribution of parison thickness ratios for each layer was located for a large deformation near the pinch-off region. The melt viscosity for each layer has an influence on the melt flow in the pinch-off region. In a comparison with an experimental data of parison thickness ratios, the simulation results are larger than the experimental data. These simulation results obtained are in good agreement with the experimental data in consideration of the standard deviations. POLYM. ENG. SCI., 50:1476-1484, 2010. (C) 2010 Society of Plastics Engineers