Modelling and simulation of nonlinear chromatographic separation processes: a comparison of di!erent modelling approaches
|
|
- Alban Wheeler
- 5 years ago
- Views:
Transcription
1 Chemical Engineering Science 55 (2000) 373}380 Modelling and simulation of nonlinear chromatographic separation processes: a comparison of di!erent modelling approaches G. DuK nnebier, K.-U. Klatt* Process Control Laboratory, Department of Chemical Engineering,University of Dortmund, Dortmund, Germany Received 31 March 1999; accepted 1 April 1999 Abstract Chromatographic processes provide a powerful tool for the separation of multicomponent mixtures in which the components have di!erent adsorption a$nities. As an alternative to conventional batch chromatography as a discontinuous process, the simulated moving bed (SMB) process as a continuous chromatographic separation process gained more and more impact recently. Mathematical modelling and simulation of batch chromatographic processes is widely documented in the literature. For the modelling of SMB chromatographic processes two di!erent strategies exist. The "rst one considers the corresponding true moving bed process neglecting the cyclic switching and thus the dominating dynmics of the process. The second alternative which is considered within this framework is to connect dynamic models of single chromatographic columns while considering the cyclic port switching in order to represent the dynamics of the real process correctly. We here describe three di!erent approaches for the modelling and simulation of chromatographic separation processes with nonlinear adsorption thermodynamics described by the Langmuir isotherm. The "rst one is a rigorous modelling and numerical simulation approach, the second one is based on the equilibrium theory and proceeds from the analytical solution for binary separation processes based on nonlinear wave propagation. The third approach only includes nonidealities caused by the limited mass transfer rates. All approaches are analyzed and compared with respect to complexity, accuracy and computational time Elsevier Science Ltd. All rights reserved. Keywords: Chromatography; Simulated moving bed; Process modelling; Dynamic simulation; Nonlinear wave theory 1. Introduction Chromatographic processes provide a powerful tool for the separation of multicomponent mixtures in which the components have di!erent adsorption a$nities, especially when the components show separation factors near one and high resolutions, yields and purities are required. Typical applications are found in the pharmaceutical industry and in the food industry where standard thermal unit operations like distillation are not suitable. Conventional batch chromatography is a discontinuous process resulting in highly diluted products. To increase the separating power of the system, a continuous countercurrent operation is desirable but the real countercurrent of the solid phase leads to serious operating problems. In this context, the simulated moving bed (SMB) technology is * Corresponding author. Tel.: # address: karsten@ast.chemietechnik.uni-dortmund.de (K.-U. Klatt) becoming an important technique for large-scale continuous chromatographic separation processes. The SMB process is realized by connecting several single chromatographic columns in series. The countercurrent movement is approximately realized by a cyclic switching of the feed stream and the inlet and outlet ports in the direction of the #uid stream. The main application in the food industry is the sugar separation, especially the separation of fructose and glucose, but also sucrose from molasses, where the SMB technology has growing economic impact (Rearick, Kearney & Costesso, 1997). In the pharmaceutical industry, the separation of enantiomers of chiral drugs on the preparative scale has recently been discovered as an e!ective application for the SMB technology Francotte & Richert, 1997; Schulte, Kinkel, Nicoud & Charton 1996; Strube, AltenhoK ner, Meurer & Schulte, 1997; Pais, Loureiro & Rodrigues, 1997). From the mathematical modelling point of view, it is useful to distinguish chromatographic processes by the type of adsorption isotherms. The type of isotherm /00/$ - see front matter 1999 Elsevier Science Ltd. All rights reserved. PII: S ( 9 9 )
2 374 G. Du( nnebier, K.-U. Klatt / Chemical Engineering Science 55 (2000) 373}380 describing the thermodynamic behavior of the modelled system in#uences the structure of the resulting mathematical problem substantially. Processes with linear adsorption isotherms (DuK nnebier, Weirich & Klatt, 1998; DuK nnebier, Engell, Klatt, Schmidt-Traub, Strube & Weirich, 1998a), as most sugar separation processes, can be modelled and simulated much easier than those with nonlinear isotherms, e.g. competitive Langmuir or Bi-Langmuir isotherms. The separation of enantiomers is a process of the last type, and this paper focusses on di!erent modelling and simulation approaches for nonlinear chromatographic processes. The rigorous dynamic modelling and simulation of chromatographic processes using "nite di!erence, "nite element or orthogonal collocation methods is widely documented in the literature and leads in general to very complex models with an enormous need for computational power (Ching & Lu, 1998; Strube & Schmidt- Traub, 1996; Kaczmarski & Antos, 1996; Lu & Ching, 1997; Ma & Guiochon, 1991). The successful design and operation of SMB separation units depend on the correct choice of the operating conditions, particularly of the #ow rates in each section. Because of the complex dynamics of the process, the choice of the operating parameters is not straightforward. The known design approaches are all based on a steady-state model of the corresponding true moving bed process, e.g. (Ma & Wang, 1997; Mazzotti, Storti & Morbidelli, 1997). Due to the simplicity of the modelling approaches used here, these design strategies only give initial guesses for the optimal operating conditions. For this task, a detailed and reliable dynamic model of the process is necessary which includes the continuous dynamics of the single columns as well as the discrete events resulting from the cyclic operation. Most dynamic modelling and simulation approaches from the literature have only been applied to batch chromatography before, and their application to the SMB process with its switching pattern and the complex initial and boundary conditions is the main and novel scope of this paper. In the sequel, this paper describes the SMB process and the dynamic modelling approach of the process chosen in Sections 2 and 3, respectively. The dynamic modelling approach of the process is based on dynamic models of single chromatographic columns. Those models can be classi"ed by the way of modelling the kinetic and dispersive e!ects which is directly linked to the model complexity. The modelling approaches are presented here in the order of descending complexity. In Section 4, a computationally e$cient solution strategy for a rigorous modelling approach is presented. In Section 5, we show the limitations of a simpli"ed approach lumping all kinetic e!ects in a single parameter. Finally, in Section 6, a solution approach for an ideal model, neglecting all kinetic e!ects, is presented. 2. The simulated moving bed process In most areas of chemical engineering, countercurrent unit operations are preferable for thermodynamic reasons. The standard concept of a countercurrent unit operation in chromatography implies a moving solid phase. However, a moving solid phase is not realizable e!ectively in practice due to technical limitations. An alternative is a simulated moving bed unit which consists of several chromatographic columns connected in series (Ruthven & Ching, 1989). The inlet and outlet ports are located between the columns. The counter current movement of the liquid and solid phases is achieved by sequentially switching the valves of the interconnected columns in the direction of the liquid #ow. In the limit of an in"nite number of columns and in"nitely short switching periods, we obtain the true moving bed process. The process is shown schematically in Fig. 1. The feed and desorbent inlets and extract and ra$nate outlets are switched one column downstream after each switching period. For economic reasons, a small number of columns is desirable. Here we consider a plant with eight columns which can be divided into four sections of two columns each. Each section has a speci"c role within the operation. A section is de"ned by its position relative to the inlets and outlets rather than by a speci"c pair of columns. The stationary regime of this process is a cyclic steady state, in which in each section an identical transient during each period between two valve switches takes place. The cyclic steady state is practically reached after a certain number of valve switches, but the system states are still varying over time because of the periodic movement of the inlet and outlet ports along the columns. In general, this dynamic behavior of the process cannot be approximated by a model of the equivalent countercurrent process. This is especially true for an economically viable SMB process with a rather small Fig. 1. Scheme of a simulated moving bed unit.
3 G. Du( nnebier, K.-U. Klatt / Chemical Engineering Science 55 (2000) 373} number of columns and SMB pocesses with nonlinear adsorption isotherms. We therefore choose a dynamic modelling approach for the purpose of system analysis and process optimization. Our model is assembled of a number of connected dynamic models of single chromatographic columns. The di!erent modelling and simulation approaches investigated are presented below. 3. Modelling the SMB process A dynamic model of the SMB process can be obtained by connecting the dynamic models of the single chromatographic columns while considering the cyclic port switching. The #uid velocities and the inlet concentrations in the di!erent sections can be calculated by mass balances around the inlet and outlet nodes (node model, Ruthven & Ching, 1989). In these equations, Q,Q,Q and Q are the #ow rates through the corresponding process sections, Q is the desorbent #owrate, Q the feed #ow rate, Q the extract #ow rate and Q the ra$nate #ow rate: Desorbent node: Q #Q "Q, (1) c Q "c Q, i"a,b, (2) Extract node: Q!Q "Q, (3) Feednode: Q #Q "Q, (4) c Q #c Q "c Q, i"a,b, (5) Ra$nate node: Q!Q "Q. (6) Dynamic modelling approaches for single chromatographic columns are presented in the following sections. 4. The general rate model Most rigorous modelling approaches found in the literature can be seen as di!erent versions of a general rate model. It takes into account all the important e!ects of nonidealities in the column, which is axial dispersion, pore di!usion and the mass transfer between liquid and solid phase. Only di!erent ways of describing these effects distinguish many modelling approaches. We consider a single chromatographic column with void fraction ε and length through which a #uid mixture containing two di!erent solutes #ows with the velocity u. We assume that the solid phase consists of porous, uniform and spherical particles with void fraction ε and that local equilibrium is established within the pores. Let c represent the concentration of the solute i in the #uid phase and q its concentration in the solid phase, respectively. D is the axial dispersion coe$cient, c the equi- librium concentration, R the particle radius and r the axial coordinate of the particle, k the "lm mass transfer resistance. The di!usion coe$cient within the particle pores is denoted with D, the concentration in the pores of the particle with c. Assuming no radial distribution of u and c, the material balance for each solute component between planes at distances x and x#δx from the entrance of the bed over a time period t to t#δt in the limit yields to the following set of coupled partial di!erential equations: c t #(1!ε)3k c (c!c (r"r ))!D εr x #uc x "0, (7) (1!ε ) εq εt #ε c t!ε D 1 r r rc "0. (8) r The most common nonlinear adsorption equilibrium relation is the competitive Langmuir isotherm with N, k being equilibrium constants: q " N k c 1# k c. (9) This isotherm takes into account both e!ects of saturation and interaction between di!erent components. A numerical solution strategy for this model is given by Gu (1995). To be solved is a coupled PDE system with partial di!erential equations for both #uid and solid phase. In the original work, a "nite element formulation is used to discretize the #uid phase and orthogonal collocation for the solid phase (Finlayson, 1980). The resulting ODEs can be solved using a conventional solver. The use of this model formulation and solution approach for the dynamic simulation of SMB chromatographic processes is a novelty and has not yet been reported. (Ching & Lu, 1998) compared the prediction of a linear driving force (LDF) model with a general rate model for the SMB case and concluded that both approaches are valid for SMB processes with nonlinear adsorption isotherms. The LDF model is based on the assumption of a parabolic concentration pro"le in the particle. Using a single interior point collocation in combination with the symmetry properties of a spheric particle, we can exactly model a parabolic pro"le. We did therefore only use one interior collocation point in our simulation study, resulting in the same number of ODEs to be solved for the more complex general rate model than for the LDF approach. Due to the complexity of the modelling approach, there is a need for simpli"ed modelling and simulation approaches. Furthermore, especially for on-line optimization and process control purposes, the complex models do not allow deep structural insight in the process dynamics.
4 376 G. Du( nnebier, K.-U. Klatt / Chemical Engineering Science 55 (2000) 373} Kinetic approach It has been shown by Van Deemter, Zuiderweg and Klinkenberg (1956), that the nonideal e!ects of axial dispersion and mass transfer resistance and particle di!usion are additive in the case of a linear isotherm under a few assumptions. By adjusting the corresponding parameters, one can therefore simplify the model by considering only one of the e!ects. This relation does not strictly hold in case of a nonlinear isotherm. In case of a Langmuir isotherm it has been shown by Golshan-Shirazi and Guiochon (1992) that, while adjusting the model parameters of several kinetic models to experimental data, the quality of the prediction of simpler models only considering mass transfer resistance does not di!er significantly from those of models with both axial dispersion and mass transfer resistance. The conclusion can be drawn that the assumption of Van Deemter is at least a useful approximation of the reality even in case of a Langmuir Isotherm. Using this simpli"cation, the general rate model reduces to the following set of partial di!erential equations; c t #1!ε ε q t "!uc x, (10) q t "k (c!c ). (11) Tan (1997) proposed a solution approach for this model equations for the case of the multicomponent competitive Langmuir Isotherm (see Eq. (9)) based on a nonlinear transformation combined with an explicit Euler method. The numerical solution of an implicit correlation for just a single reference species is required, which makes the approach very promising with respect to computational cost especially for multicomponent batch chromatography. The main drawback of this approach is the dependency of the extention of the dimensionless domain for which the system has to be solved on the numerical values of the parameters k. The values of these parameters describing the nonideality of the system tend towards in"nity if the system approaches the properties of an ideal system. In this case, the extension of the dimensionless domain tends to in"nity as well. Though the solution of the system for a single gridpoint can be obtained very fast, the number of gridpoints has to be increased in order to ensure stability and desired accuracy. Technical desirable chromatography systems should have a possibly high e$ciency, which is re#ected in relatively large numerical values of the parameter k. The advantages of this solution approach reported, for example systems by Tan (1997), do not hold in this case. This basically re#ects the di$culties commonly encountered when applying classical numerical schemes to hyperbolic partial di!erential equations. To be able to approximate speci"c features of the solutions (e.g. dicontinuous solutions, shock waves) using a conventional discretization scheme a very "ne or automatically adapted grid is needed. This basic limitation can not be circumvented by the approach presented by Tan. A solution approach especially developed for hyperbolic PDEs is presented in the following section. 6. The ideal model Neglecting all nonideal e!ects and assuming overall equilibrium between liquid and solid phase, the ideal model of chromatography is obtained. It is usually formulated after de"ning a set of dimensionless variables with being a characteristic length of the column: x"z/, τ"εut/. (12) A di!erential material balance for all components leads in the limit to the following set of quasilinear partial di!erential and algebraic equations: c x #f "0, (13) τ εc #(1!ε)q "f. (14) With the Langmuir isotherm as in Eq. (9), this gives a system of coupled hyperbolic "rst-order partial di!erential equations which can either be solved numerically or the system might be solved in closed form for some cases of initial and boundary conditions using the theory of nonlinear wave propagation. The details of this solution procedure have already been presented by DuK n- nebier and Klatt (1998) and are brie#y summarized in the sequel Summary of solution procedure The solution procedure can be summarized as follows: In the xrst step, the system is transformed with a nonlinear transformation (h-transformation) to decouple the equations. The solution for a Riemann problem (step input and constant initial conditions) can then be found analytically using the method of characteristics (Hel-!erich & Whitley, 1996; Hel!erich, 1997; Rhee, Aris & Amundsen, 1989,1970). This methodology can, in its standard form, only be applied to systems with competitive Langmuir adsorption isotherms (see Eq. (9)), but an extension to more complex isotherms can be found in (ZenhaK usern & Rippin, 1998). The solutions of these problems are travelling waves, and with the theory of wave interaction the solution of several Riemann problems can be superposed. The straightforward application of this approach is limited, since more complex
5 G. Du( nnebier, K.-U. Klatt / Chemical Engineering Science 55 (2000) 373} interaction patterns might require the solution of a partial di!erential equation themselves. Since the usual initial and boundary conditions in most applications as the SMB processes are more complex than one single step, in the second step, the initial and boundary conditions have to be discretized as a series of steps using front tracking (Wendro!, 1992). Assuming a su$cient discretization, the whole solution can now be calculated analytically in the transformed space using the rules of wave propagation and wave interaction of shock waves. In the third step, the solution has to be transformed back into the original concentration space. The approach has previously been only applied to batch chromatography by (Wendro! (1993). The implementation to study a simulated moving bed process as a novel application is a challenging task due to the far more complex initial and boundary conditions in this case. 7. Simulation example for SMB processes In this simulation example we consider a SMB unit with one column each in the "rst and fourth section and three columns in the second and third section, which re#ects an industrial separation of cyclic hydrocarbons. The adsorption equilibrium of the system can be described with a competitive Langmuir Isotherm as given in Eq. (9). The purpose of this example is to illustrate the capabilities of the modelling and simulation approaches described above. The simulation results of the general rate model are compared with the predictions of the ideal model, simulated using the method of characteristics and front tracking. For the reasons presented in Section 5, the kinetic approach has not been implemented for simulation. The isotherm parameters N have to be adjusted for the two modelling approaches. In the rigorous model, the adsorbed solid concentration q is expressed per volume of solid. In the ideal model, the adsorbed solid concentration q is expressed per particle volume. With q "q (1!ε )#c ε (15) the following relation holds: N"N(1!ε )# ε k k c. (16) The system parameters of the example system chosen are N"0.8763, N"0.6812, k "1.85 cm/g, k "5.94 c m/g, all columns have a void fraction of ε"0.4, a diameter of d"5 cm, a length of "25 cm. The parameters N have been estimated for the interesting concentration range using Eq. (16). The switching is performed after each t "3500 s, and the feed and outlet streams are chosen as follows: Q " cm/s with a concentration of c "0.23 g/cm, Q "0.828 cm/s, Q " cm/s and Q " cm/s General rate model For the discretization of the #uid phase, 12 cubic "nite elements are used for each column of Length "25 cm. This results in computation times of approx. 30 cpu seconds on a PC PII266 for one switching period Ideal model The sampling parameter for the front tracking is chosen as 0.01 (see DuK nnebier & Klatt, 1998). This results in computation times of approx. 12 cpu seconds on a PC PII266 for one switching period Comparison of simulation results In Fig. 2, the axial pro"les of the central columns of the SMB unit at di!erent instances during the startup period are displayed. Starting with at the top of the "gure, a comparison of the axial pro"les at the end of the "rst switching period leads to the following conclusions: The solution of the ideal model forms a kind of skeleton for the prediction of the general rate model. If the equilibrium is reached despite the kinetic e!ects (as in the "rst column after the feed inlet), both pro"les are nearly identical. This is only the case if the extension of the concentration plateaus is su$ciently large. If the equilibrium is not reached due to the kinetic e!ects (as for component B at the frontal adsorption wave at the beginning of the second column after the feed inlet), the plateau level predicted by the two models is di!erent. Since the concentration level a!ects the speed of propagation in case of a nonlinear isotherm, the speed of propagation of the concentration waves is not predicted correctly by the ideal model in this case. This is a speci"c feature of the nonlinear adsorption isotherm, a very good agreement between the ideal and the general rate model can be obtained in case of a linear adsorption isotherm (DuK nnebier et al., 1998a,b). This deviation is particulary large in Section III (between feed inlet and ra$nate outlet) and grows during the simulation due to the integration, as can be seen in the two lower parts of the "gure showing the axial pro"les at the end of the second and "fth switching periods, respectively. The axial pro"le of the column in cyclic steady state, which is approximately reached after the 70th switching period, as predicted by both modelling approaches is
6 378 G. Du( nnebier, K.-U. Klatt / Chemical Engineering Science 55 (2000) 373}380 Fig. 2. Development of axial pro"les during startup procedure. displayed in Fig. 3. One can consider the solution of the ideal model as a skeleton for the solution of the general rate model, and the predicted pro"les do agree closely in most parts of the process.the small error in the calculated speed of propagation of the concentration waves leads to an error in the prediction of the pro"les in steady state due to the integrating e!ect of the long time horizon. The concentration levels of the constant states are slightly di!erent, and some e!ects, like the kinks predicted by the ideal model in sections two and three, are predicted more smoothly by the general rate model due to the e!ects of mass transfer resistance and axial dispersion. The predicted location of the desorption front of component A at the extract outlet and the adsorption front of component B at the ra$nate outlet of both models show good agreement and the chosen operating conditions lead to the desired separation in both cases. As a "rst initial guess, a design of a SMB unit using the ideal model might therefore be feasible. But the di!erent prediction of the product #ows of the impurities leads to signi"cantly di!erent product purities in both cases: Purity [%] Extract Ra$nate Ideal model General rate model Due to the unrealistic assumptions of ideality, the accuracy of prediction of this modelling approach is not su$cient for design, control or optimization purposes. Those deviations between the two models are not due to the front tracking procedure, since the same e!ect can also be found in the case of batch chromatography, where analytical solutions for the ideal model are available. Great care has to be taken when using design approaches based on ideal models, as the one by
7 G. Du( nnebier, K.-U. Klatt / Chemical Engineering Science 55 (2000) 373} Fig. 3. Comparison of axial pro"les of SMB unit for cyclic steady state at the end of a switching period. Mazzotti et al. (1997). Besides the questionable assumption of ideality, these approaches use a (steady-state) model of the corresponding true moving bed process introducing an additional error. The physical plane portrait (showing the propagation directions of the di!erent waves in a τ vs. x graph) as one of the results of the ideal model solution might be used to obtain an understanding of the complex and interacting process dynamics. By following the shock paths through space and time, one can examine the complex interaction patterns. An example is given in DuK nnebier and Klatt (1998). Therefore, the ideal model can be useful for the analysis of the qualitative behavior of the process and the in#uence of di!erent process parameters. Especially the possible insight in the dynamics which can be gained by examination of the shock paths helps to understand and analyse the process dynamics and the startup procedure. The ideal model should not be chosen in any application where precise quantitative predictions are needed since those can not be given due to the unrealistic modelling assumptions. 8. Conclusions Three di!erent modelling and simulation approaches for chromatographic separation process have been presented and analyzed, starting with the rigorous general rate model considering various kinetics e!ects. Two simpli"ed modelling approaches have been presented in the sequel, one lumping the kinetic e!ects into a single parameter, the second (being the ideal model) neglecting all kinetic nonidealities. The suitable modelling and simulation approach has to be chosen depending on the speci"c application. This paper presented a case study extending and applying two modelling and simulation approaches to the complex simulated moving bed (SMB) process as a novel application. The case study was used to illustrate the capabilities of the proposed approaches. The ideal model and the solution approach based on the method of characteristics and front tracking allows speci"c qualitative insight into the complex interacting dynamics of this process especially during the startup procedure but lacks of a su$cient accuracy for quantitative predictions due to the unrealistic modelling assumptions. The approach proposed based on a kinetic model seems to be very elegant and e!ective at "rst glance but is limited to columns with low e$ciency which are very rarely used in practical applications today. The general rate model and the numerical solution approach presented here allow for a quantitative analysis of the process. As to be expected, those calculations tend to be more computational expensive than the solution of the ideal model, though far less as reported otherwise in the literature. Acknowledgements The authors are indebted to S. Engell for his support and valuable input to this work. The "nancial support of the Bundesministerium fuk r Bildung und Forschung under the grant number 03D0062B0 is very gratefully acknowledged. We also kindly acknowledge the help and support of T. Gu. Notation Variables c c c D concentration of component i in bulk phase, g/cm equilibrium concentration of component i, g/cm concentration of component i in particle pores, g/cm axial dispersion coe$cient, cm/s
8 380 G. Du( nnebier, K.-U. Klatt / Chemical Engineering Science 55 (2000) 373}380 D particle di!usivity coe$cient of component i, cm/s ε void fraction of the bed, dimensionless ε void fraction of the particle, dimensionless k "lm mass transfer resistance of component i, cm/s k lumped mass transfer resistance of component i, 1/s k Langmuir Isotherm coe$cient of component i, cm/g column length, cm N Langmuir isotherm constant, g/cm q concentration of component i in the solid phase, g/cm Q volumetric #ow rate, cm/s R particle radius, cm τ dimensionless time, dimensionless t valve switching interval for SMB process, s u #uid velocity, cm/s x axial coordinate, cm Z dimensionless axial coordinate, dimensionless Indices i in j out R E F D References component i"a,b column inlet section j"i,ii,iii,iv column outlet ra f f inate extract feed desorbate Ching, C., & Lu, Z. (1998). Parabolic intraparticle concentration pro"le assumption in modeling and simulation of nonlinear simulated moving-bed separation processes. Chemical Engineering Science, 53, 1311}1315. DuK nnebier, G., Engell, S., Klatt, K. -U., Schmidt-Traub, H., Strube, J., & Weirich, I. (1998a). Modelling of simulated moving bed chromatographic processes with regard to process control design. Computers in Chemical Engineering, 22, S855}S858. DuK nnebier, G., & Klatt, K. (1998). Modelling of chromatographic separation processes using nonlinear wave theory. In C. Georgakis, DYCOPS-5 xfth IFAC symposium on dynamics and control of process systems, Corfu, Greece (pp. 521}526). DuK nnebier, G., Weirich, I., & Klatt, K.-U. (1998). Computationally e$cient dynamic modelling and simulation of simulated moving bed chromatographic processes with linear isotherms. Chemical Engineering Science, 53, 2537}2546. Finlayson, B. (1980). Nonlinear analysis in chemical engineering. New York: McGraw-Hill. Francotte, E., & Richert, P. (1997). Applications of simulated movingbed chromatography to the separation of the enantiomers of chiral drugs. Journal of Chromatography A, 769, 101}107. Golshan-Shirazi, S., & Guiochon, G. (1992). Comparison of the various kinetic models of non-linear chromatography. Journal of Chromatography A, 603, 1}11. Gu, T. (1995). Mathematical modelling and scale up of liquid chromatography. New York: Springer. Hel!erich, F. (1997). Non-linear waves in chromatography III. Multicomponent langmuir and langmuir-like systems. Journal of Chromatography A, 768, 169}205. Hel!erich, F., & Whitley, R. (1996). Non-linear waves in chromatography II. Wave interference and coherence in multicomponent systems. Journal of Chromatography A, 734, 7}47. Kaczmarski, K., & Antos, D. (1996). Fast "nite di!erence method for solving multicomponent adsorption-chromatography models. Computers in Chemical Engineering, 20(11), 1271}1276. Lu, Z., & Ching, C. (1997). Dynamics of simulated moving-bed adsorption separation processes. Separation Science and Technology, 32, 1993}2010. Ma, Z., & Guiochon, G. (1991). Application of orthogonal collocation on "nite elements in the simulation of non-linear chromatography. Computers in Chemical Engineering, 15(6), 415}426. Ma, Z., & Wang, N.-H. (1997). Standing wave analysis of SMBchromatography: Linear systems. A.I.Ch.E. Journal, 43, 2488}2508. Mazzotti, M., Storti, G., & Morbidelli, M. (1997). Optimal operation of simulated moving bed units for nonlinear chromatographic separations. Journal of Chromatography A, 769, 3}24. Pais, L., Loureiro, J., & Rodrigues, A. (1997). Modeling, simulation and operation of a simulated moving bed for continuos chromatographic separation of 1,1'-bi-2-naphtol enantiomers. Journal of Chromatography A, 769, 25}35. Rearick, D., Kearney, M., & Costesso, D. (1997). Simulated moving-bed technology in the sweetener industry. Chemtech, 27(9), 36}40. Rhee, H., Aris, R., & Amundsen, N. (1989). First order partial diwerential equations, volume II, Theory and application of hyperbolic systems of quasilinear equations. New York: Prentice Hall. Rhee, H.-K., Aris, R., & Amundsen, N. (1970). On the theory of multicomponent chromatography. Philosophical Transactions of the Royal Society of London A, 267, 419}454. Ruthven, D., & Ching, C. (1989). Counter-current and simulated moving bed adsorption separation processes. Chemical Engineering Science, 44, 1011}1038. Schulte, M., Kinkel, J., Nicoud, R., & Charton, F. (1996). Simulierte Gegenstromchromatographie } eine e$ziente Technik zur Herstellung optisch aktiver Verbindungen im industriellen Ma{stab. Chemie Ingenieur Technik, 68, 670}683. Strube, J., AltenhoK ner, U., Meurer, M., Schmidt-Traub, H., & Schulte, M. (1997). Dynamic simulation of simulated moving bed chromatographic processes for the optimization of chiral separations. Journal of Chromatography A, 769, 81}92. Strube, J., & Schmidt-Traub, H. (1996). Dynamic simulation of simulated-moving-bed chromatographic processes. Computers Chemical Engineering, 20(suppl), 641}646. Tan, H. (1997). New algorithms for the computation of column dynamics of multicomponent liquid phase adsorption. Adsorption, 3, 137}149. Van Deemter, J., Zuiderweg, F., & Klinkenberg, A. (1956). Longitudinal di!usion and resistance to mass transfer as causes of nonideality in chromatography. Chemical Engineering Science, 5, 271}280. Wendro!, B. (1992). Composite #ows. Physica D, 60, 208}215. Wendro!, B. (1993). An analysis of front tracking for chromatography. Acta Applicandae Mathematicae, 30, 265}285. ZenhaK usern, R., & Rippin, D. (1998). Modelling and simulation of multicomponent nonlinear chromatography. Computers Chemical Engineering, 22, 259}281.
MODEL-BASED OPTIMIZATION AND CONTROL OF CHROMATOGRAPHIC PROCESSES
MODEL-BASED OPTIMIZATION AND CONTROL OF CHROMATOGRAPHIC PROCESSES Karsten-Ulrich Klatt Λ;1 Felix Hanisch Λ Guido Dünnebier Λ Sebastian Engell Λ Λ Process Control Laboratory, Department of Chemical Engineering,
More informationSUPERCRITICAL FLUID CHROMATOGRAPHY PROCESS OPTIMISATION OF THE SEPARATION OF TOCOPHEROL HOMOLOGUES
SUPERCRITICAL FLUID CHROMATOGRAPHY PROCESS OPTIMISATION OF THE SEPARATION OF TOCOPHEROL HOMOLOGUES S. Peper, S. Cammerer, M. Johannsen, G. Brunner Technical University Hamburg-Harburg, Thermal Separation
More informationNOTICE: this is the author s version of a work that was accepted for publication in Computers and Chemical Engineering. Changes resulting from the
NOTICE: this is the author s version of a work that was accepted for publication in Computers and Chemical Engineering. Changes resulting from the publishing process, such as peer review, editing, corrections,
More informationChiral Separation of Propranolol Hydrochloride Through an SMB Process Integrated with Crystallization
J. Ind. Eng. Chem., Vol. 12, No. 6, (2006) 868-876 Chiral Separation of Propranolol Hydrochloride Through an SMB Process Integrated with Crystallization Xin Wang, Yue Liu, Hong Wei Yu, and Chi Bun Ching
More informationPreparative Chromatography
Edited by Henner Andreas Seidel-Morgenstern Michael Schulte, and Preparative Chromatography Completely Revised and Updated Edition VCH Verlag GmbH & Co. KGaA About List of List of Notations XIX XXI 1 Introduction
More informationHeat and Mass Transfer over Cooled Horizontal Tubes 333 x-component of the velocity: y2 u = g sin x y : (4) r 2 The y-component of the velocity eld is
Scientia Iranica, Vol., No. 4, pp 332{338 c Sharif University of Technology, October 2004 Vapor Absorption into Liquid Films Flowing over a Column of Cooled Horizontal Tubes B. Farhanieh and F. Babadi
More informationSIMULATION ASSESSMENT OF CONTINUOUS SIMULATING MOVING BED CHROMATOGRAPHY PROCESS WITH PARTIAL FEED AND NEW STRATEGY WITH PARTIAL FEED
Brazilian Journal of Chemical Engineering ISSN 14-6632 Printed in Brazil www.abeq.org.br/bjche Vol. 26, No. 3, pp. 595-61, July - September, 29 SIMULATION ASSESSMENT OF CONTINUOUS SIMULATING MOVING BED
More informationControl of incomplete separation in simulated moving bed chromatographic processes
Preprint, 11th IFAC Symposium on Dynamics and Control of Process Systems, including Biosystems Control of incomplete separation in simulated moving bed chromatographic processes Paul Suvarov*, Alain Vande
More informationMODELING OF INDUSTRIAL PSA PROCESS FOR FUEL ETHANOL PRODUCTION. Abstract
MODELING OF INDUSTRIAL PSA PROCESS FOR FUEL ETHANOL PRODUCTION Marian Simo, University at Buffalo, Buffalo, NY, USA Christopher J. Brown, Thermal Kinetics Systems, Buffalo, NY, USA and Vladimir Hlavacek,
More informationKonrad-Zuse-Zentrum f r Informationstechnik Berlin
Konrad-Zuse-Zentrum f r Informationstechnik Berlin J. Frauhammer 1 H. Klein 2 G. Eigenberger 1 U. Nowak Solving moving boundary problems with an adaptive moving grid method: Rotary heat exchangers with
More informationStart-up and wrong-way behavior in a tubular reactor: dilution e!ect of the catalytic bed
Chemical Engineering Science 55 (2000) 3885}3897 Start-up and wrong-way behavior in a tubular reactor: dilution e!ect of the catalytic bed Margarida M. J. Quina, Rosa M. Quinta Ferreira* Department of
More informationAmmonia Gas Absorption
Ammonia Gas Absorption by Oscar D. Crisalle Professor Chemical Engineering Department University of Florida crisalle@che.ufl.edu Revision 12: September 24, 2013 Ammonia Absorption Rev 08-04/15/2013 Page
More informationModeling and Simulation of Fixed bed Adsorption column: Effect of Velocity Variation
Modeling and Simulation of Fixed bed Adsorption column: Effect of Velocity Variation B.V. Babu * and Suresh Gupta Department of Chemical Engineering Birla Institute of Technology & Science Pilani 333031
More informationTheoretical model of electrostatic precipitator performance for collecting polydisperse particles
Journal of Electrostatics 50 (2001) 177}190 Theoretical model of electrostatic precipitator performance for collecting polydisperse particles S.H. Kim, H.S. Park, K.W. Lee* Department of Environmental
More informationSimulated Moving Bed Chromatography (SMB) F DIN EN ISO 9001 E D Q U A L T N E T Y M A N E M A G SMB. The next step in chromatography
Simulated Moving Bed Chromatography (SMB) C E R T I F DIN EN ISO 9001 Q U A L I T Y M A N A G I E D E M T N E SMB The next step in chromatography Simulated Moving Bed Chromatography (SMB) SMB Pilot Unit
More informationUniversity of Groningen
University of Groningen Modelling of gas-liquid reactors - stability and dynamic behaviour of a hydroformylation reactor, influence of mass transfer in the kinetics controlled regime Elk, E.P. van; Borman,
More informationIsopropyl Alcohol Vapour Removal from Diluted Gaseous Stream by Adsorption: Experimental Results and Dynamic Model
451 A publication of CHEMICAL ENGINEERING TRANSACTIONS VOL. 47, 16 Guest Editors: Angelo Chianese, Luca Di Palma, Elisabetta Petrucci, Marco Stoller Copyright 16, AIDIC Servizi S.r.l., ISBN 978-88-9568-38-9;
More informationMathematical modeling of pulp washing on rotary drums and their numerical solution for various adsorption isotherms
ISSN 746-7233, England, UK World Journal of Modelling and Simulation Vol. 6 (2) No. 3, pp. 24-222 Mathematical modeling of pulp washing on rotary drums and their numerical solution for various adsorption
More informationUniversity of Groningen
University of Groningen Modelling of gas-liquid reactors - stability and dynamic behaviour of gas-liquid mass transfer accompanied by irreversible reaction Elk, E.P. van; Borman, P.C.; Kuipers, J.A.M.;
More informationn i,j+1/2 q i,j * qi+1,j * S i+1/2,j
Helsinki University of Technology CFD-group/ The Laboratory of Applied Thermodynamics MEMO No CFD/TERMO-5-97 DATE: December 9,997 TITLE A comparison of complete vs. simplied viscous terms in boundary layer
More informationDirect numerical simulation of compressible multi-phase ow with a pressure-based method
Seventh International Conference on Computational Fluid Dynamics (ICCFD7), Big Island, Hawaii, July 9-13, 2012 ICCFD7-2902 Direct numerical simulation of compressible multi-phase ow with a pressure-based
More informationBuoyancy-induced #ows and phase-change heat transfer in a vertical capillary structure with symmetric heating
Chemical Engineering Science 55 (2000) 2653}2661 Buoyancy-induced #ows and phase-change heat transfer in a vertical capillary structure with symmetric heating T. S. Zhao*, P. Cheng, C. Y. Wang Department
More informationTHE IRANIAN JAM PETROCHEMICAL S H 2 -PSA ENHANCEMENT USING A NEW STEPS SEQUENCE TABLE
Petroleum & Coal ISSN 1337-707 Available online at www.vurup.sk/petroleum-coal Petroleum & Coal 56 (1) 13-18, 014 THE IRANIAN JAM PETROCHEMICAL S H -PSA ENHANCEMENT USING A NEW STEPS SEQUENCE TABLE Ehsan
More informationComparison of Adsorbent Behavior in Glucose/Fructose Separation by Simulated Moving Bed (SMB) Chromatography
University of Twente Faculty of Chemical Technology Comparison of Adsorbent Behavior in Glucose/Fructose Separation by Simulated Moving Bed (SMB) Chromatography Authors: J. Blignaut K. Albataineh Y. Banat
More informationAn analysis of crystal dissolution fronts in ows through porous. media. Part 1: Compatible boundary conditions. P. Knabner
An analysis of crystal dissolution fronts in ows through porous media Part 1: Compatible boundary conditions 1 2 P. Knabner University of Erlangen{Nurnberg, Institute of Applied Mathematics, Martensstrae
More informationSEPARATION OF CITRIC AND LACTIC ACID FROM FERMENTATION LIQUID BY THE CHROMATOGRAPHIC METHOD MODELING AND PARAMETERS ESTIMATION
SEPARATION OF CITRIC AND LACTIC ACID FROM FERMENTATION LIQUID BY THE CHROMATOGRAPHIC METHOD MODELING AND PARAMETERS ESTIMATION Paweł GLUSZCZ and Jerzy PETERA Faculty of Process and Environmental Engineering,
More informationDesign of simulated moving bed and Varicol processes for preparative separations with a low number of columns
Journal of Chromatography A, 1006 (2003) 33 44 www.elsevier.com/ locate/ chroma Design of simulated moving bed and Varicol processes for preparative separations with a low number of columns Luıs S. Pais,
More informationA New Mass Transfer Model for Cyclic Adsorption and Desorption
Korean J. Chem. Eng., 63), 40-405 999) SHORT COMMUNICATION A New Mass Transfer Model for Cyclic Adsorption and Desorption Yong-Min Kim* and Sung-Sup Suh Department of Chemical Engineering, Hong-Ik University,
More information6.1 Steady, One-Dimensional Rectilinear Flows Steady, Spherically Symmetric Radial Flows 42
Contents 6 UNIDIRECTIONAL FLOWS 1 6.1 Steady, One-Dimensional Rectilinear Flows 6. Steady, Axisymmetric Rectilinear Flows 19 6.3 Steady, Axisymmetric Torsional Flows 8 6.4 Steady, Axisymmetric Radial Flows
More informationStudy of a Simulated Moving Bed (SMB) Preparative Liquid Chromatographic Appliance and Technology
Study of a Simulated Moving Bed (SMB) Preparative Liquid Chromatographic Appliance and Technology 1 Melinda Nagy, 1 Zoltán Molnár, 1 László Hanák, 1 János Argyelán, 1 Tibor Szánya, 2 Antal Aranyi, 2 Krisztina
More informationModeling chromatographic columns Non-equilibrium packed-bed adsorption with non-linear adsorption isotherms
Journal of Chromatography A, 1041 (2004) 77 85 Modeling chromatographic columns Non-equilibrium packed-bed adsorption with non-linear adsorption isotherms Ahmet R. Özdural a,, Aslı Alkan a, Piet J.A.M.
More informationCalibration and Optimization of Hydrophobic Interaction Chromatography
Calibration and Optimization of Hydrophobic Interaction Chromatography Abstract Alex Olsson Department of Chemical Engineering, Lund University, Sweden 2012-11-01 One of the most commonly used method of
More informationOptimal Synthesis and Design of Advanced Chromatographic Process Concepts
Optimal Synthesis and Design of Advanced Chromatographic Process Concepts Habilitationsschrift von Dr.-Ing. Malte Kaspereit geboren am 31. Juli 1973 in Merseburg zur Verleihung des akademischen Grades
More informationInlet Grid Block! Standard Grid Blocks Outlet Grid Point r Center of Rotation Equilibrium Grid Point Outlet Grid Block Figure 1 The Finite Dierence Gr
THE CONSTANTLY-ACCELERATING CENTRIFUGE REVISITED by Douglas Ruth University ofmanitoba Winnipeg, Manitoba, Canada Abstract SCA9812 In 1963, Homan described the use of a constantly-accelerating centrifuge
More informationLecture 7. Sorption-Separation Equipment
Lecture 7. Sorption-Separation Equipment Adsorption - Stirred-tank, slurry operation - Cyclic fixed-bed batch operation - Thermal (temperature)-swing adsorption - Fluidizing bed for adsorption and moving
More informationExplicit approximations to estimate the perturbative diusivity in the presence of convectivity and damping (Part 1): Semi-innite slab approximations
Explicit approximations to estimate the perturbative diusivity in the presence of convectivity and damping (Part 1): Semi-innite slab approximations M. van Berkel 1,2,3,4, H.J. Zwart 5,6, G.M.D. Hogeweij
More informationSELECTIVE REMOVAL OF CARBON DIOXIDE FROM AQUEOUS AMMONIA SOLUTIONS
Distillation Absorption 2010 A.B. de Haan, H. Kooijman and A. Górak (Editors) All rights reserved by authors as per DA2010 copyright notice SELECTIVE REMOVAL OF CARBON DIOXIDE FROM AQUEOUS AMMONIA SOLUTIONS
More informationDynamics Of Double Pipe Heat Exchangers: Explicit Time Domain Solutions
Dynamics Of Double Pipe Heat Exchangers: Explicit Time Domain Solutions Franco Evangelista Department of Chemistry, Chemical Engineering and Materials University of L Aquila Italy The dynamics of double
More informationThe use of cubic spline and far-side boundary condition for the collocation solution of a transient convection-diffusion problem
Korean J Chem Eng, 24(2), 204-208 (2007) SHORT COMMUNICATION The use of cubic spline and far-side boundary condition for the collocation solution of a transient convection-diffusion problem Woohyun Yun
More informationresults for temperatures, heat uxes, and Nusselt skin subjected to this localized thermal load. cooling are examined in a parametric study that
Abstract A computational study of the eectiveness of tangential convective cooling for the protection of structures from intense, localized thermal loads is presented. Two-dimensional numerical models
More informationquantum semiconductor devices. The model is valid to all orders of h and to rst order in the classical potential energy. The smooth QHD equations have
Numerical Simulation of the Smooth Quantum Hydrodynamic Model for Semiconductor Devices Carl L. Gardner and Christian Ringhofer y Department of Mathematics Arizona State University Tempe, AZ 8587-84 Abstract
More informationOPERATOR SPLITTING METHODS FOR SYSTEMS OF CONVECTION-DIFFUSION EQUATIONS: NONLINEAR ERROR MECHANISMS AND CORRECTION STRATEGIES
OPERATOR SPLITTING METHODS FOR SYSTEMS OF CONVECTION-DIFFUSION EQUATIONS: NONLINEAR ERROR MECHANISMS AND CORRECTION STRATEGIES K. HVISTENDAHL KARLSEN a, K.{A. LIE b;c, J. R. NATVIG c, H. F. NORDHAUG a,
More informationDextran and fructose separation on an SMB continuous chromatographic unit
Biochemical Engineering Journal 12 (2002) 215 221 Dextran and fructose separation on an SMB continuous chromatographic unit Mariana S. Coelho a,1, Diana C.S. Azevedo a,2, José A. Teixeira b,alírio Rodrigues
More informationVALIDATION OF AN ACTUAL SIMULATED MOVING BED UNIT TOWARDS ITS AUTOMATION CONTROL AIMING TO SEPARATE PRAZIQUANTEL ENANTIOMERS
VALIDATION OF AN ACTUAL SIMULATED MOVING BED UNIT TOWARDS ITS AUTOMATION CONTROL AIMING TO SEPARATE PRAZIQUANTEL ENANTIOMERS F. C. Cunha 1,a ; A. G. Barreto Jr 2 ; M. B. de Souza Jr 2 ; A. R. Secchi 1
More informationOptimal startup operation of simulated moving bed chromatographic processes
9th International Symposium on Dynamics and Control of Process Systems (DYCOPS 2010) Optimal startup operation of simulated moving bed chromatographic processes Suzhou Li Yoshiaki Kawairi Jörg Raisch,
More informationSOE3213/4: CFD Lecture 3
CFD { SOE323/4: CFD Lecture 3 @u x @t @u y @t @u z @t r:u = 0 () + r:(uu x ) = + r:(uu y ) = + r:(uu z ) = @x @y @z + r 2 u x (2) + r 2 u y (3) + r 2 u z (4) Transport equation form, with source @x Two
More informationA modication of the lm. model. Heikki Haario. Department of Mathematics. University of Helsinki. Hallituskatu 15, Helsinki.
A modication of the lm model Heikki Haario Department of Mathematics University of Helsinki Hallituskatu 15, 00100 Helsinki Finland e-mail: hhaario@convex.csc.i Thomas I. Seidman Department of Mathematics
More informationGravitational potential energy *
OpenStax-CNX module: m15090 1 Gravitational potential energy * Sunil Kumar Singh This work is produced by OpenStax-CNX and licensed under the Creative Commons Attribution License 2.0 The concept of potential
More informationRidge analysis of mixture response surfaces
Statistics & Probability Letters 48 (2000) 3 40 Ridge analysis of mixture response surfaces Norman R. Draper a;, Friedrich Pukelsheim b a Department of Statistics, University of Wisconsin, 20 West Dayton
More informationAnalytical solution for facilitated transport across a membrane
Chemical Engineering Science 57 () 4817 489 www.elsevier.com/locate/ces Analytical solution for facilitated transport across a membrane Mohamed Hassan Al-Marzouqi a, ees J. A. Hogendoorn b, Geert F. Versteeg
More informationUncertainty in radon measurements with CR39 detector due to unknown deposition of Po
Nuclear Instruments and Methods in Physics Research A 450 (2000) 568} 572 Uncertainty in radon measurements with CR39 detector due to unknown deposition of Po D. NikezicH, K.N. Yu* Department of Physics
More informationMultigrid Approaches to Non-linear Diffusion Problems on Unstructured Meshes
NASA/CR-2001-210660 ICASE Report No. 2001-3 Multigrid Approaches to Non-linear Diffusion Problems on Unstructured Meshes Dimitri J. Mavriplis ICASE, Hampton, Virginia ICASE NASA Langley Research Center
More informationEngineering of Online Optimizing Control - A Case Study: Reactive SMB Chromatography
Proceedings of the 7th World Congress The International Federation of Automatic Control Engineering of Online Optimizing Control - A Case Study: Reactive SMB Chromatography Achim Küpper, Sebastian Engell
More informationAdsorption dynamics and effects of carbon to zeolite ratio of layered beds for multicomponent gas adsorption
Korean J. Chem. Eng., 28(2),583-590 (2011) DOI: 10.1007/s11814-010-0399-9 INVITED REVIEW PAPER Adsorption dynamics and effects of carbon to zeolite ratio of layered beds for multicomponent gas adsorption
More informationGas hold-up and bubble diameters in a gassed oscillatory baed column
Chemical Engineering Science 56 (2001) 6143 6148 www.elsevier.com/locate/ces Gas hold-up and bubble diameters in a gassed oscillatory baed column M. S. N. Oliveira, X. Ni Department of Mechanical and Chemical
More informationNew Jersey Institute of Technology, University Heights. Newark, NJ A thin cylindrical ceramic sample is placed in a single mode microwave
Microwave Heating and Joining of Ceramic Cylinders: A Mathematical Model Michael R. Booty and Gregory A. Kriegsmann Department of Mathematics, Center for Applied Mathematics and Statistics New Jersey Institute
More informationConjugate lter approach for shock capturing
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2003; 19:99 110 (DOI: 10.1002/cnm.573) Conjugate lter approach for shock capturing Yun Gu 1 1; 2; ; and G. W. Wei 1 Department
More informationEvaluation of diffusional resistances in the process of glucose isomerization to fructose by immobilized glucose isomerase
Enzyme and Microbial Technology 28 (2001) 246 252 www.elsevier.com/locate/enzmictec Evaluation of diffusional resistances in the process of glucose isomerization to fructose by immobilized glucose isomerase
More informationMultistage Schemes With Multigrid for Euler and Navier-Stokes Equations
NASA Technical Paper 363 Multistage Schemes With Multigrid for Euler and Navier-Stokes Equations Components and Analysis R. C. Swanson Langley Research Center Hampton, Virginia Eli Turkel Tel-Aviv University
More informationAdsorption (Ch 12) - mass transfer to an interface
Adsorption (Ch 12) - mass transfer to an interface (Absorption - mass transfer to another phase) Gas or liquid adsorption (molecular) onto solid surface Porous solids provide high surface area per weight
More informationHeuristic design of pressure swing adsorption: a preliminary study
Heuristic design of pressure swing adsorption: a preliminary study S. Jain a, A.S. Moharir a,p. Li b, *, G. Wozny b a Department of Chemical Engineering, Indian Institute of Technology, Bombay 400 076,
More informationt x 0.25
Journal of ELECTRICAL ENGINEERING, VOL. 52, NO. /s, 2, 48{52 COMPARISON OF BROYDEN AND NEWTON METHODS FOR SOLVING NONLINEAR PARABOLIC EQUATIONS Ivan Cimrak If the time discretization of a nonlinear parabolic
More informationNOTICE WARNING CONCERNING COPYRIGHT RESTRICTIONS: The copyright law of the United States (title 17, U.S. Code) governs the making of photocopies or
NOTICE WARNING CONCERNING COPYRIGHT RESTRICTIONS: The copyright law of the United States (title 17, U.S. Code) governs the making of photocopies or other reproductions of copyrighted material. Any copying
More informationPressure Swing Adsorption: A Gas Separation & Purification Process
Pressure Swing Adsorption: A Gas Separation & Purification Process Pressure swing adsorption is an adsorption-based process that has been used for various gas separation and purification purposes. Separation
More informationChemical Reaction Engineering Prof. Jayant Modak Department of Chemical Engineering Indian Institute of Science, Bangalore
Chemical Reaction Engineering Prof. Jayant Modak Department of Chemical Engineering Indian Institute of Science, Bangalore Lecture No. # 26 Problem solving : Heterogeneous reactions Friends, in last few
More informationANALYSIS OF A COMPOSITIONAL MODEL FOR FLUID FLOW IN POROUS MEDIA ZHANGXIN CHEN, GUAN QIN, AND RICHARD E. EWING
ANALYSIS OF A COMPOSITIONAL MODEL FOR FLUID FLOW IN POROUS MEDIA ZHANGXIN CHEN, GUAN QIN, AND RICHARD E. EWING Abstract. In this paper we consider a compositional model for three-phase multicomponent uid
More informationPERCOLATION PARAMETERS OF BIOMOLECULES ADSORPTION ON TWO-DIMENSIONAL NETWORK SURFACES
Proceedings of COBEM 2009 Copyright 2009 by ABCM 20th International Congress of Mechanical Engineering November 15-20, 2009, Gramado, RS, Brazil PERCOLATION PARAMETERS OF BIOMOLECULES ADSORPTION ON TWO-DIMENSIONAL
More informationMathematics Research Report No. MRR 003{96, HIGH RESOLUTION POTENTIAL FLOW METHODS IN OIL EXPLORATION Stephen Roberts 1 and Stephan Matthai 2 3rd Febr
HIGH RESOLUTION POTENTIAL FLOW METHODS IN OIL EXPLORATION Stephen Roberts and Stephan Matthai Mathematics Research Report No. MRR 003{96, Mathematics Research Report No. MRR 003{96, HIGH RESOLUTION POTENTIAL
More informationNovel determination of dierential-equation solutions: universal approximation method
Journal of Computational and Applied Mathematics 146 (2002) 443 457 www.elsevier.com/locate/cam Novel determination of dierential-equation solutions: universal approximation method Thananchai Leephakpreeda
More informationError in the Estimation of Effective Diffusion Coefficients from Sorption Measurements*
Error in the Estimation of Effective Diffusion Coefficients from Sorption Measurements* D. BOBOK and E. BESEDOVÁ Department of Chemical and Biochemical Engineering, Faculty of Chemical Technology, Slovak
More informationShort-cut distillation columns based on orthogonal polynomials of a discrete variable
Short-cut distillation columns based on orthogonal polynomials of a discrete variable Peter Lory Institut für Wirtschaftsinformatik, Universität Regensburg, D 93040 Regensburg, Germany Dedicated to Professor
More informationon a Stochastic Current Waveform Urbana, Illinois Dallas, Texas Abstract
Electromigration Median Time-to-Failure based on a Stochastic Current Waveform by Farid Najm y, Ibrahim Hajj y, and Ping Yang z y Coordinated Science Laboratory z VLSI Design Laboratory University of Illinois
More informationSUPERCRITICAL FLUID CHROMATOGRAPHY AS SUCCESSFUL SEPARATION TOOL IN CHEMICAL AND PHARMACEUTICAL INDUSTRY
9th Meeting on Supercritical Fluids Trieste (Italy), 13-16 June 2004 SUPERCRIAL FLUID CHROMATOGRAPHY AS SUCCESSFUL SEPARATION TOOL IN CHEMICAL AND PHARMACEUAL INDUSTRY M. Johannsen *, S. Peper, G. Brunner
More informationTime (s)
TO: Joseph M. Powers FROM: Christian D. Hanson DATE: 18 April 1997 RE: ME 334 Project (Part II) 1 Problem Description AFortune 500 hundred company, Gooey Foods, Inc., wants to develop a piping system which
More informationA general theory of discrete ltering. for LES in complex geometry. By Oleg V. Vasilyev AND Thomas S. Lund
Center for Turbulence Research Annual Research Briefs 997 67 A general theory of discrete ltering for ES in complex geometry By Oleg V. Vasilyev AND Thomas S. und. Motivation and objectives In large eddy
More informationP 1 P * 1 T P * 1 T 1 T * 1. s 1 P 1
ME 131B Fluid Mechanics Solutions to Week Three Problem Session: Isentropic Flow II (1/26/98) 1. From an energy view point, (a) a nozzle is a device that converts static enthalpy into kinetic energy. (b)
More informationAsymptotic solution of the Boltzmann equation for the shear ow of smooth inelastic disks
Physica A 275 (2000) 483 504 www.elsevier.com/locate/physa Asymptotic solution of the Boltzmann equation for the shear ow of smooth inelastic disks V. Kumaran Department of Chemical Engineering, Indian
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT 18-05 Efficient and robust Schur complement approximations in the augmented Lagrangian preconditioner for high Reynolds number laminar flows X. He and C. Vuik ISSN
More informationA SYSTEMATIC APPROACH TO DISTILLATION COLUMN CONTROL. Sigurd Skogestad. Manfred Morari. California Institute of Technology. Pasadena, CA 91125
A SYSTEMATIC APPROACH TO DISTILLATION COLUMN CONTROL Sigurd Skogestad Manfred Morari California Institute of Technology Chemical Engineering, 206-41 Pasadena, CA 91125 Present aliation: Chem.Eng., Norwegian
More informationDetermination of single-particle drying kinetics in an acoustic levitator
Chemical Engineering Journal 86 (2002) 217 222 Determination of single-particle drying kinetics in an acoustic levitator Carola Groenewold 1, Carsten Möser, Hans Groenewold, Evangelos Tsotsas Institute
More informationAn improved MPS method for numerical simulations of convective heat transfer problems
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2006; 51:31 47 Published online 15 November 2005 in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/d.1106
More informationArticle published by EDP Sciences and available at or
http://www.emath.fr/proc/vol.1/ Axisymmetric Vortex Breakdown in a Pipe S. Wang and Z. Rusak Department of Mechanical Engineering, Aeronautical Engineering and Mechanics Rensselaer Polytechnic Institute,
More informationModeling, Simulation and Optimization of a Cross Flow Molten Carbonate Fuel Cell
Modeling, Simulation and Optimization of a Cross Flow Molten Carbonate Fuel Cell P. Heidebrecht 1, M. Mangold 2, M. Gundermann 1, A. Kienle 2,3, K. Sundmacher 1,2 1 Otto-von-Guericke-University Magdeburg,
More informationMonte Carlo Simulation of Long-Range Self-Diffusion in Model Porous Membranes and Catalysts
Monte Carlo Simulation of Long-Range Self-Diffusion in Model Porous Membranes and Catalysts Brian DeCost and Dr. Sergey Vasenkov College of Engineering, University of Florida Industrial processes involving
More informationA First Course on Kinetics and Reaction Engineering Unit 12. Performing Kinetics Experiments
Unit 12. Performing Kinetics Experiments Overview Generating a valid rate expression for a reaction requires both a reactor and and an accurate mathematical model for that reactor. Unit 11 introduced the
More informationDynamic Process Models
Dr. Simulation & Optimization Team Interdisciplinary Center for Scientific Computing (IWR) University of Heidelberg Short course Nonlinear Parameter Estimation and Optimum Experimental Design June 22 &
More informationNonlinear adaptive control for multivariable chemical processes
Chemical Engineering Science 56 (1) 6781 6791 www.elsevier.com/locate/ces Nonlinear adaptive control for multivariable chemical processes Bingjun Guo, Aiping Jiang, Xiangming Hua, Arthur Jutan Department
More informationOptimization of transport-reaction processes using nonlinear model reduction
Chemical Engineering Science 55 (2000) 4349}4366 Optimization of transport-reaction processes using nonlinear model reduction Eugene Bendersky Panagiotis D. Christo"des* Department of Chemical Engineering
More informationOutline. Definition and mechanism Theory of diffusion Molecular diffusion in gases Molecular diffusion in liquid Mass transfer
Diffusion 051333 Unit operation in gro-industry III Department of Biotechnology, Faculty of gro-industry Kasetsart University Lecturer: Kittipong Rattanaporn 1 Outline Definition and mechanism Theory of
More informationTarget Simulations. Roman Samulyak in collaboration with Y. Prykarpatskyy, T. Lu
Muon Collider/Neutrino Factory Collaboration Meeting May 26 28, CERN, Geneva U.S. Department of Energy Target Simulations Roman Samulyak in collaboration with Y. Prykarpatskyy, T. Lu Center for Data Intensive
More informationON SCALAR CONSERVATION LAWS WITH POINT SOURCE AND STEFAN DIEHL
ON SCALAR CONSERVATION LAWS WITH POINT SOURCE AND DISCONTINUOUS FLUX FUNCTION STEFAN DIEHL Abstract. The conservation law studied is @u(x;t) + @ F u(x; t); x @t @x = s(t)(x), where u is a f(u); x > 0 concentration,
More informationUltraltration Membrane Separation of Coee
Ultraltration Membrane Separation of Coee Ole Håvik Bjørkedal olehb@stud.ntnu.no Therese Bache theresba@stud.ntnu.no October 31, 2013 Abstract This experiment was performed as part of the TKP4105/TKP4110
More informationflow Cross-sectional Area, A(x) Bubble Radius, R(x,t)
1997 ASME Fluids Engineering Division Summer Meeting FEDSM'97 June 2226, 1997 FEDSM973532 ONEDIMENSIONAL BUBBLY CAVITATING FLOWS THROUGH A CONVERGINGDIVERGING NOZZLE YiChun Wang Department of Mechanical
More informationLEARNING OBJECTIVES CHEM 212: SEPARATION SCIENCE CHROMATOGRAPHY UNIT. Thomas Wenzel, Bates College. In-class Problem Set Extraction.
LEARNING OBJECTIVES CHEM 212: SEPARATION SCIENCE CHROMATOGRAPHY UNIT Thomas Wenzel, Bates College In-class Problem Set Extraction Problem #1 1. Devise a scheme to be able to isolate organic acids, bases
More informationSolution of the Two-Dimensional Steady State Heat Conduction using the Finite Volume Method
Ninth International Conference on Computational Fluid Dynamics (ICCFD9), Istanbul, Turkey, July 11-15, 2016 ICCFD9-0113 Solution of the Two-Dimensional Steady State Heat Conduction using the Finite Volume
More informationA Novel Dissolution-Diffusion Model for Investigation of Drug Release from Polymeric Microspheres
European Symposium on Computer Arded Aided Process Engineering 15. Puigjaner and A. Espuña (Editors) 2005 Elsevier Science B.V. All rights reserved. A Novel Dissolution-Diffusion Model for Investigation
More informationSimulation of Plume-Spacecraft Interaction
Simulation of Plume-Spacecraft Interaction MATÍAS WARTELSKI Master of Science Thesis in Space Physics Examiner: Prof. Lars Blomberg Space and Plasma Physics School of Electrical Engineering Royal Institute
More informationA memory-ecient nite element method for systems of reaction diusion equations with non-smooth forcing
Available online at www.sciencedirect.com Journal of Computational and Applied Mathematics 169 (2004) 431 458 www.elsevier.com/locate/cam A memory-ecient nite element method for systems of reaction diusion
More informationUsing Model Order Reduction to Accelerate Optimization of Multi-Stage Linear Dynamical Systems
Using Model Order Reduction to Accelerate Optimization of Multi-Stage Linear Dynamical Systems Yao Yue, Suzhou Li, Lihong Feng, Andreas Seidel-Morgenstern, Peter Benner Max Planck Institute for Dynamics
More informationDiffusion and Adsorption in porous media. Ali Ahmadpour Chemical Eng. Dept. Ferdowsi University of Mashhad
Diffusion and Adsorption in porous media Ali Ahmadpour Chemical Eng. Dept. Ferdowsi University of Mashhad Contents Introduction Devices used to Measure Diffusion in Porous Solids Modes of transport in
More information