[Rachit]

Hi everyone!

I have been trying to mesh flow over the busemann's biplane.

http://en.wikipedia.org/wiki/Busemann%27s_Biplane

The blockMeshDict file I have written is as follows:

vertices
(
(-1.0 0.192 -0.01)
(-0.5 0.192 -0.01)
(0.5 0.192 -0.01)
(1.0 0.192 -0.01)
(-1.0 0.096 -0.01)
(-0.5 0.096 -0.01)
(0.5 0.096 -0.01)
(1.0 0.096 -0.01)
(0 0.0842 -0.01)
(0 -0.0842 -0.01)
(-1.0 -0.096 -0.01)
(-0.5 -0.096 -0.01)
(0.5 -0.096 -0.01)
(1.0 -0.096 -0.01)
(-1.0 -0.192 -0.01)
(-0.5 -0.192 -0.01)
(0.5 -0.192 -0.01)
(1.0 -0.192 -0.01)
(-1.0 0.192 0.01)
(-0.5 0.192 0.01)
(0.5 0.192 0.01)
(1.0 0.192 0.01)
(-1.0 0.096 0.01)
(-0.5 0.096 0.01)
(0.5 0.096 0.01)
(1.0 0.096 0.01)
(0 0.0842 0.01)
(0 -0.0842 0.01)
(-1.0 -0.096 0.01)
(-0.5 -0.096 0.01)
(0.5 -0.096 0.01)
(1.0 -0.096 0.01)
(-1.0 -0.192 0.01)
(-0.5 -0.192 0.01)
(0.5 -0.192 0.01)
(1.0 -0.192 0.01)
);

blocks
(
hex (22 23 5 4 18 19 1 0) (200 1 152) simpleGrading (1 1 1)
hex (23 24 6 5 19 20 2 1) (200 1 152) simpleGrading (1 1 1)
hex (24 25 7 6 20 21 3 2) (200 1 152) simpleGrading (1 1 1)
hex (28 29 11 10 22 23 5 4) (200 1 152) simpleGrading (1 1 1)
hex (29 27 9 11 23 26 8 5) (200 1 152) simpleGrading (1 1 1)
hex (27 30 12 9 26 24 6 8) (200 1 152) simpleGrading (1 1 1)
hex (26 8 5 23 24 6 6 24) (200 1 152) simpleGrading (1 1 1)
hex (27 9 11 29 30 12 12 30) (200 1 152) simpleGrading (1 1 1)
hex (30 31 13 12 24 25 7 6) (200 1 152) simpleGrading (1 1 1)
hex (32 33 15 14 28 29 11 10) (200 1 152) simpleGrading (1 1 1)
hex (33 34 16 15 29 30 12 11) (200 1 152) simpleGrading (1 1 1)
hex (34 35 17 16 30 31 13 12) (200 1 152) simpleGrading (1 1 1)
);

edges
(
);

patches
(

);

mergePatchPairs
(
);

However, I get an error of zero or negative pyramid volume at several faces. I do understand that probably there has been mistake on my part in defining the geometry. However, it is a tedious job going through such blockMesh files this long or even longer again and again.

Are there any other ways or techniques to spot such errors such as badly defined block or face, etc?

Thanks in advance,

Rachit

[wyldckat]

Hi Rachit,

If you are using OpenFOAM 2.0 with the official reader for ParaView, not just the internal one, then try the following command when inside your case folder:

paraFoam -block

This blew my mind just now... I thought this only showed the points, but when I saw the results for the tutorial "stressAnalysis/solidDisplacementFoam/plateHole"... See the image attached!! Notice that I picked the option "Wireframe" on that drop-down box on the top center-left of the window.

Best regards,
Bruno

[wiedangel]

thank you for this option ... it is really f a big help for me. makes it easy to detect the errors