Numerical viscoelastic simulations were carried out using a K-BKZ type of separable integral constitutive equation. Both reversible and irreversible models were tried for several types of damping functions to calculate the annular extrudate behavior of high-density polyethylene (HDPE). There are two aims in this study; first, to clarify the properties of these dumping functions, and second, to investigate the influence of rheological characteristics on annular extrudate swell. In these numerical simulations, relaxation spectrum and shear viscosity were fixed, and the other characteristics were varied. The reversional response of the damping function mainly has an effect on the magnitude of the area swell even if the die is straight. The irreversible model expresses the experimental results of annular extrudate swell better than the reversible model. The accurate fitting of N1 by the damping model is important for predicting it. The magnitude of N1 predicted from the Wagner exponential model is lower than that of the PSM model, and the area swell shows the same tendency as N1. A modified PSM model that allows the N1 curve to shift can fit the magnitude of area swell. The relationship between the diameter and thickness of the extrudate depends on N2/N1, and it was estimated by simple linear elasticity of solids. The time-dependent viscosity varies with the type of damping function, and it influences the time-dependent swell.