A three-dimensional (3D) irregular and unstructured pore network was built using local topological and geometrical properties of an isometric bead pack imaged by means of a high-resolution X-ray computed micro-tomography technique. A pore network model was developed to analyze the 3D laminar/inertial (non-Darcy) flows at the mesoscopic (pore level) and macroscopic (after ensemble-averaging) levels. The non-linear laminar flow signatures were captured at the mesoscale on the basis of analogies with contraction and expansion friction losses. The model provided remarkably good predictions of macroscopic frictional loss gradient in Darcy and non-Darcy regimes with clear-cut demarcation using channel-based Reynolds number statistics. It was also able to differentiate contributions due to pore and channel linear losses, and contraction/expansion quadratic losses. Macroscopic mechanical dispersion was analyzed in terms of retroflow channels, and transverse and longitudinal Peclet numbers. The model qualitatively retrieved the Peclet-Reynolds scaling law expected for heterogeneous networks with predominance of mechanical dispersion. Advocated in watermark is the potential of pore network modeling to build a posteriori constitutive relations for the closures of the more conventional macroscopic Euler approaches to capture more realistically single-phase flow phenomena in fixed-bed reactor applications in chemical engineering. (C) 2013 Elsevier B.V. All rights reserved.