CO 2 Adsorption : Experimental Investigation and CFD Reactor Model Validation

The National Energy Technology Laboratory is investigating a new process for CO 2 capture from large sources such as utility power generation facilities as an alternative to liquid amine based adsorption processes. Many of these advanced dry processes are based upon sorbents composed of supported polyamines. In this analysis, experiments have been conducted in a laboratory-scale fluidized bed reactor and compared to CFD reactor predictions using kinetics obtained from TGA tests. Batch experiments were conducted by flowing a mixture of CO 2 , H 2 O, and N 2 (simulated flue gas) through a fluidized bed of sorbent material. The exit gas composition time series data is compared to CFD simulations using a 3-dimensional nonisothermal reacting multiphase flow model. The effects of the gas flow rate, distributor design, and particle size are explored through the CFD simulations. It is shown that the time duration for CO 2 adsorption decreased for an increase in the gas flow. Fluid bed hydrodynamics indicated that there were regions in the reactor where the inert FCC particles segregated and defluidized; without adversely affecting the capacity of the sorbent to adsorb CO 2 . The details of the experimental facility and the model as well as the comparative analysis between the data and the simulation results are discussed.


Introduction
Over the past two to three decades, there has been increasing concern over the importance of carbon dioxide emissions to the environment and the possible effect these emissions have on global climate shifting.Historical atmospheric CO 2 levels have been greater and less than the preindustrialized values shown in Figure 1 [1].Also shown in Figure 1 are the CO 2 concentration levels over the past 30 years [2].These values start in 1980 with a concentration of 338 ppm, already outside the normal range over the past 400,000 years, and increase to a value of 392 ppm in 2012, the last year that data exist.Based on these CO 2 numbers there is little doubt that CO 2 is higher now than in the recent past.Extrapolating the trend in the data taken over the past 30 years backwards (Figure 1) this trend intersects with that of the historical data in the time frame of the beginning of the industrial age.This implies with strong confidence that the increased values of CO 2 are likely caused from human activity.The recent divergence between the temperature data and the CO 2 concentration data refutes assertions that global warming can be attributed to increased CO 2 concentration as does the trend in the temperature variation between 1940 and the late 1970s [3].Even so, the buildup of CO 2 in the atmosphere, whether leading to climate change or not, is attributed to human activity (combustion of fossil fuels associated with industrial and power production) and as such it is prudent of us to investigate mitigation approaches such as the CO 2 capture technology discussed below.
The U.S. Department of Energy has set goals for carbon capture systems at 90% carbon dioxide (CO 2 ) capture with less than a 35% increase in cost of electricity [4].In this approach, the fossil energy system is designed or retrofitted such that a relative pure stream of CO 2 can be separated from the exhaust.The exhaust can then be utilized or sequestered.There are four main categories of these technologies: precombustion, postcombustion, oxy-fuel combustion, and chemical looping [5].[12] 13X Zeolite PSA and VSA N/A Chaffee et al. [13] 13X Zeolite VSA N/A Schell et al. [14] Activated Carbon PSA Fixed bed Wang et al. [15] 13X Zeolite and activated carbon VPSA Fixed bed Pirngruber and Leinekugel-Le-Cocq [16] Supported polyamine PSA N/A Siriwardane et al. [17] NaOH/CaO TSA TGA and fixed bed Fauth et al. [18] Supported mixed amines TSA Fixed bed Gray et al. [19] Supported tertiary amine TSA TGA and fixed bed Gray et al. [20] Supported amine TSA TGA Lee et al. [21] Supported amine TSA TGA Sjostrom et al. [22] Supported amine TSA CFB Sjostrom and Krutka [23] Supported amine TSA Fixed bed Monazam et al. [24] Supported amine TSA Fluidized bed Postcombustion capture technologies include liquid solvent and solid sorbent systems [6].Since there is a large existing fleet of coal fired power plants, these technologies are particularly beneficial as they are well suited for retrofit installations.Furthermore, since they are added to the backend of the power plant they clean the exhaust of CO 2 without necessarily change the normal plant operation.The solid sorbent systems can be further divided into low and high temperature systems [7][8][9][10][11].
The absolute pressure of flue gas exiting a conventional pulverized coal-fired combustor is relatively low as well as the CO 2 partial pressure.Given this fact, the temperature swing option for adsorption/desorption appears to be the more viable option when amine-based solid sorbents are used [27].Also, the use of a solid substrate impregnated with and amine compound is a superior alternative to that of a water based solution of amine due to the lower specific heat of the solid substrate compared to an aqueous carrier.By moving the sorbent between regions of adsorption and desorption a more efficient system can be achieved.
Numerous reactor types have been employed and proposed for regenerable dry CO 2 capture systems.Even though investigated extensively at small scale, fixed bed reactors have considerable design issues when employed at the utility scale.The most significant of these being reactor heating and cooling as well gas ducting, piping, and valving.In contrast, fluidized bed systems are known for good heat transfer allowing for smaller reactors and have no issues associated with gas ducting, piping, or valving.The system investigated here is a circulating fluidized bed.This system has the advantage over a fixed bed system in that the sorbent can be moved from the adsorption region where exothermic reactions occur and heat must be removed to a regeneration region where the reaction is endothermic and heat must be added.These regions can be maintained at a fixed temperature in a circulating bed operation whereas in a fixed bed concept considerable energy must be added and removed from containment vessels.In a circulating bed, the only energy transfer required is for heating or cooling of the substrate, chemical reaction, and environmental losses.
Fluidized bed reactors have been used to evaluate adsorption and regeneration using these promising polyamine sorbents [22,24,28].Pilot studies using coupled fluidized beds reported a number of operational difficulties that were thought to be associated with desorption kinetics [22].Unfortunately, fluid bed desorption studies have been limited  to a very few tests over a narrow range of temperatures, CO 2 compositions, and levels of conversion [28].More information is available on the adsorption process.Kinetic analysis for fluid bed adsorption in these sorbents was consistent with the TGA data [24,29]; however, the proposed nucleation and growth rate expressions are too numerically unstable to incorporate into computationally intensive CFD models.
As a preliminary step in the investigation of a computational method to simulate the experimental system only the adsorption phase was considered in a batch phase.When a satisfactory result is obtained the simulation can be extended to include regeneration and full loop operation.
This work focuses on computational fluid dynamics (CFD) simulations and experiments of adsorption and regeneration using amine-based solid sorbent in a fluidized bed.This is part of ongoing investigations at NETL on CO 2 capture using solid-sorbents experiments and simulation [30][31][32][33].In this study, the simpler kinetic expression was employed to describe the CO 2 adsorption process.This study was undertaken to determine whether the sorption process could be adequately simulated using this simplified expression and how the reactor internals influenced mixing and the resulting CO 2 breakthrough profiles.

Solid Sorbent Carbon Capture System
A photograph of the experimental facility is shown along with a line drawing in Figure 2. The facility consists of an adsorber and riser on the right and the regenerator on the left.Sorbent particles under normal operation travel up the riser from the adsorber across the crossover and are separated from the CO 2 lean flue gas in the cyclone.The CO 2 loaded sorbent exits the cyclone bottom passing through the diverter valve and loop seal to the regenerator.The regenerated sorbent exits the bottom of the regenerator, passing through another loop seal and is returned to the adsorber.The diverter valve is used for short periods to measure the solids circulation rate.The specific batch experiments were conducted in the regenerator, the details of which are presented below.
The reactor (regenerator vessel) is comprised of a section of high density polyethylene 16.5 cm OD, 13.7 cm ID, 1.04 m long.A plenum was attached below the reactor where gases are introduced.The plenum as shown in Figure 3 is separated from the reactor by a sintered metal distributor plate which allows gas flow into the reactor.A 5.1 cm ID tube connects the regenerator to the bottom loop seal.This allows the solids to be withdrawn from the center of gas distributor when the unit is run in the full loop mode.For semicontinuous or batch tests this tube was sealed off by closing the valve below the plenum.Both N 2 and CO 2 flows were controlled using Alicat 0-50 slpm mass flow controllers.Humidification was accomplished using a Bronkhorst controlled evaporator mixer.Temperatures were measured using thermocouples at positions 1.0, 15.0, and 35 cm above the distributor plate.CO 2 concentrations were determined using a PP systems WMA-4 CO 2 analyzer (0-10,000 ppm range).The low range of the analyzer required the introduction of dilution gas (air) to the stream which was also controlled using an Alicat mass flow controller.
Heating/cooling is accomplished by flowing oil through copper coils immersed in the bed.When heating, the oil is passed through Cromolox heater where the oil temperature is controlled to yield the desired bed temperature.During cooling phases, the oil is passed through coils immersed in a cooled water bath then to the bed coils.There are two sets of coils, the outer coil consists of 1.27 cm OD copper tubing, 4.4 meters in length.The outer coil is a 1.57 cm diameter helix with vertical separation of 2.54 cm between the coils center to center.The inner coils are comprised of 0.95 cm OD copper tubing, 3.0 meters in length.The inner coil is a 5.08 cm diameter helix with vertical separation of 1.9 cm between the coils center to center.The sorbent mixture  is pretreated to remove any previously adsorbed CO 2 by exposing the bed to N 2 flow of 30 slpm and heating the bed to 110 ∘ C. When the CO 2 concentration dropped to negligible levels the bed was cooled to 70 ∘ C under humidified N 2 flow.After the desired bed temperature was achieved the desired mixture of CO 2 /N 2 flow was introduced.This flow was continued until the CO 2 flow at the outlet matched the inlet flow indicating the sorbent had reached its capacity.
The sorbent was then heated to 110 ∘ C under humidified N 2 flow and maintained at that temperature until the outlet CO 2 concentration was negligible.The bed was then cooled to the adsorption temperature and another adsorption condition was run.A summary of the experimental test conditions is presented in Tables 2 and 3.The flow rate information is presented in Table 1 while the adsorption/regeneration cycle information is presented in Table 2.

Particle Properties and Sorbent
Chemistry/Kinetics  silica substrate.A photomicrograph of this material is shown in Figure 4 using a Leica M205 FA microscope under white reflected light.FCC particles were smooth round particles of various shades of gray as a result of the coking in the petroleum refining process, while the sorbent particles were asymmetric larger white particles.The active ingredient in the sorbent is a polymer of 1,2-ethanediamine and aziridine which was dissolved into methanol, added to the silica substrate, and then the solvent was evaporated off.The substrate was a high surface area mesoporous silica (PQ Inc. 2129) with N 2 BET surface area of 312 m 2 /g [21].Initial experiments with 100% sorbent particles and relatively deep beds produced breakthrough times in excess of 1 hour, a time scale nearly impossible to simulate for CFD, since the simulations typically provide about 30 seconds of data per day.In this case, the simulations would need to run for more than six months to simulate a single case.Engineering analysis of the breakthrough data revealed that experiments could be designed to provide breakthrough times of less than 600 seconds by correctly choosing the sorbent bed material.The sorbent bed in these experiments and simulations consists of 80% FCC and 20% sorbent.The material properties are presented in Table 4.
The particle size was measured using QICPIC particle imaging analyzer (Sympatec GmbH, Model-QP0104).The particle size distribution for the mixture was Gaussian as displayed in Figure 5.The Sauter mean size of the FCC, sorbent AX and 80 : 20 mixture are presented in Table 4 along with other relevant hydrodynamic and CO 2 adsorption properties.The particle size of the FCC was slightly greater than that of the sorbent.This compensated for the FCC's higher density to produce similar hydrodynamic properties.
The minimum fluidization velocity was measured using a 5 cm diameter fluidized bed with porous sintered metal gas distributor.This velocity represented the intersection of two nominally linear regions (packed bed and fluid bed) in the plot of pressure drop across the granular bed against gas velocity (Figure 6).The pressure drop was normalized by the weight of the fluid bed.These materials were classified as Geldart Group A materials based upon their behavior and the relative densities and particle sizes.The sorbent CO 2 uptake behavior was measured in fluid bed [26].Sorbent capacities for CO 2 uptake were between 86 and 122 mg/g sorbent depending on temperature and, to a lesser extent, CO 2 concentration.The 20 : 80 mixture of sorbent and FCC resulted in a molar capacity of 0.54 moles CO 2 /kg sorbent which along with short sorbent beds facilitated timely CFD simulations.
3.2.Chemistry/Kinetics.Abbasi, Choi, Fauth, Gray, and Sjostrom [18-20, 22, 23, 26, 34] have all considered mesoporous materials for sorbents in carbon capture systems The work presented in this paper is based on a silica-supported amine sorbent material developed at NETL.This material is very similar to the materials described by both Gray et al. [19,20] and Fauth et al. [18].The CO 2 adsorption process is modeled utilizing mechanism that is similar to the one proposed by Lee et al. [21] for a similar sorbent for a dry reaction.The CO 2 absorption process is expressed as the reaction of gaseous CO 2 with two amine (R 2 NH) sites to produce a carbamate ion (R 2 NCO 2 − ) and protonated amine ion (R 2 NH 2 + ): This work utilizes the "volume-average" reaction model in Barracuda for the simulations presented in this paper.The adsorption (forward) and desorption (reverse) rates (kmol-CO 2 /m 3 -s) take the form where   is the total mass amine in a given cell and   is the mass of the carbamate ion.It was assumed in the derivation of the above expressions that the molar concentration of the protonated amine ion is equal to the concentration of the carbamate ion.The values of the coefficients used in the simulations are provided in the rate expressions are provided Breault and Huckaby [32].The values of  and  0 were tuned so that the predicted CO 2 in the exhaust matched the experimental measurement at the baseline condition as reported by Sjostrom et al. [22] while holding the  values equal to those presented by Lee et al. [21] after conversion of the units.Lee et al. [21] describe the chemical kinetics of CO 2 adsorption on PEI impregnated mesoporous silica sorbents in more detail.The silica substrate is considered an open sorbent with micropores and a network of tubular shaped mesopores which consist of silica-PEI composite structures depicted in the right hand close up of the PEI loaded mesopores.Monazam et al. [29] found that a mesoporous PEI based sorbent exhibits 2-dimensional nucleation growth diffusion controlled rates within this composite structure.The adsorption rate is controlled by the transport of CO 2 into the deeper PEI layers in the mesopores.In addition, it was demonstrated that the bulk and macroporous diffusion is fast compared to the sorption rate within the mesopores.However, when analyzing fluidized beds CO 2 breakthrough results Monazam et al. [24] found that the effective rate could be adequately modeled using first order Arrhenius kinetics with a negative temperature dependence providing confidence that the above kinetic models can reasonably be used in CFD simulations to predict CO 2 capture.

Numerical Model
Barracuda [35] is a CFD code developed by CPFD Software which was used to perform the simulations in this study.Barracuda is based upon the multiphase particle-in-cell (MP-PIC) model developed by Andrews and O'Rourke [36] and Snider [37].Barracuda is a discrete particle simulation tool (or Eulerian-Lagrangian) tool [33], where the particles are treated in a Lagrangian manner while the fluid conservation equations are solved using Eulerian equations.The effects of particle-particle collisions in MP-PIC are approximated on the gas continuum by utilizing a pressure proportional to the particle concentration.Furthermore, the parcel or cloud concept is introduced to allow the simulation of systems with large numbers of particles.Parcels or clouds are groups of particles which are assumed to have the same position, velocity, and properties.Therefore, the dynamics of any particle in the group can be derived from the equations of one particle in the parcel or cloud.Parcels or clouds are also referred to as computational particles and computational parcels.
O'Rourke and Snider [38] made some additional improvements to the momentum balance of the particle phase by adding a particle-particle "drag" term.This term is in addition to the typical "gas-solid" drag and the granular pressure terms.The form of the mass, species mass, momentum, and energy equations for the fluid-phase has been given by Snider et al. [39] and Snider and Banerjee [40].These equations were derived from Anderson and Jackson [41] and Jackson [42] and are discussed below.
CPFD's Barracuda code solves the fluid and particle equations in three dimensions using the averaged Navier-Stokes equations.Strong coupling of the particle phase with the fluid phase is obtained in the code utilizing the concepts developed by Snider et al. [39] and Snider and Banerjee [40].The particle momentum equation uses the multiphase particle-in-cell (MP-PIC) concept as initially developed by Andrews and O'Rourke [36] including modifications to it by Snider [37].That work was further refined by O'Rourke and Snider [38] to include a "relaxation-to-the-mean" term to represent damping of particle velocity fluctuations due to particle collisions.Mass, momentum, and energy of the twophase mixture are conserved by exchange terms in the gas phase mass, momentum, and energy equations, respectively.Therefore, in summary, the mass and momentum equations for the fluid-phase are averaged forms of the detailed fluidphase mass and momentum equations [39] and [40] and are presented in (3)-( 18) as follows.
Equations Used in Barracuda.Consider In these equations,   is the fluid volume fraction,   is fluid density,   is the fluid velocity vector,  ṁ is the gas mass production rate per volume from gas-particle chemistry,  is the mean flow gas thermodynamic pressure,  is the inter-phase momentum transfer rate per unit volume,  is the gravitational acceleration, and   is the fluid stress tensor.
The equation for the fluid stress tensor,   , is presented in (3)- (18) where  and  are direction and  is shear viscosity, which is the sum of the laminar shear viscosity and a turbulence viscosity.Large eddies are calculated with the unresolved subgrid turbulence being modeled with an eddyviscosity also presented in ( 3)-( 18) and identified as (6).In this equation,  is the Smagorinsky [43] coefficient (0.01) and Δ is the subgrid length scale which is the cube root of the sum of the product of orthogonal distances across a calculation cell.
The transport equations are solved for each gas species with the total fluid phase properties calculated from the mass fractions  , of the gas species .Mass is transferred between gas species through the breaking and forming of chemical bonds with a chemical source terms  ṁ,chem in the individual gas species transport equations (see (3)-( 18)) where  is the turbulent mass diffusivity and is related to the viscosity by the definition of the Schmidt number.
The mass balance and momentum balances for the gas mixture are presented in (3) and (4), respectively.The gas mixture properties are based on the mass fractions of the gas species, solved for using (7).Moreover, the flow is compressible, and the gas phase pressure, density, and mass fractions are related through the ideal gas equation of state.Using the ideal gas equation, the partial pressure of gas species  is presented as (9) where  is the universal gas constant,   is the gas mixture temperature, and Mw i is the molecular weight of gas species  and the total mean flow gas thermodynamic pressure is obtained by (10).
In Barracuda, the dynamics of the particle phase are predicted by solving a transport equation for the particle distribution function (PDF), .This function is based upon the work of [36] in their development of the MP-PIC methodology.Barracuda assumes that the particle distribution function is a function of particle spatial location (  ), particle velocity (  ), particle mass (  ), and particle temperature (  ), and time () is represented by (11).The average number of particles per unit volume,   , is obtained from the velocities in the interval (  ,   + du p ), masses in the interval (  ,   + dm p ), and temperature in the interval (  ,   + dT p ) as shown in (12).
In CPFD's Barracuda, the fluid mass source term from ( 3) is presented as (13) where the time-rate-of-change of particle mass dm p /dt is the rate of change of the particle mass producing gasses through chemistry.Furthermore, the acceleration on a particle is obtained from (14) where   is the solids volume fraction;   is the solid material mass density;   is the drag function which depends on the particle size, velocity, position, and time;   is the solids contact stress, which depends on spatial location;   is the local mass-averaged particle velocity; and   is a particle collision damping time from [39].
The solids volume fraction is related to the particle distribution function, , as shown in (15) where the fluid volume fraction is given by (16).
The interphase momentum transfer rate per unit volume, , in (4) is found from (17) and, finally, (18) gives the solid displacement.Collectively, these equations and relationships provide the backbone to CPFD's Barracuda code that was used in this work.

Experimental Results.
Five experiments were conducted with the sorbent-FCC mixture.The results from these tests are shown in Figure 7, where the instantaneous outlet concentration normalized by the inlet concentration is plotted against the elapsed time.The breakthrough curves are "classic" S-shaped sigmoidal curves.As expected, the breakthrough time increased with decreasing experimental flow rate.The tests at 30 slpm and 15 slpm were duplicated.There was little difference between the two test results.Using this later to define the test reproducibility, the average absolute percent error (AA%E) is 3.41 and the correlation coefficient between these two tests gives an R 2 value of 0.936.
The measured pressure drop across the fluid beds was 0.53 ± 0.009 kPa and did not vary significantly as the flow varied over the tests described in Table 2.This reflects the fact that all of these experiments were conducted at conditions well over minimum fluidization.The pressure drop represented less than 60% of the weight of the FCC-sorbent loaded into the bed.This result indicated that either channeling or defluidization occurred which caused a significant portion of the mixture to be supported by the internals, that is, the heat exchange coils, or defluidized in the central region of the bed above the stagnant solids underflow piping.The CO 2 sorption capacity of the sorbent was not adversely affected by the channeling or defluidization suggesting that the FCC was the preferentially segregated.The sorbent capacities were found to be 2.61 + 0.22/−0.31moles/kg of sorbent by integrating the CO 2 adsorbed over the duration of each test.The sorbent capacity did not vary substantially with flow.

Baseline Studies.
The simulation results are shown in Figures , and 8(c) along with the respective experimental data for the three test conditions.In Figure 8(a), CO 2 breakthrough was complete by about 900 seconds for the 15 slpm tests.The AA%E and the correlation coefficient R 2 between simulation and the average experimental value is 8.88 and 0.941, respectively.Recall from above that the error between experimental repeats had an AA%E of 3.41 and an R 2 of 0.99.The error between the simulation and the experiment is smaller than the error between the experimental repeatsthat is the best that can be expected.These error values are summarized in Table 5.
The experimental data for the cases using a flow rate of 30 slpm is compared to the simulation data for that condition in Figure 8(b).For this condition, the simulation results slightly overpredicted the experimental performance at times less than about 300 seconds and under predicted the experimental performance at times greater than about 300 seconds.This over-and underprediction tends to average out.The AA%E and the R 2 values are both better than the experimental error.Figure 8(c) shows the same information for an inlet flow rate of 60 slpm.For this case, the simulation under predicted the experimental performance across the entire time domain.However, the error between the experiment and the prediction is still better than the experimental error as can be seen in ( 3)- (18).The systematic deviation indicated that the kinetic model included a rate expression that was somewhat too high; however, the comparable duration to breakthrough indicated that the CO 2 adsorption capacity was quite accurate.Rather than use the CFD model to adjust the imperfect rate expression, it was decided to conduct a parametric analysis on the rates and use the best rate parameters to evaluate the hydrodynamics and mixing.
Another way to compare the model is to calculate and compare the total CO 2 adsorbed for the simulations and the experiments.This quantity is calculated by subtracting the values in Figure 8 from a value of 1 and then integrating to determine the area under the curves.These areas are shown in Figure 9 as a function of the feed gas flow rate.For the flow rates of 15 and 30 slpm, the simulation results were bounded by the experimental values.The deviation between the experimental value and the simulation value is within the error range as defined by the data spread for the 15 slpm case.Converting these values to a capacity and averaging the experimental values and the simulation values gives an experimental capacity of 2.61 + 0.22/−0.31moles/kg of sorbent and a simulation capacity of 2.61 + 0.11/−0.07moles/kg of sorbent.These values compare with the theoretical material property of 2.7 ± 0.1 moles/kg of sorbent.

Parameter Sensitivity.
In this discussion, three different aspects of the CFD model are explored to better understand the behavior of the model: reaction rate, particle size, and distributor plate design.This will be useful in future studies CFD simulations were performed using a value of the kinetic rate preexponential value.The results of these simulations are presented in Figure 10.In this figure, the same information presented in Figure 8(b) is accompanied by simulations where the Arrhenius factor was modified by first doubling it halving it.Halving this factor reduced the to point variance from an absolute average value of 5.3% to 3.7%.Doubling the preexponential factor increased the point to point variance an absolute average value of 6.1%.Since halving the rate improved the point point variance, a simulation was run at one tenth the rate to a minimum boundary.The point point average variance for the simulation case with the rate taken as one tenth the fitted rate is 4.95%.The rate is too low as the curve at low conversions was to the right of the data and with this slower rate is to left the experimental data as can be seen in Figure 10.Looking at the simulation results and the experimental data in Figure 10, it that the rate obtained from the ADA performance [22] is between four and six times too fast for the sorbent used in these studies.This difference may partially be to external mass transfer effects.
The effect of particle has been explored and the results presented Figure 11 with respect to the baseline.In this figure, the normalized exit CO 2 in the stream is plotted against time as was done in previous plots (Figures 8 and  10).The CO 2 breakthrough for the simulated case using 300 m particles is compared experimental data and the base case simulation results using the same symbol and pattern as used in Figure 10.The 300 m simulation results produce a steeper curve.According to the analysis presented in Figure 10, the rate of removal for the larger particles faster than that for the smaller base case particles.This is because the reaction rate is independent of the size, thus the difference in the bed state.Hydrodynamics results are presented in Figures 12 13.These figures from left to right present the fraction, the gas velocity, and the species and without the diluents, The base case with smaller particles (Figure 12) has a bed that is about 20% more expanded than that for larger particles (Figure 13).This is consistent with the excess gas velocity.The in Figure 12 was obtained at a level / value of about 16 whereas the data in Figure 14 was obtained at a / mf value of 1.5.The base case produced higher and nonuniform gas velocity, while the larger particles exhibited lower and more uniform gas velocity.Also, the distribution of the reacting particles for the larger particles is more uniform within the region of particles.In summary, these results that the hydrodynamics for larger lead to a delayed breakthrough of the CO with a more compressed adsorption zone which less CO 2 reactivity in both the leading and trailing edges of that zone.
In the unit, the fluidizing gas is not uniformly injected across the entire cross section. is the of the solids withdrawal when operating the reactor in a continuous To assess the effect this might have on the simulated results, a simulation was where the gas was injected into the annular comparable to the experimental setup.The results from this simulation are in Figure along with the corresponding experimental data and the base simulation.The rate reaction for the modified distributor is than that for the base simulation.This is primarily due to the hydrodynamics in the modified simulation.Compared to the base the simulation with the modified distributor plate has more uniform distribution of the reactive species as seen by 15 and 12 together.Also, the gas flow distributed over the annular region where solids are located.These combined effects lead to a narrower reaction zone with less bypassing in the leading and trailing regions prior to and following breakthrough.The measured pressure drop across the fluidized bed was 0.63 kPa compared to the simulated value of 0.80 kPa.This represents only 58 and 70% of the weight of the mixture in the fluid beds.The simulation was able to capture some of the causes for this.In Figure 12 it is apparent that the FCC collected in the center stagnant region while the sorbent was concentrated in the annular regions of cylindrical bed.This type of segregation explains why capacity of the sorbent did not suffer from the stagnation in the center.The sorbent was apparently always available to absorb CO 2 from the flowing gas.The simulation not include the heat exchange coils and this may explain why the pressure-drop for the hydrodynamically supported solids was higher than that measured experimentally.

Conclusions
It has been that the Barracuda code by CPFD software can used to simulate the performance of a fluidized reacting system for the adsorption of CO 2 on supported polyamine sorbent material.The code accurately predicted the breakthrough time and the sorbent capacity within the error as defined by repeat experiments.The Barracuda model was able to capture the channeling and segregation of mixture a pressure drop that represented 70% of the weight of the bed, while the experimeasured a pressure drop of about 60% of the weight of the bed.However, the segregation did affect the capacity of the sorbent experimentally or in the CFD model indicating that the sorbent remained in the flow path.The theoretical capacity of the sorbent material was 2.7 ± 0.1 moles/kg of The experiments an average result within that experimental error of 2.61 + 0.22/−0.31.
The simulation results averaged to the same value as the experimental results, but with tighter deviation as defined by the range for the three conditions with the simulation results being 2.61 + 0.11/−0.07moles/kg of sorbent.Application of the model to other cases revealed that the predicted bed reactivity (or effective reaction rate) is about four to six times higher than the observed experimental rates.The parameter sensitivity studies suggest a number of paths for further improvement of the model, including the reformulation and subsequent calibration of a particle reaction model.This work has been started and should be concluded early next year.It was also shown that a nonuniform solids distribution and contacting gas velocity distribution results in a decrease in the size of the adsorption zone while reducing CO 2 bypassing.In comparing the experimental pressure drop across the reactor to the predicted from the simulations, it is concluded that the model should be expanded to explicitly resolve the internal heat transfer internal surfaces.The interaction of the heat transfer surfaces with the nonuniform fluidization will provide a channeling path for the gas, decreasing the overall bed reactivity.It is anticipated that further improvements to the overall predictive capability will be provided by including local temperature variations which will affect the reaction rate and CO 2 -absorption capacity.

Disclaimer
The U.S. Department of Energy, NETL contributions to this report were prepared as an account of work sponsored by an agency of the United States Government.Neither the United States Government nor any agency thereof, nor any of their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights.Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof.The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof.

Figure 3 :
Figure 3: Details of adsorption reactor used in tests.

Table 1 :
Summary of literature capture concepts.

Table 5 :
Summary of AA%E and  2 values.