September 2, 2011, 09:45
|
No convergence because of underdeterminedCells?
|
#1
|
New Member
Join Date: Aug 2011
Posts: 4
Rep Power: 14
|
Hi,
I have a problem running a case with SimpleFoam. The geometry is quite complex. I'm running the case with turbulence off. I made the mesh with netgen (1.4M cells, no layers). As you can see in the attached File the case doesn't converge although the solution looks physical to me.
My checkMesh complains only about some non-orthogonal Faces and a lot of under-determined Cells.
Could these cells effect the convergence?
Or what can I change in fvSolution / fvSchemes so that the case converge?
Thanks in advance.
Quote:
Checking geometry...
Overall domain bounding box (-876.472 -48 -63.9999) (884.528 48 273)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (2.12129e-17 -1.9514e-16 1.29017e-16) OK.
Max cell openness = 3.51972e-16 OK.
Max aspect ratio = 33.8371 OK.
Minumum face area = 0.146015. Maximum face area = 287.294. Face area magnitudes OK.
Min volume = 0.0386411. Max volume = 1541.7. Total volume = 3.29777e+07. Cell volumes OK.
Mesh non-orthogonality Max: 84.8534 average: 16.5116
*Number of severely non-orthogonal faces: 42.
Non-orthogonality check OK.
<<Writing 42 non-orthogonal faces to set nonOrthoFaces
Face pyramids OK.
Max skewness = 1.45798 OK.
Min/max edge length = 0.224133 28.1059 OK.
All angles in faces OK.
All face flatness OK.
Cell determinant (wellposedness) : minimum: 0 average: 1.36798
***Cells with small determinant found, number of cells: 10537
<<Writing 10537 under-determined cells to set underdeterminedCells
Failed 1 mesh checks.
End
|
Quote:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.0.0 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ddtSchemes
{
default steadyState;
}
gradSchemes
{
default Gauss linear;
}
divSchemes
{
default Gauss linear;
div(phi,U) Gauss linearUpwindV Gauss linear
div(phi,k) Gauss upwind;
div(phi,omega) Gauss upwind;
div((nuEff*dev(T(grad(U))))) Gauss linear;
div((nuEff*dev(grad(U).T()))) Gauss linear;
}
laplacianSchemes
{
default Gauss linear corrected;
}
interpolationSchemes
{
default linear;
}
snGradSchemes
{
default corrected;
}
fluxRequired
{
default no;
p;
}
// ************************************************** *********************** //
|
Quote:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.0.0 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
p
{
solver GAMG;
tolerance 1e-7;
relTol 0.1;
smoother GaussSeidel;
nPreSweeps 0;
nPostSweeps 2; //2
cacheAgglomeration on;
agglomerator faceAreaPair;
nCellsInCoarsestLevel 10;
mergeLevels 1;
}
U
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-8;
relTol 0.1;
nSweeps 1; //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;
}
potentialFlow
{
nNonOrthogonalCorrectors 10;
}
relaxationFactors
{
p 0.3;
U 0.7;
k 0.7;
omega 0.7;
}
cache
{
grad(U);
}
// ************************************************** *********************** //
|
|
|
|