How face force is calculated in saturne

[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

[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