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

Problems with cyclicBC

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 16, 2010, 03:49
Default Problems with cyclicBC
  #1
Member
 
Michael
Join Date: Mar 2009
Posts: 48
Rep Power: 8
farbfilm is on a distinguished road
Hi,

I'm trying to simulate one channel of a turbine! For this I need the cyclicBC.
I've done a periodic mesh with Icem, converted it with FluentMeshToFoam and want to use the MRFSimpleFoam-Solver.

Whatever I'm doing to setup the cyclic BC, I'm always getting this error:

--> FOAM Serious Error :
From function void cyclicPolyPatch::calcTransforms()
in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 279
Transformation tensor is not constant for the cyclic patch. Please reconsider your setup and definition of cyclic boundaries.

I checked my setup and searched the forum, but I don't get an idea what's wrong!

Can somebody help me???


Thanks, Michael
farbfilm is offline   Reply With Quote

Old   November 16, 2010, 04:51
Default
  #2
Senior Member
 
Steven van Haren
Join Date: Aug 2010
Location: The Netherlands
Posts: 149
Rep Power: 7
stevenvanharen is on a distinguished road
I assume you get this error when running createPatch?

It looks like OpenFOAM thinks the two faces are rotated with respect to eachother, are they in your geometry? And in that case did you specify a rotation axis in your dictionary?
stevenvanharen is offline   Reply With Quote

Old   November 16, 2010, 05:10
Default
  #3
Member
 
Michael
Join Date: Mar 2009
Posts: 48
Rep Power: 8
farbfilm is on a distinguished road
Hi,

thank your for your answer!

Yes, I get this error when I run createPatch, but also when I start the calculation (the calculation starts in spite of this error)!

Sorry, I don't understand your first question! What so you mean with 'are they in your geometry'??
I just specified an axis in the MRF-Zones-File! Where else do I have to specify a rotation axis in my case??

Thank you again!!

michael
farbfilm is offline   Reply With Quote

Old   November 16, 2010, 05:15
Default
  #4
Senior Member
 
Steven van Haren
Join Date: Aug 2010
Location: The Netherlands
Posts: 149
Rep Power: 7
stevenvanharen is on a distinguished road
If the two faces from which you want to create the cyclic interface are rotated with respect to each other you have to specify a rotation axis in the createPatchDict file.

Take a look in the folder where createPatch.C is located (the source), you will find and example for the createPatchDict there.
stevenvanharen is offline   Reply With Quote

Old   November 16, 2010, 05:50
Default
  #5
Member
 
Michael
Join Date: Mar 2009
Posts: 48
Rep Power: 8
farbfilm is on a distinguished road
I checked the createPatchDict from /OpenFOAM/OpenFOAM-1.5-dev/applications/utilities/mesh/manipulation/createPatch (see attachment) but couldn't find an entry where I can define a rotation axis!
Attached Files
File Type: txt createPatchDict.txt (1.8 KB, 36 views)
farbfilm is offline   Reply With Quote

Old   November 16, 2010, 06:02
Default
  #6
Senior Member
 
Steven van Haren
Join Date: Aug 2010
Location: The Netherlands
Posts: 149
Rep Power: 7
stevenvanharen is on a distinguished road
Ah, mmmm.

It is possible that this is not supported in OF 1.5, however I don't know.

This is the dictionary from 1.7. You can try if it works with 1.5
Attached Files
File Type: txt createPatchDict.txt (3.1 KB, 35 views)
stevenvanharen is offline   Reply With Quote

Old   November 16, 2010, 06:19
Default
  #7
Member
 
Michael
Join Date: Mar 2009
Posts: 48
Rep Power: 8
farbfilm is on a distinguished road
Unfortunately, it doesn't work with 1.5!

Is there another solution??
(I need 1.5-dev, cause I want to use the GGI and the turbomachinery tools later)


Thanks!!!

michael
farbfilm is offline   Reply With Quote

Old   November 16, 2010, 06:23
Default
  #8
Senior Member
 
Steven van Haren
Join Date: Aug 2010
Location: The Netherlands
Posts: 149
Rep Power: 7
stevenvanharen is on a distinguished road
Sorry, in that case I can not help you any further.
stevenvanharen is offline   Reply With Quote

Old   November 16, 2010, 06:31
Default
  #9
Member
 
Michael
Join Date: Mar 2009
Posts: 48
Rep Power: 8
farbfilm is on a distinguished road
Okay, however: Thank you very much!!

Is there anybody else who could help me??
Maybe there's somebody who is dealing with turbomachinery as well, simulates just one channel instead of the whole rotor and can help me with the general procedure of such applications (I'm a lost beginner...)??

Thank you!!

michael
farbfilm is offline   Reply With Quote

Old   November 17, 2010, 06:01
Default
  #10
Senior Member
 
maddalena's Avatar
 
maddalena
Join Date: Mar 2009
Posts: 436
Rep Power: 12
maddalena is on a distinguished road
Hi Michael,
not sure about your problem (I am using 1.7), but this:
Quote:
Originally Posted by farbfilm View Post
Transformation tensor is not constant for the cyclic patch. Please reconsider your setup and definition of cyclic boundaries.
let me think that maybe the two patches defining the cyclic are not flat, thus the transformation tensor differs from one point to the other...
I also have noticed that you use the same patches twice in your createPatchDict: you have the half0 half1 couple both when defining leftRight0 and bottom patches! This may lead the utility to some errors. Try to define only one of the two.
Hope this help,
cheers

mad
maddalena is offline   Reply With Quote

Old   November 17, 2010, 06:31
Default
  #11
Member
 
Michael
Join Date: Mar 2009
Posts: 48
Rep Power: 8
farbfilm is on a distinguished road
Hi Mad,
thank you for your answer!

The createPatchDict I posted earlier is just an example from the source of OpenFOAM 1.5-dev.
The createPatchDict that I used is attached now.

In the picture you see that my cyclic patches are flat. So, my geometry shouldn't be a problem...

michael
Attached Images
File Type: jpg mesh.jpg (41.1 KB, 51 views)
Attached Files
File Type: txt createPatchDict.txt (1.2 KB, 26 views)
farbfilm is offline   Reply With Quote

Old   November 17, 2010, 08:26
Default
  #12
Senior Member
 
maddalena's Avatar
 
maddalena
Join Date: Mar 2009
Posts: 436
Rep Power: 12
maddalena is on a distinguished road
Hi,
Quote:
Originally Posted by farbfilm View Post
The createPatchDict that I used is attached now.
That seems ok.
Quote:
Originally Posted by farbfilm View Post
In the picture you see that my cyclic patches are flat. So, my geometry shouldn't be a problem...
Ok, the looks flat, but are they really flat? What I want to say is that there can be some sort of distorted cells that does not belong to the plane... Only an idea, though.

mad
maddalena is offline   Reply With Quote

Old   November 17, 2010, 10:49
Default
  #13
New Member
 
Join Date: Feb 2010
Posts: 26
Rep Power: 7
Pekka is on a distinguished road
Hi Michael,

Are you try run renumber mesh on ICEMCfd? I think that node number must be in a consecutive order on same patch.

BR/Pekka
Pekka is offline   Reply With Quote

Old   November 18, 2010, 03:22
Default
  #14
Member
 
Michael
Join Date: Mar 2009
Posts: 48
Rep Power: 8
farbfilm is on a distinguished road
Unfortunately renumbering doesn't help, too! But thanks for that hint!

When I change the type of my periodic patches from 'patch' to 'cyclic' before using createPatch, I'm getting the following error when running createPatch:

Code:
Create time

Reading createPatchDict.

Using relative tolerance 1e-06 to match up faces and points

Create polyMesh for time = 0



face 4215 area does not match neighbour 8431 by 11.4852% -- possible face ordering problem.
patch:PERIODIC2 my area:1.59926e-07 neighbour area:1.79413e-07 matching tolerance:0.001
Mesh face:3840485 vertices:4((0.148545 0.13087 0.00101004) (0.148056 0.130485 0.00101092) (0.148056 0.130485 0.00126784) (0.148545 0.13087 0.00126674))
Neighbour face:3844701 vertices:4((0.149363 0.131512 0.834) (0.149358 0.131508 0.834) (0.149354 0.131506 0.864) (0.149359 0.131509 0.864))
Other errors also exist, only the largest is reported. Please rerun with cyclic debug flag set for more information.

    From function cyclicPolyPatch::calcTransforms()
    in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 197.

FOAM exiting
When I'm using the procedure Fabian (fs82), here Problem with cyclic boundaries in Openfoam 1.5 ,I'm getting the error I mentioned above!

I'm going round in circles...

michael
farbfilm is offline   Reply With Quote

Old   November 18, 2010, 04:47
Default
  #15
Member
 
Rasoul
Join Date: Feb 2010
Posts: 32
Rep Power: 7
aut_iut is on a distinguished road
Hi Michael,
if you define your mesh by blockMeshDict, i think you need to change some thing in it.
put all of patches which have cyclic B.C together. It should fix your problem.
for example, for my case i need to apply cyclic B.C on the front and back sides.
So, under the patches name i define BackandFront and put all of the faces of these sides together.
you need to modify your boundary file accordingly.
I hope it could be useful for you.

Rasoul
aut_iut is offline   Reply With Quote

Old   November 18, 2010, 06:38
Default
  #16
New Member
 
Join Date: Feb 2010
Posts: 26
Rep Power: 7
Pekka is on a distinguished road
Are you run on Icem "check mesh" and "beriodic boblems" is switched on? Then set the wall boundary conditions to the cyclic pair so that you have two different wall patch. Then run renumbering. Export mesh to blah.msh. Try both mesh converter "fluentMeshToFoam" and "fluent3DMeshToFoam". Then "createPatch". If it's works, test following data your boundary file under to cyclic patch and try solve. Change the angle and axis to same than is yours model. Try also "pointSync true;" and different "matchTolerance" on createPatchDict.

Code:
        type            cyclic;
        transform       rotational;
        rotationAxis    (0 0 1);
        rotationCentre  (0 0 0);
        rotationAngle    6;
of

BR/Pekka
Pekka 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
Needed Benchmark Problems for FSI Mechstud Main CFD Forum 4 July 26, 2011 12:13
Two-phase air water flow problems by activating Wall Lubrication Force challenger85 CFX 5 November 5, 2009 06:44
Some problems with Star CD Micha CD-adapco 0 August 6, 2003 13:55
unstructured grid sreekanth Main CFD Forum 1 August 6, 2001 15:09
Inverse problems Aleksey Alekseev Main CFD Forum 0 May 12, 1999 15:38


All times are GMT -4. The time now is 00:47.