[sfepy] Navier - Stokes at high viscosity

I am trying to model Stokes flow at very high viscosities (~10^20), from my results it looks like I cannot achieve this and am better suited using viscosities of ~10^5 or so. I know there is this older post [1] which states: "A word of warning though: SfePy does not have a dedicated Navier-Stokes solver (the Newton (+ and direct) solver use that is in the navier_stokes.py example is rather naive, and works for small problems with only, with viscosity "high enough")"

Has there been a different solver implemented since this was posted? If not, what do you recommend to work with such large viscosities?

Thank you

[1] "Problem with the mesh while trying to solve the Navier Stokes equations with sfepy"