Version 2022-09-27 - T.J.J.M. van Overveld
#---------------------#
# General information #
#---------------------#
Title:
The effect of the Stokes boundary layer on the dynamics of particle pairs in an oscillatory flow
Authors:
T.J.J.M. van Overveld, W.-P. Breugem, H.J.H. Clercx, M. Duran-Matute
Corresponding author:
M. Duran-Matute (m.duran.matute@tue.nl)
Fluids and Flows group and J.M. Burgers Center for Fluid Dynamics
Department of Applied Physics, Eindhoven University of Technology
P.O. Box 513, 5600 MB Eindhoven
The Netherlands
#----------------#
# Data structure #
#----------------#
Simulation naming convention:
Simulations are performed in two systems, namely the 'Oscillating box' and the 'Oscillating channel flow', which are sorted into the directories 'osc' and 'fix', respectively.
For each system, simulations are performed at different values of the particle Reynolds number (A \omega D / \nu), the reciprocal dimensionless boundary layer thickness (D / \delta), and density ratio (rho_s/rho_f). Each directory is numbered and contains the data of a single simulation, with parameters defined in the file 'key.txt'
Particle data:
pchar0001.out & pchar0002.out each corresponding to one of the two particles in a simulation
The files contain 17 columns corresponding to:
Time Step | Time | X Y Z Position | X Y Z Velocity | X Y Z Angular Velocity | X Y Z Hydrodynamic Force | X Y Z Collision Force
Flow field data in the xy plane:
xy_TimeStep_Height.txt (e.g. xy_0010000_008.txt)
with 'TimeStep' the time step in 7 digits at which the field is written
and 'Height' the height in cells in 3 digits at which the slice is taken.
The files contain 8 columns corresponding to:
X Y Cell Indices | X Y Z Velocity | Pressure | Solidity Fraction | Vorticity inside particle (for correction)
Flow field data in 3 dimensions:
The files for the 3D visualization of the flow fields are given in the directory 'Data/3DParaview', where the flowfields and lambda2-criterion visualization can be loaded in ParaView using the command 'paraview l2c.pvsm'.
#------------#
# Processing #
#------------#
All data processing is done in Jupyter-Lab. The notebooks are discussed per file hereunder.
Functions.ipynb
Contains regularly used functions
Preprocessing_{osc/fix}.ipynb
Converts the particle trajectories to derived quantities (Ap, Ar, AR, phase, Ag, Bg) which can then be stored in .txt files
Preprocessing_Circulation_{osc/fix}.ipynb
Calculates the 2D flow fields and particle trajectories into total and mean amounts of circulation (used in Figure7.ipynb).
Preprocessing_Rotation_fix.ipynb
Calculates the maximum rotation of the particles along the x-axis (used in Figure11.ipynb).
FigureX.ipynb
Uses the preprocessed data and functions to generate Figure X from the original paper.