My simple view.
You could assign each slot a voltage vector, expressed in complex notation as:
V = exp(i*n* (2*Pi)*pp/N)
where
n = stator slot number, sequentially numbered by physical position from 1 to N
N = number of stator slots
pp = number of pole pairs (pp=1 for 2-pole motor).
Then add up the voltage vectors for all slots of a given phase using a spreadsheet.
Repeat for each phase.
Now you have three vectors, one per phase.
Perform symmetrical components transform on the three vectors.
I can post example spreadsheet of what I have in mind if you're interested.
But I'll bet you already knew that much. Understanding and interpreting the the finer points of the results requires more knowledge than I have.
fwiw, I'd think 3rd harmonic zero-sequence components (present even in balanced winding) may play a large role in evaluating circulating currents in delta configuration also.
=====================================
(2B)+(2B)' ?