Classroom Analysis of a Packed Bed Reactor with and without the Use of Effectiveness Factor

Classroom analysis of a packed bed reactor with and without the use of effectiveness factor (η) is discussed in this article. The use of (η) in the calculation of an isothermal packed bed reactor with convection and first-order reaction rate is analyzed. An analytical model has been developed with and without (η). The model without (η) was presented in a 50-min lecture and later used in a 50-min computational laboratory session with various groups of two students to calculate the length and catalyst mass in a packed bed reactor. In another class, the students calculated the model using (η). Four or five class activities were tested every month in pairs of students in the place of traditional exams. A questionnaire was prepared to verify the students’ level of difficulty in understanding two theoretical topics taught in the last two classes, difficulty in doing the requested calculations, as well as among calculations using the effectiveness factor and without it. We have included a question in the questionnaire on how the students felt about the elimination of one traditional monthly exam in Chemical Reactor Design. In the two semesters when the questionnaire was applied to the Chemical Reactor Design class, on average, 86% of the students preferred various class group activities as opposed to traditional exams, 80% of the students reported that they had little difficulty in understanding the two theoretical sections in the last two classes, 57% stated that they did not have any difficulty in doing the requested calculations and 51% reported that the calculation methods with the use of (η) and without its use were somewhat similar.

neering programs and addresses many concepts that are applicable in various areas (Mahecha-Botero, Reaume, Grace, & Ellis, 2011).Computer-aided design has become extremely popular, and its use in classes may be very helpful by adding to the analysis capacity in all engineering areas (Cartaxo, Silvino, & Fernandes, 2014).
Academic laboratories have traditionally been used to complement and reinforce theoretical instruction that is conveyed through lectures in a practical way.However, data processing and model evaluation tasks are time consuming and do not add much to students' learning experiences because they reduce the time that is available for result analysis, critical thinking and report writing skill development (Stammitti, 2013).Stammitti (2013) selected three Transport Phenomena Laboratory experiments, namely metallic bar temperature profiles, transient heat conduction and fixed and fluidized bed performance, and developed a spreadsheet for each one.These spreadsheets comprise analytical and numerical solutions of different models, as well as correlations available in the literature.Fifteen chemical engineering students who tested the spreadsheets were surveyed, showing that spreadsheets were considered useful for reducing workload and boosting the quality of the analyses because students had the new possibility of quickly practicing with diverse correlations and models (Stammitti, 2013).
According to Glassey, Novakovic, & Parr (2013), traditional curriculum delivery in higher education has long been considered ineffective in promoting deep learning.Enquiry Based Learning (EBL) provides an opportunity to develop important professional attributes within subject-specific content.Computer-aided learning packages can be a useful tool in supporting the development of these skills because they enable students to explore and gain experience in new software environments in a subject-specific context.Glassey, Novakovic, & Parr (2013) developed and delivered case studies in Stages 1 and 2 of a chemical engineering degree program by utilizing a number of different software packages.The Stage 2 case studies concentrate on separation processes and the reactor engineering aspect while requiring application of the knowledge of statistics and design of experiments, respectively.A turning point questionnaire was administered to Stage 2 students (n = 66) following the completion of the case study containing seven questions about the effectiveness of the case study and students' perceptions of its impact.Anonymous comments and additional suggestions were also collected as free-text comments (Glassey, Novakovic, & Parr, 2013).
According to Regalado-Méndez, Cid-Rodríguez, & Báez-González (2010), their work is focused on a project that integrates curriculum such as thermodynamics, chemical reactor engineering, linear algebra, differential equations and computer programming.
The purpose is for students to implement most of their knowledge and tools to analyze a stirred tank chemical reactor as a simple dynamic system.The students showed a bigger interest in this practice because they worked in groups.The most important fact is that the percentage of failure among students was 10% (Regalado-Méndez, Cid-Rodríguez, & Báez-González, 2010).Dahm & Hesketh (2008) describe two chemical reaction experiments developed for Rowan University's introductory course on chemical reaction engineering, i.e., an este-rification reaction carried out in a packed bed and a competitive reaction in which the kinetics were influenced by micromixing.Dahm & Hesketh's (2008) paper describes the use of the experiments in the classroom and compares the performance of the 2005 students to the 2004 cohort for whom the course included no wet labs at all.According to Luzi, Bressa, Mazza, & Barreto (2014), many undergraduate chemical engineering curricula do not include non-catalytic or homogeneously catalyzed Gas-Liquid Reactions (GLRs) and heterogeneous catalysis.Available textbooks dealing with GLRs describe and separately formulate the different processes affecting the rate of chemical absorption.From the analogy with the effectiveness factor concept (η), a liquid utilization factor is defined in gas-liquid reactions (Froment, Bischoff, & Wilde, 2011: p. 329).The Hatta number is very similar to the Thiele modulus used in the (η) (Froment, Bischoff, & Wilde, 2011: p. 328) approach of the present article.This dimensionless group provides a convenient measure of the promoting effect of a chemical reaction on the rate of absorption (Bird, Stewart, & Lightfoot, 2002: p. 696).There are several similarities between the gas-liquid reaction modeling presented by Luzi, Bressa, Mazza, & Barreto (2014) and the subject covered in this article.
In the course on chemical reactor design in the Chemical Engineering Department of the Federal University of the Rio Grande do Norte, Brazil, the laboratory modules are carried out concomitantly with or in parallel to theoretical classes.Subsequent to the lectures, topics such as Chemical Reactor Design basic equations are developed on Excel spreadsheets and Fortran programs for the numerical simulation of reactors.The overall quality of the resulting student work is very good, and the students generally do not struggle to discuss their experimental investigation results (Ganley, 2015).
Three years ago we began to teach the Chemical Reactor Design discipline, we were basing the teaching of the Chemical Reactors design discipline on our 15-year experience in the teaching of Transport Phenomena.The traditional monthly exams were replaced with class activities.The activities were performed in pairs using Excel spreadsheets and were developed and employed in Fortran simulation programs.
The students were able to develop spreadsheets for calculations in Sections 2 and 3 of the present article.Details on the important concepts are considered as a prerequisite at this point and are taught in the mass transfer course at our university.Nevertheless, the chemical reactor design lecturer discussed the effectiveness factor (η) concept to facilitate the students' task.
In 1939, Thiele introduced the concept of an internal effectiveness factor (η internal ) in the well-known paper titled, "Relation between catalytic activity and size of particle" (Thiele, 1939).The (η internal ) was defined as being the rate of the reaction with pore diffusion resistance divided by the rate of the reaction on surface conditions R A(surf) [kg•mol/(m 3 •s)] in all catalyst particle active sites (Froment, Bischoff, & Wilde, 2011: p. 194).The observed reaction rate R A(obs) in a particle can be expressed by R A(surf) multiplied by (η internal ).
The overall effectiveness factor (η overall ) is defined as being the actual R A(obs) divided by the rate R A(reference) that would result if the entire surface of all the active sites of the catalyst particle was exposed to the bulk concentration.The R A(obs) in a particle can be ex-pressed by R A(ref) multiplied by (η overall ).Various authors (Bird, Stewart, & Lightfoot, 2002;Fogler, 1992Fogler, , 2006;;Froment, Bischoff, & Wilde, 2011;Thiele, 1939;Yamané, Araki, & Sada, 1981) have presented similar definitions of (η internal ) and (η overall ).
In this research, the use of (η) in the calculation of an isothermal packed bed reactor with convection and reaction is discussed, and it is shown that the bed length [L (m)] and catalyst mass [M p (kg)] in a packed bed can be calculated with or without (η).A questionnaire was prepared to verify the students' level of difficulty in understanding two theoretical topics taught in the last two classes, in doing the requested calculations, as well as the similarity they found between using the calculation methods with the use of the effectiveness factor and without it.We also investigated how the students felt about the elimination of traditional monthly exams in Chemical Reactor Design.

Packed Bed Reactor without Effectiveness Factor
The reaction in a packed bed is assumed to take place in many particles, rather than in a single particle.C Af (kg-mol/m 3 ) is the bulk fluid-phase concentration of (A) at any point along the length of the packed bed (Fogler, 1992).
Various spherical particles are present in a catalytic reactor and are submerged in a fluid phase containing the mixture of reactant (A) and product (B).We assume that each catalyst particle is surrounded by a fluid in movement, through which component (A) has to move to reach the surface of the catalyst.On the catalyst's surface, the reaction A → B is thought to occur on the catalyst surface and the reaction product then diffuses back out through the mixture composed of reactants and products (Bird, Stewart, & Lightfoot, 2002: pp. 551, 565).
Any radial variations in concentration are neglected and the packed bed is assumed to operate at a steady state.The mass balance for component (A) based on bulk concentration (C Af ) is solved for the case in which the flow rate through the packed bed is very large, i.e., a plug flow model is adopted and the axial dispersion can be neglected (Fogler, 1992: p. 630;Froment, Bischoff, & Wilde, 2011).A mass balance on (A) of a slice of thickness (dz), cross section ( 2 π b R ) of the packed bed of volume , where (R b ) is the packed bed radius and v sup (m/s) is the fluid superficial velocity, using a one-dimensional model and an isothermal packed bed reactor Volume of all particles in cross section m d Assuming (V p ) as the volume of all particles in the considered cross section, the fluid volume (V f ) in the same section, the total volume total and the porosity of the bed Volume of all particles in cross section m 1 Volume of considered cross section m The flux N Ap [kg•mol/(m 2 •s)] of the component (A) in a particle with external mass transfer coefficient k m [m 3 /(m 2 •s)] and the concentration of (A) on the external catalyst surface (C ARp ) for particle surface r = R p is Multiplying the surface area (A p ) of the sphere, and by dividing its volume (V p ) by the flux (N Ap ), the observed reaction rate in a spherical particle is obtained Equation ( 5) and Equation ( 6) are valid for one catalyst particle.They were adapted for a packed bed reactor by multiplying the volume fraction occupied by all catalyst particles in the bed as a function of the bed porosity ε b (dimensionless).Specifically, multiplying (1 b ε − ) by the observed reaction rate in a particle R A(obs) or using the parameter the observed reaction rate in the bed R A(obs. in bed) with volume Internal diffusive resistance in the particle 1 External mass transfer resistance in the fluid film where F 1 (dimensionless) is a parameter, φ 1 (dimensionless) is the Thiele modulus, 1 k′′ [m 3 /(m 2 •s)] is the kinetic coefficient of the superficial reaction rate, ρ p (kg/m 3 ) is the particle density, S g (m 2 /kg) is the available catalytic surface per unit mass, and Bi m (dimensionless) is the Biot number for mass transfer.
Integrating Equation ( 8) from z = 0 and (C Af0 ) to z = z (m) and (C Af ) gives ( ) The length of the packed bed (L) is obtained by reorganizing Equation ( 12) The mass of catalyst particles M p (kg) in a packed bed of volume

Example of a Packed Bed Reactor without the Effectiveness Factor
This section presents the solution performed by the students for an example of using the equations for a packed bed with reaction and convection but without effectiveness factor.
To calculate the length of the packed bed reactor (L) and weight of the porous solid in bed (M p ), it is necessary to reduce the NO concentration to a level of 0.004%, which is below the Environmental Protection Agency limit (Fogler, 2006).

Packed Bed Reactor with Effectiveness Factor
In Section 3.1, a model for a catalytic particle in which a first-order reaction is consi-dered.In Section 3.2, an equation valid for one catalyst particle is adapted for a packed bed reactor by multiplying the volume fraction occupied by all catalyst particles in the packed bed, and in Section 3.3 an example of calculation using the effectiveness factor is presented.

Diffusion and Reaction in a Single Spherical Catalytic Particle
Let us now consider a model for a catalytic particle in which a first-order reaction is carried out.It is admitted that a fluid in movement surrounds each catalyst particle and component (A) has to move through the catalyst particles to reach the catalyst surface (Bird, Stewart, & Lightfoot, 2002: p. 565).
The mass balance to describe the diffusion and chemical reaction R A [kg⋅mol/(m 3 ⋅s)] within a single porous spherical particle with radius (R p ), (C Ap ) concentration in the catalyst and effective diffusivity of the component (A) in the particle D ef (m 2 /s) is The situation where species (A) becomes, according to a first-order chemical reaction , is considered and Equation ( 22) becomes Reordering this equation results in Naturally, the internal diffusive resistance, ( ) This boundary condition is reordered ( ) It is obvious that the internal diffusive resistance ( ) ( ) includes the effect of the internal diffusional resistance, external mass transfer resistance and reaction rate resistance in active sites.
If r = R p in Equation ( 28), the concentration of (A) on the external catalyst surface (C ARp ) with parameter (F 1 ) in Equation ( 9) is ( ) Substituting Equation ( 29) into (26) gives the flux on surface (N Ap ) The observed reaction rate in a catalyst particle R A(obs) is obtained from Equation ( 6).
Replacing Equation ( 29) into ( 6) When the concentrations everywhere are C Ap = C Af in all of the active sites inside the particle, the reference reaction rate in the catalyst is Taking the ratio of the two last equations and using the (η overall ) definition Similarly to (C Ap ), the function  28), dC Ap /dr gives the flux using Equation ( 26) or Equation ( 27) and N Ap = −D ef dC Ap /dr at r = R p , which, multiplied by the outside area (A p ) per volume (V p ) of spherical particle, results in ( ) ( ) The result of Equation ( 34) is equal to Equation (31).Substituting Equation ( 34) into (33), using (φ 1 ) and (F 1 ), results in a general solution Reordering Equation ( 35) gives the analytical equation of (η overall ) for a first-order reaction and incorporates the influences of the resistances into the process in a spherical geometry

Convection and Reaction in a Packed Bed Reactor with Effectiveness Factor
Equation ( 4) is a general solution and is valid for one catalyst particle.Equation ( 4) can be adapted for a packed bed reactor by multiplying the volume fraction occupied by all catalyst particles in a packed bed, i.e., (1 b ε − ) by the actual observed reaction rate in a particle R A(obs) .Substituting this parameter in Equation ( 1), using a one-dimensional model, an isothermal packed bed reactor and assuming that the packed bed operates at a steady state ( ) ( ) Equation ( 37) is similar to Equation (11-69) from Fogler (1992: p. 630) and Equation (12-69) from Fogler (2006: p. 844) Integrating Equation ( 37) from z = 0 and C Af = C Af0 to z = L overall and C Af give a concentration of (C Af ) and the length of the packed bed based on (η overall ) and L overall ( ) ( ) Equation ( 39) is similar to Equation (E11-3.2) in Fogler (1992: p. 630) and Equation (12-70) in Fogler (2006: p. 844).
The mass of the catalyst particles in the packed bed (M p ) in volume

Example of a Packed Bed Reactor with Effectiveness Factor
The parameters were calculated using Equations ( 35), ( 39) and ( 40), with φ 1 = 18.01 and F 1 = 0.368, based on data from Fogler (2006) ( ) ( )( ) The calculated values for the reactor length and mass of the catalyst are the same for both methods.The results presented in this work are similar to the values of the calculations performed by Fogler (2006) with a small difference resulting from simplification of the equation of (η overall ).
The "mathematical trick" used in the present paper was the substitution of the concentration of component (A) on surface (C ARp ) with C ARp = F 1 C Af from Equation (29).
For example, in the flux in Equation ( 26) or Equation ( 27) the calculation of (L) and (M p ) can be completed without the necessity of using (η overall ).

Results and Discussion
The model presented without the use of (η overall ) was employed by reactor engineering students in two classes: one for a theoretical lecture by the professor and another for calculation by the students.The model without (η overall ) was presented in a 50-min class, and in another 50-min class, various groups of two students calculated the length of the reactor and the mass of the catalyst in the bed, which is presented in Section 2. In another class and in a similar way, the students calculated the model with (η overall ), which is presented in Section 3.However, we are in the whole reactor design discipline replace traditional exams for activities in pairs, fully realized in the classroom.
We have used a methodology similar to that employed by Cartaxo, Silvino, & Fernandes (2014) in the application of a questionnaire for students because privacy is a common concern for students when responding to surveys regarding pedagogical methods because they usually assume that their answers may affect their grades negatively.
Thus, to promote student engagement, no identification data were requested from the respondents and appropriate care was taken during the questionnaire application (Cartaxo, Silvino, & Fernandes, 2014).reported having little difficulty in understanding the two theoretical sections in the last two classes and 29% reported having no difficulty at all.Fifty-seven percent said that they did not have any difficulty in doing the calculations requested.Forty-three percent reported that the calculation methods with the use of the effectiveness factor and without its use were similar, and 51% found the methods somewhat similar.Stammitti (2013) concludes that the goal of updating the teaching and learning experience of chemical engineering students in the laboratory through computers was fulfilled by the introduction of spreadsheets.In the present article, our students created their own Excel spreadsheets for calculations with and without the effectiveness factor that were simpler than Stammitti's (2013).Writing Fortran programs was not necessary because the activities were simple.Only one student described another type of evaluation, which as follows: The in-class activities are quite productive; however, I suggest using a more detailed chemical reactor design for evaluation.

Exam and activity.
When necessary, we develop our own computer programs instead of using readily or commercially available simulators.One of the reasons is that the students will most likely retain their knowledge of Fortran programming even after some simulators are out of the market.Sometimes we compare our results to those of commercial simulators available in books and articles.Question correct answers; when they used the effectiveness factor, all the answers were correct.
Student feedback regarding the homogeneous kinetics laboratory module has been generally very positive (Ganley, 2015).Student attitudes and responses toward the module and its delivery indicated greater student satisfaction because the experimental hardware and the instructional documentation were each successively refined over the three semesters that the experiment has been in use.
The students used the two methodologies developed in class and considered the results that were obtained to be natural and positive for their training in chemical engineering.Student feedback regarding the calculations and studies performed in the computer laboratory were generally very positive.

Conclusions
This article describes two methods of calculation of a packed bed reactor involving reaction diffusion and convection concepts.Classroom analysis of a packed bed reactor with and without the use of the effectiveness factor (η) is discussed.The use of (η) in the calculation of an isothermal packed bed reactor with convection and first-order reaction rate is analyzed.
An analytical model has been developed with and without (η).The model without (η) was presented in a 50-min lecture and later used in a 50-min computational laboratory session with various pairs of students to calculate the length and catalyst mass in a packed bed reactor.In another class, the students calculated the model using (η).
A questionnaire was prepared to verify the students' levels of difficulty in understanding the two theoretical topics taught in the last two classes, in doing the requested calculations, as well as the similarity they found between the calculation methods with the use of the effectiveness factor and without it.However, when consulted about their satisfaction with the elimination of traditional exams in Chemical Reactor Design, in an average of two semesters, 86% of the students preferred various class group activities in comparison to traditional exams, 80% reported that they had little difficulty in understanding the two theoretical sections in the last two classes, 57% stated that they did not have any difficulty in doing the requested calculations and 51% reported that the calculation methods with the use of (η) and without its use were somewhat similar.
The calculation results obtained with or without (η) were identical; i.e., (η) can be used but it is not mandatory.The analysis presented in this research provides motiva-

(
C ARp ) is obtained from the solution of the mass balance of component (A) in an individual catalytic spherical particle with diffusion and reaction.(C ARp ) is obtained and discussed a posteriori in Equation (29) in Section 3 p. 845) (L) and (M p ) are obtained

1 1
k′′ , are included in the Thiele modulus (φ 1 ) in Equation (24).This equation is solved with a finite concentration of C Ap at r = 0 and a flux outside the surface of the particle at r = R p k m are included in the Biot number for mass transfer (Bi m ).The solution of Equation (24) using Equation (25) and Equation (26) is(Bird, Stewart, &   Lightfoot, 2002: p. 853, Equation (C.1.6a)) internal diffusional resistance, the external mass transfer resistance and the reaction rate resistance.(η overall ) can be written in another form based on the diffusive flux (−D ef dC Ap /dr) on the external surface at r = R p .The derivative of C Ap in Equation ( 0.5 4.42 10 m s 2.8 10 g m 530 m g 2.8 10 g m π 0.0254 m 0.162 m 456.8 g (7) presented by Regalado-Méndez, Cid-Rodríguez, & Báez-González (2010) is related to the present article and is presented next.I understand all the concepts, with possible responses Yes and No.According to Regalado-Méndez, Cid-Rodríguez, & Báez-González (2010), 80% of the students understood the concepts.In the present article, 80% of the students reported having little difficulty in understanding the two theoretical sections in the last two classes.Dahm & Hesketh (2008) used specific student learning performance indicators from 2004 and 2005 such as indicator (1); the problem solution is usually correct and always reasonable, with values of 4 = best and 1 = worst.In the results presented by Dahm & Hesketh (2008) from the spring of 2004, the students obtained a mean of 3.0, while in the spring of 2005, 12 students obtained a mean of 4.0.The 2004 class had some instances of students giving physically unrealistic answers such as temperature increases from endothermic reactions.Such errors were non-existent in 2005.In the present paper, 35 students performed the activity without the effectiveness factor and 34 obtained

Table 1
presents the questionnaire answered by the students.Out of the 23 students enrolled in the discipline in semester 2015.1, 19 were present and performed the calculations and 15 answered the questionnaire.Out of the 41 students enrolled in the discipline in semester 2016.1, 22 were present and performed the calculations and 20 answered the questionnaire.In an average of two semesters, 86% of the students reported feeling motivated by the varied class activities in comparison to traditional exams, 80%

Table 1 .
Questionnaire on the two last classes and the activities performed in the discipline.