Calculation of shear modulus for 2D, dense, bi-disperse, soft sphere granular mixture

On making an affine shear transformation of a square cell of granular, bi-disperse particles at equilibrium we get non-affine displacements to preserve equilibrium.

I divided the domain into many pieces. Below is an example of the domain divided into 4 pieces and their individual non-affine corrections plotted together for convenience.



In the following picture, the domain is divided into 400 pieces and the non-affine correction to the shear modulus is computed for each sub-domain and plotted.


NOTES: (i) these are computed using the formulas using the formulas in Craig's PRE of 2006. (ii) there were 2 critical bugs in the coding the potential & I had confused the diameter as the radius making the system super dense. (iii) the code that does this stuff is sitting on research-1 at:
/second_drive/asad/granularShearData/