Remove up from URIS and Darcy/Brinkman residual contribution. Convergence improves and we think it is wrong to include it. But need to justify it theoretically

Author: aabrown100-git

Code/Source/solver/fluid.cpp

@@ -2238,12 +2238,12 @@ void fluid_3d_m(ComMod& com_mod, const int vmsFlag, const int eNoNw, const int e
   // Residual contribution Birkman term 
   // Local residue
   for (int a = 0; a < eNoNw; a++) {
-      lR(0,a) = lR(0,a) + mu*K_inverse_darcy_permeability*w*Nw(a)*(u[0]+up[0])
-                        + Res*DDir*w*Nw(a)*(u[0]+up[0]);
-      lR(1,a) = lR(1,a) + mu*K_inverse_darcy_permeability*w*Nw(a)*(u[1]+up[1])
-                        + Res*DDir*w*Nw(a)*(u[1]+up[1]);
-      lR(2,a) = lR(2,a) + mu*K_inverse_darcy_permeability*w*Nw(a)*(u[2]+up[2])
-                        + Res*DDir*w*Nw(a)*(u[2]+up[2]);
+      lR(0,a) = lR(0,a) + mu*K_inverse_darcy_permeability*w*Nw(a)*(u[0])
+                        + Res*DDir*w*Nw(a)*(u[0]);
+      lR(1,a) = lR(1,a) + mu*K_inverse_darcy_permeability*w*Nw(a)*(u[1])
+                        + Res*DDir*w*Nw(a)*(u[1]);
+      lR(2,a) = lR(2,a) + mu*K_inverse_darcy_permeability*w*Nw(a)*(u[2])
+                        + Res*DDir*w*Nw(a)*(u[2]);
   }
 
 }