The steel industry is one of the technological pillars of the modern world. In 2024, more than 1.74 Gt (gigatonnes) of finished steel products were produced worldwide [1]. The steel industry is energy-intensive, and most of its growth has been based on fossil fuels, resulting in 1.92 t of carbon dioxide (CO2) emissions per ton of steel produced [2]. However, as governments and businesses commit to a net-zero horizon that requires avoiding most emissions of climate-change-inducing greenhouse gases, it is apparent that new energy systems are required to power the production and processing of steel.
High-temperature process heat is among the most challenging industrial end-uses to decarbonize because it combines high thermal power density, tight temperature–time specifications, and product-quality sensitivities that depend on furnace atmosphere and heat-transfer factors. This is particularly true for steel reheating furnaces in hot rolling mills (e.g., walking-beam and pusher furnaces), which consume 2.0~2.4 GJ/t of finished product [3], and where slabs are typically heated to > 1100 °C (often ~1250 °C) with strict uniformity requirements to ensure stable rolling, avoid defects, and minimize yield losses through oxidation-related scaling [4]. Data compiled by the European Commission highlights that this downstream step remains a nontrivial fraction of total CO2 emissions at around 0.13 tCO2 per ton of product (cf. 1.29 tCO2/t for blast furnaces), and it motivates the need to evaluate viable decarbonization alternatives at the process level, not only at the level of upstream steelmaking routes [5].
These high temperatures make it infeasible to decarbonize reheating furnaces through technologies like high-temperature heat pumps (HTHPs), which are commonly proposed in other energy-intensive, but lower-temperature processes [4]. Recent literature converges on two primary pathways for reheating-furnace decarbonization: (i) direct electrification and (ii) fuel substitution with green hydrogen (i.e. hydrogen produced from renewable electricity) [6]. A third pathway, (iii) a closed carbon loop with carbon capture and synthetic natural gas (SNG) (produced from captured CO2 and green H2 via methanation) offers process compatibility with existing furnace infrastructure but adds conversion steps and auxiliary units.
However, hydrogen shifts the decarbonization challenge upstream: the overall emissions benefit depends on the electricity mix and the availability of renewable power for electrolysis, purchase of externally produced hydrogen (via truck or pipeline), or both. In a detailed hot-rolling reheating-furnace case study comparing natural gas, electrification, hydrogen/air, and hydrogen/oxygen concepts under different European grid mixes, the CO₂ reduction potential was shown to be highly sensitive to the carbon intensity of electricity supply and the pace of grid greening [6]. As a result, hydrogen-based decarbonization cannot be assessed independently of power-system boundary conditions - particularly in contexts where industrial electrification competes for limited renewable electricity [4], [9].
Fuel substitution of hydrogen in place of natural gas causes changes in the reheating furnace atmosphere, with a notable reduction in CO2 and increase in water vapour (H2O) [10]. Although industrial-scale tests of hydrogen-combustion atmospheres are still scarce, several laboratory-scale studies have been carried out to investigate the potential effects that the different chemical composition and heat transfer characteristics may have on product quality. Scaling is a major product quality concern in hot rolled products as it causes surface defects, reduces the process yield and complicates downstream surface treatments [11]. A recent study found a moderate increase in scaling in the hydrogen-air atmosphere [12]. An earlier study showed that the differences in scaling are dependent on the specific grade of steel being reheated, however [13]. Although experiments are ongoing, a CO2-poor atmosphere could also induce decarburization issues, with implications for surface hardness and long-term fatigue [14]. The combination of changes in scaling and decarburization effects implies that a pure hydrogen fuel might not be a one-size-fits-all solution and that other alternatives should be investigated for at least some steel grades.
Nevertheless, these loops add substantial CAPEX and OPEX (capture, methanation and additional hydrogen production), and they incur an efficiency penalty through additional conversion steps. In a dedicated MILP optimization of a steel hot rolling mill, SNG produced via CO2 capture and methanation was consistently more expensive than direct hydrogen across electricity-price scenarios, even though SNG offers storage and infrastructure compatibility advantages [17]. Still, its high compatibility with existing production processes makes it a worthwhile option to consider in feasibility studies.
Because technology choice depends on discrete design decisions and time-dependent operation, recent work increasingly uses optimization-based approaches to compare decarbonization pathways under realistic constraints.
Mixed-integer linear programming (MILP) enables joint optimization of (i) discrete technology selection, (ii) capacity sizing, and (iii) time-resolved dispatch under variable production demands, energy prices and infrastructure constraints. In the hot rolling context, a recent MILP study comparing on-site hydrogen and on-site SNG under dynamic electricity pricing demonstrates that exploiting price variability via flexible operation and storage can substantially reduce levelized costs of energy (LCoE) relative to flat-price assumptions; it also shows that optimal electrolyser selection shifts with electricity-price level due to efficiency–CAPEX trade-offs [17].
Studies from other sectors that also require high-temperature process heat show that deep decarbonization is often achieved via hybrid solutions, and that the last step toward near-zero emissions can be disproportionately expensive. In glass-furnace techno-economic optimization, partial electrification combined with on-site renewables and storage can achieve moderate abatement at modest cost increases, whereas pushing toward very high abatement levels requires substantially larger investments in storage or hydrogen infrastructure [18]. Multi-stage MILP optimization in cement similarly indicates that combining multiple measures (e.g., alternative energy supply, capture, and integration) can deliver major emissions reductions while highlighting the importance of phased investments and policy incentives for economic viability [19]. Cluster-level MILP studies in high-temperature chemical sectors show that the availability of CO2 transport/storage and power infrastructure can dominate outcomes and that integration across units can reduce costs and emissions compared to isolated decarbonization measures [20]. Collectively, these findings support the expectation that hot rolling mills will benefit from evaluating portfolios of decarbonization options and from explicitly representing time-varying operation and infrastructure constraints, in order to best quantify the business impact of deep decarbonization targets.
Much of the peer-reviewed literature for decarbonization in the iron and steel sector is dedicated to the production of raw steel, while the studies on downstream processing of steel products remain comparatively sparse. The most relevant techno-economic assessment study is the doctoral dissertation by Löffler [21], recently expanded upon in a paper by Zabik, Birkelbach et al. [17]. They apply a MILP methodology to optimize the technology choice, capacity sizing and operation of a decarbonized energy supply system for a hot rolling mill in Austria, and carry out a detailed comparison of the LCoE of a hydrogen- vs SNG-based system under a variety of different energy price scenarios. There are however three additional elements not explored in their case study, and which our present work includes to provide a fuller understanding. These are namely different electrical grid capacity scenarios, which influence the planning of major investments; inclusion of onsite PV and wind power generation in the mix; explicit modelling of carbon capture subsystems; and hybrid solutions combining the advantages of more than one technology.
On the one hand, the strong increase in electricity demand expected from onsite production of green hydrogen calls for an evaluation of different grid connection regimes, especially given the structural bottlenecks that the Austrian power grid faces in the near future [22]. If a substantial part of the additional demand is satisfied through onsite renewables, grid feed-in dynamic of excess power in times of high generation and/or low process demand also need to be considered. The development of plans for investments in the grid is in turn influenced by the demands expected from electricity consumers in general and energy-intensive industries in particular. Consequently, the present work can provide insights into the interrelation of industrial planning and grid conditions and thus be helpful to decisionmakers on both sides.
On the other hand, the previous study only considers electricity supply from the general grid. The present work also introduces the possibility of building solar PV and wind power generation connected directly to the plant, which potentially decreases dependency on the grid but also requires consideration of seasonal variations in generation potential which do not necessarily correlate to variations in process heat demand. The optimal electrical mix for each scenario is calculated trough a representative period-based approach.
In addition, their study simplifies the carbon looping subsystem as a per-unit carbon capture cost; an explicit modelling of the available technologies for carbon capture would enable an optimization of the design of the carbon capture subsystem as well as of the heat integration between its energy heat demands and the waste heat from the furnaces and the methanation reaction on the other.
Finally, their model presents a conventional binary choice between electrolysis technologies, while this work allows for more general (albeit potentially less feasible in practice) hybrid solutions in the electrolysis, carbon capture and methanation subsystems that implement more than one technology in order to take advantage of each one’s relative strengths.
By examining the optimal decarbonized system under different grid connection regimes, explicitly modelling and optimizing the closed carbon loop subsystem and allowing for multi-technology solutions while simplifying the energy price conditions, the present work aims to complement previous case studies and provide a fuller picture of the challenges and opportunities the steel processing industry faces on its road to the net-zero economy.
In this section, a MILP model of a steel hot rolling mill in Austria is developed as a collection of subsystems, with the combined superstructure shown in Figure 1. The objective function and the constraints are defined as linear combinations of non-negative real and binary decision variables, which model both the design and operation of the decarbonized energy system. The superstructure allows for hybrid solutions, possibly combining different technologies in the same subsystem.
The most general superstructure of the system, showing the furnaces, the utility connections and the energy conversion units, as well as all the relevant flows of energy and mass, as described through the “Model and Methods” section. Red: the six furnaces with their fuel inputs, energy outputs and carbon outputs. Blue: the power subsystem consisting of grid connections (purchase and feed-in), onsite renewables (hydro, wind and/or solar PV) and a possible battery storage. Green: the hydrogen subsystem, which exists both in the H2 and CC+SNG routes, and consists of the external connections (purchase and feed-in), the available electrolysers (AWE/PEM and/or SOEC), and a possible H2 tank storage. Gray: the carbon capture subsystem, consisting of the available CC technologies (AS, ASHP and PSA) for both the A/B and C groups of furnaces, internal connections for CO2 utilization, and accounting of direct emissions. Orange and yellow: the methanation subsystem with the two available reactor technologies (OSM and/or TSM), as well as the natural gas and biomethane purchase connections. The latter two subsystems are only needed in the CC+SNG route.
The MILP approach was chosen because it is a well-established method for energy system expansion planning and allows the simultaneous optimization of design and operation. It remains computationally efficient even for complex system configurations and guarantees a global optimum as long as the problem is formulated linearly. This makes it particularly suitable for system-level planning studies. Its main limitation is that nonlinear effects, especially part-load behaviour of individual technologies, can only be represented through piecewise-linear approximations. Moreover, nonconvex functions would require introduction of additional binary variables and thus would make the model much harder to trace. Thus, in this work we aimed for simple linear constraints and refrained from piecewise approximations.
The available routes for decarbonization were narrowed down by a preliminary feasibility review. The industrial partner ruled out direct electrification for a variety of reasons: resistive heating might pose difficulties with achieving the desired temperatures in the existing furnace design due to low heating density; the thickness of some pieces makes induction unsuitable for this case study; and plasma burners are not at a sufficiently high TRL yet. Furthermore, the internal tests on the product quality impact of the hydrogen combustion atmosphere are still inconclusive, and therefore the potentially more expensive, but production-wise less problematic CC+SNG route is also investigated as an alternative to the pure-H2 route.
For computational tractability, the whole year is modelled by taking a randomly selected sample 𝒫 of 22 days 𝑝 (one-sixteenth of the year), each composed of 24 one hour (1 h) time steps t:
An analysis of the trade-off between computational tractability and robustness of results was conducted (see Appendix) and yielded the chosen sample size of 22 days as a number that produces < 5% uncertainty in the objective function value across 13 realizations while keeping computation time low enough to enable quick iteration over different scenarios and parameter values. The choice of 1-hour steps, meanwhile, is based on the time resolution of the data for process heat and electricity demand that the industrial partner made available for the study.
The plant consists of several reheating furnaces, of which six are considered in this study. They will be referred to as furnaces furn = A1, A2, B1, B2, C1 and C2. Furnace A1 is used for reheating steel billets from room temperature to forming temperature, and furnace A2 for keeping them hot after rolling. Subsequently, the steel pieces are further processed in one of two heat treatment lines: either at furnaces B1 and B2, which are co-located with furnaces A1 and A2; or at furnaces C1 and C2, which are located in a separate building.
Time series of temperature and load were provided by the plant for all furnaces, with a duration of one year and resolution of one hour or better. These data were used in the model but are not shown in detail here for anonymization reasons; instead a summary is provided in Table 2.
Statistics of process heat demands. The reference load MWq is the average thermal demand of the sum of all the furnaces.
Furnace |
A1 |
A2 |
B1 |
B2 |
C1 |
C2 |
Total |
|---|---|---|---|---|---|---|---|
Peak load [MW/MWq] |
1.28 |
0.25 |
0.37 |
0.19 |
0.19 |
0.07 |
1.93 |
Mean load [MW/MWq] |
0.62 |
0.06 |
0.17 |
0.07 |
0.05 |
0.02 |
1.00 |
Mean/peak load [%] |
49 |
24 |
46 |
39 |
28 |
28 |
52 |
σ [MW/MWq] |
0.32 |
0.04 |
0.12 |
0.05 |
0.05 |
0.05 |
0.47 |
σ / Mean load [%] |
51 |
73 |
69 |
64 |
91 |
87 |
47 |
The thermal energy that goes into the furnace at each time-step is modelled as the sum of the Lower Heating Values (LHV) of the fuels
where
with the furnace-specific
The efficiency of each furnace when burning natural gas was inferred from yearly-averaged energy flow data including fuel inputs, product heat, flue gas heat and cooling loads. For simplicity, the efficiencies for other fuels were extrapolated by taking into account the different heating values and stoichiometries but not possible changes in heat transfer dynamics.
Parametrization of furnaces in the MILP model
Furnace |
Prod. eff. |
Fuel efficiency |
Waste heat >150ºC |
||||
|---|---|---|---|---|---|---|---|
A1 |
0.660 |
0.868 |
0.886 |
0.869 |
0.112 |
0.112 |
0.112 |
A2 |
0.440 |
0.758 |
0.773 |
0.759 |
0 |
0 |
0 |
B1 |
0.733 |
0.833 |
0.850 |
0.834 |
0 |
0 |
0 |
B2 |
0.500 |
0.762 |
0.777 |
0.763 |
0 |
0 |
0 |
C1 |
0.708 |
0.837 |
0.854 |
0.838 |
0.031 |
0.031 |
0.031 |
C2 |
0.575 |
0.870 |
0.887 |
0.871 |
0.029 |
0.029 |
0.029 |
In the following the external connections available to the plant are described. These comprise:
The OPEX associated with gas and biomethane are:
According to data from the International Energy Agency, during the decade and a half from 2010 to 2024 the average spot price for electricity in Austria was 76.47 EUR/MWh in nominal terms [26]. As detailed forecasting of future energy market dynamics falls outside the scope of this work, for the purposes of the MILP model a constant nominal price of
As of 2026, the current price of CO2 emissions certificates in Austria is set by law at 55 EUR/t [29]. However, from 2027 on, the prices will be set dynamically in the EU-wide emissions trading system, with a recent forecast giving a predicted price of 145 EUR/ton by 2030 in the baseline scenario [30]. Thus, as a middle point for the near future, it seems reasonable to set the model’s emissions price at
Parametrization of utility prices in the MILP model
Utility |
Parameter |
Symbol |
Value |
Reference |
|---|---|---|---|---|
Natural gas |
Capacity |
Unlimited |
see text |
|
Price (buy) |
40 EUR/MWh |
[23] |
||
Biomethane |
Capacity |
Unlimited |
see text |
|
Price (buy) |
120 EUR/MWh |
|||
Emissions factor |
2.743 tCO2/t |
see text |
||
Electricity |
Capacity |
Varies |
scenario-dependent |
|
Price (draw) |
80 EUR/MWh |
[26] |
||
Price (feed-in) |
-55 EUR/MWh |
see text |
||
H2 |
Capacity |
Unlimited |
see text |
|
Price (draw) |
150 EUR/MWh |
|||
Price (feed-in) |
-60 EUR/MWh |
see text |
||
CO2 emissions |
Cost of emissions |
100 EUR/tCO2 |
For utility-scale PV, the Fraunhofer report cites a CAPEX range of 700~900 EUR/kWp, which amounts to 140~180 EUR/m2 under 20% efficiency and the standard 1 kW/m2 [32]. Such value for the efficiency is in the 18~25% range observed for commercial PV modules in another Fraunhofer study [34]. Here a more conservative value of
The investment cost for renewable generation is:
Statistics of renewable generation profiles
Technology |
Parameter |
Symbol |
Peak |
Mean |
σ |
H when =0 |
|---|---|---|---|---|---|---|
Wind |
Load factor |
1.00 |
0.27 |
0.30 |
15 % |
|
Solar PV |
Irradiance |
886 W/m2 |
133 W/m2 |
213 W/m2 |
47 % |
|
Hydro |
Load factor |
1.00 |
0.44 |
0.23 |
0 % |
Parametrization of renewable energy sources in the MILP model
Technology |
Parameter |
Symbol |
Value |
Reference |
|---|---|---|---|---|
Wind |
Specific CAPEX |
3,614 kEUR/MW |
60 EUR/MWh [32] |
|
PV |
Efficiency |
20 % |
[34] |
|
Maximum area |
717 m2/MWq |
Site measurements |
||
Specific CAPEX |
250 EUR/m2 |
Computed from [32] |
Green hydrogen is produced by electrolysis of purified water using green electricity, either directly from renewable power or with power drawn from the electrical grid with the appropriate green generation certificates. In some cases, a heat input is also required to maintain the temperature of the reaction.
In general, the electrolyser stack is modelled as a linear relation between the power input
where
There is a maximum fraction
and also a technology-dependent minimal part load, such that if the nominal capacity is
The logic of (16) is implemented in the MILP framework by using a binary time series commitment variable paired with the appropriate constraints.
If
Parametrization of electrolysis technologies in the MILP model
Technology |
Parameter |
Symbol |
Value |
Reference |
|---|---|---|---|---|
AWE/PEM |
Efficiency |
62 % |
Calculated from [36] |
|
Minimal load |
5 % |
[39] |
||
Maximal heat ratio |
0 % |
[39] |
||
Specific CAPEX |
850 kEUR/MW |
Industry source, supported by [36] |
||
SOEC |
Efficiency |
75 % |
[40] |
|
Minimal load |
50 % |
[40] |
||
Maximal heat ratio |
20 % |
[40] |
||
Specific CAPEX |
3,750 kEUR/MW |
In the reheating furnace context, post-combustion carbon capture consists generally of collecting the flue gas of furnaces and passing it through a device that separates the CO2 from the rest of the flue gas, so that it can be either stored away or utilized in some other process. In this case, the purpose of the CC is to produce SNG that is burned again in the furnaces, in a quasi-closed loop. The carbon capture process requires electricity and/or heat to separate the CO2 from the rest of the species in the flue gas, with the specific values varying by technology. Assuming constant efficiency over the operating range, the linearised input-output relations for a CC unit with recovery rate
where the subindices "in", "cap" and "slip" refer to the mass flows of CO2 that comes in with the flue gas, that is captured, and that slips from the CC unit without being captured, respectively.
A multitude of different CC technologies exist, divided into several families and at different technology readiness levels (TRL) [41]. In this work, two possibilities from the absorption family and one from the adsorption family are considered.
As the CC units should ideally be located as close as possible to the source of the flue gas, the system is envisioned as two separate parts: a choice of AS, ASHP and/or PSA for the group of furnaces A/B and another for the pair of furnaces C, with access to recovered waste heat shared between all. This is illustrated on the right side of Figure 1. If the nominal capacity of each unit in terms of captured mass per unit of time is
Parametrization of carbon capture technologies in the MILP model
Technology |
Parameter |
Symbol |
Value |
Reference |
|---|---|---|---|---|
AS |
Separation rate |
95.0 % |
[46] |
|
Heat demand |
3.79 MJ/kg |
[46] | ||
Electricity demand |
0.1744 MJ/kg |
[46] | ||
Specific CAPEX |
225 kEUR/(t/h) |
Industry sources | ||
ASHP |
Separation rate |
95.0 % |
[42] | |
Heat demand |
- |
[42] | ||
Electricity demand |
1.62 MJ/kg |
[42] | ||
Specific CAPEX |
1.697 kEUR/(t/h) |
[42] | ||
PSA |
Separation rate |
98.2 % |
[43] | |
Heat demand |
- |
[43] | ||
Electricity demand |
1.7 MJ/kg |
[43] | ||
Specific CAPEX |
1,111 kEUR/(t/h) |
Hydrogen (H2) and carbon dioxide (CO2) can react exothermically to produce methane (CH4). In the ideal reaction, the so-called Sabatier process can be summarised [47] in the following stoichiometric formula:
In a methanation reactor, this reaction can be carried out in a controlled manner to achieve specific compositions and conversion rates, depending on the desired properties of the resulting gas mixture. The fraction of residual carbon dioxide and hydrogen left in the resulting gas depend on the actual stochiometric ratios in the initial mixture and on the respective conversion rates
The investment costs associated with the methanation reactor(s) are:
Parametrization of methanation technologies in the MILP model
Technology |
Parameter |
Symbol |
Value |
Reference |
|---|---|---|---|---|
OSM |
CO2 conversion rate |
94 % |
Industry source |
|
H2 conversion rate |
93 % |
Industry source |
||
LHV of SNG |
44.52 MJ/kg |
Calculated from industry source |
||
Specific CAPEX |
820 kEUR/MW |
Industry source |
||
TSM |
CO2 conversion rate |
100 % |
Industry source |
|
H2 conversion rate |
99 % |
Industry source |
||
LHV of SNG |
50.36 MJ/kg |
Calculated from industry source |
||
Specific CAPEX |
1,172 kEUR/MW |
Industry source |
In case of load and or/cost variations in the energy supply (due to the intermittency of renewable sources) or in the energy demand (due to production requirements), energy storage can provide flexibility by smoothing out those fluctuations. In this work two possibilities have been considered: electro-chemical storage in batteries and chemical storage in H2 tanks. Both are modelled using a linearized approach with a charging/discharging efficiency parameter η and a standing state-of-charge loss λ.
If each storage has a maximum state-of-charge (SoC) capacity
Charging/discharging losses for the battery storage were (somewhat optimistically) considered negligible for the purposes of this study, but auxiliary losses were accounted for at 1%/h as a typical value for multi-MWh-scale, fast-cycling storages [49].The investment cost was set at
Parametrization of energy storage technologies in the MILP model
Technology |
Parameter |
Symbol |
Value |
Reference |
|---|---|---|---|---|
Battery |
SoC loss |
1 %/h |
[49] |
|
Roundtrip eff. |
100 % |
Simplification |
||
Specific CAPEX |
100 kEUR/MWh |
[50] |
||
H2 tank |
SoC loss |
1 %/h |
[49] |
|
Roundtrip eff. |
80 % |
[51] |
||
Specific CAPEX |
400 kEUR/MWh |
Assumption |
Bringing together eqs. (6), (7), (8) and (9) yields the total annual OPEX:
while the total CAPEX are obtained from eqs. (13), (17), (20), (24) and (25):
bringing the total annualized costs to:
with an interest rate i =3% and amortization period n =20 years yielding a CAPEX annualization factor of 6.7%.
It is expected that the grid will be strained in the future due to increased demand both in the plant and in the broader region as both the steel industry and other sectors electrify their processes. Although quantitative forecasts of future grid availability are outside the scope of this work, it was deemed relevant to establish three exploratory scenarios regarding the grid connection regime. The three scenarios respectively represent optimistic, baseline and pessimistic possibilities for the medium-term development of the power grid in the region and in Austria at large.
This section consists of three parts. First, the energy costs for both decarbonization routes are studied, as well as the associated optimal design of the energy system. Secondly, a more detailed look is given at the optimal operation of the electrolysis and CC subsystems. Thirdly, a parameter variation analysis is conducted to identify the influence of some of the major model parameters.
Figure 2 shows a breakdown of the levelized costs of heat (LCoH), in terms of EUR/MWh of product heat. The most apparent result is that in all scenarios, switching to a decarbonized energy system incurs a stark increase in the LCoH. If hydrogen can be burned directly in the furnaces, the LCoH is shown to increase by between 124% (if the electrical grid is unconstrained) and 166% (if it is strongly constrained). If the CC+SNG subsystem is required, the costs increase further: by 187% in the unconstrained grid (28% more expensive than only hydrogen) and by 231% in the strongly constrained grid (24% more expensive), a more than three-fold increase. There are two main factors for this difference: on the one hand, the additional CAPEX and OPEX associated with the CC+SNG subsystem, which are not present in the pure-H2 path; on the other, the additional CAPEX and OPEX induced in the electricity and hydrogen supplies, due to the conversion losses in the methanation step and to the lower combustion efficiency of SNG relative to hydrogen. The results are consistent with the lower end of the 25~43% bracket that Zabik, Birkelbach et al. found for the SNG cost premium [17].
Waterfall plot for the cost breakdown in both the H2 and the SNG path, including expenses (green), incomes (red) and intermediate totals (blue). For each trio of bars, the top series is for the unconstrained grid case; the middle one, for the partially constrained grid; and the one in the bottom, for the strongly constrained grid. The dashed line marks the BAU case cost of 109.80 EUR/MWh, consisting of 73.80 EUR/MWh of natural gas costs and 36.00 EUR/MWh of emissions certificates.
Upon a closer look at the individual grid scenarios, some interesting differences arise, which can be explained by the Sankey diagrams of energy flows in Figs. Figure 3 and Figure 4:
Sankey diagrams of the optimal systems for the H2 route (a: unconstrained grid; b: partially constrained grid; c: strongly constrained grid). All quantities are peak capacities per MW of product heat; the quoted percentages are each unit's average load factor. The diagrams are not to scale with each other; the heat demand is the same in all cases.
Sankey diagrams of the optimal systems for the CC+SNG route (top: unconstrained grid; middle: partially constrained grid; bottom: strongly constrained grid). All quantities are peak capacities per MW of product heat; the quoted percentages are each unit 's average load factor. CO2 slip and emissions compensation from biomethane are not shown. In the CC units, the upper row is for A/B and the lower row for C.
When different options are available whose OPEX and CAPEX form a Pareto front, optimization can show which of them, or which combination of them, offers the lowest total cost under a set of boundary conditions. The optimal operating patterns optimized simultaneously with the capacity decisions shown in Figure 3 and Figure 4 are highlighted in the following.
Optimal operation of the electrolysis subsystem for both routes and under the three grid scenarios, during one of the representative days.
Looking more closely at the three grid scenarios, some extra insights can be gained. When the grid is unconstrained, the electrolyser capacity is dominated by AWE/PEM in the H2 route and by SOEC in the CC+SNG route. This is due to the highly exothermal nature of the methanation process increasing the availability of waste heat, which can be reused in the SOEC unit to further lower its operational costs. Meanwhile, in the strongly constrained scenario, the intermittency of renewable generation causes the AWE/PEM to dominate the optimal solution - a higher capacity of SOEC would require a higher base load supply of electricity, in turn necessitating additional investment in battery storage. The role of storage is in any case quite limited, in line with the static price scenarios from Zabik, Birkelbach et al., and would presumably be more relevant if fluctuating prices were introduced into the model [17].
Optimal operation of the carbon capture subsystem under the three grid scenarios, during the same representative day.
In scenarios with unconstrained grid access, the model favours technologies with higher investment cost but lower operating cost and/or higher efficiency, because sufficient electricity is available to utilize them at high load factors. This leads to steadier operation of units such as SOEC and PSA. By contrast, under constrained grid conditions the model shifts toward more flexible technologies and supply options, including AWE/PEM, hydrogen imports, wind power and storage, because electricity scarcity increases the value of operational flexibility. Thus, the optimal design and the optimal operation should be interpreted jointly: the installed capacities reflect the expected operating regime, while the operating profiles explain why those capacities are economically preferred.
It follows from Figure 2 that the CAPEX of wind power, the cost of purchasing green hydrogen from the gas grid and the CAPEX of the electrolysers are some of the largest contributors to the total energy costs. However, there is considerable uncertainty in the values of those parameters, which may cause significant deviations in the actual energy costs. Therefore, a sensitivity analysis was conducted to study if and to what extent the possible changes in those parameters affect the LCoH.
Considering all the possible combinations of parameter values, decarbonization routes and grid capacities, a total of 270 separate optimization runs were evaluated. It should be noted that, although the parameter variations were treated independently in this study, some of these varied parameters may be correlated in practice. In particular, the cost of grid-purchased hydrogen can partly be shaped by upstream electricity prices, electrolyser investment costs and wind farm CAPEX. In the present study, these parameters were varied independently in order to isolate their directional influence on the optimal solution. The main results, again using LCoH as the main metric, are shown in Figure 7.
Distribution of LCoH generated by the 45 variations, for both routes (H2 in blue, CC+SNG in orange). Each dots represents an individual variation. The dotted lines mark the respective baseline cost.
The most obvious result is that, for the price of purchasing hydrogen from the grid, there is a threshold under which the optimal solution consists of forgoing any on-site electrolysis capacity simply supplying all the H2 through the grid, regardless of the values of the other cost parameters. The exception is the strongly constrained grid scenario, where some electrolysis capacity is still required to process excess renewable energy.
The variations in wind power CAPEX do not have a large influence on the LCoH, as they induce a reduction in wind power capacity and an increased electricity purchase from the grid, again with the exception of the strongly constrained grid. In that case, higher investment costs translate into a higher LCoH.
Finally, the variations in the specific investment cost of electrolysers seem not to have a large impact on average, but cause a larger spread in the distribution of the LCoH generated by the variations of the other parameters. However, higher investment costs do translate to an increase in LCoH while they are still below the grid purchase costs. This can be observed especially in the case of an unconstrained electricity grid connection, where increasing electrolyser costs are directly linked to an increase in LCoH.
The sensitivity analysis reveals that the response of the system is not purely gradual, but exhibits threshold-like behaviour. Most importantly, the attractiveness of on-site hydrogen production depends on the interaction between hydrogen purchase price, grid constraint and costs for renewable electricity. When hydrogen from the grid is sufficiently cheap, the model tends to forgo on-site electrolysis in favour of external supply. However, under strongly constrained grid conditions, on-site electrolysis remains valuable even at low hydrogen purchase prices because it provides a sink for excess renewable electricity. Hence, grid constraints determine whether the system shifts towards a purchase-based or production-based hydrogen supply strategy.
Expressed in absolute terms, the modelled LCoH rises from 109.8 EUR/MWh in the BAU case to approximately 246–292 EUR/MWh for the H2 route and 315–363 EUR/MWh for the CC+SNG route, depending on grid conditions. If product quality and heating process considerations allow for a mix of pure hydrogen and SNG, then the LCoH for the H2 and the CC+SNG pathway can be interpreted respectively as the lower and upper ends of the cost band for such a hybrid solution.
From an industrial perspective, these increases in levelized heat costs are not merely percentage differences but cost levels that would be difficult to absorb in a highly competitive steel market if borne by the plant alone. In other words, even the most techno-economically favourable decarbonized configurations identified here (i.e. pure-hydrogen fuel without major electricity grid constraints) do not appear to be business cases driven by direct heat-cost savings; rather, their implementation would likely depend on sufficiently strong policy support measures such as carbon pricing, CAPEX financing, OPEX subsidies, contracts for difference, or green-steel price premiums.
The results should therefore be interpreted primarily as showing the relative economic attractiveness of the modelled pathways under different infrastructure constraints, while the broader investment decision remains contingent on whether these absolute cost levels can be recovered in the market or offset by policy support. This is particularly relevant because industrial decision-makers do not choose only between hydrogen and CC+SNG, but also between these routes, partial decarbonization steps, efficiency measures, and in other cases electrification-based concepts that were not investigated in this study.
In this work, a model of a steel hot rolling mill in Austria was developed, with (i) direct combustion of green hydrogen and (ii) post-combustion carbon capture with methanation as the possible routes for decarbonizing its energy supply, and several technologies available for each subsystem. The design and operation of the energy system were simultaneously optimized using a mixed integer linear problem (MILP) formulation. Even after optimizing for minimal annualized costs, a major increase in the levelized cost of heat (LCoH) was predicted.
However, the increase in LCoH was found to be noticeably dependent on the constraints of the system. If direct combustion of green hydrogen is allowed, the cost is predicted to increase to 2.5 - 2.9 times the BAU LCoH; if synthetic natural gas (SNG) is required, then the increase is 3.2 - 3.6 times the BAU LCoH, mainly due to the extra investment costs and conversion losses of the carbon capture and methanation steps. The strictness of the limitations on the electrical grid to provide flexibility not only has a large impact on the headline LCoH, but also on the underlying configuration of the energy system.
Finally, the parameter variation analysis showed that the capital costs of wind power generation and the purchase costs of hydrogen from the gas grid also have a large influence on the LCoH and the configuration of the system, with the investment costs of electrolysers having a comparatively minor influence. Together, the results highlight the importance of accurate modelling of the boundary conditions when predicting the costs of the energy transition, but show that optimization methods can help in revealing the cost-optimal strategies for minimizing said costs.
From a business perspective, the absolute LCoH levels suggest that, under the assumptions of this study, neither route is likely to be adopted on the basis of direct cost competitiveness alone. Rather, deployment would depend on favourable infrastructure conditions together with policy support and/or market mechanisms that reward low-carbon steel production.
The insights obtained from the present work suggest a range of refinements of the model that could prove to be interesting to investigate in future research. On the one hand, (i) the outsize impact of grid boundary conditions on the design and operation of the system calls for a more detailed modelling of the grid dynamics and the electricity and hydrogen markets, including time-dependent prices and capacities, with price fluctuations in particular having a potentially major effect on the sizing and operation of energy storages. On the other hand, (ii) more detailed modelling of operational conditions of the technologies involved in the electrolysis, carbon capture and methanation subsystems could prove useful in the context of broader feasibility studies, with more emphasis on the efficiency changes over the operating range of each technology, as well as more realistic models of energy storage technologies. Finally, the model could be expanded with (iii) additional alternative fuel possibilities (ammonia, methanol, enrichment of combustion air with oxygen produced as by-product of electrolysis) as well as with (iv) the search for optimal decarbonization pathways, i.e. intermediate solutions with partial decarbonization steps and the associated optimization of the timing of investment decisions.
This project is supported with the funds from the Climate and Energy Fund and implemented in the framework of the RTI-initiative “Flagship region Energy”.
- SteelFacts
- World Steel in Figures 2025
- Tracking industrial energy efficiency and CO2 emissions: In support of the G8 plan of action
- ,
Electrification or Hydrogen? The Challenge of Decarbonizing Industrial (High-Temperature) Process Heat ,Multidisciplinary Scientific Journal , Vol. 7 (4), :439-456, https://doi.org/10.3390/j7040026 - Prospective scenarios on energy efficiency and CO2 emissions in the EU iron & steel industry: Re-edition
- ,
Towards CO2-neutral process heat generation for continuous reheating furnaces in steel hot rolling mills – A case study ,Energy , Vol. 224 , , https://doi.org/10.1016/j.energy.2021.120155 - , Electrification of reheating furnaces: state of the art and future research needs
- ,
Rolling mill decarbonization: Tenova SmartBurners with 100% hydrogen ,Matériaux & Techniques , Vol. 109 (3-4), , https://doi.org/10.1051/mattech/2022012 - ,
A Review of the Role of Hydrogen in the Heat Decarbonization of Future Energy Systems: Insights and Perspectives ,Energies , Vol. 17 (7), , https://doi.org/10.3390/en17071688 - ,
Utilization of hydrogen fuel in reheating furnace and its effect on oxide scale formation of low-carbon steels ,International Journal of Hydrogen Energy , Vol. 140 , :1212-1220, https://doi.org/10.1016/j.ijhydene.2024.11.230 - ,
Scale Formation Minimization on Internal Surface of Seamless Hot-Rolled Pipes ,Steel Transl. , Vol. 50 (5), :340-346, https://doi.org/10.3103/S0967091220050046 - ,
Oxidation of Carbon Steels in Novel Reheating Conditions: Changes to Oxidation Kinetics ,Metall Mater Trans B , Vol. 56 (4), :3762-3773, https://doi.org/10.1007/s11663-025-03588-7 - ,
Oxide Scale Formation on Low-Carbon Steels in Future Reheating Conditions ,Metals , Vol. 14 (2), , https://doi.org/10.3390/met14020189 - , Investigation of Decarburization of Spring Steel During Hot Rolling
- ,
Demonstration of a high-efficiency carbon dioxide capture and methanation system with heat/material integration for power-to-gas and zero carbon dioxide emissions in flue gasses ,International Journal of Greenhouse Gas Control , Vol. 114 , , https://doi.org/10.1016/j.ijggc.2022.103584 - ,
Synthetic natural gas via integrated high-temperature electrolysis and methanation: Part I—Energy performance ,Journal of Energy Storage , Vol. 1 , :22-37, https://doi.org/10.1016/j.est.2015.04.002 - ,
Hydrogen and Synthetic Natural Gas for decarbonizing steel hot rolling mills: Economic viability under dynamic electricity pricing ,Energy Conversion and Management , Vol. 341 , , https://doi.org/10.1016/j.enconman.2025.119912 - ,
Techno-Economic Assessment of Electrification and Hydrogen Pathways for Optimal Solar Integration in the Glass Industry ,Solar , Vol. 5 (3), , https://doi.org/10.3390/solar5030035 - ,
Modeling and Multi-Stage Planning of Cement-IIES Considering Carbon-Green Certificate Trading ,Processes , Vol. 11 (4), , https://doi.org/10.3390/pr11041219 - ,
Optimizing Emissions Reduction in Ammonia–Ethylene Chemical Clusters: Synergistic Integration of Electrification, Carbon Capture, and Hydrogen ,Ind. Eng. Chem. Res. , Vol. 64 (8), :4479-4497, https://doi.org/10.1021/acs.iecr.4c03817 - , , Transformation of a Steel Rolling Plant: Identification of Decarbonisation Pathways using MILP
- , APG’s Future Grid Puts Supply-Secure Energy Transition on the Fast Track, Austrian Power Grid
- , Energy prices, taxes, Statistik Austria
- Marktbericht 2023 (accessed )
- Potenzialstudie Grüne Gase: Analyse und Bewertung der Potenziale Grüner Gase in der Innovationsregion Mitteldeutschland (accessed )
- , Real-Time Electricity Tracker, International Energy Agency
- Importmöglichkeiten für erneuerbaren Wasserstoff (accessed )
- , , Ministerialentwurf Erneuerbares-Gas-Gesetz
- , Emissionspreisindex, Austrian Energy Agency
- , Scenarios shaping EU ETS prices
- IEA Wind TCP Austria 2022
- Study: Levelized Cost of Electricity- Renewable Energy Technologies
- National Survey Report of PV Power Applications in Austria 2023 (accessed )
- Photovoltaics Report (accessed )
- ,
Assessment of the three most developed water electrolysis technologies: Alkaline Water Electrolysis, Proton Exchange Membrane and Solid-Oxide Electrolysis ,Materials Today: Proceedings , Vol. 66 , :140-145, https://doi.org/10.1016/j.matpr.2022.04.264 - Cost Forecast for Low-Temperature Electrolysis - Technology Driven Bottom-Up Prognosis for PEM and Alkaline Water Electrolysis Systems
- Electrolyzers for Hydrogen Production: Technical and Economic Characteristics Executive Summary
- Next Level Solid Oxide Electrolysis: Upscaling potential and techno-economical evaluation for 3 industrial use cases
- , Silyzer 300: The next paradigm of PEM electrolysis, Siemens Energy
- , Sunfire-HyLink SOEC, Sunfire SE
- ,
Technology readiness level assessment of carbon capture and storage technologies ,Renewable and Sustainable Energy Reviews , Vol. 200 , , https://doi.org/10.1016/j.rser.2024.114578 - , Improving Energy Efficiency of Carbon Capture Processes with Heat Pumps
- ,
Overview on Pressure Swing Adsorption (PSA) as CO2 Capture Technology: State-of-the-Art, Limits and Potentials ,Energy Procedia , Vol. 114 , :2390-2400, https://doi.org/10.1016/j.egypro.2017.03.1385 - ,
Energy and cost estimates for capturing CO2 from a dry flue gas using pressure/vacuum swing adsorption ,Chemical Engineering Research and Design , Vol. 102 , :354-367, https://doi.org/10.1016/j.cherd.2015.06.033 - ,
How much can novel solid sorbents reduce the cost of post-combustion CO2 capture? A techno-economic investigation on the cost limits of pressure–vacuum swing adsorption ,Applied Energy , Vol. 306 , , https://doi.org/10.1016/j.apenergy.2021.117955 - ,
Techno Economic Evaluation of Amine based CO2 Capture: Impact of CO2 Concentration and Steam Supply ,Energy Procedia , Vol. 23 , :381-390, https://doi.org/10.1016/j.egypro.2012.06.053 - ,
Review on methanation – From fundamentals to current projects ,Fuel , Vol. 166 , :276-296, https://doi.org/10.1016/j.fuel.2015.10.111 - ,
Advancements in CO2 methanation: A comprehensive review of catalysis, reactor design and process optimization ,Chemical Engineering Research and Design , Vol. 201 , :457-482, https://doi.org/10.1016/j.cherd.2023.11.060 - , 5 MWh Battery Energy Storage System for North America, Chint Power Systems
- Battery technology in the European Union: 2024 Status report on technology development, trends, value chains and markets
- DOE Hydrogen and Fuel Cells Program Record 9013: Energy requirements for hydrogen gas compression and liquefaction as related to vehicle storage needs

