The reduction of emissions from large industrial furnaces critically relies on insights gained from numerical models of turbulent non-premixed combustion. In the article Mitigating Thermal NOx by Changing the Secondary Air Injection Channel: A Case Study in the Cement Industry, the authors present the use of the open-source OpenFoam software environment for the modeling of the combustion of Dutch natural gas in a cement kiln operated by our industrial partner. In this paper, various model enhancements are discussed. The steady-state Reynolds-Averaged Navier-Stokes formulation is replaced by an unsteady variant to capture the time variation of the averaged quantities. The infinitely fast eddy-dissipation combustion model is exchanged with the eddy-dissipation concept for combustion to account for the finite-rate chemistry of the combustion reactions. The injection of the gaseous fuel through the nozzles occurs at such a high velocity that a comprehensive flow formulation is required. Unlike in Mitigating Thermal NOx by Changing the Secondary Air Injection Channel: A Case Study in the Cement Industry, wave transmissive boundary conditions are imposed to avoid spurious reflections from the outlet patch. These model enhancements result in stable convergence of the time-stepping iteration. This in turn increases the resolution of the flow, combustion, and radiative heat transfer in the kiln. This resolution allows for a more accurate assessment of the thermal NO-formation in the kiln. Results of a test case of academic interest are presented. In this test case, the combustion air is injected at a low-mass flow rate. Numerical results show that the flow in the vicinity of the hot end of the kiln is unsteady. A vortex intermittently transports a fraction of methane into the air stream and a spurious reaction front is formed. This front causes a transient peak in the top wall temperature. The simulated combustion process is fuel-rich. All the oxygen is depleted after traveling a few diameters into the kiln. The thermal nitric oxide is formed near the burner and diluted before reaching the outlet. At the outlet, the simulated thermal NO concentration is equal to 1 ppm. The model is shown to be sufficiently mature to capture a more realistic mass inflow rate in the next stage of the work.@en

No alternatives are currently available to operate industrial furnaces, except for hydrocarbon fuels. Plant managers, therefore, face at least two challenges. First, environmental legislation demands emission reduction. Second, changes in the origin of the fuel might cause unforeseen changes in the heat release. This paper develops the hypothesis for the detailed control of the combustion process using computational fluid dynamic models. A full-scale mock-up of a rotary cement kiln is selected as a case study. The kiln is fired by the non-premixed combustion of Dutch natural gas. The gas is injected at Mach (Formula presented.) via a multi-nozzle burner located at the outlet of an axially mounted fuel pipe. The preheated combustion air is fed in (co-flow) through a rectangular inlet situated above the attachment of the fuel pipe. The multi-jet nozzle burner enhances the entrainment of the air in the fuel jet. A diffusion flame is formed by thin reaction zones where the fuel and oxidizer meet. The heat formed is transported through the freeboard, mainly via radiation in a participating medium. This turbulent combustion process is modeled using unsteady Favre-averaged compressible Navier–Stokes equations. The standard k- (Formula presented.) equations and standard wall functions close the turbulent flow description. The eddy dissipation concept model is used to describe the combustion process. Here, only the presence of methane in the composition of the fuel is accounted for. Furthermore, the single-step reaction mechanism is chosen. The heat released radiates throughout the freeboard space. This process is described using a P1-radiation model with a constant thermal absorption coefficient. The flow, combustion, and radiative heat transfer are solved numerically using the OpenFoam simulation software. The equations for flow, combustion, and radiant heat transfer are discretized on a mesh locally refined near the burner outlet and solved numerically using the OpenFoam simulation software. The main results are as follows. The meticulously crafted mesh combined with the outlet condition that avoids pressure reflections cause the solver to converge in a stable manner. Predictions for velocity, pressure, temperature, and species distribution are now closer to manufacturing conditions. Computed temperate and species values are key to deducing the flame length and shape. The radiative heat flux to the wall peaks at the tip of the flame. This should allow us to measure the flame length indirectly from exterior wall temperature values. The amount of thermal nitric oxide formed in the flame is quantified. The main implication of this study is that the numerical model developed in this paper reveals valuable information on the combustion process in the kiln that otherwise would not be available. This information can be used to increase fuel efficiency, reduce spurious peak temperatures, and reduce pollutant emissions. The impact of the unsteady nature of the flow on the chemical species concentration and temperature distribution is illustrated in an accompanying video.

@en