View Issue Details

IDProjectCategoryView StatusLast Update
0003119OpenFOAMBugpublic2018-12-05 00:19
ReporterweyburneAssigned Towyldckat 
PrioritynormalSeveritymajorReproducibilityalways
Status closedResolutionno change required 
PlatformblueCFD-Core-2017OSwin10OS Version
Product Version5.x 
Fixed in Version 
Summary0003119: paraFoam reporting an incorrect pressure and velocity value at the multi-grading boundary
DescriptionThe problem occurs with paraFoam (paraView 5.4.1, Win 10) when requesting pressure and velocity slices through a mult-graded region. In the slices, the one data point at the the boundary between the multi-grading areas is slightly incorrect. It maybe an interpolation error, not sure.
Steps To ReproduceThe problem can be seen using a simple 2-D channel, one block geometry simulation. BlockMesh is used with the multi-grading option used in both the flow direction, x-plane, and normal to the channel walls, y-plane. Air flow at RT is used with the channel length and the air velocity chosen to keep the flow laminar (Rex=5e5 at exit). SimpleFoam with the laminar option is run until P and U residuals are less than 2x10^7. paraFoam is then run and center line slices in both flow direction and wall normal are obtained. Attached files include all relevant OpenFoam files as well as two jpg graphs of the pressure illustrating the problem, In simple pressure or velocity plots, the problem is difficult to see unless the boundary region is blown up. It is very noticeable if one takes the derivative of the pressure or velocity.
TagsNo tags attached.

Activities

weyburne

2018-12-04 16:26

reporter  

U (1,372 bytes)
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.4.1-dev                             |
|   \\  /    A nd           | Web:      http://www.openfoam.org               |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/

FoamFile
{
    version 2.0;
    format ascii;

    root "/home/mjaworsk/OpenFOAM/mjaworsk-1.4.1-dev/run/conjugateTray2D";
    case "lithium-5mm";
    instance "0";
    local "";

    class volVectorField;
    object U;
}

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

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0.9375 0 0);

boundaryField
{
    
    inlet
    {
        type            fixedValue;
	      value           uniform (0.9375 0 0);
    }
    outlet
    {
        type            zeroGradient;
    }
    plate
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    

    frontAndBack
    {
        type            empty;
    }
}


// ************************************************************************* //
U (1,372 bytes)
fvSolution (2,137 bytes)
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  5                                     |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    p
    {
        solver          GAMG;
        tolerance       1e-010;
        relTol          0.1;
        smoother        GaussSeidel;
        nPreSweeps      1;
        nPostSweeps     3;
        nFinestSweeps   3;
        scaleCorrection true;
        directSolveCoarsest off;
        catchAgglomeration   on;
        nCellsInCoarstLevel  30;
        agglomerator     faceAreaPair;
        mergeLevels     1;
    }
    pFinal
    {
        $p;
        relTol          0;
    }
    
    "(U|k|epsilon|omega|f|v2)"
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-08;
        relTol          0.1;
    }
}

SIMPLE
{
    nNonOrthogonalCorrectors 0;
    consistent      yes;

    residualControl
    {
        p               1e-10;
        U               1e-8;
        "(k|epsilon|omega|f|v2)" 1e-8;
    }
}

relaxationFactors
{
    equations
    {
        P               0.8
        U               0.95; // 0.9 is more stable but 0.95 more convergent
        ".*"            0.8; // 0.9 is more stable but 0.95 more convergent
    }
}
PISO
{
    nCorrectors     2;
    nNonOrthogonalCorrectors 0;
    pRefCell        0;
    pRefValue       0;
}

potentialFlow
{
    nNonOrthogonalCorrectors 10;
}

// ************************************************************************* //
fvSolution (2,137 bytes)
turbulenceProperties (1,093 bytes)
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  5                                     |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      turbulenceProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

simulationType laminar;

RAS
{
    // Tested with kEpsilon, realizableKE, kOmega, kOmegaSST, v2f,
    // ShihQuadraticKE, LienCubicKE.
    RASModel        kEpsilon;

    turbulence      off;

    printCoeffs     off;
}


// ************************************************************************* //
turbulenceProperties (1,093 bytes)
transportProperties (1,035 bytes)
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.3                                   |
|   \\  /    A nd           | Web:      http://www.openfoam.org               |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/

FoamFile
{
    version         2.0;
    format          ascii;

    root            "";
    case            "";
    instance        "";
    local           "";

    class           dictionary;
    object          transportProperties;
}

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

transportModel  Newtonian;

nu              nu [0 2 -1 0 0 0 0] 1.5e-5;




// ************************************************************************* //
transportProperties (1,035 bytes)
fvSchemes (1,555 bytes)
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  5                                     |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default         steadyState;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default         none;
    div(phi,U)      bounded Gauss linearUpwind grad(U);
    div(phi,k)      bounded Gauss limitedLinear 1;
    div(phi,epsilon) bounded Gauss limitedLinear 1;
    div(phi,omega)  bounded Gauss limitedLinear 1;
    div(phi,v2)     bounded Gauss limitedLinear 1;
    div((nuEff*dev2(T(grad(U))))) Gauss linear;
    div(nonlinearStress) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

wallDist
{
    method meshWave;
}


// ************************************************************************* //
fvSchemes (1,555 bytes)
controlDict (1,359 bytes)
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.3                                   |
|   \\  /    A nd           | Web:      http://www.openfoam.org               |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/

FoamFile
{
    version         2.0;
    format          ascii;

    root            "";
    case            "";
    instance        "";
    local           "";

    class           dictionary;
    object          controlDict;
}

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

application     simpleFoam;

startFrom       startTime;

startTime       0;

stopAt          endTime;

endTime         12.0;

deltaT          0.002;

writeControl    runTime;

writeInterval   3.0;

purgeWrite      0;

writeFormat     ascii;

writePrecision  15;

writeCompression uncompressed;

timeFormat      general;

timePrecision   6;

runTimeModifiable yes;

functions
{
     #includeFunc residuals
}

// ************************************************************************* //
controlDict (1,359 bytes)
blockMeshDict (1,746 bytes)
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.3                                   |
|   \\  /    A nd           | Web:      http://www.openfoam.org               |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/

FoamFile
{
    version         2.0;
    format          ascii;

    root            "";
    case            "";
    instance        "";
    local           "";

    class           dictionary;
    object          blockMeshDict;
}

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


convertToMeters 0.001;

vertices        
(
    // Block 1
    (0 0 0)   //0
    (8000 0 0)     //1
    (8000 1000 0)   //2
    (0 1000 0) //3
    (0 0 1)   //4
    (8000 0 1)     //5
    (8000 1000 1)   //6
    (0 1000 1) //7
    
  
);

blocks          
(
    // Block 1
    hex (0 1 2 3 4 5 6 7) (1000 600 1)    
    simpleGrading
        (
            ((0.05 0.3 10) (0.9 0.4 1) (0.05 0.3 0.1))
            ((0.05 0.3 10) (0.9 0.4 1) (0.05 0.3 0.1))
            1
         )
    
  
);

edges           
(
);

patches         
(

    wall inlet
    (
        (0 4 7 3)
    )
    wall outlet
    (
        (1 5 6 2)
 
    )
    wall plate
    (    
	      (7 6 2 3)
	      (0 1 5 4)
        
    )
    
    empty frontAndBack 
    (
        (0 1 2 3)
        (4 5 6 7)
	    
 
    )
);




// ************************************************************************* //
blockMeshDict (1,746 bytes)
p (1,351 bytes)
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.3                                   |
|   \\  /    A nd           | Web:      http://www.openfoam.org               |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/

FoamFile
{
    version         2.0;
    format          ascii;

    root "/home/mjaworsk/OpenFOAM/mjaworsk-1.4.1-dev/run/conjugateTray2D";
    case "lithium-5mm";
    instance "0";
    local "";

    class           volScalarField;
    object          p;
}

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


dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{

    inlet    
    {
        type            zeroGradient;
    }

    outlet
    {
        type            fixedValue;
	      value		        uniform 0;
    }
    plate
    {
        type            zeroGradient;
    }
    
    frontAndBack    
    {
        type            empty;
    }
}


// ************************************************************************* //
p (1,351 bytes)
GraphPy4.jpg (149,405 bytes)
GraphPy4.jpg (149,405 bytes)

wyldckat

2018-12-05 00:19

updater   ~0010204

The blueCFD-Core project has its own issue tracker: https://github.com/blueCFD/Core/issues

Please do not use the OpenFOAM Foundation's issue tracker to report problems that occur with independent projects, unless of course the same issue was reproduced with a version that is provided by the OpenFOAM Foundation.

Issue History

Date Modified Username Field Change
2018-12-04 16:26 weyburne New Issue
2018-12-04 16:26 weyburne File Added: U
2018-12-04 16:26 weyburne File Added: fvSolution
2018-12-04 16:26 weyburne File Added: turbulenceProperties
2018-12-04 16:26 weyburne File Added: transportProperties
2018-12-04 16:26 weyburne File Added: fvSchemes
2018-12-04 16:26 weyburne File Added: controlDict
2018-12-04 16:26 weyburne File Added: blockMeshDict
2018-12-04 16:26 weyburne File Added: p
2018-12-04 16:26 weyburne File Added: GraphPy4.jpg
2018-12-05 00:19 wyldckat Note Added: 0010204
2018-12-05 00:19 wyldckat Assigned To => wyldckat
2018-12-05 00:19 wyldckat Status new => closed
2018-12-05 00:19 wyldckat Resolution open => no change required