Techno-Economic Assessment of Green Hydrogen, Carbon Capture and Synthetic Natural Gas as Decarbonization Options for a Steel Hot Rolling Mill in Austria

Original scientific paper

Journal of Sustainable Development of Smart Energy Networks
Volume 2, Issue 1, pp 1-30
DOI: https://doi.org/10.13044/j.sdi.d3.0733 (registered soon)
Mikel Olaciregui-Segura1 , Anton Beck2, Christoph Zauner2
1 AIT Austrian Institute of Technology GmbH, WIen, Austria
2 AIT Austrian Institute of Technology GmbH, Wien, Austria

Abstract

Hot rolling mills for steel processing are an energy-intensive and notably hard-to-abate industry. Previous studies have not explicitly modelled the interaction between the different technologies that potentially make up the decarbonized energy supply system of a steel rolling mill, nor considered potential limitations of electrical grid infrastructure. The combination of additional capital and operational expenses together with infrastructural limitations potentially makes deep decarbonization of such an industrial energy system economically unfeasible. In this paper, we present a complete model of a hot rolling mill in Austria with two decarbonization routes, namely green hydrogen and synthetic natural gas with carbon capture, and we use Mixed-Integer Linear Programming optimization to derive the least-cost design and operation of the energy system under three power grid capacity scenarios and a choice of different technologies. We find a 124~166% increase in the levelized cost of heat for the hydrogen route and of 187~231% for the synthetic natural gas route, depending on power grid capacity. The results highlight the crucial role of the electrical grid in providing flexibility and confirm that, under current and near-future techno-economic conditions, a major reduction in carbon emissions from the steel processing sector is only possible with significant policy support.

Keywords: Steel rolling, Reheating furnace, Decarbonization, Green hydrogen, Carbon capture, Synthetic natural gas, Mixed-integer linear optimization

Creative Commons License
Views (in 2026): 227 | Downloads (in 2026): 246
Total views: 227 | Total downloads: 246
Introduction

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 in Hot Rolling

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].

Decarbonization Options for Steel Reheating and Hot Rolling

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.

Direct electrification. Direct electrification eliminates combustion emissions at the point of use and can be highly efficient, but its application to large, continuous reheating furnaces is subject to important practical constraints. Replacing fuel burners with electric heating removes the flue-gas stream that supports heat recovery and affects convective mixing, meaning that furnace flow and heat-integration concepts may need to be substantially redesigned to maintain productivity and temperature uniformity [4]. Major technologies for direct electrification include resistive heaters, induction coils and plasma torches [7]. In the present study, resistive heating was examined but found unsuitable for the larger furnaces because its heat density is insufficient. Induction heating was also excluded, as it is not compatible with the site-specific products due to spatial constraints and their inhomogeneous shapes. Plasma torches, by contrast, offer a burner-like electrification concept with high power density and a hot gaseous heat-transfer medium, making them a promising option for future applications. However, because they currently remain at a relatively low technology readiness level and still face challenges related to efficiency, heat transfer, NOx formation, product quality, and maintenance, they were not included in the present study either.

Fuel substitution with green hydrogen. Hydrogen (H2) combustion is often positioned as a relatively straightforward retrofit option because it can retain established burner–furnace architectures. High-temperature process-heat reviews note that steel reheating furnaces can technically be converted to hydrogen and that pilot/industrial deployments are approaching high readiness levels (often framed as TRL ~8–9), while also emphasizing that feasibility hinges on site constraints and available infrastructure [4]. Industrial burner development for reheating furnaces has demonstrated stable operation up to 100% H2 with low-NOx concepts at typical furnace temperatures, supporting the practical plausibility of hydrogen retrofits [8].

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.

Fuel substitution with carbon capture and Synthetic Natural Gas. An alternative that requires minimal process modifications is to continue using methane gas (CH4) as fuel while making the system near-carbon-neutral via CO2 capture and utilization: CO2 is captured from the flue gas and methanated with green hydrogen to form SNG, which is recycled as furnace fuel. Pilot demonstrations show that high capture and conversion performance are technically achievable, particularly when heat and material integration are leveraged to reduce auxiliary energy demand [15]. Process-integration studies of power-to-SNG chains further quantify how overall efficiency depends strongly on electrolysis type (steam electrolysis vs co-electrolysis), operating pressure, and heat integration between (exothermic) methanation and heat demands [16].

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.

Techno-Economic Assessment Studies Using Mixed-Integer Linear Programming

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.

Research Gaps and Motivation for this Study

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.

Model and Methods

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:

tT=0,,23 pP=31,34,47,71,136,151,153,164,183,185,187,227,245,257,262,268,272,283,300,304,347,350

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.

Industrial Process

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 fFthat are being burned at the time, weighted by a fuel type- and furnace-dependent efficiency parameter ηffurn that accounts for the heat that is lost in the flue gas:

qt,pfurn=fFηffurn×LHVf×m˙t,pf,furn

where m˙t,pf,furn is the mass flow of each fuel being burned at each time-step. The heat that goes into the product itself is calculated as:

qpr,t,pfurn=ηprfurn×qt,pfurn

with the furnace-specific ηprfurn  accounting for heat lost via radiation, cooling, etc. Similarly, the energy that can be recovered as high temperature waste heat (defined here as temperatures higher that 150 °C) is given by another furnace- and fuel-specific parameter r>150furn:

wt,pfurn=fFr>150furn×ηffurn×LHVf×m˙t,pf,furn

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

ηprfurn

ηNGfurn

ηH2furn

ηlgSNGfurn

rNGfurn

rH2furn

rlgSNGfurn

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

Utilities

In the following the external connections available to the plant are described. These comprise:

Gas grid. Most of the energy demand in the business-as-usual (BAU) case is supplied via natural gas pipeline. Since the gas prices are specific to every mill and usually extremely confidential, in this study a fixed price cbuyNG  = 40 EUR/MWh is used, based on the historical average of 38.97 EUR/MWh (in nominal terms) between 2010 and 2024 [23]. The capacity of the gas pipeline used for supplying the mill is not a relevant constraint, as it covers all the current demand and, in any decarbonization scenario, the use of natural gas is unlikely to increase significantly – quite the opposite.

Biomethane. A small amount of carbon emissions, caused by carbon slip in the CC subsystem, may be offset by purchasing biomethane at a price cbuyBM, which is has an emissions factor efBM of embedded biogenic carbon. This quantity is approximated as the stoichiometric carbon content of the fuel, such that a zero emissions solution is feasible. The capacity of the biomethane delivery system (be it pipeline or trucks) is not considered a significant constraint, as the demand of biomethane allowed in the decarbonized system is considerably lower than the current demand of natural gas. An Austrian government survey of several existing biomethane production sites quotes a production cost between 102 and 182 EUR/MWh [24]. Additionally, a German study on development potential for green gases puts the costs of waste-sourced biomethane at 118 EUR/MWh in the baseline scenario [25]. Thus, in the present model a future price of cbuyBM = 120 EUR/MWh is assumed.

The OPEX associated with gas and biomethane are:

OPEXgas=pPtTcbuyNG×m˙CH4,t,pNG+cbuyBM×m˙CH4,t,pBM

Power grid. The power grid is currently the main electricity supply infrastructure. The cost is approximated as a constant cbuyp, with the option to feed-in excess renewable production to the grid at a price csellp. The operating costs associated to the power grid are:

OPEX power = pP tT c buy p × p t,p buy + c sell p × p t,p sell

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 cbuyp = 80 EUR/MWh is used as a proxy for year-round averages. For excess power, the feed-in income is set at csellp = -55 EUR/MWh. This value is set somewhat below the levelized cost of electricity for wind power, so that the grid is allowed to provide some flexibility in times of high generation without the sale of renewable electricity becoming the main business case by itself.

Hydrogen grid. It is expected that in the future it will be possible to import a volume m˙H2,t,pbuy of green hydrogen through the existing gas pipeline infrastructure at a price cbuyH2. The chosen price reflects a midpoint of the 93~171 EUR/MWh projected for pipeline imports into Austria plus some costs for end-user distribution [27], and also takes into account the future production costs targeted by the penalty-and-investment mechanism envisioned in the Austrian ministry’s draft Renewable Gases Act [28]. The pipeline in the model also allows to export excess green hydrogen at a rate mH2,t,psell and price csellH2 when it is not needed on-site. The income is set at a much lower value of 55 EUR/MWh as a deliberate modelling choice to reflect the assumption that feed-in of hydrogen will happen at times of high supply, and thus low market price, across the region. The hydrogen grid-associated OPEX are thus calculated as:

OPEXH2=pPtTcbuyH2×LHVH2×m˙H2,t,pbuy+csellH2×LHVH2×m˙H2,t,psell

CO2 Emissions. For calculation of the business-as-usual (BAU) scenario, the CO2 emissions certificates need to be taken into account at price cemitCO2:

OPEXCO2=pPtTcemitCO2×m˙CO2,emit,t,pefBM×m˙BM,emit,t,p

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 cemitCO2 = 100 EUR/t.

Parametrization of utility prices in the MILP model

Utility

Parameter

Symbol

Value

Reference

Natural gas

Capacity

m˙CH4,peak

Unlimited

see text

Price (buy)

cbuyNG

40 EUR/MWh

[23]

Biomethane

Capacity

m˙BM,peak

Unlimited

see text

Price (buy)

cbuyBM

120 EUR/MWh

[24], [25]

Emissions factor

efBM

2.743 tCO2/t

see text

Electricity

Capacity

ppeak

Varies

scenario-dependent

Price (draw)

cbuyp

80 EUR/MWh

[26]

Price (feed-in)

csellp

-55 EUR/MWh

see text

H2

Capacity

m˙H2,peak

Unlimited

see text

Price (draw)

cbuyH2

150 EUR/MWh

[27], [28]

Price (feed-in)

csellH2

-60 EUR/MWh

see text

CO2 emissions

Cost of emissions

cemitCO2

100 EUR/tCO2

Interpolation from [29] and [30]

Wind power. Onshore wind has had a rapid development in Austria in recent years, with an average net yearly increase of 215 MW in installed capacity [31]. There are already operational wind parks in the studied mill's region with more planned, which enabled the use of historical generation profile data (Table 4). In this work, it is assumed that onshore wind will be the main avenue for expansion of low-emission electricity supply, with wind parks located close enough to the factory that they can be connected directly without relying on the already constrained electrical grid. The specific CAPEX are set such that with amortization and interest they yield a levelized cost of electricity (LCoE) of 60 EUR/MWh, in line with the Fraunhofer Institute’s 53~68 EUR/MWh estimate of current LCoE for onshore wind in regions with 2500 full load hours [32], which is very close to the 27% average load factor shown in Table 4. If the time-dependent load factor is LFwind, then the supplied electricity is:

pt,pwind=LFt,pwind×ppeakwind

Solar PV. As of 2023 Austria had nearly 6.4 GW of solar PV installed capacity, representing a ten-fold increase over the previous decade [33]. In the studied steel mill, a modest surface area was considered suitable for onsite PV installation, potentially acting as a complement to wind power. The generated electricity is proportional to the solar irradiation IPV through the efficiency ηPV, and the surface area that can be built is capped at amaxPV:

pt,pPV=It,pPV×ηPV×aPV,  aPVamaxPV

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 cinvPV = 250 EUR/m2 is used to account for lower irradiance, land acquisition costs, and taxes.

Hydropower. A small capacity of hydropower is currently installed that supplies part of the baseline electrical demand. However, it was considered not feasible to expand beyond the current capacity. The load factor is LFt,phydro, thus:

pt,phydro=LFt,phydro×ppeakhydro

The investment cost for renewable generation is:

CAPEXren=cinvwind×ppeakwind+cinvPV×aPV

Statistics of renewable generation profiles

Technology

Parameter

Symbol

Peak

Mean

σ

H when =0

Wind

Load factor

LFt,pwind

1.00

0.27

0.30

15 %

Solar PV

Irradiance

It,pPV

886 W/m2

133 W/m2

213 W/m2

47 %

Hydro

Load factor

LFt,phydro

1.00

0.44

0.23

0 %

Parametrization of renewable energy sources in the MILP model

Technology

Parameter

Symbol

Value

Reference

Wind

Specific CAPEX

cinvwind

3,614 kEUR/MW

60 EUR/MWh [32]

PV

Efficiency

ηPV

20 %

[34]

Maximum area

amaxPV

717 m2/MWq

Site measurements

Specific CAPEX

cinvPV

250 EUR/m2

Computed from [32]

Electrolysis Subsystem

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 pt,pelec and heat input qt,pelec on one side and the mass flow output of hydrogen m˙H2,t,pelec  on the other:

m˙H2,t,pelec=ηelecLHVH2×pt,pelec+qt,pelec

where ηelec  is the efficiency of the electrolyser stack and LHVH2 =119.96 MJ/kg is the lower heating value of hydrogen. For the capacity-expansion model used in this work, AWE/PEM and SOEC units are represented with constant efficiency over their admissible operating range. This simplification is consistent with the planning focus of the model, which prioritizes long-term investment, sizing, and infrastructure interactions over detailed stack electrochemistry. This simplified model improves the computational tractability at the cost of slightly overestimating efficiency at low part load.

There is a maximum fraction rqmaxelec of energy that can be supplied directly as heat:

qt,pelecrqmaxelec×pt,pelec+qt,pelec

and also a technology-dependent minimal part load, such that if the nominal capacity is mH2,peakelec  and the minimal load is limminelec, then depending on whether the electrolyser is active:

m˙H2,t,pelec[limminelec×m˙H2,peakelec,m˙H2,peakelec],if electrolyser is running=0,if electrolyser is shut down

The logic of (16) is implemented in the MILP framework by using a binary time series commitment variable paired with the appropriate constraints.

Alkaline Water Electrolysis (AWE) and Proton Exchange Membrane electrolysis (PEM). In AWE, OH- anions move from cathode to anode through a liquid electrolyte, releasing H2 gas at the cathode. In PEM, water is supplied at the anode and H+ cations traverse a polymer membrane to the cathode, forming H2 gas at the cathode. While both types of electrolysis have their own specificities, in the linearized approach of this work they are grouped together due to their generally lower investment cost, lower efficiency, room-temperature operation and lower minimal part load, the latter being relevant when dealing with intermittent renewable energy sources [35]. The investment cost of cinvAWE/PEM = 850 kEUR/MW was quoted by an industry insider source, but it lines up with the upper band of the forecasted range laid out at the beginning of the decade by the Fraunhofer institute, which is a system cost of 450~950 kEUR/MW for AWE and 500~980 kEUR/MW for PEM, with the upper limit corresponding to 5 MW-scale in 2020 and the lower one for 100 MW-scale in 2030 [36]. The constant part load efficiency of 62%, assumed for both technologies (Eq. 14) in line with a 52.8~55.1 kWh/kg system electricity demand [36], is valid for PEM electrolysers in the range of around 40% to 100% load. Below 40% a drop in efficiency can be expected.

Solid Oxide Electrolysis Cell (SOEC). In SOEC, high-temperature steam is supplied at the cathode and O2- anions are transported from cathode to anode through a solid ceramic membrane. This type of technology has a higher investment cost, higher efficiency, heat for steam generation as part of the energy input (which can be supplied as waste heat from other processes) and a more restrictive minimal part load, which complicates its use with intermittent renewable sources [35]. The investment costs of cinvSOEC = 3,750 kEUR/MW were provided by industry insider sources, including not just the electrolysis cells but also heat integration and the associated infrastructure. This lines up reasonably with costs reported in the literature – the World Bank reported 5,000~5,800 kUSD/MW for current pilot projects and a potential for cost reduction to 2,000 kUSD/MW in upscaled systems [37], while another report in the European context projected 4,800 kEUR/MW in the base case for 2020 and 1,200 kEUR/MW for the upscaled system in 2030 [38].

If m˙H2,peakAWE/PEM  and m˙H2,peakSOEC are the nominal capacities selected for the respective technologies, then the investment costs for the electrolysis subsystem are:

CAPEXelec=cinvAWE/PEM×LHVH2×m˙H2,peakAWE/PEM+cinvSOEC×LHVH2×m˙H2,peakSOEC

Parametrization of electrolysis technologies in the MILP model

Technology

Parameter

Symbol

Value

Reference

AWE/PEM

Efficiency

ηAWE/PEM

62 %

Calculated from [36]

Minimal load

limminAWE/PEM

5 %

[39]

Maximal heat ratio

rqmaxAWE/PEM

0 %

[39]

Specific CAPEX

cinvAWE/PEM

850 kEUR/MW

Industry source, supported by [36]

SOEC

Efficiency

ηSOEC

75 %

[40]

Minimal load

limminSOEC

50 %

[40]

Maximal heat ratio

rqmaxSOEC

20 %

[40]

Specific CAPEX

cinvSOEC

3,750 kEUR/MW

Industry source, supported by [37], [38]

Carbon Capture Subsystem

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 RCC, power demand dpCC and heat demand dqCC are:

m˙CO2,cap,t,pCC=RCC×m˙CO2,in,t,pCC,  m˙CO2,slip,t,pCC=1R×m˙CO2,in,t,pCC pt,pCC=dpCC×m˙CO2,in,t,pCC,  qt,pCC=dqCC×m˙CO2,in,t,pCC

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.

Amine Scrubber (AS). This method works by passing the flue gas through an amine bath in a continuous process. The amine is then regenerated and the CO2 stored. The process has a low power demand but higher heat requirements, and is currently at a very high TRL [41]. The investment costs are taken to be cinvAS= 225 kEUR/(t/h), based on quotes from suppliers.

Amine Scrubber with Heat Pump (ASHP). The AS process can be made more efficient by integrating a heat pump for internal heat recovery. This eliminates the need for external heat input and replaces it with a smaller electricity input, at the cost of higher installation costs [42]. Those costs cinvASHP=  1.697 kEUR/(t/h) are based on the costs for the AS itself plus the cost premium for the heat pump provided in the aforementioned paper.

Pressure Swing Adsorption (PSA). This method separates CO2 by feeding flue gas into a pressurized tank and then greatly reducing the pressure of the tank, which causes the CO2 to be adsorbed into a purposefully designed solid coating (typically zeolite-based) while the other components of the gas remain free. This batch process requires no heat input but instead is powered by electricity, and although its TRL for CO2 capture is lower, it is considered in the literature to promise low operation and investment costs [41], [43]. As exact costs figures for this type of system are hard to find in the literature, an estimate of cinvPSA  = 1,111 kEUR/(t/h) has been calculated from two studies that provide partial cost figures [44], [45].

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 CO2,peakCC, then the associated investment costs are:

CAPEXCC=cinvAS×m˙CO2,peakAS,A/B+m˙CO2,peakAS,C+cinvASHP×m˙CO2,peakASHP,A/B+m˙CO2,peakASHP,C+ cinvPSA×m˙CO2,peakPSA,A/B+m˙CO2,peakPSA,C

Parametrization of carbon capture technologies in the MILP model

Technology

Parameter

Symbol

Value

Reference

AS

Separation rate

RAS

95.0 %

[46]

Heat demand

dqAS

3.79 MJ/kg

[46]

Electricity demand

dpAS

0.1744 MJ/kg

[46]

Specific CAPEX

cinvAS

225 kEUR/(t/h)

Industry sources

ASHP

Separation rate

RASHP

95.0 %

[42]

Heat demand

dqASHP

-

[42]

Electricity demand

dpASHP

1.62 MJ/kg

[42]

Specific CAPEX

cinvASHP

1.697 kEUR/(t/h)

[42]

PSA

Separation rate

RPSA

98.2 %

[43]

Heat demand

dqPSA

-

[43]

Electricity demand

dpPSA

1.7 MJ/kg

[43]

Specific CAPEX

cinvPSA

1,111 kEUR/(t/h)

Estimate from [44] and [45]

Synthetic Natural Gas Subsystem

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:

CO2+4H2CH4+2H2O+164.0 kJ/mol

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 χCO2meth and χH2meth. In the linearized model it is assumed that the reactor is operated at a fixed stochiometric ratio of feedstocks and thus at fixed conversion ratios and waste heat generation Then the linearized input-output relations for the reactor model, relating the CO2 and H2 inputs to the CH4, waste heat, CO2 and H2 outputs are, respectively:

m˙CH4,out,t,pmeth=χCO2meth×m˙CO2,in,t,pmeth,  wht,pmeth=10.2 MW/(kg/s)×m˙CH4,out,t,pmeth m˙CO2,out,t,pmeth= 1χCO2meth×m˙CO2,in,t,pmeth,    m˙H2,out,t,pmeth= 1χH2meth×m˙H2,in,t,pmeth

One-Stage Methanation (OSM). In a conventional methanation reactor, not all the CO2 is converted into methane; in addition to a large concentration of CH4, also some amount of residual H2 and CO2 are present. This changes the LHV and flame properties of the fuel with respect to the near-pure methane of natural gas; consequently in this work said mixture is referred to as low-grade synthetic natural gas (lgSNG). All the parameters for this type of reactor were collected from an industry supplier source.

Two-Stage Methanation (TSM). By cooling the lgSNG and passing it through another methanation reactor, a much higher purity of methane can be achieved, making it effectively identical to fossil-sourced natural gas, i.e. [48]. All the parameters for this type of reactor were collected from an industry supplier source as well.

The investment costs associated with the methanation reactor(s) are:

CAPEXmeth=cinvOSM×m˙lgSNG,peakOSM+cinvTSM×m˙CH4,peakOSM

Parametrization of methanation technologies in the MILP model

Technology

Parameter

Symbol

Value

Reference

OSM

CO2 conversion rate

χCO2OSM

94 %

Industry source

H2 conversion rate

χH2OSM

93 %

Industry source

LHV of SNG

LHVlgSNG

44.52 MJ/kg

Calculated from industry source

Specific CAPEX

cinvOSM

820 kEUR/MW

Industry source

TSM

CO2 conversion rate

χCO2TSM

100 %

Industry source

H2 conversion rate

χH2TSM

99 %

Industry source

LHV of SNG

LHVSNG

50.36 MJ/kg

Calculated from industry source

Specific CAPEX

cinvTSM

1,172 kEUR/MW

Industry source

Energy Storage

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 SoCpeak , the investment cost for storages is:

CAPEXstorage=cinvbat×SoCpeakbat+cinvH2 tank×SoCpeakH2 tank

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 cinvbat =100 kEUR/MWh, in line with European Commission expectations of battery packs falling below 100 USD/kWh in 2027 [50]. For the hydrogen tank storage the round-trip efficiency was set at 80% as an intermediate point between the 60~70% range of liquefied storage and the (pressure-dependent) 80~95% range of compressed storage [51], the auxiliary demand equal to the battery, and the investment cost assumed to be a higher cinvH2 tank = 400 kEUR/MWh, with these less optimistic parametrization likely favoring batteries.

Parametrization of energy storage technologies in the MILP model

Technology

Parameter

Symbol

Value

Reference

Battery

SoC loss

λbat

1 %/h

[49]

Roundtrip eff.

ηbat

100 %

Simplification

Specific CAPEX

cinvbat

100 kEUR/MWh

[50]

H2 tank

SoC loss

λH2 tank

1 %/h

[49]

Roundtrip eff.

ηH2 tank

80 %

[51]

Specific CAPEX

cinvH2 tank

400 kEUR/MWh

Assumption

Objective Function

Bringing together eqs. (6), (7), (8) and (9) yields the total annual OPEX:

OPEX=OPEXgas+OPEXpower+OPEXH2+OPEXCO2

while the total CAPEX are obtained from eqs. (13), (17), (20), (24) and (25):

CAPEX=CAPEXren+CAPEXelec+CAPEXCC+CAPEXmeth+CAPEXstorage

bringing the total annualized costs to:

annual cost=OPEX+i1+in1+in1CAPEX

with an interest rate i =3% and amortization period n =20 years yielding a CAPEX annualization factor of 6.7%.

Scenario Definition

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.

Unconstrained grid regime. With a peak capacity of around 7.2 MW/MWq, it represents a scenario where all the investments in grid capacity are realized and the mill can draw, in practice, as much power as it needs from the grid. This boundary condition potentially tilts the optimal solution towards a more interconnected system.

Partially constrained grid regime. With a 0.72 MW/MWq peak capacity, this roughly corresponds to the current situation before decarbonization and electrification measures are undertaken.

Strongly constrained grid regime. With only 0.24 MW/MWq of peak capacity available, which equals around one-third of today’s maximum theoretical capacity, it represents a pessimistic scenario in which grid expansion investments to alleviate structural bottlenecks lag behind industry decarbonization commitments and thus place limits on the mill’s ability to draw electrical power from the grid, potentially pushing the optimal solution towards more electricity-efficient technologies and/or a more self-contained system with expanded onsite renewable generation.

Results and Discussion

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.

Energy Costs and Optimal Design

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.

Unconstrained grid. In the unconstrained grid scenario (top bar of each trio in Figure 2), most of the electricity is supplied directly from the grid, with only a minor contribution from wind power, and lesser amounts from hydro and PV, with no battery storage required. Hydrogen is produced largely onsite with a smaller contribution from grid purchases.

Partially constrained grid. In the partially constrained grid scenario (middle bar), around half of the electricity demand is satisfied by the grid and the other half by renewables (mostly wind) and some battery storage. However, the total electricity consumption is low, because the optimal solution involves purchasing most hydrogen from the grid.

Strongly constrained grid. In the scenario where the grid is most strongly constrained (bottom bar), it can only supply part of the necessary electrical power, so a large capacity of wind power is installed and there is also substantial buying of hydrogen from the grid. Flexibility to offset the intermittency of wind generation is provided through three mechanisms: exporting excess electricity to the grid, charging or discharging the battery storage depending on the supply and demand at each time step, and producing excess hydrogen and exporting it.

CAPEX-OPEX Trade-Offs and Optimal Operation

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.

Electrolysis subsystem. Figure 5 shows the operation of the electrolysis subsystem during one of the 24-hour representative periods (specifically, October 10th). It is apparent from the plots that, regardless of the decarbonization route and grid capacity the three methods of supplying green H2 fulfil specific roles. First, the high-CAPEX, low-OPEX SOEC acts as a base load supply, operating at a near-constant, very high load factor which takes advantage of its high efficiency to offset the larger investment cost. Second, the low-CAPEX, high-OPEX AWE/PEM follows the demand profile more closely whenever abundant wind power is available (typically at night), taking advantage of its more flexible operation possibilities. Finally, the option to directly purchase green hydrogen from the grid is used for peak-of-the-peak supply, thus avoiding extra electrolyser capacity that would otherwise go largely unused; it is also used to fulfil variable loads when the availability of wind power is low, particularly in the partially constrained grid scenario.

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.

CC subsystem. In the carbon capture subsystem, meanwhile, the optimal operation is strongly related with the availability of electricity. In the unconstrained grid scenario, the PSA (using abundant electricity) is used to fulfil the baseline demand, while the AS (using waste heat for which it competes with the SOEC) is activated as a backup when the load of the furnaces is high. In the two constrained grid scenarios, on the other hand, the PSA provides a small, uniform capacity while the AS handles most of the carbon capture demand. In these two latter scenarios, due to the lower capture rate, a larger carbon slip is present which must be offset by a larger purchase of biomethane. These results again highlight the crucial role that the grid conditions have on the optimal decarbonization strategy.

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.

Parameter Variations

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.

Wind CAPEX. The 60 EUR/MWh LCoE assumed for the onshore wind production is realistic but somewhat on the optimistic side, for two reasons. On the one hand, one could expect a "low-hanging fruit" effect by which the best locations for wind parks are built out first, and future expansions will have a lower average load factor and thus higher LCoE. On the other hand, transmission lines required to connect new wind parks directly to the steel mill could have an extra cost. Therefore, the cost variations studied are 1×, 1.25× and 1.5× the baseline value of cinvwind.

Hydrogen purchase price. The 150 EUR/MWh assumed as baseline, corresponding to about 4.5 EUR/kg, is within the lower range of current prices for renewable H2 currently, but it is expected that prices will go down in the future thanks to economies of scale and government subsidies (the "hydrogen economy"). Thus, the studied variations are 1×, 0.8× and 0.6× the baseline cost cbuyH2.

Electrolysers CAPEX. The investment costs for electrolysis systems have uncertainties in both directions. On the one hand, it is likely that technological innovations and economies of scale will bring down the production cost of electrolysis cells. However, it is also possible that ancillary infrastructure and adaption of gas supply systems in the furnaces (for the H2) will induce additional investment costs. To account for all the possibilities, the studied variations are 0.5×, 0.75×, 1×, 1.5× and 2× the baseline cinvAWE/PEM and cinvSOEC.

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.

Implications for Business

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.

Conclusion

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.

Future Work

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.

Acknowledgment
Acknowledgement

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”.

REFERENCES
  1. SteelFacts
  2. World Steel in Figures 2025
  3. Tracking industrial energy efficiency and CO2 emissions: In support of the G8 plan of action
  4. Leicher J., Giese A., Wieland C., 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
  5. Pardo N., Moya J., Vatopoulos K. Prospective scenarios on energy efficiency and CO2 emissions in the EU iron & steel industry: Re-edition
  6. Schmitz N., Sankowski L., Kaiser F., Schwotzer C., Echterhof T., Pfeifer H., 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
  7. Häggström G., Electrification of reheating furnaces: state of the art and future research needs
  8. Della Rocca A., Astesiano D., Malfa E., Rolling mill decarbonization: Tenova SmartBurners with 100% hydrogen, Matériaux & Techniques, Vol. 109 (3-4), , https://doi.org/10.1051/mattech/2022012
  9. Ameli H., Strbac G., Pudjianto D., Ameli M., 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
  10. Airaksinen S., Haapakangas J., Gyakwaa F., Heikkinen E., Fabritius T., 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
  11. Emel’yanov A., Lube I., Kuznetsov V., Levchenko D., Scale Formation Minimization on Internal Surface of Seamless Hot-Rolled Pipes, Steel Transl., Vol. 50 (5), :340-346, https://doi.org/10.3103/S0967091220050046
  12. Haapakangas J., Airaksinen S., Heikkinen E., Fabritius T., 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
  13. Haapakangas J., Riikonen S., Airaksinen S., Heikkinen E., Fabritius T., Oxide Scale Formation on Low-Carbon Steels in Future Reheating Conditions, Metals, Vol. 14 (2), , https://doi.org/10.3390/met14020189
  14. Şahin A., Işikgül A., Işikgül O., Imer S., Yaşacan D., Investigation of Decarburization of Spring Steel During Hot Rolling
  15. Yamamoto S., Sayama S., Kunitomi S., Ogasawara K., Baba N., 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
  16. Giglio E., Lanzini A., Santarelli M., Leone P., 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
  17. Zabik G., Birkelbach F., Hofmann R., 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
  18. Miserocchi L., Franco A., 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
  19. Guo Z., Zhou S., Modeling and Multi-Stage Planning of Cement-IIES Considering Carbon-Green Certificate Trading, Processes, Vol. 11 (4), , https://doi.org/10.3390/pr11041219
  20. Tiggeloven J., Faaij A., Kramer G., Gazzani M., 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
  21. Löffler S., , Transformation of a Steel Rolling Plant: Identification of Decarbonisation Pathways using MILP
  22. , APG’s Future Grid Puts Supply-Secure Energy Transition on the Fast Track, Austrian Power Grid
  23. , Energy prices, taxes, Statistik Austria
  24. Wlcek B. Marktbericht 2023 (accessed )
  25. Kutz C. Potenzialstudie Grüne Gase: Analyse und Bewertung der Potenziale Grüner Gase in der Innovationsregion Mitteldeutschland (accessed )
  26. , Real-Time Electricity Tracker, International Energy Agency
  27. Kathan J. Importmöglichkeiten für erneuerbaren Wasserstoff (accessed )
  28. , , Ministerialentwurf Erneuerbares-Gas-Gesetz
  29. , Emissionspreisindex, Austrian Energy Agency
  30. Altaghlibi M., Gentile G., Scenarios shaping EU ETS prices
  31. Jaksch-Fliegenschnee M., Wonisch P., Krenn A. IEA Wind TCP Austria 2022
  32. Kost C., Müller P., Sepúlveda Schweiger J., Fluri V., Thomsen J. Study: Levelized Cost of Electricity- Renewable Energy Technologies
  33. Fechner H. National Survey Report of PV Power Applications in Austria 2023 (accessed )
  34. Philipps S., Warmuth W. Photovoltaics Report (accessed )
  35. Sebbahi S., Nabil N., Alaoui-Belghiti A., Laasri S., Rachidi S., Hajjaji A., 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
  36. Holst M., Aschbrenner S., Smolinka T., Voglstätter C., Grimm G. Cost Forecast for Low-Temperature Electrolysis - Technology Driven Bottom-Up Prognosis for PEM and Alkaline Water Electrolysis Systems
  37. Electrolyzers for Hydrogen Production: Technical and Economic Characteristics Executive Summary
  38. van 't Noordende H., van Berkel F., Stodolny M. Next Level Solid Oxide Electrolysis: Upscaling potential and techno-economical evaluation for 3 industrial use cases
  39. , Silyzer 300: The next paradigm of PEM electrolysis, Siemens Energy
  40. , Sunfire-HyLink SOEC, Sunfire SE
  41. Bukar A., Asif M., 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
  42. Wilk V., Leibetseder D., Zauner C., Rath A., Schwaiger M., Improving Energy Efficiency of Carbon Capture Processes with Heat Pumps
  43. Riboldi L., Bolland O., 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
  44. Susarla N., 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
  45. Subraveti S., Roussanaly S., Anantharaman R., Riboldi L., Rajendran A., 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
  46. Husebye J., Brunsvold A., Roussanaly S., Zhang X., 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
  47. Rönsch S., Review on methanation – From fundamentals to current projects, Fuel, Vol. 166 , :276-296, https://doi.org/10.1016/j.fuel.2015.10.111
  48. Tommasi M., Degerli S., Ramis G., Rossetti I., 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
  49. , 5 MWh Battery Energy Storage System for North America, Chint Power Systems
  50. Bielewski M. Battery technology in the European Union: 2024 Status report on technology development, trends, value chains and markets
  51. Gardiner M. DOE Hydrogen and Fuel Cells Program Record 9013: Energy requirements for hydrogen gas compression and liquefaction as related to vehicle storage needs

DBG