A PDF propagation approach to model turbulent dispersion in swirling flows

Abstract

Listen

Translate article

A computationally efficient approach that solves for the spatial covariance matrix along the dense particle ensemble-averaged trajectory has been successfully applied to describe turbulent dispersion in swirling flows. The procedure to solve for the spatial covariance matrix is based on turbulence isotropy assumption, and it is analogous to Taylor's approach for turbulent dispersion. Unlike stochastic dispersion models, this approach does not involve computing a large number of individual particle trajectories in order to adequately represent the particle phase; a few representative particle ensembles are sufficient to describe turbulent dispersion. The particle Lagrangian properties required in this method are based on a previous study (Shirolkar and McQuay, 1998). The fluid phase information available from practical turbulence models is sufficient to estimate the time and length scales in the model. In this study, two different turbulence models are used to solve for the fluid phase – the standard k– ε model, and a multiple-time-scale (MTS) model. The models developed here are evaluated with the experiments of Sommerfeld and Qiu (1991). A direct comparison between the dispersion model developed in this study and a stochastic dispersion model based on the eddy lifetime concept is also provided. Estimates for the Reynolds stresses required in the stochastic model are obtained from a set of second-order algebraic relations. The results presented in the study demonstrate the computational efficiency of the present dispersion modeling approach. The results also show that the MTS model provides improved single-phase results in comparison to the k– ε model. The particle statistics, which are computed based on the fundamentals of the present approach, compare favorably with the experimental data. Furthermore, these statistics closely compare to those obtained using a stochastic dispersion model. Finally, the results indicate that the particle predictions are relatively unaffected by whether the Reynolds stresses are based on algebraic relations or on the turbulence isotropy assumption.