Heterogeneous aperture distributions are an intrinsic characteristic of natural fractures. The presence of highly heterogeneous aperture distributions can lead to flow channeling, thus influencing the macroscopic behavior of the fluid flow. High-fidelity numerical simulation tools are needed for realistic simulation of fracture flow when such features are present. Here, focus is set on the role of mechanical fracture closure for fluid flow and appropriate simulation by a fully hydro-mechanically (HM) coupled numerical model. In a laboratory experiment, an artificial fracture in a granodiorite sample is created. During different sequential loading cycles, the development of fracture closure, contact area and contact stress are examined. Constant fluid flow rate injection into the center of the rough fracture is modelled to investigate the impact of fracture closure on the flow field and injection pressure. Results show that the numerical framework for heterogeneous fracture surfaces allows for reproducing experimental data of dry, mechanical tests at the laboratory scale, and it may offer advanced understanding and prediction of the behavior of reservoirs that are subject to high-pressure fluid injections.