Subsurface flows are influenced by the presence of faults and large fractures which actas preferential paths or barriers for the flow. In literature models were proposed tohandle fractures in a porous medium as objects of codimension 1. In this work we considerthe case of a network of intersecting fractures, with the aim of deriving physicallyconsistent and effective interface conditions to impose at the intersection betweenfractures. This new model accounts for the angle between fractures at the intersectionsand allows for jumps of pressure across intersections. This fact permits to describe theflow when fractures are characterized by different properties more accurately with respectto other models that impose pressure continuity. The main mathematical properties of themodel, derived in the two-dimensional setting, are analyzed. As concerns the numericaldiscretization we allow the grids of the fractures to be independent, thus in generalnon-matching at the intersection, by means of the extended finite element method(XFEM). This increases the flexibility of the method in the case of complexgeometries characterized by a high number of fractures.