We report our recent theoretical calculations for $b\bar{b}$ dijet production in high-energy nuclear collisions. The NLO+parton shower (PS) event generator SHERPA has been employed to provide the pp baseline of $b\bar{b}$ dijet production. A framework which combines the Langevin transport equation to describe the evolution of heavy quarks and the higher-twist scheme to consider the inelastic energy loss of both light and heavy partons has been implemented. Within this framework, we present the theoretical results for the transverse momentum imbalance $x_J=p_{T2}/p_{T1}$ both for inclusive dijets and $b\bar{b}$ dijets in Pb+Pb collisions at $\sqrt{s_{NN}}=5.02$~TeV. The energy loss of b-jets is expected to shift $x_J$ to smaller values relative to the pp reference which is consistent with the CMS data. In addition, we show the medium modification for angular correlation of $b\bar{b}$ dijets in A+A collisions at $\sqrt{s_{NN}}=5.02$~TeV. We observe a stronger suppression in the small $\Delta \phi=|\phi_{b1}-\phi_{b2}|$ region where the gluon splitting processes dominate relative to the large $\Delta \phi$ region. The difference leads to a modest suppression on the near-side ($\Delta \phi\sim 0$) and enhancement on the away-side ($\Delta\phi\sim\pi$).