Parameters of equations (1) - (4), (Gogová & Hanika (2009,b) ).
1. Introduction
Oxidations of organic compounds in liquid phase by oxygen have been applied for years in many important industrial and waste water treatment processes. There are several variants of technical design of these processes – ranging from homogeneous through heterogeneous to biotechnological routes. As an example, the process of gluconic acid production by glucose oxidation can be arranged in all of these variants.
Bioprocesses are usually carried out in aqueous media at ambient temperature and atmospheric pressure in the presence of living microorganisms and their enzymatic apparatus (e.g. Aspergillus niger), or by using pure enzymes (glucose oxidase and catalase), (Sikula, et al. (2006), (2007) ). In the former case, the biomass represents a solid phase in the reaction system, whereas in the latter case, the reaction system is homogeneous – liquid. Some drawbacks are also inherent with the bioprocesses – e.g. strong sensitivity of microorganisms to impurities present in the reaction system, losses of the substrate transformed to carbon dioxide or utilized for the microorganisms growth, low solubility of oxygen in the reaction system owing to the presence of ionic salts and nutrients (e.g. glucose) with a high rate of oxygen consumption by the microorganisms on the other hand, frequent occurrence of non-newtonian hydrodynamic properties of biomass suspension, foam formation, etc. For such bioprocesses, gas-lift reactors (known as air-lift reactors) (GLRs) are often used for their capability of delivering oxygen to the growing culture at a sufficient rate, while maintaining low shear stress.
Heterogeneous catalyst application to the oxidation process is advantageous compared to homogeneous catalytic systems with respect to simpler separation of a catalyst from the reaction mixture by filtration. At the heterogeneous variant however, a catalyst selection and its optimization is one of the crucial points to be considered. Another one is the reactor type selection and its design (estimation of the geometry and the size of the reactor selected). Neither the authors’ experience, nor literature search provide many generalizations on selection of a reactor type for which a counter-example could not be thought up. An open mind and good ideas are probably more important here than any generalization. Furthermore, because of complexity in scale dependency of various reactor selection criteria, the authors incline to agree with the statement of Bisio & Kabel (1985) :
Heterogeneously catalyzed wet air oxidation of glucose (Glc) to gluconic acid (Glcac) in aqueous alkaline solution serves as a model reaction. Palladium on activated carbon commercial catalyst enables to run the reaction selectively at ambient conditions. On an industrial scale, biotechnological routes of Glcac production currently prevail over the catalytic one. This is mainly because of the Glcac broad utilization in the food industry. The other reason is a problem with activity of the catalysts used. Pt-group catalysts suffer from gradual reversible deactivation due to an action of oxygen during the reaction course.
One way to overcome the problem leads through the catalyst optimization. Recently, good activity, selectivity and long-term stability were reported for supported gold catalysts (Biella, et al. (2002), Comotti, et al. (2006), Thielecke, et al. (2007) ).
Another approach to solve the problems with the catalyst activity deals with the process and/or reactor optimization. It is based on correct choice of a reactor type for a chemical process, the reactor well-suited design and on setting an appropriate mode of the reactor operation. These are the crucial aspects for maximizing the technological output.
The text is focused on solving the problem with the catalyst unstable activity through the reactor / reaction step optimization. Optimization of continuous stirred tank reactor (CSTR) and gas-lift reactor (GLR) productivity through the gas feed modulation is attempted. For any input operational conditions the task is to find conditions of the highest possible productivity of the reactors, i.e. to find conditions where the reaction and reactivation times are shared optimally, so that neither any time is wasted in prolonged activation process, nor is an insufficient activation time provided.
Beneficial effect of composition modulation on a CSTR performance is demonstrated. The catalyst activity can be maintained long-term steady by periodically alternating the gas feed composition. In the case of CSTR, period length and the period split represent independent variables. They both can be varied independently within one CSTR unit of a given construction. It is demonstrated here that for any period length always a split value exists, where the maximum reactor productivity is achieved. By connecting the points of optimal split for every period length, trajectory of the maximal CSTR productivity is obtained.
GLR was selected as the reactor type suitable to carry out the model reaction in. A GLR natural operation enables the catalyst to be periodically exposed to reaction and activation conditions, in the riser and downcomer sections of the GLR, respectively. Such reactors are in their nature multifunctional. In a GLR both, the period length and the split value are bound with geometry of the GLR given. Therefore only one geometrical optimum exists for given set of input operational conditions. The maximal GLR productivity is guaranteed only in this geometrical optimum, because the residence time in riser (reaction time) and the residence time in downcomer (activation time) are only here shared optimally.
1.1. The model reaction kinetics
Glucose (Glc) wet air oxidation over palladium on activated carbon catalyst is used as the model reaction. It takes place in three-phase medium. Advantage of the catalyst is its ability to catalyze the reaction selectively towards gluconic acid (Glcac) at mild conditions. Its drawback is fast deactivation during the oxidation reaction, when oxygen forms surface oxides or penetrates from the topmost Pd layer to subsurface layer forming subsurface oxide (Simmons, et al. (1991), Lundgren, et al. (2002), Ketteler, et al. (2005) ). These Pd-O phases are less active compared to chemisorbed oxygen (oxygen adsorbed above the first metallic layer), and will be referred to as those responsible for change in the catalyst activity. This deactivation was proved reversible (Vleeming, et al. (1997)).
Activation and reactivation of the catalyst is based on reduction of the catalyst active sites. Glucose is good reduction agent to pre-reduce the catalyst. Existence of optimal activation times depending on the reaction mixture composition was observed (see Gogová & Hanika (2009,a) for details).
Equations (1) and (2) form the kinetic model of the model reaction (Gogová & Hanika (2009,b) ). They describe mathematically processes of the main surface reaction, the catalyst deactivation and the catalyst reactivation. Change in the catalyst activity is described through a change in fractional coverage by inactive oxygen species,
This kinetic model applies all the time except when oxygen concentration in the liquid phase approaches zero. Then the reaction and reactivation mechanism changes: inactive oxygen species (responsible for the catalyst deactivation) take over the function of the chemisorbed oxygen in the main surface reaction. When the mechanism changes, its mathematical description also changes, and equations (3) and (4) apply as the kinetics model instead of equations (1) and (2).
Expressions of the lumped rate and adsorption constants of equations (1) - (4) are listed in Table 1.
Rate and adso rption constant | Dimension | Value |
k w K Glc K O K Glcac k D k A
|
[m 4.5 kg -1 mol -0.5 s -1] [m 3 mol -1] [m 1.5 mol -0.5] [m 3 mol -1] [m 1.5 mol -0.5 kg -1 s -1] [kg -1 s -1] [m 3 kg -1 s -1] | 0.00313 0.0169 4.50 0.384 0.00612 0.00518 5.47 10 -5 |
In the kinetic model (equations (1) and (2) or alternatively (3) and (4)), the change in the catalyst activity is expressed through the change in the fractional coverage by inactive oxygen species,
Relative activity of the catalyst at time
Thus the relative activity is useful parameter that characterizes changes in the reaction rate as the catalyst deactivates, and it is obtained conveniently from the experimental results. The equation (5) applies to all deactivation processes, no matter if the rate equation is separable or not according to the concept of separability (Szépe & Levenspiel (1970), Butt & Petersen (1988) ).
A rate equation is separable if it can be expressed as a product of two terms – the reaction rate on the fresh catalyst and the catalyst activity in the following form:
Active fraction
The rate equation (1) is separable. Combination of the equations (1), (6) and (7) gives the following relation between the catalyst activity and the active fraction:
Thus the activity depends only on the amount of inactive oxygen species.
1.2. Reversible deactivation of the Pd/C catalyst
Figure 1 provides an insight into the findings made during the model reaction kinetics study. Several regions can be recognized there. Before each experiment in semi-continuous stirred tank reactor (SSTR), the catalyst was activated in the reactor by its reduction with Glc as a component of the reaction mixture in inert atmosphere. The reaction was started up replacing nitrogen flow by flow of nitrogen/oxygen mixture with the desired partial pressure of oxygen. Each experiment consisted of one or more consecutive oxidation runs. Between these oxidation runs the catalyst was reactivated in inert atmosphere with Glc.
The reaction rate at the beginning of the second reaction cycle in SSTR (full circles in Figure 1) is lower because of change in the reaction mixture composition during the first reaction cycle in the batch system (SSTR). To prove full reversibility of the Pd/C catalyst deactivation, the primary experimental data in Figure 1 were corrected for the reaction mixture composition change (empty circles in Figure 1). To serve this purpose, the Glc and Glcac concentrations at the reaction start-up were applied in the kinetic model (equations (1) and (2) with parameters of Table 1).
Figure 1 indicates possibility of improving the reactor performance by periodically exchanging the reaction and reactivation cycles. In the text below, this approach is analyzed in process of selection and optimization of a target reactor suitable to carry out the model reaction in.
1.3. Strategy for elimination of the catalyst deactivation
The advantage of the Pd/C catalyst is its ability to catalyze the model reaction efficiently at mild conditions maintaining high selectivity towards Glcac. Its drawback is fast deactivation during the oxidation reaction. In one hour the catalyst activity can drop to less than 40% of its original value depending on the reaction conditions. Although reversible, the deactivation rate presents a crucial problem for industrial implementation of the process.
This text is devoted to one of many strategies aimed to eliminate the problems with unstable activity of the catalyst. It leads through the process and/or reactor optimization. This approach deals with correct choice of a reactor type for a chemical process, the reactor well-suited design and with setting an appropriate mode of its operation. These are the crucial aspects for maximizing the technological output.
For process similar to the model one, Markusse, et al. (2001) found that the catalyst activity can be maintained steady by periodically switching between oxygen and nitrogen flow to a CSTR. In general, the term “periodic operation” refers to operation regimes in which one or more reactor parameters vary in time. Modulation of mostly composition and/or feed flow rate was researched by e.g.: Boelhouwer, et al. (2002), Silveston & Hanika (2002), Tukač, et al. (2003), Silveston & Hanika (2004), Liu, et al. (2008) etc., with the aim to improve chemical reactors performance through forcing the reactor to operate under transient rather than steady-state conditions. Silveston (1998) in his monograph pays attention to several catalytic processes operated in this way.
For the reaction of Glc oxidation by air over reversibly deactivating Pd catalyst, it was indicated in Figure 1 that CSTR productivity can be enhanced by the gas feed flow modulation. Beneficial effects of composition modulation on a CSTR performance were studied in Gogová & Hanika (2009,b) with the model reaction. The task was to share optimally the reaction and reactivation times within given period length, so that neither any time is wasted in prolonged activation process, nor an insufficient activation time is provided.
2. Dynamic operation of CSTR with feed periodic modulation
Possibility of improving the reactor performance by exchanging the reaction and reactivation cycles is indicated in Figure 1. For deeper insight into the model system behaviour under periodic mode of operation, the kinetic model (equations (1) and (2); or (3) and (4)) was implemented in mathematical model of a CSTR (equations 9-12) operating at constant Glc and Glcac concentrations in time and with varying volumetric flow rate of the liquid feed stream according to the value of immediate reaction rate,
with initial conditions:
where
The value of
Oxygen saturation concentration in the liquid phase,
The kinetics model (equations (1) and (2) or alternatively (3) and (4)) embedded in model of CSTR enables to separate the effect of the reagents concentrations from the effect of the change in the catalyst activity itself on deviation of the reaction rate in time. The CSTR model makes it possible to express the catalyst activity directly through the ratio of the immediate to the initial reaction rate (i.e. by using equation (5)), and reveals directly the progress in the extent of the catalyst deactivation, see Figure 3. Figure 2 illustrates the effect of varying oxygen partial pressure in the gas feed stream on performance of the CSTR operated under O2/N2 periodic mode.
Figure 2a shows the rate of Glcac formation in time and Figure 2b reveals the progress in the catalyst fractional coverage by the inactive oxygen species, which stand behind the catalyst deactivation. It can be seen that in addition to the observed inhibition effect of oxygen on the Glc oxidation rate (Vleeming, et al. (1997), Gogová & Hanika (2009,a), (2009,b) ), also the rate and the extent of the catalyst deactivation are influenced by oxygen concentration in the liquid phase. With its raise, both the
2.1. Optimization of the CSTR under forced periodic operation
The CSTR operation was optimized in the conditions outlined above with equations (9) – (12), i.e. at constant Glc and Glcac concentrations in time and varying volumetric flow rate of the liquid feed stream according to the immediate reaction rate. The reactor is now operated under on-off periodic mode, i.e. with alternating cycles of switching on and off air feed stream. The catalyst reactivation in this case takes over in the oxygen-free intervals.
The period length is defined as a sum of reaction and reactivation times:
The split of period is the time the reaction takes in relation to the entire period length:
The reactor productivity is represented by cycle-time-averaged reaction rate,
In case of CSTR, period length and the period split value represent independent variables. They both can be varied independently within one CSTR unit of a given (and constant) construction. Beneficial effect of the gas feed modulation on the CSTR performance is showed in Figure 4.
The catalyst activity can be maintained long-term steady by periodically alternating the reaction and activation periods of the catalyst operation. As can be seen in Figure 4, for any period length always a split value exists, where the maximal reaction rate is achieved. The CSTR optimization task was to find conditions that guarantee the highest reactor productivity at any period given. In other words, the optimal reaction-reactivation time-share had to be found within any given period length. The maximal CSTR productivity is only guaranteed in the split optimum, where no time is wasted in prolonged activation process, neither an insufficient activation time is provided. Figure 4 maps the CSTR performance for selected input conditions. By connecting the points of optimal split for every period length, trajectory of the maximal CSTR productivity is obtained. For illustration, the trajectory is highlighted in Figure 4.
In the direction of decreasing the values of
approaches value of steady-state reaction rate of such system. It should be noted that the model solution in Figure 4 ought to be taken as an approximation in the region near
3. Multifunctional gas-lift reactor (GLR) employment
3.1. Characteristics of gas-lift reactors
What makes gas-lift reactors attractive for chemical and biotechnological applications is their relatively simple construction with possible segregation into various reaction zones, low and homogeneously distributed shear forces, good (and cheap) mixing with elimination of backmixing, and whole lots of possible design modifications.
Gas-lift reactor (GLR) consists of four main sections (see Figure 5): riser, gas-liquid separator, downcomer and bottom of the reactor. Operation of the GLR is relatively simple. It is based on spontaneous circulation of the reaction mixture along these four sections of the reactor as a result of difference in apparent density of the media present in riser and downcomer of the GLR. Successful application of these reactors to a specific (bio)chemical process is closely related with proper design of the reactor and on optimization of the mode of its operation.
However, in a GLR of a given construction (geometry), superficial gas velocity is the only independent variable that affects the entire hydrodynamics within the reactor - see the scheme in Figure 5, which documents the complexity of phenomena that occur in a GLR. Understanding of such hydrodynamic phenomena as gas hold-up, flow regimes or circulation velocity leads to more insight into the resulting mixing, heat and mass transfer.
When a GLR is employed for Glc production with living cultures, then quite contrary to the catalytic route, oxygen has to be present in riser, as well as in downcomer, to ensure living conditions all over the GLR circulation loop. But since the difference in the gas hold-up in riser and that in downcomer (
The geometric design of GLR (scale of reactor, separator design, slightness of the reactor, ratio of cross-sectional areas of the downcomer and the riser etc.), the superficial gas velocity, pressure drop (friction) along the flow path and physical properties of the liquid phase have a strong influence on both the gas hold-up and the liquid velocity.
In a gas-lift (air-lift) reactor the hydrodynamics, transport and mixing properties, gas hold-up, interfacial areas and interphase mass transfer coefficients depend strongly on the prevailing flow regime. Following regimes occur in direction of increasing flow rate and can be deducted from both visual observation and gas hold-up:
3.2. Natural periodic operation of GLR
Gas-lift reactor (GLR) was selected as the reactor type suitable to carry out the model reaction in. GLR natural operation resembles the above mentioned forced periodic operation of the CSTR as follows: In GLR, the model reaction, as well as the catalyst reactivation proceeds within one multifunctional reactor unit. In principle, GLR operation is based on spontaneous circulation of the catalyst dispersion in liquid, which results from difference in apparent density of the reaction mixture present in riser and downcomer sections of the reactor. Therefore, if complete separation of the gas phase is ensured after the reaction media passes the riser where the main reaction (and the catalyst deactivation) takes part, the downcomer serves as reactivation zone of the reactor.
Kluytmans, et al. (2003) proposed application of GLR for reaction similar to the model one. To design the target three-phase GLR that would suit the reaction of Glc oxidation over the Pd/C catalyst, we constructed GLR mathematical model (see Gogová & Hanika (2009,c) ) that is presented and applied in the following text. For any input operational conditions maximal productivity of the reactor is only guaranteed in the point of the target GLR optimal geometry. The GLR model proposed employs tanks-in-series mixing model sketched in Figure 6 to combine hydrodynamics of a real GLR (Bello, et al. (1984), Blažej, et al. (2004), Blažej, et al. (2004), Juraščík, et al. (2006), Sikula, et al. (2007), Sikula (2008), Sikula & Markoš (2008) ) and the kinetics of the model reaction from Section 1.1. This GLR model is much simpler and a bit more realistic than that of Kluytmans, et al. (2003) who employed axial dispersion mixing model with hydrodynamics measured in small-scale 2D bubble-column reactor and considered intra-particle diffusion.
3.3. Modelling and optimization of the GLR productivity
Optimization of GLR productivity is attempted in the following couple of sections. In the case of GLR it is impossible to move along a trajectory of maximal productivity as was the case with the CSTR of Section 2.1, without reconstruction of the GLR itself, as both, the period length and the split value, are bound with geometry of the GLR given. Since both, the reaction and the reactivation times in GLR are set by the reaction mixture residence times in riser and downcomer of the GLR, respectively, only one geometrical optimum exists for given set of input conditions. The maximal GLR productivity is guaranteed only in this geometrical optimum. The optimization task in this case therefore is to find these geometrical optima for any set of input conditions.
GLR mathematical model was derived and applied to aid the target reactor design. The GLR model consists of two main parts (Figure 7). In the first one (hydrodynamics cycle), hydrodynamics and optimal geometry of the reactor is iteratively calculated. The second part (reactor performance cycle) uses the results of the first part, links them up with the model reaction kinetics and iteratively calculates the actual GLR steady-state performance.
Tanks-in-series model (as sketched in Figure 6) is employed in the second part of the GLR mathematical model to grasp the way of mixing in real GLR. Every tank of the tanks-in-series model is described by a set of nonlinear algebraic equations (NAE) linking hydrodynamics and kinetics that apply in the given section of GLR. Tanks-in-series and axial dispersion models are the most frequently used mixing models for GLRs. Both were applied recently for simulation of the biotechnological equivalent of the model reaction run in GLR (Znad, et al. (2004), Sikula, et al. (2006), (2007), Sikula & Markoš (2008) ).
The following set of assumptions applies with the derived GLR mathematical model:
The GLR operates in steady state in the region of homogeneous bubbly flow.
The downcomer gas hold-up is zero; this is achievable by correct design of separator (Gogová, et al. (2002) ).
The number of tanks that form riser and downcomer corresponds to the extent of axial dispersion in these sections. Plus, there are two tanks for each – the bottom and the separator. In the latter two tandems, the sizes of the individual twin tanks vary depending on the respective volumes taken up by either liquid (then they become a functional part of downcomer) or G-L dispersion (and then act as a part of riser).
Each of the tanks in series operates isothermally and is perfectly mixed.
The mathematical model combines description of the GLR hydrodynamics with the kinetics of the glucose oxidation reaction and the catalyst deactivation and reactivation.
The reactor operates at low Glc conversion, up to 10%. Then, according to Kunz & Recker (1995) assumption of 100% selectivity towards Glcac can be taken.
The Pd/C catalyst particles (45 µm) are assumed to be homogeneously dispersed in liquid phase. Then it is justifiable, for the mathematical description purpose, to accept a concept of pseudo-homogeneous reaction phase. However, apart from the liquid phase, the catalyst does not leave the reactor.
Constant liquid volumetric flow rate within the reactor sections is assumed, i.e. the liquid density is constant – independent on the conversion degree.
3.3.1. Algorithm for the GLR mathematical model solution
The GLR model is used to firstly aid the target GLR design for any given input conditions, and secondly, to predict steady state characteristics of the target reactor. The simulations presented were run by using commercial software Matlab.
Algorithm of the GLR mathematical model is sketched in Figure 7. The GLR model input adjustable parameters are: degree of Glc conversion, all the feed (and start-up) concentra-tions, the catalyst concentration, kinetics parameters, temperature (fixed, 30 C), atmospheric pressure (fixed), superficial gas velocity, number of tanks and basic geometrical parameters. The basic geometrical parameters (the reactor type – internal-loop gas-lift reactor (ILGLR), the GLR volume, its separator volume, its bottom height, the liquid height, the outer column diameter and the riser wall thickness) were adopted from experimental GLR, volume 40L (see Blažej, et al. (2004), Blažej, et al. (2004), Juraščík, et al. (2006), Sikula, et al. (2007), Sikula (2008), Sikula & Markoš (2008) ). This particular reactor has been a source of many hydrodynamic correlations employed in the GLR model. Non-adjustable parameters have complex dependencies on the input adjustable ones. The procedure pictured in Figure 7 is applied to calculate them. See Gogová & Hanika (2009,c) for more details on the GLR mathematical model.
The one-tank model is an integral part of the GLR mathematical model (see Figure 7). It uses concept of one CSTR, which for the period of
The one-tank model serves to estimate a first approximation of the optimal split value (defined in Eq.(15)) for the target gas-lift reactor. At the optimal split, the target GLR productivity reaches its maximum (given by the cycle-time-averaged reaction rate). The
3.4. Optimal design of the target GLR
To optimize GLR productivity, optimal residence times in the riser and downcomer sections have to be found. Optimal value of split,
Case |
|
|
A D /A R [-] | D R [m] | D D,Eqv. [m] |
a) |
0.00286 | 0.0658 | 1.6989 | 0.0909 | 0.1184 |
b) |
0.00349 | 0.1332 | 0.8588 | 0.1084 | 0.1005 |
c) |
0.00287 | 0.2380 | 0.4065 | 0.1236 | 0.0788 |
Figure 8 shows effect of superficial gas velocity
It can be seen in Fig. 8, that a change in
It was found during the model reaction kinetics study that the reaction rate is inhibited by oxygen. Moreover, oxygen affects the extent of the catalyst deactivation. The variation range of oxygen content in gas feed stream is therefore limited. Figure 8 covers major part of the target GLR operational window in terms of
For the input operational parameters given by the points on the
As proved above, depending on the input conditions the optimum reaction-reactivation time-share varies and so does the optimum in the target GLR geometry, which is also reflected in the profiles of the reaction rate and the concentrations along the GLR. In Figure 10, calculated profiles of the actual concentrations and reaction rates along the circulation loop are presented for selected set of input parameters. Each symbol in Figure 10 represents one tank of the tanks-in-series within the loop (their abbreviations are for illustration marked at the top of the Figures 10a, b). The space time in Figure 10 is defined as follows:
The target GLR operates continuously. In the profiles, inlet and outlet points in the reactor are visible (compare with Figure 6). For the input conditions listed with Figure 10, the calculated reactor productivity (cycle-time-averaged reaction rate) is 3.49 mmol Glcac per 1 kg of the catalyst per second.
In Figures 11 and 12 calculated profiles of reaction rates and molar fraction of oxygen in the gas phase are shown along the GLR circulation loop. The arrows indicate the trends as the
Similarly to Figure 10, the maximum immediate reaction rate is achieved in separator(D) tank for every
Simulation results in Figure 12 show profiles of the immediate reaction rates and molar fraction of oxygen in the gas phase for various
3.5. Practical aspects
The derived GLR mathematical model was used for computer aided design and optimization of a target multifunctional gas-lift reactor with the aim to solve the main problem of the model reaction - the catalyst unstable activity. Examples of such processes are e.g. wet air oxidation of waste waters or syntheses of chemical specialties. Glucose oxidation (in alkaline aqueous solution) with oxygen in the presence of a palladium catalyst was used as the case study. Application of GLR to this reaction system enabled simulta-neous reaction in riser and the catalyst reactivation in downcomer section of the reactor.
The GLR model assumptions lean on the kinetics of the reaction and the catalyst deacti-vation. All the simulations assume homogeneous bubbly flow of the reaction mixture in riser, zero gas hold-up in downcomer and predict the system steady state operation. The isothermal tanks-in-series mixing model describes axial dispersion of the reaction mixture in the oxidation and the catalyst reactivation sections of the reactor. Hydrodynamic parameters of the GLR model were taken from pilot plant data (Blažej, et al. (2004), Blažej, et al. (2004), Juraščík, et al. (2006), Sikula, et al. (2007), Sikula (2008), Sikula & Markoš (2008)).
The GLR mathematical model proposed helps to overcome the problem with the catalyst unstable activity by appropriate calculation of the target GLR geometry for any given input operational conditions. Moreover, the model is capable of predicting optimal geometry of the target GLR, its maximal productivity and other steady state characteristics for reactions similar to the model one, i.e. for G-L-S oxidations with reversible deactivation of a catalyst due to an action of any substance present in the gas phase. The limitations are only given by meeting the ranges of GLR operational window as explained with Figures 8 and 9.
The GLR model derived is not limited to the model reaction only. It can easily be extended for other G-L-S oxidations with reversible deactivation of a catalyst due to an action of a substance present in the gas phase. The limitation here is given by at least partly overlapping the ranges of a GLR applicability (see the reasoning given with Figures 8 and 9 about the ranges of the model reaction operational window) with the new reaction requirements (set by the new reaction kinetics).
The proposed GLR model can also be extended to non-isothermal process conditions. In future, the area of the GLR model employment may be broadened for process scale-up and the reactor safe control. But, the experimental validation of the model solutions remains a challenge for the future research.
For the biotechnological routes of Glcac production, GLRs are also of interest due to several advantages that they offer over alternative bioreactors. However, a GLR design for the biotechnological applications is based on different policy, compared to the catalytic application. In the catalytic process the condition of zero gas hold-up in downcomer was directive for the reactor design. On the contrary, in the biotechnological application this condition is no longer relevant as the living conditions for the cell cultures involved have to be ensured all over the GLR circulation loop, i.e. oxygen has to be present in downcomer. Therefore, GLR correctly designed for a bioprocess should provide sufficient G-L mass transfer, and it should operate at an optimal gas hold-up in downcomer.
4. Conclusion
The model reaction chosen appears to be interesting, not only due to its versatility for the industrial applications, but also because it raises many chemico–engineering problems to be solved in unconventional ways. The main problem of the model reaction is fast but fully reversible deactivation of Pd catalyst due to an action of oxygen during the reaction course. Various approaches can be taken to get over the problem. In the work described in this chapter, the problem is tackled through a tailored selection, design and optimization of the catalytic reactor. For reasons explained below, gas-lift reactor (GLR) was selected as the target reactor suitable to carry out the model reaction in. It allows the reaction along with the catalyst reactivation proceed within one reactor unit. Such reactors are in their nature multifunctional.
Deeper insight into the catalyst deactivation is made by analyzing the model reaction behaviour under conditions of a continuous stirred tank reactor (CSTR) operation. Optimization of the CSTR productivity through the gas feed stream composition / flow modulation was attempted. Dynamic mathematical model of CSTR operating under on-off periodic mode was used to aid this task. In the periodic operation, enhancement of the overall reaction rate results from forcing the catalyst to operate under transient conditions. In CSTR, period length, as well as the period split value represent independent variables and can be varied independently within one CSTR unit of a given (and constant) construction. For any period length always a split value exists, where maximum reaction rate is achieved (see Figure 4). The optimization task was to find conditions that guarantee the highest reactor productivity at any period given. In other words, the optimal reaction-reactivation time-share had to be found within a given period length, where the maximum CSTR productivity is guaranteed, because in the period split optimum no time is wasted in prolonged activation process, neither insufficient activation time is provided. Figure 4 maps the CSTR performance for selected input conditions. Benefits of running the model reaction under conditions of periodic exposure to oxidative and reductive environment were explored and a trajectory of maximal CSTR productivity was defined. A real production plant operational requirements and limitations decide about the position on this trajectory.
GLR (gas-lift reactor) natural operation resembles the above mentioned forced periodic operation of the CSTR as follows: In GLR, the main reaction, as well as the catalyst reactivation proceeds within one multifunctional reactor unit. In principle, GLR operation is based on spontaneous circulation of the catalyst dispersion in liquid, which results from difference in apparent density of the reaction mixture present in riser and downcomer sections of the reactor. Therefore, if complete separation of the gas phase is ensured after the reaction media passes the riser where the main reaction (and the catalyst deactivation) takes part, the downcomer serves as reactivation zone of the reactor. If economic aspects were considered, GLR natural periodic operation would be cheaper than the forced periodic operation of CSTR.
In GLR it is impossible to move along a similar trajectory of the highest reactor productivity (as was the case with CSTR) without reconstruction of the GLR itself. The period value is given by liquid circulation velocity, i.e. the time one circulation loop takes; and the split value is set by the given GLR geometry. Only one geometrical optimum exists for given set of input operational conditions. The maximal GLR productivity is only guaranteed in this geometrical optimum, because the residence time in riser (reaction time) and the residence time in downcomer (activation time) are only here shared optimally. The optimization task for this GLR case was to find the optimal geometry for any set of input conditions. The derived GLR mathematical model was used to aid design of target reactor for reaction of heterogeneously catalyzed glucose oxidation. The model helps to overcome the problem of the catalyst’s fast reversible deactivation by appropriate calculation of the target GLR.
In this chapter presented theoretical analysis of the target reactor optimal design procedure also offers several extensions to the future research. It should firstly focus on experimental validation of the presented GLR mathematical model and after that on the model extension for non-isothermal reactions. This would make the GLR model a useful tool for a process scale-up and its safety control.
Another direction of subsequent research efforts might be a critical comparison of chemical and biotechnological oxidation routes for syntheses of chemical specialties. An objective confrontation would be valuable from the viewpoint of technical arrangement of the processes as well as from the viewpoint of the two processes economics.
5. Nomenclature
Across-sectional area (m2)
0initial, at
Σ
Sepseparator
SSTRsemi-continuous stirred tant reactor
References
- 1.
Bello R. A. Robinson C. W. Moo-Young M. 1984 Liquid circulation and mixing characteristics of airlift contactors .62 573 577 . - 2.
Biella S. Prati L. Rossi M. 2002 Selective oxidation of D-glucose on gold catalyst .206 242 247 ;doi:10.1006/jcat.2001.3497. - 3.
Bisio A. Kabel R. L. 1985 Scaleup of chemical processes; Conversion from laboratory scale tests to successful commercial size design. John Wiley & Sons. - 4.
Blažej M. 2004 , Department of Chem. and Biochem. Engng., Slovak University of Technology in Bratislava. - 5.
Blažej M. Juraščík M. Annus J. Markoš J. 2004 Measurent of mass transfer coefficient in an airlift reactor with internal loop using coalescent and non-coalescent liquid media.79 1405 1411 ;doi: 10.1002/jctb.1144 . - 6.
Blažej M. Kiša M. Markoš J. 2004 Scale influence on the hydrodynamics of an internal loop airlift reactor .43 1519 1527 ;doi: 10.1016/j.cep.2004.02.003 . - 7.
Boelhouwer J. G. Piepers H. W. Drinkenburg A. A. H. 2002 Advantages of forced non-steady operated trickle-bed reactors .25 647 650 . - 8.
Butt J. B. Petersen E. E. 1988 Activation, deactivation and poisoning of catalysts. Academic Press, Inc. - 9.
Comotti M. Della Pina C. Falletta E. Rossi M. 2006 Aerobic oxidation of glucose with gold catalyst: Hydrogen peroxide as intermediate and reagent . 348(3),313 316 . - 10.
Eya H. Mishima K. Nagatani M. Iwai Y. Arai Y. 1994 Measurements and correlation of solubilities of oxygen in aqueous solutions containing glucose, sucrose and maltose.94 201 209 . - 11.
Godó Š. Klein J. Polakovič M. Báleš V. 1999 Periodical changes of input air flowrate- a possible way of improvement of oxygen transfer and liquid circulation in airlift bioreactors.54 4937 4943 . - 12.
Gogová Z. Čamaj V. Hronec M. Stanček F. 2002 , PatentWO2004047980 EP1569747 (SK appl. No.: 1676/2002) - 13.
Gogová Z. Hanika J. 2009a Reactivation of a palladium catalyst during glucose oxidation by molecular oxygen .635 520 526 ;doi: 10.2478/s11696-009-0053-3 . - 14.
Gogová Z. Hanika J. 2009b Dynamic modelling of glucose oxidation with palladium catalyst deactivation in multifunctional CSTR; Benefits of periodic operation .1501 223 230 ;doi: 10.1016/j.cej.2009.02.020 . - 15.
Gogová Z. Hanika J. 2009c Model aided design of three-phase gas-lift reactor for oxidation accompanied with catalyst reversible deactivation. Chem. Eng. Technol.3212 1929 1940 ;doi:10.1002/ceat.200900191 . - 16.
Juraščík M. Blažej M. Annus J. Markoš J. 2006 Experimental measurements of the volumetric mass transfer coefficient by the dynamic pressure-step method in an internal loop airlift reactors of different scale.125 81 87 ;doi: 10.1016/j.cej.2006.08.013 . - 17.
Ketteler G. Ogletree D. F. Bluhm H. Liu H. Hebenstreit E. L. D. Salmeron M. 2005 In Situ Spectroscopic Study of the Oxidation and Reduction of Pd(111).127 18269 18273 ;doi:10.1021/ja055754y. - 18.
Kluytmans J. H. J. van Wachem B. G. M. Kuster B. F. M. Schouten J. C. 2003 Design of an Industrial-Size Airlift Loop Redox Cycle (ALRC) Reactor for Catalytic Alcohol Oxidation and Catalyst Reactivation .42 4174 4185 ;doi: 10.1021/ie020916+ . - 19.
Kunz M. Recker C. 1995 A new continuous oxidation process for carbohydrates.13 11 15 . - 20.
Liu G. Zhang X. Wang L. Zhang S. Mi Z. 2008 Unsteady-state operation of trickle-bed reactor for dicyclopentadiene hydrogenation .63 4991 5002 . - 21.
Lundgren E. Kresse G. Klein C. Borg M. Andersen J. N. De Santis M. Gauthier Y. Konvicka C. Schmid M. Varga P. 2002 Two-Dimensional Oxide on Pd(111). 88(24);doi:10.1103/PhysRevLett.88.246103. - 22.
Markusse A. P. Kuster B. F. M. Schouten J. C. 2001 Platinum catalysed aqueous methyl-a-D-glucopyranoside oxidation in a multiphase redox-cycle reactor .66 191 197 . - 23.
Shah Y. T. Kelkar B. G. Godbole S. P. Deckwer W. D. 1982 Design parameters estimation for bubble column reactors. 28(3),353 379 . - 24.
Sikula I. 2008 , Department of Chem. and Biochem. Engng., Slovak University of Technology in Bratislava. - 25.
Sikula I. Juraščík M. Markoš J. 2006 Modeling of enzymatic reaction in an internal loop airlift bioreactor. 60(6), 446-453;doi: 10.2478/s11696-006-0081-1. - 26.
Sikula I. Juraščík M. Markoš J. 2007 Modeling of fermentation in an internal loop airlift bioreactor.62 5216 5221 ;doi: 10.1016/j.ces.2007.01.050. - 27.
Sikula I. Markoš J. 2008 Modeling of enzymatic reaction in an airlift reactor using an axial dispersion model . 62(1),10 17 ;doi:10.2478/s11696-007-0073-9. - 28.
Silveston P. L. 1998 Composition modulation of catalytic reactors . Gordon and Breach Science Publishers. - 29.
Silveston P. L. Hanika J. 2002 Challenges for the periodic operation of trickle-bed catalytic reactors .57 3373 3385 . - 30.
Silveston P. L. Hanika J. 2004 Periodic operation of three-phase catalytic reactors. 82 1105 1142 . - 31.
Simmons G. W. Wang Y. Marcos J. Klier K. 1991 Oxygen Adsorption on Pd(100) Surface: Phase Transformations and Surface Reconstruction.95 4522 4528 . - 32.
Szépe S. Levenspiel O. 1970 Catalyst deactivation. . Oxford. - 33.
Thielecke N. Vorlop K. D. Prusse U. 2007 Long-term stability of an Au/Al2O3 catalyst prepared by incipient wetness in continuous-flow glucose oxidation .122 266 269 ;doi:10.1016/j.cattod.2007.02.008. - 34.
Tukač V. Hanika J. Chyba V. 2003 Periodic state of wet oxidation in trickle-bed reactor. 79-80, 427-431. - 35.
Vleeming J. H. Kuster B. F. M. Marin G. B. 1997 Selective Oxidation of Methyl a-D-Glucopyranoside with Oxygen over Supported Platinum: Kinetic Modeling in the Presence of Deactivation by Overoxidation of the Catalyst .36 3541 3553 . - 36.
Znad H. Báleš V. Markoš J. Kawase Y. 2004 Modeling and simulation of airlift bioreactors .21 73 81 ;doi: 10.1016/j.bej.2004.05.005 .