Direct numerical simulation and subgrid analysis of a transitional droplet laden mixing layer
- 1 March 2000
- journal article
- research article
- Published by AIP Publishing in Physics of Fluids
- Vol. 12 (3), 650-671
- https://doi.org/10.1063/1.870271
Abstract
Direct numerical simulations of a temporally developing, droplet laden mixing layer undergoing transition to mixing turbulence are conducted. The formulation includes complete two-way couplings of mass, momentum, and energy. As many as 18×106 grid points are used to discretize the Eulerian gas phase equations and up to 5.7×106 initially polydisperse evaporating droplets are tracked in the Lagrangian reference frame. The complete transition to mixing turbulence is captured for several of the higher Reynolds number simulations and it is observed that increasing the droplet mass loading ratio results in a more “natural” turbulence characterized by increased rotational energy and less influence of the initial forcing perturbations. An increased mass loading also results in increased droplet organization within the layer. An a priori subgrid analysis is then conducted which shows that neglecting subgrid velocity fluctuations in the context of large eddy simulations may result in significant errors in predicting the droplet drag force for Stokes numbers St∼1 (with the flow time scale based on the mean velocity difference and initial vorticity thickness). Similar possible errors of lesser magnitude are also observed for the droplet heat flux and evaporation rate when thermodynamic subgrid fluctuations are neglected. An extension of the eddy interaction model commonly used in Reynolds-averaged simulations is then proposed in order to account for the missing subgrid information. Probability density functions (PDFs) of the subgrid fluctuations calculated across homogeneous planes are shown to be highly intermittent, particularly near the laminar–turbulent boundaries of the mixing layer. However, the actual subgrid PDFs calculated locally are much less intermittent and may be adequately modeled by the Gaussian distribution throughout the majority of the mixing layer. A scale similarity model is then employed to predict both the velocity and thermodynamic subgrid variances. The similarity model is well correlated with the actual subgrid variances and shows good agreement in predicting the local fluctuation intensities when a filter width-dependent model constant is used. The subgrid fluctuation variances acting on the droplets are then shown to be well modeled if the Eulerian subgrid variance model is interpolated to the droplet locations.This publication has 40 references indexed in Scilit:
- Direct numerical simulation of turbulence modulation by particles in isotropic turbulenceJournal of Fluid Mechanics, 1998
- Direct numerical simulations of evaporating droplet dispersion in forced low Mach number turbulenceInternational Journal of Heat and Mass Transfer, 1998
- Turbulent dispersion of particles using eddy interaction modelsInternational Journal of Multiphase Flow, 1996
- The structure and the small-scale intermittency of passive scalars in homogeneous turbulenceJournal of Scientific Computing, 1995
- Effect of Selective Modification of Turbulence on Two-Equation Models for Particle-Laden Turbulent FlowsJournal of Fluids Engineering, 1994
- Preferential concentration of particles by turbulenceInternational Journal of Multiphase Flow, 1994
- Several new numerical methods for compressible shear-layer simulationsApplied Numerical Mathematics, 1994
- On the two-way interaction between homogeneous turbulence and dispersed solid particles. I: Turbulence modificationPhysics of Fluids A: Fluid Dynamics, 1993
- Johnson-Edgeworth Translation for Probability Modeling of Binary Scalar Mixing in Turbulent FlowsCombustion Science and Technology, 1993
- Particle Dispersion by Vortex Structures in Plane Mixing LayersJournal of Fluids Engineering, 1992