SfePy does not handle units at all, so the units of the results are determined by units of the input parameters (material constants, loads, etc.). For >example, SI units on input -> SI units on output. Note that all the input units need to be consistent. Some consistent unit sets: http://www.eng->tips.com/viewthread.cfm?qid=355753
I do not understand what you mean by the strain rate, because the solution of your problem does not depend on time. Maybe you just need the >ev_cauchy_strain applied to the velocity u?
Try using stiffness_from_lame() (see examples/linear_elasticity/linear_elastic.py) - in linear isotropic elasticity, the stress is, given the Lame parameters >lambda, mu:
sigma = lambda tr(e) I + 2 mu e
and in incompressible fluid:
sigma = -p I + 2 mu e(velocity)
so I guess you could set lambda = 0, and mu to the viscosity to get the shear stress part using ev_cauchy_stress.
I guess my main question regarding strain rate is weather there are equations in place to solve for non-newtonian flow?