/TaylorGreenParticles

Simulate motion of particles in Taylor Green Vortex with different Stokes number

Primary LanguagePython

TaylorGreenParticles

particles in taylor green vortex

Motions of particles in Taylor Green vortex with

TaylorGreenParticles is a single python script to simulate motions of particles in Taylor Green vortex. The animations can help to understand the phenomenon of preferential concentration of inertial particles in turbulence, see Preferential concentration. The fortran code is added by Guo Chen.


Solve the fluid phase

The analytical solution for 2D Taylor Green vortex is given by:


where u,v,\omega are respectively the x-velocity,y-velocity and vorticity. In this simulation, we set both x and y in range [0,2*pi]. And they are discretized into a PERIODIC mesh with shape of 64*64 (change it as you wish).

Solve the particle phase

The particle dispersed in fluid experiences a hydrodynamic force, and this force will drive the particle's motion following the Newton's Second Law. The procedure to caculate the force lists as follows.

  • First, caculate the particle Reynolds number re_p:


  • Next, caculate the drag coefficient drag_coff :


  • If we define the particle relaxation time as tau_p, we can get:


  • Finally, the particle translation motion (location) is governed by :


  • In practice, the particle velocity and location is updated by:


Results

The particle relaxation time (or particle inertia) can be changed by adjusting either dp or rho_p. Here we give five typical results.

  • If tau_p = 5.56, it means particle has small inertia, so it will exactly follow the motion of fluid.

particles in taylor green vortex

  • If we increase tau_p to 555.56, the particle will have similar inertia with the fluid. There will be a strong competition between the particle and fluid vortex. It turns out that the particle's power is much smaller than the vortex core, so they choose to stay in the edge of the vortexs.

particles in taylor green vortex

  • Further increase tau_p to 55555.56, the particles now can cross the outer part of the vortex, but they can still not penetrate into the vortex core.

particles in taylor green vortex

  • Finally, tau_p is increased to 277777.78, now they have much higher inertia than the vortex, so they just freely move in the whole domain. Particles win!

particles in taylor green vortex

Contact

Send me an email at zjuwangzhuo@zju.edu.cn (Zhuo Wang) or 11827039@zju.edu.cn (Guo Chen) if you have any suggestion.