Bitumen residues in the oil sand tailings can be a threat to the environment that separating them from tailings before disposal is crucial. This study establishes an Eulerian-Eulerian computational fluid dynamics (CFD) model for an industrial-scale oil sand tailings pipeline. A comprehensive sensitivity analysis was conducted on the selection of carrier-solid and solid-bitumen drag models. The combination of small and large particle sizes (i.e., 75 & 700 µm) and bitumen droplet size (i.e., 400 µm) provided good agreement with field data in velocity profiles and pressure drop. The validated model was subsequently extended to investigate the influence of the secondary phase (i.e., bitumen droplets and bubbles) on flow characteristics in a tailing pipeline. The investigation covered a range of bitumen droplet size (100-400 µm), bitumen fraction (0.0025-0.1), bubble size (5-1000 µm), and bubble fraction (0.0025-0.3) and their influences on the velocity, solids, and bitumen distribution are revealed. For an optimum bubble size of 500 µm, a maximum recovery of 59% from the top 50 % and 83 % from the top 75 % of the pipe cross-section was obtained. The present study demonstrates the preferential distribution of bitumen and provides valuable insight on bitumen recovery from an industrial-scale tailings pipeline.