In this paper, the Eulerian stochastic fields (ESF) method in a LES framework is applied to a generic selective catalytic reduction (SCR) configuration in order to retrieve seamlessly the effect of turbulence–chemistry interaction on the NH3-conversion and AdBlue film formation. The ESF method is based on the transport equation of a joint scalar filtered density function. The injection of AdBlue and relevant spray dynamics are modeled by the Lagrange-particle description where a computational parcel represents a cluster of real AdBlue droplets. Owing to the low injection pressure, weak atomization of AdBlue leads to intense droplet wall-impingement and film formation on the SCR duct wall and/or on the mixer elements, if present. A 2-D thin film approach is adopted to model the film formation and dynamics in combination with multi-regimes droplet–wall interaction model. After verification tasks against the reference data for simple gas phase reactive cases, the verified ESF method is coupled with a one-equation based LES turbulence model together with the Lagrange particle-tracking which is combined with the 2-D thin film approach to simulate the AdBlue injection and the film formation in the generic SCR configuration. At first, the assessment of the adopted LES mesh is carried out in terms of the so-called LES quality of index to determine the optimal mesh resolution. Next, to ensure the convergence with respect to the number of required ESF and mesh resolution, SCR simulations are performed by using various (2, 4, 6, 8, 12 and 16) stochastic fields and also with refined mesh. Thereby, the sensitivity of the predictive capability of the numerical tool is evaluated for the production of NH3 and HNCO with and without ESFs. This clearly indicates the importance of an accurate description of the turbulence–chemistry interaction to retrieve reliably the conversion of NH3. Finally, the film dynamics especially the evolution of the film thickness described with six required ESF and without ESF is compared with experimental data from the generic SCR configurations. To further demonstrate the potential of suggested numerical approach, a detailed numerical analysis is provided in terms of spray-impingement dynamics, scalar uniformity index, droplet life time and characteristics of the droplets prone to form solid deposits in the monolith channel.
Read full abstract