On random walk models for simulation of particle-laden turbulent flows

Abstract

In this investigation, the accuracy of the discrete and continuous random walk (DRW, CRW) stochastic models for simulation of fluid (material) point particle, as well as inertial and Brownian particles, was studied. The corresponding dispersion, concentration, and deposition of suspended micro- and nano-particles in turbulent flows were analyzed. First, the DRW model used in the ANSYS-Fluent commercial CFD code for generating instantaneous flow fluctuations in inhomogeneous turbulent flows was evaluated. For this purpose, turbulent flows in a channel were simulated using a Reynolds-averaged Navier–Stokes (RANS) approach in conjunction with the Reynolds Stress Transport turbulence model (RSTM). Then spherical particles with diameters in the range of 30 µm to 10 nm were introduced uniformly in the channel. Under the assumption of one-way coupling, ensembles of particle trajectories for different sizes were generated by solving the particle equation of motion, including the drag and Brownian forces. The DRW stochastic turbulence model of the software was used to include the effects of instantaneous velocity fluctuations on particle motion, and the steady state concentration distribution and deposition velocity of particles of various sizes were evaluated. In addition, the improved CRW model based on the normalized Langevin equation was used in an in-house Matlab code. Comparisons of the predicted results of the DRW model of ANSYS-Fluent with the available experimental data and the DNS simulation results and empirical predictions showed that this model is not able to accurately predict the flow fluctuations seen by the particles in that it leads to unreasonable concentration profiles and time-varying deposition velocities. However, the predictions of the improved CRW model were in good agreement with the experimental data and the DNS results. Possible reasons causing the discrepancies between the DRW predictions and the experimental data were discussed. The improved CRW model was also implemented through user-defined functions into the ANSYS-Fluent code, which resulted in accurate concentration distribution and deposition velocity for different size particles.