pressure corrector turbFoam
Dear all,
I have a question concerning the pressure correction in turbFoam: Code:
fvScalarMatrix pEqn This confuses me a little bit as I expected the right hand side to be fvc::div(H/A)=sumOverFaces(Hf/Af.Sf) and with phi=Uf.Sf and Uf=H/A at this point in the algorithm, I assume the code does the following: fvc::div(phi)=fvc::div(Hf/Af.Sf)=sumOverFaces(Hf/Af.Sf.Sf) with one Sf too many I probably overlook something, can someone tell me what? Thanks in advance, Sara |
thanks
Henrik,
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, Sara |
All times are GMT -4. The time now is 00:10. |