In the present study, a numerical calculation procedure based on the finite volume method was employed to simulate flow fields in double-branched tubes. The configuration was a tube with two vertical branches; the two branches were either on the same side or on the opposite side. The study focused on the baseline flow fields and the possible flow interaction between the two branches. The branching ratio and the branch /main tube diameter ratio were fixed in this study. The results showed that when the two branches were on the same side, the low/oscillating shear regions were found on the ventral walls of the branches and on the dorsal wall of the main tube distal to the branches. The flow field proximal to each branch was similar to that in a single-branched tube when the two branches were distant. When the branches were on the opposite side with the staggering distance S=0 (symmetric case), the low/oscillating shear regions were found on the lateral walls of the main tube. As S increased, the interaction between the two branches weakened, the low/oscillating shear regions were found on the lateral walls of the main tube to the side of the second branch. The flow field near the branch was significantly different from that of a single-branched tube. Care should be taken on localization of plaques in multi-branched vessels due to the flow pattern change. The numerical results were qualitatively consistent with what observed experimentally, by other investigators.