Brownian and turbulent diffusion
The phenomenon of small scale chaotic motion of particles through a fluid, due to the sub‑continuum interactions with individual fluid particles, is known as Brownian motion. In the case of a turbulent flow, chaotic motion of particles is also observed. It is due to the instantaneous forces arising from the turbulent, continuum-scale fluctuations of the velocity and pressure fields. This chaotic motion is estimated within the boundaries of an element during a particle time step.
A particle moves an average distance | r2 | in Δt according to , where Dp is the particle diffusivity:
- DT is the turbulent diffusivity associated with the unresolved turbulent motions.
- DB is the Brownian diffusivity associated with the sub-continuum effects.
The turbulent diffusivity is computed based on the turbulent eddy diffusivity, as follows:
The Brownian diffusivity is given by:
In these equations:
- μt is the eddy viscosity based on the turbulent model.
- ρf is the fluid density.
- μf is the fluid viscosity.
- kB is the Boltzmann constant.
- T is the fluid temperature.
- d is the diameter of the particle.
Therefore the perturbations of velocity field due to Brownian and turbulent forces is estimated as:
The presence of the particle inside the element is estimated as:
Substituting expressions for Brownian and turbulent diffusivities along with the expression for the perturbation of flow field, into the expression for the drag force, becomes:
where the unit vector that represents the direction of the velocity perturbation is generated randomly, and Ve is the element volume of the mesh.
By default, the flow solver computes the Brownian and turbulent diffusivity term . When you deactivated this term, the flow solver uses . Usually, you neglect this term when you apply the particle slip correction. For more information, see Particle slip correction factor.