CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Post-Processing (https://www.cfd-online.com/Forums/openfoam-post-processing/)
-   -   Problem with circle sampling type, of Sample utility (https://www.cfd-online.com/Forums/openfoam-post-processing/118446-problem-circle-sampling-type-sample-utility.html)

nik May 28, 2013 04:17

Problem with circle sampling type, of Sample utility
 
Hi,

I am trying to use the sample utility to extract data along a circular line. To do so, I am using the "circle" sampling type. The problem that I have is that in the first column of the output file, where theta should be, I always get zeros. This "circle" type is not listed in the users guide, so the following dictionary that I am using is a result of trial and error, and by looking at the source sampleSet.H. So I am not sure if I am using all the right parameters.

This is the sampleDict that I am using:

Code:

interpolationScheme cell;
setFormat          raw;

sets
(
    alongCircWall
    {
        type        circle;
        origin      (0 0 0);    // Origin (x, y, z) in global cartesian co-ordinates
        circleAxis  (1 0 0);    // Axis of the circle
        startPoint  (0 0 0.6);  // Defines start point on circle (x, y, z) in global cartesian co-ordinates
        dTheta      2;          // Sampling interval in degrees about the origin
        axis        x;          // This I don't know why is needed, but needs to be defined as well?
    }
);

fields              ( T );

By the way, I am using OpenFOAM 2.2.x. Any help would be appreciated.

Regards,
Nicolas

nik May 30, 2013 09:14

Figured it out
 
I have figured it out after all.

The problem was with the "axis x;" parameter. This parameter defines which coordinates are to be written in the output file for each point. You can select, x, y, z, for a single coordinate, or xyz for all coordinates. Since writing the angle theta is not an option, I used xyz to get the coordinates of each point along the circle, and some awk commands to calculate the angle.

It is also worth notting that this sampling type writes out the data along the circle in a weird order. So before plotting the values of a field versus the angle theta, you need to take a close look in the order in which the points are written.

Nicolas

daniel54431 February 23, 2015 12:01

Circle sampling on boundary patch
 
Hello all,

i know this thread is a bit older but my question fits in this thread very well.
I try to use the circle sample to get circular distributed values from
a boundary patch (e.g. wallheatflux) from a pipe.
The problem is when I set dTheta=1 for example I do not receive
360 values from 360 points, I just get around 50 values from 50 points.
For the remaining points I get the message:

Code:

--> FOAM Warning :
    From function void circleSet::calcSamples(DynamicList<point>&, DynamicList<label>&, DynamicList<label>&, DynamicList<label>&, DynamicList<scalar>&) const
    in file sampledSet/circle/circleSet.C at line 128
    Unable to find cell at point id 45 at location (0.0248137 0.00304673 0.5)

The points are definitely in the domain.
So how could it be, that OF finds no cells for interpolation
at most points?


Thanks for help.

Regards,

Daniel

Cenk20 May 5, 2016 09:01

Quote:

Originally Posted by daniel54431 (Post 533060)
Hello all,

i know this thread is a bit older but my question fits in this thread very well.
I try to use the circle sample to get circular distributed values from
a boundary patch (e.g. wallheatflux) from a pipe.
The problem is when I set dTheta=1 for example I do not receive
360 values from 360 points, I just get around 50 values from 50 points.
For the remaining points I get the message:

Code:

--> FOAM Warning :
    From function void circleSet::calcSamples(DynamicList<point>&, DynamicList<label>&, DynamicList<label>&, DynamicList<label>&, DynamicList<scalar>&) const
    in file sampledSet/circle/circleSet.C at line 128
    Unable to find cell at point id 45 at location (0.0248137 0.00304673 0.5)

The points are definitely in the domain.
So how could it be, that OF finds no cells for interpolation
at most points?


Thanks for help.

Regards,

Daniel

Hi, do you find a solution for your problem. I have the same problem now...

Regards,
Cenk

niran May 29, 2017 06:27

circularDistribution
 
surfaceElevation1
{
type waveGauge;

pointDistribution circularDistribution;
N 36;

centre (-0.5657 -0.2 0);
add (0 0.4 0);
radius 0.23;
axis y;
}

it is working in waves2Foam
depth in Y direction
N number of WaveGauge


All times are GMT -4. The time now is 21:40.