An extended method of semidiscrete high-resolution finite volume is used in this paper to obtain numerical solutions for a formulated nonlinear lumped kinetic model of liquid chromatographic process to examine the effect of chromatographic column overloading gradient elution considering core-shell particles. The model constitutes linear solvent strength (LSS), Henry’s constant, coefficient of nonlinearity, and coefficient of axial dispersion. The effects of modulator concentration changes for the elution of single and two components are analyzed. The advantages of introducing gradient elution against isocratic elution in terms of core radius fraction are investigated intensively. Numerical temporal moments are generated from the solutions obtained for a more in-depth examination of the considered model. Moreover, multiple forms of a single- and two-component mixture are generated to analyze the influences of core radius fractions on gradient elution. For example, the obtained results are utilized to investigate the effects of the slope of gradient, concentration of modulator, solvent strength parameter, coefficient of nonlinearity, coefficient of mass transfer, and coefficient of axial dispersion on the profiles of concentration in order to improve the process performance using optimal core radius fraction parameter values.