Numerical analysis of seepage law for radial fluid flow in a single fracture: Comparison between smooth and rough fractures

This paper investigated the differences for the hydraulic characteristics in a single fracture between using the Navier–Stokes (N–S) equation and Darcy's law, which would be benefit to understand the seepage mechanism in the fracture. A numerical model of the radial flow was established considering the aperture size and water injection flow rate. Some conclusion could be given. First, the Darcy's law only described the seepage characteristics when the flow rate was small when the flow rate and pressure response have a linear relationship. While the N–S equation could describe the linear and nonlinear seepage characteristics, resulting in a better model of the actual fracture seepage flow. Second, the aperture size had a limited influence on the water pressure and seepage velocity inside the fracture when the flow rate was small. It began to have a significant impact influence on the seepage characteristics inside the fracture with the aperture increased. Third, the flow–pressure response conformed to the Forchheimer equation in the fracture. The critical Reynold number would decrease from 1.2 to 0.0116 when fracture aperture decreased from 3 to 0.5 mm using the N–S equation. The degree of nonlinearity of the fluid flow increases with fracture roughness increasing. This work gave a guidance to the difference in the two seepage theories and correction for the result by Darcy law, which was widely used in the engineering calculation.