Progress in roughness research, mapping any given roughness geometry to its fluid dynamic behaviour, has been hampered by the lack of accurate and direct measurements of skin-friction drag, especially in open systems. The Taylor-Couette (TC) system has the benefit of being a closed system, but its potential for characterizing irregular, realistic, 3-D roughness has not been previously considered in depth.Here, we present direct numerical simulations (DNSs) of TC turbulence with sand grain roughness mounted on the inner cylinder. The model proposed by Scotti (Phys. Fluids, vol. 18, 031701, 2006) has been improved to simulate a random rough surface of monodisperse sand grains, which is characterized by the equivalent sand grain height k s . Taylor numbers range from T a = 1.0 × 10 7 (corresponding to Re τ = 82) to T a = 1.0 × 10 9 (Re τ = 635). We focus on the influence of the roughness height k + s in the transitionally rough regime, through simulations of TC with rough surfaces, ranging from k + s = 5 up to k + s = 92, where the superscript '+' indicates non-dimensionalization in viscous units.We analyse the global response of the system, expressed both by the dimensionless angular velocity transport N u ω and by the friction factor C f . An increase in friction with increasing roughness height is accompanied with enhanced plume ejection from the inner cylinder. Subsequently, we investigate the local response of the fluid flow over the rough surface. The equivalent sand grain roughness k + s is calculated to be 1.33k, where k is the size of the sand grains. We find that the downwards shift of the logarithmic layer, due to transitionally rough sand grains exhibits remarkably similar behavior to that of the Nikuradse (VDI-Forschungsheft 361, 1933) data of sand grain roughness in pipe flow, regardless of the Taylor number dependent constants of the logarithmic layer.