CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Meshing & Mesh Conversion

[Other] cellLevel field not readable by topoSet

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

LinkBack Thread Tools Search this Thread Display Modes
Old   June 30, 2015, 05:35
Default cellLevel field not readable by topoSet
Join Date: Mar 2009
Location: Germany
Posts: 83
Rep Power: 15
blaise is on a distinguished road
Dear All,

i am facing a problem with topoSet when i try to use the field cellLevel to select cells. topoSet cannot find the cellLevel file. I need this because I want to do mesh refinement in dependence of the previous refinement level. This worked in OpenFOAM 2.2.x. Currently i am using build 2.3.x-2d875bf8e478.

Here is my topoSetDict:

        name    FS;
        type    cellSet;
        action  new;
        source boxToCell;
            box (-1000 -1000 -1.7) (1000 1000 1.7);
    // values of field within certain range
        name    FS;
        type    cellSet;
        action  subset;    
        source  fieldToCell;
            fieldName   cellLevel;   
            min         4;
            max         4;
This is the error message:

Create polyMesh for time = 66

Reading topoSetDict

Time = 66
    mesh not changed.
Created cellSet FS
    Applying source boxToCell
    Adding cells with center within boxes 1((-1000 -1000 -1.7) (1000 1000 1.7))
    cellSet FS now size 406636
Read set cellSet FS with size 406636
    Applying source fieldToCell
--> FOAM Warning : 
    From function fieldToCell::applyToSet(const topoSetSource::setAction, topoSet& set)
    in file sets/cellSources/fieldToCell/fieldToCell.C at line 203
    Cannot read field cellLevel from time 66
    cellSet FS now size 0

There is a cellLevel file in 66/polyMesh/, but i suppose this is not the place it should be, at least in comparison to version 2.2.x. In 2.2.x this file would sit at the top level of the respective time directory, and could even be read by paraFoam for display of cellLevels.

I thought this issue could depend on the debug setting in snappyHexMeshDict. Setting this to 2 is meant to give:
2 : write volScalarField with cellLevel for postprocessing
but gives strangely numbered time directories (like 66, probably intermediate mesh states).

So, how can topoSet get access to cellLevel information again?

Thank you!
blaise is offline   Reply With Quote


Thread Tools Search this Thread
Search this Thread:

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

Similar Threads
Thread Thread Starter Forum Replies Last Post
simpleDyMFoam on v1806 gamemakerh OpenFOAM Programming & Development 0 November 8, 2018 09:15
potential flows, helmholtz decomposition and other stuffs pigna Main CFD Forum 1 October 26, 2017 09:34
Access to field which is evaluated at the moment Tobi OpenFOAM Programming & Development 6 April 19, 2017 14:09
[General] How to create an additional vector with {Field 4, Field 5, Field 6} Bombacar ParaView 1 August 15, 2015 19:05
Demagnetization field jan0385 ANSYS 0 April 12, 2012 10:09

All times are GMT -4. The time now is 01:42.