|November 24, 2010, 10:41||
Bug in wall distance calculation
Join Date: Mar 2009
Location: Zürich, Switzerland
Posts: 63Rep Power: 10
I found a problem with the wall distance calculation in conjunction with cyclic patches in OpenFOAM 1.7.x and 1.5-dev. I will describe it for 1.5-dev as the rotationaAngle in 1.5-dev could help us to find the root of the problem.
Unfortunately I couldn't reproduce this bug with simpler meshes. The geometry is the inlet section (straight pipe) of my case with rotational periodicity (360°/7). I converted the mesh with fluentMeshToFoam from a fluent mesh file. It doesn't help to use fluent3DMeshTofoam. checkMesh doesn't report any problem. The settings of the cyclic patch are:
rotationCentre (0 0 0);
rotationAxis (0 0 1);
This setup gives me a "correct" Urel, k, omega and pressure but the result of the wall distance is wrong. I get the same results in 1.7.x and 1.5-dev.
If I change the rotationAngle to 51.42857, the simulation does not converge due to the wrong transformation of the velocity vectors at the cyclic patches. The big surprise is that the wall distance seems to be correct in this case. The attached image shows the wall distance y for a positive and negative rotation angle.
I could imagine that the problem can be found somewhere in FaceCellWave.C. I will try to find it but maybe somebody is faster than me. If anybody is interested I can send him the mesh.
|March 15, 2011, 10:32||
Join Date: Jan 2011
Posts: 73Rep Power: 8
I am using the files in yPlus.tar.gz from the thread http://www.cfd-online.com/Forums/ope...tml#post277386. This uses the nearWallDist.H file you are talking about in here. I had a look at the values of this d (nearWallDistance). I expected them to be the same all over the airfoil I am computing (since I am using a mesh with constant 1st cell height of 2.75e-5). To be sure that this is constant after having converted it to OF, I checked that in a post processor and this 1st cell height is constant, so there are no "conversion issues"). :S
What I get is that this value of nearWallDistance oscillates between 9e-6 and 1.3e-05. I don´t like this because the value of the computed yPlus then is not perfect.
Can you tell me anything about this?
|Thread||Thread Starter||Forum||Replies||Last Post|
|Wall Shadow interfering with Nusselt Number calculation?||ryzd||FLUENT||1||March 3, 2010 07:16|
|calculation of wall heat flux in CHT||richard||CFX||1||July 23, 2009 09:00|
|cell to wall distance||azmir||CD-adapco||0||March 27, 2007 01:03|
|MACRO FOR DIMENSIONLESS WALL DISTANCE Y+ ?||Wael||FLUENT||7||October 6, 2005 04:24|
|Wall functions||Abhijit Tilak||Main CFD Forum||6||February 5, 1999 02:16|