CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Bugs

Discrepancy ensight slice and full data for parallel decomposed case

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   August 27, 2013, 04:47
Default Discrepancy ensight slice and full data for parallel decomposed case
  #1
New Member
 
Jasper
Join Date: May 2013
Location: Netherlands
Posts: 1
Rep Power: 0
Japsor is on a distinguished road
Dear all,

Yesterday I encountered some unwanted behaviour when analyizing only a slice of data, generated by an addition to system/controldict. This occures for a decomposed case run in parallel.

When comparing results from the full dataset, by opening the .foam file with paraview4.1 the results look fine but when loading an ensight .case file used to generate only a slice, the slice shows sudden high velocity regions where it is both physically unrealistic and also does not match the full data. I reproduced this problem with the simple tut/icofoam/cavity-case.

I work with OpenFoam-2.2.0, viewed the full .foam file with paraview4.1 and generated an additional .foam file for each processor (in this case 3).
This was done to check if the boundaries of the processor-domains match the edges of the anomalies. It turns out they do.
I attached some figures to illustrate the problem. Plotted are the seperated domains, using a slice at the centre, just like the generated .case slice-file. The controldict and decomposePardict are below this message.

This problem was not found in OpenFoam-2.1.0. Has anyone else encountered a problem like this and know how to solvethis? As I can work with OpenFoam2.1 as well this means that it problably will not pose much of a problem for me, but I would at least like to report the bug and hopefully the problem can be solved. I hope you can help. Lastly, I hope this whole post is in an acceptable layout, as it is only my first post.

gr. Jasper

Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
libs ( 
       "libOpenFOAM.so"
//     "libdynamicSmagorinskyModel.so"
       "libmyEnsightSampling.so"
      );

application     PimpleFoam;
 
startFrom       startTime;

startTime       0;

stopAt          endTime;

endTime         20; 

deltaT          0.001;

writeControl    timeStep;

writeInterval   500;

purgeWrite      0;

writeFormat     binary;

writePrecision  6;

writeCompression on;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;


// ************************************************************************* /

functions
{
    sliceCentre
    {
        type            surfaces;
        functionObjectLibs
        (   
            "libsampling.so" 
        );  
        outputControl   outputTime;
        outputInterval  1; 
//      surfaceFormat   reducedEnsight;
    surfaceFormat   ensight;   
    fields
        (   
            U
            p
        nuSgsi
//        Psi
//        J
        );  
        surfaces
        (   
            sliceCentre
            {   
                type            plane;
                basePoint       (0.0 0.4 0.0);
                normalVector    (0.0 0.0 1.0);
                interpolate     true;
            }   
        );  
        interpolationScheme cell;
    }

    timeAverage
    {
        type    fieldAverage;
        enable  on;
        functionObjectLibs ("libfieldFunctionObjects.so");
        cleanRestart    false;
        resetOnOutput   false;
        outputControl   outputTime;
        outputInterval  1;
        timeStart       20; 
        timeEnd         $endTime;
        fields
        (
            U
            {
                mean        on;
                prime2Mean  on;
                base        time;
            }
        );
    }
}
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    note        "mesh decomposition control dictionary";
    object      decomposeParDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

numberOfSubdomains  3;

method          scotch;

scotchCoeffs
{
    //processorWeights
    //(
    //    1
    //    1
    //    1
    //    1
    //);
    //writeGraph  true;
    //strategy "b";
}


// ************************************************************************* //
Attached Images
File Type: png slice.png (71.6 KB, 2 views)
File Type: png processor0.png (17.3 KB, 2 views)
File Type: png processor1.png (10.8 KB, 2 views)
File Type: png processor2.png (21.2 KB, 1 views)

Last edited by Japsor; August 27, 2013 at 07:15. Reason: clarification
Japsor is offline   Reply With Quote

Reply

Tags
decomposed, ensight, parallel

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On



All times are GMT -4. The time now is 02:17.