A systematic investigation of algorithm impact in preparative chromatography with experimental verifications.

Computer-assisted optimization of chromatographic separations requires finding the numerical solution of the Equilibrium-Dispersive (ED) mass balance equation. Furthermore, the competitive adsorption isotherms needed for optimization are often estimated numerically using the inverse method that also solves the ED equations. This means that the accuracy of the estimated adsorption isotherm parameters explicitly depends on the numerical accuracy of the algorithm that is used to solve the ED equations. The fast and commonly used algorithm for this purpose, the Rouchon Finite Difference (RFD) algorithm, has often been reported not to be able to accurately solve the ED equations for all practical preparative experimental conditions, but its limitations has never been completely and systematically investigated. In this study, we thoroughly investigate three different algorithms used to solve the ED equations: the RFD algorithm, the Orthogonal Collocation on Finite Elements (OCFE) method and a Central Difference Method (CDM) algorithm, both for increased theoretical understanding and for real cases of industrial interest. We identified discrepancies between the conventional RFD algorithm and the more accurate OCFE and CDM algorithms for several conditions, such as low efficiency, increasing number of simulated components and components present at different concentrations. Given high enough efficiency, we experimentally demonstrate good prediction of experimental data of a quaternary separation problem using either algorithm, but better prediction using OCFE/CDM for a binary low efficiency separation problem or separations when the compounds have different efficiency. Our conclusion is to use the RFD algorithm with caution when such conditions are present and that the rule of thumb that the number of theoretical plates should be greater than 1000 for application of the RFD algorithm is underestimated in many cases.

[1]  Georges Guiochon,et al.  Application of orthogonal collocation on finite elements in the simulation of non-linear chromatography , 1991 .

[2]  Torgny Fornstedt,et al.  An improved algorithm for solving inverse problems in liquid chromatography , 2006, Comput. Chem. Eng..

[3]  Concentration dependence of lumped mass transfer coefficients linear versus non-linear chromatography and isocratic versus gradient operation. , 2003, Journal of chromatography. A.

[4]  Andreas Seidel-Morgenstern,et al.  Experimental determination of single solute and competitive adsorption isotherms. , 2004, Journal of chromatography. A.

[5]  Anita M. Katti,et al.  Fundamentals of Preparative and Nonlinear Chromatography , 1994 .

[6]  Georges Guiochon,et al.  Numerical determination of the competitive isotherm of enantiomers. , 2003, Journal of chromatography. A.

[7]  M. Morbidelli,et al.  Modeling fixed-bed adsorption columns through orthogonal collocations on moving finite elements , 1997 .

[8]  G. Guiochon,et al.  Comparative modeling of breakthrough curves of bovine serum albumin in anion-exchange chromatography. , 2001, Journal of chromatography. A.

[9]  Krzysztof Kaczmarski,et al.  Estimation of adsorption isotherm parameters with inverse method--possible problems. , 2007, Journal of chromatography. A.

[10]  S. Rohani,et al.  Numerical determination of competitive adsorption isotherm of mandelic acid enantiomers on cellulose-based chiral stationary phase. , 2008, Journal of chromatography. A.

[11]  A. K. Ray,et al.  Determination of competitive adsorption isotherm parameters of pindolol enantiomers on alpha1-acid glycoprotein chiral stationary phase. , 2006, Journal of chromatography. A.

[12]  Pierre Rouchon,et al.  Numerical Simulation of Band Propagation in Nonlinear Chromatography , 1987, Preparative-Scale Chromatography.

[13]  G. D. Byrne,et al.  VODE: a variable-coefficient ODE solver , 1989 .

[14]  T. Fornstedt,et al.  Validation of the accuracy of the perturbation peak method for determination of multicomponent adsorption isotherm parameters in LC. , 2004, Analytical chemistry.

[15]  J. Villadsen,et al.  Solution of differential equation models by polynomial approximation , 1978 .

[16]  T. Fornstedt,et al.  A quest for the optimal additive in chiral preparative chromatography. , 2009, Journal of chromatography. A.