The Boundary element approximation for the parabolic variational inequalities of the second kind is discussed.

The quality of approximation is as good as polygonal boundary reduction, and it can process the data online in a constant buffer.

The approximation order of the interpolation with the appropriate boundary conditions and constraints is analyzed.