Methanol Production via CO2 Hydrogenation: Sensitivity Analysis and Simulation—Based Optimization

El dióxido de carbono (CO2) es uno de los gases de efecto invernadero, que puede causar el calentamiento global. Uno de los estudios para mitigar las emisiones de CO2 a la atmósfera es convertir el CO2 en productos valiosos (es decir, metanol). Para hacer que la producción de metanol a través de la hidrogenación de CO2 sea un proceso competitivo, se deben considerar las condiciones operativas óptimas con un coste de producción mínimo. Este artículo estudió una aplicación de la metodología de superficie de respuesta en la optimización de la producción de metanol a través de la hidrogenación de CO2. El objetivo de esta optimización fue minimizar el coste de producción de metanol por tonelada producida de metanol. El análisis de sensibilidad se realizó para determinar los parámetros que muestran impactos significativos en el coste de producción de metanol. La metodología de superficie de respuesta junto con el solucionador de programación no lineal se utilizaron como herramienta de optimización. Los resultados mostraron que la metodología de superficie de respuesta se aplicó con éxito a la producción de metanol a través del proceso de hidrogenación de CO2. El costo mínimo de producción de metanol obtenido fue de $ 565.54 por tonelada de metanol producido con las condiciones operativas óptimas de la siguiente manera. Presión de entrada al primer reactor: 57,8 bar, temperatura de entrada al primer reactor: 183,6 ◦C, presión de entrada al segundo reactor: 102,6 bar, temperatura de salida del enfriador de corriente líquida después del segundo reactor: 63,5 ◦C, temperatura de entrada a la primera columna de destilación: 51,8 ◦C.

section_content: Most of energy in the world is currently from combustion of carbonaceous fuels, which are coal, oil, and natural gas. CO 2 emission from this combustion is considered as the second contribution to the greenhouse effect (9-26%), after the water vapor and clouds (36-72%). The recent attempts predict that CO 2 will show stronger greenhouse effect when its amount in the atmosphere is double (Jaworowski et al., 1992). The net CO 2 emissions could increase at around 5. 4% over the next few decades (Radhi, 2009). Due to this concern, many applications and researches on CO 2 conversion and utilization have been studied to control amount of CO 2 releasing to the atmosphere. CO 2 can either be used directly or as feedstock to produce useful chemicals and materials. For the direct use, CO 2 is utilized in many different applications such as food preservation, beverage carbonation, fire extinguisher, supercritical extraction, dry ice, etc (Song et al., 2002). For conversion of CO 2 to other products, urea synthesis is the largest production while the methanol synthesis is the second largest production in this area (Naims, 2016). 

Methanol is commonly used as both solvent and reactant in chemical industry. Uses of methanol is found in many household products, including paints, varnishes, cleaning products (Dasgupta and Klein, 2014). Moreover, methanol can use as a motor fuel or gasoline blending component (Ingamells and Lindquist, 1975). The methanol-fueled vehicles use a blend of 85 percent methanol with 15 percent unleaded gasoline (M85) (Ingamells and Lindquist, 1975; EPA, 2002; Bukhtiyarova et al., 2017). From laboratory and road tests indicate that adding of 10% methanol can raise octane 2-3 numbers. In the research of Gabele (1990), the results showed that increase of methanol content does not affect the emission rate of exhaust gas. 

In the conventional methanol production, methanol is produced from petroleum product (synthesis gas) via hydrogenation of CO and CO 2, and reversed water-gas shift reaction (María et al., 2013). The commercial methanol productions from synthesis gas use CuO/ZnO/Al 2 O 3 as catalyst (Sayah et al., 2010; Jadhav et al., 2014). The commercial CuO/ZnO/Al 2 O 3 coupled with a zeolite membrane reactor can provide higher CO 2 conversion, methanol yield, and selectivity compared with a traditional reactor (Gallucci et al., 2004). 

Copper based catalyst is mostly used for CO 2 hydrogenation to methanol due to its cheap and higher catalyst activity (Ali et al., 2015). Saito and Murata studied of Al 2 O 3 supported Cu-based catalyst with Co-precipitation technique. The catalyst shows high activity and the stability was improved by adding colloidal silica (Saito and Murata, 2004). 

Several developments of catalysts used in catalytic conversion to improve the performance of catalyst were studied. Mg and Mn was promoted on CuZnZr catalyst. Adding of MgO and MnO lead to increase in catalytic activity of the catalyst (Sloczyński et al., 2003). Doping of Mn to Cu/Zn/Zr catalyst can increase the methanol production rate. The zirconium indicates the advantageous influence on the catalyst activity (Lachowska and Skrzypek, 2004). The Pd/ZnO catalysts over multi-walled carbon nanotubes was found to have the turnover frequency reached 0. 015 per second (Liang et al., 2009). 

Recently, the production of methanol from direct CO 2 hydrogenation is of interest (Samimi et al., 2017; Marlin et al., 2018). The process has potential to mitigate the CO 2 emission to the atmosphere. The methanol production from direct CO 2 (using pure sources of CO 2 and H 2 ) has several advantages over the conventional process-it results in significantly less byproducts, and requires less energy in product purification (Marlin et al., 2018). However, the methanol production cost via direct CO 2 hydrogenation is 2-2. 5 times higher than the cost of conventional process (Atsonics et al., 2015). The process of methanol production via CO 2 hydrogenation consumes more utilities than the conventional methanol production (Machado et al., 2014). 

There are some studies on simulation-based optimization of methanol synthesis to increase the methanol production rate (Hoseiny et al., 2016; Leonzio, 2017). For methanol production from synthesis gas, the effects of changes in operating conditions on the production rate were studied to maximize the production rate. The studied operating parameters were feed flow rate, pressure and temperature of feed, and the cooling water temperature. The results showed that the methanol production rate can be increased by 7% with higher feed pressure and lower feed temperature (Hoseiny et al., 2016). For methanol synthesis via direct CO 2 hydrogenation, Grazia Leonzio developed mathematical model of the reactor used in methanol production. The impacts of reaction temperature, reaction pressure, H 2 /CO 2 ratio, and the recycle factor on methanol production rate and reactor volume were studied (Leonzio, 2017). 

However, optimization of methanol production via CO 2 hydrogenation, which considers all possible operating parameters, for the minimum production cost has not been studied. This paper studied an application of response surface methodology (RSM) in optimization of methanol production via CO 2 hydrogenation process. To be able to apply the RSM in optimization, the response surface has to be in the form of the second order model. The sensitivity analysis was used to determine the significant operating parameters. The RSM coupled with non-linear solver was employed to obtain the local optimal operating conditions and the minimum methanol production cost. 

This paper is organized as follows: section Process simulation and economic evaluation gives details of process simulation and economic evaluation. Section Methodology provides methodology used in this work including details of sensitivity analysis and simulation-optimization framework. Section Results and discussion discusses results of process optimization. Section Conclusion provides conclusion.
section_name: INTRODUCTION

section_content: In the study of methanol production process via CO 2 hydrogenation, the process simulation was combined with the economic analysis to evaluate the objective function (the methanol production cost) corresponding to the decision variables. In what follows, details of process simulation and economic evaluation are explained.
section_name: PROCESS SIMULATION AND ECONOMIC EVALUATION

section_content: The methanol production process via CO 2 hydrogenation was simulated using Aspen Hysys version 8. 8 process simulator. Peng-Robinson was used as the thermodynamic property package. The physical properties were predicted using the thermodynamics based equation as shown in Equations ( 1 ) and (2). Where H ID is the Ideal Gas Enthalpy and S ID is the Ideal Gas Entropy. 

In the process simulation, two reactors were employed due to low conversion of the CO 2 hydrogenation reaction. represents the process flow diagram of methanol production via CO 2 hydrogenation (Wiesberg et al., 2016). In this process, the feed of 1,000 kmoles per hour of carbon dioxide at 40 • C and 20 bar was mixed with the 3,000 kmoles per hour of hydrogen (at the same conditions). The mixture was then compressed, heated and sent to the first equilibrium reactor. 

The first reactor partially converted CO 2 to methanol as liquid product, as shown in Equations ( 3 )-( 5 ) (Tidona et al., 2013). The unreacted CO 2 and H 2 then entered the second equilibrium reactor to produce more methanol product. The pressure of gas phase leaving the second reactor was reduced to recover methanol as liquid phase. All liquid methanol products were sent to the first distillation column, where the light components (CO, CO 2, and H 2 ) leaved at the top of the column. The mixture of methanol and water leaved the column at the bottom and entered the second distillation column, where the methanol product with purity of 99. 5%mole was obtained at the top of the column. 

For equipment specifications, both reactors were simulated using equilibrium model. The sets of reaction used in each reactors are Equations ( 3 )-( 5 ). The efficiency of the pump was assumed at 75% (adiabatic). The efficiencies of all compressors were assumed at 75% (adiabatic). The specifications of column T-100 were condenser temperature at 40. 1 • C and reflux ratio at 0. 5. The specifications of column T-101 were 99. 5%mole methanol product purity, and reboiler temperature at 143. 2 • C. The detailed process simulation of methanol production via CO 2 hydrogenation is provided in Appendix A.
section_name: Process Simulation

section_content: In an economic analysis, capital, and operating costs were included in calculation of methanol production cost. The capital cost involved all major equipment, except pump, and piping. The capital cost was calculated using equations and data from the capital equipment-costing program (Turton et al., 2003). The data was adjusted for inflation from year 2001-2017 by using values of the Chemical Engineering Plant Cost Index, CEPCI. The CEPCI value in 2001 is 297, and the CEPCI value in 2017 is 541. 7 (Jenkins, 2018). 

The equipment costs were estimated base on total module costs (C TM ), shown in Equation ( 6 ), where n represents the total number of pieces of equipment. C BM is the bare module cost, which can be estimated from Equation ( 7 )

Where C P o is the purchased cost for base conditions, which can be determined from Equation (8) F P is the pressure factor F M is the material factor

In Equation ( 8 ), parameters A is the capacity or size parameter for the equipment, K 1, K 2, K 3 are the maximum and minimum values used in the correlation. The description of major equipment, and K 1, K 2, and K 3 used in this studied are shown in Table 1. 

For cost modeling of reactors, the dimension of reactors were fixed at height of 5. 8674 m and diameter of 1. 068 m. The values of K 1 is 3. 4974, K 2 is 0. 4485, K 3 is 0. 1074, B 1 is 2. 25, B 2 is 1. 82, and F M is 3. 1. The pressure factors (Fp) for the reactors were determine using Equation ( 9 ), where P is in barg. 

Assumption used in the economic analysis are listed below. 

-The plant operates for 8,400 h per year. 

-The plant is expected to have a 10-year plant life with on salvage value. -The labor cost is $40,000 per operator per year. 

-The working capital is 15% of fixed capital investment. 

-Total capital investment including fixed capital investment and working capital. -The maintenance and repairs are 5% of fixed capital investment. -The operating supplies are 15% of maintenance and repairs. 

-The laboratory charge is 15% of labor. 

-The administrative expenses is 50% of labor cost. 

-The maintenance and repairs, the laboratory charge and the administrative expenses increase annually 3%. -The local taxes and insurance are 4% of fixed capital investment. -The plant overhead cost is 60% of labor. 

-The cost of carbon dioxide is $12. 10 per ton (Wiesberg et al., 2016). -The cost of hydrogen is $1,250 per ton (Wiesberg et al., 2016) -The cost of process water at 25 • C is $0. 0259 per cubic meter. 

-The cost of saturated steam at 6. 9 bar is $55 per ton. 

-The cost of electricity is $0. 127 per kW h.
section_name: Economic Evaluation

section_content: The methodology in this work consists of two sections. The first section is "Sensitivity analysis" to determine the parameters that show significant impacts on the methanol production cost. The second section is "Simulation-Optimization" to optimize the significant operating parameters of methanol production process for the minimum methanol production cost. Details of Sensitivity analysis and Simulation-Optimization are described in the following sections.
section_name: METHODOLOGY

section_content: Sensitivity analysis is an efficient tool, which can be used to determine the significant parameters. In this section, the impacts of eight parameters on the methanol production cost were determined using sensitivity analysis. An increase in pressure and decrease in temperature leads to increase of the equilibrium conversion of CO 2 to methanol (Witoon et al., 2015). The eight studied parameters are (1 and 2) inlet pressure and temperature to the first reactor, ERV-100, for CO 2 and H 2 conversion conditions, (3) inlet temperature to the first separator for methanol product separation condition, (4) inlet pressure to the second reactor, ERV-101, for CO 2 and H 2 conversion conditions, (5 and 6) outlet temperature of both coolers, located after the second reactor, for methanol product separation condition from gas and liquid mixture, (7) outlet pressure at the valve VLV-102 for recovering methanol from gas phase mixture, (8) Inlet temperature to the first distillation column for separation between liquid methanol and water products, and gaseous unreacted reactant separation. For all parameters, except inlet temperature to the first reactor, we performed the sensitivity analysis by varying value of interested parameter within the range ± 25% of its base value, and fixing other parameters at their base values. For inlet temperature to the first reactor parameter, we varied the value of temperature within the range ± 10% of its base value. Table 2 shows the base value used in sensitivity analysis of each parameter.
section_name: Sensitivity Analysis

section_content: The results from sensitivity analysis show that there are five parameters that show significant impacts on methanol production cost. Figure 2 shows inversely proportional relationship between operating parameters and methanol production cost (the methanol production cost decreases with an increase in the value of parameter). In the figure, the inlet pressure to the first reactor and the inlet temperature to the first distillation column show significant inversely proportional relationship with the methanol production cost. 

On the other hand, Figure 3 shows proportional relationship between operating parameters and methanol production cost. In the figure, the inlet temperature to the first reactor, the inlet pressure to the second reactor, and the outlet temperature of the liquid stream cooler after the second reactor show significant proportional relationship to the methanol production cost. 

Five significant variables obtained from this sensitivity analysis will be used as decision variables to the optimization problem.
section_name: Sensitivity Analysis: Results

section_content: In this section, RSM (Nuchitprasittichai and Cremaschi, 2011) coupled with non-linear solver were employed to optimize the methanol production via CO 2 hydrogenation process. The objective is to minimize the methanol production cost per tons produced methanol ($/tons produced methanol). Response surface methodology is a statistical tool used to represent the relationship between independent variables and dependent variable(s). In this work, RSM was employed to represent the relationship between significant operating parameters and the methanol production cost per tons produced methanol. The significant operating parameters (decision variables) are the significant parameters obtained from sensitivity analysis (in section Sensitivity analysis) which are x 1 : inlet pressure to the 1st reactor, x 2 : inlet temperature to the 1st reactor, x 3 : inlet pressure to the 2nd reactor, x 4 : outlet temperature of the liquid stream cooler after the 2nd reactor, x 5 : inlet temperature to the 1st distillation column. 

Figure 4 represents simulation-optimization algorithm. In the algorithm, the range of each decision variables was first determined. In this work, the ranges of decision variables were obtained from literature (Shen et al., 2000; Witoon et al., 2015). Factorial design with two levels (2 k Factorial design, section 2 k Factorial design) and a center point were used as a data set of the decision variables. Then, process simulation, coupled with cost analysis, was run corresponding to operating conditions in the data set to obtained methanol production costs. The data of operating conditions and the corresponding methanol production costs were combined and fit the first order model by regression analysis (using Design Expert 11 software). If the data fit the first order model (R-squared > 0. 7), it means that the region of operating conditions and methanol production costs is in linear region. The steepest descent (section Steepest descent) was performed to move operating conditions to the region of lower methanol production cost. On the other hand, if the data do not fit the first order model, it means that the data are in nonlinear region. The Box-Behnken Design [BBD, section Box-Behnken design (BBD)] was performed to collect more data in non-linear region. All data were then combined to fit the second order (non-linear) model. The non-linear model was then used as the objective function in optimization problem. Microsoft Excel (non-linear programing, NLP) solver was used to solve the non-linear optimization problem (section Optimization formulation) for the minimum methanol production cost. The optimal solution obtained from this RSM is considered as a local optimal solution.
section_name: Simulation-Optimization

section_content: The 2 k factorial design is a design for k factors. Each factor is at only two levels, which are values at the lower and upper bounds of the parameter range. The design is generated in coded variables. The coded variable -1 represents the lower bound value, and the coded variable 1 represents the upper bound value. The design is widely used in factor screening experiments. The design observes at each corner of the cube (corner points). Figure 5A demonstrates 2 k factorial design for three independent factors.
section_name: k Factorial Design

section_content: The steepest descent was performed when the data obtained from 2 k factorial design (including one center point) fits the first order model. The steepest descent moves the independent variables to the direction of maximum decrease in the response (Douglas, 2013). To perform the steepest descent, the independent variables will be coded to the (-1, 1) interval as shown in Equation (10). Where x i is the coded variable, and ξ i is the natural variable. 

The regression analysis was performed to fit the data set of 2 k factorial design and a center point with the first order model. Equation ( 11 ) represents example of the first order model with five independent factors

A step size ( x i ) in each independent variable was determined by Equation ( 12 ). Where b i, highest is the highest value among all coefficients b i, except b o. 

The path of steepest descent for x i was determined by using Equation ( 13 ). Where j is point along the path. 

x i condition along the path = x i, mid point + j x i 1 2 Range of x i ;

Experiments were conducted along the path of steepest descent until no further decrease in response.
section_name: Steepest Descent

section_content: The Box-Behnken design is the three-level design for fitting response surface (the second order model). The design consists of the midpoint of each edge of the space, and a center point. Figure 5B demonstrates Box-Benhnken design for three independent factors.
section_name: Box-Behnken Design (BBD)

section_content: This section gives detailed optimization formulation. In the optimization formulation, a second order model was used to represent the response surface. The form of a full second-order model is shown in Equation ( 14 ). All five decision variables studied in this work are continuous variables. The decision variables were coded in the ranges between -1 and 1. Therefore, all the constraints are in the boundary of coded variables. 

Where y is the methanol production cost ($ per tons produced methanol) Decision variables:x 1 : inlet pressure to the 1st reactor x 2 : inlet temperature to the 1st reactor x 3 : inlet pressure to the 2nd reactor x 4 : outlet temperature of the liquid stream cooler after the 2nd reactor

x 5 : inlet temperature to the 1st distillation column Subject to : -1 ≤ x i ≤ 1; i = 1, 2, 3, 4, 5
section_name: Optimization Formulation

section_content: We analyzed the impacts of eight operating parameters on the methanol production cost. Then, optimized the value of significant parameters for the minimum methanol production
section_name: RESULTS AND DISCUSSION

section_content: Data set 1 Data set 2 cost. The results obtained from the sensitivity analysis, the steepest descent, and the optimization are discussed in this section.
section_name: Decision variables Ranges

section_content: In performing simulation-optimization algorithm, the first data set (data set 1) from 2 k factorial design with initial range of each decision variable (shown in Table 3, Data set 1 column) was fitted with the first order regression model. The obtained R-squared value is 0. 9995. This means the data still fit the linear model. Therefore, we performed steepest descent to search for the region of lower methanol production cost. Figure 6 shows steepest descent in moving the operating conditions to lower methanol production cost region. 

From the steepest descent results, for the first reactor, an increase in inlet pressure from 56. 0 to 58. 0 bar, and a decrease in inlet temperature from 192. 0 to 183. 0 • C yielded an increase in methanol product at equilibrium condition (%conversion increased from 40 to 45%). For the second reactor, a decrease in inlet pressure from 118. 0 to 102. 0 bar yielded an increase in methanol product in the reactor. The outlet temperature of the liquid stream cooler after the 2nd reactor decreased from 74 • C to around 65 • C yielded an increase in liquid methanol separated from the separator. The inlet temperature to the first distillation column increased from 48. 0 to 52. 0 • C yielded higher efficiency to separate most of remaining CO 2 and lighter components from methanol and water. 

The steepest descent moves the operating conditions to the new conditions as follows: inlet pressure to the first reactor: 58. 84 bar, inlet temperature to the first reactor: 184. 0 • C, inlet pressure to the second reactor: 106. 64 bar, outlet temperature of the liquid stream cooler after the second reactor: 65. 54 • C, and inlet temperature to the first distillation column 52. 75 • C. We then used these new operating conditions as the middle value of the range of decision variables. The ranges of the new operating conditions are shown in Table 3, Data set 2 column. 

The new data set (Data set 2) was then constructed with 2 k factorial design. The corresponding methanol production costs were collected, as shown in Figure 6. Most of the data in data set 2 shows lower methanol production cost than data set 1, except the data at low inlet pressure to the first reactor. The results agree with the results from sensitivity analysis that the inlet pressure to the first reactor has inversely proportional relationship with the methanol production cost. 

By fitting data set 2 with the first-order regression model, the data did not appropriately fit the model as Rsquared value is 0. 6353 (<0. 7). This means that there is high possibility to find the non-linear region at this step. We then used theses ranges of operating condition with BBD in collecting more data for constructing the nonlinear model.
section_name: Simulation-Optimization: Steepest Descent

section_content: Optimal values

x 1 : Inlet pressure to the 1st reactor (bar) 57. 8

x 2 : Inlet temperature to the 1st reactor ( • C) 183. 6

x 3 : Inlet pressure to the 2nd reactor (bar) 102. 6

x 4 : Outlet temperature of the liquid stream cooler after the 2nd reactor ( • C) 63. 5

x 5 : Inlet temperature to the 1st distillation column ( • C) 51. 8

Predicted methanol production cost ($/ton) 559. 59

Actual methanol production cost ($/ton) 565. 54

after the second reactor (x 4 ), and inlet temperature to the first distillation column (x 5 ), do not show significant impacts on the production cost in the non-linear model (Equation 15 ), the values of these three parameters were set at their lower bounds. 

The prediction accuracy of the model was determined by the percent error. The percent error of the model was estimated by comparing between the predicted methanol production cost with the actual production cost. We obtained the actual production cost by running the process simulation with the corresponding optimal operating conditions. The percent error of this non-linear model is 1. 05. For regression analysis of this model, the R-squared value is 0. 8705, and the adjusted Rsquared value is 0. 8610. Therefore, with the small number of percent error, it can be concluded that this methanol production via CO 2 hydrogenation process can be represented by the second order model. The process is successfully optimized using RSM. 

The optimization formulation was solved in the Intel Pentium 4 processor, core i3. The local optimal solutions were obtained in 2 s. The obtained minimum production cost ($ per tons produced methanol) is an offset between the methanol production cost and the amount of produced methanol. The methanol production cost includes energy consumption, utilities, and capital cost costs. For the optimal solution, the amount of produced methanol increased so that the minimum production cost per tons of produced methanol was obtained. The optimal methanol production rate is 964 kmoles per hours, which produced from 1,000 kmoles per hour of feed carbon dioxide and 3,000 kmoles per hour of feed hydrogen.
section_name: Parameters

section_name: ACKNOWLEDGMENTS

section_content: The financial support from Suranaree University of Technology is greatly acknowledged.
section_name: ACKNOWLEDGMENTS

section_name: DATA AVAILABILITY

section_content: All datasets generated and analyzed for this study are included in the manuscript/Supplementary Files.
section_name: DATA AVAILABILITY

section_name: The Optimal Conditions

section_content: The reduced non-linear model (the model with only significant terms) was constructed with 73 sample points (32 sample points from 2 k factorial design, and 41 sample points from BBD) as shown in Equation ( 15 ). The model represents the relationship between the operating conditions and the methanol production cost. 

Where y is the methanol production cost per tons produced methanol ($ per ton produced methanol), x 1 is inlet pressure to the first reactor, and x 2 is inlet temperature to the first reactor. Please be noted that the obtained values of x 1 and x 2 in Equation ( 15 ) are in code variables (-1 to +1). 

The obtained results have to be converged to the actual operating condition values by using the ranges in Table 3, data set 2 column and Equation ( 10 ), where -1 represents the lowest value and +1 represents the highest values in the range. 

Table 4 shows the local optimal operating conditions (x i ) of methanol production via CO 2 hydrogenation with the minimum production cost per tons produced methanol (y). Since three parameters, which are inlet pressure to the 2nd reactor (x 3 ), outlet temperature of the liquid stream cooler
section_name: The Optimal Conditions

section_content: In the optimization of methanol production via CO 2 hydrogenation, the sensitivity analysis coupled with RSM was successfully represent the relationship between methanol production operating conditions and the methanol production cost per tons produced methanol ($ per ton produced methanol). The non-linear solver was applied to optimize the non-linear relationship model for the minimum methanol production cost. Two reactors were employed in the process. In the optimal region of study, the inlet pressure and temperature to the first reactor show significant impacts on the methanol production cost. The optimal operating conditions of inlet pressure to the first reactor is 57. 8 bar, inlet temperature to the first reactor is 183. 6 • C, and the other three insignificant parameters were set at their lower bounds. The obtained minimum methanol production cost is $565. 54 per ton produced methanol.
section_name: CONCLUSION

section_content: PB simulated the process simulation, collected data, analyzed results, and wrote the first draft of the manuscript. AN designed the methodology of the work, and revised the manuscript for the final version.
section_name: AUTHOR CONTRIBUTIONS

section_content: Suranaree University of Technology funding on Aspen Hysys simulator.
section_name: FUNDING

section_content: The Supplementary Material for this article can be found online at: https://www. frontiersin. org/articles/10. 3389/fenrg. 2019. 00081/full#supplementary-material
section_name: SUPPLEMENTARY MATERIAL

section_content: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
section_name: Conflict of Interest Statement: