SnappyHexMesh Refinement nCellsBetweenLevels
 

Hi all,
I am trying to formulate effect of nCellsBetweenLevels parameter. I have meshed one case with different value.

I was thinking that when nCellsBetweenLevels=10 ( which i thought that is a very high value) there will be 10 cells between 2 levels of refinement in each direction. But as u can see from the screenshots its not like that at all. When i increase the value to 50, i get hardly 6 cells between two levels.

Does anyone can explain me what does the value of nCellsBetweenLevels exactly refers to?

nCellsBetweenLevels =10
Attachment 9028
nCellsBetweenLevels =50
Attachment 9029

Greetings to all!

@laymelek: Sorry for the reaaaally late reply. That looks like a bug that existed on that older version of OpenFOAM :(.

I'm answering to this thread so that questions related to this parameter "nCellsBetweenLevels" this can be discussed in a single thread.

At the request of the forum member CFDUser_, I've added a new section "nCellsBetweenLevels" at openfoamwiki.net in the page respective to snappyHexMesh: http://openfoamwiki.net/index.php/Sn...sBetweenLevels

If anyone has questions about this, feel free to ask!

Best regards,
Bruno


PS: In case anyone is wondering: the CFD-Online forums prefer that discussions on a particular topic are kept in a single thread.

Dear Bruno,

What happened when you set nCellsBetweenLayers to 0? can you upload that picture too?
I have some suggestions if you may consider them as valid point.
1. When you consider a case by default please specify OpenFOAM version so that one can understand which version can reproduce the same results.
2. Since I have observed the definition of some parameters changed with different versions, please give us Dict files so that one can have clear idea about other parameters too.

Thanks a lot for your time and efforts. :)

Best Regards,
CFDUser_

Dear CFDUser_,

Thanks for the feedback!

Ooops, I forgot to document the valid values. I've updated the subsection with an entry "Valid values": http://openfoamwiki.net/index.php/Sn...sBetweenLevels

Good reminder, I forgot to indicate which versions that analysis applies to. I've updated the subsection accordingly: http://openfoamwiki.net/index.php/Sn...sBetweenLevels

Strange, I thought that the wiki page already had an indication of where examples could be found. I've added a new section "Examples" on how they can be found: http://openfoamwiki.net/index.php/Sn...xMesh#Examples

Beyond this, feel free to edit the wiki page and improve it!

Best regards,
Bruno

Dear Bruno,

Good day :). Thanks once again for your explanations.
what is gapLevelIncrement in refine surfaces sub-dict?


Thanks & Regards,
CFDUser_

Hi CFDUser_,

If you document in that wiki page one or more parameter(s) with which you are already familiar with, I can try next weekend to figure out and document "gapLevelIncrement".

Best regards,
Bruno

Dear Bruno,

Coming 2-3 months are like hectic time for me. I will upload some of the valid cases in feature. My analysis are problem specific, so its not good to update in wiki pages. I will find out better way to post as I cant post my current cases as well.

Thanks & Regards,
CFDUser_

Hi CFDUser_,

Well, the reason why I used one of the cases from the snappyWiki was because it was a very simple case in which certain features can be tested.

It's very likely that if you extract only the problematic location of your geometry and create a base mesh only for that location, you'll be able to test a lot faster and a lot more possible workarounds for the problems you are getting with the whole geometry on that particular location.
Furthermore, those small test cases will then be so unrelated to the real geometry, that it will probably be almost impossible to deduce what it was from.

Beyond this, I'm quoting some feedback you sent me over PM back on August 26th, because this way it'll be easier for me to keep track of the things I still want to do in the near future, and because this way it's kept a public record of who discovered what:
It's the end of the week for me and my brain isn't very cooperative today :( I'll keep a note to look into this as soon as my brain is rested and able to figure this out ;)

Best regards,
Bruno

Hello everyone,
I am new in this forum, but I would still like to know, if there is a fix for this bug? I am having the exact same trouble meshing with SHM.

Quick request: Please be more specific to what bug you are referring to.I've taken a quick look at the posts above and there is no clear indication of any bug existing here.

Sorry. I seem to have the same problem as the OP.

I am using snappyHexMesh to mesh a simple geometry. When I set up the value nCellsBetweenLevels to 10, my resulting mesh does not have 10 cells between each level. When I set it up to 100 I get a significant change, but still not 100 cells, but someting like 12.
Could you please give me some advice?


It is almost like what OP has written:

Quick answer: OK... what about the instructions I wrote on the aforementioned wiki page section:
Did you check those settings?

First of all, thanks for your reply.
I read the instructions and checked my SHM-file. Neither point 2 and 3 is the reason for the error. I did not really understand what is meant by point 1
My case is a simple flow around a rectangular cylinder. The interflow is 5*D (where D is the height of the rectangular shape). The backflow is 15*D and the domainheight is 14*D. I have also a refinementBox and a refinementCylinder. The level of my rectangular surface is set to 4 and the refinementRegion is set to 2.


Here is a screenshot of the case. Although I set nCellsBetweenLevels to 100, the cells between levels are less than 20.
Attachment 66755

Quick answer: Sorry for the late reply, but here is what I've gotten with OpenFOAM 6 (up commit d8da4bd09212, Dec 4 21:32:25):

1. The attached "with20.png" image shows what with 20 cells per level got me.
2. The attached "with100.png" image shows what with 100 cells per level got me.
3. Attached is the case "testcase.tar.gz" that I used for the last run, i.e. 100.

What I saw in practice:

1. In the case with 20 cells, I think I saw only 19 cells, at least in one direction.
2. In the case with 100 cells, I used the ruler filter to measure that it had 15 metres along Z, which at 0.2m width, gave me only 75 cells and not 100.
   • However, it seemed to me that the limitation was triggered because it has stopped seeing more cells along Y and X, hence it stopping at 75 and not 100.

So, what can be done as a workaround:

• Use the distance based algorithm instead, e.g.:
  
  Code:

  ---

      refinementRegions
    {
        refinementBox
        {
            mode distance;
            levels ((1.0 2) (20.0 1));
        }
    }

  ---

  which is then shown on the attached image "distance_based.png". This sets a level 2 for 1m distance from the surface of the refinement box, along with a level 1 for a distance up to 20m.
• You can use any geometry as a reference for the refinement regions.