I'm a little confused when it comes to the calculation of strain and shear rate. I know, what the matlab codes does
Gradient of u-velocity: [px1 py1] = gradient(u, hx, hy);
Gradient of v-velocity: [qx2 qy2] = gradient(v, hx, hy);
which gives the data for
px1 = du/dx, py1 = du/dy
qx2 = dv/dx, qy2 = dv/dy
To get the shear rate, pivlab calculates du/dy+dv/dx and for the strain rate du/dx-dv/dy. If i comapre the calculation of a counter-rotating vortex pair with the literature (eg. Steinberg and Driscoll, 2009 http://www.sciencedirect.com/science/article/pii/S0010218009001850), I get the same results with artificial image pairs only if i calculate du/dx - dv/dy. But they say, the strain rate is calculated using du/dy + dv/dx, which would be the shear rate in pivlab and gives totally different results. Where is my mental error?
If I compare the calculated vorticity via the curl-function and the gradient dv/dx-du/dy, I get the same results, as it should be the case, but strain and shear totally confuse me. Maybe that's why chemists shouldn't touch the field of fluid dynamics, but still... can anybody enlighten me?