[fabianpk]

I have been looking into the implementation of TVD schemes in OpenFOAM, and I found something that I'm not entirely sure I understand.

When the function for the limiter variable is calculated in LimitedScheme.C (so before it is actually calculated) it goes into LimitFuncs.C where it does magSqr if the variable that's being interpolated is a vectorField. So, later on when the gradient is calculated, it is not the gradient of the field but the gradient of the mag-squared field.

The idea I have is that this is done because the later calculations of the limiter etc requires a scalar. And since it doesn't really matter if we use magSqr or mag of the field, the one that's cheaper to calculate is used.

Is my theory to why magsqr is used correct? Or is it some entirely different reason?

/Fabian