Presence of fracture roughness and occurrence of nonlinear flow complicate fluid flow through rock fractures. This paper presents a qualitative and quantitative study on the effects of fracture wall surface roughness on flow behavior using direct flow simulation on artificial fractures. Previous studies have highlighted the importance of roughness on linear and nonlinear flow through rock fractures. Therefore, considering fracture roughness to propose models for the linear and nonlinear flow parameters seems to be necessary. In the current report, lattice Boltzmann method is used to numerically simulate fluid flow through different fracture realizations. Flow simulations are conducted over a wide range of pressure gradients through each fracture. It is observed that creeping flow at lower pressure gradients can be described using Darcy's law, while transition to inertial flow occurs at higher pressure gradients. By detecting the onset of inertial flow and regression analysis on the simulation results with Forchheimer equation, inertial resistance coefficients are determined for each fracture. Fracture permeability values are also determined from Darcy flow as well. According to simulation results through different fractures, two parametric expressions are proposed for permeability and inertial resistance coefficient. The proposed models are validated using 3D numerical simulations and experimental results. The results obtained from these two proposed models are further compared with those obtained from the conventional models. The calculated average absolute relative errors and correlation coefficients indicate that the proposed models, despite their simplicity, present acceptable outcomes; the models are also more accurate compared to the available methods in the literature.