We present calculations of charm and beauty production at Tevatron within the framework of k _T -factorization, using the unintegrated gluon distributions as obtained from the Linked Dipole Chain model. The analysis covers transverse momentum and rapidity distributions and the azimuthal correlations between b and b̅ quarks (or rather muons from their decay) which are powerful tests for the different unintegrated gluon distributions. We compare the theoretical results with recent experimental data taken by DØ and CDF collaborations at the Tevatron Run I and II.