View Issue Details

Jump to NotesJump to History
IDProjectCategoryView StatusDate SubmittedLast Update
0001183OpenFOAMBugpublic2014-02-25 02:152014-02-25 09:57
Reporteruser851Assigned Tohenry  
PrioritynormalSeveritymajorReproducibilityalways
Status resolvedResolutionno change required 
PlatformLinuxOSOpenSuseOS Version11.3
Summary0001183: phiHbyA in pEqn.H of bouyantPimpleFoam
DescriptionAccording to report from Japanese user of me, "surfaceScalarField phiHbyA" seems to have issue.
In the "pEqn.H", it is as follow.

    surfaceScalarField phiHbyA
    (
        "phiHbyA",
        (
            (fvc::interpolate(rho*U) & mesh.Sf())
          + rAUf*fvc::ddtCorr(rho, U, phi)
        )
      + phig
    );

But, Japanese user says that it should be as follow.

    surfaceScalarField phiHbyA
    (
        "phiHbyA",
        (
            (fvc::interpolate(rho*HbyA) & mesh.Sf())
          + rAUf*fvc::ddtCorr(rho, U, phi)
        )
Additional InformationFor checking it, he compared original and modification "pEqn.H" by hotRoom.
Here are residuals of these. You can see betterment on "p_rgh".
TagsNo tags attached.

Activities

user851

2014-02-25 02:15

 

residual.zip (20,422 bytes)

henry

2014-02-25 08:30

manager   ~0002893

It is not clear why Japanese user is suggenting the removal of the buoyancy term from the flux; this is an essential part of the algorithm.

henry

2014-02-25 09:57

manager   ~0002894

Sorry, I was confused by the missing phig and did not notice the change from U to HbyA.

Resolved by commit b766ecd3f957d8d31643963a28a06bec7100081d

Issue History

Date Modified Username Field Change
2014-02-25 02:15 user851 New Issue
2014-02-25 02:15 user851 File Added: residual.zip
2014-02-25 08:30 henry Note Added: 0002893
2014-02-25 08:30 henry Status new => closed
2014-02-25 08:30 henry Assigned To => henry
2014-02-25 08:30 henry Resolution open => no change required
2014-02-25 09:57 henry Note Added: 0002894
2014-02-25 09:57 henry Status closed => resolved