View Issue Details

IDProjectCategoryView StatusLast Update
0000519OpenFOAM[All Projects] Bugpublic2012-04-23 21:36
Reporteruser158Assigned Tohenry 
PrioritynormalSeveritytrivialReproducibilityalways
Status resolvedResolutionfixed 
PlatformLinuxOSOpenSuseOS Version12.1
Product Version 
Fixed in Version 
Summary0000519: code missing in nonNewtonianIcoFoam.C
Descriptionsmall code is missing in nonNewtonianIcoFoam.C

in nonNewtonianIcoFoam.C, replace...

        fvVectorMatrix UEqn
        (
            fvm::ddt(U)
          + fvm::div(phi, U)
          - fvm::laplacian(fluid.nu(), U)
        );

...with...

        fvVectorMatrix UEqn
        (
            fvm::ddt(U)
          + fvm::div(phi, U)
          - fvm::laplacian(fluid.nu(), U)
          - (fvc::grad(U) & fvc::grad(fluid.nu()))
        );

That is, the code "- (fvc::grad(U) & fvc::grad(fluid.nu()))" is missing in UEqn.
See for example interFoam solver.
Additional Information"- (fvc::grad(U) & fvc::grad(fluid.nu()))" must be added to get correct results
Have done several simulations for coaxial cylinders viscometer using Bingham viscoplastic fluid. The only time I can reproduce the correct result is when I have this term added. The correct analytical solution is attained by the Reiner-Riwlin equation.
TagsNo tags attached.

Activities

henry

2012-04-23 21:36

manager   ~0001283

Thanks for the bug-report.
Resolved by commit ab294fde8152d65bc572622692e0c55168692cab

Issue History

Date Modified Username Field Change
2012-04-23 18:28 user158 New Issue
2012-04-23 21:36 henry Note Added: 0001283
2012-04-23 21:36 henry Status new => resolved
2012-04-23 21:36 henry Resolution open => fixed
2012-04-23 21:36 henry Assigned To => henry