[sharonyue]

Hi

In icoFoam we have:

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

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

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

Thanks

[sharonyue]

div(phiHbyA)?


Im confused about this divergence of a scalar thing.

[sharonyue]

I changed it to

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

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

Whats the problem..

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]

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.