Periodicity¶
The domain above has a periodicity on the side walls and from block 1 to block 2 on the suction side. Plot3D library can find periodic surfaces. Periodic surfaces can only be found by rotation about the x-axis. The user also needs to specify which direction is periodic either i, j, or k. Combinations of outer faces are searched to find the ones that are periodic when rotated about the x-axis and in the “k” direction (default option). Below is an example on how to find periodic surfaces.
In this example we will use the file PahtCascade-ASCII
1from plot3d import write_plot3D, read_plot3D, periodicity, connectivity
2
3blocks = read_plot3D('PahtCascade-ASCII.xyz', binary = False)
4face_matches, outer_faces_formatted = connectivity(blocks)
5periodic_surfaces, outer_faces_to_keep = periodicity(blocks,outer_faces,face_matches,periodic_direction='k',rotation_axis='x',nblades=55)
6# Append periodic surfaces to face_matches
7face_matches.extend(periodic_surfaces)
