0001365: twoPhaseEulerFoam: Wall lubrication force: Incorrect wall normals

ID	Project	Category	View Status	Date Submitted	Last Update

0001365	OpenFOAM	Bug	public	2014-07-31 09:38	2014-08-14 16:45

Reporter	Juho	Assigned To	will
Priority	normal	Severity	major	Reproducibility	always
Status	resolved	Resolution	fixed

Summary	0001365: twoPhaseEulerFoam: Wall lubrication force: Incorrect wall normals
Description	The 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.
Tags	No tags attached.

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.

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