|
[Sponsors] |
[Commercial meshers] Handling cyclic BC from gambit to openfoam for a cascade airfoil problem - OF 1.6 |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
April 13, 2010, 22:32 |
Handling cyclic BC from gambit to openfoam for a cascade airfoil problem - OF 1.6
|
#1 |
New Member
maverick
Join Date: Feb 2010
Posts: 9
Rep Power: 16 |
Dear all,
after successfully installing openfoam and after doing some tutorials, i am doing a problem on cascade flow analysis of a transonic turbine . I am having problem in creating the cyclic patch. i have given below the procedure i have followed.There were 3-4 threads i have referred and i have implemented the ideas given in them.still i have problem when i use Paraview. Problem description : My problem consists of blade centered domain in which there is airfoil in center, inlet to left of airfoil,outlet to the right , one periodic surface on top and other periodic surface on bottom and airfoil itself is the wall 1. I have created the mesh in gambit, i have linked the upper periodic edge and lower periodic edge meshes(from thread).While creating BC i have given them as " wall " with names wall_pressure and wall_suction respectively.exported the 2d mesh. 2.imported the mesh in openfoam with fluentMeshToFoam and before importing i have set the writeprecision in controldict to 12 (from thread ).after importing i gave checkmesh and found ok.Then i viewed the mesh in paraview after making the changes in p and u files . 3.Then i made a createpatchdict in system folder with the following lines in it (from thread) Code:
matchTolerance 0.2; pointSync true; patches ( { name wall_pressure; type cyclic; constructFrom patches; patches (wall_pressure wall_suction); } ); 4.I ran createPatch in terminal with matchtolerance set at 0.2 (from thread).I got the following result Code:
Create time Reading createPatchDict. Using relative tolerance 0.2 to match up faces and points Create polyMesh for time = 0 Moving faces from patch wall_pressure to patch 1 Moving faces from patch wall_suction to patch 1 Doing topology modification to order faces. Synchronising points. Synchronising points. Points changed by average:0 max:0 Removing patches with no faces in them. Removing empty patch wall_suction at position 2 Removing patches. Writing repatched mesh to 4e-08 End 6.reran the createPatch with the following changes Code:
matchTolerance 0.2; pointSync true; patches ( { name wall_pressure; type cyclic; constructFrom patches; patches (wall_pressure); } ); Code:
Create time Reading createPatchDict. Using relative tolerance 0.2 to match up faces and points Create polyMesh for time = 4e-08 Dumping cycPatch.name() half0 faces to "initial_wall_pressure_half0.obj" Dumping cycPatch.name() half1 faces to "initial_wall_pressure_half1.obj" Dumping cyclic match as lines between face centres to "initial_wall_pressure_match.obj" Moving faces from patch wall_pressure to patch 1 Doing topology modification to order faces. cyclicPolyPatch::order : Writing half0 faces to OBJ file "wall_pressure_half0_faces.obj" cyclicPolyPatch::order : Writing half1 faces to OBJ file "wall_pressure_half1_faces.obj" cyclicPolyPatch::order : Dumping currently found cyclic match as lines between corresponding face centres to file "/home/kaka/openfoam_validation/tutorials/compressible/sonicFoam/ras/nacaAirfoil/wall_pressure_faceCentres.obj" --> FOAM Serious Error : From function cyclicPolyPatch::order(const primitivePatch&, labelList&, labelList&) const in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 1547 Patch:wall_pressure : Cannot match vectors to faces on both sides of patch Perhaps your faces do not match? The obj files written contain the current match. Continuing with incorrect face ordering from now on! Dumping cycPatch.name() half0 faces to "coupled_wall_pressure_half0.obj" Dumping cycPatch.name() half1 faces to "coupled_wall_pressure_half1.obj" Dumping cyclic match as lines between face centres to "coupled_wall_pressure_match.obj" Synchronising points. On coupled patch wall_pressure forcing uniform rotation of (1 -0 0 0 1 0 0 0 1) On coupled patch wall_pressure forcing uniform rotation of 1((1 -0 0 0 1 0 0 0 1)) Synchronising points. --> FOAM Warning : From function syncPoints(const polyMesh&, pointField&, const CombineOp&, const point&) in file createPatch.C at line 483 There are decomposed cyclics in this mesh with transformations. This is not supported. The result will be incorrect Points changed by average:0.410037731332 max:105.851509145 face 2 area does not match neighbour 108 by 200% -- possible face ordering problem. patch:wall_pressure my area:1.02186e-17 neighbour area:63.8042 matching tolerance:0.2 Mesh face:28987 vertices:4((52.6199 -49.349 1.05888) (52.6199 -49.349 -1.05888) (51.6153 -49.349 -1.05888) (52.6199 -49.349 -1.05888)) Neighbour face:29093 vertices:4((53.6245 -49.349 -1.05888) (-29.0738 14.1267 1.05888) (51.6153 -49.349 -1.05888) (52.6199 -49.349 -1.05888)) Rerun with cyclic debug flag set for more information. From function cyclicPolyPatch::calcTransforms() in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 180. FOAM exiting any help will be greatly appreciated |
|
July 7, 2010, 07:42 |
|
#2 |
New Member
Irin Sun
Join Date: Jan 2010
Posts: 5
Rep Power: 16 |
I have the same problem
|
|
June 18, 2011, 04:36 |
|
#3 |
New Member
Join Date: Oct 2010
Posts: 24
Rep Power: 15 |
Hi,maverick
Have you solved this problem? I have the same problem |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[Other] Cyclic patch problem! | SMartella | OpenFOAM Meshing & Mesh Conversion | 7 | May 6, 2024 08:03 |
Problem with airfoil shape optimization | robyTKD | SU2 Shape Design | 7 | March 7, 2022 16:18 |
Problem with using oscillatingFixedvalue in OpenFoam 1.6 | markong | OpenFOAM | 3 | November 5, 2020 02:56 |
Problem with Acoustic Pulse Simulation | Weng-Ji | OpenFOAM Running, Solving & CFD | 0 | November 6, 2016 04:40 |
High drag for airfoil compared to XFOIL and wind tunnel data | Ry10 | SU2 | 15 | October 30, 2016 17:27 |