# Is that a checker-boarding pattern? How to avoid it?

 Register Blogs Members List Search Today's Posts Mark Forums Read

December 10, 2009, 07:02
Is that a checker-boarding pattern? How to avoid it?
#1
Member

Robert Castilla
Join Date: Apr 2009
Location: Spain
Posts: 94
Rep Power: 10
Hi,

I am simulating internal turbulent flow in complex geometries with simpleFoam and k-Omega SST, and I have found this strange behavior (see attached picture).

The fvSolution and fvSchemes files are:

solvers
{
p
{
solver GAMG;
tolerance 1e-7;
relTol 0.1;
// preconditioner DIC
smoother GaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
cacheAgglomeration on;
agglomerator faceAreaPair;
nCellsInCoarsestLevel 10;
mergeLevels 1;
};

U
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-8;
relTol 0.1;
nSweeps 1;
};

k
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-8;
relTol 0.1;
nSweeps 1;
};

omega
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-8;
relTol 0.1;
nSweeps 1;
};
}

SIMPLE
{
nNonOrthogonalCorrectors 0;
}

relaxationFactors
{
p 0.3;
U 0.7;
k 0.7;
epsilon 0.7;
}

-----------------------------------------------------------------------

ddtSchemes
{
}

{
default Gauss linear;
}

divSchemes
{
default none;
div(phi,U) Gauss linearUpwindV cellLimited Gauss linear 1;
div(phi,k) Gauss upwind;
div(phi,omega) Gauss upwind;
}

laplacianSchemes
{
default Gauss linear corrected;
// default Gauss linear limited 0.5;
// default Gauss linear limited 0.333;
}

interpolationSchemes
{
// default upwind phi;
// interpolate(U) upwind phi;
default linear;
interpolation linear;
}

{
default corrected;
}

fluxRequired
{
default no;
p;
}
-----------------------------------------------------------------------------------

I am not sure if it is the checker-boarding pattern. If it is, it should be avoided with SIMPLE, isn't it? Maybe have I to tune the SIMPLE parameters? Can somebody give me some hint?

Thank you

Robert
Attached Images
 press.jpg (36.6 KB, 55 views)

 December 10, 2009, 09:45 #2 Member   Robert Castilla Join Date: Apr 2009 Location: Spain Posts: 94 Rep Power: 10 I am trying with pisoFoam (transient solver) and results are not better.

 December 11, 2009, 12:50 #3 Member   Robert Castilla Join Date: Apr 2009 Location: Spain Posts: 94 Rep Power: 10 It seems so be related to mesh. It has been meshed with snappyHexMesh. If there are sharp transitions of cell size, this kind of patterns arise. I don't know the reason, but I have tried with a regular (not grading) hexahedral mesh, it went fine. Maybe it is something that only happens when meshing with snappyHexMesh inan internal flow.

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post sean Main CFD Forum 0 February 2, 2009 05:36 bala Siemens 5 March 6, 2008 05:23 Jan CFX 2 June 11, 2005 23:38 Selina Tracy Main CFD Forum 17 March 30, 2001 04:15 lee ji hyun Main CFD Forum 1 January 26, 2000 09:35