View Issue Details

IDProjectCategoryView StatusLast Update
0001365OpenFOAMBugpublic2014-08-14 16:45
ReporterJuho Assigned Towill  
PrioritynormalSeveritymajorReproducibilityalways
Status resolvedResolutionfixed 
Summary0001365: twoPhaseEulerFoam: Wall lubrication force: Incorrect wall normals
DescriptionThe wall normal calculation employed in the twoPhaseEulerFoam wall lubrication force models produces incorrect wall normals near inlets, outlets and processor boundaries.

    volVectorField nWall(- fvc::grad(yWall_));
    nWall /= mag(nWall) + SMALL

See the attached "wallNormal.png" for an illustration of the issue and the effect it has on the simulated flow.
TagsNo tags attached.

Activities

Juho

2014-07-31 09:38

reporter  

wallNormals.png (35,749 bytes)   
wallNormals.png (35,749 bytes)   

Juho

2014-07-31 10:45

reporter   ~0003192

This issue can be fixed by replacing the gradient based wall normal calculation with a vector provided by wallDistReflection class. See the attached "modifiedWallNormals.png" for an illustration.

wallLubricationModel.H:

    // Wall normal unit vector
    wallDistReflection yr_;

Model.C:

    nWall = yr_.n();

Juho

2014-07-31 10:46

reporter  

modifiedWallNormals.png (28,260 bytes)   
modifiedWallNormals.png (28,260 bytes)   

will

2014-08-14 16:44

manager   ~0003207

Thanks for the report. This has been fixed in commit #173b252.

Issue History

Date Modified Username Field Change
2014-07-31 09:38 Juho New Issue
2014-07-31 09:38 Juho File Added: wallNormals.png
2014-07-31 10:45 Juho Note Added: 0003192
2014-07-31 10:46 Juho File Added: modifiedWallNormals.png
2014-08-14 16:44 will Note Added: 0003207
2014-08-14 16:44 will Status new => resolved
2014-08-14 16:44 will Resolution open => fixed
2014-08-14 16:44 will Assigned To => will