We have to modify a subroutine to be able to use a non Newtonian fluid. There are several things to do.
At the end, we have to modify a Fortran program, we add the following lines.
We compare the velocity gradient (gam) with a limit value (gamc). Depending on the actual value of gam, the fluid has a different viscosity.
The constants are determined using the theoretical relation given by our supervisor.