Saramito constitutive eqn in simple shear

[ktiqbal]

Hello,

I'm trying to simulate the Saramito EVP constituive eqn in simple shear (the same case as in Fig. 5a in [1]), however I am getting some discrepancy. I attach my input dicts for the problem and the solution that I get with rheoTool compared to the solution by Saramito. Please note that I am not trying to reproduce exactly the result in Fig. 5a in [1] as I am aware that Saramito did not include the factor of 1/sqrt(2) when computing the contraction of the deviator as rheoTool does; instead I am solving Eq. 13 in [1] with my own solver using the 2nd invariant of the deviator for the effective stress (as rheoTool does) and comparing with rheoTool's solution.
0.zip
constant.zip
system.zip

For reference, following the same procedure as in section 3.2 in [1], I get excellement agreement with rheoTool for other viscoelastic constitutive equations like Oldroyd-B and Giesekus (see below), but some discrepancy with the Saramito.. I have also solved the sysem of ODEs that the Saramito constitutive equation reduces to for simple shear with an imposed flow solution (Eq. 13 in [1]) (using Saramito's contraction norm) and can reproduce the solution presented in Fig. 5a.

Is the formulation in rheoTool by chance different to Saramito's formulation in [1] (other than the contraction vs 2nd invariant for the effective stress)? I looked through the equations for Saramito in the user guide and the source code Saramito.C and for n = 1 and k = etap, the equation seems should simplify to the same as in [1], but I'm stumped.

One specific question I have regarding this:

When calculating the deviator tau_d = tau - tr(tau)/N * I, am I correct in saying that N = 2 in rheoTool for 2D problems (and N = 3 for 3D problems)?

I'd be grateful for any help with this!

[1] Saramito, Pierre. "A new constitutive equation for elastoviscoplastic fluid flows." Journal of Non-Newtonian Fluid Mechanics 145.1 (2007)