Model non-constant viscosity with fem::cell_quan