[Xulia]

Dear Foamers:

I have been using OF-2.1.x for the past two months, forgive me if my problem is very easy/obvious to fix. I am running simulations on open channel so my BCs are cyclic in the spanwise direction. I am generating turbulence (using LES) that I will save as initial conditions for the next simulation: turbine performance. I want to save a plane from the middle of the channel (normal to streamwise direction) containing data of U, p, U_0, etc during different times, that I will use as inlet for the turbine simulation. I found that cuttingPlane can do that and I imitated the case incompressible/simpleFoam/motorbike.

This is my cuttingPlane function:
cuttingPlane
{
type surfaces;
functionObjectLibs ("libsampling.so");
outputControl outputTime;

surfaceFormat vtk;
fields ( p U U_0 nuSgs);

interpolationScheme cellPoint;

surfaces
(
yNormal
{
type cuttingPlane;
planeType pointAndNormal;
pointAndNormalDict
{
basePoint (2 1 1);
normalVector (0 1 0);
}
interpolate true;
}
);
}


When I run the solver (channelFoam), it gives me this error:

--> FOAM FATAL ERROR:
Unhandled coupledPolyPatch type cyclic

From function isoSurface::collocatedFaces(const coupledPolyPatch&) const
in file sampledSurface/isoSurface/isoSurface.C at line 94.

FOAM aborting

I can't avoid the cyclic BC, since it represents the width of my channel and I need it to be that way... I am confused because I went through the cuttingPlane case (/src/sampling/cuttingPlane) and read all the functions and I don't know when isoSurface.C is called... I can't find that line....
Has anybody succeed in using cuttingPlane involving cyclic BCs?
Is there an alternative way to save a plane during simulations?

Many thanks!!

[vishwakarma]

Originally Posted by Xulia

Dear Foamers:

I have been using OF-2.1.x for the past two months, forgive me if my problem is very easy/obvious to fix. I am running simulations on open channel so my BCs are cyclic in the spanwise direction. I am generating turbulence (using LES) that I will save as initial conditions for the next simulation: turbine performance. I want to save a plane from the middle of the channel (normal to streamwise direction) containing data of U, p, U_0, etc during different times, that I will use as inlet for the turbine simulation. I found that cuttingPlane can do that and I imitated the case incompressible/simpleFoam/motorbike.

This is my cuttingPlane function:
cuttingPlane
{
type surfaces;
functionObjectLibs ("libsampling.so");
outputControl outputTime;

surfaceFormat vtk;
fields ( p U U_0 nuSgs);

interpolationScheme cellPoint;

surfaces
(
yNormal
{
type cuttingPlane;
planeType pointAndNormal;
pointAndNormalDict
{
basePoint (2 1 1);
normalVector (0 1 0);
}
interpolate true;
}
);
}


When I run the solver (channelFoam), it gives me this error:

--> FOAM FATAL ERROR:
Unhandled coupledPolyPatch type cyclic

From function isoSurface::collocatedFaces(const coupledPolyPatch&) const
in file sampledSurface/isoSurface/isoSurface.C at line 94.

FOAM aborting

I can't avoid the cyclic BC, since it represents the width of my channel and I need it to be that way... I am confused because I went through the cuttingPlane case (/src/sampling/cuttingPlane) and read all the functions and I don't know when isoSurface.C is called... I can't find that line....
Has anybody succeed in using cuttingPlane involving cyclic BCs?
Is there an alternative way to save a plane during simulations?

Many thanks!!

Did you get the answer to your question? I am facing the same problem and not able to resolve it.