[PanPeter]

Hi,

I am sampling values on a half-plane using swak. The example beneath works fine but I would like to know if there is a possibility to also export the defined surface as a stl.

Or would it be possible to pre-define a surface that I can supply my swak expression with? Since my mesh is moving, it should not be a faceSet or faceZone, as these are moving with the mesh.

Thanks

name
    {
    	type swakExpression;
    	valueType surface;
    	surfaceName name2;
    	surface {
    	    type plane;
    	    surfaceType searchablePlate;
    	    span         (100 -100 0);
    	    origin       (0 1 0);
    	    normalVector    (0 0 1);
    	    basePoint    (0 1 0);
    	    // interpolate false;
    	}
    	verbose true;
    	expression "mag(U)";
    	ignore_unimplemented_simpleFunctionObject::movePoints true;
    	accumulations (
    	    min
    	    max
    	    );
    }

[PanPeter]

I can't define a half plane as described above:

When changing the expression to

mag(Sf())

the surface are of my plane should be returned. It shows the surface area of the whole cut plane, not of my half-plane only.

Has anybody defined a searchablePlate inside a swakExpression and could give me a hint how this could be solved?

[PanPeter]

Hi,
I've tried to find solutions to define the sampling over a half plane and couldn't find anything but the definition of external stl files.

So I thought of this work-around: A plane is defined, instead of a half plane, and the swak expression only takes values of cells, whose dependent coordinate is smaller than some threshold (y0 in this case).

underScreenSurface{
	y0 // define it here
    	type swakExpression;
    	valueType surface;
	surfaceName normalZ;
	surface{
	    type plane;
	    basePoint		(0 0 0);
	    normalVector    	(0 0 1);
	}
    	verbose true;
    	expression "(pos().y <= $y0) ? mag(Sf()) : 0";
    	accumulations (
	    sum
    	    );
    }

[gschaider]

Originally Posted by PanPeter

Hi,
I've tried to find solutions to define the sampling over a half plane and couldn't find anything but the definition of external stl files.

So I thought of this work-around: A plane is defined, instead of a half plane, and the swak expression only takes values of cells, whose dependent coordinate is smaller than some threshold (y0 in this case).

underScreenSurface{
    y0 // define it here
        type swakExpression;
        valueType surface;
    surfaceName normalZ;
    surface{
        type plane;
        basePoint        (0 0 0);
        normalVector        (0 0 1);
    }
        verbose true;
        expression "(pos().y <= $y0) ? mag(Sf()) : 0";
        accumulations (
        sum
            );
    }

Try adding

        writeSurfaceOnConstruction true;
        autoWriteSurface true;
        surfaceFormat vtk;

You should get a directory surfacesRepository with the vtk-files in each time-step that is written (only the geometry. No values)