1. Introduction
To meet the CO
2 reduction targets of the Paris Agreement on climate change, countries around the world are developing strategies to produce hydrogen and apply it in the industrial, mobility, and energy sectors [
1,
2]. Therefore, the USA is aiming to halve its green house as emissions by 2030 compared to 2005 levels and achieve complete climate neutrality by 2050 [
3]. Furthermore, all countries within the EU are legally obligated to reduce their carbon emissions to 55% of 1990 levels by 2030 and to achieve climate neutrality by 2050 [
4]. In 2021, the German government set even more ambitious targets to reduce greenhouse gas emissions by 65% by 2030 and 88% by 2040 compared to 1990 levels and to achieve climate neutrality by 2045 [
5]. Increasing energy efficiency is seen as a fundamental pillar of the energy transition along with the expansion of renewable energies, as it requires less energy infrastructure to operate production processes [
6] (p. 7) [
7].
With an annual energy consumption of 4.2 TWh and estimated total emissions of 1.8 million tons of CO
2-equivalent, non-ferrous metal foundries (abbreviated as NF-metal foundries) are among the particularly energy- and emission-intensive industries, with the larger share of 2.4 TWh being provided by fossil fuels. Additionally, 50% of the total energy consumption is attributed to the melting process [
8] (p. 1) [
9]. Green hydrogen enables the partial or complete reduction of greenhouse gas emissions through the blending or complete substitution of fossil fuels.
However, integrating hydrogen into energy systems presents several challenges. In Ref. [
10], the authors evaluated over 59 field projects and 33 stakeholder dialogues, identifying five barrier domains (technical, infrastructural, socioeconomic, environmental, and institutional) which need to be overcome for the Netherlands. Focusing on the technical domain, challenges in efficiency and durability on the component level are mentioned. In the case of fuel-fired thermal processing plants, this may require adjustments on the burner unit or combustion chamber, as hydrogen combustion is associated with higher fuel volume flows at identical firing rates, higher reactivity, and higher adiabatic flame temperatures [
11]. Therefore, challenges such as significantly increased NOx emissions and material resistance at high temperatures may occur [
12] [
13] (pp. 8–28) [
14].
Crucible furnaces are mainly used in NF-metal foundries and currently operate electrically or with fossil fuels. Switching to hydrogen thus presents significant potential for reducing CO
2 emissions in NF-metal foundries [
15]. Computational fluid dynamics (CFD) simulations can help to identify energy-saving potentials and to make design decisions, such as adapting the geometry to the changed requirements of a hydrogen flame [
16] (pp. 8–10) [
17] (p. 117) [
18] (p. 21). In Ref. [
19], a literature review on CFD simulations of fuel-fired crucible furnaces was conducted, resulting in the identification and assessment of models for fossil-fired crucible furnaces with varying modelling quality.
A comparison of key performance indicators by changing the fuel composition to hydrogen, however, has not yet been carried out, but due to substantially different combustion characteristics, significant changes on these can be expected [
13]. High energy utilization rates and even heat distributions on the crucible surface are desired to achieve uniform melt bath temperatures and minimal thermal stresses within the materials [
20]. Therefore, in this paper, two CFD models—one for a hydrogen and a second for a methane flame in a crucible furnace designed for fossil fuels—are set up, and key performance indicators at identical firing powers are compared.
To do so, at first, typical crucible furnace operating conditions are described, from which the case study is derived. Then, based on literature research for jet and swirl burners for methane and hydrogen flames, applicable numerical models are selected, and the simulation setup is placed in the context of real operating modes. Subsequently, temperature and radiation fields, radiation parameters, and heat transfer characteristics are evaluated to support hydrogen-related design decisions within the crucible furnace and burner unit.
2. Furnace Case Study and Numerical Modelling
In this section, combustion models representing methane and hydrogen flames in jet and swirl burners are researched and evaluated. Then, the case study is described, and the numerical simulation model of the crucible furnace, considering its operational characteristics, is described.
2.1. State of the Art in Modelling Natural Gas, Methane, Hydrogen, and Hydrogen-Blended Flames in Jet and Swirl Burners
A literature review was conducted on the numerical simulation of jet and swirl burners to select numerical models, capable of describing methane and hydrogen flames for the case study. Numerous publications successfully use numerical methods to describe temperature, flow, and species fields for a wide range on burners for natural gas (NG), methane, hydrogen, and hydrogen-blended flames. The Reynolds-averaged Navier–Stokes equations (RANS equations) are commonly employed using turbulence and combustion models for solving these simulations.
In Refs. [
21,
22], the authors simulate jet burners with varying hydrogen–natural gas compositions using the RANS approach and validate the results against experimental data.
In Ref. [
21], combustion is modelled using a combined finite-rate and eddy break-up model, the second one uses a flamelet approach with assumed probability density functions. The first publication employs the standard k-ε, the renormalization group k-ε (RNG), and the Reynolds stress model, finding that all turbulence models have similar axial profiles and are therefore suitable for describing the flame. Radiation is modeled using the P1 radiation model. Validation was performed on measurement data evaluating the maximum temperature of the hydrogen flame, exhaust gas temperatures for H
2 and CH
4 as well as the exhaust gas composition showing overall adequate results except for the CH
4 exhaust gas temperature which was predicted to low. The authors attribute the deviation to a slightly different location of the measurement point compared to the simulation and model inaccuracies.
The latter publication applies the standard k-ε-, the Reynolds stress, as well as a modified k-ε turbulence model. Validation was performed for the CH4 flame with three radial temperature profiles at different combustion chamber depths. The results show that all models describe the temperature profiles as qualitatively adequate, with the modified k-ε turbulence model providing the most accurate results. If and which radiation model was used is not described.
In Refs. [
23,
24,
25,
26], the combustion of variable hydrogen fractions in different swirl burners is simulated, validated, and evaluated. To describe the turbulence, the first publication applies the Realizable k-ε and the RNG models, yielding better results for the former. Furthermore, the authors state that in principle, both models are applicable to the slow-swirling flow in their case. The discrete ordinate method has been used for the radiation modelling and the eddy dissipation combustion model, with adjusted reaction mechanisms for the different fuel compositions used. The simulation data were compared with temperature and exhaust gas composition measurements at the outlet for the natural gas case yielding a temperature deviation of about 14%.
In Ref. [
24] the authors use the coupled finite-rate and eddy break-up model to describe the combustion, the standard k-ε as well as the RNG model for the turbulence, and the discrete ordinate method for radiation. In their validation, they compare the temperature fields of different hydrogen contents in a cross section, with measurement data showing that both models can describe the flame structure with more accurate results for the RNG model.
To describe the turbulence, the authors in Ref. [
25] use the RNG model, the combustion is modeled with a finite rate approach, and the radiation is represented with the P1 model. They successfully validate their results with temperature measurements at the flame center.
In Ref. [
26], a 100% natural gas flame and a mixture containing 30% hydrogen are simulated, applying different flamelet approaches as well as the realizable k-ε and the Reynolds stress turbulence models. If and which radiation model was used is not described. The simulation results were validated against flow patterns, radial temperature profiles, and exhaust gas concentrations, yielding the most accurate results with the partially premixed combustion and the Reynolds stress turbulence model.
In
Table 1, the numerical methods used for describing the flames and the maximum temperatures are summarized. In summary, the temperature and flow fields of the mentioned fuel mixtures in jet and swirl burners can be simulated using the RANS equations in combination with the coupled finite-rate eddy dissipation or flamelet approach as well as the P1 or discrete ordinates radiation model. For swirl burners, models such as the Realizable k-ε or Reynolds stress model, which deviate from the standard k-ε Model, provide more accurate results in swirling flows. In all publications, an increase in hydrogen fraction is observed to result in higher combustion chamber temperatures and NOx emissions, while CO
2 emissions decrease.
2.2. Case Study of the Fuel Fired Crucible Furnace
Figure 1 illustrates the internal contours of the considered crucible furnace in the front and side views, including inlets, outlets, burner positioning, and used burner heads to generate the computational domain for both simulations. The geometry was obtained by simplifying the CAD of a conventional fossil fired crucible. Additionally, the inlet and outlet surfaces on which the boundary conditions were applied are marked. As can be seen in both cases, the same inlet areas are used for air and gas. On the remaining surface, wall boundaries were applied. Such a furnace can be used for melting or holding metals warm and is operated in batch mode for aluminum. During the melting process, the furnace melts the metal at high power and overheats it to 680 °C before keeping it warm at a lower power or removing it for further processing by tilting or ladling. Due to the many time-varying parameters, the furnace operation is considered transient. The burner unit can regulate the target temperature through pulsation or power modulation [
15] [
27] (p. 321) [
28]. As shown in
Figure 1, the fuel gas flows separately from the air through the burner tube and is distributed into the combustion chamber through the nozzle head. The air is injected around the burner head and through corresponding holes. Therefore, fuel and oxidizer are mixed near the nozzle head and downstream until the ignition limit is reached, triggering the combustion reaction and flame formation. To avoid excessive temperatures in the vicinity of the nozzle head, the air holes for mixing were removed in the hydrogen case.
To compare both flames, identical steady-state operating points with the same heating power were selected at a melt bath temperature of 680 °C. This approach neglects time-dependent factors such as thermal masses in the furnace, including walls, burners, and other components. Thus, the simulation represents approximately a snapshot at the end of the melting process when the temperature in the melt bath reaches 680 °C. At this point, the burner operates at its maximum heating power of 150 kW at λ = 1.1 using either methane or hydrogen before transitioning to the holding operation mode or removing the molten material.
2.3. Numerical Setup
Two CFD models were developed and simulated to compare both combustion cases; the corresponding computational domains are shown in
Figure 1. The RANS equations were discretized using the finite volume method and solved using the commercial flow solver ANSYS Fluent
®. As shown in
Section 2.2, this simulation methodology is widely used in combustion simulations, including for melting furnaces [
21,
23,
26,
29,
30,
31].
It is necessary to consider not only the turbulent flow field but also the chemical processes to predict turbulent non-premixed flames formed by the mixing of a fuel and air stream through flow and turbulence. The large number of involved species and interactions between turbulent fluctuations in the flow field and chemical kinetics significantly increase the computational effort, which is why generally simplified models with limited validity are used for calculations [
32] (pp. 343–416) [
33] (p. 38). As described in
Section 2.2, the flamelet model, among others, has been successfully used in the literature to describe hydrogen and methane flames. This approach assumes that the turbulent non-premixed flame can be represented as an ensemble of local laminar flames called flamelets. These, in turn, can be efficiently calculated from the conservation equations of a one-dimensional counterflow flame in a pre-processing step and tabulated, taking extensive reaction mechanisms into account. The chemistry–turbulence interaction in the reaction rate is modelled with assumed probability density functions. Due to its computational efficiency and the ability to consider chemistry–turbulence interaction via reaction mechanisms, this model was selected [
22,
26,
34]. In ANSYS Fluent, the non-premixed combustion model with a diabatic steady-state diffusion flame was accordingly chosen. For the methane flame, the internal GRI 3.0 mechanism was used, and for the hydrogen flame, the hydrogen reaction mechanism was used.
To close the Reynolds stress tensor, the standard k-ε turbulence model was chosen, and the enhanced wall treatment was selected to describe the flow near walls. This is a two-layer model that maps the boundary layer using empirical wall functions depending on the dimensionless wall distance y+ or resolves it explicitly. This allows the (turbulent) viscosity, the (turbulent) dissipation rate ε, and the velocity gradients to be explicitly calculated near the wall, which, for example, allows for a more accurate prediction of heat transfer [
35] (pp. 369–389).
Radiation is modelled using the discrete ordinates model, which transforms the radiative transfer equation into a transport equation for radiation intensity and solves it for a finite number of radiation directions. Additionally, the combustion products H
2O and CO
2 participate in the radiation exchange through absorption and emission and must be considered in the calculations. For this purpose, the weighted sum of gray gases model (Abr. WSGGM) is used, which has established itself in combustion simulations as a compromise between computational time and accuracy [
32] (p. 442).
Spatially, the gradients were discretized using the Green–Gauss cell-based method, pressure with the PRESTO! scheme, and all other variables with the second-order upwind scheme. The pressure-based solver was used for solving, and the SIMPLE scheme was employed for the pressure-velocity coupling. The convergence criterion was set to 10
−3 for the mass, momentum, turbulence, and mixture fraction equations and to 10
−6 for the energy and radiation equations [
36] (pp. 952+954).
Below the equations for species, momentum, energy and radiation are depicted, which are solved with above described modelling:
where
is the density,
the species fraction,
the velocity vector,
the diffusion flux,
the net chemical rate of production of species
i,
p the pressure,
the stress tensor,
the gravitation,
external body forces,
e the internal energy,
h the enthalpy
the effective conductivity,
T the temperature, and
the volumetric source term. Furthermore,
represents radiation intensity,
position,
direction,
scattering direction vector, s as path length,
a as absorption,
as scattering coefficient,
n as refractive index,
as Stefan–Boltzmann constant,
as phase function, and
as solid angle [
36] (pp. 3+160+174+224).
2.4. Boundary Conditions
To solve the posed boundary value problem, all boundary conditions must be defined. Therefore, the exterior walls and the crucible inner walls in contact with the melt were chosen. The explicitly modelled domain describes the inner volume of the furnace, while the external walls or the wall thickness of the crucible up to the melt are taken into account through internal shell elements of the program. Material properties and wall thicknesses can be assigned to these, with which the steady-state heat flows through the walls are described. The corresponding boundary conditions, wall thicknesses, and materials are given in
Table 2 and
Table 3, which describe the 150 kW heating cases at 680 °C melt temperature. In addition, the melt bath is not modelled explicitly but is represented by a constant inner wall temperature. This neglects temperature fluctuations and free convection effects in the melt, which is assumed to have a minor impact compared to the comparison of hydrogen and methane flames.
2.5. Mesh Independence Study
Applying a mesh independence study for the hydrogen case with three uniformly refined meshes, the dependency of the solution to the mesh was evaluated. Calculations were performed on the high-performance computer system of Baden-Württemberg’s Universities and Universities of Applied Sciences bwUniCluster 2.0. In
Figure 2, the volume and surface mesh of the outlet plane of the second refinement level are depicted, and they consist of about 2.6 × 10
7 nodes and 7.5 × 10
6 cells, respectively. In addition, on the lower right side, the contours of the burner with the air and gas inlets are shown. For uniformly coarser meshes below 1*10
7, the convergence limits could not be reached. The mesh resolution is considered sufficiently accurate if no change can be seen in the solution upon further refinement. The total heat transfer into the crucible and the average temperature of the outlet plane were chosen as control quantities and are depicted in
Figure 3. As can be observed, only minor changes for both maximum quantities—1% and below to the finest mesh—can be stated, and therefore, all meshes identify as sufficiently accurate to describe the control quantities. As a compromise between solution accuracy and stability, as well as computational resources, the second mesh is used in the subsequent investigation.
3. Evaluation
In the following, the combustion characteristics of the methane and hydrogen flame with previously described numerical setup, and operating conditions are compared by evaluating temperature and radiation fields as well as radiation parameters. Furthermore, the results are plausibilized with the literature in
Section 2.1. Then, the energy utilization rates and heat distributions are analyzed, identifying as one of the key performance indicators for the furnace operation. Therefore, the total convective and radiative heat fluxes, as well as the heat distribution into the crucible, were compared and related to the changed combustion characteristics. Finally, the results are discussed in the context of furnace design and burner adjustments.
3.1. Characterization and Comparison of the Temperature Fields the Methane and Hydrogen Flames
Figure 4 depicts the temperature fields of the inlet and outlet cross sections for the hydrogen and methane flames as well as the volume-averaged temperatures. Additionally, regions of the crucible furnace are divided into the areas inlet, outlet, front side, and back side. To evaluate both simulations, in
Table 4, the temperatures were initially classified as either medium or high and then compared qualitatively.
Overall, a similar combustion characteristic of both flames can be observed. The reactants begin to mix in the inlet area, forming a reaction zone with strong temperature gradients that continues around the back side of the crucible. Further downstream on the outlet side, the highest temperatures and homogenization of the temperature profile across the cross-section can be observed. At this point, a small part of the flow directly exits the furnace space, leading to very high temperatures in the outlet area and indicating possible approaches for optimizing the furnace geometry. The larger flow portion continues around the crucible, where further homogenization and a slight decrease in temperatures can be seen. Subsequently, the flow leaves the combustion chamber via the outlet.
Comparing both flames, similar temperature fields with slightly higher temperatures and broader reaction zone of the hydrogen flame can be observed in the inlet area and on the back side. Further downstream around the crucible on the outlet side, the temperature rises by about 200 K above that of the methane flame due to the higher flame temperatures of hydrogen. Overall, an increase in the average combustion chamber temperature of 61 K from 1434 K in the case of methane to 1495 K for the hydrogen flame can be observed, which is largely attributable to the outlet and front side downstream.
To verify the simulation setup, the maximum flame temperatures and the combustion characteristics are compared with the literature in
Table 1. As a result, similar maximum temperatures for the hydrogen and methane flames were observed, ranging between the lowest and highest maximum flame temperatures. In addition, the maximum temperature increases by about 258 K when switching to hydrogen, which is also within the range attested in the literature.
Comparing the flame characteristic by switching to hydrogen, all authors claim an increase in temperature due to the higher adiabatic flame temperature, which is also observed in
Figure 4. Ref. [
23] describes a longer, broader flame in more detail as well as increased radial temperatures with increasing hydrogen content, which we also observe (
Figure 4).
As a result, comparing the maximum temperatures and the changed combustion characteristics, both flames presented are within a plausible range. Deviations within that range seem plausible and may occur due to different numerical setups, boundary conditions, burner or combustion chamber geometry.
3.2. Characterization and Comparison of the Radiation Characteristics the Methane and Hydrogen Flame
Looking at the net radiation power of the gases in
Figure 5, patterns similar to the temperatures can be seen. In areas of low temperatures, low or negative net radiation powers can be observed, indicating a decrease in radiation intensity, whereas in areas of high temperatures, high net radiation powers of the gases can be observed, which is reflecting the strong temperature dependence of radiative heat exchange.
Figure 5 shows the corresponding absolute absorption, emission, and net radiation powers across the temperature quantiles. For both flames, up to about the 50% temperature quantile, the absorption rates exceed the emission rates, indicating a decrease in radiation intensity. Consequently, in these cells, net radiation heat is absorbed by the gas. Beyond that, the emission rates dominate due to the significantly stronger temperature than absorption dependency, resulting in a sharply increased radiation power with rising temperature. Thus, the regions of the inlet area and the back side can be associated with comparatively low net heat radiation, with the highest found in the outlet and front side area.
Comparing the emission and absorption behavior of both flames in
Figure 6 (top), higher emission and absorption rates can be observed for the methane flame at the respective temperature quantiles. With increasing temperature quantiles, both ratios decrease continuously, and at the 85% quantile, the emission rates of the hydrogen flame exceed those of the methane flame. Overall, up to the 0.8% temperature quantile, the different radiation characteristics do not lead to any change in the net radiation behavior. From there, the hydrogen flame tends to emit more net power at each quantile.
According to the radiative transfer equation and by neglecting scattering effects, the change in radiation intensity depends on temperature as well as emission and absorption coefficients. Therefore, temperature quantiles and the associated volume are shown in
Figure 6 (bottom). As seen, the respective temperature quantiles of both flames occupy the same volumes, with an over-proportional rise in temperatures for the hydrogen flame from approximately the 60% temperature quantile, indicating smaller emission and absorption coefficients for the hydrogen case at the same temperature. Due to the similar temperature fields of both flames, the temperature ranges occupy approximately the same regions. Between the 0.6–0.8 quantiles, higher temperatures are needed for the same net emission, which then turns into higher net emission rates for further rising temperatures. Therefore, higher net radiation powers occur mainly in areas of high temperatures above ~1750 K in the hydrogen case. The difference in net radiation powers of both flames correlates with the temperatures, suggesting a linear relationship. Consequently, the increased net radiation can be almost exclusively related to higher temperatures than the methane flame.
3.3. Characterization and Comparison of the Crucible Heat Input of the Methane and Hydrogen Flames
Apart from a high total heat transfer rate into the crucible, which describes the energy utilization at the simulated operating point, an even heat distribution on the crucible surface is desired since a uniform melt bath temperature and minimal thermal stresses are promoted. Therefore, the following evaluation of energy utilization is based on the transferred absolute radiation and convective heat flux (
Figure 7, bottom7). For the even heat distribution, the spatial distribution of the heat flux density, its statistical moments, and the associated histogram are analyzed (
Figure 7, top).
At the bottom of
Figure 7, an overall increase of the total heat transfer by 19.0% for the hydrogen flame can be observed, indicating a significantly enhanced energy utilization of the fuel, which is almost exclusively due to an increased radiant heat share. The convective portion remains nearly constant (hydrogen: 4.5 kW, methane: 4.6 kW) and plays a minor role in heat transfer for both flames (hydrogen: 12.1%, methane: 14.8%).
The heat transfer on the crucible surface is similar for both flames (cf.
Figure 8), with low heat transfer rates on the inlet and back side and higher rates on the outlet side, where the hydrogen flame reaches higher transfer rates across the entire crucible surface (cf.
Figure 8, top and bottom). On the outlet side, the transfer rates can exceed those of the front side, leading to high standard deviations of 28 (hydrogen) and 22 kW/m
2 (methane), observed in the distribution function. Compared to the increase in the absolute amount of heat transferred by 19.0%, the standard deviation in the hydrogen case rises disproportionately by 27.3%. This describes a stronger scattering of the heat flux on the crucible surface around its mean value and thus a more uneven heat flux distribution, mainly due to the very high maximum transfer rates on the outlet side. Therefore, in the hydrogen case, more uneven melt bath temperatures and higher thermal stresses can be expected.
In summary, a significantly increased fuel utilization (19.0%) and thus higher energy efficiency can be observed during hydrogen operation, which could almost exclusively be attributed to an increased radiant heat share. The higher heat transfer rates could clearly be assigned to the outlet side, resulting in a significantly worse heat flux distribution. Furthermore, the net radiation power of the gases in both simulations is identified as the main influencing factor on the behavior of heat transfer to the crucible. In turn, this strongly depends on the temperature as well as the absorption and emission coefficients, which differ significantly for both flames (radiation H
2O and CO
2). In
Section 3 (cf.
Figure 5 and
Figure 6), the increased net radiation of the hydrogen case was related to the higher temperatures compared to the methane flame, which is also identified as the main cause of the increased energy utilization rate and poorer heat distribution.
With higher and more uneven heat transfer rates into the crucible, thermal stresses in the material increase, and more uneven melt bath temperatures are to be expected. Furthermore, higher temperatures can lead to an exponential increase in NOx emissions, which may exceed legal limits. For robust and compliant operation, measures such as increased air numbers or exhaust gas recirculation may be necessary to reduce combustion temperatures or further adjustments. How energy efficiency behaves with the implementation of measures aimed at reducing combustion temperatures remains to be evaluated. The results presented here suggest that radiation heat transfer rate into the crucible and thus energy efficiency should decrease with reduced temperatures.
4. Outlook
In this publication, a hydrogen and a methane flame with identical firing power were simulated in a crucible furnace at a steady-state operating point at 680 °C, describing the end of a melting run. Switching the fuel from methane to hydrogen significantly changes the temperature and flow fields as well as the heat transfer behavior into the crucible. Increased furnace chamber temperatures for the hydrogen case were determined, which could predominantly be attributed to the front and outlet sides. Furthermore, maximum temperatures and flame characteristics could be successfully corroborated with results from the literature.
An increased fuel utilization rate could be determined, which led to a 19.0% higher absolute heat flux into the crucible, which could allow for shorter heating times or a reduction in power. The higher heat flux could predominantly be attributed to the radiative heat part and the over-proportional increase in temperatures compared to methane starting at about 1550 K (cf.
Figure 6 and
Figure 7). Simultaneously, a disproportional worsening of the spatially distributed heat transfer rate was determined (increase of the standard deviation by 27.3%).
Furthermore, the evaluation results indicate that measures to reduce combustion temperatures—such as exhaust gas recirculation, leaner operation to reduce NOx emissions, or compliance with temperature limits of materials—should result in a decreasing fuel utilization rate and more evenly distributed heat flux into the crucible. Additionally, constructive adjustments to the combustion chamber geometry, such as repositioning the outlet or extending the burner introduction, could increase the residence time of the hot exhaust gases in the furnace and thereby increase fuel utilization.
To further examine the conclusions of this publication, future research could focus on a more extensive parameter study, systematically investigating the aforementioned influencing parameters (e.g., firing rate, air number, burner and furnace geometry) and comparing them to the original fossil-fuel-fired configuration. Additionally, the model could be developed to be fully transient to be able to describe the full transient heating up case, aiming to quantify the fuel utilization rate over the whole heating process.