Regular Article
Design of adsorption column for reclamation of methyldiethanolamine using homogeneous surface diffusion model
Department of Chemical Engineering, Khalifa University, P.O. Box: 127788, Abu Dhabi, United Arab Emirates
^{*} Corresponding authors: pravin.kannan@ku.ac.ae; fawzi.banat@ku.ac.ae
Received:
7
April
2020
Accepted:
14
September
2020
A predictive simulation model was applied to design a fixedbed adsorber for studying the removal of Total Organic Acid (TOA) anions from lean Methyldiethanolamine (MDEA) solution using Calcium Alginate Bentonite (CAB) clay hybrid composite adsorbent. Unlike other conventional techniques typically used for packed bed design, the predictive Homogeneous Surface Diffusion Model (HSDM) does not require any test column breakthrough curves a priori. Mass transfer coefficients and isotherm model parameters are provided as input data to HSDM for simulating column breakthrough curves. Various isotherm models were fitted to batch equilibrium data for TOA adsorption on CAB composite adsorbent. Based on Akaike Information Criterion (AIC), Freundlich isotherm was selected and the model parameters were obtained by nonlinear regression. Film transfer coefficients and surface diffusivities were determined using appropriate empirical correlations available in the literature. HSDM predictions were first validated using labscale column adsorption data generated at lower residence times. The effects of dimensionless numbers (Biot and Stanton) on breakthrough times were investigated using the dimensionless HSDM system and a suitable scaleup regime (Bi& ~& 1 and St& >& 10) was established wherein the sensitivity of mass transfer parameters would be minimal. Using similitude rules on key design parameters, a pilotscale adsorption column was designed and breakthrough curves were generated using the validated HSDM. The appropriateness of the design technique was verified by comparing the estimated breakthrough data and column design parameters with conventional scaleup and kinetic approaches.
© P. Kannan et al., published by IFP Energies nouvelles, 2020
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Almost all natural gas has H_{2}S, CO_{2} or both that needs to be removed before the gas is pumped through transmission pipelines. The sweetening process is carried out using aqueous Methyldiethanolamine (MDEA, 45–50& wt.%) in a regeneration column where heat is applied to strip the acid gas components and recover lean/aqueous MDEA solution (Keewan et al., 2018; Mehassouel et al., 2018; Younas and Banat, 2014). However, contamination of industrial lean MDEA from heat stable salts such as total organic acids (produced by the reaction between aerial oxygen and CO_{2}/H_{2}S) and heavy metal ions like chromium, lead, etc. (produced from the makeup water or due to the metal corrosion or erosion caused during the continuous running of the plant) always remain a challenge to the gas industry (Cummings et al., 2007; Pal et al., 2015).
The acidic heat stable salts play a vital role in the regeneration column, acting as enhancers for the stripping process. On the other hand, the presence of HSS (Heat Stable Salts) in a lean amine solution is detrimental to the absorption process intended for amine enrichment (Verma and Verma, 2009; Weiland, 2008). Hence, partial removal of Total Organic Acid (TOA) anions from aqueous amines is crucial for avoiding some operational issues encountered during natural gas sweetening process. Different methods have been used for the removal of Heat Stable Salts (HSS) from amine solvents in natural gas sweetening units. Currently, vacuum distillation, electrodialysis, ion exchange, and adsorption are used for the removal of HSS from lean MDEA solutions. Adsorption is widely used among others and has been in practice for several years. The technique is efficient, easy to operate, and requires low maintenance cost (Edathil et al., 2020). The removal has mainly been facilitated through the development of novel adsorbents, including Calcium Alginate Bentonite (CAB) clay composites that serve to remove TOA and metal ions from lean amine solvents (Pal et al., 2013). The oxidative degradation of MDEA produces high concentrations of organic acid anions that remain as contaminants in the solvent. MDEA contains both amine and hydroxyl groups that adhere strongly to these contaminants, thereby making separation by any available technique largely challenging. Calcium alginate biopolymer contains two functional groups, namely carboxyl and hydroxyl groups that serve as adsorption sites and aid in the removal of organic acid anions from a lean MDEA solution. Further, reinforcing bentonite into the alginate matrix increases the mechanical properties, density, and adsorptive efficiency (Edathil et al., 2018; Pal et al., 2019).
However, testing the efficiency of the adsorbent in a broad range of process conditions is sometimes not feasible due to the sensitive nature of MDEA that would otherwise alter the physical properties. For instance, it is a challenging task to generate adsorption equilibrium data as any dilution to lean amine solution would result in significant changes in pH, viscosity and other molecular properties that would significantly affect the stripping tendency of the solvent. Due to these limitations in adsorbent performance testing, reliable simulation capable of accurately modeling the breakthrough behavior under real plant conditions is imperative. Such models could be beneficial for scaleup studies by identifying critical design parameters, interpreting labscale test results, and designing fullscale adsorber.
For the design of fixed bed adsorber, many methods have been developed that could be categorized into shortcut or scoping methods and rigorous methods (Crittenden and Thomas, 1998; Patel, 2019). Design of fixedbed adsorption column by shortcut methods including Mass Transfer Zone Length (MTZL) model, Length of Unused Bed (LUB), Empty Bed Contact Time (EBCT) method, Bed Depth Service Time (BDST) method, Transfer Unit Approach (NTU and HTU) and breakpoint capacity methods requires labscale, pilotscale or plantscale experimentally determined Breakthrough Curve (BTC) data. On the other hand, rigorous methods are based on complex solutions of the conservation, transport and equilibrium relationships. Under certain assumptions, a simplified analytical solution (equilibrium and reaction rate models) could be obtained for BTC prediction that is, however, limited for specific cases.
Predictive models, on the other hand, consider all the subprocesses of film transfer, intraparticle diffusion, and adsorption kinetics to generate the complete sinusoidal shaped BTC. This is achieved by the development of advanced numerical techniques that aid in solving the complex set of partial differential equations used to describe the column behavior. A brief review of literature indicates that predictive models have been developed and employed for ion exchange and fixedbed adsorption operations in wastewater/water treatment (Chowdhury et al., 2015; Crittenden et al., 1986a; Hand et al., 1984, 1989; Hudaya and Rachmat, 2019; Srivastava et al., 2008). The most commonly used general rate models are usually considered as complete models since they predict the entire breakthrough curves of the adsorption process (Xu et al., 2013). Among these, the Homogeneous Surface Diffusion Model (HSDM) has been successfully proved and thus, is widely used for predicting the fixed bed adsorber dynamics of many adsorbate/adsorbent systems. Generally, a complete BTC model would constitute the material balance equation, the equilibrium relationship and a set of equations that describe the external and internal mass transfer processes. Since dispersion is negligible in typical adsorption process flow conditions, mass balance equation based on plug flow type was considered. The superficial velocity was assumed to be constant along the bed as removal of trace TOA components from bulk MDEA solution will pose a negligible effect on material balance. It was also assumed that surface diffusion is the predominant intraparticle mass transfer mechanism and is not a function of local liquid concentration. Many complete BTC models that encompass all three mass transfer processes (film, pore, and surface diffusion) have been reported in the literature (Crittenden et al., 1987). These models differ in the assumptions about the flow processes and mass transfer mechanisms used to describe the column behavior. Homogeneous Surface Diffusion Model (HSDM) has been successfully employed to predict the fixedbed dynamics for many adsorbateadsorbent systems reported earlier (Hand et al., 1984). Many simulation packages [FAST, AdDesignS^{TM}] that have been developed utilize HSDM (dimension and dimensionless forms), along with a powerful numerical solver for simultaneously solving the complex transportreaction equations and nonlinear adsorption isotherms to predict the breakthrough behavior of the column satisfactorily. Like other predictive models, HSDM requires minimal equilibrium data on isotherms and mass transfer characteristics that could be sourced from literature or through labscale experiments.
For the masstransfer controlled process, column design could be solely based on techniques that employ similitudes in mass transfer characteristics. A frequently employed method for fullscale performance prediction is the Rapid SmallScale Column Test (RSSCT) that employs scaling relationships of the design and operational parameters between the smallscale and largescale adsorber. By maintaining perfect similarity and using a relatively smaller adsorbent particle, the RSSCT would exhibit identical breakthrough profiles as the fullscale unit (Crittenden et al., 1991). Few other researchers also accomplished similarity between scales by maintaining the dimensional parameters namely, solute distribution parameter (D_{g}), the surface diffusion modulus (Ed = St/Bi), Stanton number (St) and Reynolds number (Re) (Crittenden et al., 1986b; Hand et al., 1983). Equating the dimensionless parameters of the smallscale and fullscale column, mathematical equations describing the relationships between the critical design parameters were established. However, the derived relationships required varying adsorbent particle diameters between the scales to maintain similitudes. Predictive models are more advantageous when constant particle size needs to be maintained between the scales to minimize any operational and hydraulic problems in fullscale unit. This is particularly true for some adsorbent materials where it is practically impossible to engineer uniform particle size during synthesis. Hence, in the scaleup of fixedbed adsorption columns, typically particle size and hydraulic loading rate or superficial velocity remain unchanged between the scales. However, it is not feasible to keep the velocity fixed during scaleup due to geometrical design limitations in L/D ratio of the labscale column (Inglezakis and Poulopoulos, 2006). In fact, when higher velocities were used at larger scales, it offered increased advantages in terms of driving the system to diffusion control rather than film control. In other words, it resulted in higher breakthrough times as film resistances are minimized, and uniform flow distribution inside the column. Any issues associated with higher velocities, like attrition and fluidization, could be easily resolved by switching to downflow operation in largescale unit.
Inglezakis and Poulopoulos (2006) have summarized the critical design parameters and the influence it may have on the performance of fixedbed operations. As it can be seen from Table 1, apart from particle size, contact time is another critical parameter that has a significant effect on breakthrough times. By maintaining similar particle size and contact times, results from labscale studies could be directly transferred to plantscale making the scaleup procedure more simple and precise. However, contact times and the resulting breakthrough times encountered in labscale are too small and impractical to be employed in plantscale operations where higher breakthrough times are desired. Therefore, apart from other similitude rules mentioned in Table 1, particle size and superficial velocity are typically maintained constant in the conventional scaleup techniques.
Critical design parameters reproduced from Inglezakis and Poulopoulos (2006).
As it can be seen, the optimum scale and design depends on several key variables associated with the adsorbent characteristics and process variables. It can be challenging to determine the appropriate pilotscale residence time and loading rates with an experimental method that is based only on column breakthrough studies. Numerous experiments may be required to conclude on the design values, including flow characteristics making the approach costly and timeconsuming. Hence, an efficient predictive model in conjunction with dimensionless number(s) that define the mass transfer behavior would be necessary to ensure similitudes during the design process. For example, column adsorption process has been successfully modeled by HSDM and characterized using Bi number by several researchers (Hand et al., 1983; Lee et al., 1983; Smith, 1997; Traegner and Suidan, 1989; Wolborska, 1999). However, in an application towards adsorption of arsenate and other contaminants using granular ferric hydroxide filters, Sperlich et al. (2008) concluded that characterization based on Bi number alone was not sufficient to completely characterize the HSDM model. It was suggested to complement Bi with St number in order to predict the BTC satisfactorily and no dimensionless similarities were required. In the current work, a similar approach has been demonstrated to design a pilotscale fixedbed adsorption column for TOA adsorption using CAB particles.
Simulation and column design of TOA adsorption from lean MDEA solutions has never been reported in literature till date. The objective of the present work is to utilize HSDM for designing a fixed bed adsorber removing TOA present in lean MDEA solution. A series of labscale breakthrough experiments would be conducted at different hydraulic loading rates using different sized columns packed with CAB adsorbent. These experiments would serve to validate the HSDM results and assess the qualitative features of the model in terms of breakthrough curves. The optimum design range would be established based on sensitivity analysis of Bi and St in conjunction with certain critical design parameter guidelines available in literature. Finally, the appropriateness of the design technique would be verified by comparing the scaleindependent HSDM predictions with conventional scaleup and kinetic approaches. Figure 1 illustrates the overall design methodology employed in this work.
Fig. 1 HSDM based fixedbed adsorber design technique employed in this study. 
2 Method description
2.1 Selection of isotherm model
Equilibrium experiments were performed using a batch system to generate TOACAB adsorbent isotherm data, and the results have been published in the literature (Edathil et al., 2020). In this study, the batch experimental data were evaluated with the isotherm models of Langmuir, Freundlich, Jovanovic, Twostep Langmuir, LangmuirFreundlich, FritzSchlunder and Li, and the constants were obtained by nonlinear regression. A standard procedure in choosing an appropriate model equation is based on the value of regression coefficient (R^{2}). However, regression coefficients and other error functions are prone to error when comparing equations with different degrees of freedom (DávilaJiménez et al., 2014). To overcome this problem, Akaike Information Criterion (AIC) was used to rank the isotherm models since it is more sensitive to model deviations and takes into consideration the number of parameters in an equilibrium isotherm model, unlike other error functions. The AIC values for all the models mentioned above were calculated using the following expression:
(1)where N is the number of isotherm data points, SSE is the sum of squared residuals and N_{p} is the number of fitted parameters. As can be seen, AIC takes both accuracy and model complexity into account while regression coefficients consider only accuracy as the determining factor. For a given experimental dependent variable response (equilibrium solid concentration) on the independent variable (influent concentration), the abovementioned models with a different number of parameters could be rated based on the values of model AIC’s. The particular isotherm model that exhibits the lowest AIC value would best describe the equilibrium between the adsorbent and adsorbate under consideration.
2.2 Labscale column experiments
Adsorbent material properties of CAB adsorbent including porosity and density were reported in previous works (Edathil et al., 2018). Figure 2 presents a schematic of the fixedbed adsorption setup used in this work for the removal of TOA from lean MDEA solution. The adsorption system was designed to treat industrial lean MDEA solution and provide clean MDEA (without TOA) at the outlet of the column. The system was designed in such a way that the same column could be utilized for both adsorption and regeneration. Adsorption studies were conducted using different borosilicate glass columns (BUCHI, Switzerland) of varying dimensions (d (cm) × h (cm): 1.5& ×& 10, 2.6& ×& 10, 2.6& ×& 23 and 4.6& ×& 23). The column was filled with a known quantity of 2% CAB adsorbent and then lean MDEA solution of known TOA concentration was pumped through the column using a peristaltic pump at the desired flow rate in an upflow mode. Treated MDEA effluent samples were collected from the outlet of the column at definite time intervals and the concentration of TOA ions in the effluent was measured using a UVvis spectrophotometer. The operation of the column was stopped once the concentration of TOA ions in effluent samples reached the influent concentration. All adsorption experiments were performed at room temperature and an influent pH& ~& 10.5. For the column dimensions and flow rates considered in this study, the Empty Bed Contact Time (EBCT) varied from 23.5& min to 382& min.
Fig. 2 Schematic diagram representing the fixedbed column adsorption setup. 
2.3 BTC predictive model
In this work, a complete Breakthrough Curve (BTC) model considering both adsorption equilibrium and kinetics was used to model the real Sshaped BTCs for TOA adsorption on CAB adsorbents. HSDM is simultaneously represented by two partial differential equations that describe the macroscale liquid phase fluid movement over a bed of adsorbent particles (Eq. (2)) and the unsteady state surface diffusion into the spherical adsorbent particle (Eq. (3)):
It would be shown later in this study from the AIC criterion analysis, Freundlich isotherm would be the most suitable model to describe the isothermal equilibrium between the diluted fluid mass and adsorbed mass in solid phase. The nonlinear temperatureindependent Freundlich isotherm equation couples the two partial differential equations through the adsorption term of the transport equation. With the nonlinear adsorption isotherm embedded, it becomes highly challenging to find an analytical solution. Solutions to the set of PDEs along with the boundary and initial conditions were obtained using finite differences method as published in literature (Sperlich et al., 2008). All calculations were performed using the software FAST 2.0 Beta (Fixedbed Adsorption Simulation Tool) that was developed originally for water treatment applications.
2.4 Calculation of mass transfer coefficients
In order to estimate the dimensionless numbers, several model parameters need to be determined either through experiments or through empirical correlations. The column and adsorbent geometric parameters such as particle size, particle density (ρ_{a}), particle porosity (ε_{p}), superficial velocity, empty bed contact time, bed porosity (ε), influent initial concentration (C_{0}) were directly measured. Film diffusion coefficients,k_{L} and surface diffusion coefficients, D_{s} were estimated using the Gnielinski correlation which is a function of Reynolds and Schmidt numbers and Sontheimer correlation, respectively (Sperlich et al., 2008):
In the above equations, D represents molecular diffusivity and is calculated as:
(6)where, η is kinematic viscosity, V_{b} is normal molar volume, τ_{p} is tortuosity (set to 1), q_{0} is equilibrium loading, and SPDFR is the surface to pore diffusion flux ratio.
2.5 HSDM model validation
HSDM can be satisfactorily used to predict breakthrough curves for different adsorption column sizes and flow conditions, without any experimental data a priori. However, in order to establish reliability and confidence over the chosen isotherm model and mass transfer correlations, HSDM predictions were first compared to labscale experimental data generated using different sized adsorption columns and flow rates with CAB adsorbent of 1& mm particle size. Five different adsorption column experiments were designed representing varying residence times. The bed porosity was maintained between 0.35 and 0.5 by adjusting the mass of adsorbent loaded into the column. Experiments were run using lean MDEA with an initial concentration of 3250–3500& ppm and breakthrough data were recorded as a function of residence time. The corresponding Biot and Stanton numbers have also been shown in Table 2 along with other process conditions. The same feed and process conditions were provided as input to the HSDM model and simulations were performed to predict the complete breakthrough curves. Instead of comparing C/C_{o} at certain points, HSDM simulation provides a way to compare the entire BTC which is significant in establishing similarities in flow pattern between scales.
Column parameters used in HSDM validation study.
2.6 Sensitivity analysis
To investigate the effect of these dimensionless numbers on breakthrough curves, PérezFoguet et al. developed a dimensionless analysis of HSDM (PérezFoguet et al., 2013). The macroscale liquid phase mass transport equation and the intraparticle diffusion equation were transformed into a dimensionless form using dimensionless variables for the liquid and solid phase concentration, axial position and contact time. The derivation of the dimensionless model equations and the corresponding boundary and initial conditions could be found elsewhere (Sperlich et al., 2008). The two characteristic partial differential equations of the HSDM were transformed into a system of two ordinary differential equations coupled with the macroscale transportreaction PDE. Discretization was accomplished using a discontinuous Galerkin scheme and the overall system evolution was integrated with a timemarching scheme based on the forward Euler method. The overall system was successfully used to simulate the adsorption of different adsorbates on granular ferric hydroxide.
The detailed dimensionless analysis of HSDM presented by PérezFoguet et al. (2013) facilitates establishing limit behavior of the model to determine the values of Bi and St number (PérezFoguet et al., 2013). The influence of dimensionless numbers, including Bi, St, Ed, D_{g} and n on the breakthrough curves was thoroughly assessed using the dimensionless HSDM. Excerpts from this work have been presented in Section 3.3 and the results have been directly applied to determine the working range of Bi and St numbers.
Design parameters from this technique were compared with direct scaleup and kinetic approaches. Details of these two conventional techniques, including key equations have been summarized in the Supplementary Section.
3. Results and discussion
3.1 Equilibrium studies
Different isotherm models were fitted to batch equilibrium data and the model parameters were determined by nonlinear regression. Equilibrium isotherms for TOA adsorption on CAB adsorbents have been presented in the Supplementary Section, see Figure S1. Figure 3 shows the comparison of simulated data from various models with batch equilibrium data obtained by varying adsorbent mass at 23& °C and Table 3 shows the regressed model parameters and the corresponding AIC values for different isotherm models. It can be noticed that the Freundlich model exhibited the least AIC value indicating the best fit with experimental data. Hence, Freundlich model was incorporated into the HSDM model to describe the equilibrium behavior between the solute present in the liquid and solid phases.
Fig. 3 Comparison of batch equilibrium data with various isotherm models (legend details: Exp – experimental, JC – Jovononic, FR – Freundlich, LF – LangmuirFreundlich, 2sL – 2 stepLangmuir, FS – FritzSchlunder, Lang – Langmuir). 
Summary of isotherm model equations and parameters.
3.2 HSDM model validation
Figure 4 shows a comparison plot of normalized concentration against adsorption time between labscale experimental data and HSDM predictions for different residence times. An inset has been provided in order to better visualize the dynamics at early adsorption times. As shown in Figure 4, HSDM is able to capture the overall process dynamics described by the “S” shaped curve for all the HRT’s considered. Also, the general trend in BTC timeshift with varying HRTs was predicted reasonably well by the HSDM. On comparison of simulated BTCs with experimental data in Figure 4, some discrepancies are evident, especially for the case of HRT& =& 137& min. Many different factors influence the shape of the BTCs, primarily the values of mass transfer coefficients employed in the HSDM. The sensitivity of these coefficients has been analyzed via two dimensionless parameters, namely Bi and St numbers. As outlined in Section 3.3 and in the range of Bi values considered, St remains a significant parameter in determining BTC dynamics. However, in order to minimize the discrepancies between experimental data and simulation, the coefficients should be determined accurately from controlled batch experiments rather than from empirical correlations. Other known possible factors that influence BTC dynamics include HSDM assumptions concerning negligible pore diffusion and unaccounted experimental factors like backpressure and influent pumping issues encountered during column experiments. Nevertheless, as shown in this study, HSDM could be employed for the rapid design of the adsorption column where the application of other techniques is severely challenged by the inherent dependency on experimental data.
Fig. 4 Validation of HSDM with labscale BTC data (lines represent simulation and marker represents lab data, inset depicts the goodness of fit at lower times). 
3.3 Sensitivity analysis
Sensitivity analysis on the dimensionless numbers provides knowledge on the variation of column breakthrough behavior in different scales. Breakthrough curves were predicted for ranges of St and Bi numbers to understand the mass transfer behavior at different process conditions. The analysis is based on a comparison of normalized concentration C/C_{0}versus dimensionless time T, defined as the ratio of operation time to ideal (stoichiometric) breakthrough time.
The values of dimensionless parameters of the labscale experiment used for HSDM validation study are presented in Table 2. Since solute distribution parameter (D_{g}) and Freundlich exponent (n) were determined from batch experiments in this study, they were treated as constants and not included in the sensitivity analysis. Further, since Ed& =& St/Bi, the influence of Ed is dependent on the individual values of Bi and St numbers. Detailed results on the dimensionless analysis could be found elsewhere (PérezFoguet et al., 2013) and only key findings have been discussed in this section. It was shown by PérezFoguet et al. that BTC shapes and limit behaviors for both linear and nonlinear cases (n& ≤& 1) are similar with minor differences in the sharpness of the wavefront. This was also verified in this study by performing a dimensionless analysis on “n” varying from 0.25 to 1 (equivalent to an initial concentration of 7.5–3250& mg/L) while keeping Bi and St constant at 0.9 and 10, respectively. As it could be seen from Figure 5, the effect of n is not so significant on the BTC dynamics and thus extrapolation of Freundlich isotherm towards zero adsorbate concentration may be assumed to impart negligible discrepancies in the HSDM predictions. Hence for all subsequent sensitivity and scaleup studies, the value of n was held constant at the regressed value of 0.8.
Fig. 5 Effect of Freundlich exponent on breakthrough curve dynamics 
A representative plot illustrating the effect of St at low and high Bi numbers, see Figures 6a and 6b, was reproduced and assessed using the dimensionless HSDM system. Variation of St at fixed Bi would be attained by varying the flow velocity of the solvent flowing through the packed bed column. Breakthrough curves shown in both the plots were generated by varying Stanton numbers between 0 and 10^{4} for fixed Bi values of 0.9 and 100 while other dimensionless parameters, D_{g} and n were held constant. At low Bi number, film diffusion dominates and thus acts as the controlling mechanism. Since differences between the results obtained at smaller pairs of Bi number were negligible (PérezFoguet et al., 2013), only Bi& =& 0.9 and Bi& =& 100 were considered for the analysis. Also, since limit behavior was found at these values for Bi, the predicted breakthrough curves are strongly dependent on St number as evidenced in Figure 6. At Bi number close to 1 and St& >& 10, refer Figure 6a, Mass Transfer Zone (MTZ) pattern is fully established as indicated by the sinusoidal BTC shaped curves. This regime (St& >& 10) could be thought of instantaneous adsorption regime where film transfer dominates advection. However, at St& <& 10, the MTZ is still developing and shows a varying trend, indicating relatively slower adsorption and at St& =& 0, the condition changes to no adsorption. Also, since Bi number is small and close to 1, the BTCs are purely dependent on St number within the range from 0 to 10. However, at higher Bi number and thus higher Ed, see Figure 6b, when intraparticle diffusion is rate controlling, BTC shows a varying trend till the value of St is around 1000. The sensitiveness of breakthrough curves at low Stanton number (St& <& 10) was evident in the labscale column experiments reported in previous sections. As a comparison, Test 2 and Test 4 were run at varying Stanton numbers (Test 2: St& =& 2.4, Test 4: St& =& 5.6) and at constant Biot number of 0.76. It can be noticed from Figure 4 that the dynamics of the two breakthrough curves are different with significant difference in breakthrough times.
Fig. 6 Effect of Stanton number on the breakthrough curves at (a) Bi& =& 0.9, (b) Bi& =& 100. 
Thus, columns of different sizes would exhibit similar mass transfer behavior or similar controlling mechanism if Bi was chosen close to 1 and St& ≥& 10. This analysis provided a working range of influent flow rates and loading rates or superficial velocity that could serve as a guideline while designing adsorption columns of larger scales.
HSDM assumes that solid phase mass transfer occurs only by surface diffusion and hence tortuosity and SPDFR are not considered significant. The sensitivity of intraparticle diffusion coefficient (D_{s}) was tested on the BTC characteristics in a wide range from 10^{−11} to 10^{−9}& m^{2}/s, see Figure 7. The effect of D_{s} on the BTC dynamics could be better understood through analysis of the dimensionless Biot number. Since Bi and D_{s} are inversely proportional, increasing D_{s} by two orders of magnitude from 10^{−11} to 10^{−9} would proportionally decrease Bi that would eventually lead to faster adsorption characterized by sharp wavefront. Figure 8 illustrates the dimensionless BTC for different Biot numbers and a fixed St of 10. This implies a case of constant film transfer rate and flow velocity but varying surface diffusion efficient. As it can be noticed form sharpness of the wavefront, adsorption rate increases with decreasing Bi or increasing D_{s}. However, the effect of Bi on BTC shape beyond Bi = 0.1 (or 10^{−1}) is not significant and practically remains unchanged, thereby setting the limits for surface diffusion coefficient.
Fig. 7 Effect of surface diffusion coefficient on BTC dynamics. 
Fig. 8 Effect of Biot number on Dimensionless BTC dynamics. 
By performing a sensitivity analysis on the Freundlich exponent, the impact of extrapolation of model parameters could be analyzed. BTC shapes and limit behaviors for both linear and nonlinear cases (n& ≤& 1) are similar with minor differences in the sharpness of the wavefront (PérezFoguet et al., 2013). This was also verified in this study by performing a dimensionless analysis on “n” varying from 0.25 to 1 (equivalent to an initial concentration of 7.5 to 3250& mg/L) while keeping Bi and St constant at 0.8 and 10, respectively. As it could be seen from Figure 9, the effect of Freundlich exponent is not so significant on the BTC dynamics and thus extrapolation of Freundlich isotherm towards zero adsorbate concentration can be assumed to have negligible influence on the HSDM predictions. Hence for all sensitivity and scaleup studies, the value of n was held constant at the regressed value of 0.8.
Fig. 9 Effect of Freundlich exponent on breakthrough curve dynamics. 
3.4 Design scale assessment
From sensitivity analysis, the range of Bi and St numbers that would exhibit similar mass transfer phenomena for TOA adsorption were determined. Based on these results, a pilotscale adsorption column was designed using the critical design parameter rules recommended by Inglezakis and Poulopoulos, as shown in Table 1. For all subsequent analysis, labscale test column Case 3, shown in Tables 2 and 4, was chosen as the representative labscale design and compared with the pilotscale design parameters.
Summary of labscale and pilotscale adsorption column design parameters.
Figure 10 illustrates the complete breakthrough curves predicted for different scales using the validated HSDM. The scales from the lab to pilot differ in residence time; however, the critical design parameters remain the same. For a labscale column with a residence time of 34& min, the estimated breakthrough time was ca. 26 min that is in good agreement with the experimentally determined value of 28.5& min. Similarly, pilotscale unit with a residence time of 205& min resulted in relatively higher breakthrough times of ca. 330& min.
Fig. 10 Comparison of BTC curves between different scales as predicted by HSDM. 
As expected, breakthrough time to attain a normalized effluent concentration of 0.1 increases with increasing residence time. It should be noted that mass transfer coefficients were calculated by empirical correlations from the literature because they are not scalable from batchscale studies to pilotscale studies due to differences in flow pattern in the reactor. The accuracy of HSDM predictions depends strongly on the appropriateness of these correlations and the estimated equilibrium parameters from batch studies.
Additionally, the column design was carried out using the conventional packedbed scaleup procedure and simplified kinetic models (detailed in the Supplementary Section) and compared with the HSDM design parameters. It is worth noting that both the scaleup and the kinetic approach depends entirely on the breakthrough data generated using test column, either laboratory or pilot plant. In the scaleup approach using the Length of Unused Bed (LUB) model, the loading rate and the unused bed length for both the labscale and pilotscale units had to be maintained constant to obtain similar mass transfer characteristics. In this study, the labscale test column of 4.3& cm diameter and 23& cm height was used with a filtration or loading rate of 0.256& cm^{3}/min& cm^{2} and an EBCT of 96& min. The length of unused bed corresponding to the fraction of bed unused was calculated as 13.1& cm. As per the LUB model, this length should remain the same in scaleup, and thus the length of unused bed in pilotplant design was taken as 13.1& cm and the corresponding breakthrough time was estimated as 345& min.
The kinetic approach was based on the simplified Bohart and Adams model that employed a kinetic rate equation to determine the reaction constant and maximum solid phase loading. However, this approach necessitates a breakthrough volume or time to be specified in the design equation. A breakthrough time of 335& min (taken from HSDM) was used to calculate other design parameters including the mass of adsorbent required and breakthrough volume. Table 5 shows a comparison of various design parameters calculated from different techniques. It can be noted that all three techniques yield similar results for the design of a pilotscale adsorber, however with different complexities. HSDM could be used to design adsorption column of various scales with minimal input about equilibrium data and mass transfer coefficients. However, the other two techniques, scaleup, and kinetic approaches are dependent on the accuracy of the supplied breakthrough data.
Comparison of pilotscale design parameters between various techniques.
4 Conclusion
The objective of the present work was to utilize the predictive homogeneous surface diffusion model for designing a fixedbed adsorber to remove TOA present in lean MDEA solution using CAB adsorbent. As part of the design technique, various isotherm models were fitted to batch equilibrium data, and based on Akaike Information Criterion (AIC), it was found that the equilibrium isotherms were best described by the Freundlich equation. Subsequently, Freundlich isotherm parameters and appropriate mass transfer correlations for film and surface diffusivities were incorporated into the HSDM equation matrix. Numerical solution to the simultaneous transportreaction equations and nonlinear Freundlich equation was executed through commercial software. Simulation results were compared with labscale experimental data collected at lower residence times (<140 min), and it was evident that HSDM could predict breakthrough curves with reasonable accuracy. Dimensionless HSDM equations were employed to describe the limit behavior of the model based on dimensionless numbers, Bi and St. Sensitivity analysis on the two parameters established the operating range for the design units as Bi& ~& 1 and St& >& 10. Further, in order to preserve the flow pattern during scaleup process, key column parameters and similitude rules from literature were reviewed and integrated into the column design. Based on the design guidelines, 30& cm by 150& cm fixedbed adsorber with a continuous throughput of 11.1& L/h was considered suitable. HSDM predicted a 330 min column operating time with an equivalent lean amine treatment capacity of 60 L based on a 10% breakthrough limit for the designed column. Accuracy of the HSDM based design technique was evaluated by comparing with conventional scaleup and kinetic approaches and was found to be in good agreement. Results demonstrated the rapid, easeofuse and accuracy of the HSDM technique for the design of fixedbed adsorption columns for complex systems.
Acknowledgments
The authors would like to acknowledge the support provided by the Gas Research Center (GRC) at Khalifa University under research grant GRC11006.
References
 Chowdhury Z.Z., Hamid S.B., Zain S.M. (2015) Evaluating design parameters for breakthrough curve analysis and kinetics of fixed bed columns for Cu(II) cations using lignocellulosic wastes, BioResources 10, 1, 732–749. [Google Scholar]
 Crittenden B., Thomas W.J. (1998) Adsorption technology and design, ButterworthHeinemann, Woburn, MA, USA. [Google Scholar]
 Crittenden J.C., Berrigan J.K., Hand D.W. (1986a) Design of rapid smallscale adsorption tests for a constant diffusivity, J. Water Pollut. Control Fed. 58, 4, 312–319. [Google Scholar]
 Crittenden J.C., Hand D.W., Arora H., Lykins B.W. (1987) Design considerations for GAC treatment of organic chemicals, J. Am. Water Works Ass. 79, 1, 74–82. [CrossRef] [Google Scholar]
 Crittenden J.C., Hutzler N.J., Geyer D.G., Oravitz J.L., Friedman G. (1986b) Transport of organic compounds with saturated groundwater flow: Model development and parameter sensitivity, Water Resour. Res. 22, 3, 271–284. [CrossRef] [Google Scholar]
 Crittenden J.C., Reddy P.S., Arora H., Trynoski J. (1991) Predicting GAC performance with Rapid SmallScale Column Tests, J. Am. Water Works Ass. 83, 1, 77–87. [CrossRef] [Google Scholar]
 Cummings A.L., Smith G.D., Nelson D.K. (2007) Advances in amine reclaiming: Why there is no excuse for operating a dirty amine system, in: Laurance Reid Gas Conditioning Conference, Dickinson TX, USA, pp. 227–244. [Google Scholar]
 DávilaJiménez M.M., ElizaldeGonzález M.P., GarcíaDíaz E., GonzálezPerea M., GuevaraVilla M.R.G. (2014) Using akaike information criterion to select the optimal isotherm equation for adsorption from solution, Adsorpt. Sci. Technol. 32, 7, 605–622. [CrossRef] [Google Scholar]
 Edathil A.A., Pal P., Banat F. (2018) Alginate clay hybrid composite adsorbents for the reclamation of industrial lean methyldiethanolamine solutions, Appl. Clay Sci. 156, 213–223. [CrossRef] [Google Scholar]
 Edathil A.A., Pal P., Kannan P., Banat F. (2020) Total organic acid adsorption using alginate/clay hybrid composite for industrial lean amine reclamation using fixedbed: Parametric study coupled with foaming, Int. J. Greenh. Gas. Con. 94, 102907. [CrossRef] [Google Scholar]
 Hand D.W., Crittenden J.C., Arora H., Miller J.M., Lykins B.W. (1989) Designing fixedbed adsorbers to remove mixtures of organics, J. Am. Water Works Ass. 81, 1, 67–77. [CrossRef] [Google Scholar]
 Hand D.W., Crittenden J.C., Thacker W.E. (1983) Useroriented batch reactor solutions to the homogeneous surface diffusion model, J. Environ. Eng. 109, 1, 82–101. [CrossRef] [Google Scholar]
 Hand D.W., Crittenden J.C., Thacker W.E. (1984) Simplified models for design of fixedbed adsorption systems, J. Environ. Eng. 110, 2, 440–456. [CrossRef] [Google Scholar]
 Hudaya T., Rachmat V. (2019) Activated carbon fixedbed adsorber design for treating chromium hexavalent wastewater, Makara J. Technol. 22, 3, 135–141. [CrossRef] [Google Scholar]
 Inglezakis V.J., Poulopoulos S.G. (2006) Adsorption, Ion Exchange and Catalysis, Elsevier, Amsterdam, The Netherlands. [Google Scholar]
 Keewan M., Banat F., Pal P., Zain J., Alhseinat E. (2018) Foaming of industrial lean methyldiethanolamine solution in the presence of hydrocarbon and fatty acid based corrosion inhibitors, Oil Gas Sci. Technol.  Rev. IFP Energies nouvelles 73, 76, 1–7. [CrossRef] [Google Scholar]
 Lee M.C., Crittenden J.C., Snoeyink V.L., Ari M. (1983) Design of carbon beds to remove humic substances, J. Environ. Eng. 109, 3, 631–645. [CrossRef] [Google Scholar]
 Mehassouel A., Derriche R., Bouallou C. (2018) Kinetics study and simulation of CO_{2} absorption into mixed aqueous solutions of methyldiethanolamine and hexylamine, Oil Gas Sci. Technol.  Rev. IFP Energies nouvelles 73, 19, 1–10. [CrossRef] [Google Scholar]
 Pal P., AbuKashabeh A., AlAsheh S., Banat F. (2015) Role of aqueous methyldiethanolamine (MDEA) as solvent in natural gas sweetening unit and process contaminants with probable reaction pathway, J. Nat. Gas Sci. Eng. 24, 124–131. [CrossRef] [Google Scholar]
 Pal P., Banat F., AlShoaibi A. (2013) Adsorptive removal of heat stable salt anions from industrial lean amine solvent using anion exchange resins from gas sweetening unit, J. Nat. Gas Sci. Eng. 15, 14–21. [CrossRef] [Google Scholar]
 Pal P., Edathil A.A., Banat F. (2019) Calcium alginate gel and hard beads for the removal of total organic acid anions and heavy metal ions from industrial lean methyldiethanolamine solvent, Polym. Bull. 76, 1, 103–118. [CrossRef] [Google Scholar]
 Patel H. (2019) Fixedbed column adsorption study: a comprehensive review, Appl. Water Sci. 9, 45, 1–17. [CrossRef] [Google Scholar]
 PérezFoguet A., Casoni E., Huerta A. (2013) Dimensionless analysis of HSDM and application to simulation of breakthrough curves of highly adsorbent porous media, J. Environ. Eng. 139, 5, 667–676. [CrossRef] [Google Scholar]
 Smith E.H. (1997) Wave front analysis for design of fixedbed adsorbers, Chem. Eng. Commun. 159, 1, 17–37. [CrossRef] [Google Scholar]
 Sperlich A., Schimmelpfennig S., Baumgarten B., Genz A., Amy G., Worch E., Jekel M. (2008) Predicting anion breakthrough in Granular Ferric Hydroxide (GFH) adsorption filters, Water Res. 42, 8–9, 2073–2082. [CrossRef] [PubMed] [Google Scholar]
 Srivastava V.C., Prasad B., Mishra I.M., Mall I.D., Swamy M.M. (2008) Prediction of breakthrough curves for sorptive removal of phenol by bagasse fly ash packed bed, Ind. Eng. Chem. Res. 47, 5, 1603–1613. [CrossRef] [Google Scholar]
 Traegner U.K., Suidan M.T. (1989) Parameter evaluation for carbon adsorption, J. Environ. Eng. 115, 1, 109–128. [CrossRef] [Google Scholar]
 Verma N., Verma A. (2009) Amine system problems arising from heat stable salts and solutions to improve system performance, Fuel Process. Technol. 90, 4, 483–489. [CrossRef] [Google Scholar]
 Weiland R. (2008) Heat stable salts and amine unit performance, Hydrocarbon World 3, 1, 30–33. [Google Scholar]
 Wolborska A. (1999) External film control of the fixed bed adsorption, Chem. Eng. J. 73, 2, 85–92. [CrossRef] [Google Scholar]
 Xu Z., Cai J., Pan B. (2013) Mathematically modeling fixedbed adsorption in aqueous systems, J. Zhejiang Univ. Sci. A 14, 3, 155–176. [CrossRef] [Google Scholar]
 Younas O., Banat F. (2014) Parametric sensitivity analysis on a GASCO’s acid gas removal plant using ProMax simulator, J. Nat. Gas Sci. Eng. 18, 247–253. [CrossRef] [Google Scholar]
Supplementary Material
Details of these two conventional techniques (scaleup approach, kinetic approach), including key equations have been summarized in the supplementary section.
Fig. S1. Plot of ln ((C_{0}/C_{t})& −& 1) vs. time.
Table S1. Kinetic parameters obtained by linear regression
Fig. S2. Adsorption equilibrium curve, q_{e}vs. C_{e} for the adsorptive removal of TOA using CAB composites.
Supplementary References.
(Access here)All Tables
All Figures
Fig. 1 HSDM based fixedbed adsorber design technique employed in this study. 

In the text 
Fig. 2 Schematic diagram representing the fixedbed column adsorption setup. 

In the text 
Fig. 3 Comparison of batch equilibrium data with various isotherm models (legend details: Exp – experimental, JC – Jovononic, FR – Freundlich, LF – LangmuirFreundlich, 2sL – 2 stepLangmuir, FS – FritzSchlunder, Lang – Langmuir). 

In the text 
Fig. 4 Validation of HSDM with labscale BTC data (lines represent simulation and marker represents lab data, inset depicts the goodness of fit at lower times). 

In the text 
Fig. 5 Effect of Freundlich exponent on breakthrough curve dynamics 

In the text 
Fig. 6 Effect of Stanton number on the breakthrough curves at (a) Bi& =& 0.9, (b) Bi& =& 100. 

In the text 
Fig. 7 Effect of surface diffusion coefficient on BTC dynamics. 

In the text 
Fig. 8 Effect of Biot number on Dimensionless BTC dynamics. 

In the text 
Fig. 9 Effect of Freundlich exponent on breakthrough curve dynamics. 

In the text 
Fig. 10 Comparison of BTC curves between different scales as predicted by HSDM. 

In the text 