Abstract A Taylor–Galerkin finite element time marching scheme is derived to numerically simulate the flow of a compressible and nonisothermal viscoelastic liquid between eccentrically rotating cylinders. Numerical approximations governing constitutive equations are computed over custom refined unstructured grid piecewise linear Galerkin elements. An original extension DEVSS formulation for flu...