  The study of anisotropic harmonic flow coefficients $ v_{n}$(n=2,3,4) is
performed in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV under Monte Carlo
HYDJET++ model (HYDrodynamics plus JETs) framework. Anisotropic flow of
identified particles and correlation between the azimuthal harmonic flow
amplitudes is presented. Here, we have considered body-body and tip-tip type of
geometrical configurations for Xe-Xe collision systems. The kinematic ranges
$|\eta|<0.8$, $0<p_{T}<5.0$ GeV/c, and $|\delta / \eta|> 2$ are considered. The
results have been shown for seven classes of centrality and compared with the
ALICE experimental data. The anisotropic flow of identified charged particles
show a strong centrality dependence. Mass ordering is observed for
$v_{2},v_{3}$ and $v_{4}$. Mass ordering is different for different ranges of
transverse momentum $p_{T}$. Strong correlation is observed between
$v_{3}-v_{2}$, $v_{4}-v_{2}$, and $v_{4}-v_{3}$. Such correlation is centrality
dependent and is different in different centrality windows. The anisotropic
flow coefficients show a clear dependence on the total charged particle
multiplicity. HYDJET++ model justifies experimental data well enough.
