CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (http://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   icoFoam's pressure possion equation. (http://www.cfd-online.com/Forums/openfoam-programming-development/120786-icofoams-pressure-possion-equation.html)

sharonyue July 15, 2013 03:49

icoFoam's pressure possion equation.
 
1 Attachment(s)
Hi

In icoFoam we have:
Code:

fvScalarMatrix pEqn
                (
                    fvm::laplacian(rAU, p) == fvc::div(phiHbyA)
                );

But via my equation(see the image), Shouldnt it be

Code:

fvScalarMatrix pEqn
                (
                    fvm::laplacian(rAU, p) == fvc::div(HbyA)
                );

Thanks

sharonyue July 16, 2013 03:30

div(phiHbyA)?
\nabla \cdot a

Im confused about this divergence of a scalar thing.

sharonyue July 16, 2013 21:34

I changed it to
Code:

fvScalarMatrix pEqn
                (
                    fvm::laplacian(rAU, p) == fvc::div(HbyA)
                );

It compiles successfully and works in cavity and elbow case!And the solution is almost the same with
Code:

fvScalarMatrix pEqn
                (
                    fvm::laplacian(rAU, p) == fvc::div(PhiHbyA)
                );

Whats the problem..

Code:

fvc::div(HbyA):dimensions      [0 0 -1 0 0 0 0];

internalField  nonuniform List<scalar> 9(0.157469 -0.0674757 -0.487085 -0.960691 0.789192 0.715649 30.1762 0.258981 -30.5822);

fvc::div(phiHbyA):dimensions      [0 0 -1 0 0 0 0];

internalField  nonuniform List<scalar> 9(-0.435688 -0.0982547 0.152804 -0.0448476 0.720155 -0.152663 28.2072 0.236062 -28.5848);


santiagomarquezd July 26, 2013 12:39

Hi, check the p-U coupling section (2.7) in my Ph.D thesis. It will also give you some pointers to good material on this topic.

Link to the PhD thesis is provided in my user page at openfoamwiki.net: http://openfoamwiki.net/index.php/User:Santiagomarquezd

Regards.


All times are GMT -4. The time now is 21:44.