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/)
-   -   wrong result using "interpolate" on postProcess cuttingPlane (https://www.cfd-online.com/Forums/openfoam-post-processing/220152-wrong-result-using-interpolate-postprocess-cuttingplane.html)

Za-ck August 26, 2019 19:42

wrong result using "interpolate" on postProcess cuttingPlane
 
2 Attachment(s)
Hi,

I would like to plot the alpha field on a cutting plane using the runTimePostProcessing tool. Everything looks right without interpolation, but with activated interpolation the values seems to get mixed up. Same error on an isoSurface. No problems using a plane instead of a cuttingPlane.

Any ideas?

Cheers,
Za-ck

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  v1906                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.com                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
FoamFile
{
        version    1906;
        format      ascii;
        class      dictionary;
        location    "system";
        object      visualization;
}


functions
{
    planes
    {
        type            surfaces;
        libs            ("libsampling.so");

        writeControl    writeTime;
        writeFields    true;

        fields          (alpha.water);

        surfaceFormat  none;

        store          true;

        //cell cellPatchConstrained cellPoint cellPointFace cellPointWallModified pointMVC
        interpolationScheme cellPoint;

       

        surfaces
        {
            xPlane
            {
                type            cuttingPlane; //plane; //
                planeType      pointAndNormal;
                interpolate    true;
//                regularise      false;
                pointAndNormalDict
                {
                    point  (0 0 0);
                    normal  (1 0 0);
                }
            }           
               
        };

    }



    postPro1
    {
        #includeEtc "caseDicts/postProcessing/visualization/runTimePostPro.cfg"

        // Time control etc
        writeControl    writeTime;

        output
        {
            name        image;
            width      1920;
            height      1200;
        }
       
       
       
        camera
        {
            nFrameTotal    1;
            parallelProjection  no;
            position        (-20 0 0);
            focalPoint      (-1 0 0.0);

            up              (0.0 0.0 1.0);
        }

        // Default colours
        colours
        {
            ${..colourScheme.paraview};
        }

        surfaces
        {
            //A cutting plane from sampled surfaces:
            stored1
            {
                type            functionObjectSurface;
                functionObject  planes.xPlane; //planes.freeSurface; //planes.test2; //
                colourMap      coolToWarm;
                representation  surface;
                visible        yes;
                featureEdges    no;
                colourBy        field;
                field          alpha.water;
                range          (0.0 1.0);
                opacity        1;

            }
 

        }
       
        text
        {
           
        }

    }

}


tecmul August 28, 2019 11:08

I'm confused about what the "interpolate" switch does. Doesn't sampling on a plane always need some form of interpolation?


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