A nodal high-order discontinuous Galerkin finite element (DG-FE) method is presented to solve the equilibrium-dispersive model of chromatography with arbitrary high-order accuracy in space. The method can be considered a high-order extension to the total variation diminishing (TVD) framework used by Javeed et al. (2011a,b, 2013) with an efficient quadrature-free implementation. The framework is used to simulate linear and non-linear multicomponent chromatographic systems. The results confirm arbitrary high-order accuracy and demonstrate the potential for accuracy and speed-up gains obtainable by switching from low-order methods to high-order methods. The results reproduce an analytical solution and are in excellent agreement with numerical reference solutions already published in the literature.