[kbr]

Hi Khaled,

I think your fvm::Sp(rho, alpha)-Term should just look like 'rho*alpha' due to its explicit character. 'fvm' is used for source terms like fvm::Sp(rho,U) in any UEqn.

[booz]

Hey Kerstin,

my goal is to solve the steady-state navier-stokes-equation in simpleFoam with a volume sourceterm rho*alpha.

You mean it is enough to write

fvm::div(phi, U)
+ rho*alpha
+ turbulence->divR(U) );

?

Khaled

[kbr]

Hi Khaled,

yes I think so but I'm not sure. In my opinion the only way to discretize this term implicitly is

fvm::Sp(rho * alpha / U, U)

but I'm not sure if this is a good idea. Have a look in the source code of some turbulence models like kOmegaSST (epsEqn).

Kerstin

[booz]

Ok Kerstin,

thanks for the answer.

But last question: what means the term

fvm::Sp(rho * alpha / U, U)

exactly or how openFoam discretize this term?
I am a beginner and don´t understand openfoam really good.

Regards,
Khaled