A numerical method has been proposed to predict 3-dimensional flow in a duct system with multiple outlets. For the duct system, it is supposed that the pressure values are given at multiple outlets while the velocity profile is given at a inlet. To maintain the continuity of pressure distribution between main and branch duct, present method allows that the pressure value taken from analysis of branch duct can be converted to the main duct analysis. The result from present method which can handle the pressure boundary condition closely coincided with that from regular method which can handle the velocity boundary condition only. Furthermore the flow distribution from present method showed good agreement with that from the single block method. From the comparison of the present method with the total pressure method used for engineering duct design, 13% of discrepancy in pressure loss was shown between the main duct inlet and the branch duct outlet.