 |
 |
 |
|
|
|
[Sponsors] | ||||
June 9, 2022, 08:40
|
reconstructPar problem
|
#1 |
|
New Member
Afu Lp
Join Date: Mar 2022
Posts: 4
Rep Power: 4  |
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
Any help will be welcome. Best, Dong Yan |
|
|
|
|
|
June 10, 2022, 00:09
|
|
#2 |
|
New Member
Afu Lp
Join Date: Mar 2022
Posts: 4
Rep Power: 4 |
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? |
|
|
|
|
|
|
June 10, 2022, 00:18
|
|
#3 |
|
New Member
Afu Lp
Join Date: Mar 2022
Posts: 4
Rep Power: 4 |
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;
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
|
|
|
|
|
|
July 1, 2022, 00:09
|
|
#4 |
|
New Member
Afu Lp
Join Date: Mar 2022
Posts: 4
Rep Power: 4 |
It seems only occur in ESI version. This error can be repeated in v2112 and v2206, but not in OpenFOAM8.
|
|
|
|
|
|
 |
| Tags |
| reconstructpar |
|
|
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 03:31 |
| BuoyantBoussinesqSimpleFoam_Facing problem | Mondal131211 | OpenFOAM Running, Solving & CFD | 1 | April 10, 2019 20:41 |
| Mesh& steptime independant: conduction-convection problem | Fati1 | Main CFD Forum | 1 | October 28, 2018 14:52 |
| natural convection problem for a CHT problem | Se-Hee | CFX | 2 | June 10, 2007 07:29 |
| Adiabatic and Rotating wall (Convection problem) | ParodDav | CFX | 5 | April 29, 2007 20:13 |
Linear Mode
