# A Population Balance Model for Large Eddy Simulation of Polydisperse Droplet Evolution

## TL;DR

This paper develops a population balance model integrated with LES to predict the evolution of polydisperse droplet size distributions in turbulent flows, accounting for breakup dynamics and validated against experimental data.

## Contribution

It introduces a novel LES-compatible population balance model for droplet breakup, extending kernels to include viscous effects and fitting parameters with laboratory experiments.

## Key findings

- Good agreement with experimental droplet size distributions
- Droplet size distribution variability is highly non-Gaussian
- Model captures intermittent behavior of key droplet properties

## Abstract

In the context of many applications of turbulent multi-phase flows, knowledge of the dispersed phase size distribution and its evolution is critical to predicting important macroscopic features. We use a population dynamics model for polydisperse droplet distributions specifically adapted to a LES framework including a model for droplet breakup due to turbulence, neglecting coalescence. Following earlier methods used in the Reynolds averaged Navier--Stokes framework, the droplet breakup due to turbulent fluctuations is modelled by treating droplet-eddy collisions as in kinetic theory of gases. In order to also model smaller droplets comparable to or smaller than the Kolmogorov scale we extend the breakup kernels using a structure function model that smoothly transitions from the inertial to the viscous range. The model includes a dimensionless coefficient that is fitted by comparing predictions in a one-dimensional version of the model with a laboratory experiment of oil droplet breakup below breaking waves. After initial comparisons of the one-dimensional model to measurements of oil droplets in an axisymmetric jet, it is then applied in a three-dimensional LES of a jet in crossflow with large oil droplets of a single size being released at the source of the jet. We model the concentration fields using $N_d =15$ bins of discrete droplet sizes and solve scalar transport equations for each bin. The resulting droplet size distributions are compared with published experimental data, and good agreement for the relative size distribution is obtained. The LES results also enable us to quantify size distribution variability. We find that the probability distribution functions of key quantities such as the total surface area and the Sauter mean diameter of oil droplets are highly variable, some displaying strong non-Gaussian intermittent behavior.