[pbo]

one question about the way the viscous force is computed in turbulentForceAndTorqueFunctionObject.C:

// Viscous component
viscousForces[patchI] =
gSum
(
-nuEff.boundaryField()[curPatch]*
U.boundaryField()[curPatch].snGrad()*
mesh.magSf().boundaryField()[curPatch]
);

This is only valid in regions where the velocity field is more or less linear (snGrad being approximated by a one sided finite difference if I am not mistaken). In other words, for models with wall functions (log law), the above formulation will give incorrect viscous force estimates at the wall. Am I correct?