Three-dimensional Reynolds-averaged Navier–Stokes (RANS) solutions from CFX were utilized to investigate the leakage flow characteristics in the labyrinth honeycomb seal of steam turbines. At first, the accuracy and reliability of the utilized RANS approach was demonstrated using the published experimental data of the honeycomb seal. It showed that the utilized numerical method has sufficient precision to predict the leakage performance in seals. Then a range of sealing clearances, cell diameters, cell depths, rotation speeds, and pressure ratios were investigated to determine how these factors affect the leakage flow rate of the labyrinth honeycomb seal. The computed leakage flow rate increased with increasing sealing clearance and pressure ratios. Furthermore, the results show that the studied labyrinth honeycomb seal has the optimum sealing performance in the case of honeycomb cell diameter equals labyrinth step width, and the ratio of the honeycomb cell depth to honeycomb cell diameter is 0.93 under the designed condition. The flow pattern of each case is also illustrated to describe the leakage flow characteristics in labyrinth honeycomb seals.