Carbon Emissions Constrained Energy Planning for Aluminum Products

: The production of primary aluminum is an energy-intensive industry which produces large amounts of direct and indirect greenhouse gas emissions, especially from electricity consumption. Carbon Emissions Constrained Energy Planning proved to be an e ﬃ cient tool for reducing energy-related greenhouse gas emissions. This study focuses on energy planning constrained by CO 2 emissions and determines the required amount of CO 2 emissions from electricity sources in order to meet speciﬁed CO 2 emission benchmark. The study is demonstrated on and applied to speciﬁc aluminum products, aluminum slugs and aluminum evaporator panels. Three di ﬀ erent approaches of energy planning are considered: (i) an insight-based, graphical targeting approach, (ii) an algebraic targeting approach of cascade analysis, and (iii) an optimization-based approach, using a transportation model. The results of the three approaches show that approximately 2.15 MWh of fossil energy source should be replaced with a zero-carbon or 2.22 MWh with a low-carbon energy source to satisfy the benchmark of CO 2 emissions to produce 1 t of aluminum slug; however, this substitution results in higher costs. This study is the ﬁrst of its kind demonstrated on and applied to speciﬁc aluminum products, and represents a step forward in the development of more sustainable practices in this ﬁeld.


Introduction
The current increase in global greenhouse gas (GHG) emissions should be substantially reversed to mitigate climate change and prevent global temperatures from rising beyond the 2 • C target [1]. Currently, global GHG emissions are attributed mainly to energy (more than 70% of GHG emissions), while the remaining are due to agriculture, land-use change and forestry, and industrial processes and waste [2]. In the European Union (EU), a binding target has been set to reduce GHG emissions by at least 40% by 2030 compared to 1990 levels [3].
In the EU, energy-intensive industries including iron, steel, cement and aluminum either have to pay carbon taxes or, if they are included in the EU's Emissions Trading System (EU ETS), they need to have sufficient allowances to cover all emissions produced [4]. Under EU ETS, installations could receive free allocation allowances every year and/or buy (remaining) emission allowances on the market. Every year, aluminum producers in the EU receive some allowances allocated for free, owing to international competitiveness of carbon leakage [5]. The number of free allocations is being decreased Sinha and Chaturvedi [40] reviewed carbon reduction technologies in industries with a focus on carbon emission limits planning.
Based on the literature review, it was found that none of the studies focused on CCEP applied to the aluminum industry and to specific aluminum products. It is also worth noting that most of the CCEP studies considered only one approach to energy planning for specific applications, while this study considers three different approaches: the first is the graphical approach of CEPA, the second is the algebraic approach implementing a cascade analysis technique, and the third is the optimization-based approach, applying a transportation model [13].
The aim of this study is the determination of the required amount of CO 2 emissions from electricity sources to meet specified CO 2 emission limits, considering the aforementioned approaches. All three approaches are systematic, which provides decision support in determining an appropriate energy mix to achieve a specific target, by considering the characteristics of energy sources and demands. CCEP techniques provide support in planning for more sustainable production of industrial products. They can also be used to optimize the energy supply mix based on carbon emission constraints subject to economic and other environmental constraints, except for carbon footprint and the like. This study focuses on the production of more sustainable aluminum products, and is the first CCEP work applied to the aluminum industry. CCEP approaches are applied to one aluminum product, i.e., aluminum slugs (graphical and numerical approaches), and two aluminum products, i.e., slugs and evaporator panels (all three approaches). All three approaches are finally compared and discussed.

Carbon Emission-Constrained Electricity Planning for Aluminum Products
Aluminum is widely available and is the second most widely used metal worldwide. Aluminum is used in several sectors, including transportation, packaging, building and construction, furniture, equipment and machinery. However, the aluminum industry is highly energy-consuming, especially in the production of primary aluminum. Aluminum production consumes about 3.5% of electricity globally [41], while in China [42] and the US [43], that figure is around 5% and in Montenegro it is more than 25% [44]. According to data from the Slovenian Ministry of Infrastructure [45] and the data collected by the authors, aluminum production in Slovenia consumed about 8.9% of the total electricity produced in 2017. Primary production of aluminum starts with bauxite mining and further refining of alumina from bauxite. Alumina is further processed into aluminum via electrolysis, known also as the Hall-Héroult process. Subsequently, aluminum is processed into a range of products by casting, rolling, extrusion and other operations [46]. The highest energy consumption and CO 2 emissions are due to electrolysis [47]. The production of secondary or recycled aluminum, on the other hand, requires less than 5% of the energy required to produce aluminum from ore [48].
The aim of this study is to minimize CO 2 emissions in order to meet specified CO 2 emission limits, which could be set at the national or regional level. The emission limit (target) for this study was set at 0.376 t CO 2 /MWh [5]. This value was set by the European Commission in delegated regulation (EU) 2019/331, where transitional rules for the free allocation of emission allowances were set for the period 2021-2030. This emission factor is used for the determination of indirect emissions, where a given factor (0.376 t CO 2 /MWh) is multiplied by the electricity consumption required to produce specific products.
In this study, two aluminum products are considered, aluminum slugs and evaporator panels. The production of aluminum slugs and evaporator panels is described in Figure 1. Both products are made from molten aluminum produced by the electrolysis process, while slugs are also made from recycled aluminum produced as process waste material from the further processing of slugs. This study performs CCEP of aluminum products for a Slovenian aluminum company. In the company, aluminum slugs are produced from around 0.67 t of electrolysis aluminum per t of slugs; the remaining material used for their production is secondary aluminum. Aluminum evaporator panels are produced only from electrolysis aluminum, where from 1 t of molten aluminum, about 1 t of evaporator panels are produced. panels according to the needs of the customer. Both production processes follow the same value chain up to slug and evaporator panel production. For the sake of simplification, the same amount of the primary aluminum used in the production of aluminum slugs and evaporator panels is assumed; 0.67 t of primary aluminum is required for the production of 1 t of aluminum slugs, and the same amount of primary aluminum is used for 0.67 t of aluminum evaporator panels. To produce 0.67 t of primary aluminum, 0.063 MWh of electricity are required for the anode production and 9.12 MWh for the electrolysis. This consumption is the same for both products, since they follow the same value chain up to the final steps. They have similar electricity consumption, although the production of evaporator panels requires slightly more electricity. For slug production, an additional 0.263 MWh of electricity is required for 1 t of the product, while for evaporator panels, an additional 1.016 MWh of electricity is required for 0.67 t of the product. Total electricity consumption for 1 t of slugs is 9.446 MWh, and for 0.67 t of evaporator panels is 10.199 MWh.
The electricity delivered to a company is generated from a mix of fossil, nuclear and renewable energy sources. Electricity sources of fossil origin are coal, lignite, natural gas and oil. Hydro, wind, photovoltaics, geothermal, biomass, biogas and biodiesel are considered types of renewable energy. For the sake of simplification, energy sources are grouped into three main categories. Fractions of energy sources in the electricity mix, and emission factors of the different sources are summarized in Table 1. Average values for prices and emission factors are used for each energy source. The prices for each energy source are regarded as confidential and thus they are not presented. Three different energy planning approaches (graphical, algebraic and optimization-based approaches) will be demonstrated in the following sections. The first two approaches are presented and illustrated on one aluminum product (aluminum slugs) and are shown for zero-and low-carbon energy sources (in the current work, the low-carbon energy source is renewable energy). The third approach deals with two aluminum products (aluminum slugs and evaporator panels) and considers two cases: the case where electricity is pinched for each product, and the hypothetical case where only one Pinch exists for both products. In the case of optimization, the trade-off solution is obtained between emissions and the cost of energy sources to achieve a given benchmark.
Graphical and numerical approaches consist of several steps and follow similar algorithms whose detailed steps are described in related previous papers. First, energy supply and demand are defined; then, they are arranged into emission intervals (segments) by increasing the emission factors Aluminum slug production (see Figure 1) consists of casting, the production of the narrow strip, hot and cold rolling, stamping, annealing, surface treatment and packaging [47]. The production process of Roll-Bond evaporator panels also consists of several steps: casting of the wide strip, rolling it into the band for evaporator panels, roll-bonding, recrystallization annealing, inflation of imprinted channels, final dimension cut or stamping and packaging, and final assembly of the evaporator panels according to the needs of the customer. Both production processes follow the same value chain up to slug and evaporator panel production.
For the sake of simplification, the same amount of the primary aluminum used in the production of aluminum slugs and evaporator panels is assumed; 0.67 t of primary aluminum is required for the production of 1 t of aluminum slugs, and the same amount of primary aluminum is used for 0.67 t of aluminum evaporator panels. To produce 0.67 t of primary aluminum, 0.063 MWh of electricity are required for the anode production and 9.12 MWh for the electrolysis. This consumption is the same for both products, since they follow the same value chain up to the final steps. They have similar electricity consumption, although the production of evaporator panels requires slightly more electricity. For slug production, an additional 0.263 MWh of electricity is required for 1 t of the product, while for evaporator panels, an additional 1.016 MWh of electricity is required for 0.67 t of the product. Total electricity consumption for 1 t of slugs is 9.446 MWh, and for 0.67 t of evaporator panels is 10.199 MWh.
The electricity delivered to a company is generated from a mix of fossil, nuclear and renewable energy sources. Electricity sources of fossil origin are coal, lignite, natural gas and oil. Hydro, wind, photovoltaics, geothermal, biomass, biogas and biodiesel are considered types of renewable energy. For the sake of simplification, energy sources are grouped into three main categories. Fractions of energy sources in the electricity mix, and emission factors of the different sources are summarized in Table 1. Average values for prices and emission factors are used for each energy source. The prices for each energy source are regarded as confidential and thus they are not presented. Three different energy planning approaches (graphical, algebraic and optimization-based approaches) will be demonstrated in the following sections. The first two approaches are presented and illustrated on one aluminum product (aluminum slugs) and are shown for zero-and low-carbon energy sources (in the current work, the low-carbon energy source is renewable energy). The third approach deals with two aluminum products (aluminum slugs and evaporator panels) and considers Energies 2020, 13, 2753 5 of 18 two cases: the case where electricity is pinched for each product, and the hypothetical case where only one Pinch exists for both products. In the case of optimization, the trade-off solution is obtained between emissions and the cost of energy sources to achieve a given benchmark.
Graphical and numerical approaches consist of several steps and follow similar algorithms whose detailed steps are described in related previous papers. First, energy supply and demand are defined; then, they are arranged into emission intervals (segments) by increasing the emission factors from the lowest to the highest. By the graphical approach, the energy source and demand curves are then plotted on the energy consumption-emissions diagram. The energy source curve is further shifted horizontally until it intersects with the energy demand curve (at the "Pinch Point") [12]. By the numerical approach, the energy surplus/deficit is cascaded through the intervals, and the cumulative net energy consumption should be non-negative. The deficit should then be added to the zero-or low-carbon energy source [13]. The optimization approach, on the other hand, involves formulating a working model whose solution is determined numerically via established solution techniques (i.e., Simplex algorithm) that is embedded in standard optimization software. The data and variables should be defined, and together with the constraints, they form a feasible region. The optimal solution is then defined based on the given objective, which is minimized or maximized. The procedure for each approach is described in more detail, together with the case studies, in the following section.

Graphical Targeting Approach
The first technique for CCEP is the insight-based graphical targeting approach, called CEPA, originally developed in 2007 by Tan and Foo [12]. CEPA uses a graphical approach based on the principles of conventional PA [21]; it enables the identification of the minimum amount of low-or zero-carbon energy needed to satisfy demand-side emission constraints (target) [49]. Note that the term "zero-carbon" is usually applied to sources with very low CO 2 intensity compared to fossil fuels, even if the actual magnitude is nonzero. Various applications and extensions of the CEPA methodology have been developed recently, such as applications to energy sectors in different countries (Ireland [50], New Zealand [51], USA [52], China [53], Nigeria [54], the Baltic States [55], the EU [56] and others), alternative metrics and footprints (land [13], water [15]), and other sectors, such as transport systems [16], economic systems [17] and others. Conceptual or graphical techniques applied in PA are useful tools in the preliminary stages of energy planning, energy policy [57], in graphical representation, in step-by-step user control and in the verification of results. The graphical display also provides visualization that is useful for the analysis of a problem and subsequent communication of the results. On the other hand, graphical methods are limited to relatively simple problems [12]. Another limitation could be the accuracy of the results, which depends on the quality of the graphical display [12].
All the steps of the CEPA algorithm are presented in detail by Tan and Foo [12]. The first step in the CEPA procedure is the preparation of the energy source and demand data, which are shown in Table 2. Energy sources are sorted in order of increasing emission factors (see also Table 1). The consumption described in Table 2 (in MWh/t) was obtained by multiplying the fraction of the energy source by the total consumption of electricity for the production of aluminum slugs (9.446 MWh/t). As stated previously, only one product (aluminum slugs) was considered for the graphical and algebraic approaches. Emissions (in t CO 2 /t) were obtained by multiplying the emission factor and total electricity consumption for 1 t of slugs. Energy demand represents the total consumption of electricity for 1 t of slugs with target emission constraints (benchmark of 0.376 t CO 2 /MWh [5]). Emissions (3.552 t CO 2 /t) were obtained by multiplying the benchmark emission factor and total electricity consumption.
The second step in the CEPA procedure is the generation of source and demand composite curves (CCs). Source curves for energy supply were obtained as CO 2 emissions vs. electricity supply plots, and demand curves as CO 2 emissions vs. electricity demand plots. The horizontal axis represents electricity supply or demand, while CO 2 emissions are plotted as the vertical axis. The slope of each electricity source is equal to its emission factor, while the slope of electricity demand is equal to Energies 2020, 13, 2753 6 of 18 benchmark emission factor. Fossil energy has the highest emission factor (1.015 t CO 2 /MWh; see also Table 1) of all electricity sources and a much steeper slope compared to renewable and nuclear energy. The next step is to construct source and demand CCs (see Figures 2 and 3). The source curve consists of three segments, i.e., nuclear, renewable, and fossil segments, and is constructed by arranging the segments end-to-end in ascending order of emission factors. The resulting source CC curls upwards [12]. The demand curve consists of only one segment, i.e., target CO 2 emissions for 1 t of aluminum slugs. The demand curve is thus linear over the entire interval.
To supply suitable energy sources that meet the CO 2 emissions targets, the source CC is then shifted to the right. The Pinch Point is where the demand curve touches the shifted source curve. On the lower end of the shifted source CC (minimum distance between origin and shifted CC) is the minimum amount of zero-carbon energy sources needed to meet the CO 2 emission limit, while at the upper end (distance between the Pinch Point and the end of the shifted source CC) is excess energy from fossil sources. At the lower end, it is generally desirable to maximize the use of zeroor low-carbon sources [58], while at the upper end, it is desirable to minimize the use of the most carbon-intensive energy source. Figure 2 shows a Pinch diagram for 1 t of aluminum slug production while minimizing the zero-carbon energy source.
The minimum amount of zero-carbon energy ( Figure 2) in this study is the same as the excess of the fossil source because of the same amount of electricity supply and demand. The graphical approach shows that for 1 t of aluminum slugs, about 2.14 MWh of the fossil source should be replaced with zero-carbon energy to achieve an emission limit of 3.552 t CO 2 . As is shown in the left-upper window of Figure 2, the zero-carbon source has a zero slope, and nuclear energy, which is the lowest carbon emission source, is shifted for the minimum zero-carbon energy source to the right.
Energies 2020, 13, x FOR PEER REVIEW 6 of 18 The second step in the CEPA procedure is the generation of source and demand composite curves (CCs). Source curves for energy supply were obtained as CO2 emissions vs. electricity supply plots, and demand curves as CO2 emissions vs. electricity demand plots. The horizontal axis represents electricity supply or demand, while CO2 emissions are plotted as the vertical axis. The slope of each electricity source is equal to its emission factor, while the slope of electricity demand is equal to benchmark emission factor. Fossil energy has the highest emission factor (1.015 t CO2/MWh; see also Table 1) of all electricity sources and a much steeper slope compared to renewable and nuclear energy.
The next step is to construct source and demand CCs (see Figures 2 and 3). The source curve consists of three segments, i.e., nuclear, renewable, and fossil segments, and is constructed by arranging the segments end-to-end in ascending order of emission factors. The resulting source CC curls upwards [12]. The demand curve consists of only one segment, i.e., target CO2 emissions for 1 t of aluminum slugs. The demand curve is thus linear over the entire interval.
To supply suitable energy sources that meet the CO2 emissions targets, the source CC is then shifted to the right. The Pinch Point is where the demand curve touches the shifted source curve. On the lower end of the shifted source CC (minimum distance between origin and shifted CC) is the minimum amount of zero-carbon energy sources needed to meet the CO2 emission limit, while at the upper end (distance between the Pinch Point and the end of the shifted source CC) is excess energy from fossil sources. At the lower end, it is generally desirable to maximize the use of zero-or lowcarbon sources [58], while at the upper end, it is desirable to minimize the use of the most carbonintensive energy source. Figure 2 shows a Pinch diagram for 1 t of aluminum slug production while minimizing the zero-carbon energy source. The minimum amount of zero-carbon energy ( Figure 2) in this study is the same as the excess of the fossil source because of the same amount of electricity supply and demand. The graphical approach shows that for 1 t of aluminum slugs, about 2.14 MWh of the fossil source should be replaced with zero-carbon energy to achieve an emission limit of 3.552 t CO2. As is shown in the left- low emission factor, with advantages such as renewability. Figure 3 shows a Pinch diagram for 1 t of aluminum slugs with a minimum low-carbon energy source (renewable energy). In the low-carbon case, about 2.22 MWh is the minimum low-carbon energy source, and the same amount of energy is the excess of the fossil source. Owing to the addition of the low-carbon source, slightly more excess energy from the fossil source with the highest emission factor is obtained compared to the zerocarbon example. Both cases (zero-and low-carbon energy source) will also be presented using the algebraic targeting approach.

Algebraic targeting approach (cascade analysis)
The algebraic targeting approach with the cascade analysis technique was originally developed for resource conservation networks [59], and has its roots in the Problem Table Algorithm and the Heat Cascade developed for Maximum Heat Recovery networks [21]. Cascade analysis for CCEP is an algebraic targeting approach for the identification and determination of the minimum amount of zero-and/or low-carbon sources to achieve emission limits. The method was introduced by Foo et al. [13] and could overcome the inaccuracy problems of the graphical approach.
Cascade analysis for determining the minimum amount of zero-or low-carbon sources consists of various steps [13], all of which are presented in detail by Foo et al. [13]. In the following section, the steps of the algorithm are demonstrated on an example which targets electricity demand in the production of 1 t of aluminum slugs.
In the first column of the cascade analysis table (see Tables 3-6 However, it should be noted that it is virtually impossible to produce electricity with no CO 2 emissions (zero-carbon electricity). Even renewable energy sources can only approach CO 2 neutrality, but they do not reach it [58]. In many cases, it is more desirable to minimize the low-carbon energy source. Renewable energy is assumed to be a low-carbon energy source in this study, as it exhibits a low emission factor, with advantages such as renewability. Figure 3 shows a Pinch diagram for 1 t of aluminum slugs with a minimum low-carbon energy source (renewable energy). In the low-carbon case, about 2.22 MWh is the minimum low-carbon energy source, and the same amount of energy is the excess of the fossil source. Owing to the addition of the low-carbon source, slightly more excess energy from the fossil source with the highest emission factor is obtained compared to the zero-carbon example.
Both cases (zero-and low-carbon energy source) will also be presented using the algebraic targeting approach.

Algebraic Targeting Approach (Cascade Analysis)
The algebraic targeting approach with the cascade analysis technique was originally developed for resource conservation networks [59], and has its roots in the Problem Table Algorithm and the Heat Cascade developed for Maximum Heat Recovery networks [21]. Cascade analysis for CCEP is an algebraic targeting approach for the identification and determination of the minimum amount of zeroand/or low-carbon sources to achieve emission limits. The method was introduced by Foo et al. [13] and could overcome the inaccuracy problems of the graphical approach.
Cascade analysis for determining the minimum amount of zero-or low-carbon sources consists of various steps [13], all of which are presented in detail by Foo et al. [13]. In the following section, the steps of the algorithm are demonstrated on an example which targets electricity demand in the production of 1 t of aluminum slugs.
In the first column of the cascade analysis table (see Tables 3-6), emission factors are sorted in ascending order in k intervals (c k ). Their increments between successive rows (∆c k ) are presented in the second column of the cascade analysis table. The following columns represent energy demand (D j ) and energy source (S i ) at the interval of the corresponding emission factor, where the number of intervals k is equal to the number of energy demand (j) and source types (i). The difference between supply and demand (S i − D j ) gives the net surplus or deficit of energy at each interval k. A cumulative energy surplus or deficit (F k = ∆(S i − D j )) is the sum of increments between successive rows, starting from either a zero-carbon or a low-carbon source. In the next column, the CO 2 emission load at interval k (E k ) is calculated as: Next, the cumulative CO 2 load (∆E k ) between successive rows of E k is calculated. If ∆E k has a negative value, the CO 2 limit has been exceeded, and the cascade is infeasible. In the last column, clean source demand at interval k (F CS,k ) is shown, which is calculated as: where c CS represents an emission factor or the clean (zero-or low-carbon) energy source. In the following section, both cases will be shown, i.e., the minimum demand for the zero-and low-carbon energy sources. The cascade analysis approach was implemented in Excel by using the actual (not rounded) numbers for calculations.

Minimum Demand for the Zero-Carbon Source
First, the demand for the zero-carbon source was determined to achieve the specified emission limits. The data used in the study are as shown in Table 1. Cascade analysis was applied to one aluminum product, i.e., aluminum slugs, and is shown in Table 3.
The analysis in Table 3 shows an infeasible cascade for a minimum zero-carbon energy source. A feasible cascade should have a non-negative value for the cumulative CO 2 load in the last interval. A feasible cascade is shown in Table 4, where almost all the steps remain the same as before. The exception is that the largest negative value of F CS,k is considered as the quantity of the zero-carbon energy source, and is taken as an absolute positive value.  The negative value of F CS,k in the last interval from Table 3 (−2.137 MWh) is the amount of required zero-carbon energy source S i in the first interval in Table 4 (where emission factor c k is zero). In this way, the cascade is adjusted to become feasible. From Table 4, it can also be seen that the excess energy from the fossil source is the same as the minimum amount of zero-carbon energy (2.137 MWh; last interval of F k ). The Pinch Point is obtained at interval k, where ∆E k is equal to zero, which is at 1.015 t CO 2 /MWh.

Minimum Demand for the Low-Carbon Source
The algebraic targeting approach is further extended to determine the amount of the low-carbon energy source to achieve the specified emission limits. All steps in the cascade analysis remain the same as for the zero-carbon example, except that S i for the renewable energy source is now the sum of the previous amount of renewable energy source (1.425 MWh) plus the minimum amount of zero-carbon energy obtained above (2.137 MWh, see Table 4). F CS,k was obtained by considering a c CS of 0.038 (see also Table 1). The cascade analysis determining the minimum amount of the low-carbon energy source is shown in Table 5. Table 5. Cascade analysis for the low-carbon energy source where not all CO 2 loads are non-negative.  The cascade analysis in Table 5 gives an infeasible solution for the low-carbon source owing to the negative value of ∆E k in the last interval. A feasible solution to the cascade analysis is shown in Table 6, where again, the most negative value of F CS,k (−0.083 MWh) is added to the low-energy source (S i for renewable energy source with c k of 0.038). This is the amount (in absolute numbers) of the required addition to the low-carbon supply shown in Table 5 (3.562 MWh, increased to 3.645).
As obtained also by the graphical approach, the amount of excess energy from fossil sources is increased to 2.22 MWh. The Pinch Point is again obtained at 1.015 t CO 2 /MWh.
In the following section, the optimization-based approach is used to minimize the cost of electricity sources while achieving CO 2 emission targets. Two products are considered for optimization, i.e., aluminum slugs and aluminum evaporator panels.

Optimization-Based Approach
Mathematical Programming (MP) is another broad category of PI methodologies [60]. MP enables the search to be automated, while finding the optimal solution among many design alternatives [61]. The simplest representation of the MP model is defined as a deterministic model (the model where all the parameters are constant and not uncertain or random), and is written in the form of a Linear Program (LP). One advantage over the previous two approaches is the flexibility to change the objective function to suit the context of the problem. For PI problems, one of the first problem formulations of MP was the transportation model [62].
The transportation model (its variants are also known as the source-sink model or resource-allocation problem) is based on the principle that a product is transported from a number of sources to a number of sinks (destinations) based on a given objective (minimum cost, maximum profit, etc.). Dealing with both mass and energy for CCEP problems as extensive properties, in contrast to temperature in Heat Integration, which is an intensive property, the transportation model is a suitable approach for solving CCEP problems. The transportation model, which is written in a generic way and could be applied to any transportation or planning problem, is shown in Equations (3)-(7). The sets, data and variables applied in the transportation model are described on the left.
Sets I source (fossil, nuclear, renewable) with elements i ∈ I J demand (slug, evaporator panels) with elements j ∈ J Parameters a i supply or availability of energy source i (MWh/y) b j demand for electricity to produce product j (MWh) F i emission factor for electricity source i (t CO 2 /MWh) P i price of electricity source i (EUR/MWh)

Positive variable
x ij amount of electricity from source i to demand j (MWh/y) w ij fraction of electricity source i to produce product j (-) E j CO 2 emissions for each product j (t CO 2 /t products) Equation (3) is calculated for each source i, and implies that the sum of specific electricity source i consumed for all j products (x ij ) should be greater than or equal to the availability of that electricity source i (a i ). Equation (4) is calculated for each demand for energy j and requires that the amount of electricity from all the sources satisfy the demand for electricity to produce each product j (b j ). CO 2 emissions originating from electricity consumption for the production of each product j (E j ) are calculated with Equation (5), where the amounts of electricity consumed x ij are multiplied by the emission factors for energy source i (F i ). Equation (6) calculates the fraction of electricity source i (w ij ) used to produce each product j. Finally, Equation (7) represents the objective, which is to minimize the cost of the electricity supply mix, which is calculated by multiplying the amount of electricity consumed (x ij ) by the price of each electricity source i (P i ).
Three different cases are studied: • Current case: The electricity mix for producing aluminum slugs and evaporator panels is fixed to the current shares of each energy source in the mix. To simulate the current case, CO 2 emissions are also fixed to their current emission values. To calculate the cost for the supply of electricity mix, Equations (3)-(7) are applied. • Case 1: It is assumed that the electricity mix delivered to the company allows optimal fractions of the main energy sources to be in the electricity supply mix. Also, emission levels for each aluminum product should be achieved for every product. The upper emission limit is set to 3.55 t CO 2 /t slugs and to 3.835 t CO 2 per 0.67 t of evaporator panels. As explained previously, the emission limit is set to 0.376 t CO 2 /MWh of electricity, while 9.446 MWh are required for 1 t of slugs and 10.199 MWh for 0.67 t of evaporator panels. To calculate the cost for this case, Equations (3)- (7) are applied. where Equation (5) is reformulated so that emissions from electricity consumed do not exceed the upper emission limits for each product j (3.55 t and 3.835 t CO 2 ): • Case 2: This case is similar to Case 1, but here, it is assumed that the total emission level should be achieved for both aluminum products together. The upper emission limits are the same as those for the first case. Equation (5) is again reformulated to consider the upper limit on emissions for both products combined: The models for the three cases determine the minimum cost of the electricity supply mix considering emission limits. The main results of the transportation model for the three cases are summarized in Table 7, where the fractions of the electricity sources used to produce each product and the total cost of electricity supply are shown. It should be noted that, for reasons of confidentiality, the values of cost are normalized, where a value of 1 represents the cost for current electricity mix. Table 7 shows that the cost of the current electricity supply mix is lower by 26% compared to cases considering emission limits (Case 1 and 2). Currently, there is a relatively high share of fossil sources in the electricity mix, which is cheaper compared to nuclear and renewable sources. With the current electricity mix, however, emission levels cannot be achieved because of the higher emission factor of fossil energy (see also Table 1). For both the current case and Case 1, the fractions of electricity sources in the mix are the same for both products. Case 2 is a hypothetical case, where the fractions of sources in the electricity mix for slugs and evaporator panels are different. Renewable energy is not selected in either Case 1 and 2, owing to the higher emission factor and lower cost compared to nuclear energy.
The Pinch diagrams in the next sections show where an electricity consumption Pinch occurs for each product (Case 1), and when only one Pinch for all the products combined should be achieved (Case 2).

Electricity Consumption Pinch for Each Product
Case 1 presents the case whereby for each product, the emission limit should not be exceeded. Electricity consumption is pinched for each product, as shown in Figure 4. The source CC represents the current case, where emissions and the composition of the electricity mix for both products are fixed. The demand curve consists of two segments, one for each product. However, both segments have the same slope, since the CO 2 emission benchmark is the same for both products (0.376 t CO 2 /MWh [5], as stated previously). For Case 1, it can be seen that both products use nuclear and fossil energy to reach the target, as presented in Table 7. Renewable energy is not selected, as it has a higher emission factor and is more expensive than nuclear.
Energies 2020, 13, x FOR PEER REVIEW 12 of 18 The Pinch diagrams in the next sections show where an electricity consumption Pinch occurs for each product (Case 1), and when only one Pinch for all the products combined should be achieved (Case 2).

Electricity Consumption Pinch for Each Product
Case 1 presents the case whereby for each product, the emission limit should not be exceeded. Electricity consumption is pinched for each product, as shown in Figure 4. The source CC represents the current case, where emissions and the composition of the electricity mix for both products are fixed. The demand curve consists of two segments, one for each product. However, both segments have the same slope, since the CO2 emission benchmark is the same for both products (0.376 t CO2/MWh [5], as stated previously). For Case 1, it can be seen that both products use nuclear and fossil energy to reach the target, as presented in Table 7. Renewable energy is not selected, as it has a higher emission factor and is more expensive than nuclear. A sensitivity analysis was further performed by fixing the fraction of renewable energy in the electricity supply mix. It is expected that the share from renewable sources will rise in future, while the expansion of nuclear in the future is less certain. The fraction of renewable energy is increased from 0 % to 100 %. The results of the sensitivity analysis are shown in Figure 5 for fractions of fossil and nuclear sources (left) and cost (right).
It can be seen that the fraction of nuclear energy in the electricity mix decreases almost linearly with the increase in the renewable energy share in the electricity mix up to renewable w reaching about 0.7. Subsequently, nuclear energy is no longer selected because of cheaper fossil electricity, while the emission limit is satisfied. The cost of the electricity supply mix is shown in Figure 5 on the right. For reasons of confidentiality, the values of cost, on the vertical axis, are hidden. The lowest cost was obtained with no renewable energy in the supply mix (as shown in Table 7). It should be noted that renewable energy is the most expensive of all the sources (see Table 1), and the cost of electricity is constantly increasing with increasing the share of renewable w . The cost increases slightly for renewable w < 0.6, where renewable energy replaces nuclear. It begins to increase more significantly for renewable w > 0.6, when renewable replaces fossil energy, owing to the significant differences in price. A sensitivity analysis was further performed by fixing the fraction of renewable energy in the electricity supply mix. It is expected that the share from renewable sources will rise in future, while the expansion of nuclear in the future is less certain. The fraction of renewable energy is increased from 0% to 100%. The results of the sensitivity analysis are shown in Figure 5 for fractions of fossil and nuclear sources (left) and cost (right).
It can be seen that the fraction of nuclear energy in the electricity mix decreases almost linearly with the increase in the renewable energy share in the electricity mix up to w renewable reaching about 0.7. Subsequently, nuclear energy is no longer selected because of cheaper fossil electricity, while the emission limit is satisfied. Fossil source, on the other hand, remains almost constant, between w renewable 0 and about 0.7. At higher values of w renewable , the fraction of w fossil starts decreasing almost linearly with the increase in w renewable . It should be noted that w renewable + w fossil + w nuclear should always be equal to 1.
The cost of the electricity supply mix is shown in Figure 5 on the right. For reasons of confidentiality, the values of cost, on the vertical axis, are hidden. The lowest cost was obtained with no renewable energy in the supply mix (as shown in Table 7). It should be noted that renewable energy is the most expensive of all the sources (see Table 1), and the cost of electricity is constantly increasing with increasing the share of w renewable . The cost increases slightly for w renewable < 0.6, where renewable energy replaces nuclear. It begins to increase more significantly for w renewable > 0.6, when renewable replaces fossil energy, owing to the significant differences in price.
Energies 2020, 13, x FOR PEER REVIEW 13 of 18 Figure 5. Composition of energy mix (left) and cost (right) while increasing the fraction of renewable energy in the electricity mix.

Electricity Consumption Pinch for all the Products Considered
Case 2 presents a hypothetical case where for all products, the overall emission limit should not be exceeded. Only one electricity Pinch is obtained for the total electricity demand, as shown in Figure  6. The demand curve is the same as in Case 1, while the source curves for the Current case and Case 2 are now different. Both source curves are arranged in ascending order of emission factors for all the products combined. It can be seen that for slugs, only nuclear energy is used, while for evaporator panels, first nuclear energy is used, and later, also fossil energy to satisfy benchmark CO2 emissions. It should be noted that even here, fossil energy is the preferred source because it has the lowest cost; thus Cases 1 and 2 both exhibit the same cost, as shown in Table 7.

Conclusions
In this work, an extension of CEPA was developed and applied to the energy-and carbonintensive aluminum industry. Graphical, algebraic and optimization-based approaches for CCEP were applied for the production of specific aluminum products, i.e., aluminum slugs and evaporator panels. Graphical and algebraic approaches determined the minimum amount of zero-and lowcarbon energy sources required to achieve the specified emission limit. The transportation model was

Electricity Consumption Pinch for All the Products Considered
Case 2 presents a hypothetical case where for all products, the overall emission limit should not be exceeded. Only one electricity Pinch is obtained for the total electricity demand, as shown in Figure 6. The demand curve is the same as in Case 1, while the source curves for the Current case and Case 2 are now different. Both source curves are arranged in ascending order of emission factors for all the products combined. It can be seen that for slugs, only nuclear energy is used, while for evaporator panels, first nuclear energy is used, and later, also fossil energy to satisfy benchmark CO 2 emissions. It should be noted that even here, fossil energy is the preferred source because it has the lowest cost; thus Cases 1 and 2 both exhibit the same cost, as shown in Table 7.

Electricity Consumption Pinch for all the Products Considered
Case 2 presents a hypothetical case where for all products, the overall emission limit should not be exceeded. Only one electricity Pinch is obtained for the total electricity demand, as shown in Figure  6. The demand curve is the same as in Case 1, while the source curves for the Current case and Case 2 are now different. Both source curves are arranged in ascending order of emission factors for all the products combined. It can be seen that for slugs, only nuclear energy is used, while for evaporator panels, first nuclear energy is used, and later, also fossil energy to satisfy benchmark CO2 emissions. It should be noted that even here, fossil energy is the preferred source because it has the lowest cost; thus Cases 1 and 2 both exhibit the same cost, as shown in Table 7.

Conclusions
In this work, an extension of CEPA was developed and applied to the energy-and carbonintensive aluminum industry. Graphical, algebraic and optimization-based approaches for CCEP were applied for the production of specific aluminum products, i.e., aluminum slugs and evaporator panels. Graphical and algebraic approaches determined the minimum amount of zero-and lowcarbon energy sources required to achieve the specified emission limit. The transportation model was

Conclusions
In this work, an extension of CEPA was developed and applied to the energy-and carbon-intensive aluminum industry. Graphical, algebraic and optimization-based approaches for CCEP were applied for the production of specific aluminum products, i.e., aluminum slugs and evaporator panels. Graphical and algebraic approaches determined the minimum amount of zero-and low-carbon energy sources required to achieve the specified emission limit. The transportation model was further applied to energy planning for two specific products to calculate the optimal electricity mix to reach the specified emission limit with minimum cost.
All the approaches yielded similar results; the first two showed that about 2.14 MWh of zero-carbon or 2.22 MWh of low-carbon, and the same amounts of excess energy (2.14 MWh or 2.22 MWh), are obtained for 1 t of aluminum slugs. The third approach was based on cost optimization and showed that 26% higher electricity cost would be required to achieve the specified CO 2 emission target compared to the current case. In principle, this incremental cost can potentially be covered by economic incentives such as an emissions tax. CCEP can thus be used in the future to estimate the amount of such incentives needed to drive industrial fossil-based GHG emissions abatement.
Each approach showed its strengths and weaknesses, as summarized in Table 8. The graphical and algebraic approaches proved to be powerful tools for energy planning of industrial products, since they are intuitive and provide better insights into the problem for both analysis and communication.
Optimization-based results are less intuitive; however, optimization automates the procedure and could include more details, including cost optimization. Hybrid approaches are suggested which combine different methods to yield the synergistic advantages of each separate method [63]. Table 8. Comparison of the approaches.

Approach Strengths Weaknesses
Graphical -useful approach for preliminary planning -enables graphical (visual) presentation -useful for communication of results -enables user control -enables verification of results -limited to simple problems -limited accuracy of results -size and unit dependence -not flexible -procedure is not automated -cannot include details and objectives, such as cost optimization Algebraic -application to both simple and complex problems -precision of results -more systematic approach -provides insight into the problem -useful for communication of results -size and unit dependence -not flexible -procedure is not automated -cannot include details and objectives, such as cost optimization Optimization -application to both simple and complex problems -precision of results -flexible -automated procedure -equations can be modified -could include more details and various objectives -relatively simple input/code determination -requires background of mathematical modeling -more difficult communication of results -verification of results is less straightforward This work demonstrated the usefulness of CCEP approaches in the move towards a more sustainable aluminum industry. Since aluminum is widely used in various applications and is highly energy intensive, optimizing the electricity supply mix used in its production could significantly contribute to conserving resources and decreasing the global carbon footprint. However, regarding the optimization of the electricity mix, more detailed exploration regarding the electricity sector should be performed in terms of the cost, emissions, sustainability and availability of low-, zero-and negative-carbon emission sources.
In the future, CCEP analyses will be extended for targeting GHG (CO 2 Equation emissions) and other footprints. Various options for carbon capture could be incorporated, including technologies to achieve GHG emissions reductions. Additionally, CO 2 emissions reductions could be analyzed by including mass integration to minimize emissions and the use of materials. Analyses could also be expanded to more aluminum products and to include more details regarding types of energy sources (renewable energy could be separated to hydro, photovoltaics, geothermal, wind, biomass, etc.), and their availability and pricing (e.g., hourly based), including future prices. Multiscale application of PI techniques could be used to improve material and energy efficiency to reduce energy consumption and emissions at the process, TS and supply chain network levels [33]. The methodology will be extended to account for variations in the ratio of primary and secondary aluminum within a circular economy framework and to account for negative emission technologies, and will include the combined use of all three proposed approaches.

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