Event-by-event viscous hydrodynamics is combined with heavy quark energy loss models to compute heavy flavor flow cumulants $v_2{2}$, $v_3{2}$, and $v_2{4}$ as well as the nuclear modification factors of $D^0$ and $B^0$ mesons in PbPb collisions at 2.76 TeV. Our results indicate that bottom quarks can flow as much as charm quarks in the $p_T$ range 8--30 GeV.