The Inﬂuence of Variable CO 2 Emission Tax Rate on Flexible Chemical Process Synthesis

: The emission of greenhouse gasses is a major environmental problem, and efforts are being made worldwide in various ways to encourage producers to reduce their emissions. There is a need to incorporate environmental measures into process design and synthesis, as pollution prevention is a higher priority than waste management, and in this way, more sustainable solutions can also be achieved. One possibility is to introduce a CO 2 tax, the value of which is very uncertain in the future. This paper demonstrates how the CO 2 tax affects the optimal results of synthesizing chemical processes using mixed-integer nonlinear programming (MINLP). It was found that the tax increase promotes the use of better-quality raw materials and more efﬁcient process units. Energy consumption and emissions are reduced and economic performance deteriorates. A multi-period, two-stage stochastic approach with recourse is suitable to incorporate the uncertainty of the CO 2 tax in the MINLP process synthesis and gives better results than a simpler deterministic approach. In the case of the heat exchanger network synthesis, the costs obtained with the stochastic approach were 5% lower, and the emissions 7% lower than with a deterministic approach.


Introduction
One of the most troubling environmental threats is CO 2 emissions [1] and the related problem of climate change [2]. Countries around the world are struggling to reduce CO 2 emissions, which account for a large part of greenhouse gas (GHG) emissions [3]. According to the IEA, emissions in 2020 decreased by 4.8%, compared to the previous year, due to the COVID-19 pandemic [4]. Europe was the fourth-largest producer of GHG emissions in 2019 [5]. It was preceded by China, the USA, and India.
The European Union Emissions Trading System was introduced in 2005 to promote the reduction in greenhouse gas emissions in a cost-effective and economically efficient manner [6]. Taxes on greenhouse gas emissions have been shown to be beneficial in promoting the reduction of pollutant emissions [7]. With appropriate environmental instruments, governments can encourage investment in carbon capture, utilization, and storage technologies [8]. Imposing a carbon tax on coal and oil would increase the cost of electricity compared to cleaner electricity generation [9]. A tax policy that charges higher emissions encourages generators to reduce emissions [10] by adopting more modern technologies [11]. On the other hand, the introduction of energy efficiency measures in countries with high potential for efficiency improvements reduces emissions more than a carbon tax [12]. Compliance with GHG emission tax leads to sustainable processes and changes the way energy is supplied and used [13]. Energy consumption and emissions can be further reduced by integrating heat and power generation in the process [14].
Fossil fuels, which are used to generate heat and electricity and contribute most to global CO 2 emissions [15], are expected to be replaced in the future by renewable energy

Methods
In this work, a mixed-integer nonlinear programming (MINLP) synthesis approach was used for process synthesis under variable CO 2 emission tax rate. A general mathematical model (1) [31] allowed simultaneous optimization of process unit selection (topology), process unit sizes, and operating and control variables. In this problem, there was also an uncertain parameter-namely, the CO 2 emission tax. The variables in the model (1) were divided into the following two groups: • The first-stage variables (topology and size of the process) were the same for all CO 2 taxes and were determined in advance; • The second-stage variables (operating and control variables) were adjusted later when the values of the CO 2 tax were known.
where Z objective variable; y first-stage binary variables for selection of topological process alternatives; d first-stage design variables for sizes of process units; x second-stage control and operating variables; A matrix of constants; a vector of constants; B matrices of constants; c T fixed cost; f variable cost function; g inequality constraints; h equality constraints; θ CO 2 uncertain parameter-CO 2 tax.
A characteristic of this model is that a parameter θ CO 2 , which is a real number, can take infinitely many values between specified lower and upper bounds. Therefore, this problem cannot be solved directly but requires discretization of the uncertain parameter values. This means selecting a finite number of values for the CO 2 emission tax, estimating their probabilities, and solving the corresponding mathematical problem for multiple scenarios. Four approaches to incorporate CO 2 emission tax in process synthesis were established and tested on two chemical engineering case studies: (i) deterministic approach, (ii) modified deterministic approach, (iii) stochastic approach considering different CO 2 taxes in one period, and (iv) stochastic approach considering different taxes on CO 2 emissions over several periods.

Deterministic Approach
In the deterministic approach, the first-stage variables were determined using the expected value of the CO 2 tax, as shown in the one-scenario problem (2): where θ CO 2 the expected value of CO 2 emission tax (EUR/t CO 2 equivalent). The obtained optimal values of the first-stage variables were then fixed, and the secondstage variables were determined by solving the one-scenario problem (3) at different values of the CO 2 tax.
where y opt optimal values of binary variables determined with problem (2); d opt optimal values of design variables determined with problem (2); t index of scenarios, i.e., different values of CO 2 tax rates; T set of CO 2 tax values; T = {t; t = 1, 2, . . . , N t }; Z t objective variable at scenario t; N t the number of CO 2 tax values selected.
The values of the objective variable obtained at different CO 2 tax rates (Z t ) are multiplied by their probabilities, and the deterministic expected value of the objective function is obtained, as shown in Equation (4).
where E(Z deter ) deterministic expected value of an objective variable; pr t probability of emission tax value t.

Deterministic Approach with Possible Correction of the First-Stage Variables
The approach described in the previous section is conservative in the sense that it does not allow for any adjustment of the first-stage variables with respect to the actual realization of the CO 2 tax value in the future. In practice, however, some changes and process modifications are usually possible, which are designed as process reconstruction or retrofit. A modified approach was, therefore, proposed in which the values of the first-stage variables are determined at the initial value of the CO 2 tax using problem (2), and then the values of the design variables at the specific CO 2 tax rate can be modified (∆d t ) using problem (5). In this problem, the design variables of those process units that were determined in the previous step with problem (2) can be increased. New units can also be added, which are determined by optimizing the binary variables in (5). The increases in the design variables are limited by the upper bound since in practice retrofitting is often only possible within the limits of certain physical dimensions.
where d d values of design variables obtained at the initial (expected) value of CO 2 tax; ∆d t increase of design variables at a given value of CO 2 tax; U upper bound for the changes in design variables.

A Stochastic Approach with Different CO 2 Emission Taxes within a Single Time Period
Deterministic approaches consider different emission tax rates sequentially, which usually does not allow to establish optimal synergies between first-and second-stage variables, and the obtained results may not be optimal in the long run. Therefore, a stochastic approach was introduced in which the scenarios are considered simultaneously during the optimization. The hypothesis was that the simultaneous approach would provide better optimal results by establishing an appropriate trade-off between the firstand second-stage variables. The stochastic MINLP synthesis is formulated as a two-stage multi-scenario problem with recourse (6), where the index t represents the CO 2 emission tax rates that determine several possible scenarios within the same time period. Tax rates (scenarios) can be determined using, e.g., Gaussian quadrature points.
The first-stage variables y and d are the same for all CO 2 tax values and are therefore not indexed over the set of scenarios. The second-stage variables are adjusted for each CO 2 tax value and are therefore indexed over the set t∈T. The expected value of the objective variable is obtained simultaneously at different tax values by multiplying the part of the objective function f by their probabilities, pr t . If the CO 2 tax is a single uncertain parameter, the number of scenarios would be moderate, and the multi-scenario problem (6) can usually be solved simultaneously with manageable computer effort.

A Stochastic Approach with Varying CO 2 Emission Taxes over Multiple Time Periods
Tax rates on CO 2 emissions are expected to be in phases, increasing over the long run. The goal was, therefore, to synthesize an optimal process under variable CO 2 tax rates over several time periods in which different distributions of CO 2 tax rates are expected. To accomplish this task, a multi-period multi-scenario two-stage stochastic problem with recourse was formulated (7).
where p index of time periods; P set of time periods; P = {p; p = 1, 2, . . . , N p }; N p number of time periods.
The synthesis of the process with the problem (7) is performed simultaneously for different values of the CO 2 emission tax in several time periods. The second-stage variables are indexed by the index of CO 2 tax rates, t, and by the time periods, p. The size of problem (7), measured by the number of conditional constraints and variables, is about t × p times that of problem (2).

Results
The objectives of the study were i) to investigate how the tax on CO 2 emissions affects the MINLP synthesis of processes and ii) to develop a methodology for the synthesis of processes that takes into account the uncertainty in the values of this tax. Case studies of a chemical production process and a heat exchanger network synthesis were used to illustrate different approaches. Both examples presented in this paper represent an extended version of the authors' earlier work [32].

Influence of CO 2 Emission Tax on the Process Synthesis
The process superstructure for the production of chemical B is shown in Figure 1. There were two feeds to choose from: the more expensive Feed 1 contained 80% of reactant A and 20% of inert component C, while Feed 2 was cheaper and contained 70% of reactant A and 30% of inert component C. The selected feed was mixed with the recycled stream and heated in unit C1 before entering the reactor. Two alternative reactors were available in which reactant A was converted to product B. Reactor RCT1 was cheaper since the conversion of reactant A was 10%, while reactor RCT2 was more expensive, with a 20% conversion. The effluent from the reactor was cooled in unit H1 and sent to a separator (FLASH). Product B was obtained in the liquid phase, and unreacted reactant A and inert component C in the gaseous phase. Part of this stream was recycled, while the other part was purged to prevent the accumulation of the inert component in the process. The proportion of the purge stream was optimized. A process synthesis for profit maximization under different rates of CO2 emission tax was performed. Two variants were investigated: (i) without heat integration, in which all heating and cooling demand was met by utilities and (ii) with heat integration of the process, in which maximum heat exchange between process streams was sought [21]. For comparison, a synthesis with minimization of GHG emissions was also performed.
The results of the synthesis of the thermally nonintegrated process are presented in Table 1. In terms of topological choices, the more expensive feedstock and reactor were chosen for all CO2 tax values.
In the case of profit maximization, the investment, the overall conversion, the recycle stream, and the utility consumption increased as the CO2 emission tax increased, while the purge stream, the total GHG emissions, and the profit decreased. In the case of GHG emissions, a trade-off was established between the consumption of hot utility and the amount of GHG component in the purge stream. The impact of the latter was greater in this case, as the purge stream was reduced when the CO2 emission tax was increased, while the consumption of hot utility increased.
At a very high CO2 tax rate, the total GHG emission approached the value obtained by minimizing emissions. The profit became negative when the value of the CO2 emission tax exceeded 14 EUR/t.  A greenhouse effect was assumed for the chemicals in the feed stream, the purge stream, and the product. Heaters C1, C2, and C3 use a hot utility, which also contributes to the total greenhouse gas (GHG) emission of the process. The process synthesis was modeled as an MINLP problem. The profit was maximized, which is the difference between the revenue from product sales and the cost of raw materials, utilities, depreciation, and CO 2 emission tax.
A process synthesis for profit maximization under different rates of CO 2 emission tax was performed. Two variants were investigated: (i) without heat integration, in which all heating and cooling demand was met by utilities and (ii) with heat integration of the process, in which maximum heat exchange between process streams was sought [21]. For comparison, a synthesis with minimization of GHG emissions was also performed.
The results of the synthesis of the thermally nonintegrated process are presented in Table 1. In terms of topological choices, the more expensive feedstock and reactor were chosen for all CO 2 tax values. In the case of profit maximization, the investment, the overall conversion, the recycle stream, and the utility consumption increased as the CO 2 emission tax increased, while the purge stream, the total GHG emissions, and the profit decreased. In the case of GHG emissions, a trade-off was established between the consumption of hot utility and the amount of GHG component in the purge stream. The impact of the latter was greater in this case, as the purge stream was reduced when the CO 2 emission tax was increased, while the consumption of hot utility increased.
At a very high CO 2 tax rate, the total GHG emission approached the value obtained by minimizing emissions. The profit became negative when the value of the CO 2 emission tax exceeded 14 EUR/t.
In the synthesis of the heat-integrated process (Table 2), an integrated model of Duran and Grossmann [33] was used. At a low value of CO 2 emission tax (0 and 10 EUR/t), a cheaper Feed 2 with a lower amount of reactant A and a more expensive reactor (RCT2) with a higher conversion were chosen. As the emission tax increased, a more expensive Feed 1 was selected. When minimizing the total GHG emissions, the highest investment and conversion were obtained, while the lowest consumption of hot utility and raw material, as well as the total GHG emissions, were achieved. The obtained results are shown graphically in Figures 2-6. Figure 2 shows the total CO 2 emissions under different CO 2 tax rates. The heat-integrated process has significantly lower CO 2 emissions than the nonintegrated process. For the integrated process, the GHG emissions are high only at low CO 2 tax values because the cheaper raw material with lower content of reactant A was selected. At a CO 2 tax value higher than 7 EUR/t, GHG emissions are significantly reduced. The obtained results are shown graphically in Figures 2-6. Figure 2 shows the total CO2 emissions under different CO2 tax rates. The heat-integrated process has significantly lower CO2 emissions than the nonintegrated process. For the integrated process, the GHG emissions are high only at low CO2 tax values because the cheaper raw material with lower content of reactant A was selected. At a CO2 tax value higher than 7 EUR/t, GHG emissions are significantly reduced. The percentage of the stream that is purged at different CO2 emission tax rates is shown in Figure 3. A larger fraction of the purge occurs in the nonintegrated process. The difference between the purge fraction in the nonintegrated and integrated process at higher CO2 tax rates is about 5 percentage points. At lower CO2 tax rates, the difference is smaller due to the higher emissions caused by choosing a cheaper raw material.  The total conversion of reactant A increases with the CO2 emission tax, as shown in Figure 4. The conversion in the process with heat integration is 13% higher than in the process without heat integration. The flow rate of the feed stream at different CO2 tax rates is shown in Figure 5. The heat-integrated process requires less raw material than the nonintegrated process, except at very low CO2 tax rates, when a cheaper feedstock is chosen in the integrated process. Higher CO2 tax rates encourage lower feedstock consumption and more efficient reactant conversion. The total conversion of reactant A increases with the CO2 emission tax, as shown in Figure 4. The conversion in the process with heat integration is 13% higher than in the process without heat integration. The flow rate of the feed stream at different CO2 tax rates is shown in Figure 5. The heat-integrated process requires less raw material than the nonintegrated process, except at very low CO2 tax rates, when a cheaper feedstock is chosen in the integrated process. Higher CO2 tax rates encourage lower feedstock consumption and more efficient reactant conversion. The percentage of the stream that is purged at different CO 2 emission tax rates is shown in Figure 3. A larger fraction of the purge occurs in the nonintegrated process. The difference between the purge fraction in the nonintegrated and integrated process at higher CO 2 tax rates is about 5 percentage points. At lower CO 2 tax rates, the difference is smaller due to the higher emissions caused by choosing a cheaper raw material.
The total conversion of reactant A increases with the CO 2 emission tax, as shown in Figure 4. The conversion in the process with heat integration is 13% higher than in the process without heat integration. The flow rate of the feed stream at different CO2 tax rates is shown in Figure 5. The heat-integrated process requires less raw material than the nonintegrated process, except at very low CO2 tax rates, when a cheaper feedstock is chosen in the integrated process. Higher CO2 tax rates encourage lower feedstock consumption and more efficient reactant conversion.  Processes 2021, 9, 1720 10 of 16 Figure 6 shows the investment under different CO2 tax rates. The investment of the heat-integrated process is higher than that of the nonintegrated process and increases with the increase of the CO2 emission tax. It is also higher when a cheaper feed stream with a lower content of raw material A is chosen.

Process Synthesis under Variable CO2 Emission Tax Rate
The first part of the study has shown that the values of the CO2 tax rate have a significant impact on the topological and operational variables in optimal process solutions. Since the value of the CO2 emission tax is uncertain in the future, the next part of the study tested different approaches to account for this uncertainty in the process synthesis. The process synthesis problem was extended by the following strategies: 1. Considering multiple values of CO2 tax rates within one period; 2. Introducing multiple periods with specific tax values to synthesize an optimal process over a longer period; 3. Performing a deterministic synthesis of the process for comparison.

One-Period Process Synthesis under Variable CO2 Tax Rate
The aim of this part of the study was to set up a one-period process model with simultaneous consideration of several CO2 taxes (scenarios) and to test how different probabilities of these scenarios affect the optimal result. Three CO2 tax values, 10 EUR/t, 60 EUR/t, and 120 EUR/t, and three combinations of probabilities within a period were assumed, as shown in Table 3, to investigate the impact of the probability distribution on the results of the multi-scenario process synthesis. The MINLP synthesis model for the heat-integrated process was solved for each combination simultaneously at three emission tax values. Expected profits fell enormously as the probability of a higher CO2 tax value increased.  The flow rate of the feed stream at different CO 2 tax rates is shown in Figure 5. The heat-integrated process requires less raw material than the nonintegrated process, except at very low CO 2 tax rates, when a cheaper feedstock is chosen in the integrated process. Higher CO 2 tax rates encourage lower feedstock consumption and more efficient reactant conversion. Figure 6 shows the investment under different CO 2 tax rates. The investment of the heat-integrated process is higher than that of the nonintegrated process and increases with the increase of the CO 2 emission tax. It is also higher when a cheaper feed stream with a lower content of raw material A is chosen.

Process Synthesis under Variable CO 2 Emission Tax Rate
The first part of the study has shown that the values of the CO 2 tax rate have a significant impact on the topological and operational variables in optimal process solutions. Since the value of the CO 2 emission tax is uncertain in the future, the next part of the study tested different approaches to account for this uncertainty in the process synthesis. The process synthesis problem was extended by the following strategies:

1.
Considering multiple values of CO 2 tax rates within one period; 2.
Introducing multiple periods with specific tax values to synthesize an optimal process over a longer period; 3.
Performing a deterministic synthesis of the process for comparison.

One-Period Process Synthesis under Variable CO 2 Tax Rate
The aim of this part of the study was to set up a one-period process model with simultaneous consideration of several CO 2 taxes (scenarios) and to test how different probabilities of these scenarios affect the optimal result. Three CO 2 tax values, 10 EUR/t, 60 EUR/t, and 120 EUR/t, and three combinations of probabilities within a period were assumed, as shown in Table 3, to investigate the impact of the probability distribution on the results of the multi-scenario process synthesis. The MINLP synthesis model for the heat-integrated process was solved for each combination simultaneously at three emission tax values. Expected profits fell enormously as the probability of a higher CO 2 tax value increased. Since emissions and associated taxes occur over a longer period of time, it may be useful to treat multiple periods in the model. The duration of the period can be arbitrarily determined as the period over which a particular distribution of the CO 2 tax rate is expected. In this study, we introduced three periods, and in each of these periods, three different values of the CO 2 tax were assumed with corresponding probabilities (Table 4). Lower CO 2 tax values would tend to occur in the early periods, while higher CO 2 tax rates would take effect later. The process synthesis problem was solved simultaneously for nine scenarios. The expected profit (loss) of the nonintegrated process was −57,324 k EUR/y (Table 5), the investment was 1610 kEUR, and the more expensive raw material and the reactor with higher conversion were chosen. The same choice was made in the case of the heat-integrated process, but in this case, the investment and the expected profit were significantly higher.

Deterministic Synthesis under a Variable CO 2 Tax Rate
The average value of the CO 2 emission tax over the three periods, calculated based on the probabilities in Table 4, is 56 EUR/t. For this value of the CO 2 tax, the first-stage variables, i.e., the choice of the feed stream, reactor type, and its size, were determined. Then, for the fixed values of the first-stage variables and the investment obtained, the second-stage variables (flow rates of the streams, hot and cold utility consumptions) were calculated for all individual values of the tax rate. The problems were solved in the Processes 2021, 9, 1720 11 of 16 framework of the General Algebraic Modeling System (GAMS) [34], using the DICOPT solver as the MINLP solution algorithm. The results are shown in Tables 6 and 7.   The profits obtained at different values of CO 2 tax rates were multiplied by the total probability of an individual CO 2 tax value over all three periods, yielding deterministic expected values of the profit.
Synthesis of nonintegrated process: The expected loss calculated by the deterministic approach, in the case without heat integration, (−57,482 kEUR/y) is slightly more negative than the expected loss of the stochastic approach (−57,324 kEUR/y) reported in Table 5. The expected profit of the optimal heat-integrated process obtained by the deterministic approach (50,999 kEUR/y) is lower than the expected profit of the process obtained by the stochastic approach (51,253 kEUR/y, Table 5).
In both cases, the stochastic approach gave a slightly better result than the deterministic approach, indicating the advantage and importance of stochastic optimization in process synthesis under uncertain CO 2 tax rates.

Synthesis of Heat Exchanger Network under Variable CO 2 Emission Tax Rate
The effect of the tax rate on CO 2 emissions has also been studied in the synthesis of the heat exchanger network (HEN), where very good trade-offs between operating costs and investment are established in an MINLP synthesis model by Yee and Grossmann [25]. The studied HEN consisted of two hot streams and three cold streams [20]. Inlet temperatures (T s ), outlet temperatures (T t ), heat capacity flow rates (CF), and heat transfer coefficients (α) are shown in Table 8. The objective function was to minimize the total annual cost (TAC), which consisted of the annual investment cost (depreciation), operating cost, and CO 2 emission tax. The GHG emissions came from the consumption of the hot utility and from the manufacturing of the heat exchangers, the latter being estimated according to the investment of HEN [35].
The values of CO 2 emission tax in the three periods were provided in the range from 10 EUR/t to 120 EUR/t. Five values (15.2, 35.4, 65.0, 94.8, and 114.8) EUR/tCO 2 that correspond to zeros of the Legendre polynomial of the fifth order were used [36]. It was assumed that lower CO 2 tax values would be more likely at the beginning, while higher CO 2 tax values are more likely in the later period. For the probabilities of the CO 2 tax values, a beta distribution was assumed defined by the parameters (α, β) for periods 1, 2, and 3 as follows: (2, 6), (3,3), and (7,3). The problems were solved in GAMS [34] using an SBB solver.

Synthesis of HEN for Equal Time Periods
In the first step, we assumed three time periods of equal length, which means that each period contributes the same proportion (1/3) to the expected value. The results obtained with the four approaches, described in Section 2, are presented in Table 9. The flexible heat exchanger network determined using the deterministic approach ( Figure 7) has a total area of 1083 m 2 , the expected hot utility consumption is 4354 kW, and the expected cost over three periods is 1,730,303 EUR/y. The network consists of seven process units, of which four are heat exchangers, two are coolers, and one is a heater. The same network topology was obtained using a modified deterministic approach, but the total area was increased by 2.3%, or 25 m 2 . The flexible heat exchanger network determined using the deterministic approach ( Figure 7) has a total area of 1083 m 2 , the expected hot utility consumption is 4354 kW, and the expected cost over three periods is 1,730,303 EUR/y. The network consists of seven process units, of which four are heat exchangers, two are coolers, and one is a heater. The same network topology was obtained using a modified deterministic approach, but the total area was increased by 2.3%, or 25 m 2 .  A flexible network obtained using a multi-scenario approach with Gaussian quadrature ( Figure 8) has a total area of 985 m 2 , the expected hot utility consumption is 4018 kW, and the total expected cost is 1,646,609 EUR/y. A flexible network obtained using a multi-scenario approach with Gaussian quadrature ( Figure 8) has a total area of 985 m 2 , the expected hot utility consumption is 4018 kW, and the total expected cost is 1,646,609 EUR/y. A network obtained with a stochastic approach has one more unit than with a deterministic approach but a lower area and utility consumption. The total GHG emissions of the stochastic approach are 7.4% lower than the GHG emissions of the deterministic approach.

Synthesis of HEN for Different Time Periods
In the second step, we assumed that the three periods are of different lengths, i.e., the first period is the longest and the last the shortest, since it is realistic to assume that higher taxes will not be enforced very soon. We assumed that the duration of the periods is in the ratio 2.333:1.666:1. The results are shown in Table 10.  Tables 9 and 10 shows that the expected TAC values are lower in the case in which a longer duration of the first period with a lower value of the CO2 tax was assumed. In the case of a stochastic multi-period approach, the topology of the network (Figure 9) is different from Figure 8. The network consists of eight process units, of which three are heat exchangers, two are coolers, and three are heaters. The total area is slightly lower. In the case of the deterministic approach, the network topology is A network obtained with a stochastic approach has one more unit than with a deterministic approach but a lower area and utility consumption. The total GHG emissions of the stochastic approach are 7.4% lower than the GHG emissions of the deterministic approach.

Synthesis of HEN for Different Time Periods
In the second step, we assumed that the three periods are of different lengths, i.e., the first period is the longest and the last the shortest, since it is realistic to assume that higher taxes will not be enforced very soon. We assumed that the duration of the periods is in the ratio 2.333:1.666:1. The results are shown in Table 10. A comparison of the results in Tables 9 and 10 shows that the expected TAC values are lower in the case in which a longer duration of the first period with a lower value of the CO 2 tax was assumed. In the case of a stochastic multi-period approach, the topology of the network (Figure 9) is different from Figure 8. The network consists of eight process units, of which three are heat exchangers, two are coolers, and three are heaters. The total area is slightly lower. In the case of the deterministic approach, the network topology is the same as in Figure 7.

Variables
Variables One Period Multi-Period A comparison of the results in Tables 9 and 10 shows that the expected TAC values are lower in the case in which a longer duration of the first period with a lower value of the CO2 tax was assumed. In the case of a stochastic multi-period approach, the topology of the network (Figure 9) is different from Figure 8. The network consists of eight process units, of which three are heat exchangers, two are coolers, and three are heaters. The total area is slightly lower. In the case of the deterministic approach, the network topology is the same as in Figure 7. Both variants of this case study confirmed that a simultaneous multi-scenario approach leads to lower total cost and lower emissions than a deterministic approach. It is, therefore, important to use stochastic methods to integrate the CO 2 tax into the long-term synthesis of optimal processes.

Conclusions
In this paper, the influence of the CO 2 emission tax rate on the synthesis of chemical processes and other process systems was presented using the MINLP approach. Methods for synthesizing long-term optimal and flexible processes considering the uncertain CO 2 tax rate were also demonstrated, ranging from a deterministic approach to a single-period stochastic approach, to a multi-period multi-scenario stochastic approach.
The results show that an increase in the CO 2 tax increases the conversion of reactants, reduces the consumption of raw materials, increases investment, and, consequently, worsens the economic performance of the process. A higher tax rate encourages the use of more efficient process units, e.g., reactors with higher conversion, raw materials with lower levels of impurities, networks with a higher degree of integration.
The synthesis of a flexible process under an uncertain CO 2 tax rate was formulated as a deterministic and stochastic problem, as well as a multi-scenario and multi-period problem that takes into account the longer time period in which the CO 2 tax changes.
The result of this approach is an optimal flexible process that can adjust over time to the variable value of the CO 2 tax. The study showed the advantage of the stochastic approach over the deterministic one, as the stochastic approach achieves better-expected values of the optimization criterion. Greenhouse gas emissions were also noticeably lower, by 7% in the case of the HEN synthesis.

Conflicts of Interest:
The authors declare no conflict of interest.