Profile of the Effectiveness Factor under Optimal Operating Conditions for the Conversion of Ortho-Xylene to Phthalic Anhydride in a Fixed-Bed Tubular Reactor

: The objective of this research is to find the effectiveness factor of the catalyst particles for the most favorable conditions of the phthalic anhydride production in a fixed bed reactor, with the aim of achieving the highest rate of phthalic anhydride production compared to other secondary products and analyzing the areas of lower effectiveness for the modification of the reactor design. Initially, the material and the energy balances in the catalytic bed are solved to obtain the concentration and temperature profiles based on the radius and length of the reactor, using polymath software(Polymath ® v6.2 Software Minitab 19 Matlab 2019) with the data from literature. Once the profiles reproducibility was verified using the initial data (inlet temperature, pressure in the reactor, reactor wall temperature, reactor radius and mass flow rate) the experimental design 3 5 carry out, which generates 243 “experiments”, whose response variable (phthalic anhydride concentration) was obtained using Matlab. Subsequently, the variables were analyzed using the Minitab 18 ® that, through the response surface analysis method, allowed us to obtain the optimal values of the tested variables. Then, Subsequently, material and energy balances coupled with Fourier and Fick’s laws, along with the effectiveness factor equation, were applied, resulting in the generation of 9 coupled differential equations. Upon implementing the finite difference method, this yielded 90 nonlinear algebraic equations, which were solved using the Polymath software. A total of 78 particles were preselected based on their radial and axial positions to determine the effectiveness factor profile, with values ranging from 0.83 to near unity. The lower values correspond to the points with higher temperature, as evidenced by the calculations performed.


Introduction
Packed beds of catalyst particles are the most widely used for gas phase reagents in chemical production.Studies about fixed-bed tubular reactors can be found referring to different internal designs, different applications, microreactors, among others.Researchers have been studying the phenomenon that occurs in them, since the last century [1][2][3].Fixed bed tubular reactors are used preferentially for exothermic catalytic processes [4], mainly bidimensional models with radial diffusion.
The reaction of the organic compound, phthalic anhydride, is carried out in a fixed bed reactor.This product is of industrial interest, as it serves as a raw material in the production of plastics, solvents, resins, cellulose, agricultural fungicides, plastics, amines, among others [5].One way to obtain this product is by conversion of o-xylene.In the literature, there are different conversion routes for this process.Historically, the raw material for phthalic anhydride production was the oxidation of naphthalene, which has gradually been replaced by o-xylene over the past 50 years.Considering the quantities produced on an industrial scale [6,7], even a 1% increase in selectivity has a substantial economic impact [8].An important limitation in optimization processes is having experimental reactor data, operating conditions, and physical property information; considering this aspect, as most data related to the partial oxidation of ortho-xylene is available, we proceeded with its optimization, supplementing missing data with typical values of physical properties.Ivanovskaya and Sembayev [9] comment that there are discrepant conversion routes of o-xylene to phthalic anhydride.However, the reaction is still the focus of different research, both theoretical and experimental as shown below.
Chandrasekharan and Calderbank [10] studied the production of phthalic anhydride from the oxidation reaction of o-xylene in a tubular reactor.The results showed that the incorporation of TiO 2 into V 2 O 5 , alters the kinetic parameters so that the phthalic anhydride precursors are less evident than in the case of V 2 O 5 alone.Gimeno et al. [11] studied the kinetics of the partial oxidation process of ortho-xylene in a fluidized bed reactor using V 2 O 5 /SiO 2 catalyst, indicating that the reactor allowed obtaining reliable kinetic data despite the strong exothermic nature of the reaction.A tool to predict or estimate parameters is the mathematical model, which with well-founded hypotheses leads to good results, which can be validated experimentally.
Nikolov and Anastasov [12] researched by means of a mathematical model the effect of the inlet temperature both on the temperature regime in the reactor and the selectivity with respect to the desired product, managing to verify the results experimentally.These results were used to validate the vapor phase oxidation process of o-xylene in phthalic anhydride with air.
The experimental study performed by Papageorgiou et al. [13] in an integral isothermal fixed-bed reactor exposed that the selectivity's of products some of the catalytic oxidations of o-xylene were strongly influenced by the partial pressure of oxygen.The product is o-tolualdehyde, phthalide and phthalic anhydride.Also, the authors developed a kinetic model for analyzed the dependence of the selectivity's on the composition of the reaction mixture and the results in both are in excellent agreement.Andrigo et al. [3] carried out an analyze of the phenomena that occur in fixed bed reactors from the reaction with a single pellet, where reaction and diffusion competing, to reactions where dispersion and heat transfer play a role important.
The deactivation in the oxidation process of o-xylene into phthalic anhydride was studied in a fixed bed industrial reactor by Anastasov [14].The study was carried out experimentally to obtain the temperature profiles along the length of the reactor.And a mathematical model was used to determine the degree of catalyst deactivation along the bed as a function of time.
Castillo-Araiza and López-Isunza [15] analyzed the performance of hydrodynamics in the packed-bed catalytic reactor.In addition, they used the partial oxidation of o-xylene to study the behavior of redox dynamics of the catalyst surface together with the use of a catalyst activity profile.The catalyst used was vanadium oxide.The same authors presented some perceptions of the fundamental role of catalytic and deactivation kinetics to partial oxidation o-xylene on a commercial V 2 O 5 /TiO 2 catalyst to both in a steady state and in a transient state in an industrial-scale reactor with packed bed single-tube [16].Sethapokin et al. [17] developed a two-dimensional heterogeneous mathematical model to predict the molar flux rate of the by-subproduct of phthalic anhydride and maleic anhydride in an industrial fixed-bed reactor.They estimated the kinetic parameters of the catalytic bed from a design of experiments.The simulated results were validated with production data from a case study using V 2 O 5 /TiO 2 catalyst achieving an accuracy of 1.8% deviation from the production data.
Also, some researchers used computational fluid dynamics (CFD) and ANSYS fluent were used to analyze the conversion behavior of o-xylene to phthalic anhydride or other reactions in a fixed-bed catalytic reactor [18][19][20].Thus, the aim of this study is to determine the effectiveness factor profile of catalyst particles in a fixed-bed tubular reactor, in order to identify the regions of highest catalyst activity based on the radial and axial position within the reactor.To achieve this, operational conditions within the reactor are optimized, using ortho-xylene concentration as the dependent variable and inlet temperature, tube wall temperature, pressure within the reactor bed, reactor diameter, and mass flow rate as independent variables.Carrasco et al. developed a method for calculating the effectiveness factor profile applied to cyclohexanol dehydrogenation, where the model deduction and calculation procedure are illustrated using numerical techniques and statistical tools [21].This procedure serves to identify regions of highest catalytic activity, which are primarily dependent on temperature, enabling appropriate precautions to be taken before conducting the respective experiment.Additionally, it allows for the indirect determination of the minimum required reactor length for the given conditions.The mathematical modeling considers the oxidation of o-xylene over a vanadium oxide catalytic medium, utilizing the triangular scheme presented in Figure 1.
ChemEngineering 2024, 8, x FOR PEER REVIEW 3 of 24 presented some perceptions of the fundamental role of catalytic and deactivation kinetics to partial oxidation o-xylene on a commercial V2O5/TiO2 catalyst to both in a steady state and in a transient state in an industrial-scale reactor with packed bed single-tube [16].Sethapokin et al. [17] developed a two-dimensional heterogeneous mathematical model to predict the molar flux rate of the by-subproduct of phthalic anhydride and maleic anhydride in an industrial fixed-bed reactor.They estimated the kinetic parameters of the catalytic bed from a design of experiments.The simulated results were validated with production data from a case study using V2O5/TiO2 catalyst achieving an accuracy of 1.8% deviation from the production data.
Also, some researchers used computational fluid dynamics (CFD) and ANSYS fluent were used to analyze the conversion behavior of o-xylene to phthalic anhydride or other reactions in a fixed-bed catalytic reactor [18][19][20].Thus, the aim of this study is to determine the effectiveness factor profile of catalyst particles in a fixed-bed tubular reactor, in order to identify the regions of highest catalyst activity based on the radial and axial position within the reactor.To achieve this, operational conditions within the reactor are optimized, using ortho-xylene concentration as the dependent variable and inlet temperature, tube wall temperature, pressure within the reactor bed, reactor diameter, and mass flow rate as independent variables.Carrasco et al. developed a method for calculating the effectiveness factor profile applied to cyclohexanol dehydrogenation, where the model deduction and calculation procedure are illustrated using numerical techniques and statistical tools [21].This procedure serves to identify regions of highest catalytic activity, which are primarily dependent on temperature, enabling appropriate precautions to be taken before conducting the respective experiment.Additionally, it allows for the indirect determination of the minimum required reactor length for the given conditions.The mathematical modeling considers the oxidation of o-xylene over a vanadium oxide catalytic medium, utilizing the triangular scheme presented in Figure 1.This scheme was modified from the literature, that it is found in generalized form diagrams or reaction network for o-xylene oxidation [22][23][24][25][26].The corresponding reaction mechanism is: This scheme was modified from the literature, that it is found in generalized form diagrams or reaction network for o-xylene oxidation [22][23][24][25][26].The corresponding reaction mechanism is:

Reactor Model
In this work, the modeling is based on the scheme presented in Figure 2, where material and energy balances are considered in the reactor bed with radial diffusion.From this, radial and axial profiles of the concentration of each component, denoted as Cj(r, z), and temperature, denoted as T(r, z), are obtained.Subsequently, material balances for each component are performed within the spherical particle of radius r, denoted as Cj(r), and temperature T(r).This must be coupled with Fick's and Fourier's laws, resulting in concentration profiles for each component, molar flux density of each component, temperature profile, heat flux density profile, and the effectiveness factor of each particle along the bed using the implicit finite difference method.and temperature, denoted as T(r, z), are obtained.Subsequently, material balances fo each component are performed within the spherical particle of radius r, denoted as Cj(r and temperature T(r).This must be coupled with Fick's and Fourier's laws, resulting i concentration profiles for each component, molar flux density of each component, tem perature profile, heat flux density profile, and the effectiveness factor of each particl along the bed using the implicit finite difference method.It is known that the general equation of mass transport is given by: whose modified form for the equimolar counter-diffusion heterogeneous system is pre sented as follows: It is known that the general equation of mass transport is given by: whose modified form for the equimolar counter-diffusion heterogeneous system is presented as follows: the energy balance general form is given by: a suitable simplification of Equation (3) can be written as: when the reaction rate is given by: The heat generation rate, by chemical reaction, is given by: the necessary boundary conditions are: When z = 0, then When r = 0, then ∂C j ∂r When r = R, then From Equation (2), is obtained: The mass transfer equation can be written in terms of the Peclet number for mass transfer given by: with which, Equation (12) becomes: The following variable changes are made: and substituting (15) and ( 16) in (14): Equation ( 17) is used for all points on the radial position except the center.When r = 0, Equation ( 17) is indeterminate, therefore, the L'Hôpital's rule is applied, resulting in: this equation is valid only in the tube center.Discretizing the Equations ( 17) and (18) we get: From Equation (4) it is obtained: this equation can be written in terms of the Peclet number for heat transfer shown in: by substituting (22) in (21) we get: Following equations indicate variable changes: substituting both equations in (23), is obtained: L'Hôpital's rule is applied to solve the previous equation when r = 0, resultant in: this equation is valid only in the tube center.Discretizing the Equations ( 26) and ( 27) we get: Additionally, it is required the use of the Biot number in the wall as boundary conditions:

Reactor Simulation
Garcia-Ochoa et al.
[27] used the orthogonal collocation method to solve the material and energy balances.In the present work, the partial differential equations system was solved using the discretization with method of line (MOL) by polymath (v6.2), this with the purpose of validating the model and to identify the factors that influence the optimal conversion of o-xylene to phthalic anhydride.The data used to simulate the process on the catalytic bed and in the catalyst particle are presented in Tables 1 and 2. properties Adapted from [27,28].
For this work, the following hypotheses was considered: Effective diffusivity of the mixture: this data fluctuates between 10 −8 and 10 −10 m 2 /s.In this case, the value of 9 × 10 −9 m 2 /s has been proposed.
The concentration profiles of o-xylene, phthalic anhydride, and by-products oxidation produced in the catalytic bed for the ten placement points are shown in Figure 3.The phthalic anhydride concentration increases as the reaction progresses through the reactor.
For this work, the following hypotheses was considered: Effective diffusivity of the mixture: this data fluctuates between 10 −8 and 10 −10 m 2 /s.In this case, the value of 9 × 10 −9 m 2 /s has been proposed.
The concentration profiles of o-xylene, phthalic anhydride, and by-products oxidation produced in the catalytic bed for the ten placement points are shown in Figure 3.The phthalic anhydride concentration increases as the reaction progresses through the reactor.The temperature profiles produced in the catalytic bed as a function of the reactor length are presented in Figure 4.Note that the reaction temperature increases rapidly near the reactor inlet reaching the maximum temperature.This behavior is the same as that reported in literature [27] for the same conditions of T0 = Tw = 625.15K.The validation of the fixed bed reactor model was obtained with the reproduction of the profiles (Figures 3 and 4) of Garcia-Ochoa [27].Note that, after 3 m of length, the temperature variation is minimal.The factors selected for optimization by statistical treatment were the feed temperature, the reactor radius, the temperature of the reactor wall, the total pressure, and the mass flow rate and as a response it is proposed to obtain the highest concentration of phthalic anhydride.

Factorial Design
The conversion process of o-xylene in phthalic anhydride was analyzed in a fixedbed tubular reactor in order to find the best concentration at the reactor outlet.The research was carried out using the design of experiments (DOE), which allows studying the relationship between the factors (input variables) and the responses (output variables).
The Matlab ® software (R2015a) was used for the data generation through the resolution of the differential equations system, and Minitab was used for optimization.The results of interest were the concentration profiles.In this sense, for an adequate management of the information, the indicator "average value of phthalic anhydride concentration at the reactor outlet" has been defined as a response variable.The flowchart shown in Figure The validation of the fixed bed reactor model was obtained with the reproduction of the profiles (Figures 3 and 4) of Garcia-Ochoa [27].Note that, after 3 m of length, the temperature variation is minimal.The factors selected for optimization by statistical treatment were the feed temperature, the reactor radius, the temperature of the reactor wall, the total pressure, and the mass flow rate and as a response it is proposed to obtain the highest concentration of phthalic anhydride.

Factorial Design
The conversion process of o-xylene in phthalic anhydride was analyzed in a fixed-bed tubular reactor in order to find the best concentration at the reactor outlet.The research was carried out using the design of experiments (DOE), which allows studying the relationship between the factors (input variables) and the responses (output variables).
The Matlab ® software (R2015a) was used for the data generation through the resolution of the differential equations system, and Minitab was used for optimization.The results of interest were the concentration profiles.In this sense, for an adequate management of the information, the indicator "average value of phthalic anhydride concentration at the reactor outlet" has been defined as a response variable.The flowchart shown in Figure 5 summarizes the process.The analyses for the experimental design use data involving five input variables and three levels, which are shown in Table 3.For the experimentation development, the factorial design table has been generated using the Minitab 18 ® software.The number of experimental runs corresponds to 3 5 = 243, according to Table 4.

Source
Degrees of Free-F-Value p-Value The analyses for the experimental design use data involving five input variables and three levels, which are shown in Table 3.For the experimentation development, the factorial design table has been generated using the Minitab 18 ® software.The number of experimental runs corresponds to 3 5 = 243, according to Table 4.  Thus, Figure 6 shows the graph of main effects for the concentration of phthalic anhydride in the reactor outlet.The three points in each graph represent low, medium, and high values for each factor; the medium values were taken from Table 1 (inlet temperature, external reactor wall temperature, operating pressure, reactor radius, and mass flow rate of ortho-xylene), which likely have the greatest influence on the dependent variable (CB); while the high and low values were assigned with reference to the central point.And it is observed that for the condition T w = 625.15,radius = 0.0105 and pressure = 1.5, values that maximize the response are obtained.
It is important to note that the mass flow and inlet temperature factors do not represent a significant variation in terms of the values in the response since they have a p value > 0.05 (see Table 5), therefore, it is determined that the flow factors mass and inlet temperature are non-significant variables in the search process for the maximum conversion value.Furthermore, the interactions T w × T in , Radius × flow, Pressure × Flow, Pressure × T in and Flow × T in are statistically insignificant.The factors T w , Radius, Pressure and the interactions T w × Radius, T w × Pressure, T w × Pressure and Radius × Pressure are statistically significant.Thus, Figure 6 shows the graph of main effects for the concentration of phthalic anhydride in the reactor outlet.The three points in each graph represent low, medium, and high values for each factor; the medium values were taken from Table 1 (inlet temperature, external reactor wall temperature, operating pressure, reactor radius, and mass flow rate of ortho-xylene), which likely have the greatest influence on the dependent variable (CB); while the high and low values were assigned with reference to the central point.And it is observed that for the condition Tw = 625.15,radius = 0.0105 and pressure = 1.5, values that maximize the response are obtained.
It is important to note that the mass flow and inlet temperature factors do not represent a significant variation in terms of the values in the response since they have a p value > 0.05 (see Table 5), therefore, it is determined that the flow factors mass and inlet temperature are non-significant variables in the search process for the maximum conversion value.Furthermore, the interactions Tw × Tin, Radius × flow, Pressure × Flow, Pressure × Tin and Flow × Tin are statistically insignificant.The factors Tw, Radius, Pressure and the interactions Tw × Radius, Tw × Pressure, Tw × Pressure and Radius × Pressure are statistically significant.Since there are non-significant variables in the process, their value can be managed with other criteria, such as the cost of cooling the gases from a highly exothermic process.Since there are non-significant variables in the process, their value can be managed with other criteria, such as the cost of cooling the gases from a highly exothermic process.
Taking this premise into account, the assigned values for the non-significant factors are set at Flow = 4684 kg/m 2 /h and T in = 625.15K to continue with the analysis.

Response Surface Methodology
After carrying out the factorial design and finding the best configuration of the experimentation levels, the optimal point was search, now applying a response surface design.In Figure 7 is shown for the process.
ChemEngineering 2024, 8, x FOR PEER REVIEW 11 of 24 Taking this premise into account, the assigned values for the non-significant factors are set at Flow = 4684 kg/m 2 /h and Tin = 625.15K to continue with the analysis.

Response Surface Methodology
After carrying out the factorial design and finding the best configuration of the experimentation levels, the optimal point was search, now applying a response surface design.In Figure 7 is shown for the process.Firstly, the steepest ascent method was applied to find the optimum region, which was possible with the linear regression equation (LR1, Equation (31)) derived from the factorial design, where only the significant factors are considered, wall temperature, radius, and pressure.
from this equation, the signs of the coefficients were mainly used and a step size proportional to the magnitude of each factor was set (see Table 6), in order to explore the region before and next to the optimal configuration point.Firstly, the steepest ascent method was applied to find the optimum region, which was possible with the linear regression equation (LR1, Equation (31)) derived from the factorial design, where only the significant factors are considered, wall temperature, radius, and pressure.
from this equation, the signs of the coefficients were mainly used and a step size proportional to the magnitude of each factor was set (see Table 6), in order to explore the region before and next to the optimal configuration point.Then the exploration was carried out with the step size assumed according to Table 7 and considering the fixed mass flow and inlet temperature.The results presented in Figure 8 show the behavior of C B along the steps.Then the exploration was carried out with the step size assumed according to Table 7 and considering the fixed mass flow and inlet temperature.The results presented in Figure 8 show the behavior of CB along the steps.Note that in Figure 8, in the vicinity of step 2 there is a region of local maximum, so we worked with the coordinates of the points close to step 2 to create the response surface design table, according to the Box-Behnken design, in Minitab 18 ® , from which 15 runs were obtained, presented in Table 8, which were simulated with the Polymath ® program.Note that in Figure 8, in the vicinity of step 2 there is a region of local maximum, so we worked with the coordinates of the points close to step 2 to create the response surface design table, according to the Box-Behnken design, in Minitab 18 ® , from which 15 runs were obtained, presented in Table 8, which were simulated with the Polymath ® program.The analysis of the response surface design applied to the phthalic anhydride concentration process gave the results show in Table 9: This analysis was made in Minitab 18 ® to evaluate the significance of the values.From the analysis of variance, it is inferred that the curvature (with a value p =< 0.05) is significant and the factors involved are significant for the design of the response surface.
With the help of the following regression equation of the response surface (LR2, Equation (32)), considering only the factors that vary, the graphs of Figures 9 and 10 are generated keeping one factor constant at a time, which represent the second order interaction of the factors.From the contour and response surface graphs (Figures 9 and 10) it can be that there is no optimal point; but there is an optimal region for the process, w concentration of phthalic anhydride is greater than 0.0001405 kmol m −3 .This me the experimenter has an optimal operating range provided by the modeling and tion, which is reflected in the contour graphs.But within this optimal region, the s performs the calculation and suggests the follow optimal point shown in Figure 1 Reactor  From the contour and response surface graphs (Figures 9 and 10) it can be inferred that there is no optimal point; but there is an optimal region for the process, where the concentration of phthalic anhydride is greater than 0.0001405 kmol m −3 .This means that the experimenter has an optimal operating range provided by the modeling and simulation, which is reflected in the contour graphs.But within this optimal region, the software performs the calculation and suggests the follow optimal point shown in Figure 11.With these optimal results, the parameters of Tables 1 and 2 and the differential equa tions of the reactor presented in Section 2.1, the new concentration profiles of o-xylene phthalic anhydride and oxidation by-products produced in the catalytic bed are deter mined for the ten selected points of the reactor, according to Figures 12 and 13  Then, the optimal results obtained are: With these optimal results, the parameters of Tables 1 and 2 and the differential equations of the reactor presented in Section 2.1, the new concentration profiles of o-xylene, phthalic anhydride and oxidation by-products produced in the catalytic bed are determined for the ten selected points of the reactor, according to Figures 12 and 13.
ChemEngineering 2024, 8, x FOR PEER REVIEW 15 With these optimal results, the parameters of Tables 1 and 2 and the differential e tions of the reactor presented in Section 2.1, the new concentration profiles of o-xy phthalic anhydride and oxidation by-products produced in the catalytic bed are d mined for the ten selected points of the reactor, according to Figures 12 and 13    Note that the optimal phthalic anhydride concentration, both in Equation (32) or LR2 and in the system of differential Equations ( 19), ( 20) and ( 28), ( 29) is 0.00014033 kmol/m 3 and 0.0001405 kmol/m 3 , respectively, which means an error of 0.12%.

Modeling of the Processes
Since the mass transfer process by molecular diffusion in the spherical particle is in the radial direction, and in a steady state; Equation ( 1) is reduced to: Equation ( 2) applies to each of the component's A, B and C respectively, that is: Then Fick's Law is applying: The molecular diffusion heat transfer equation applied to the catalyst particle by simplifying Equation (4) becomes: Now, we use Fourier's law given by: Note that the optimal phthalic anhydride concentration, both in Equation (32) or LR2 and in the system of differential Equations ( 19), ( 20) and ( 28), ( 29) is 0.00014033 kmol/m 3 and 0.0001405 kmol/m 3 , respectively, which means an error of 0.12%.

Studies on the Catalyst Particle 4.1. Modeling of the Processes
Since the mass transfer process by molecular diffusion in the spherical particle is in the radial direction, and in a steady state; Equation ( 1) is reduced to: Equation ( 2) applies to each of the component's A, B and C respectively, that is: Then Fick's Law is applying: The molecular diffusion heat transfer equation applied to the catalyst particle by simplifying Equation (4) becomes: Equations ( 37)-( 39), ( 42), ( 45)-( 47), ( 52) and (54) must be solved simultaneously with the help of Equations ( 48)-(53), with the following boundary conditions: From the Equations ( 19)-( 21): From the Equations ( 37)-(39): From the Equation (54): From the Equation (42):

Simulation for Calculating the Effectiveness Factor in Each Particle
To generate the concentration and temperature profile at each point in the reactor, a run was made with the Polymath using the optimized values of the wall temperature, the pressure, and the radius of the reactor; beside the inlet temperature and feed flow and the other data in Table 1.This allows knowing the external conditions of temperature and concentration of each particle located at a point (r,z) of the reactor.A specific point of the reactor had the following conditions used (T = 625.161K, C A = 1.7486 × 10 −3 kmol/m 3 , C B = 0.2809 × 10 −3 kmol/m 3 , C C = 1.3606 × 10 −3 kmol/m 3 ).With these conditions we proceeded to calculate the effectiveness factor of that specific particle.For the other particles, the point (r,z) is again identified, along with their respective concentrations and temperature, and the effectiveness factor is calculated.It is important to highlight that to find the effectiveness factor, the set of Equations ( 58)-(66) were used for each one of the 78 selected catalyst particles within the reactor bed.
The molar flux density profiles inside the catalyst particle are presented in Figure 14.Note that N A takes negative values since the flow of this component goes in the opposite direction to the reference axis located in the geometric center of the spherical particle.The o-xylene molar flux density is higher at the surface of the particle, compared to what occurs in the center.The N C has a slightly higher value at the surface of the particle, relative to the center.In Figure 15 is shown the molar concentration of components A, B and C inside the catalyst particle in function radius.In this case, the center of the catalyst particle contains more phthalic anhydride than on the surface.The opposite occurs with o-xylene.It is observed that there is a very small temperature gradient between the temperature in the center and the surface of the catalytic particle, according to Figure 16.This is due to the very small radius of the particle.In Figure 15 is shown the molar concentration of components A, B and C inside the catalyst particle in function radius.In this case, the center of the catalyst particle contains more phthalic anhydride than on the surface.The opposite occurs with o-xylene.In Figure 15 is shown the molar concentration of components A, B and C inside the catalyst particle in function radius.In this case, the center of the catalyst particle contains more phthalic anhydride than on the surface.The opposite occurs with o-xylene.It is observed that there is a very small temperature gradient between the temperature in the center and the surface of the catalytic particle, according to Figure 16.This is due to the very small radius of the particle.It is observed that there is a very small temperature gradient between the temperature in the center and the surface of the catalytic particle, according to Figure 16.This is due to the very small radius of the particle.There is an increase in the heat flux density from the center to the particle outside, reaching its maximum value at the end, as observed in Figure 17.There is an increase in the heat flux density from the center to the particle outside, reaching its maximum value at the end, as observed in Figure 17.There is an increase in the heat flux density from the center to the particle outside, reaching its maximum value at the end, as observed in Figure 17.An analysis of the axial profiles of effectiveness factor for the different radial coordinates of the reactor, shows that the factor of effectiveness falls rapidly near the entrance of the reactor and subsequently increases until reaching the maximum point.For example, for the center of the reactor (r0) at the reactor position (z = 0.30), the effectiveness factor is 0.84135 and for the radius of the reactor (r10) for the same position (z = 0.30) η is equal to 0.9342.Already for half of the reactor (z = 1.5) in the center of the reactor (r0) η is equal to 0.97154 and for the radius of the reactor (r10) the effectiveness factor is 0.9763, as shown in Figure 19.Note that the initial drop in the effectiveness factors is due to the increase in temperature (see Figure 4), that is, they are inversely proportional.For a better visualization, only some profiles are shown.

Conclusions
With the application of material and energy balances in the catalytic reactor bed, using operating conditions reported in the literature (shown in Table 1), the production rate of phthalic anhydride has been optimized in terms of concentration as the dependent variable, considering the inlet temperature of the reactive mixture, reactor radius, external reactor wall temperature, operating pressure, and mass flow rate of ortho-xylene as independent variables.An analysis of the axial profiles of effectiveness factor for the different radial coordinates of the reactor, shows that the factor of effectiveness falls rapidly near the entrance of the reactor and subsequently increases until reaching the maximum point.For example, for the center of the reactor (r0) at the reactor position (z = 0.30), the effectiveness factor is 0.84135 and for the radius of the reactor (r10) for the same position (z = 0.30) η is equal to 0.9342.Already for half of the reactor (z = 1.5) in the center of the reactor (r0) η is equal to 0.97154 and for the radius of the reactor (r10) the effectiveness factor is 0.9763, as shown in Figure 19.Note that the initial drop in the effectiveness factors is due to the increase in temperature (see Figure 4), that is, they are inversely proportional.For a better visualization, only some profiles are shown.An analysis of the axial profiles of effectiveness factor for the different radial coordinates of the reactor, shows that the factor of effectiveness falls rapidly near the entrance of the reactor and subsequently increases until reaching the maximum point.For example, for the center of the reactor (r0) at the reactor position (z = 0.30), the effectiveness factor is 0.84135 and for the radius of the reactor (r10) for the same position (z = 0.30) η is equal to 0.9342.Already for half of the reactor (z = 1.5) in the center of the reactor (r0) η is equal to 0.97154 and for the radius of the reactor (r10) the effectiveness factor is 0.9763, as shown in Figure 19.Note that the initial drop in the effectiveness factors is due to the increase in temperature (see Figure 4), that is, they are inversely proportional.For a better visualization, only some profiles are shown.

Conclusions
With the application of material and energy balances in the catalytic reactor bed, using operating conditions reported in the literature (shown in Table 1), the production rate of phthalic anhydride has been optimized in terms of concentration as the dependent variable, considering the inlet temperature of the reactive mixture, reactor radius, external reactor wall temperature, operating pressure, and mass flow rate of ortho-xylene as independent variables.

Conclusions
With the application of material and energy balances in the catalytic reactor bed, using operating conditions reported in the literature (shown in Table 1), the production rate of phthalic anhydride has been optimized in terms of concentration as the dependent variable, considering the inlet temperature of the reactive mixture, reactor radius, external reactor wall temperature, operating pressure, and mass flow rate of ortho-xylene as independent variables.This approach also enabled the determination of the length at which the maximum concentration of phthalic anhydride is achieved, surpassing 3 m (see Figure 3).To enhance the conversion of ortho-xylene with a shorter reactor length, the o-xylene oxidation process was analyzed using DOE with five factors and three levels for each factor.The analysis identified significant variables (factors)-wall temperature, reactor radius, and operating pressure, with a p-value < 0.05 (for a 95% confidence level)-crucial for achieving the maximum conversion value in the reactor.
Utilizing the response surface methodology, the optimal region for the ortho-xylene oxidation process was determined.The best configuration included an outer wall temperature of 630.8 K, reactor radius of 0.011 m, operating pressure of 1.65 atm, mass flow rate of 4684 kg/m 2 /h, and feed temperature of 625.15 K.Under these conditions, the reactor length required to achieve the maximum concentration of phthalic anhydride was found to be 1.5 m (see Figure 12).
The effectiveness factor profiles were obtained for individual catalytic particles at discrete points in the reactor.This allowed establishing a relationship with the temperature profile within the catalytic bed, indicating an inverse correlation between the effectiveness factor and temperature.Given the highly exothermic nature of the reaction, these findings facilitate the appropriate sizing of reactor diameter and length.This, coupled with heat transfer processes, establishes the most favorable conditions for the production of phthalic anhydride from the partial oxidation of ortho-xylene.
The developed method can also be extended to other processes, given the prior knowledge of the kinetic model, operating conditions, physical properties, and geometric dimensions, which must be incorporated into the input data line of the program.Upon execution, it primarily provides the radial and axial profiles of temperature and composition required for calculating the effectiveness factor profile by solving the system of nonlinear algebraic equations.Effective thermal conductivity in the catalytic particle (kJ m −2 s −1 K −1 ) k er : Effective radial conductivity in the catalytic bed (kJ m −1 s −1 K −1 ) k w : Thermal conductivity in the tube wall (kJ m −1 s −1 K −1 ) L: Reactor length (m) N j : Molar flux density based on stationary component coordinates j (kmol m −2 s −1 ) N p: Maximum number of subdivisions in the reactor or maximum number of subdivisions in the particle NR: Number of reactions p j : Partial pressure of each component j (atm) p O 2 : Partial pressure of oxygen (atm) P eMr : Radial mass transfer Peclet number P eHr : Radial heat transfer Peclet number q r : Radial heat flux density (kJ m −2 s −1 ) R g : Universal gas constant (0.082 L atm mol −1 K −1 ) R: Reactor radius or particle radius (m) r: Radial coordinate (m) R i : Reaction velocity per unit volume of catalyst (kmol m −3 cat s −1 ) r i : Reaction velocity per unit weight of catalyst (kmol kg −1 cat s −1 ) t: Time (s) T: Temperature (K) T i : Temperature at each placement point (K) T W : Reactor wall temperature (K) V: Gas superficial velocity (m s −1 ) W i : Arrhenius equation parameter (K) Z: Axial coordinate ∆H Ri : Reaction enthalpy of i-reaction (kJ kmol −1 ) Greek Letters Catalyst Surface w: Wall

Figure 2 .
Figure 2. Scheme of the fixed bed catalytic reactor.

Figure 2 .
Figure 2. Scheme of the fixed bed catalytic reactor.
Feed composition 1% mol O-xylene * 99% mol of dry air Total pressure of feed gases P = 1 atm * Inlet temperature of reaction gases to the reactor T 0 = 625.15K * Inner reactor wall temperature T W = 625.15K * Flow properties and gas phase

Figure 3 .
Figure 3. Axial concentration profile for ten points.Figure 3. Axial concentration profile for ten points.

Figure 3 .
Figure 3. Axial concentration profile for ten points.Figure 3. Axial concentration profile for ten points.The temperature profiles produced in the catalytic bed as a function of the reactor length are presented in Figure4.Note that the reaction temperature increases rapidly near the reactor inlet reaching the maximum temperature.This behavior is the same as that reported in literature [27] for the same conditions of T 0 = T w = 625.15K.

Figure 4 .
Figure 4. Axial profile of temperature for ten points.

Figure 4 .
Figure 4. Axial profile of temperature for ten points.

Figure 6 .
Figure 6.Main effects of experimental factors at low, medium, and high levels on phthalic anhydride concentration.

Figure 6 .
Figure 6.Main effects of experimental factors at low, medium, and high levels on phthalic anhydride concentration.

Figure 7 .
Figure 7. Flowchart for the response surface design.

Figure 7 .
Figure 7. Flowchart for the response surface design.

Figure 12 .
Figure 12.Concentrations profile at the optimum point.

Figure 12 .
Figure 12.Concentrations profile at the optimum point.

Figure 13 .
Figure 13.Temperature profile at the optimum point.

Figure 13 .
Figure 13.Temperature profile at the optimum point.

Figure 14 .
Figure 14.Molar flux density profiles in a catalyst particle.

Figure 15 .
Figure 15.Molar concentration profiles in a catalyst particle.

Figure 14 .
Figure 14.Molar flux density profiles in a catalyst particle.

Figure 15 .
Figure 15.Molar concentration profiles in a catalyst particle.

Figure 15 .
Figure 15.Molar concentration profiles in a catalyst particle.

Figure 16 .
Figure 16.Temperature profiles inside a catalyst particle.

Figure 17 .
Figure 17.Heat flux density inside a catalyst particle.

Figure 18
Figure18shows the particle effectiveness factor profiles as a function of the reactor radius increment for reactor axial positions different.η has only small variations for Z greater than 1.5, so it is not shown in the figure.

Figure
Figure Temperature profiles inside a catalyst particle.

Figure 17 .
Figure 17.Heat flux density inside a catalyst particle.

Figure 18
Figure18shows the particle effectiveness factor profiles as a function of the reactor radius increment for reactor axial positions different.η has only small variations for Z greater than 1.5, so it is not shown in the figure.

Figure 17 .
Figure 17.Heat flux density inside a catalyst particle.

Figure 18
Figure18shows the particle effectiveness factor profiles as a function of the reactor radius increment for reactor axial positions different.η has only small variations for Z greater than 1.5, so it is not shown in the figure.

Figure 18 .
Figure 18.Effectiveness factor profiles as a function of radius for each reactor length.

Figure 19 .
Figure 19.Effectiveness factor profiles as a function of the length for each reactor radio.

Figure 18 .
Figure 18.Effectiveness factor profiles as a function of radius for each reactor length.

Figure 18 .
Figure 18.Effectiveness factor profiles as a function of radius for each reactor length.

Figure 19 .
Figure 19.Effectiveness factor profiles as a function of the length for each reactor radio.

Figure 19 .
Figure 19.Effectiveness factor profiles as a function of the length for each reactor radio.
Gradient vector (1/m) ρ b : Catalyst bed density (kg m −3 ) ρ g : Density of the gas mixture (kg m −3 ) ϕ: Auxiliary variable v ij : Stoichiometric coefficient of component j in reaction i Subscripts e f : Effective i: Number of reactions: i = 1, 2, 3 in: Inlet j: Component: A, B and C k: Number of subdivisions in the reactor or particle: 1, 2, 3, . . .Np s:

Table 1 .
Catalyst data and operating conditions.

Table 2 .
Kinetic and thermodynamic data.

Table 3 .
Factors and levels of experimentation.
Subsequently, the variance analysis was carried out for the data recorded according to the experimental design according to the Table5obtained by Minitab ® 18.

Table 3 .
Factors and levels of experimentation.

Table 4 .
Experimentation block.Subsequently, the variance analysis was carried out for the data recorded according to the experimental design according to the Table5obtained by Minitab ® 18.

Table 6 .
Optimum setting point and step size.

Table 6 .
Optimum setting point and step size.

Table 9 .
Variance analysis of response surface design.
DF = Degrees of freedom.* = Values that could not be calculated.