Direct numerical simulation (DNSs) are used to systematically investigate applicability of minimal channel approach for characterization of roughness-induced drag in irregular rough surfaces. Roughness is generated mathematically using a random algorithm, in which the power spectrum (PS) and probability density function (PDF) of surface height function can be prescribed. 12 different combinations of PS and PDF are examined and both transitionally and fully rough regimes are investigated (roughness heights varies in the range $k^+$ = 25 -- 100).
It is demonstrated that both roughness function ($\Delta U^+$) and zero-plane displacement can be predicted within $\pm5%$ accuracy using DNS in properly sized minimal channels. Notably, the predictions do not deteriorate when a limited range of large horizontal roughness scales are filtered out due to the small channel size (here up to 10% of original roughness height spectral energy based on 2D PS). Additionally, examining the results obtained from different random realizations of roughness shows that a certain combination of PDF and PS leads to a nearly unique $\Delta U^+$ for deterministically different surface topographies.
In addition to the global flow properties, the distribution of time-averaged surface force exerted by the roughness onto the fluid is calculated and compared for different cases.
It is shown that patterns of surface force distribution over irregular rough surfaces can be well captured when the sheltering effect is taken into account. This is made possible applying the sheltering model proposed by Yang et al. to each specific roughness topography.
Furthermore, an analysis of the coherence function between roughness height and surface force distributions reveals that the coherence drops at larger streamwise wavelengths, which can be an indication that very large horizontal scales are less dominant in contributing to the skin friction drag.
Finally, some existing roughness correlations are assessed using the present roughness dataset, and it is shown that the correlation predictions for the values of equivalent sand-grain roughness mainly lie within $\pm30%$ error in comparison to the DNS results.