When the change rate of the stress tensor is approximated by a time derivative, the constitutive equation is reduced to the linear function of a deformation rate tensor. The derived pseudo-viscosity decreases as relaxation time increases. To enhance the pseudo-viscosity matrix, a relaxation method is also presented. By the proposed formulation, the stress overshoot effect of a transient Couette flow is definitively simulated.