[Thomas Baumann]

Hallo all,

I built an axissymmetric mesh with the blockmesh utility. It is a pipe with an inlet and outlet, on top a wall, down is the axis and left and right are wedges.
The radius is about one meter.

Running with icofoam works, but looking at the velocity at the wall there are nonzero values.

I plotted the velocity-values in paraview.
I also used the sampledict with different interpolation schemes (cell, cellPoint, cellPointFace).
Plotting the wall surface shows again zero values.

Has anybody an idea?

Thanks, Thomas

[sega]

Originally Posted by Thomas Baumann

I plotted the velocity-values in paraview.

Looks like ParaView didn't plot the value directly at the wall.
Can you show us the result with with sample?

[Thomas Baumann]

Hallo,

here are the velocity-profils made by the sampledict with the three interpolation schemes: cell, cellPoint and cellPointFace.

The set was made by a uniform line with 40 values from r =0.9R until r = R.

[sega]

Instead of using a line with 40 points you can try using a midPointAndFace type, so you will get values at points corresponding to your discretization.
In this case you can even omit the nPoints-entry!

Something like this:

interpolationScheme cellPointFace;

setFormat     raw;

surfaceFormat    foamFile;

sets      
(
    line_horizontal 
    {
        type            midPointAndFace;

        axis            z;
        start           (0.8838 0.8838 0);
        end             (0.8838 0.8838 5);
    }
);

[Thomas Baumann]

I have used midPointAndFace, too, but nothing changed. There is still nonzero velocity at the wall.

[sega]

I have played with the case but could not find a solution to this problem, at first glance.

Anyone got further ideas?

[Thomas Baumann]

I think the interpolation sheme is the problem. I changed in the sampledict the interpolationScheme to cell (constant over cells) and the type to midPoint (one Point per cell) and the results are much better. Plotting the velocity graphs with gnuplot show good results.