[fisch]

Hi,

i want to add a (nondependent) source term to the UEqn of my simpleFoam solver.
I tried different ways where to place the term and the problem is that it lead to two different results and i don't know what is wrong...

way a) was:
tmp<fvVectorMatrix> UEqn
(
fvm::div(phi, U)
- fvm::laplacian(nu, U)
-sourceterm
);
eqnResidual = solve
(
UEqn() == -fvc::grad(p)
).initialResidual();


way b) was:
tmp<fvVectorMatrix> UEqn
(
fvm::div(phi, U)
- fvm::laplacian(nu, U)
);
eqnResidual = solve
(
UEqn() == -fvc::grad(p) + sourceterm
).initialResidual();


Can anybody help me please? What is the difference between these two ways?
Thanks a lot,
rupert

[l_r_mcglashan]

Can you post a complete case up? How do they differ?