Periodic boundaries for particle tracking

Understand how the flow solver relates the aproach and departure velocities for rotational and translational periodic boundaries.

In the flow solver, a particle that reaches a periodic boundary is presumed to pass through to the corresponding point on the matching periodic face. If is the location across the periodic boundary for any , then for the approach point obtained from integration of the particle motion equations, denoted by , the departure point is given by:

In addition, at a rotationally periodic boundary defined by a rotation of an angle Θ, about an axis directed along the unit vector ω and with the center of rotation x0, the departure velocity is obtained by a rotation imposing conservation of the cylindrical components of the velocity vector across the periodic boundary, according to:

where δij is the Kronecker delta tensor and εijk is the permutation tensor. In the case of a translational periodicity condition, the following condition is imposed:

The integration of the particle motion equations is then continued from the modified state defined by the imposed location with the imposed velocity .