  Using the HYDJET++ model, we measure the higher-order flow harmonics
$v_{n}$(n=2,3,4) of (multi-) strange hadrons in Au+Au collisions at
$\sqrt{s_{NN}}$ = 200 GeV and Pb+Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.
We have compared our model results and experimental data at RHIC and LHC
energies. The model reproduces the higher flow harmonics of (multi-) strange
hadrons as a function of $p_{T}$, centrality, and quark content. We have
studied and discussed mass ordering of flows $v_{n}$(n=2,3,4) among
$\pi^{+}+\pi^{-}$, $K^{-}+K^{+}$, $K_{s}^{0}$, $p+\overline{p}$,
$\Lambda+\overline{\Lambda}$, $\Xi^{-}+\overline{\Xi}^{+}$,
$\Omega^{-}+\overline{\Omega}^{+}$ at low $p_{T}$ and the baryon-meson grouping
at intermediate $p_{T}$. NCQ scaling using the HYDJET++ model and its
comparison with experimental results at RHIC and LHC energies have also been
discussed. Study of flow harmonics behaviour in RHIC and LHC energy regimes
will provide an additional insight into the dynamics of anisotropic flow and
the effect of radial flow expansion in the system.
