Rothe method and numerical analysis for a new class of fractional differential hemivariational inequality with an application