A two-dimensional general rate model of liquid chromatography incorporating slow rates of adsorption-desorption kinetics, axial and radial dispersions, and core-shell particles is formulated. Radial concentration gradients are generated inside the column by considering different regions of injection at the inlet. Analytical solutions are obtained for a single-component linear model by simultaneously utilizing the Laplace and Hankel transformations for the considered two sets of boundary conditions. These linear solutions are useful for simulating liquid-chromatographic columns with diluted or small-volume samples and those in which radial concentration gradients are significant. To gain further insight into the process, analytical moments are also deduced from the Laplace-Hankel-domain solutions. For situations of concentrated and large-volume samples, which are not solvable analytically, formulation of nonlinear models is necessary. In this study, a semidiscrete, high-resolution, finite-volume scheme is extended to approximate the resulting nonlinear-model equations for multicomponent mixtures. The performance of the column is analyzed by implementing a specified criterion of performance. A few numerical case studies are conducted to inspect the effects of the model parameters on the elution profiles.