The method of volume averaging is applied to estimate the Taylor-Aris dispersion tensor of solute advected in columns consisting of ordered pillar arrays with wall retention of the type used in chromatographic separation. The appropriate closure equations are derived and solved in a unit cell with periodic boundary conditions to obtain the dispersion tensor (or the reduced plate height) as a function of the Peclet number (reduced velocity); pillar pattern, shape and size; partition coefficient; and resistance to mass transfer. The contributions of the velocity profile, the wall adsorption, and the mass transfer resistance to the dispersion tensor are identified and delineated. The model is verified by comparing its predictions and obtaining favorable agreement with results of direct numerical simulations and with experimental data for columns containing ordered pillars. The model is then used to study the effect of pillars' shape and pattern on the longitudinal dispersion coefficient (plate height).