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

reconstructPar problem

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 9, 2022, 07:40
Default reconstructPar problem
  #1
New Member
 
Afu Lp
Join Date: Mar 2022
Posts: 4
Rep Power: 3
Dong Yan is on a distinguished road
Hi, foamers.

I performed a simulation of curved pipe in parallel and it looks fine with type of "decomposed case" in paraview. However, I restructure the data by "reconstructPar" and read the data with type of "reconstructed case" in paraview. The flow field looks very strange. The attachment has a picture of my grid, it should not be a problem with the grid. I think it should be that the data is not being properly restructured, but I don't know how to fix it.

There is no error when running reconstructPar.
the log is :
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time



Reconstructing fields
region=region0

Time = 3800

Reconstructing FV fields

    Reconstructing volScalarFields

        epsilon
        k
        nut
        p
        turbulenceProperties:I
        turbulenceProperties:L
        turbulenceProperties:epsilon
        turbulenceProperties:k
        yPlus

    Reconstructing volVectorFields

        U
        wallShearStress

    Reconstructing volSymmTensorFields

        turbulenceProperties:R

    Reconstructing surfaceScalarFields

        phi

Reconstructing point fields

No point fields

No lagrangian fields

No FA fields


End
I'm using v2112.
Any help will be welcome.

Best,
Dong Yan
Attached Images
File Type: png 11.png (133.7 KB, 13 views)
File Type: jpg 12.jpg (35.2 KB, 12 views)
File Type: jpg 13.jpg (125.1 KB, 9 views)
File Type: jpg 14.jpg (202.7 KB, 10 views)
Dong Yan is offline   Reply With Quote

Old   June 9, 2022, 23:09
Default
  #2
New Member
 
Afu Lp
Join Date: Mar 2022
Posts: 4
Rep Power: 3
Dong Yan is on a distinguished road
I have tried both the mesh generated by ICEM and the mesh generated by blockMesh without extrudeMesh. The above error did not occur.

Then I tried to convert the mesh with the above error to “.msh” format using foamMeshToFluent, and then using fluentMeshToFoam to convert it back to the polymesh format, and the above error did not occur.

Therefore, I guess it may be caused by extrudeMesh. Maybe there is a problem with overlapping points when extruding the mesh with extrudeMesh?
Dong Yan is offline   Reply With Quote

Old   June 9, 2022, 23:18
Default
  #3
New Member
 
Afu Lp
Join Date: Mar 2022
Posts: 4
Rep Power: 3
Dong Yan is on a distinguished road
my blockMeshDict and extrudeMeshDicts are below:

Code:
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  8
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

//convertToMeters 0.001;  // minimeter 
scale           0.001;

//Auxiliary settings
rPipe                       20;
dPipe                       40;
lengthUpstream              #calc "3 * $dPipe"; // length of upstream pipe = 3d , 6d for period pipe flow
lengthDownsteam             #calc "3 * $dPipe"; // length of downstream pipe = 3d

ogridFactor                 0.68; // 12/(12+5.5) JFM 2016 
rPipeOgrid                  #calc "$rPipe * $ogridFactor";
rPipeOgrid2                 #calc "$rPipe * $ogridFactor / pow(2,0.5) * 1.18"; // 1.18

/*
curvature                   #calc "2 * $dPipe"; //radius of curvature = 2d
    // so, the central point of bendouelet is located in ( " $lengthUpstream + $curvature " " $curvature " 0)
xCentralPointOfBendOutlet   #calc "$lengthUpstream + $curvature";
yCentralPointOfBendOutlet   #calc "$curvature";
zCentralPointOfBendOutlet   0;

xCentralPointOfOutlet       #calc "$lengthUpstream + $curvature";
yCentralPointOfOutlet       #calc "$lengthDownstream + $curvature";
zCentralPointOfOutlet       0;
*/

xCellsUpsteam               128;  // streamwise, when xcell=1000, the dleta x plus = 43.6 
//xCellsBend                  0;
//xCellDownsteam              0;

//xcells                      100;
ycells                      50;
zcells                      42;

ogridCells                  42;

yPlus1                      0;
//stretch                     0.25;
//Auxiliary settings

vertices
(
    //upsteam pipe
    // block 0
    (0                  #calc "$rPipeOgrid / pow(2,0.5)"                 #calc "-1 * $rPipeOgrid / pow(2,0.5)")               //0
    ($lengthUpstream    #calc "$rPipeOgrid / pow(2,0.5)"                 #calc "-1 * $rPipeOgrid / pow(2,0.5)")               //1
    ($lengthUpstream    #calc "$rPipe / pow(2,0.5)"                      #calc "-1 * $rPipe / pow(2,0.5)")                    //2
    (0                  #calc "$rPipe / pow(2,0.5)"                      #calc "-1 * $rPipe / pow(2,0.5)")                    //3
    (0                  #calc "$rPipeOgrid / pow(2,0.5)"                 #calc "$rPipeOgrid / pow(2,0.5)")                    //4
    ($lengthUpstream    #calc "$rPipeOgrid / pow(2,0.5)"                 #calc "$rPipeOgrid / pow(2,0.5)" )                   //5
    ($lengthUpstream    #calc "$rPipe / pow(2,0.5)"                      #calc "$rPipe / pow(2,0.5)")                         //6
    (0                  #calc "$rPipe / pow(2,0.5)"                      #calc "$rPipe / pow(2,0.5)")                         //7


    //block 1 new vertex
    (0                  #calc "-1 * $rPipeOgrid / pow(2,0.5)"            #calc "$rPipeOgrid / pow(2,0.5)")                    //8
    (0                  #calc "-1 * $rPipe / pow(2,0.5)"                  #calc "$rPipe / pow(2,0.5)")                         //9
    ($lengthUpstream    #calc "-1 * $rPipeOgrid / pow(2,0.5)"            #calc "$rPipeOgrid / pow(2,0.5)")                    //10
    ($lengthUpstream    #calc "-1 * $rPipe / pow(2,0.5)"                  #calc "$rPipe / pow(2,0.5)")                         //11

    //block 2 new vertex
    (0                  #calc "-1 * $rPipeOgrid / pow(2,0.5)"            #calc "-1 * $rPipeOgrid / pow(2,0.5)")               //12
    (0                  #calc "-1 * $rPipe / pow(2,0.5)"                 #calc "-1 * $rPipe / pow(2,0.5)")                    //13
    ($lengthUpstream    #calc "-1 * $rPipeOgrid / pow(2,0.5)"            #calc "-1 * $rPipeOgrid / pow(2,0.5)")               //14
    ($lengthUpstream    #calc "-1 * $rPipe / pow(2,0.5)"                 #calc "-1 * $rPipe / pow(2,0.5)")                    //15

    //bend pipe
    //block 5 new vertex
    
);

blocks
(
    //block 0
    hex (0 1 2 3 4 5 6 7)       ($xCellsUpsteam $ycells $ogridCells) simpleGrading (
                                                                          1                   // x
                                                                            (
                                                                        (0.8 0.45   0.5)        // y 20%, cellNumber 30%, expension ratio=0.25
                                                                        (0.2 0.55 0.071)        //   
                                                                                )
                                                                                        1     // z
                                                                                        )
    //block 1
    hex (4 5 6 7 8 10 11 9)     ($xCellsUpsteam $ycells $ogridCells) simpleGrading (
                                                                          1                   // x
                                                                            (
                                                                        (0.8 0.45   0.5)        // y 20%, cellNumber 30%, expension ratio=0.25
                                                                        (0.2 0.55 0.071)        //   
                                                                                )
                                                                                        1     // z
                                                                                        )
    //block 2
    hex (8 10 11 9 12 14 15 13) ($xCellsUpsteam $ycells $ogridCells) simpleGrading (
                                                                          1                   // x
                                                                            (
                                                                        (0.8 0.45   0.5)        // y 20%, cellNumber 30%, expension ratio=0.25
                                                                        (0.2 0.55 0.071)        //   
                                                                                )
                                                                                        1     // z
                                                                                        )
    //block 3
    hex (12 14 15 13 0 1 2 3)   ($xCellsUpsteam $ycells $ogridCells) simpleGrading (
                                                                          1                   // x
                                                                            (
                                                                        (0.8 0.45   0.5)        // y 20%, cellNumber 30%, expension ratio=0.25
                                                                        (0.2 0.55 0.071)        //   
                                                                                )
                                                                                        1     // z
                                                                                        )
    //block 4
    hex (0 1 5 4 12 14 10 8)    ($xCellsUpsteam $ogridCells $ogridCells) simpleGrading (1 1 1)
    //block 5

);

edges           
(
    //block 0 arc
    arc     3   7   (0                  $rPipe      0)
    arc     2   6   ($lengthUpstream    $rPipe      0)

    //block 1
    arc     7   9   (0                  0           $rPipe)
    arc     6   11  ($lengthUpstream    0           $rPipe)

    //block 2
    arc     9   13  (0                  #calc "-1 * $rPipe"     0)
    arc     11  15  ($lengthUpstream    #calc "-1 * $rPipe"     0)

    //block 3
    arc     13  3   (0                  0                       #calc "-1 * $rPipe")
    arc     15  2   ($lengthUpstream    0                       #calc "-1 * $rPipe")

    //block 4
    arc     0   4   (0                  $rPipeOgrid2                0)
    arc     4   8   (0                  0                           $rPipeOgrid2)
    arc     8   12  (0                  #calc "-1 * $rPipeOgrid2"   0)
    arc     12  0   (0                  0                           #calc "-1 * $rPipeOgrid2")

    arc     1   5   ($lengthUpstream    $rPipeOgrid2                0)
    arc     5   10  ($lengthUpstream    0                           $rPipeOgrid2)
    arc     10  14  ($lengthUpstream    #calc "-1 * $rPipeOgrid2"   0)
    arc     14  1   ($lengthUpstream    0                           #calc "-1 * $rPipeOgrid2")
);

boundary
(
    inlet 
    {
        //type            cyclic;
        //neighbourPatch  outlet;
        type            patch;
        faces
        (
            (3 7 4 0)
            (4 7 9 8)
            (8 9 13 12)
            (12 13 3 0)
            (0 4 8 12)
        );
    }

    outlet 
    {
        //type            cyclic;
        //neighbourPatch  inlet;
        type            patch;
        faces
        (
            (2 6 5 1)
            (5 6 11 10)
            (10 11 15 14)
            (14 15 2 1)
            (1 5 10 14)
        );
    }

    walls 
    {
        type            wall;
        faces
        (
            (3 7 6 2)
            (7 9 11 6)
            (9 13 15 11)
            (13 3 2 15)
        );
    }
);

mergePatchPairs
(

);

// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2112                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      extrudeMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //



constructFrom   mesh;

sourceCase      ".";
sourcePatches   (outlet);

flipNormals     false;
exposedPatchName    outlet;

//extrudeModel    linearDirection;//polyline;
extrudeModel    sector;//polyline;

nLayers         128;

expansionRatio  1.0;

/*
polylineCoeffs
{
    vertices
    (
        (0.24 0 0)
        (0.32 0.08 0)
    );

    edges
    (
        arc 0 1  (0.24 0.08 0)
    );
}
*/

sectorCoeffs    //<- Also used for wedge
{
    point       (0.12 0.08 0);
    axis        (0 0 1);
    angle       90;  // For nLayers=1 assume symmetry so angle/2 on each side
}

mergeFaces false;

mergeTol 0;
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2112                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      extrudeMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //



constructFrom   mesh;

sourceCase      ".";
sourcePatches   (outlet);

flipNormals     false;
exposedPatchName    outlet;

//extrudeModel    linearDirection;//polyline;
//extrudeModel    sector;//polyline;
extrudeModel    linearDirection;

nLayers         62;

expansionRatio  1.0;

/*
polylineCoeffs
{
    vertices
    (
        (0.24 0 0)
        (0.32 0.08 0)
    );

    edges
    (
        arc 0 1  (0.24 0.08 0)
    );
}
*/

/*
sectorCoeffs    //<- Also used for wedge
{
    point       (0.24 0.08 0);
    axis        (0 0 1);
    angle       90;  // For nLayers=1 assume symmetry so angle/2 on each side
}
*/

linearDirectionCoeffs
{
    direction       (0 1 0);
    thickness       0.06;
}

mergeFaces false;

mergeTol 0;
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Dong Yan is offline   Reply With Quote

Old   June 30, 2022, 23:09
Default
  #4
New Member
 
Afu Lp
Join Date: Mar 2022
Posts: 4
Rep Power: 3
Dong Yan is on a distinguished road
It seems only occur in ESI version. This error can be repeated in v2112 and v2206, but not in OpenFOAM8.
Dong Yan is offline   Reply With Quote

Reply

Tags
reconstructpar

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Is COMSOL Multi Physics is suitable to solve complex flow problem? steve lee COMSOL 8 January 5, 2023 02:31
BuoyantBoussinesqSimpleFoam_Facing problem Mondal131211 OpenFOAM Running, Solving & CFD 1 April 10, 2019 19:41
Mesh& steptime independant: conduction-convection problem Fati1 Main CFD Forum 1 October 28, 2018 13:52
natural convection problem for a CHT problem Se-Hee CFX 2 June 10, 2007 06:29
Adiabatic and Rotating wall (Convection problem) ParodDav CFX 5 April 29, 2007 19:13


All times are GMT -4. The time now is 22:47.