[matthias]

Hi Folks,

I have splitted the convection into different terms

U*U = U*U + B*Udelta*U + B*U*Udelta + B*B*uDelta*uDelta

where B is a volScalarField and uDelta a volVectorField. Both fields are calculated explicitly in a previous step.

A possible discretization could be

fvm::div(phi, U) + fvc::div( B*Udelta*U + B*U*Udelta + B*B*uDelta*uDelta)

Is there a way to discretize the remaining terms also implicitly?

(If I use fvm::div or fvm::Sp I'm getting an error reporting incompatible fields for operation)

Best regards

Matthias