We conducted direct numerical simulations (DNSs) of turbulent flow over three-dimensional sinusoidal roughness in a channel. A passive scalar is present in the flow with Prandtl number P r = 0.7, to study heat transfer by forced convection over this rough surface. The minimal-span channel is used to circumvent the high cost of simulating high Reynolds number flows, which enables a range of rough surfaces to be efficiently simulated. The near-wall temperature profile in the minimal-span channel agrees well with that of the conventional full-span channel, indicating it can be readily used for heattransfer studies at a much reduced cost compared to conventional DNS. As the roughness Reynolds number, k + , is increased, the Hama roughness function, ∆U + , increases in the transitionally rough regime before tending towards the fully rough asymptote of κ −1 m log(k + ) + C, where C is a constant that depends on the particular roughness geometry and κ m ≈ 0.4 is the von Kármán constant. In this fully rough regime, the skin-friction coefficient is constant with bulk Reynolds number, Re b . Meanwhile, the temperature difference between smooth-and rough-wall flows, ∆Θ + , appears to tend towards a constant value, ∆Θ + F R . This corresponds to the Stanton number (the temperature analogue of the skin-friction coefficient) monotonically decreasing with Re b in the fully rough regime. Using shifted logarithmic velocity and temperature profiles, the heat transfer law as described by the Stanton number in the fully rough regime can be derived once both the equivalent sand-grain roughness k s /k and the temperature difference ∆Θ + F R are known. In meteorology, this corresponds to the ratio of momentum and heat transfer roughness lengths, z 0m /z 0h , being linearly proportional to the inner-normalised momentum roughness length, z + 0m , where the constant of proportionality is related to ∆Θ + F R . While Reynolds analogy, or similarity between momentum and heat transfer, breaks down for the bulk skin-friction and heat-transfer coefficients, similar distribution patterns between the heat flux and viscous component of the wall shear stress are observed. Instantaneous visualisations of the temperature field show a thin thermal diffusive sublayer following the roughness geometry in the fully rough regime, resembling the viscous sublayer of a contorted smooth wall.