|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
Join Date: Jun 2009
Location: Belgium
Posts: 3
Rep Power: 15 ![]() |
Dear all,
I have a question concerning the pressure correction in turbFoam: Code:
fvScalarMatrix pEqn ( fvm::laplacian(rUA, p) == fvc::div(phi) ); 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 Last edited by Sara D; July 9, 2009 at 04:26. |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Senior Member
Henrik Rusche
Join Date: Mar 2009
Location: Göttingen, Niedersachsen, Germany
Posts: 280
Rep Power: 16 ![]() |
Sara,
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: Quote:
Henrik |
|
![]() |
![]() |
![]() |
![]() |
#3 |
New Member
Join Date: Jun 2009
Location: Belgium
Posts: 3
Rep Power: 15 ![]() |
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 |
|
![]() |
![]() |
![]() |
Thread Tools | Search this Thread |
Display Modes | |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Neumann pressure BC and velocity field | Antech | Main CFD Forum | 0 | April 25, 2006 03:15 |
Gas pressure question | Dan Moskal | Main CFD Forum | 0 | October 24, 2002 23:02 |
what the result is negatif pressure at inlet | chong chee nan | FLUENT | 0 | December 29, 2001 06:13 |
Hydrostatic pressure in 2-phase flow modeling (CFX4.2) | HB &DS | CFX | 0 | January 9, 2000 14:19 |
Hydrostatic pressure in 2-phase flow modeling (long) | DS & HB | Main CFD Forum | 0 | January 8, 2000 16:00 |