Abstract
A new singular finite element method for the viscoelastic flow problems was investigated. As for a small core region near singularity, we adopted a special γ^<-0.5> type element for elastic stress and pressure, and a special γ^<0.5> type element for velocity, where we intentionally removed a nodal point from the position of singularity, and then each element contained 12 nodes. It was found that the new method can predict successfully the die swell flow of viscoelastic fluids at a higher Weissenberg number over We=70 with less spurious oscillations near the die-lip.