Minimum Entropy Generation Rate and Maximum Yield Optimization of Sulfuric Acid Decomposition Process Using NSGA-II.

Based on the theory of finite-time thermodynamics (FTT), the effects of three design parameters, that is, inlet temperature, inlet pressure, and inlet total mole flow rate, of a tubular plug-flow sulfuric acid decomposition reactor on the total entropy generation rate (EGR) and SO2 yield are analyzed firstly. One can find that when the three design parameters are taken as optimization variables, the minimum total EGR and the maximum SO2 yield of the reference reactor restrict each other, i.e., the two different performance objectives cannot achieve the corresponding extremum values at the same time. Then, the second-generation non-dominated solution sequencing genetic algorithm (NSGA-II) is further used to pursue the minimum total EGR and the maximum SO2 yield of the reference reactor by taking the three parameters as optimization design variables. After the multi-objective optimization, the reference reactor can be Pareto improved, and the total EGR can be reduced by 9% and the SO2 yield can be increased by 14% compared to those of the reference reactor. The obtained results could provide certain theoretical guidance for the optimal design of actual sulfuric acid decomposition reactors.

1. Introduction

At present, the Hybrid-Sulphur (H-S) thermochemical cycle and the Sulphur-Iodine (S-I) thermochemical cycle are considered to be the two most promising recycling methods in the preparation of hydrogen from water by thermochemical cycles [1], and the schematic diagram of S-I thermochemical cycle is shown in Figure 1. Both the H-S and the S-I cycles contain the sulfuric acid decomposition process. Therefore, it is important and necessary to improve the performance of the sulfuric acid decomposition process.

Figure 1

The schematic diagram of S-I thermochemical cycle.

The S-I thermochemical cycle consists of three main chemical reactions: (1) the endothermic decomposition of hydrogen iodide in gas phase; (2) the spontaneous absorption of sulfur dioxide in liquid phase; (3) the sulfuric acid decomposition reaction. The corresponding reaction equations are given as follows:

Reaction type (I) is the spontaneous decomposition of H2SO4 into SO3 and H2O at 400–500 °C. Reaction type (II) is the reaction of SO3 over 750 °C to produce SO2 and O2 under the action of a catalyst. In this process, a great deal of heat is consumed, which is also the main energy consumption process in the S-I and H-S thermochemical cycles.

In the aspect of thermodynamic analysis and optimization of sulfuric acid decomposition, Van der ham et al. [1] assumed that the reaction mixture satisfies the ideal gas equation of state, established the physical model of sulfuric acid decomposition reaction, and analyzed the minimization of entropy generation rate (EGR) of a sulfuric acid decomposition reactor by using the optimal control theory. Kuchi et al. [2] carried out a numerical simulation of a high-temperature shell and tube heat exchanger and decomposer, investigated the fluid flow, heat transfer, and chemical reaction processes in the decomposer by using the porous media method, and established a two-dimensional axisymmetric tubular plug-flow reactor model. Ponyavin et al. [3] studied the sulfuric acid decomposer process in a high-temperature ceramic heat exchanger and established a three-dimensional calculation model of the reactor. Van der ham et al. [4] further compared two methods to improve the efficiency of sulfuric acid decomposition reactor and proposed two design schemes to improve the efficiency of the reactor. On the basis of Ref. [1], Wang et al. [5,6] optimized the decomposition of sulfuric acid in the tubular plug-flow reactor with the goal of maximum yield [5], further analyzed the influences of the design parameters of the reactor on the SO2 yield and specific EGRs [6], and obtained the optimal parameters corresponding to the minimum specific EGRs.

Many scholars have optimized other types of thermochemical reaction processes by using the theory and method of finite-time thermodynamics (FTT) [7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22]. For example, Wang et al. [23] investigated the isotherm chemical reaction A⇔B⇔C and obtained the best concentration configuration of the reaction. Johannessen and Kjelstrup [24] studied the EGR minimization of sulfur dioxide oxidation process. The second-generation non-dominated solution sequencing genetic algorithm (NSGA-II) has been widely used in multi-objective optimization of various engineering problems [25,26,27,28,29,30].

On the basis of Refs. [1,5,6], this paper will further analyze the effects of reactant inlet temperature, pressure, and total molar flow rate on total EGR and SO2 yield, and perform the multi-objective optimization of the process by using the NSGA-II algorithm by applying FTT.

2. Modeling of the Sulfuric Acid Decomposition Process

A reference reactor used in the performance analysis and optimization as well as the kinetics and thermodynamics models will be introduced in this section.

2.1. Reference Reactor

The model of a tubular plug-flow reactor for sulfuric acid decomposition is shown in Figure 2. It is assumed that the temperature (Tw) of the outer wall of tubular plug-flow reactor does not change with time and its distribution is linear along the axial direction of the reactor. The distribution follows (K). The reaction mixture in the reactor is regarded as an ideal gas and only flows along the axial direction of the reactor. The radial concentration gradient and temperature gradient of the reaction mixture in the reactor are ignored without both radial diffusion and back-mixing. The total molar flow rate and velocity of the reaction mixture at the cross-section of the reactor are as follows: where is the molar flow rate of reaction component , i.e., H2SO4, SO3, H2O, SO2 and O2; is the radial cross section area of the reactor, and is the universal gas constant.

Figure 2

Schematic of tubular plug-flow reactor.

The data of catalyst selection, reactor structure, and thermodynamic parameters of the reaction mixture are determined according to Ref. [1], as listed in Table 1.

Table 1

Parameters of the reference reactor.

Parameter	Symbol	Value
Overall heat transfer coefficient/J/(K·m^2·s)	U	170
Dynamic viscosity/kg/(m·s)	η	4 × 10−5
Catalyst bed porosity	ε	0.45
Catalyst pellet density/kg/m^3	ρ_p	4200
Catalyst pellet diameter/m	D_p	0.003
Inner diameter of reactor/m	D	0.030
Length of reactor/m	L	3.090
Inlet temperature/K	T in	800
Inlet pressure/bar	P in	7.1
Inlet total molar flow rate	F tot,in	0.034
Molar fraction of inlet H2SO4	F_H_2SO_4,in	0.094
Molar fraction of inlet SO3	F_SO_3,in	0.425
Molar fraction of inlet H2O	F_H_2O,in	0.481
Molar fraction of inlet SO2	F_SO_2,in	0.000
Molar fraction of inlet O2	F_O_2,in	0.000

2.2. Models of Kinetics and Thermodynamics

The fluid flow, heat transfer, and chemical reaction of the reaction mixture in a tubular plug-flow reactor follow momentum, energy, and mass conservation equations, respectively, which are given by: where and are the density and flow velocity of the reaction mixture on the entrance section, respectively; subscript = 1, 2 represents the reaction types (I) and (II); is the reaction rate of mass per unit catalyst, and they are and ; and are the component molar constant-pressure heat capacity and the reaction enthalpy of the reaction type , and their expressions are given in the Appendix A.

The heat transfer from the heat source outside the tube to the reaction mixture inside the tube follows Newtonian heat transfer law:

For different reaction conditions and mechanisms, the driving force in the kinetic equation could be written as different mathematical forms, and the corresponding coefficients in the kinetic equation should be determined by experiments and also be different for different choices of the driving force. According to Ref. [1], the condition that the chemical reaction occurred at the vicinity of the equilibrium is assumed to be satisfied, and all components are assumed to have stoichiometric reaction order, so the chemical reaction rates of reaction types (I) and (II) are as follows: where and are the reaction rate constants of reaction types (I) and (II), according to Ref. [1], , ; P represents the partial pressure of the corresponding component; is the equilibrium constant of the chemical reaction type ; is the standard Gibbs free enthalpy of the reaction type , and the expression is given in the Appendix A. The driving force in the kinetic Equation (12) is written as , and effects of the different forms of the driving force on the optimization results will be considered in another paper in the future.

The SO2 yield of the tubular plug-flow reactor is as follows:

The local EGR of the tubular plug-flow reactor is as follows: where subscripts ht, f, and cr represent the local EGRs of heat transfer, fluid flow, and chemical reaction, respectively.

The total EGR is obtained by integrating the local EGR, i.e.,

3. Parameter Analyses of Sulfuric Acid Decomposition Reactor

By changing the inlet parameters of the reference reactor, including the inlet temperature Tin, pressure Pin and the total molar flow rate Ftot,in, the total EGR and the SO2 yield of the reference reactor are analyzed, and the influences of the initial inlet conditions on the two performance objectives can be obtained. The variation ranges of the initial inlet parameters are: 750 K ≤ Tin ≤ 900 K, 4 MPa ≤ Pin ≤ 9.5 MPa, and 0.0027 mol/s ≤ Ftot,in ≤ 0.1 mol/s.

Figure 3 shows the effects of the temperature Tin of the reaction mixture on the total EGR and the SO2 yield. It can be seen that the total EGR decreases nonlinearly with the increase of the temperature Tin, and the decreasing trend is fast firstly and then slow; when the temperature Tin increases from 750 °C to 900 °C, the total EGR decreases from 0.331 W/K to 0.189 W/K, i.e., decreases by 43%. The main reason is that with the temperature Tin of the reaction mixture increases, the heat transfer temperature difference between the reaction mixture and the external heat source decreases, which reduces the local EGR of heat transfer and the total EGR. The SO2 yield increases very slowly with the increase of the temperature Tin, and when the temperature Tin increases from 750 °C to 900 °C, the SO2 yield increases by only 0.4%. It can be seen that the total EGR can be reduced by increasing the temperature Tin of the reaction mixture, i.e., the irreversibility of the sulfuric acid decomposition process could be reduced by increasing the Tin of the reaction mixture. However, it is not significant to increase the SO2 yield by increasing the temperature Tin of the reaction mixture.

Figure 3

The effects of on the total EGR and the SO2 yield.

Figure 4 shows the effects of the pressure Pin of the reaction mixture on the total EGR and the SO2 yield. It can be seen that the curve of the total EGR is concave and parabolic-like with the increase of the pressure Pin, and the minimum value is 0.224 W/K when the pressure Pin is about 0.85 MPa. The SO2 yield decreases linearly with the increase of the pressure Pin. When the pressure Pin increases from 0.4 MPa to 1 MPa, the SO2 yield decreases from 0.0118 mol/s to 0.0105 mol/s, i.e., decreases by 11.02%.

Figure 4

Effects of on the total EGR and the SO2 yield.

Figure 5 shows the effects of the molar flow rate Ftot,in of the reaction mixture on the total EGR and the SO2 yield. It can be seen that the total EGR and the SO2 yield increase with the increase of the molar flow rate Ftot,in, and the minimum total EGR and the maximum SO2 yield are mutually restricted. When the molar flow rate Ftot,in increases from 0.027 mol/s to 0.10 mol/s, the total EGR and the SO2 yield increases by 4.8 times and 1.8 times, respectively.

Figure 5

Effects of on the total EGR and the SO2 yield.

4. Multi-Objective Optimization and Result Analyses

From the analyses in Section 3, when the three inlet parameters are chosen as optimization variables, and the minimum total EGR and the maximum SO2 yield are taken as optimization objectives, respectively, there is no optimal solution to achieve the extremum values of the total EGR and SO2 yield at the same time. Therefore, how to select the appropriate initial inlet conditions to achieve the relative optimal total EGR and SO2 yield is very important. The NSGA-II algorithm is one of the excellent algorithms to solve multi-objective optimization problems, and can give a series of non-inferior solutions (solutions that cannot be optimized for arbitrary objectives without making other objectives worse) of multi-objective problems. The corresponding improvement process is called Pareto improvement, the corresponding set of non-inferior solutions is called the Pareto-optimal solution set, and the corresponding objective function solution is called the Pareto-optimal front.

Figure 6 shows the flow chart of the NSGA-II algorithm. In this section, all of the Tin, Pin and Ftot,in are taken as the optimization variables to minimize the total EGR and maximize the SO2 yield. The optimization intervals of the variables are consistent with the previous single-variable analysis.

Figure 6

Basic flow chart of NSGA-II algorithm.

Figure 7 is Pareto optimal frontier of a reference reactor based on the objective of minimizing total EGR and maximizing SO2 yield, where points A and B represent the solution of the maximum SO2 yield and the minimum total EGR, respectively. At point A, the weighting coefficient of SO2 yield in multi-objective optimization is 1, and the weighting coefficient of total EGR is 0, it is also the solution of maximizing the SO2 yield. Similarly, point B is the solution of minimizing the total EGR. From Figure 7, it can be seen that the minimum total EGR and the maximum SO2 yield are mutually constrained, and they cannot achieve the extremum values at the same time. Only the relative optimal solutions of the two objectives under different weighting coefficients can be found, that is, the non-inferior solution. One can select the appropriate optimal solution from the Pareto-optimal solution set according to different needs to meet the different demands of decision-making purposes. Commonly used multi-objective decision-making methods are Shannon, LINMAP, and TOPSIS, but in the actual decision-making process, decision-making is usually based on actual engineering experience and personal preferences of decision-makers, there is no universal way to make decisions.

Figure 7

Pareto optimal frontiers of reference reactor.

In this paper, in order to facilitate the comparison with the reference reactor, a suitable multi-objective decision point (point C) is selected for comparison. Because the solution of the minimum specific EGRs is the solution of the total EGR and the yield under a certain ratio, the decision point of the minimum specific EGR must be on the Pareto-optimal front, which can be used as an important basis to verify the accuracy of the NSGA-II algorithm results.

Figure 8 is the bar chart of the target value of the reference reactor under optimization and non-optimization. Table 2 lists the results of each optimization target condition. It can be seen that compared with the reference reactor, the SO2 yield of the reactor with the maximum yield increases by 118%, but the total EGR increases by 222%; the total EGR of the minimum EGR reactor decreases by 40%, and the corresponding SO2 yield also decreased by 22%; the total EGR and the SO2 yield of the reactor with the minimum specific EGR decrease by 38% and 16%, respectively. From Figure 7, it can be easily concluded that the reference reactor is not located at the Pareto optimal frontier, so the reference reactor can be optimized by Pareto improvement. A non-inferior solution (point C) is obtained by the multi-objective optimization method, in which the total EGR of the reactor decreases by 9% and the SO2 yield of the reactor increases by 14% compared to the reference reactor. Also, from Figure 7, it can be seen that a series of non-inferior solutions located at the upper left of the decision point (point E) of the reactor have good properties of reducing the total EGR and increasing the SO2 yield.

Figure 8

Comparison of total EGR and the yield of optimized objectives.

Table 2

Calculation results of each target.

	Reactor Inlet Parameters			EGR		SO2 Yield
	Temperature Tin(K)	PressureP_in(1×10^5Pa)	Molar RateF_tot,in(mol·s^−1)	Σ_tot/W·K^−1		ΔF_SO_2/mol·s^−1
Reference reactor	800	7.10	0.034	0.2316	——	0.01100	——
Maximum yield	896	8.97	0.010	0.7450	↑222%	0.02395	↑118%
Minimum EGR	893	8.69	0.027	0.1388	↓40%	0.00862	↓22%
Specific EGR	900	8.62	0.030	0.1446	↓38%	0.00930	↓16%
Multi-objective optimization	894	9.18	0.041	0.2111	↓9%	0.01256	↑14%

Figure 9, Figure 10 and Figure 11 show the distribution of the Tin, Pin and Ftot,in in Pareto-optimal fronts, and the black and white dots in the figures represent the total EGR and the SO2 yield, respectively, which exist in pairs. As seen from Figure 9 and Figure 10, the Tin and Pin of the reaction mixture in Pareto-optimal fronts are mainly distributed in high-temperature (892–896 K) and high-pressure (9.0–9.2 bar) area, so increasing the Tin and Pin of the reaction mixture is an important means for Pareto improvement. Figure 11 shows that the Ftot,in of the reaction mixture in Pareto-optimal fronts distributes uniformly in its optimal range, which indicates that adjusting the Ftot,in of the reaction mixture in Pareto-optimal fronts is an important means to reconcile the contradiction between the minimum total EGR and the maximum SO2 yield.

Figure 9

Distribution of inlet temperature in Pareto-optimal fronts.

Figure 10

Distribution of inlet pressure in Pareto-optimal fronts.

Figure 11

Distribution of total inlet molar flow rate in Pareto-optimal fronts.

5. Conclusions

In this paper, the effects of reaction mixture inlet parameters on the total EGR and SO2 yield of the tubular plug-flow sulfuric acid decomposition reactor are analyzed, and the multi-objective optimization for the two performance objectives are carried out by using FTT. The results show that:

When the Tin increases from 750 °C to 900 °C, the total EGR decreases by 43% and the SO2 yield increases by 0.4%. When the Pin increases from 0.4 MPa to 1 MPa, the curve of the total EGR versus the Pin is a concave parabolic-like, the minimum value of the total EGR is 0.224 W/K when the Pin equals to 0.85 MPa, and the corresponding SO2 yield decreases by 11%. When the Ftot,in increases from 0.027mol/s to 0.10mol/s, the total EGR and the SO2 yield increase by 4.8 times and 1.8 times, respectively.

The reference reactor can be Pareto improvement, one of the non-inferior solutions can reduce the total EGR by 9% and increase the SO2 yield by 14% compared to those of the reference reactor.

FTT is a powerful theoretical tool for the performance analysis and optimization of tubular plug-flow sulfuric acid decomposition reactor. The NSGA-II algorithm is an effective mathematical tool for the multi-objective optimization of tubular plug-flow sulfuric acid decomposition reactor. The Pareto-optimal fronts obtained in this paper has a certain theoretical guiding significance for the optimal designs of the actual sulfuric acid decomposition reactors.

Table A1

Thermodynamic coefficients.

Gas	M_W,i/kg·mol^−1	A_i	B_i	C_i	D_i	E_i	F_i	G_i	T_min/K	T_max/K
SO2	6.40×10^−2	21.430	74.351	−57.752	16.355	0.087	−305.769	254.887	298	1200
O2	3.20×10^−2	29.659	6.137	−1.187	0.096	−0.220	−9.861	237.948	298	6000
SO3	8.01×10^-2	24.025	119.461	−94.387	26.926	−0.118	−407.853	253.51	298	1200
H2O	1.80×10^−2	30.092	6.833	6.793	−2.534	0.082	−250.881	223.397	500	1700
H2SO4	9.81×10^−2	47.289	190.331	−148.123	43.868	−0.740	−758.953	301.296	298	1200