Paper revision: add plots for shear-thinning case

[pkottmanjr]

For the revised article, we have decided to:

• remove different dependencies in the modeling argument
• remove power-law-like simulation
• add shear-thinning case

[pkottmanjr]

For reference: the profiles we want for numerical experiments - (formula resp. implementation):

• linear: $$f\Delta(u) = \overline f\Delta \cdot u$$ resp.

  viscosity[eps_, r_] := viscosity1 + viscosity2 * hEps[eps, r-R];

• power law illustrative: $$f\Delta(u) = \overline f\Delta \cdot \dfrac{\mathrm d u}{\mathrm d r}$$ resp.

  viscosity[eps_, r_] := viscosity1 + viscosity2 * D[hEps[eps, r-R], r];

• full power law: $$\eta_g = \frac{\tilde{\eta}C}{r^{2q}}\left|-\frac{3}{r\rho} + \frac{\rm{d}}{{\rm{d}}r}\left(\frac{1}{\rho}\right)\right|^q$$ resp.

  viscosity[eps_, r_] := viscosity1 * c/r^(2*q) * Abs[ -3/(r * specVol[eps, r]) + D[specVol[eps, r], r]]^q;

For the full power law, we want 2 separate cases:

• shear-thickening ($q > 1$)
• shear-thinning ($q < 1$)