-
in mpfa schemes, k*gradh is recovered on neumann boundaries by dividing the neumann flux to be applied by the ratio density/viscosity. On neumann boundaries we thus have to make sure that this ratio is multiplied onto the calculated fluxes in the flux variables.
a6104213