Abstract
The simulated moving bed (SMB) technology has attracted considerable attention for its efficiency as a chromatographic adsorptive separation. It has been increasingly applied to the separation of binary mixtures with low separation factors, namely to separate isomers. Although quite a vast amount of information has been published concerning the simulation and design of operating conditions of existing SMB plants, fewer works have addressed the question of design and optimisation of geometric parameters and operating conditions of a new adsorber, especially when mass transfer resistances are significant. The present work extends an algorithm developed elsewhere to design SMB equipment and optimize its operating conditions and applies it to the case of fructoseglucose separation using a cationexchange resin as stationary phase in order to obtain nearly pure fructose in the extract and glucose in the raffinate. The constraints were set as 99% purity for both products. The objective function was chosen to be the adsorbent productivity. The algorithm attempted to find the minimum column lengths for increasing throughputs, which met the required purity constraint. Then, the best construction parameters and operating conditions were chosen as those for which the adsorbent productivity was maximum. The effects of the safety margins applied on the velocity ratios in sections 1 and 4 were examined and a heuristic rule for optimum eluent flowrate was derived. The effect of the purity requirements was also investigated. Finally, the calculated optimal operating points, in terms of flowrate ratios in SMB sections 2 and 3, were analysed in the frame of the equilibrium theory. Sound coherence was verified, which confirmed the accuracy and adequacy of the extended algorithm for the design and optimisation of a SMB adsorber with strong mass transfer effects.
SMB; Design; Optimisation
PROCESS SYSTEM ENGINEERING
Design and optimization of new simulated moving bed plants
D. C. S. AzevedoI, ^{*} * To whom correspondence should be addressed ; A. E. Rodrigues^{II}
^{I}Grupo de Pesquisa em Separações por Adsorção (GPSA), Departamento de Engenharia Química, Universidade Federal do Ceará, Fone: +(55) (85) 40089611, Fax: +(55) (85) 40089610, Campus Universitário do Pici, Bl. 709, CEP 60.455760 Fortaleza  CE, Brasil. Email: diana@gpsa.ufc.br
^{II}Laboratory of Separation and Reaction Engineering (LSRE), Faculdade de Engenharia, University of Porto Rua Dr. Roberto Frias s/n, 4200465 Porto, Portugal
ABSTRACT
The simulated moving bed (SMB) technology has attracted considerable attention for its efficiency as a chromatographic adsorptive separation. It has been increasingly applied to the separation of binary mixtures with low separation factors, namely to separate isomers. Although quite a vast amount of information has been published concerning the simulation and design of operating conditions of existing SMB plants, fewer works have addressed the question of design and optimisation of geometric parameters and operating conditions of a new adsorber, especially when mass transfer resistances are significant. The present work extends an algorithm developed elsewhere to design SMB equipment and optimize its operating conditions and applies it to the case of fructoseglucose separation using a cationexchange resin as stationary phase in order to obtain nearly pure fructose in the extract and glucose in the raffinate. The constraints were set as 99% purity for both products. The objective function was chosen to be the adsorbent productivity. The algorithm attempted to find the minimum column lengths for increasing throughputs, which met the required purity constraint. Then, the best construction parameters and operating conditions were chosen as those for which the adsorbent productivity was maximum. The effects of the safety margins applied on the velocity ratios in sections 1 and 4 were examined and a heuristic rule for optimum eluent flowrate was derived. The effect of the purity requirements was also investigated. Finally, the calculated optimal operating points, in terms of flowrate ratios in SMB sections 2 and 3, were analysed in the frame of the equilibrium theory. Sound coherence was verified, which confirmed the accuracy and adequacy of the extended algorithm for the design and optimisation of a SMB adsorber with strong mass transfer effects.
Keywords: SMB; Design; Optimisation
INTRODUCTION
The simulated moving bed technology has been very popular for the last four decades as a continuous chromatographybased operation to separate binary mixtures. The most commonly used configuration is depicted in Figure 1. It consists of a set of fixed beds connected in a closed circuit with two inlet and two outlet streams. The inlet streams are the mixture to be separated (the load or feed) and a suitable solvent to displace the chemical species being separated (the eluent or desorbent). The outlet streams collect the purified chemical species from the initial mixture: the extract contains the most strongly adsorbed species and the raffinate, the weakly adsorbed one. From time to time (called a switching period, t*), these streams are simultaneously shifted one column ahead in the direction of fluid flow. This mode of operation provides an ingenious way of simulating the countercurrent motion of adsorbent, which enhances the driving force for mass transfer. As a result, relatively inefficient columns may be used in SMB yielding satisfactory separation. Recently, more sophisticated modes of operation have been proposed, such as the Varicol (Subramani et al., 2003), ModiCon (Schramm et al., 2002) and Powerfeed (Zhang et al., 2003), which essentially derive from the conventional SMB just described, providing higher resolution and general superior performance.
Due to the relative complexity of its physical implementation, the definition of the number and dimensions of columns and operating conditions (flowrates and switching time) of a SMB adsorber is not a straightforward task. The flowrates and switching time must be such in order to allow the strongly retained species to concentrate around the extract port and be majorily adsorbed around the raffinate port. The reverse must be true for the weakly adsorbed species: it should concentrate in the bulk liquid phase surrounding the raffinate port.
Since it is a fruitless effort to search for suitable SMB operating conditions by trialand error procedures, the design of SMB adsorbers has attracted significant attention from the academic community. The simplest case may be formulated for systems with linear uncoupled adsorption isotherms. The analysis of the equivalent representation of a true moving bed under an equilibrium model leads to explicit inequality relations between solid and liquid flowrates in the four TMB sections (see Table 1). Nicoud (1992) as well as Ching and Ruthven (1985) changed these inequality relations into design equations by introducing a common safety margin b. Therefore, for linear systems, the SMB system is exclusively a flowcontrolled process. That is to say, its design does not depend on the feed concentration in the frame of equilibrium theory.
Nonlinear behaviour of the isotherms adds to the dependency of the flowrates the feed concentration as an influencing factor. The equilibrium theory applied to nonlinear equilibria was originally addressed by Storti et al (1993) and Mazzotti et al (1996), who proposed explicit relationships to locate the boundaries of the separation region in a m_{2}xm_{3} plane for isotherms of Langmuir and modified Langmuir type. The main finding provided by the equilibrium theory is certainly the explicit definition of the boundaries of separations regions in terms of the solid and fluid flowrates in SMB sections. Furthermore, it is always a useful check and source of comparison for design strategies that incorporate nonideal effects. However, equations from the equilibrium theory refer to flowrate conditions under which a SMB adsorber may achieve 100% purity for both extract and raffinate products, independently on column number and dimensions. In practice, the results are only applicable for columns with a sufficiently high number of theoretical stages (low HETP). Fortunately, for most applications, fairly inefficient columns may safely be used in a SMB operated under the conditions stated by the equilibrium theory (Charton and Nicoud, 1995). Alternatively, adsorbent particle size and column dimensions may also be tailored to overcome such effects as axial mixing and mass transfer resistances, which are main sources of deviation from the results of the equilibrium theory. Nevertheless, some authors have attempted to refine the scope of the equilibrium theory, so as to include those nonideality effects into the design of operating conditions of SMB adsorbers. Pais et al (1997) shaped the regions of separation under nonideal effects by using a detailed model with intraparticle mass transfer being described with a simple LDF approximation. It was shown that the set of values of fluid/solid flowrate ratios in the m_{2}xm_{3} is considerably reduced when mass transfer effects are present, even for a constraint of 99% on product purities. Similar results were obtained by Migliorini et al (1999). By using a complete detailed model, they defined the regions of separation in the m_{2}xm_{3} plane for decreasing purity requirements (99%, 95%, , 60%). For SMB columns having the number of theoretical stages above a threshold value (40), the regions of separation would enlarge with decreasing purity requirements as compared to the region defined by the equilibrium theory. On the other hand, for less efficient columns, the regions of separation will virtually "shrink" in comparison with the ideal region and they may eventually not exist whether the constraint on purity is too strict or the columns have a plate number far below the limiting value. Both of these works attributed values to the flowrate ratios in sections 1 and 4 in accordance with the explicit relations defined in the equilibrium theory with a given safety margin. Azevedo and Rodrigues (1999) observed the impact of varying the values of the constraints for these sections on the resulting separation regions obtained for lowefficient columns. Instead of a 2dimensional parameter space, a three dimensional separation volume is used to present the obtained separation regions as a function of the velocity ratio in section 1.
The theory underlying the equilibrium models is a very useful starting point when it comes to the design of SMB construction parameters. More or less sophisticated process models are used in order to assess the performance criteria on which the optimization and design steps are based. Charton and Nicoud (1995) published a pioneer work on which the authors describe the steps to be followed in designing not only the flowrates but also the construction parameters such as column dimensions and adsorbent particle size. The proposed design procedure makes use of the following pieces of information: adsorption isotherms; optimal flowrate ratios as defined from the equilibrium theory; relationship giving the HETP as a function of mobile phase interstitial velocity (van Deemter); and a relationship for the calculation of pressure drop across packed beds. For a given bead diameter, optimum column length and mobile phase velocity are derived. Following a similar approach, other optimization strategies have been reported (Biressi et al, 2000; LudemannHombourger et al, 2000), which basically elected the adsorbent productivity as an objective function to be maximized. The present work applies a strategy to find design construction parameters and operating conditions required to achieve a desired product purity to the case of fructoseglucose separation using exchange resins. These parameters are then optimised by electing the adsorbent productivity as objective function. The algorithm used is based on the work of Biressi et al (2000); however, instead of using a staged equilibrium model, this work incorporates a detailed steady state process model proposed in a previous publication (Azevedo and Rodrigues, 1999). The model is based on the equivalent representation of a true moving bed (TMB) due to the much lesser computational effort required in its solution as compared to dynamic SMB models. The results are presented for different safety margins b applied on the regeneration sections (sections 1 and 4) so that an optimal value is found. Optimal operating points, in terms of the velocity ratios in sections 2 and 3, are plotted in a parameter space and sound coherence is observed when compared to the predictions from the equilibrium theory.
THEORY UNDERLYING THE DESIGN ALGORITHM
This design strategy essentially determines the minimum column dimensions (length and section) and operating conditions that allow treating a required throughput and achieving a desired purity in the outlet streams without exceeding given pressure drop limits imposed by the packing material. The number of degrees of freedom in the design problem is, therefore seven: column length and diameter, switching time and the four section velocity ratios g_{j}. The algorithm to be applied is an extension of that developed by Biressi et al (2000). However, instead of using a model of equilibrium stages, we have used a detailed TMB steadystate model, as described in a previous publication (Azevedo and Rodrigues, 1999), since this system (fructoseglucose on ionexchange resin) involves considerable mass transfer resistance, which is adequately predicted with a rigorous model. Another novel aspect is the search for the optimal safety margin b to be applied in the velocity ratio constraints in sections 1 and 4, using the adsorbent productivity as objective function.
The flow constraints as stated from the equilibrium model, equations (2), (3) and (4), have been used as general guidelines to provide an initial guess for g_{j}. The other pieces of information to be used in the design/optimisation algorithm are:
A correlation to estimate pressure drop in packed beds, such as the KozenyKárman equation:
A detailed and computationally fast process model, such as that presented by Azevedo and Rodrigues (1999).
According to the guidelines proposed by Biressi et al (2000), the following two assumptions are made:
iA given SMB recycle flowrate (Q1) is assumed and all other flowrates and column cross section are obtained based on that value. After the algorithm calculations are finished, a scale parameter W is calculated as the ratio between the desired feed flowrate and the theoretical feed flowrate, which is that obtained for the recycle flowrate assumed a priori. All other calculated flowrates and section area may be multiplied by the scale factor to obtain the values necessary to process the desired feed flowrate. Calculated g_{j} and L_{c} values remain unchanged.
iiSince pressure drop is proportional to the throughput of the plant, we expect that the productivity will be the highest when pressure drop is the highest possible in the plant. Therefore, pressure drop in section 1, which is where the fluid velocity is maximum, should be equal to the allowable upper limit, hence fixing the value of the u_{1}xL_{c} product from equation (6).
From the adsorption isotherms and using a certain safety margin , values for g1 and g4 are set so as to ensure the proper behaviour of sections 1 and 4, that is, complete regeneration of the adsorbent and eluent, respectively. b is an input parameter to the algorithm so that and Design results for different values of b are compared.
At this point, , g1 and g4 are defined. Three from the initially seven degrees of freedom (g1,g2,g3,g4,t*,L_{c},A) of the design problem are left to be defined: the column length and the velocity ratios g2 and g3. To define the most adequate values for these variables, heuristic rules based on the equilibrium theory are applied having product purities as decision parameters.
The parameter h is defined as the dimensionless distance of a certain (g2,g3) pair to the optimum point given by the equilibrium . In other words:
For each value of , within a certain range defined by the user, the algorithm manages to find the minimum column length and g2, which lead to the required product purities. For a given value of h, the algorithm starts with the minimum g2 value defined by the equilibrium theory and a sufficiently small value for L_{c}. All information required by the process model simulator is available by then and process performance is calculated using a detailed TMB model. Depending on the product purities calculated, the algorithm applies a heuristic rule and makes one of the following decisions:
· Increase column length;
· Increase/decrease g2;
The flow sheet describing the decisionmaking process of the algorithm may be found in Figure 2. The inner loop shows that, for a certain h, the algorithm will stop when the minimum column length that provides the required purities is found. These data will be stored together with the corresponding operating conditions. Plots of h versus column length and adsorbent productivity are then constructed. Adsorbent productivity (kg/hr/m^{3}) is defined as the amount of fructose/glucose obtained per unit time divided by the inventory of adsorbent. The value of h which maximises productivity, for instance, defines the optimal operating point and column length for a given. b. The scale factor corresponding to this optimal condition is also obtained . By multiplying the obtained section area and flowrates by W, one is able to rescale all necessary operating conditions for any desired throughput.
The procedure described in the previous paragraphs (following assumptions i and ii) is repeated for different values of b. An optimum b value (and the respective construction and operating parameters) may be found from the maximum in adsorbent productivity.
RESULTS AND DISCUSSION
This procedure was applied to the separation of fructoseglucose mixtures under the following conditions:

Concentration of each sugar = 400 g/l;

Viscosity = 3 cP;

Specific mass = 1260 kg/m
^{3};

Maximum pressure drop = 2.5 bar/column;

Adsorbent: Dowex Monosphere 99/Ca, d
_{p}= 320 µm;

Linear isotherms with
=0.27 and
= 0.53;

T=50ºC;

12 columns;

Purity constraint = 99% for both extract and raffinate;
The results obtained from the design algorithm will be presented as plots showing the minimum column length and adsorbent productivity as a function of parameter h. This parameter is intrinsically related to the amount of feed being treated, that is, the required throughput. Figure 3 illustrates the results obtained for increasing safety margins b, so that and were assumed in each case. The degree of subdivision was considered, a priori, as 3 columns per section. The axis <PR_{i}> stands for the average adsorbent productivity calculated as follows:
It is clear that an increasing adsorbent inventory becomes necessary as the SMB adsorber is required to treat greater throughputs. This is shown by the rise in the values of L_{c} obtained, which is initially slow and then becomes steep for h>0.8. This may indicate that it is wise to avoid operating a SMB plant close to the theoretical maximum throughput (h=1) because a very large adsorbent inventory would be required to achieve the same purity. This is further confirmed by the fact that the packing is less efficiently used under this condition. In all plots, the average adsorbent productivity reaches a maximum value at h around 0.6 to 0.7 and it declines abruptly thereafter.
In regards to the effect of the safety margin b, the design algorithm results confirm that shorter columns, or smaller adsorbent inventory, are required as this parameter is increased. The adsorbent productivity also tends to increase for higher values of b. The optimal column length was selected for each value of b as that for which productivity is maximum. Figure 4 shows the optimal column lengths as a function of the average adsorbent productivity (<PR_{i}>) and the average solvent consumption (<SC_{i}>). The average solvent consumption is defined as:
The numbers in the plot stand for the respective values of b. The solvent consumption increases steadily as the safety margin is increased with a steep rise for b>1.3. The adsorbent productivity reaches a new peak at the very same value of b. Therefore, for the present case, the parameters calculated for b=1.3 define the optimal operating conditions for a SMB adsorber with 3 columns per section given a required minimum purity of 99%.
Another interesting point is raised when confronting the optimal (g2,g3) pairs obtained for each value of h. Figure 5 shows the path of points in a g2 x g3 parameter space for increasing values of the safety margin b. All "paths" converge to the point , which is the optimal operating point as predicted by the equilibrium theory. However, from the standpoint of adsorbent productivity, the design algorithm using a detailed model indicates that the optimal operating point is located at ^{h}=0.7, rather than at ^{h}=1.
For b=1.3, the optimal geometric parameters, flowrates and switching time are summarized in Table 2, assuming a recycle flowrate of 1000 ml/min in section 1. The rows corresponding to h= 0.3, 0.6 and 0.8 were selected and the respective geometric parameters and operating conditions (t*,g1,g4) were used in order to find the separation area for which both product purities would be higher than 99%. This was done by scanning the equilibrium triangle with a sufficiently fine grid and by performing successive simulations. Figure 6 shows the obtained regions of separation as compared to the optimal (g2,g3) pairs found by the design algorithm. The vertex of each of the separation "triangles" coincides with the corresponding optimal point as calculated by the design algorithm for that value of h.
The same analysis carried out previously was applied to the following SMB column subdivisions: 4332; 3432 and 3342. The purity constraint was kept at 99% and b=1.1. Little deviations were observed by changing the relative length of the sections. Figure 7 shows the results obtained for minimum column length and corresponding productivity, where the lines join the points obtained for the same value of h. Some productivity enhancement is obtained when the size of section 4 is reduced in favour to sections 1 and 3. Yet, the maximum productivity obtained for this configuration is inferior to that obtained for b=1.3 and 3 columns per section.
Finally, the effect of the requirement on product purity was investigated. For a SMB with 3 columns per section and b=1.3, the design algorithm was applied for the following minimum product purities: 99.9%, 99%, 95%, 90% and 80%. Figure 8 shows the obtained plots of the minimum column lengths (a) and respective productivity (b) as a function of h. As expected, the required column length decreases as the constraint on product purity is relaxed and should approach zero as the purity tends to 50% (no separation). Likewise, the adsorbent productivity is enhanced.
CONCLUSIONS
A design package to define construction and operating parameters of a SMB unit was applied to the case of fructoseglucose separation on ionexchange resins, which generally impose significant internal mass transfer resistance. The results were used for optimization purposes. For the linear system under study, it has been demonstrated that SMB units operate inefficiently at the vertex of the equilibrium triangle , which is the optimum operating point in the frame of the equilibrium theory. Maximum productivity is achieved at 0.6<h<0.8. A safety margin of 1.3 was found to be ideal in order to obtain optimal performance in terms of adsorbent productivity. The optimal operating points in a parameter space g_{2}×g_{3} are coherent with the predictions drawn from the equilibrium theory. Changes in the relative size of one section to another do not significantly alter the obtained results. Superior productivity may eventually be obtained by decreasing the size of section 4 in favour of section 1 or 3. For decreasing purity requirements, required column lengths become shorter and the obtained productivity is higher. The bottom line for these trends was found when the required product purity comes to as low as 90%. A lower constraint on purity does not produce either further improvement in productivity or a decrease in column length.
ACKNOWLEDGEMENTS
The authors wish to thank CAPES (Ministry for Education of Brazil) and Fundação da Ciência e Tecnologia (FCTMinistry for Science and Technology of Portugal) for funding this research.
NOMECLATURE
A
Area of a SMB column cross section
(m^{2})d_{p}
Adsorbent particle diameter
(m)K'
Linear adsorption constant for a homogeneous adsorbent particle
(dimensionless) L_{c} Length of a SMB column (m) m_{j} Ratio between liquid and solid flowrates in a TMB section j (dimensionless) Q_{j} Liquid flowrate in a SMB section j (m^{3}/s) T Temperature (ºC) t* Switching time (min) V_{c} Volume of a SMB column (m^{3})Greek Letters
b
Safety factor applied on section velocity ratios (>1)
(dimensionless)gj
Ratio between liquid and solid interstitial velocities in a TMB section j
(dimensionless)DP
Pressure drop across a packed bed
(Pa) e Void fraction between particles in a packed bed (dimensionless) h Dimensionless feed flowrate (dimensionless) µ Absolute viscosity (kg/m/s) u Ratio between solid and fluid volumes (=(1e)/e) (dimensionless) u Fluid superficial velocity (m/s) W Scale factor (dimensionless)Subscripts/Superscripts
1, ,4
SMB/TMB section
()A, B
Strongly and weakly adsorbed chemical species, respectively
()eq
In the frame of the equilibrium theory
() F Feed () FR, GL Fructose and glucose, respectively ()Received: April 15, 2005
Accepted: March 06, 2006
 Azevedo, D. C. S. and Rodrigues, A.E., Bilinear driving force approximation in the modeling of simulated moving bed using bidisperse adsorbents, Ind. Eng. Chem. Res., 38, 35193529 (1999).
 Azevedo, D. C. S. and Rodrigues, A.E., Design of a Simulated Moving Bed in the Presence of MassTransfer Resistances, AIChE Journal, 45, 956966 (1999).
 Azevedo, D. C. S. and Rodrigues, A.E., Design Methodology and Operation of a Simulated Moving Bed Reactor for the Inversion of Sucrose and GlucoseFructose Separation, Chem. Eng. Journal, 82, 95107 (2001).
 Biressi, G., LudemannHombourger, O., Mazzotti, M., Nicoud, R.M. and Morbidelli, M., Design and Optimisation of a SMB Unit: Role of Deviations from Equilibrium Theory, J. Chromatography A, 876, 315 (2000).
 Charton, F. and Nicoud, R.M, Complete Design of a Simulated Moving Bed, J. Chromatography A, 702, 9712 (1995).
 Ching, C. B. and Ruthven, D. M., Experimental Study of a Simulated CounterCurrent Adsorption System  III. Sorbex Operation, Chem. Eng. Sci.,40, 14111417 (1985).
 LudemannHombourger, O., Bailly, M., Nicoud, R.M., Design of a Simulated moving bed: Optimal Particle Size of the Stationary Phase, Sep. Sci. & Technol., 35, 12851305 (2000).
 Mazzotti, M., Storti, G., Morbidelli, M., Robust Design of Countercurrent Adsorption Separation: 3. Nonstoichiometric Systems, AIChE Journal, 42, 27842796 (1996).
 Migliorini, C., Gentilini, A., Mazzotti, M. and Morbidelli, M, Design of Simulated Moving Bed Units under Nonideal Conditions, Ind. Eng. Chem. Res., 38, 24002410 (1999).
 Nicoud, R.M., The Simulated Moving Bed: A Powerful Chromatographic Process, LCGC Intl., 5, 4347 (1992).
 Pais, L. S., Loureiro, J. M. and Rodrigues, A. E., Modeling, Simulation and Operation of a Simulated Moving Bed for Continuous Chromatographic Separation of 1,1'bi2naphtol Enantiomers, J. Chromatography A, 827, 215233 (1997).
 Schramm, H., Kaspareit, M., Kienle, A. and SeidelMorgenstern, A., Improving SMB process by cyclic modulation of the feed concentration, Chem. Eng. Technol. 25, 11511155 (2002).
 Storti, G., Mazzotti, M., Morbidelli, M. and Carrá, S., Robust Design of Binary Countercurrent Adsorption Separation Processes, AIChE J., 39, 471492 (1993).
 Subramani, H. J., Hidajat, K. and Ray, A. K. Optimization of Simulated Moving Bed and Varicol Processes for GlucoseFructose Separation, Trans IchemE, 81, part A, 549567 (2003).
 Zhang, Z., Mazzotti, M. and Morbidelli, M., Powerfeed Operation of Simulated Moving Bed Units: Changing FlowRates during the Switching Interval, J. Chromatography A, 1006, 8799 (2003).
Publication Dates

Publication in this collection
07 July 2006 
Date of issue
June 2006
History

Accepted
06 Mar 2006 
Received
15 Apr 2005