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

Possible bug in cyclicPolyPatch::order

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 26, 2009, 10:00
Default Possible bug in cyclicPolyPatch::order
  #1
Member
 
Niklas Wikstrom
Join Date: Mar 2009
Posts: 85
Rep Power: 8
wikstrom is on a distinguished road
Hi,

Running recent version of 1.6.x:

I am trying to create a (rotational) cyclic patch from two patches using createPatch. rotationAxis and rotationCentre is is given. Although several tests with different tolerances are performed, attempts fail with a Serious error in cyclicPolyPatch:rder, followed by a Fatal Error due to face area mismatch (log file attached).

A coarser mesh on the same geometry runs without errors and the problem is not easy to reproduce.

Funny thing is that the faces reported to be of different area in the FATAL ERROR, does not match the correspondingly reported vertices. The attached image show the vertex location as orange spheres and the faces location barely seen as green triangles. The faces are loaded into paraview from the faceSet "wrongFaces" generated by createPatch.

The wavefront .obj files XX_half0_faces.obj and XX_half1_faces.obj generated by createPatch does contain the correct faces. With debug flag set for cyclics the lines showing found face to face connections looks much like an inverted hedgehog (attached).

I am aware of ANSAS lack of fine tolerances, but I have made sure the patches are perfectly flat and aligned to respective coordinate planes. Moreover the surface mesh on the two patches are in perfect match, using "eye norm".

My test case is 3M compressed and I am more than happy to send a mail containing the failing case.

Best regards and warmest congratulations
Niklas
Attached Images
File Type: png wrongFaces_wrongVertices.png (69.1 KB, 24 views)
File Type: jpg hedgeHog.jpg (93.3 KB, 20 views)
Attached Files
File Type: txt createPatchLog.txt (2.8 KB, 5 views)
wikstrom is offline   Reply With Quote

Old   October 26, 2009, 13:33
Default
  #2
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 16
mattijs is on a distinguished road
Hi Niklas,

- can you switch on the debug flag for cyclic. Check the printout for 'Specified rotation' n0 and n1.
- check the obj files written for 'automatic ordering result'. Do they correspond to the two halves rotated on top of one another?

if you can create a smaller case I'm happy to have a quick look at it.
mattijs is offline   Reply With Quote

Old   October 27, 2009, 08:43
Default
  #3
Member
 
Niklas Wikstrom
Join Date: Mar 2009
Posts: 85
Rep Power: 8
wikstrom is on a distinguished road
Code:
Specified rotation : n0:(0 1 0) n1:(0 5.29112e-08 -1)
Which seems reasonable?

HOWEVER:
The problem is solved, for this test case, by setting the matchTolerance 1. It feels a bit awkward to allow such a high tolerance which is why I did not try it before. The points un-alignment is far less than the shortest face edge (~10) in this case, I am quite sure, and the point sync reports
Code:
Points changed by average:9.91435e-16 max:2.84237e-14
It remains to be seen if a usable, more complex geometry with cyclics can be created with Ansa.
Sloppy tolerances is a problem with Ansa, who on the other hand puts a nice effort into exporting directly to OpenFOAM.

Thank you Mattijs for answering! It helped me help myself :-)

Best regards!
wikstrom is offline   Reply With Quote

Old   October 27, 2009, 12:57
Default
  #4
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 16
mattijs is on a distinguished road
The problem is I guess due to an accumulation of errors. Normal is slightly off, transformed points (one half transformed on to other half) even more off and probably just beyond the tolerance (relative per face) for a few small faces.
mattijs 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
Bug in Meshing Parameters menu Spacing1 (1e+10) Karna ANSYS Meshing & Geometry 1 October 12, 2009 14:38
On Bug of Fluent 12.0 lzgwhy FLUENT 0 August 26, 2009 06:41
bug in Rcomponents utility cedric_duprat OpenFOAM Bugs 1 May 7, 2009 02:56
Bug reports Mattijs Janssens (Mattijs) OpenFOAM 0 January 10, 2005 11:05
Forum y2k Bug Jonas Larsson Main CFD Forum 1 January 5, 2000 11:22


All times are GMT -4. The time now is 21:40.