pressure corrector turbFoam
I have a question concerning the pressure correction in turbFoam:
This confuses me a little bit as I expected the right hand side to be
and with phi=Uf.Sf and Uf=H/A at this point in the algorithm, I assume the code does the following:
with one Sf too many
I probably overlook something, can someone tell me what?
Thanks in advance,
I think what you are missing is that fvc::div(phi) DOES a simple sum over all cell faces and DOES NEITHER dot-product with the face area vectors NOR interpolate to the face.
Yes, this is not intuitive, but it makes a lot of sense when you think about it.
Unfortunately, on page P-39 the Programmer's Guide says:
Thank you, now I understand.
That is indeed a pitfall and the programmers guide doesn't help in this case.
Thanks a lot for clarifying,
|All times are GMT -4. The time now is 06:51.|