|
[Sponsors] |
August 10, 2023, 14:33 |
flow over cylinder in openFoam
|
#1 |
Senior Member
Saeed Jamshidi
Join Date: Aug 2019
Posts: 214
Rep Power: 8 |
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 |
|
August 10, 2023, 14:34 |
|
#2 |
Senior Member
Saeed Jamshidi
Join Date: Aug 2019
Posts: 214
Rep Power: 8 |
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 |
|
August 11, 2023, 07:42 |
|
#3 |
Senior Member
Yann
Join Date: Apr 2012
Location: France
Posts: 1,174
Rep Power: 27 |
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] Code:
--> FOAM FATAL IO ERROR: (openfoam-2112 patch=220610) Cannot find patchField entry for top file: 0/p.boundaryField at line 25 to 46. 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 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 |
|
August 11, 2023, 15:16 |
|
#4 |
Senior Member
Saeed Jamshidi
Join Date: Aug 2019
Posts: 214
Rep Power: 8 |
Thank you dear Yann 🙏
|
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |