CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Native Meshers: blockMesh

cyclic && Non-orthogonality !!

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 27, 2013, 09:08
Angry cyclic && Non-orthogonality !!
  #1
Member
 
ali jafari
Join Date: Sep 2012
Posts: 50
Rep Power: 4
ali jafari is on a distinguished road
hi

I have flow between two cylinders their top and bottom are cyclic(see attached file). after >blockMesh "OKmesh" display on terminal , but >checkMesh display a fail mesh (nonorthogonality face ) when I see these faces in paraview I find cyclic boundary !!
I check orientations of faces and other blockmesh rules but have not effect
Note : this error occurs in OF 2.2.0 , checkmesh in other version dose not display any error!
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

//convertToMeters 0.1;

vertices
(
(0 0.5 0) //0
//inner , bottom circles
(0 0.5 0) //1
(0.5 0 0) //2
(0 -0.5 0) //3
(-0.5 0 0) //4
//outer , bottom circles
(0 0.566 0) //5
(0.566 0 0) //6
(0 -0.566 0) //7
(-0.566 0 0) //8
//inner , top circles
(0 0.5 0.66) //9
(0.5 0 0.66) //10
(0 -0.5 0.66) //11
(-0.5 0 0.66) //12
//outer , top circles
(0 0.566 0.66) //13
(0.566 0 0.66) //14
(0 -0.566 0.66) //15
(-0.566 0 0.66) //16


);

blocks
(
hex (9 13 14 10 1 5 6 2) (35 104 90) simpleGrading (1 1 1) //first 1/4
hex (9 12 16 13 1 4 8 5) (104 35 90) simpleGrading (1 1 1) //2nd 1/4
hex (12 11 15 16 4 3 7 8) (104 35 90) simpleGrading (1 1 1) //3th 1/4
hex (10 14 15 11 2 6 7 3) (35 104 90) simpleGrading (1 1 1) //4th 1/4

//hex (1 5 7 3 9 11 15 13) (4 4 4) simpleGrading (1 1 1) //right half
//hex (1 3 7 5 9 13 15 11) (4 4 4) simpleGrading (1 1 1)
);

edges
(
// in - bottom circle
arc 2 1 (0.3428 0.3639 0)
arc 2 3 (0.3428 -0.3639 0)
arc 3 4 (-0.3428 -0.3639 0)
arc 4 1 (-0.3428 0.3639 0)

// out - bottom circle
arc 5 6 (0.4174 0.3823 0)
arc 6 7 (0.4174 -0.3823 0)
arc 7 8 (-0.4174 -0.3823 0)
arc 8 5 (-0.4174 0.3823 0)

// in - top circle
arc 9 10 (0.3428 0.3639 0.66)
arc 10 11 (0.3428 -0.3639 0.66)
arc 11 12 (-0.3428 -0.3639 0.66)
arc 12 9 (-0.3428 0.3639 0.66)

// out - top circle
arc 13 14 (0.4174 0.3823 0.66)
arc 14 15 (0.4174 -0.3823 0.66)
arc 15 16 (-0.4174 -0.3823 0.66)
arc 16 13 (-0.4174 0.3823 0.66)

);

boundary
(
out
{
type wall;
faces
(
(13 14 6 5)
(13 16 8 5)
(14 15 7 6)
(15 16 8 7)
);
}

in
{
type wall;
faces
(
(1 2 10 9)
(1 4 12 9)
(10 11 3 2)
(11 12 4 3)
);
}

top
{
type cyclic;
neighbourPatch bottom;

faces
(
(16 12 9 13)
(16 15 11 12)
(11 15 14 10)
(14 13 9 10)


);

}

bottom
{
type cyclic;
neighbourPatch top;

faces
(
(8 4 1 5)
(8 7 3 4)
(3 7 6 2)
(6 5 1 2)

);
}

);

mergePatchPairs
(

);

// ************************************************** *********************** //
ali jafari is offline   Reply With Quote

Old   May 27, 2013, 09:56
Default
  #2
Member
 
Yosmcer Mocktai
Join Date: Apr 2013
Location: Behind a computer
Posts: 41
Rep Power: 8
Yosmcer will become famous soon enough
Hello,

Why did you put two points (0 and 1) at the same place?
You do not use the 0 one.

I launched your mesh with
Code:
paraFoam -block
I attached a screenshot of the result. The block are represented by trapezoidal one, this is the way the -block function work, in reality, they follow the edges.

The first 4 blocks are well defined (if the 0 is not considered), but I see a problem with the two last:
//hex (1 5 7 3 9 11 15 13) (4 4 4) simpleGrading (1 1 1) //right half
//hex (1 3 7 5 9 13 15 11) (4 4 4) simpleGrading (1 1 1)

As I can see, all these points are coplanar, so cannot constitute any block.

[EDIT]Sorry, they cannot be the problem as they are commented[/EDIT]

Yosmcer
Attached Images
File Type: jpg Blocks.jpg (34.7 KB, 28 views)

Last edited by Yosmcer; May 27, 2013 at 12:59.
Yosmcer is offline   Reply With Quote

Old   May 27, 2013, 12:06
Default
  #3
Member
 
ali jafari
Join Date: Sep 2012
Posts: 50
Rep Power: 4
ali jafari is on a distinguished road
thank you for quick reply but two last dont important because these just note , you can delete them.

please tell me your idea clearly , thank you
ali jafari is offline   Reply With Quote

Old   May 27, 2013, 12:57
Default
  #4
Member
 
Yosmcer Mocktai
Join Date: Apr 2013
Location: Behind a computer
Posts: 41
Rep Power: 8
Yosmcer will become famous soon enough
Quote:
Originally Posted by ali jafari View Post
thank you for quick reply but two last dont important because these just note , you can delete them.
Oh, yes, I was distracted, they are commented, so they cannot be the problem.
Sorry.

[EDIT]

I have the same error with checkMesh:
Code:
Checking geometry...
    Overall domain bounding box (-0.566 -0.566 0) (0.566 0.566 0.66)
    Mesh (non-empty, non-wedge) directions (1 1 1)
    Mesh (non-empty) directions (1 1 1)
    Boundary openness (-7.75343e-16 1.43147e-15 -5.31732e-17) OK.
    Max cell openness = 3.24349e-16 OK.
    Max aspect ratio = 4.38502 OK.
    Minimum face area = 1.38286e-05. Maximum face area = 6.26916e-05.  Face area magnitudes OK.
    Min volume = 1.04619e-07. Max volume = 1.18166e-07.  Total volume = 0.145991.  Cell volumes OK.
    Mesh non-orthogonality Max: 180 average: 7.01164
 ***Number of non-orthogonality errors: 29120.
  <<Writing 29120 non-orthogonal faces to set nonOrthoFaces
    Face pyramids OK.
    Max skewness = 0.019251 OK.
    Coupled point location match (average 0) OK.

Failed 1 mesh checks.
But once I defined cyclic boundary fields for the top and bottom patches in the U and p files, paraFoam runs stressfully without errors.

Code:
(...)
    top
    {
        type            cyclic;
    }

    bottom
    {
        type            cyclic;
    }
(...)
I put all the case I used to sample this blockMeshDict file in attachment.
Attached Images
File Type: jpg Tube.jpg (61.3 KB, 14 views)
Attached Files
File Type: zip Ali.zip (77.8 KB, 5 views)

Last edited by Yosmcer; May 28, 2013 at 04:02.
Yosmcer is offline   Reply With Quote

Old   May 31, 2013, 11:56
Default
  #5
Member
 
Yosmcer Mocktai
Join Date: Apr 2013
Location: Behind a computer
Posts: 41
Rep Power: 8
Yosmcer will become famous soon enough
Is your problem solved?

If I change the patch type to patch instead of cyclic, there is no more non-orthogonality errors.
If I look to the number of non-orthogonality errors, it is equal to 104*35*8 so proportional to two meshing dimensions, but not to the third (90).

However, this error do not prevent you to run the case (I tested with icoFoam, this is just slow as you have a big amount of cells).

[EDIT]

This problem seems to be already known and look solved in the next OpenFOAM version:

CheckMesh fails when cyclic boundaries are defined - OF 2.2.0

Last edited by Yosmcer; May 31, 2013 at 12:56.
Yosmcer is offline   Reply With Quote

Old   May 31, 2013, 13:49
Default
  #6
Member
 
ali jafari
Join Date: Sep 2012
Posts: 50
Rep Power: 4
ali jafari is on a distinguished road
thanks for your contribution .
yes , it seems a bug ! I couldn't solve this problem. I changed nCorrector and NnonorthogonalCorrection , but it just was useful for icoFoam and and it didn't work for steady-state solver-problem (simpleFoam) . my problem is steady-state and I'm confused....icoFoam very slow for my case and dont converge to stationary solution.

PLEASE HELP OpenFoam EXPERT !!!!!!!!!!
ali jafari is offline   Reply With Quote

Old   May 31, 2013, 14:39
Default
  #7
Member
 
Yosmcer Mocktai
Join Date: Apr 2013
Location: Behind a computer
Posts: 41
Rep Power: 8
Yosmcer will become famous soon enough
Quote:
Originally Posted by ali jafari View Post
My problem is steady-state and I'm confused....icoFoam very slow for my case and dont converge to stationary solution.

PLEASE HELP OpenFoam EXPERT !!!!!!!!!!
I'm unfortunately not an expert. I noticed that it is slow to run, but I think it is because of the number of nodes. The more node you have, the more time it will take.

Maybe you could try with a less accurate mesh (eg: (5 45 52) instead of (35 90 104)) for the tests, and put the good mesh when you want the best results.

I never used simpleFoam cause icoFoam was asked in the work I had to do, so I cannot help with it.

For icoFoam, pay attention to the Courant Number. Max Courant Number must be under 1 to converge.
You can change this number by changing the mesh resolution and the time step.
It is displayed at each iteration.
See point 2.1.6.2 in tutorial.

I think icoFoam is enough for steady-state, but simpleFoam should also work.

If you think your problem is related to the mesh (But it seams that this IS OK as the non-orthogonality errors seems to not cause problem and just a bug), the tread could be continued here.
If you think the problem is about using the solver (which I think is), you may ask for help in a solver part of the forum, you will have more expert on this.
If you do so, put a link here and say if solved.

Hope you will find the solution.
Yosmcer is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Cyclic boundary conditions concerning channel with constant flow rate pascool OpenFOAM Pre-Processing 1 November 16, 2012 14:04
Cyclic Mesh Error vbchris OpenFOAM Meshing & Mesh Conversion 4 November 15, 2012 14:04
Cyclic bc. and overlapGgi Pekka OpenFOAM 12 February 15, 2012 08:10
Problem with cyclic patches sven OpenFOAM Running, Solving & CFD 0 December 5, 2011 15:27
Pressure instability with rhoSimpleFoam daniel_mills OpenFOAM Running, Solving & CFD 44 February 17, 2011 18:08


All times are GMT -4. The time now is 02:37.