Entropy Generation Rate Minimization for Methanol Synthesis via a CO2 Hydrogenation Reactor

The methanol synthesis via CO2 hydrogenation (MSCH) reaction is a useful CO2 utilization strategy, and this synthesis path has also been widely applied commercially for many years. In this work the performance of a MSCH reactor with the minimum entropy generation rate (EGR) as the objective function is optimized by using finite time thermodynamic and optimal control theory. The exterior wall temperature (EWR) is taken as the control variable, and the fixed methanol yield and conservation equations are taken as the constraints in the optimization problem. Compared with the reference reactor with a constant EWR, the total EGR of the optimal reactor decreases by 20.5%, and the EGR caused by the heat transfer decreases by 68.8%. In the optimal reactor, the total EGRs mainly distribute in the first 30% reactor length, and the EGRs caused by the chemical reaction accounts for more than 84% of the total EGRs. The selectivity of CH3OH can be enhanced by increasing the inlet molar flow rate of CO, and the CO2 conversion rate can be enhanced by removing H2O from the reaction system. The results obtained herein are in favor of optimal designs of practical tubular MSCH reactors.


Introduction
Over the past century, a huge amount of CO 2 produced through burning fossil fuels has been released into the atmosphere, which has led to global warming. Nowadays, much effort is being put into carbon emission reduction. There are mainly three ways to reduce carbon emissions: (1) Utilizing clean energy sources, such as the solar, wind, nuclear, and tidal energy; (2) utilizing carbon sequestration and storage technology; (3) collecting and recycling CO 2 through chemical reactions. The methanol synthesis via CO 2 hydrogenation (MSCH) reaction is in fact an effective scheme for alleviating the greenhouse effect. What is more, methanol is a primary liquid petrochemical, and is widely used in the chemical and energy fields in applications such as hydrogen storage, dimethyl ether and hydrocarbon production, etc. [1][2][3][4].
where ∆ r H i is the enthalpy of reaction i.

Reactor Model
Two routes of methanol synthesis can be found in public literature, i.e., (1) CO2 shifts to CO through RWGS, and then CO convert into CH3OH; (2) CO2 convert into CH3OH through an intermediate HCOO [3]. Undergo the experimental research of many years, Skrzypel et al. [54] and Vanden Bussche and Froment [55] confirmed that CO2 is the main source of methanol synthesis, and CO and CO2 can mutual converts through the RWGS reaction [2,9]. Two reactions occurring in the reactor are as follows [9,36,44]: where  r i H is the enthalpy of reaction i .
In the chemical industrial production, two reactions generally reacts in a suit of catalyst-filled tubes, which are placed inside a tubular heat exchanger. In order to simplify the model, this paper will study one such tube, and establish a one-dimensional pseudo-homogeneous plug flow reactor model as shown in Figure 1 [15,35,[41][42][43][44][45][46][47][48][49][50][51][52][53]. According to the judgmental criterion of the assumption on the plug flow, one reactor can be taken as a plug flow reactor, in the case of ti / L d (L and dti are, respectively, the reactor length and diameter) exceeds 150 [56]. Herein, ti / L d is equal to 300, therefore, the MSCH reactor can be taken as a plug flow reactor. The mathematical model of diffusion listed in Appendix A, Appendix B, and Appendix C are consistent with recent treatments of reaction-diffusion descriptions by Niven [57] and England [58] and a recent temperature-dependent description by Ritchie [59]. The effect of the diffusion within the catalyst pellet on the chemical reaction rate can be ignored, when the effectiveness factor of chemical reaction is close to 1. In this paper, the diffusion within the catalyst pellet can be ignored, and the MSCH reactor can be taken as a pseudo-homogeneous reactor through the verification. The detailed discussion is shown in Appendices A and B. The MSCH reactor herein can be studied with one-dimensional model, since the size of MSCH reactor is not large, and radial gradient of temperature and mass can be ignored. In summary, the MSCH reactor model herein includes following assumptions: (1) the reaction mixture is not back-mixed in the axial, and mixes uniformly in the radial; (2) the radial temperature and concentration gradients are neglected.
In the MSCH reactor, both the reactants and products are taken as the ideal gas. The reactor geometry parameters, the ICI 51-2 Cu/ZnO/Al2O3 catalyst geometry parameters, the physical parameters of mixture gas and the inlet conditions derived from the reference reactor in References [44,45] are listed in Table 1. The overall heat transfer profile of the reference reactor along the dimensionless axial can be obtained based on the empirical formula proposed by According to the judgmental criterion of the assumption on the plug flow, one reactor can be taken as a plug flow reactor, in the case of L/d ti (L and d ti are, respectively, the reactor length and diameter) exceeds 150 [56]. Herein, L/d ti is equal to 300, therefore, the MSCH reactor can be taken as a plug flow reactor. The mathematical model of diffusion listed in Appendix A, Appendix B, and Appendix C are consistent with recent treatments of reaction-diffusion descriptions by Niven [57] and England [58] and a recent temperature-dependent description by Ritchie [59]. The effect of the diffusion within the catalyst pellet on the chemical reaction rate can be ignored, when the effectiveness factor of chemical reaction is close to 1. In this paper, the diffusion within the catalyst pellet can be ignored, and the MSCH reactor can be taken as a pseudo-homogeneous reactor through the verification. The detailed discussion is shown in Appendices A and B. The MSCH reactor herein can be studied with one-dimensional model, since the size of MSCH reactor is not large, and radial gradient of temperature and mass can be ignored. In summary, the MSCH reactor model herein includes following assumptions: (1) the reaction mixture is not back-mixed in the axial, and mixes uniformly in the radial; (2) the radial temperature and concentration gradients are neglected.
In the MSCH reactor, both the reactants and products are taken as the ideal gas. The reactor geometry parameters, the ICI 51-2 Cu/ZnO/Al 2 O 3 catalyst geometry parameters, the physical parameters of mixture gas and the inlet conditions derived from the reference reactor in References [44,45] are listed in Table 1. The overall heat transfer profile of the reference reactor along the dimensionless axial can be obtained based on the empirical formula proposed by Dixonge [60,61], and the calculation results indicate that the overall heat transfer coefficient is approximately equal to 60 W/ K·m 2 .

Reaction Kinetic Model
The MSCH reaction processes can be classified as the high-pressure process, 25-30 MPa, the medium-pressure process, 10-25 MPa, and the low-pressure process, 5-10 MPa [36]. Among them, there are more researches on low-pressure process, since the low-pressure of methanol synthesis process is applied widely in chemical industries [1,3,62].
In this paper, the kinetic model established by Vanden Bussche and Froment [44] is selected and utilized, since this kinetic model is based on the ICI 51-2 Cu/ZnO/Al 2 O 3 catalyst, which is applied widely in chemical industries, such as the LURGI type methanol synthesis reactor [9]. Another reason is that the kinetic model proposed by Reference [44] has been checked by experiments of lab-scale [44,55] and commercial-scale [9], what is more, the kinetic model has a wider application range. The temperature varies between 180 and 280 • C, the pressures varies between 15 and 51 bar in the experiment of Vanden Bussche and Froment [44]. The pressure from industrial data of LURGI type methanol synthesis reactor is 66.7 bar [9].
According to the kinetic model proposed by Reference [44], the reaction rates of the MSCH and RWGS reactions are: where r i is the reaction rate of reaction i, P k = (PF k )/F T is the partial pressure of component k, P and T are the pressure and temperature of reaction mixture (TRM), respectively, F k and F T are the mole flow rate of component k and reaction mixture in the axial position z, respectively, K * 1 and K * 2 are thermodynamic equilibrium constants of Reactions (1) and (2) [63], respectively, and κ j is parameter group j related to the adsorption, equilibrium, and rate constants of elementary reactions. These parameter groups can be calculated as follows [44]: where R g = 8.314 J/(mole·K) is the universal gas constant, A(i) is the frequency factor, B(i) represents either E or −∆H, and these factors are listed in Table 2 [44].
124,119 ---* k 5a and k 1 are, respectively, the rate constants of the rate-determining steps for the MSCH and RWGS reactions, K H 2 and K H 2 O are, respectively, the adsorption constants of H 2 and H 2 O, K m s (m = 2, 3,4,8,9) are the equilibrium constants of all elementary reactions except for some ignored.

Conservation Equation
The heat transfer between the reaction mixture and the exterior wall heat reservoir obeys the Newton heat transfer law, i.e., q ∝ (∆T), and the heat flux passed through the tube wall is where U is the overall heat transfer coefficient, and T a is the EWT. The energy conservation equation is: where ε c is the void fraction of the catalyst bed, ρ c is the catalyst density, C p,k is the mole heat capacity at constant pressure of the component k, and ∆ r H T,i is the standard mole enthalpy of reaction i. These parameters can be obtained as follows [64]: where υ k,i is the stoichiometric number of component k in reaction i, I i is the integration constant of reaction i that can be obtained based on the standard mole enthalpy ∆ f H 0 298.15K, k and the coefficients A, B, C, D, E, ∆ f H 0 298.15K and M k are listed in Table 3 [65]. The momentum equation utilized to describe the pressure drop along the reactor axial is generally described by the Ergun equation, in the case of Re p /(1 − ε c ) < 500 [66]. The Reynolds numbers Re p is: (15) where G is the superficial mass flow rate, G = ∑ k (F k M k )/A c , M k is the mole mass of component k, µ mix is the viscosity of reaction mixture (see Appendix C), and d p is the catalyst particle diameter. Re p /(1 − ε c ) is less than 0.005 in the case discussed in this paper, therefore, the momentum equation can be described by the Ergun equation as follows: where c g = F T R g T/(PA c ) is the superficial velocity of mixture gas. The mole balance equations are described by the yields of CH 3 OH and CO (ξ 1 and ξ 2 ), which are both defined with the inlet mole flow rate of CO 2 [15]: where subscript 'in' denotes the inlet state of variables, and the mole flow rate of component k and the total mole flow rate are defined as follows:

Entropy Generation Rate of the MSCH Reactor
Non-equilibrium thermodynamics indicates that the irreversible process always generates EGR, which is described as the product sum of the conjugate fluxes and forces [67][68][69]. In a plug flow reactor, the EGR is mainly produced by the chemical reactions, heat transfer and viscous flow [15,[67][68][69].
In term of the EGR due to chemical reactions, the driving forces are: where ∆ r G i is the Gibbs free energy of reaction i. In term of the EGR due to heat transfer, the driving force is: In term of the EGR due to viscous flow, the driving force is: The local EGR is: where the three terms on right side are the local EGR due to heat transfer, σ H , MSCH reaction, σ MS , RWGS reaction, σ RW , and viscous flow, σ F , respectively, and the total EGR, (dS/dt) TOT obtained by the integral of the local EGR along the reactor axial is as follows: where ∑ H , ∑ MS , ∑ RW and ∑ F are, respectively, the total EGR due to heat transfer, MSCH reaction, RWGS reaction and viscous flow. The optimization objective is: Equation (27) takes into account almost all phenomena generating EGR in the MSCH reactor, while Reference [43] only considered the EGR from the MSCH reaction. Therefore, the study for the optimal reactor herein is more comprehensive, and the optimal results can be favorable for the optimal design for the practical MSCH reactor.

Mathematical Description of the Optimization Problem
The optimization problem herein is to minimize the total EGR subjected to some constraints. The fixed methanol yield (ξ 1 = 1.3139, ∆F CH 3 OH = 1.32 × 10 −4 ) and the conservation equations are taken as the constraints. The control variable, i.e., the EWT can be controlled completely. The inlet PRM and the inlet mole flow rate of components in the optimal reactor are same with those in the reference reactor [43]. At the outlet, all variables except for the methanol yield are allowed to change freely. The geometry sizes and the catalyst properties are same with those of the reference reactor in Reference [43]. Optimal control theory is utilized to solve this optimization problem and to find the optimal EWT profile.

Application of Optimal Control Theory
In the early 1950s, optimal control theory had been applied in the study of the minimum time control problem. However, the original optimal control theory cannot solve an optimization problem that the admissible control belongs to a closed set. In order to solve this problem, Pontryagin created the minimum principle on 1958 [70].
The purpose of this paper is to minimize the total EGR of the MSCH reactor and find the optimal EWT profile by using optimal control theory. Therefore, the EWT is taken as the control variable, the total EGR is taken as the performance objective, and the TRM (T), PRM (P), methanol yield (ξ 1 ), and carbon monoxide yield (ξ 2 ) are taken as the state variables controlled by the conservation equations. The inlet and outlet values of state variables and the multiplier functions corresponding to the state variables are taken as the boundary condition [15].
Hamiltonian also needs to satisfy the extremum condition, i.e., ∂H/∂T a = 0. The relation between the EWT and the TRM can be derived via this extremum condition [15]: The boundary conditions of optimal control theory are obtained based on the transversal condition, i.e., x n = x n,specification or λ n = 0. The boundary conditions are listed in Table 4. The property of Hamilton [70] is often used to check the veracity of the optimal results. Hamiltonian herein keeps on a constant along the optimal trajectory, since it does not depend on explicitly the axial reactor length, i.e., Hamiltonian is autonomous. Table 4. Boundary conditions of optimal control problems *.
are, respectively, the inlet and outlet values of state variables, λ(0) and λ(z) are, respectively, the multiplier functions of the corresponding state variables at the inlet and outlet.

Numerical Calculations of Optimization Problem
Solving the optimization problem must has a reasonable initial value. The calculation results of the reference reactor with constant heat reservoir T a = 523 K [43] solved by the numerical discretization are taken as the initial values of the optimal reactor. Optimization problem herein involves eight differential equations, i.e., (Equations (29)-(36)), 1 algebraic equation, i.e., Equation (37), and eight boundary conditions listed in Table 4. Therefore, the optimal control comes down to solve a two-point boundary value problem involving the differential equation. Optimal solution can be obtained by the 'bvp4c' solver in Matlab. The solution accuracy is mainly influenced by the initial multiplier functions, i.e., λ ξ 1 , λ ξ 2 , and the number of grid points. In order to improve the solution accuracy, this paper utilizes 5000 grid points. The calculation results show that the solution error does not exceed 6.0 × 10 −8 .

Numerical Results and Discussions
The reference reactor with the constant heat reservoir T a = 523 K and the optimal reactor with minimum EGR due to the MSCH reaction [43], i.e., the optimal reactor in [43] are both utilized to compare with the optimal reactor herein, i.e., the optimal reactor. In the three reactors, the geometry sizes of reactor, the catalyst properties, and the inlet components are the same with those of the reference reactor in Reference [43]. Table 5 lists the total EGRs of the three reactors. Compared with the results of the reference reactor in [43], the total EGR of the optimal reactor decreases by 20.5%, the total EGR due to the heat transfer decreases by 68.8%, the total EGR due to two chemical reactions decreases by 3.3%, and the total EGR due to the viscous flow is almost unchanged. In the reference and optimal reactors, the total EGRs due to the MSCH reaction account for 69% and 84% of the total EGR, respectively. In the optimal reactor of [43], the total EGR due to the MSCH reaction only accounts for 7.1% of the total EGR, however, the total EGR increases by 123.9% compared with that of the reference reactor.
* The symbols T a = const, dS 1 chem /dt = min, and dS T /dt = min represent the reference reactor and optimal reactor in Reference [45], and the optimal reactor herein, respectively. Figure 2 describes the TRM profiles of the reference and optimal reactors, and the equilibrium temperature profile of the MSCH reaction in reference reactor. As shown in Figure 2, the TRM in the reference reactor (dashed line) increases from 493.2 K at the reactor inlet up to a hot spot of 552.6 K, and then it decreases toward a value of 532.4 K. The TRM in the optimal reactor (solid line) increases from 497.9 K at the reactor inlet up to a hot spot of 555.1 K, next it decreases toward a value of 533.4 K, and then it increases toward a value of 533.5 K at the reactor outlet. Except for the range near the reactor outlet, the TRM profiles in the two reactors are almost the same. The above results can be explained according to the theory of the chemical reaction equilibrium, i.e., Le Chatelier's principle. As shown in Figure 2, the equilibrium temperature of the MSCH reaction (the calculation results based on Equation (22)) decreases steeply from 783.6 K at the reactor inlet up to a turning point of 588 K at the dimensionless axial position z = 0.25, and then levels off to a more flat decrease until a value of 540.9 K. Near to the inlet of the two reactors, the difference between the TRMs and equilibrium temperature is so large that the chemical driving force is large and the reaction rate is also large. Therefore, a large amount of heat is released by the MSCH reaction, which results into the steep increase of the TRM. With the decrease of the difference between the TRM and equilibrium temperature, the EWTs (as shown in Figure 3) become the dominating contribution to the TRMs. Therefore, the TRMs show a flat decrease under the cooling effect of the EWTs. Figure 3 describes the EWT profiles of the reference and optimal reactors along the reactor axial. It is noteworthy that the optimal EWT profile (solid line) shown in Figure 3 is a significant result herein. In industrial, this ideal EWT profile can be reached approximately through assigning some heat exchanger providing different constant cooling temperature, e.g., two-stage cooling strategy (dash dot line) in this case can reduce the total EGR 10.78% with respect to the reference reactor. As shown in Figure 3, the EWT in reference reactor (dash line) shows a constant profile along the reactor axial. The EWT in the optimal reactor increases from 497.9 K up to a local maximum of 539.2 K at the dimensionless axial position z = 0.29, next it decreases to a local minimum of 526.7 K at the dimensionless axial position z = 0.78, and then it increases to a maximum of 533.6 K at the reactor outlet. 588 K at the dimensionless axial position   0.25 z , and then levels off to a more flat decrease until a value of 540.9 K. Near to the inlet of the two reactors, the difference between the TRMs and equilibrium temperature is so large that the chemical driving force is large and the reaction rate is also large. Therefore, a large amount of heat is released by the MSCH reaction, which results into the steep increase of the TRM. With the decrease of the difference between the TRM and equilibrium temperature, the EWTs (as shown in Figure 3) become the dominating contribution to the TRMs. Therefore, the TRMs show a flat decrease under the cooling effect of the EWTs.    Figure 3 describes the EWT profiles of the reference and optimal reactors along the reactor axial. It is noteworthy that the optimal EWT profile (solid line) shown in Figure 3 is a significant result herein. In industrial, this ideal EWT profile can be reached approximately through assigning some heat exchanger providing different constant cooling temperature, e.g., two-stage cooling strategy (dash dot line) in this case can reduce the total EGR 10.78% with respect to the reference reactor. As shown in Figure 3, the EWT in reference reactor (dash line) shows a constant profile along the reactor axial. The EWT in the optimal reactor increases from 497.9 K up to a local maximum of 539.2 K at the dimensionless axial position   0.29 z , next it decreases to a local minimum of 526.7 K at the dimensionless axial position   0.78 z , and then it increases to a maximum of 533.6 K at the reactor outlet. Figure 4 describes that the optimal EWT (solid line) and TRM (dashed line) profiles in the optimal (large frame) and reference reactors (small frame). As shown in large frame of Figure 4, the EWT is almost all less than the TRM. According to the Le Chatelier's principle, decreasing the TRM is favor of increasing the carbon dioxide conversion and methanol production rates for the endothermic reaction. Therefore, the EWT should be lower than the TRM. It is noteworthy that the EWTs are equal to the TRMs at the inlet and outlet in the optimal reactor. The reason for this phenomenon is that the EWT is equal to the TRM, in the case of the TRM changes freely at the boundary according to optimal control theory.  Figure 3. The EWT profiles for the reference, optimal, and two-stage cooling reactor, T a ( z). Figure 4 describes that the optimal EWT (solid line) and TRM (dashed line) profiles in the optimal (large frame) and reference reactors (small frame). As shown in large frame of Figure 4, the EWT is almost all less than the TRM. According to the Le Chatelier's principle, decreasing the TRM is favor of increasing the carbon dioxide conversion and methanol production rates for the endothermic reaction. Therefore, the EWT should be lower than the TRM. It is noteworthy that the EWTs are equal to the TRMs at the inlet and outlet in the optimal reactor. The reason for this phenomenon is that the EWT is equal to the TRM, in the case of the TRM changes freely at the boundary according to optimal control theory.
As shown in Figure 4, the temperature difference in the optimal reactor between the TRM and EWT is almost all less than that in the reference reactor. Therefore, the local EGR due to the heat transfer in the optimal reactor is less than that in the reference reactor (as shown in Figure 5). As shown in Figures 5 and 6, compared with the reference reactor, the local EGR and driving force due to the heat transfer in the optimal reactor distribute more smooth which is in accordance with the equipartition principle of entropy generation and driving force [15,[71][72][73][74].
EWT is almost all less than the TRM. According to the Le Chatelier's principle, decreasing the TRM is favor of increasing the carbon dioxide conversion and methanol production rates for the endothermic reaction. Therefore, the EWT should be lower than the TRM. It is noteworthy that the EWTs are equal to the TRMs at the inlet and outlet in the optimal reactor. The reason for this phenomenon is that the EWT is equal to the TRM, in the case of the TRM changes freely at the boundary according to optimal control theory.  As shown in Figure 4, the temperature difference in the optimal reactor between the TRM and EWT is almost all less than that in the reference reactor. Therefore, the local EGR due to the heat transfer in the optimal reactor is less than that in the reference reactor (as shown in Figure 5). As shown in Figures 5 and 6, compared with the reference reactor, the local EGR and driving force due to the heat transfer in the optimal reactor distribute more smooth which is in accordance with the equipartition principle of entropy generation and driving force [15,[71][72][73][74].   As shown in Figure 4, the temperature difference in the optimal reactor between the TRM and EWT is almost all less than that in the reference reactor. Therefore, the local EGR due to the heat transfer in the optimal reactor is less than that in the reference reactor (as shown in Figure 5). As shown in Figures 5 and 6, compared with the reference reactor, the local EGR and driving force due to the heat transfer in the optimal reactor distribute more smooth which is in accordance with the equipartition principle of entropy generation and driving force [15,[71][72][73][74].     Figure 7, more than 75% of the total EGR distributes in the first 30% length of the two reactors. The local EGR of the optimal reactor (solid line) decreases from 0.2911 W/(K·m) up to a local minimum of 0.2623 W/(K·m) at the dimensionless axial position z = 0.029, next it increases toward a value of 0.3025 W/(K·m) at the dimensionless axial position z = 0.115, next it decreases steeply toward a value of 0.028 W/(K·m) at the dimensionless axial position z = 0.33, and then it levels off to a more flat decrease until a value of 0.002 W/(K·m) at the reactor outlet. Compared with the local EGR of the reference reactor (dash line), the local EGR of the optimal reactor distribute more even, especially after z = 0.33. At the first 30% length of the two reactors, the dominating contribution for the local EGR comes from the MSCH reaction, after 30% length, the dominating contribution for the local EGR comes from the heat transfer.   Figure 8, the contribution of the MSCH reaction on local EGRs are maximal in the two reactors. After the dimensionless axial position   0.35 z , the dominating contribution for the local EGR comes from the heat transfer. The local EGR caused by the viscous flow is minimal, since the size of the lab-scale reactor herein is so small that the pressure drop is negligible. However, the local EGR due to viscous flow cannot be ignored in practical methanol synthesis reactor. It is noteworthy that compared with other profiles, the local EGR profiles due to the heat transfer change most dramatically. The reason of this phenomenon is that under the constraints of fixed methanol yield and fixed reactor, the EGR due to heat transfer can be minimized in a larger optimization potential.  Figure 9 describes the methanol yields of the two reactors. As shown in Figure 9, the profiles of   Figure 8, the contribution of the MSCH reaction on local EGRs are maximal in the two reactors. After the dimensionless axial position z = 0.35, the dominating contribution for the local EGR comes from the heat transfer. The local EGR caused by the viscous flow is minimal, since the size of the lab-scale reactor herein is so small that the pressure drop is negligible. However, the local EGR due to viscous flow cannot be ignored in practical methanol synthesis reactor. It is noteworthy that compared with other profiles, the local EGR profiles due to the heat transfer change most dramatically. The reason of this phenomenon is that under the constraints of fixed methanol yield and fixed reactor, the EGR due to heat transfer can be minimized in a larger optimization potential.   Figure 8, the contribution of the MSCH reaction on local EGRs are maximal in the two reactors. After the dimensionless axial position   0.35 z , the dominating contribution for the local EGR comes from the heat transfer. The local EGR caused by the viscous flow is minimal, since the size of the lab-scale reactor herein is so small that the pressure drop is negligible. However, the local EGR due to viscous flow cannot be ignored in practical methanol synthesis reactor. It is noteworthy that compared with other profiles, the local EGR profiles due to the heat transfer change most dramatically. The reason of this phenomenon is that under the constraints of fixed methanol yield and fixed reactor, the EGR due to heat transfer can be minimized in a larger optimization potential.  Figure 9 describes the methanol yields of the two reactors. As shown in Figure 9, the profiles of the methanol yields of the two reactors are similar, since the outlet methanol yields of the two reactors are fixed. The inlet methanol yield is slightly larger than that of the reference reactor, since the initial MSCH reaction rate in the optimal reactor is slightly larger than that in the reference reactor.  Figure 9 describes the methanol yields of the two reactors. As shown in Figure 9, the profiles of the methanol yields of the two reactors are similar, since the outlet methanol yields of the two reactors are fixed. The inlet methanol yield is slightly larger than that of the reference reactor, since the initial MSCH reaction rate in the optimal reactor is slightly larger than that in the reference reactor.  Figure 10 describes mole flow rates of components in the two reactors. The mole flow rate of N2 is not given, since N2 doesn't participate in the MSCH and RWGS reactions. As shown in Figure 10, the mole flow rate of CO decreases over the overall reactor, since the concentration of CO in reaction mixture exceeds the equilibrium concentration of CO under the TRM, which results in that the reaction direction of RWGS is changed. Undoubtedly, the reverse reaction of the RWGS reaction is beneficial for improving the reaction selectivity. The mole flow rate of H2O increases over the overall reactor. During the reaction process, with the increase of H2O, the MSCH reaction tends to reach the chemical equilibrium, which goes against the increase of methanol yield.

Conclusions
This paper studies the performance of MSCH reactor and obtains the optimal configuration of the EWT using FTT theory. The minimum EGR caused by heat transfer, viscous flow and chemical reaction is taken as the optimization objective, the fixed methanol yield, the fixed inlet pressure, and the fixed inlet components are taken as the constraints, and the completely controllable EWT is taken as the control variable. The mathematical model of the optimization problem is established using optimal control theory. The optimal results indicate that the total EGR decreases by 20.5% and the EGR due to heat transfer decreases by 68.8% compared with that of the reference reactor with a constant EWT profile. The local EGRs of the two reactors mainly distribute in the first 30% reactor  is not given, since N 2 doesn't participate in the MSCH and RWGS reactions. As shown in Figure 10, the mole flow rate of CO decreases over the overall reactor, since the concentration of CO in reaction mixture exceeds the equilibrium concentration of CO under the TRM, which results in that the reaction direction of RWGS is changed. Undoubtedly, the reverse reaction of the RWGS reaction is beneficial for improving the reaction selectivity. The mole flow rate of H 2 O increases over the overall reactor. During the reaction process, with the increase of H 2 O, the MSCH reaction tends to reach the chemical equilibrium, which goes against the increase of methanol yield.  The mole flow rate of N2 is not given, since N2 doesn't participate in the MSCH and RWGS reactions. As shown in Figure 10, the mole flow rate of CO decreases over the overall reactor, since the concentration of CO in reaction mixture exceeds the equilibrium concentration of CO under the TRM, which results in that the reaction direction of RWGS is changed. Undoubtedly, the reverse reaction of the RWGS reaction is beneficial for improving the reaction selectivity. The mole flow rate of H2O increases over the overall reactor. During the reaction process, with the increase of H2O, the MSCH reaction tends to reach the chemical equilibrium, which goes against the increase of methanol yield.

Conclusions
This paper studies the performance of MSCH reactor and obtains the optimal configuration of the EWT using FTT theory. The minimum EGR caused by heat transfer, viscous flow and chemical reaction is taken as the optimization objective, the fixed methanol yield, the fixed inlet pressure, and the fixed inlet components are taken as the constraints, and the completely controllable EWT is taken as the control variable. The mathematical model of the optimization problem is established using optimal control theory. The optimal results indicate that the total EGR decreases by 20.5% and the EGR due to heat transfer decreases by 68.8% compared with that of the reference reactor with a constant EWT profile. The local EGRs of the two reactors mainly distribute in the first 30% reactor The effectiveness factors η 1 and η 2 are taken as the value of 1. The reason of this phenomenon is that the complex MSCH reactor model involving the catalyst model, which cannot be optimized based on the optimal control. Therefore, the MSCH reactor is assumed to be a pseudo-homogeneous reactor. The reasonableness of η 1 = 1 and η 2 = 1 herein will be verified based on the MSCH reference reactor considering the catalyst pellet model.
The effectiveness factor η i of reaction i is defined as follow [75] where η i is the effectiveness factor of reaction i, R p is the radius of the catalyst pellet, and r i and r p, i are the reaction rates in the catalyst bed and pellet, respectively. The r i and r p, i depend on the partial pressure in the catalyst bed P k and pellet P p, k , respectively. The partial pressure profiles inside the catalyst pellet can be obtained based on the component balances equations as follows: where D e m,k is the effective diffusivities. The boundary conditions at the pellet surface and at the center of the catalyst are as follows: Appendix B Figure A1 describes the partial pressures of CO 2 and H 2 O along the dimensionless radial coordinate of catalyst pellet at the dimensionless axial position z = 0.1. As shown in Figure A1, the partial pressure of CO 2 inside the catalyst pellet is less than that in the catalyst bed, i.e., P p, CO 2 ≤ P CO 2 , when R ≤ 1, while the partial pressure of H 2 O inside the catalyst pellet is more  Figure A2 describes the effectiveness factors in the reference reactor considering the catalyst model. As shown in Figure A2, the effectiveness factors of two reactions are close to the constant of 1, especially the effectiveness factor of the MSCH reaction (solid line). Table A1 lists the effect of the effectiveness factors on the key parameters of the reference reactor. As shown in Table A1, the effectiveness factors  1 (z) and  2 (z) in the reference reactor can be taken as 1. Therefore, the assumption of pseudo-homogeneous reactor is reasonable. It is noteworthy that the effectiveness factor of the RWGS reaction (dash line) is more than 1 between   0 z and   0.27 z . The reason of this phenomenon is that, the reverse reaction rate of RWGS inside the catalyst pellet is more quickly than that in the catalyst bed, since the partial pressure of H2O inside the catalyst pellet is more than that in the catalyst bed.   Figure A2 describes the effectiveness factors in the reference reactor considering the catalyst model. As shown in Figure A2, the effectiveness factors of two reactions are close to the constant of 1, especially the effectiveness factor of the MSCH reaction (solid line). Table A1 lists the effect of the effectiveness factors on the key parameters of the reference reactor. As shown in Table A1, the effectiveness factors η 1 (z) and η 2 (z) in the reference reactor can be taken as 1. Therefore, the assumption of pseudo-homogeneous reactor is reasonable. It is noteworthy that the effectiveness factor of the RWGS reaction (dash line) is more than 1 between z = 0 and z = 0.27. The reason of this phenomenon is that, the reverse reaction rate of RWGS inside the catalyst pellet is more quickly than that in the catalyst bed, since the partial pressure of H 2 O inside the catalyst pellet is more than that in the catalyst bed. . The reason of this phenomenon is that the main reaction, i.e., MSCH reaction consumes CO2 and simultaneously produces H2O, and the diffusion rate of H2O is less than the formation rate of H2O inside the catalyst pellet.  Figure A2 describes the effectiveness factors in the reference reactor considering the catalyst model. As shown in Figure A2, the effectiveness factors of two reactions are close to the constant of 1, especially the effectiveness factor of the MSCH reaction (solid line). Table A1 lists the effect of the effectiveness factors on the key parameters of the reference reactor. As shown in Table A1, the effectiveness factors  1 (z) and  2 (z) in the reference reactor can be taken as 1. Therefore, the assumption of pseudo-homogeneous reactor is reasonable. It is noteworthy that the effectiveness factor of the RWGS reaction (dash line) is more than 1 between   0 z and   0.27 z . The reason of this phenomenon is that, the reverse reaction rate of RWGS inside the catalyst pellet is more quickly than that in the catalyst bed, since the partial pressure of H2O inside the catalyst pellet is more than that in the catalyst bed.   * where the symbol, η 1 = η 2 = 1 and η 1 (z), η 2 (z) represent the reference reactor model herein and the reference reactor model considering the catalyst pellet model, respectively.