Page 1 of 1

How face force is calculated in saturne

Posted: Thu Dec 21, 2017 5:20 pm
by lzhang
Hello,

I'd like to know how the face force is calculated in code_saturne (which can be obtained using field_get_val_v(iforbr, forbr)). I'm using the version 5.0.4, so where are the related subroutines located? And in my simulation, the face force is always 0 at the end of the first iteration, afterwards, it becomes non zero. Is there anything wrong?

And in saturne, the face force is calculated based on pressure, not the total pressure, am I right? Because the stress that I have in my simulation is very small.

Best regards,
Lei

Re: How face force is calculated in saturne

Posted: Fri Dec 22, 2017 12:07 am
by Yvan Fournier
Hello,

The force is updated in multiple places, to handle not only pressure but viscous effects (the pressure part being usually much stronger).

The best solution is to "grep" forbr in the source tree.

Best regards,

Yvan