Comprehension of Effects of Intraparticle Effective Diffusivity on Simulated Moving Bed

The influence of intraparticle effective diffusivity on simulated moving bed (SMB) is investigated in this work. The intraparticle effective diffusivity has significant influence on the separation process of SMB and it affects the upper and left sides of separation region of SMB in different styles. The result is explained according to the separation of binary mixture in single column in C1-C2 space. This indicates the necessity of considering the effects of intraparticle mass transfer processes for chemical engineers during their designing the operating conditions of SMB. *Corresponding author: Hong Gao, Department of Physics, Anshan Normal University, Anshan, Liaoning, 114005, PR China, Tel: +8615641280037; E-mail: gh71@163.com Received June 23, 2017; Accepted June 27, 2017; Published June 30, 2017 Citation: Gao H (2017) Comprehension of Effects of Intraparticle Effective Diffusivity on Simulated Moving Bed. J Chromatogr Sep Tech 8: 371. doi: 10.4172/21577064.1000371 Copyright: © 2017 Gao H. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.


Introduction
Significant effects of intraparticle mass transfer on separation result in single chromatographic column have been realized by many researchers [1][2][3][4][5][6][7][8][9][10]. But the effects of intraparticle mass transfer on simulated moving bed (SMB) process have seldom been reported. In practice, the flow rate varies in different zones of SMB and sometimes is high when the SMB is run, and then the influence of intraparticle mass transfer is important and becomes significant.
It is the goal pursued by every user of SMB that choose the appropriate operating conditions for SMB and then obtain the optimal separation result. And practically, the complete separation region in m 2 -m 3 plane is generally the key for the determination of the operating conditions of SMB, where m 2 and m 3 are the ratios of net fluid flow rate to net solid flow rate in second and third zone of SMB, respectively. Under ideal conditions, the complete separation region of SMB can be easily obtained according to the results in literatures [11][12][13]. While for the situations that the mass transfer resistances exist, such separation region has to be calculated based on the numerical solution of the model of SMB.
The General Rate (GR) model which involves the effect of intraparticle mass transfer processes is used as the model for single column of SMB in this work. And a modified equilibrium dispersive (MED) model of which the numerical solution is equivalent to that of GR model is used in this work to simplify the calculations. The influence of different levels of intraparticle mass transfer parameters, i.e., effective diffusivity and external mass transfer coefficient on separation region (purities above 99%) for binary mixture in four-zone SMB is investigated and then explained by the results in the C 1 -C 2 space of binary mixture separation under same conditions in single column.
For single column of SMB units, the General Rate (GR) model [1], which involves the effects of intraparticle mass transfer processes but neglects the process of adsorption/desorption, is adopted in this work, because in most cases the kinetic process of adsorption/desorption is fast enough and its contribution to the overall mass transfer processes can be neglected [1]. The model parameters are considered to be same for two components for the sake of simplification. And then the model equations can be written as follows, And the Langmuir isotherm is considered: The initial conditions are: The boundary conditions are The calculation of separation region of SMB using GR model directly as the single column model for SMB units is not easy. In fact, it is computation intensive and extremely time-consuming. And then a simplified but equivalent method for solving the problem is expected.
A modified equilibrium dispersive (MED) model has been developed [14] and can replace the GR model to simplify the numerical calculations. Therefore, the total band broadening originated from axial dispersion and intraparticle mass transfer processes can be considered equal to the band broadening resulted from an equivalent axial dispersion, which represented by axial dispersion coefficient i eq D , in the modified equilibrium dispersive (MED) model. The modified equilibrium dispersive (MED) model is as follows, ). And the initial and boundary conditions for C i and q i are same with that of GR model above.
The HETP equation for the GR model is deduced from the moment equations of the model [15], and it is written as below, A four-zone SMB which involves two columns per zone is considered in this work. And the values of some parameters needed for calculations are listed as follows. Column length L=12 cm, column diameter R=1 cm, external porosity .05 mg/mL, Langmuir isotherm parameters G 1 =2, G 2 =3, b 1 =0.5 mL/mg, b 2 =0.5 mL/mg, axial dispersion coefficient D L =0.01 cm 2 /min. The operating parameters for the first and fourth zone of SMB are considered as far away from their critical values, and the switch time of SMB is 2 min, the flow rates in the first and fourth zones are 33.9 mL/min and 18.8 mL/min, respectively.
To obtain the edge of separation region of SMB, the lattice with The modified equilibrium dispersive (MED) model above is employed to obtain the numerical solutions equivalent to that of the GR model. Firstly, Equation (3) is discretized with the method of orthogonal collocation on finite elements. And then, the discretized ordinary differential equations were solved by Runge-Kutta method of order four.
The separation region of SMB is calculated at different levels of intraparticle mass transfer parameters, i.e., effective diffusivity D eff and external mass transfer coefficient k ext . The ranges of values of these parameters were determined according to the literatures [2,[4][5][6].
The separation regions of SMB calculated at four levels of D eff are shown in Figure 1. The triangles in the figure are corresponding to the separation regions under different conditions, and are named as Triangle 1, Triangle 2, Triangle 3 and Triangle 4, respectively and in the order of decrease of D eff . It is clear that the triangle shrinks remarkably when D eff decreases from 1 × 10 -4 cm 2 /min to 5 × 10 -6 cm 2 /min. And more specifically, with decrease of D eff , the upper side of the triangle moves downwards remarkably and almost vertically. The upper sides of the triangles obtained under higher D eff turn gradually upward slightly with increase of m 3 , but for the triangles obtained under lower D eff , this gradual raise can not be observed. While for the left side of triangle, it moves rightward when D eff decreases. But the points in left side of triangle move rightward in different distances. The points near the diagonal move farther than those far from it. In summary, the decrease of D eff affects the two sides of the separation region in different ways; the upper side moves downwards farther with its slope decreases slightly; while for the left side, its slope decreases remarkably and at the same time it moves rightward but not so far as the upper side does.
The external mass transfer coefficient has similar effect on the separation region of SMB. The result is shown as Figure 2. With the value of k ext decreases from 10 cm/min to 0.1 cm/min, the corresponding calculated separation region of SMB shrinks slightly. The upper side of the triangle moves downwards and the left side of the triangle moves rightward. And it seems that the reduction of the area of the separation region is mainly due to downward movement of the upper side of the separation region rather than the rightward movement of the left one. And compared with the triangles in Figure 1, the difference between the triangles in Figure 2 is much less.
The results shown in Figures 1 and 2 indicate that slowing the intraparticle mass transfer processes would reduce the separation region of SMB. And the change of the separation region is more sensitive to the decrease of effective diffusivity than that of external mass transfer coefficient in their ranges which generally reported in literatures.
In fact, the results of the effect of intraparticle mass transfer parameters on the separation region of SMB mentioned above is easy to understand when the result of the binary mixture separation in single column in C 1 -C 2 space is considered. The sample feeding style is a wide rectangle injection of binary mixture solution which is similar with that of SMB. The effects of effective diffusivity on the separation results in C 1 -C 2 space of a wide rectangular pulse of binary mixture solution eluted in single column are presented here for the explanation of the point.  Figure 3 shows the results of different levels of intraparticle effective diffusivity when the binary mixture solutions of lower and higher concentrations are fed, respectively. The following points can be concluded according to Figure 3. At first, different from two straight lines in C 1 -C 2 space under ideal condition [16], the lines in C 1 -C 2 space curve and are tangent with the coordinate axes when they approach the axes. And with the decrease of effective diffusivity, the tangent points get closer to the coordinate origin. The shorter distances between the tangent points and the coordinate origin mean the lower concentrations of two pure components taken from the binary mixture. That is to say, with the decrease of effective diffusivity, the complete separation of binary mixture in single column gets more difficult. Secondly, the influences of the nonlinearity which caused by the increase of feed concentration on the fronts and the rears of eluted band profiles are different. When lower concentration of binary mixture being fed, with the decrease of effective diffusivity, the straight part of the line in C 1 -C 2 space corresponding to the fronts of eluted band profile shortens and disappears eventually. While with the higher sample concentration being fed, all these lines are curved. The difference between the lines under different levels of D eff is always obvious whether for the situation of lower or higher concentration of binary mixture being fed. In other words, the increase of the feed concentration can only slightly reduce the difference between the fronts of the eluted band profiles calculated under different levels of effective diffusivity. That is to say, the influence of effective diffusivity is stronger than that of concentration nonlinearity on the lines in C 1 -C 2 space corresponding to the fronts of eluted band profiles. But the situation for the lines in C 1 -C 2 space which corresponding to the rears of eluted band profiles is much different. When the feed concentration increases, the difference between these lines almost disappears completely. And this indicates that the increased concentration nonlinearity strongly diminishes the difference between the rears of eluted band profiles calculated at different levels of effective diffusivity. Now, the effects of effective diffusivity on the separation region of SMB can be interpreted based on the understanding of the results in C 1 -C 2 space for binary mixture separation in single column.
Firstly, Figure 3a and 3b shows that with decrease of effective diffusivity, the complete separation of binary mixture in single column gets more difficult. And this can be the explanation for the result shown in Figure 1 that the separation region reduces when effective diffusivity decrease.
Secondly, as we know, the left side of the triangle (separation region) for SMB indicates the left border of the region in m 2 -m 3 plane beyond which the pure more retained component (C 2 ) can not be obtained from extract port of SMB. So, the behavior of the lines in C 1 -C 2 space corresponding to the rears of eluted band profiles in single column under different levels of effective diffusivity can be employed here to understand the effect of effective diffusivity on the left boundary of separation region of SMB. On the other hand, the upper side of the triangle in m 2 -m 3 plane is the border of the separation region for SMB beyond which the pure less retained component (C 1 ) can not be collected from raffinate port of SMB. And the results discussed above about the lines in C 1 -C 2 space corresponding to the fronts of eluted band profiles under different levels of effective diffusivity should be used to explain the effect of effective diffusivity on the upper side of separation region.
To facilitate the explanation, two solid straight lines are drawn in Figures 1 and 2, which named as L 1 and L 2 , respectively. They are arbitrary lines parallel with the diagonal of m 2 -m 3 plane, but L 2 is  near the diagonal and L 1 far from it. When the SMB is operated along each one of the two lines, m 2 and m 3 increase synchronously and the feeding flow rate keeps constant. Line L 1 is corresponding to the higher feeding flow rate compared with line L 2 . The higher feeding flow rate means more quantity of binary mixture being fed into the columns in the third zone of SMB per minute. Although the same concentration is fed, compared with the lower feeding flow rate along line L 2 , the operating condition of higher feeding flow rate along line L 1 means the stronger nonlinear effect in the columns of SMB which similar with that aroused by high concentration in single column. Therefore, it is easy to understand that in Figure 1, the distances between the left sides of the triangles along line L 1 is much less than that along line L 2. For it has been stated in the discussion about Figure 3 that the increased nonlinearity strongly diminishes the difference between the rears of band profiles eluted from single column under different levels of effect diffusivity. As for the upper sides of the triangles in Figure 1, the remarkable difference between them along line L 1 is nearly same with that along line L 2 except for the slightly upward deviation of the points which obtained under the condition of D efi =1 × 10 -4 cm 2 /min. The reason for the phenomenon is as follows. Just as what mentioned above according to Figure 3, the influence of effective diffusivity is much stronger than that of concentration nonlinearity on the lines in C 1 -C 2 space corresponding to the fronts of eluted band profiles. Therefore, the increase of feed concentration only slightly reduces the difference between the fronts of the eluted band profiles calculated under different levels of effective diffusivity.
A similar interpretation as above can be applied for the effect of external mass transfer coefficient (k ext ) on the separation region of SMB shown in Figure 2.
The stronger influence of concentration nonlinearity on the left sides of separation regions of SMB under different levels of effective diffusivity is an important result for practice. Because in most cases, we hope the objective component of separation is collected in the extract port of SMB. The result above indicates that the purification process of more retained component (C 2 ) is sensitive to either nonideal or nonlinear factors. Therefore, careful determination of the operation condition of SMB is necessary and this is especially important for the separation involves slow intraparticle mass transfer process.

Conclusion
In summary, the main conclusions as follows are obtained. Firstly, the influence of effective diffusivity is more than that of external mass transfer coefficient on the SMB process in practical range of these parameters. Secondly, the separation region of SMB shrinks in the following style with the decrease of intraparticle mass transfer parameters (D eff and k ext ). The upper side of separation region moves downwards evidently and with its slope decreasing slightly; the left side of separation region moves rightward but with its slope decreasing obviously. The behaviors of the separation region of SMB above can be explained by the result in C 1 -C 2 space of binary mixture separation in single column.