[yousuf]

i have posted boundary conditions, check mesh results and div schemes.........have a look at them

[mkraposhin]

I'm running your case, and results i will obtain tomorrow.
velocity speed 73.9 m/s is for this mesh (convertion factor .0254)?
so, when you was running for convertion factor 1 and 25.4 does you change velocity speed also?

why are you using 'slip' for periodic walls? (not 'cyclic')?

[mkraposhin]

i'm trying turbFoam variant, with div settings for
div(phi,k) and div(phi,epsilon):

div(phi,k) Gauss upwind;
div(phi,epsilon) Gauss upwind;

and it converges (11 iterations), but epsilon is still bounded

because it takes approx 2-3 minutes for each iteration, i will report you about results tomorrow

[mgz1985]

hi all

i am trying to do an incompressible flow over a 2D airfoil. everything seems to be fine. my mesh check is OK.when i run icoFoam i get the following error

This mesh contains patches of type empty but is not 1D or 2D
by virtue of the fact that the number of faces of this
empty patch is not divisible by the number of cells.

i am new to openfoam and this is my first test case after having practiced thru the tutorials. I ahd set the faces in the z-direction to empty as mentioned in the user guide.

can ne1 help?

thanx

[ngj]

Hi

It would mean that you have added cells in the z-direction but still uses empty on the x-y-normal patches, which means that there are no boundary conditions for the z-direction.

/ Niels

P.S. Sometime I have experienced that the boundary file is not updated even though you make a blockMesh. Thus try to delete /case/constant/polyMesh/boundary and see if that helps.

[mgz1985]

hi neil,

i dont understand the meaning of setting the BC for the empty patches. if u do not mind can u elaborate.yeah i have just 1 cell in the z-direction as mentioned in the tutorials. i deleted the boundary file and tried but did not help. same error.

i m putting my U boundaries here.

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
movingWall
{
type fixedValue;
value uniform (-55.55 0 0);
}

fixedWall
{
type fixedValue;
value uniform (0 0 0);
}

frontAndBack
{
type empty;
}
}

[ngj]

Could you please upload your case and I'll take a look at it.

/ Niels

[mgz1985]

convertToMeters 0.001;

vertices
(
(181 -27.39 0)
(209.26 61.17 0)
(600 0.19 0)
(600 500 0)
(209.26 500 0)
(181 -500 0)
(600 -500 0)
(181 -27.39 1)
(209.26 61.17 1)
(600 0.19 1)
(600 500 1)
(209.26 500 1)
(181 -500 1)
(600 -500 1)

(2000 500 0)
(2000 -500 0)
(2000 500 1)
(2000 -500 1)

);

blocks
(
hex (0 1 4 5 7 8 11 12) (60 50 1) simpleGrading (1 1 1)
hex (0 5 6 2 7 12 13 9) (50 30 1) simpleGrading (1 1 1)
hex (1 2 3 4 8 9 10 11) (30 50 1) simpleGrading (1 1 1)
hex (3 6 15 14 10 13 17 16) (60 50 1) simpleGrading (1 1 1)
);

edges
(
polyLine 0 1 (
(151.99 -26.69 0)
(122.49 -25.57 0)
(92.87 -23.67 0)
(63.09 -20.64 0)
(48.1 -18.59 0)
(32.99 -15.91 0)
(17.63 -11.97 0)
(9.72 -8.96 0)
(6.43 -7.24 0)
(4.73 -6.13 0)
(0 0 0)
(1.27 7.88 0)
(2.56 9.7 0)
(5.28 12.72 0)
(12.37 18.52 0)
(27.01 27.09 0)
(41.9 33.62 0)
(56.91 39 0)
(87.13 47.45 0)
(117.51 53.7 0)
(148 58.21 0)
(178.58 61.17 0)
)

polyLine 7 8(
(151.99 -26.69 1)
(122.49 -25.57 1)
(92.87 -23.67 1)
(63.09 -20.64 1)
(48.1 -18.59 1)
(32.99 -15.91 1)
(17.63 -11.97 1)
(9.72 -8.96 1)
(6.43 -7.24 1)
(4.73 -6.13 1)
(0 0 1)
(1.27 7.88 1)
(2.56 9.7 1)
(5.28 12.72 1)
(12.37 18.52 1)
(27.01 27.09 1)
(41.9 33.62 1)
(56.91 39 1)
(87.13 47.45 1)
(117.51 53.7 1)
(148 58.21 1)
(178.58 61.17 1)
)

polyLine 4 5 ((-1000 500 0) (-1000 -500 0))
polyLine 11 12 ((-1000 500 1) (-1000 -500 1))

polyLine 0 2 (
(210.74 -27.35 0)
(239.87 -26.76 0)
(269.09 -25.82 0)
(298.62 -24.49 0)
(328.36 -22.81 0)
(358.25 -20.8 0)
(388.27 -18.51 0)
(408 -16.89 0)
(448.6 -13.44 0)
(490.2 -9.92 0)
(509.15 -8.37 0)
(539.44 -5.79 0)
(569.74 -0.19 0)

)

polyLine 7 9 ( (210.74 -27.35 1)
(239.87 -26.76 1)
(269.09 -25.82 1)
(298.62 -24.49 1)
(328.36 -22.81 1)
(358.25 -20.8 1)
(388.27 -18.51 1)
(408 -16.89 1)
(448.6 -13.44 1)
(490.2 -9.92 1)
(509.15 -8.37 1)
(539.44 -5.79 1)
(569.74 -0.19 1)

)

polyLine 1 2 (
(270.92 60.74 0)
(301.38 57.75 0)
(331.64 53.76 0)
(361.75 48.94 0)
(391.73 43.47 0)
(421.6 37.5 0)
(451.39 31.24 0)
(472 26.86 0)
(481.13 24.92 0)
(510.84 18.6 0)
(512.25 18.3 0)
(516 17.34 0)
(520.2 16.62 0)
(540.56 12.32 0)
(570.25 6.15 0)

)

polyLine 8 9 (
(270.92 60.74 1)
(301.38 57.75 1)
(331.64 53.76 1)
(361.75 48.94 1)
(391.73 43.47 1)
(421.6 37.5 1)
(451.39 31.24 1)
(472 26.86 1)
(481.13 24.92 1)
(510.84 18.6 1)
(512.25 18.3 1)
(516 17.34 1)
(520.2 16.62 1)
(540.56 12.32 1)
(570.25 6.15 1)

)

);

patches
(

wall movingWall
(
(0 1 8 7)
(0 2 9 7)
(1 8 9 2)
)
wall fixedWall
(
(3 6 13 10)
(3 10 16 14)
(6 13 17 15)
(3 4 11 10)
(4 5 12 11)
(5 6 13 12)
)
empty frontAndBack
(
(0 1 4 5)
(0 2 6 5)
(1 2 3 4)
(3 6 15 14)
(10 13 17 16)
(7 8 11 12)
(7 9 13 12)
(8 9 10 11)
)



);

mergePatchPairs
(
);

that is my blockMeshDict

Re is 1.67 million
U is 200kph

thanx a lot neil for ur patience

[ngj]

Please zip the ENTIRE case, as that is what I need.

/ Niels

[mgz1985]

/home/frpelikan/OpenFOAM/mayank/tutorials/icoFoam/check8.tar.gz

[mgz1985]

can u tell me how to paste a zipped folder here. I have never done it b4 and dont c an option to do it

thanx

[ngj]

Look under the /Documentation/Formatting to your left on the Message Board.

[mgz1985]

\attach

[mgz1985]

Mayank_Gupta.zip

[mgz1985]

thanx neils for teaching how to attach in a post.

[mgz1985]

hi neils

any luck . I tried after changing the cell thickness also but same error.

[ngj]

The problem is present, because you have two regions, which are not connected. One of these (or both), has not been defined in the blockMesh file, thus they are given 'empty' as it belong to defaultFaces.

These patches are between the two regions just downstream the airfoil. You will definitely not get the results you anticipate using that mesh. You need to make some connection between the two regions.

/ Niels

[mgz1985]

which patch r u xactly talking bout. I have all the patches attached.

downstream of the airfoil in block 3 i have 1 internal face and the 5 external faces are defined
(3 6 15 14)
(10 13 17 16)
(3 6 13 10) (this is the internal face and i tried without it also
(3 10 16 14)
(6 13 17 15)

(14 15 16 17 ) (the software fails to recognize this face)

all my patches are connected at the internal faces and i do not define them

did the solver run at your system

thanks

[mgz1985]

which patch r u xactly talking bout. I have all the patches attached.

downstream of the airfoil in block 3 i have 1 internal face and the 5 external faces are defined
(3 6 15 14)
(10 13 17 16)
(3 6 13 10) (this is the internal face and i tried without it also
(3 10 16 14)
(6 13 17 15)

(14 15 16 17 ) (the software fails to recognize this face)

all my patches are connected at the internal faces and i do not define them

did the solver run at your system

thanks

[mgz1985]

thanx neils

figured out what u meant. but i have a question, i was working with I-DEAS for meshing. And it allows the meshing. i was trying that as i am new to openfoam

thanx again.