Evaluation AC2-08
Premixed Methane-Air Swirl Burner (TECFLAM)
Application Challenge 2-08 © copyright ERCOFTAC 2011
In order to assess the simulation a comparison of the time averaged velocity, temperature and species distribution with experimental data is presented in this section. A more extensive evaluation can be found in the cited references. Kuenne et al.[1] .
Comparison of CFD Results with Experimental Data
RANS Simulations of the Isothemal and Reacting Flow
In Fig. 5.1 results of the isothermal flow field obtained with FASTEST and the k-epsilon model are shown in comparison with experimental data. Typically for this type of swirl flow, the velocity field starts to expand right after the nozzle exit, which produces a positive pressure gradient in the axial and radial direction, high enough to form a central recirculation zone as intended for flame stabilization. As will be outlined later on (see the Best Practice section) --- depending on the choice of the computational domain --- the inlet boundary conditions are crucial and not completely defined regarding the swirl intensity. Therefore isothermal RANS simulations have been used to identify the sensitivity related to this parameter. These simulations provide a cost effective way to gain knowledge for the more compuational expensive simulations dealing with the reacting case where additionally effects would overlap or compensate during such a parametric study.
In Fig. 5.2 results of the reacting case obtaind with ANSYS CFX are shown in comparison with experimental data. As mentioned befor the Turbulent Flame Speed Closure model has been used to treat the chemial reaction during the simulation where a strong sensitivity related to its model coefficient can be observed. Using the Zimont correlation for the turbulent flame speed with the standard value of A=0.5 overestimates the velocity a lot. This trend was also observed in further simulations where the formulation of Peters developed for the G-equation model was applied. The reason for these deviations becomes obvious if one considers the profiles of the axial velocity at 10mm. While at higher radii within the outer shear layer the agreement is promising (as observed in the isothermal simulations)a huge acceleration occurs at the flame front because of the density change of the reaction. The turbulent flame speed that determines the flame position is too large and as a result the combustion takes place too early. As shown above, in the isothermal case of this configuration, the k-epsilon model is generally able to predict the velocity field near the nozzle exit satisfactorily and so also the shear layer, which is important for the flame stabilization. Consequently, the deviations observed so far must result from the combustion model and therewith from the turbulent flame speed. In the literature strong variations in the right choice of the modeling constant entering the correlation for the turbulent flame speed have been reported depending on the configuration investigated. For that reason a series of simulations with a different model constant have been carried out and adequate agreement could be achieved with A=0.18 which is a strong deviation from its default vaulue A=0.5. The strong sensitivity indicates that this modeling approach is not well suited to help during the design process since even with the tuning of the flame position the scalar distribuation was poorly predicted.
LES Simulations of the Isothemal and Reacting Flow
In Fig. 5.3 the simulation results of the axial velocity for the isothermal and reacting case are compared with experimental data. Overall the spreading of the turbulent swirling flow and the size and intensity of the recirculation zone are well predicted and an excellent agreement with the measurements can be observed in all axial planes for both cases. Comparing the isothermal and reacting case a strong acceleration of the velocity caused by the thermal expansion through the flame front can be observed which is accurately reproduced by the simulation.
In Fig. 5.4 the temperature and methane mass fraction are compared with experimental data. Since the flame is very sharp the distribution is almost completely caused by the flame-turbulence interaction resulting in a wide turbulent flame brush. Except that a small radial offset exists and the gradient at higher axial positions is slightly over predicted the temperature is well predicted by the simulation. The mean profile of the methane mass fraction reflects its consumption by the chemical reaction on the left and the mixing with the coflowing air on the right. At the higher axial positions the methane mass fraction exceeds the measurement indicating insufficient fuel consumption in the simulation.
Comparing the efficiency functions (which is the main model parameter) the differences (marked with arrows in Fig. 5.3 and 5.4) between the two formulations are small and hence the model uncertainty will not influence the results crucially.
References
Contributors: Guido Kuenne (EKT), Andreas Dreizler (RSM), Johannes Janicka (EKT)
EKT: Institute of Energy and Power Plant Technology, Darmstadt University of Technology
RSM: Institute Reactive Flows and Diagnostics, Center of Smart Interfaces, Darmstadt University of Technology