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

flow over cylinder in openFoam

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

Like Tree1Likes
  • 1 Post By saeed jamshidi

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 10, 2023, 14:33
Default flow over cylinder in openFoam
  #1
Senior Member
 
Saeed Jamshidi
Join Date: Aug 2019
Posts: 214
Rep Power: 7
saeed jamshidi is on a distinguished road
Hi,
I,m going to simulate flow over cylinder in openfoam. Here is my procedures:
First i create background mesh.
second i create cylinder and fluid stl file and add it to triSurface folder.
finally i apply snappyHexMesh, but there is an error when i apply pimpleFoam:(!
Also, where i should apply cylinder boundary condition as a wall?
i have attached all of my files:
Would you please assist me at this matter:)?
Thank you.

*************************************************

myfoam@DESKTOP-TK3D7CI:~/OpenFOAM-sims/tutorials/incompressible/pimpleFoam/laminar/cylinder2D$ cd
myfoam@DESKTOP-TK3D7CI:~$ cd OpenFOAM-sims/tutorials/incompressible/pimpleFoam/laminar/cylinder2D
myfoam@DESKTOP-TK3D7CI:~/OpenFOAM-sims/tutorials/incompressible/pimpleFoam/laminar/cylinder2D$ blockMesh
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2112 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : _6e1fca0e-20220610 OPENFOAM=2112 patch=220610 version=2112
Arch : "LSB;label=32;scalar=64"
Exec : blockMesh
Date : Aug 10 2023
Time : 21:50:18
Host : DESKTOP-TK3D7CI
PID : 1570
I/O : uncollated
Case : /home/myfoam/OpenFOAM-sims/tutorials/incompressible/pimpleFoam/laminar/cylinder2D
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Creating block mesh from "system/blockMeshDict"
--> FOAM IOWarning :
Found [v1012] 'convertToMeters' entry instead of 'scale' in dictionary "system/blockMeshDict"

This keyword is 132 months old.

Creating block edges
No non-planar block faces defined
Creating topology blocks

Creating topology patches - from boundary section

Creating block mesh topology

Check topology

Basic statistics
Number of internal faces : 0
Number of boundary faces : 6
Number of defined boundary faces : 6
Number of undefined boundary faces : 0
Checking patch -> block consistency

Creating block offsets
Creating merge list (topological search)...
Deleting polyMesh directory "constant/polyMesh"

Creating polyMesh from blockMesh
Creating patches
Creating cells
Creating points with scale (1 1 1)
Block 0 cell size :
i : 0.01999999999999996 .. 0.02000000000000002
j : 0.02000000000000007 .. 0.01999999999999996
k : 0.02 .. 0.02


There are no merge patch pairs

Writing polyMesh with 0 cellZones
----------------
Mesh Information
----------------
boundingBox: (-0.4000000000000001 -0.4000000000000001 -0.01) (1 0.4000000000000001 0.01)
nPoints: 5822
nCells: 2800
nFaces: 11310
nInternalFaces: 5490
----------------
Patches
----------------
patch 0 (start: 5490 size: 70) name: top
patch 1 (start: 5560 size: 40) name: inlet
patch 2 (start: 5600 size: 40) name: outlet
patch 3 (start: 5640 size: 70) name: Bottom
patch 4 (start: 5710 size: 5600) name: frontAndBack

End

myfoam@DESKTOP-TK3D7CI:~/OpenFOAM-sims/tutorials/incompressible/pimpleFoam/laminar/cylinder2D$ surfaceFeatureExtract
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2112 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : _6e1fca0e-20220610 OPENFOAM=2112 patch=220610 version=2112
Arch : "LSB;label=32;scalar=64"
Exec : surfaceFeatureExtract
Date : Aug 10 2023
Time : 21:50:23
Host : DESKTOP-TK3D7CI
PID : 1571
I/O : uncollated
Case : /home/myfoam/OpenFOAM-sims/tutorials/incompressible/pimpleFoam/laminar/cylinder2D
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time


Note: Feature line extraction only valid on closed manifold surfaces

Reading surfaceFeatureExtractDict

Surfaces : cylinder.stl
Output : cylinder
Load options : offset
Write options: writeObj=off writeVTK=off

Statistics:
Triangles : 360 in 1 region(s)
Vertices : 182
Bounding Box : (-20 -19.98781585693359 0) (20 19.98781585693359 10)

Constructing feature set from included angle 150
Selecting edges based purely on geometric tests: no

Initial Feature set:
points : 0
edges : 182
of which
region edges : 0
external edges : 182
internal edges : 0

Final Feature set:
points : 182
of which
convex : 0
concave : 0
mixed : 0
non-feature : 182
edges : 182
of which
external edges : 182
internal edges : 0
flat edges : 0
open edges : 0
multiply connected : 0

Writing extendedFeatureEdgeMesh to "/home/myfoam/OpenFOAM-sims/tutorials/incompressible/pimpleFoam/laminar/cylinder2D/constant/extendedFeatureEdgeMesh/cylinder.extendedFeatureEdgeMesh"

Writing featureEdgeMesh to "/home/myfoam/OpenFOAM-sims/tutorials/incompressible/pimpleFoam/laminar/cylinder2D/constant/triSurface/cylinder.eMesh"
Surfaces : fluid.stl
Output : fluid
Load options : offset
Write options: writeObj=off writeVTK=off

Statistics:
Triangles : 380 in 1 region(s)
Vertices : 190
Bounding Box : (-400 -400 0) (1000 400 10)

Constructing feature set from included angle 150
Selecting edges based purely on geometric tests: no

Initial Feature set:
points : 8
edges : 194
of which
region edges : 0
external edges : 194
internal edges : 0

Final Feature set:
points : 190
of which
convex : 8
concave : 0
mixed : 0
non-feature : 182
edges : 194
of which
external edges : 194
internal edges : 0
flat edges : 0
open edges : 0
multiply connected : 0

Writing extendedFeatureEdgeMesh to "/home/myfoam/OpenFOAM-sims/tutorials/incompressible/pimpleFoam/laminar/cylinder2D/constant/extendedFeatureEdgeMesh/fluid.extendedFeatureEdgeMesh"

Writing featureEdgeMesh to "/home/myfoam/OpenFOAM-sims/tutorials/incompressible/pimpleFoam/laminar/cylinder2D/constant/triSurface/fluid.eMesh"
ExecutionTime = 0.03 s ClockTime = 0 s

End
Attached Files
File Type: zip cylinder.zip (167.7 KB, 7 views)
saeed jamshidi is offline   Reply With Quote

Old   August 10, 2023, 14:34
Default
  #2
Senior Member
 
Saeed Jamshidi
Join Date: Aug 2019
Posts: 214
Rep Power: 7
saeed jamshidi is on a distinguished road
myfoam@DESKTOP-TK3D7CI:~/OpenFOAM-sims/tutorials/incompressible/pimpleFoam/laminar/cylinder2D$ snappyHexMesh


\*---------------------------------------------------------------------------*/
Build : _6e1fca0e-20220610 OPENFOAM=2112 patch=220610 version=2112
Arch : "LSB;label=32;scalar=64"
Exec : snappyHexMesh
Date : Aug 10 2023
Time : 21:50:29
Host : DESKTOP-TK3D7CI
PID : 1572
I/O : uncollated
Case : /home/myfoam/OpenFOAM-sims/tutorials/incompressible/pimpleFoam/laminar/cylinder2D
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Read mesh in = 0.05 s

Overall mesh bounding box : (-0.4000000000000001 -0.4000000000000001 -0.01) (1 0.4000000000000001 0.01)
Relative tolerance : 1e-06
Absolute matching distance : 1.612575579624099e-06

Reading refinement surfaces.
--> FOAM Warning :
From Foam::refinementSurfaces::refinementSurfaces(const Foam::searchableSurfaces&, const Foam::dictionary&, Foam::label, bool)
in file refinementSurfaces/refinementSurfaces.C at line 467
Reading "system/snappyHexMeshDict.castellatedMeshControls.refineme ntSurfaces" at line 72
Not all entries in refinementSurfaces dictionary were used. The following entries were not used : 1(solid)
Read refinement surfaces in = 0 s

Reading refinement shells.
Read refinement shells in = 0 s

Setting refinement level of surface to be consistent with shells.
Checked shell refinement in = 0 s

Reading features.
Read features in = 0.01 s


Determining initial surface intersections
-----------------------------------------

Edge intersection testing:
Number of edges : 11310
Number of edges to retest : 0
Number of intersected edges : 0
Edge intersection testing:
Number of edges : 11310
Number of edges to retest : 11310
Number of intersected edges : 0
Calculated surface intersections in = 0 s

Initial mesh : cells:2800 faces:11310 points:5822
Cells per refinement level:
0 2800

Adding patches for surface regions
----------------------------------

Patch Type Region
----- ---- ------
Added patches in = 0 s

Edge intersection testing:
Number of edges : 11310
Number of edges to retest : 0
Number of intersected edges : 0
Decomposition method none [1]

Refinement phase
----------------

Found point (-0.25 0.0001 0.0001) in cell 1407 on processor 0

Surface refinement iteration 0
------------------------------

Marked for refinement due to surface intersection : 0 cells.
Determined cells to refine in = 0 s
Selected for refinement : 0 cells (out of 2800)
Stopping refining since too few cells selected.


Removing mesh beyond surface intersections
------------------------------------------

Setting cellZones according to locationsInMesh:
Location : (-0.25 0.0001 0.0001)
cellZone : none

For cellZone none found point (-0.25 0.0001 0.0001) in global region 0 out of 1 regions.
Keeping all cells containing inside points
Selected for keeping : 2800 cells.
Edge intersection testing:
Number of edges : 11310
Number of edges to retest : 0
Number of intersected edges : 0
Setting cellZones according to locationsInMesh:
Location : (-0.25 0.0001 0.0001)
cellZone : none

For cellZone none found point (-0.25 0.0001 0.0001) in global region 0 out of 1 regions.
Keeping all cells containing inside points
Selected for keeping : 2800 cells.
Edge intersection testing:
Number of edges : 11310
Number of edges to retest : 0
Number of intersected edges : 0

Big gap refinement iteration 0
------------------------------

Determined cells to refine in = 0.03 s
Selected for refinement : 0 cells (out of 2800)
Stopping refining since too few cells selected.


Shell refinement iteration 0
----------------------------

Marked for refinement due to distance to explicit features : 0 cells.
Marked for refinement due to refinement shells : 0 cells.
Determined cells to refine in = 0 s
Selected for internal refinement : 0 cells (out of 2800)
Stopping refining since too few cells selected.


Dangling coarse cells refinement iteration 0
--------------------------------------------

Determined cells to refine in = 0.01 s
Selected for refinement : 0 cells (out of 2800)
Stopping refining since too few cells selected.


Dangling coarse cells refinement iteration 0
--------------------------------------------

Determined cells to refine in = 0 s
Selected for refinement : 0 cells (out of 2800)
Stopping refining since too few cells selected.


Refinement transition refinement iteration 0
--------------------------------------------

Determined cells to refine in = 0 s
Selected for refinement : 0 cells (out of 2800)
Stopping refining since too few cells selected.


Splitting mesh at surface intersections
---------------------------------------

Introducing baffles for 0 faces that are intersected by the surface.

Setting cellZones according to locationsInMesh:
Location : (-0.25 0.0001 0.0001)
cellZone : none

For cellZone none found point (-0.25 0.0001 0.0001) in global region 0 out of 1 regions.
Created baffles in = 0 s


After introducing baffles : cells:2800 faces:11310 points:5822
Cells per refinement level:
0 2800

Introducing baffles to block off problem cells
----------------------------------------------

markFacesOnProblemCells : marked 0 additional internal faces to be converted into baffles.
Analyzed problem cells in = 0.02 s


Introducing baffles to delete problem cells.

Created baffles in = 0 s


After introducing baffles : cells:2800 faces:11310 points:5822
Cells per refinement level:
0 2800
Setting cellZones according to locationsInMesh:
Location : (-0.25 0.0001 0.0001)
cellZone : none

For cellZone none found point (-0.25 0.0001 0.0001) in global region 0 out of 1 regions.

Remove unreachable sections of mesh
-----------------------------------

Split mesh in = 0 s


After subsetting : cells:2800 faces:11310 points:5822
Cells per refinement level:
0 2800

Introducing zones for interfaces
--------------------------------

Setting cellZones according to locationsInMesh:
Location : (-0.25 0.0001 0.0001)
cellZone : none

For cellZone none found point (-0.25 0.0001 0.0001) in global region 0 out of 1 regions.
Setting faceZones according to neighbouring cellZones:

Edge intersection testing:
Number of edges : 11310
Number of edges to retest : 0
Number of intersected edges : 0

Handling cells with snap problems
---------------------------------

Introducing baffles for 0 faces that are intersected by the surface.

Setting cellZones according to locationsInMesh:
Location : (-0.25 0.0001 0.0001)
cellZone : none

For cellZone none found point (-0.25 0.0001 0.0001) in global region 0 out of 1 regions.
Created baffles in = 0 s


After introducing baffles : cells:2800 faces:11310 points:5822
Cells per refinement level:
0 2800

Introducing baffles to block off problem cells
----------------------------------------------

markFacesOnProblemCells : marked 0 additional internal faces to be converted into baffles.
Analyzed problem cells in = 0 s


Introducing baffles to delete problem cells.

Created baffles in = 0 s


After introducing baffles : cells:2800 faces:11310 points:5822
Cells per refinement level:
0 2800
Setting cellZones according to locationsInMesh:
Location : (-0.25 0.0001 0.0001)
cellZone : none

For cellZone none found point (-0.25 0.0001 0.0001) in global region 0 out of 1 regions.

Remove unreachable sections of mesh
-----------------------------------

Split mesh in = 0 s


After subsetting : cells:2800 faces:11310 points:5822
Cells per refinement level:
0 2800

Merge free-standing baffles
---------------------------

freeStandingBaffles : detected 0 free-standing baffles out of 0

Detected free-standing baffles : 0
Merged free-standing baffles in = 0 s


dupNonManifoldPoints : Found : 0 non-manifold points (out of 5822)
Detected unsplittable baffles : 0

Merge refined boundary faces
----------------------------

Merging 0 sets of faces.

No faces merged ...

Merging all points on surface that
- are used by only two boundary faces and
- make an angle with a cosine of more than 0.7071067811865476.

No straight edges simplified and no points removed ...
Refined mesh : cells:2800 faces:11310 points:5822
Cells per refinement level:
0 2800
Writing mesh to time 0.05
Wrote mesh in = 0.01 s.
Mesh refined in = 0.07000000000000001 s.

Morphing phase
--------------

Duplicating 0 points on faceZones of type boundary
Constructing mesh displacer ...
Using mesh parameters
{
maxNonOrtho 65;
maxBoundarySkewness 20;
maxInternalSkewness 4;
maxConcave 80;
minVol 1e-18;
minTetQuality 1e-09;
minArea -1;
minTwist 0.05;
minDeterminant 0.001;
minFaceWeight 0.05;
minVolRatio 0.01;
minTriangleTwist -1;
nSmoothScale 4;
errorReduction 0.75;
relaxed
{
maxNonOrtho 75;
}
}


Checking initial mesh ...
Checking faces in error :
non-orthogonality > 65 degrees : 0
faces with face pyramid volume < 1e-18 : 0
faces with face-decomposition tet quality < 1e-09 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Detected 0 illegal faces (concave, zero area or negative cell pyramid volume)
Checked initial mesh in = 0.02 s


Smoothing patch points ...
Smoothing iteration 0
Found 0 non-manifold point(s).

Scaling iteration 0
Moving mesh using displacement scaling : min:1 max:1
Correcting 2-D mesh motion ...done
Checking faces in error :
non-orthogonality > 65 degrees : 0
faces with face pyramid volume < 1e-18 : 0
faces with face-decomposition tet quality < 1e-09 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Successfully moved mesh

Smoothing iteration 1
Found 0 non-manifold point(s).

Scaling iteration 0
Moving mesh using displacement scaling : min:1 max:1
Correcting 2-D mesh motion ...done
Checking faces in error :
non-orthogonality > 65 degrees : 0
faces with face pyramid volume < 1e-18 : 0
faces with face-decomposition tet quality < 1e-09 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Successfully moved mesh

Smoothing iteration 2
Found 0 non-manifold point(s).

Scaling iteration 0
Moving mesh using displacement scaling : min:1 max:1
Correcting 2-D mesh motion ...done
Checking faces in error :
non-orthogonality > 65 degrees : 0
faces with face pyramid volume < 1e-18 : 0
faces with face-decomposition tet quality < 1e-09 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Successfully moved mesh

Patch points smoothed in = 0.24 s



Morph iteration 0
-----------------
Calculating patchDisplacement as distance to nearest surface point ...
non-zone points : attract to nearest of all non-zone surfaces
zone points : attract to zone surface only

Calculated surface displacement in = 0 s


Detecting near surfaces ...
Overriding nearest with intersection of close gaps at 0 out of 0 points.
Smoothing displacement ...
Iteration 0
Iteration 10
Iteration 20
Displacement smoothed in = 0.01 s


Moving mesh ...

Iteration 0
Moving mesh using displacement scaling : min:1 max:1
Correcting 2-D mesh motion ...done
Checking faces in error :
non-orthogonality > 65 degrees : 0
faces with face pyramid volume < 1e-18 : 0
faces with face-decomposition tet quality < 1e-09 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Successfully moved mesh
Moved mesh in = 0.08 s


Repatching faces according to nearest surface ...
Repatched 0 faces in = 0.02 s


Edge intersection testing:
Number of edges : 11310
Number of edges to retest : 11310
Number of intersected edges : 0
Merging 0 sets of faces.

No faces merged ...

Merging all points on surface that
- are used by only two boundary faces and
- make an angle with a cosine of more than 0.8660254037844387.

No straight edges simplified and no points removed ...
Snapped mesh : cells:2800 faces:11310 points:5822
Cells per refinement level:
0 2800
Writing mesh to time 0.1
Wrote mesh in = 0.04 s.
Mesh snapped in = 0.41 s.
Layer specification as
- final layer thickness ('finalLayerThickness')
- expansion ratio ('expansionRatio')
--> FOAM Warning :
From Foam::labelHashSet Foam::polyBoundaryMesh::patchSet(const Foam::UList<Foam::wordRe>&, bool, bool) const
in file meshes/polyMesh/polyBoundaryMesh/polyBoundaryMesh.C at line 915
Cannot find any patch or group names matching cylinder
--> FOAM Warning :
From Foam::layerParameters::layerParameters(const Foam::dictionary&, const Foam::polyBoundaryMesh&, bool)
in file snappyHexMeshDriver/layerParameters/layerParameters.C at line 540
Reading "system/snappyHexMeshDict.addLayersControls.layers" at line 131
Layer specification for cylinder does not match any patch.
Valid patches are 5(top inlet outlet Bottom frontAndBack)

Shrinking and layer addition phase
----------------------------------

Using mesh parameters
{
maxNonOrtho 65;
maxBoundarySkewness 20;
maxInternalSkewness 4;
maxConcave 80;
minVol 1e-18;
minTetQuality 1e-09;
minArea -1;
minTwist 0.05;
minDeterminant 0.001;
minFaceWeight 0.05;
minVolRatio 0.01;
minTriangleTwist -1;
nSmoothScale 4;
errorReduction 0.75;
relaxed
{
maxNonOrtho 75;
}
}



Merging all faces of a cell
---------------------------
- which are on the same patch
- which make an angle < 30 degrees (cos:0.8660254037844387)
- as long as the resulting face doesn't become concave by more than 90 degrees
(0=straight, 180=fully concave)

Merging 0 sets of faces.

No faces merged ...

Merging all points on surface that
- are used by only two boundary faces and
- make an angle with a cosine of more than 0.8660254037844387.

No straight edges simplified and no points removed ...

No layers to generate ...
Layer mesh : cells:2800 faces:11310 points:5822
Cells per refinement level:
0 2800
Writing mesh to time 0.15
Wrote mesh in = 0.05 s.
Layers added in = 0.05 s.
Checking final mesh ...
Checking faces in error :
non-orthogonality > 65 degrees : 0
faces with face pyramid volume < 1e-18 : 0
faces with face-decomposition tet quality < 1e-09 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Finished meshing without any errors
Finished meshing in = 0.61 s.
End

myfoam@DESKTOP-TK3D7CI:~/OpenFOAM-sims/tutorials/incompressible/pimpleFoam/laminar/cylinder2D$ pimpleFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2112 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : _6e1fca0e-20220610 OPENFOAM=2112 patch=220610 version=2112
Arch : "LSB;label=32;scalar=64"
Exec : pimpleFoam
Date : Aug 10 2023
Time : 21:50:42
Host : DESKTOP-TK3D7CI
PID : 1573
I/O : uncollated
Case : /home/myfoam/OpenFOAM-sims/tutorials/incompressible/pimpleFoam/laminar/cylinder2D
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


PIMPLE: Operating solver in PISO mode

Reading field p



--> FOAM FATAL IO ERROR: (openfoam-2112 patch=220610)
Cannot find patchField entry for top

file: 0/p.boundaryField at line 25 to 46.

From void Foam::GeometricField<Type, PatchField, GeoMesh>::Boundary::readField(const Foam::DimensionedField<TypeR, GeoMesh>&, const Foam::dictionary&) [with Type = double; PatchField = Foam::fvPatchField; GeoMesh = Foam::volMesh]
in file /usr/src/packages/BUILD/src/OpenFOAM/lnInclude/GeometricBoundaryField.C at line 172.

FOAM exiting
saeed jamshidi is offline   Reply With Quote

Old   August 11, 2023, 07:42
Default
  #3
Senior Member
 
Yann
Join Date: Apr 2012
Location: France
Posts: 1,088
Rep Power: 26
Yann will become famous soon enough
Hello Saeed,

Firstly, I would recommend to use code tags when pasting your terminal outputs here, it will make things way easier to read:

HTML Code:
[CODE]your terminal output[/CODE]
Now back to your question:

Code:
--> FOAM FATAL IO ERROR: (openfoam-2112 patch=220610)
Cannot find patchField entry for top

file: 0/p.boundaryField at line 25 to 46.
pimpleFoam complains because it cannot find a boundary condition for the patch named "top" in the file 0/p
Your mesh contains a patch named "top" and you need to define boundary conditions for it.
Actually you need to define boundary conditions for all variables in the 0 directory, for all patches in your mesh.

Before trying to fix boundary conditions, I think you have a problem with your mesh. When running snappyHexMesh, you get this:
Code:
Reading refinement surfaces.
--> FOAM Warning :
From Foam::refinementSurfaces::refinementSurfaces(const Foam::searchableSurfaces&, const Foam::dictionary&, Foam::label, bool)
in file refinementSurfaces/refinementSurfaces.C at line 467
Reading "system/snappyHexMeshDict.castellatedMeshControls.refineme ntSurfaces" at line 72
Not all entries in refinementSurfaces dictionary were used. The following entries were not used : 1(solid)
Read refinement surfaces in = 0 s
Which is basically telling you defined some parameters in refinementSurfaces which are not used by snappy. It means there is probably a mistake in your snappyHexMeshDict file.

Then when checking the rest of the snappy log, it basically does nothing to your initial mesh: blockMesh creates a 2800 cells mesh and snappy does not seem to create any cells:

Code:
Big gap refinement iteration 0
------------------------------

Determined cells to refine in = 0.03 s
Selected for refinement : 0 cells (out of 2800)
Stopping refining since too few cells selected.


Shell refinement iteration 0
----------------------------

Marked for refinement due to distance to explicit features : 0 cells.
Marked for refinement due to refinement shells : 0 cells.
Determined cells to refine in = 0 s
Selected for internal refinement : 0 cells (out of 2800)
Stopping refining since too few cells selected.


Dangling coarse cells refinement iteration 0
--------------------------------------------

Determined cells to refine in = 0.01 s
Selected for refinement : 0 cells (out of 2800)
Stopping refining since too few cells selected.

You need to go step by step: first deal with snappyHexMesh: run it, check the result in paraView, make sure everything is fine.
Once you are satisfied with the mesh, you can move on setting up boundary conditions and running pimpleFoam.

You can use user guide and tutorials to have some examples on how to use snappyHexMesh and setup a case.

Regards,
Yann
Yann is offline   Reply With Quote

Old   August 11, 2023, 15:16
Default
  #4
Senior Member
 
Saeed Jamshidi
Join Date: Aug 2019
Posts: 214
Rep Power: 7
saeed jamshidi is on a distinguished road
Thank you dear Yann 🙏
Yann likes this.
saeed jamshidi is offline   Reply With Quote

Reply

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
Simulation DNS for two phase flow in openFoam ht2017 OpenFOAM Running, Solving & CFD 0 December 8, 2020 05:50
Drag force coefficient too low for a flow past cylinder at Re= 1e05 Scabbard STAR-CCM+ 2 June 5, 2020 14:44
Flow past a cylinder at Re 1e05 using LES, drag force coefficient is to low Scabbard Main CFD Forum 21 June 19, 2018 13:58
[Other] OpenFoam Flow over a Circular Cylinder WolfgangS. OpenFOAM Meshing & Mesh Conversion 12 March 3, 2014 10:53
benchmark: flow over a circular cylinder goodegg Main CFD Forum 12 January 22, 2013 11:47


All times are GMT -4. The time now is 10:57.