CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Post-Processing

sampling on cellSet

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 22, 2013, 02:44
Default sampling on cellSet
  #1
Senior Member
 
Tarak
Join Date: Aug 2010
Location: State College, PA
Posts: 105
Rep Power: 6
Tarak is on a distinguished road
Hii,

I need to use the "sample" utility to sample field values in a cellSet. Can someone please tell me the entries needed in sampleDict to do this?

Thanks
Tarak is offline   Reply With Quote

Old   May 22, 2013, 03:47
Default
  #2
Senior Member
 
Tushar Chourushi
Join Date: Jul 2009
Location: IIT-Indore, India
Posts: 318
Blog Entries: 1
Rep Power: 9
Tushar@cfd is on a distinguished road
Quote:
Originally Posted by Tarak View Post
Hii,

I need to use the "sample" utility to sample field values in a cellSet. Can someone please tell me the entries needed in sampleDict to do this?

Thanks
I think the following link will solve your problem..

SampleDict
Tushar@cfd is offline   Reply With Quote

Old   May 22, 2013, 11:13
Default
  #3
Senior Member
 
Tarak
Join Date: Aug 2010
Location: State College, PA
Posts: 105
Rep Power: 6
Tarak is on a distinguished road
Hii Tushar,

Thanks for your reply. If I understand correctly, the link talks about sampling on a patch, which I can do by writing something like
surfaces
(
airfoil
{
type patch;
patches (airfoil);;
// interpolate false;
}
);

inside sampleDict.
But I need to sample data along a cellSet (not patch)
Tarak is offline   Reply With Quote

Old   May 23, 2013, 00:33
Default
  #4
Senior Member
 
Tushar Chourushi
Join Date: Jul 2009
Location: IIT-Indore, India
Posts: 318
Blog Entries: 1
Rep Power: 9
Tushar@cfd is on a distinguished road
Quote:
Originally Posted by Tarak View Post
Hii Tushar,

Thanks for your reply. If I understand correctly, the link talks about sampling on a patch, which I can do by writing something like
surfaces
(
airfoil
{
type patch;
patches (airfoil);;
// interpolate false;
}
);

inside sampleDict.
But I need to sample data along a cellSet (not patch)
By the cellSet do you mean to say a line/surface inside the fluid domain?

If suppose you want to plot along a line follow the below shown procedure.

Make your sampleDict like the one shown below:

interpolationScheme cellPoint;

setFormat raw;

surfaceFormat raw;

sets
(
name
{
type uniform;
axis y;
start (0 0 0.0005);
end (0 0.425 0.0005);
nPoints 800; }
);
fields (U p);

In a similar fashion you can progress for the surface with some changes like surfaces instead of sets, type patch..etc.
Tushar@cfd is offline   Reply With Quote

Old   May 23, 2013, 01:09
Default
  #5
Senior Member
 
Tarak
Join Date: Aug 2010
Location: State College, PA
Posts: 105
Rep Power: 6
Tarak is on a distinguished road
Hii Tushar,

Thanks for your reply. By cellSet I meant the list of cells which we can create based on multiple criteria (See section 5.1 of http://openfoamwiki.net/index.php/TopoSet).
Tarak is offline   Reply With Quote

Old   May 23, 2013, 01:40
Default
  #6
Senior Member
 
Tushar Chourushi
Join Date: Jul 2009
Location: IIT-Indore, India
Posts: 318
Blog Entries: 1
Rep Power: 9
Tushar@cfd is on a distinguished road
Quote:
Originally Posted by Tarak View Post
Hii Tushar,

Thanks for your reply. By cellSet I meant the list of cells which we can create based on multiple criteria (See section 5.1 of http://openfoamwiki.net/index.php/TopoSet).

Hi Tarak,

Thanks for explaining..
Try looking at the already available topoSetDict with the OF. $tutorial/incompressible/pimpleFoam/TjunctionFan/system/ .
May be bad co-ordinates could be sometimes a problem with the effecting running of utility.
Tushar@cfd is offline   Reply With Quote

Old   July 9, 2013, 16:18
Default
  #7
Senior Member
 
Francois Beaubert
Join Date: Mar 2009
Location: Lille, France
Posts: 133
Rep Power: 8
francois is on a distinguished road
Tarak, Did you find a solution to sample your data on a cellSet/faceSet ? Thanks Regards François
francois is offline   Reply With Quote

Old   July 9, 2013, 17:44
Default sampling of cellSet
  #8
New Member
 
Ashvin Chaudhari
Join Date: Aug 2011
Location: Finland
Posts: 20
Rep Power: 5
ashvinc9 is on a distinguished road
Hi,
One can first run the topoSetDict and try to make a cellSet. Below is the code to keep in topoSetDict in order to make box typed cellSet.

Code:
actions
(
    {
        name    c0; // name of the cellSet
        type    cellSet;
        action  new;
        source  boxToCell; //box type of cellSet
        sourceInfo
        {
            box (1 0.5 0) (3 1.5 0.5); //size or dimensions of the Box, you may change this 
        }
    }
);
After creating the cellSet, one can run following command to produce a VTK file of that cellSet, e.g.,

Code:
foamToVTK -cellSet c0 -latestTime
Hope this helps !

- Ashvin
ashvinc9 is offline   Reply With Quote

Old   July 11, 2013, 07:34
Default
  #9
Senior Member
 
Francois Beaubert
Join Date: Mar 2009
Location: Lille, France
Posts: 133
Rep Power: 8
francois is on a distinguished road
Thanks Ashvin for your tip.

Indeed it's possible to obtain VTK file from cellSet with foamToVTK but i'm looking for something different: I'll like to obtain a raw file of the sampled cellSet which I can load into python.

Is it possible to do it without writing a parser to convert the VTK file to raw data ?

Have a nice day
Regards

Francois
francois is offline   Reply With Quote

Old   July 14, 2013, 21:34
Default
  #10
Senior Member
 
Tarak
Join Date: Aug 2010
Location: State College, PA
Posts: 105
Rep Power: 6
Tarak is on a distinguished road
Hii Francois,

No, I couldn't find a solution. Please let me know if you have one.
Tarak is offline   Reply With Quote

Old   August 21, 2013, 10:08
Default Simple sampling using cellSet
  #11
New Member
 
Johan Roenby
Join Date: May 2011
Location: Denmark
Posts: 9
Rep Power: 6
roenby is on a distinguished road
Dear Tarak

Maybe I misunderstood something, but I think what you are trying to do is straight forward:

The cellSet you create e.g. with topoSet will be in a file, say, <case>/constant/polyMesh/sets/c0. This contains a list of cell indices, e.g. for all cells having centres within a box defined in your topoSetDict.

If the field that you want to sample is a volVectorField (such as the velocity field U) or volScalarField (such as the pressure field p) then the file <case>/<time>/<volFieldName> contains a list with N vectors or scalars, where N is the number of cells in your domain. The n'th value in the list is the field value at the n'th cell.

So if your cellSet contains the values (0 4 10) your sampling is simply a matter of picking out elements 0, 4 and 10 in the files <case>/<time>/<volFieldName>.

Of course for this to be usefull you also need the cell centres for these cells. This can be calculated with the utility writeCellCentres (type writeCellCentres -help in the terminal).

Hope this helped.

Johan

Ps: If the field you want to sample is a surfaceScalarField or a surfaceVectorField the values in the time directories are the field values on the cell faces. Then I suppose you can do exactly the same using a faceSet. The n'th value in your <case>/<time>/<surfaceField> corresponds to the n'th face in your <case>/constant/polyMesh/faces file (the numbers in this faces file refers to the elements in the constant/polyMesh/points file). Unfortunetaly, I don't think there is any writeFaceCentres utility currently so you would have to calculate the face centres yourself. For this you can extract the points for the faces in your faceSet using writeMeshObject -faceSet <faceSetFileName>.
roenby 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
sampling interval and sampling frequency beauty FLUENT 3 November 6, 2014 11:09
topoSet geometry than more than 1 box Grimoli OpenFOAM Pre-Processing 1 April 14, 2013 08:54
Creating a cellSet fom yPlus buffi OpenFOAM Post-Processing 0 December 17, 2012 05:19
cellSet - how to define a hollow cylinder RugbyGandalf OpenFOAM Meshing & Mesh Conversion 2 December 17, 2012 03:34
sampling Data during calculation JimKnopf OpenFOAM Post-Processing 0 May 27, 2011 04:10


All times are GMT -4. The time now is 20:08.