Adaptive predictive collocation with a cubic spline interpolation function for convection-dominant fixed-bed processes: Application to a fixed-bed adsorption process

Interference effects and nonlinearities that are inherent in an adsorption equilibrium add challenges to the task of modeling fixed-bed adsorption processes. In particular, the interference effect tends to produce steep transient concentration fronts, which are enhanced by the high convection to dispersion ratio that is found in normal bed operation conditions. This causes spatial stiffness and renders numerical computations difficult and inaccurate. To solve this problem, we propose a novel numerical procedure using gradient-directed adaptive predictive collocation with a cubic spline interpolation function and far-side boundary conditions. The spatial domain is divided into fixed subdomains and the number of collocation points for each subdomain is adaptively adjusted according to the present location and advancement speed of the maximum gradient in each subdomain. The proposed method was applied to a fixed-bed adsorption process capturing CO2 from a CO2/N2 gas mixture and its effectiveness was compared with that of cubic spline and orthogonal collocation methods with fixed interpolation points, respectively. The effects of kinetic rates on the simulation results were also investigated from the viewpoint of accuracy and CPU time.

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

[2]  R. T. Yang,et al.  Kinetic separation by pressure swing adsorption : method of characteristics model , 1990 .

[3]  Woohyun Yun,et al.  The use of cubic spline and far-side boundary condition for the collocation solution of a transient convection-diffusion problem , 2007 .

[4]  K. S. Knaebel,et al.  Pressure swing adsorption: Development of an equilibrium theory for binary gas mixtures with nonlinear isotherms , 1989 .

[5]  Soon-Haeng Cho,et al.  NON-ISOTHERMAL ADSORPTION OF NITROGEN-CARBON DIOXIDE MIXTURE IN A FIXED BED OF ZEOLITE-X , 1994 .

[6]  R. T. Yang,et al.  Bulk separation of multicomponent gas mixtures by pressure swing adsorption: Pore/surface diffusion and equilibrium models , 1985 .

[7]  Lorenz T. Biegler,et al.  Optimization of a Pressure-Swing Adsorption Process Using Zeolite 13X for CO2 Sequestration , 2003 .

[8]  R. T. Yang,et al.  Comparison of Activated Carbon and Zeolite 13X for CO2 Recovery from Flue-Gas by Pressure Swing Adsorption , 1995 .

[9]  Kwang Soon Lee,et al.  Reduced-order distillation model using collocation method with cubic splines , 1991 .

[10]  Alírio E. Rodrigues,et al.  Solution of partial differential equations systems by the moving finite element method , 1992 .

[11]  D. Do,et al.  Adsorption analysis : equilibria and kinetics , 1998 .

[12]  D. Kunii,et al.  Heat transfer characteristics of porous rocks , 1960 .

[13]  S. Ergun.  Fluid flow through packed columns , 1952 .

[14]  Dianne E. Wiley,et al.  Reducing the Cost of CO2 Capture from Flue Gases Using Pressure Swing Adsorption , 2008 .

[16]  R. T. Yang,et al.  Gas Separation by Adsorption Processes , 1987 .

[17]  Q. Yu,et al.  Computer simulations of the dynamics of multicomponent ion exchange and adsorption in fixed beds. Gradient-directed moving finite element method , 1989 .

[18]  Alírio E. Rodrigues,et al.  Applications of a moving finite element method , 2001 .

[19]  Keith Miller,et al.  Design and Application of a Gradient-Weighted Moving Finite Element Code I: in One Dimension , 1998, SIAM J. Sci. Comput..

[20]  L. Biegler,et al.  Simulation and Optimization of a Pressure Swing Adsorption System: Recovering Hydrogen from Methane , 2005 .

[21]  T. N. E. Greville,et al.  Theory and applications of spline functions , 1969 .

[22]  Soon-Haeng Cho,et al.  Numerical Analysis on the Power Consumption of the PSA Process for Recovering CO2 from Flue Gas , 2002 .