CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ...

Handling cyclic BC from gambit to openfoam for a cascade airfoil problem - OF 1.6

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

Like Tree1Likes
  • 1 Post By maverick

Reply
 
LinkBack Thread Tools Display Modes
Old   April 13, 2010, 22:32
Default 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: 7
maverick is on a distinguished road
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);
}
);
here wall_pressure and wall_suction are the periodic edges

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
5. i got a boundary file in the directory 4e-08/polymesh/boundary. i changed the type of wall_pressure from wall to cyclic.copied this boundary file and pasted this file in the constant/polymesh directory replacing the existing boundary file.

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);
}
);
7.Now i get the following error when i run createPatch
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
i am not much familiar with linux and openFOAM,
any help will be greatly appreciated
mm.abdollahzadeh likes this.
maverick is offline   Reply With Quote

Old   July 7, 2010, 07:42
Default
  #2
New Member
 
Irin Sun
Join Date: Jan 2010
Posts: 5
Rep Power: 7
irinsun is on a distinguished road
I have the same problem
irinsun is offline   Reply With Quote

Old   June 18, 2011, 04:36
Default
  #3
New Member
 
Join Date: Oct 2010
Posts: 24
Rep Power: 6
LijieNPIC is on a distinguished road
Hi,maverick
Have you solved this problem? I have the same problem
LijieNPIC is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Problem with using oscillatingFixedvalue in OpenFoam 1.6 markong OpenFOAM 2 September 16, 2010 06:05
kOmegaSST in openfoam 1.6 Gearb0x OpenFOAM 2 March 3, 2010 07:02
Sample Utility not working in OpenFoam 1.6 titio OpenFOAM Post-Processing 0 February 5, 2010 13:12
Cross-compiling OpenFOAM 1.6 on Linux for Windows 32 and 64bits with Mingw-w64 wyldckat OpenFOAM Announcements from Other Sources 7 January 19, 2010 16:39
GAMBIT meshing problem Gauthier Lambert Main CFD Forum 1 August 3, 2000 09:22


All times are GMT -4. The time now is 03:56.