Evaluation AC3-12

Front Page

Description

Test Data

CFD Simulations

Evaluation

Best Practice Advice

Particle-laden swirling flow

Application Challenge AC3-12   © copyright ERCOFTAC 2013

Comparison of Test Data and CFD

A rather good agreement between the experiments and predictions was obtained for gas and particle phase in both swirling cases considered. The comparison of the calculated streamlines of the gas flow with those obtained from the integration of the measured axial velocity shows that the flow field is predicted reasonably well for both conditions (Figure 10). The most obvious difference is that the axial extension of the central recirculation bubble is predicted to be larger at the top and downstream ends for both cases. The predicted width of the central recirculation bubble and the extension of the recirculation at the edge of the pipe expansion are in good agreement with the measured results.

Figure 10: Measured and calculated gas-phase streamlines (the upper parts of each figure corresponds to the calculations and the lower parts show the measurements); (a) Case 1; (b) Case 2.

The measured cross-sectional profiles of the three velocity components are compared with the calculations in Figure 11 for Case 2. The agreement is very good, except for the tangential velocity which is under-predicted in the region downstream of the location where the recirculation bubble has its largest radial extension. Although, the turbulent kinetic energy of the gas phase is considerably under- predicted in the initial mixing region between the primary and annular jets and within the recirculation at the edge of the pipe expansion (z = 52 mm), the agreement is reasonably good for the cross-sections further downstream. Similar results have been obtained for swirl Case 1 which was summarized in a previous publication (Sommerfeld et al. 1992).

Figure 11: Comparison between measurements and numerical calculations for the gas-phase in Case 2: (a) axial mean velocity; (b) radial mean velocity; (c) tangential mean velocity; (d) turbulent kinetic energy.

The measured and calculated particle mean velocities and the associated velocity fluctuations are compared in Figure 12 and 13 again for case 2. All mean velocity components are generally well-predicted, except for the radial velocity which is predicted to be positive at z = 315 mm (i.e. the particles move towards the wall) whereas the experiments show negative velocities. This implies that in the experiment the particles move on average away from the wall. Although, it might be expected that the numerical results in the near-wall region are very sensitive with respect to the modelling of the wall collision process, it was found that a variation in the normal restitution ratio and the friction coefficient in the three-dimensional inelastic collision model did not result in considerable changes in the velocity profiles. The scattering of the numerical results in the core region downstream of the location, where the central recirculation bubble has its largest radial extension is a result of the low number of collected particles, since the majority of the particles have already moved out of the core region.

The particle velocity fluctuations are more or less under-predicted for all three velocity components, which presumably is caused by the under- prediction of the gas-phase turbulent kinetic energy in the initial region of flow development (Figure 11(d)). This is demonstrated in Figure 12(b) for the fluctuation of the axial particle velocity component.

Figure 12: Comparison of measured and calculated axial particle velocity profiles: (a) mean particle velocity (m/s); (b) mean particle velocity fluctuation (m/s), Case 2.

Figure 13: Comparison of measurements and calculations: (a) radial particle velocity (m/s); (b) tangential particle velocity (m/s), Case 2.

The calculated particle mass flux shows very good agreement with the experiments in the initial region just downstream of the inlet, namely at z = 52 mm, and further downstream at z = 315 (Figure 14(a)). In the intermediate region, where the particles rebound from the wall and interact again with the central recirculation bubble (see Figure 7), some differences between the experiments and predictions are observed in the near-wall region. At z = 155 and 195 mm, the measurements show a relatively narrow layer of particles near the wall with the maximum being directly at the wall. The calculations show a slightly wider particle layer with a maximum in the particle mass flux located at a radial position of 65 mm, which might result from the assumptions made in the particle-wall collision model.

A sensitive study, by using different normal restitution ratios and wall friction coefficients in the particle-wall collision model, shows that their influence on the particle mass flux distribution at z = 155 and 195 mm is not too strong (Figure 15). A reduction in the normal restitution ratio from 0.8 to 0.6 gives a slightly narrower particle layer with a more pronounced maximum closer to the wall. When the particle lift force is switched off in the calculations, the particle layer becomes wider with a less pronounced maximum in the mass flux distribution near the wall. This also shows that the particle lift force due to rotation is only important in the near-wall region, where high rotational velocities are induced by the wall collision. Since the particles lag behind the air flow in the near-wall region, the direction of the lift force is directed towards the centreline. From these results, one may conclude that the particle motion in the near- wall region might also be affected by electrostatic forces, although the Plexiglas test section was carefully grounded.

Figure 14: Comparison of measurements and calculations: (a) stream-wise particle mass flux (kg/m2s); (b) particle number mean diameter (μm), Case 2.

In the publications of Sommerfeld et al. (1992) additional numerical results for the swirl Case 1 may be found. Details on measured and calculated particle size-velocity correlations are presented in Sommerfeld and Qiu (1993).

Figure 15: Measured and calculated axial particle mass flux at z  = 115 and 195 mm for different parameters in the wall collision model: △, with lift force, e = 0.8, μ = 0.3; ◊, with lift force, e = 0.6, μ = 0.3; ✱, without lift force, e = 0.8, μ = 0.3; □, experiment, Case 2.

---

Contributed by: Martin Sommerfeld — Martin-Luther-Universitat Halle-Wittenberg

Front Page

Description

Test Data

CFD Simulations

Evaluation

Best Practice Advice

© copyright ERCOFTAC 2013