Hi! I want to add the Darcy-
 

Hi!
I want to add the Darcy-term in the simpleFoam solver. I tried this:

tmp<fvvectormatrix> UEqn
(
fvm::div(phi, U)
+ turbulence->divR(U)+fvm::SuSp(nu*G,U)
);

UEqn().relax();

solve(UEqn() == -fvc::grad(p));

Everything works fine until i changed my Field G from a volScalarField in a volTensorField, because i have directed permeabilities. I tried nu*(G & U), but it did not give the right solution.

Hope somebody can help me.
Thanks

Nico