The code hosted on this repository generates an animation of a rigid disc that rotates passively under channel-flow turbulence.
The solver Incompact3D is used to perform a Direct Numerical Simulation of a turbulent channel fitted with a free-to-move solid disc on the wall.
A collection of two-dimensional snapshots of the shear-stress field at the wall are stored in the directory flow_snapshots and used as input to the script animation_disc.py.
The Matplotlib library is then used to render the flow field snapshots into animation frames and export the resulting video to file.
The video below is an example of the output, obtained using a larger range of flow snapshots from than those included in flow_snapshots:
More complete information on the flow physics and the numerical methods can be found in:
Olivucci P, Wise DJ, Ricco P, Reduction of turbulent skin-friction drag by passively rotating discs, J. Fluid Mech., 2021.