The operation of large industrial furnaces will continue to rely on hydrocarbon fuels in the near foreseeable future. Mathematical modeling and numerical simulation is expected to deliver key insights to implement measures to further reduce pollutant emissions. These measures include the design optimization of the burners, the dilution of oxidizer with exhaust gasses, and the mixing of natural gas with hydrogen. In this paper, we target the numerical simulation of non-premixed turbulent combustion of natural gas in a single heating section of a ring pit anode baking furnace. In previous work, we performed combustion simulations using a commercial flow simulator combined with an open-source package for the three-dimensional mesh generation. This motivates switching to a fully open-source software stack. In this paper, we develop a Reynolds-Averaged Navier-Stokes model for the turbulent flow combined with an infinitely fast mixed-is-burnt model for the non-premixed combustion and a participating media model for the radiative heat transfer in OpenFoam. The heat transfer to the refractory brick lining is taken into account by a conjugate heat transfer model. Numerical simulations provide valuable insight into the heat release and chemical species distribution in the staged combustion process using two burners. Results show that at the operating conditions implemented, higher peak temperatures are formed at the burner closest to the air inlet. This results in a larger thermal nitric-oxide concentration. The inclusion of the heat absorption in the refractory bricks results in a more uniform temperature on the symmetry plane at the center of the section. The peak in thermal nitric-oxides is reduced by a factor of four compared to the model with adiabatic walls

This work studies how non-premixed turbulent combustion in a rotary kiln depends on the
geometry of the secondary air inlet channel. We target a kiln in which temperatures can reach values above 1800 degrees Kelvin. Monitoring and possible mitigation of the thermal nitric-oxide (NOx) formation is of utmost importance. The performed reactive flow simulations result in detailed maps of the spatial distribution of the flow, thermodynamics and chemical conditions of the kiln. These maps provide valuable information to the operator of the kiln. The simulations show the difference between the existing and the newly proposed geometry of the secondary air inlet. In the existing configuration, the secondary air inlet is rectangular and located above the base of the burner pipe. The secondary air flows into the furnace from the top of the flame. The heat release by combustion is unevenly distributed throughout the flame. In the new geometry, the secondary air inlet is an annular ring placed around the burner pipe. The secondary air flows circumferentially around the burner pipe. The new secondary air inlet geometry is shown to result in a more homogeneous spatial distribution of the heat release throughout the flame. The peak temperatures of the flame and the production of thermal NOx are significantly reduced. Further research is required to resolve limitations of various choices in our modeling approach.