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

cellLevel field not readable by topoSet

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

LinkBack Thread Tools Display Modes
Old   June 30, 2015, 04:35
Default cellLevel field not readable by topoSet
Join Date: Mar 2009
Location: Germany
Posts: 47
Rep Power: 8
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
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
New BoundaryCondition: how to modify a single component of a field? look86 OpenFOAM Programming & Development 4 August 28, 2012 16:48
Demagnetization field jan0385 ANSYS 0 April 12, 2012 09:09
Porosity field in Fluent wojciech FLUENT 1 September 20, 2010 11:19
problem with internal field and nonuniform list OFU OpenFOAM Running, Solving & CFD 1 October 5, 2009 03:35
Flow field measurement in turbomachinery! Wanlai Lin Main CFD Forum 3 September 3, 1999 12:06

All times are GMT -4. The time now is 16:00.