Discrete fracture model has become one of the alternatives for the classical continuum model to simulate the irregular aspects of the fluid flow and the solute transport in fractured rocks. It is based on the assumptions that the discharge in a single fracture is proportional to the cube of the aperture and the fractured rock can be represented by the statistical assemblage of such single fractures. This study is intended to evaluate the effect of the fracture junction on the cubic law. Numerical solution of flow in junction system was obtained by using the Boundary-Fitted Coordinate System (BFCS) method. Results with different intersection angles in crossing fractures show that the geometry of the junction affects the discharge pattern under the same simulation conditions. Therefore, strict numerical and experimental examinations on this subject are required.