Chromatographic separation processes need efficient simulation methods, especially for nonlinear adsorption isotherms such as the Langmuir isotherms which imply the formation of concentration shocks. The focus of this paper is on the space−time conservation element/solution element (CE/SE) method. This is an explicit method for the solution of systems of partial differential equations. Numerical stability of this method is guaranteed when the Courant−Friedrichs−Lewy condition is satisfied. To investigate the accuracy and efficiency of this method, it is compared with the classical cell model, which corresponds to a first-order finite volume discretization using a method of lines approach (MOL). The evaluation is done for different models, including the ideal equilibrium model and a mass transfer model for different adsorption isotherms—including linear and nonlinear Langmuir isotherms—and for different chromatographic processes from single-column operation to more sophisticated simulated moving bed (SMB) processes for the separation of binary and ternary mixtures. The results clearly show that CE/SE outperforms MOL in terms of computational times for all considered cases, ranging from 11-fold for the case with linear isotherm to 350-fold for the most complicated case with ternary center-cut eight-zone SMB with Langmuir isotherms, and it could be successfully applied for the optimization and control studies of such processes.