[ehsan_vaghefi]

Hi there

I am trying to use 'SimpleFoam' for a mesh where there is a known sink in each element.

I have read the 'simplefoam' guide on line and it says :
The first term of H{U} represents the matrix coefficients of the neighbouring cells multiplied by their velocity, while the second part contains the unsteady term and all the sources except the pressure gradient.

So I have the idea of adding my sinks as VolScalarField to the Ueqn :

tmp<fvvectormatrix> UEqn
(
fvm::div(phi, U)
+ turbulence ->divDevReff(U)
);

UEqn().relax();

eqnResidual = solve
(
UEqn() == -fvc::grad(p) + SINK
).initialResidual();

maxResidual = max(eqnResidual, maxResidual);

but it seems that no matter what I do, the answers don't change as if the SINK term is not there.

How can I solve this problem ?