We study helicity polarization through the $(3+1)$-dimensional relativistic viscous hydrodynamic models at $\sqrt{{s}_{NN}}=200$ GeV $\text{Au}+\text{Au}$ collisions. Similar to the local spin polarization, we consider the helicity polarization beyond global equilibrium and investigate the contributions induced by thermal vorticity, shear viscous tensor, and the fluid acceleration. We find that the local helicity polarization induced by thermal vorticity dominates over other contributions. It also implies that, in low-energy collisions, the fluid vorticity as part of thermal vorticity may play the crucial role in the total helicity polarization. Such a finding could be useful for probing the local strength of vorticity in rotational quark gluon plasmas by measuring helicity polarization. Our simulation confirms the strict space reversal symmetry, whereas we also compare our numerical results with approximated relations derived from ideal Bjorken flow. Our studies also provide a baseline for the future investigation on local parity violation through the correlations of helicity polarization.