A numerical modeling study is presented to investigate flow–focusing phenomena in a large-scale fracture network, constructed using field data collected from the unsaturated zone of Yucca Mountain, Nevada, the proposed repository site for high-level nuclear waste. The two-dimensional fracture network for an area of 100 by 150 m contains more than 20000 fractures. Steady state unsaturated flow in the fracture network is investigated for different boundary conditions and rock properties. Simulation results indicate that flow paths are generally vertical, and that horizontal fractures mainly provide pathways between neighboring vertical paths. In addition to fracture properties, flow–focusing phenomena are also affected by rock–matrix permeability, with lower matrix permeability leading to a high degree of flow focusing. The simulation results further indicate that the average spacing between flow paths in a layered system tends to increase and flow tends to become more focused, with depth.