[skyinventorbt]

I am trying to simulate a plane jet problem. Basically started from pitzDaily tutorial.
I have created mesh using blockMesh and the blockMeshDict file is as below...

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 0.001;

vertices
(
(-1 0 0) //0
(0 0 0) //1
(0 0 0.1) //2
(-1 0 0.1) //3
(-1 0.5 0) //4
(0 0.5 0) //5
(0 0.5 0.1) //6
(-1 0.5 0.1) //7
(250 0 0) //8
(250 0 0.1) //9
(250 0.5 0) //10
(250 0.5 0.1) //11
(0 200 0) //12
(250 200 0) //13
(250 200 0.1) //14
(0 200 0.1) //15
);

blocks
(
hex (0 1 5 4 3 2 6 7) (10 20 1) simpleGrading (1 1 1)
hex (1 8 10 5 2 9 11 6) (500 20 1) simpleGrading (10 1 1)
hex (5 10 13 12 6 11 14 15) (500 400 1) simpleGrading (10 40 1)
);

edges
(
);

patches
(
patch inlet
(
(0 4 7 3)
)
patch outlet
(
(8 10 11 9)
(10 13 14 11)
)
wall nozzleup
(
(4 5 6 7)
)
wall domainleft
(
(5 12 15 6)
)
symmetryPlane sym
(
(0 1 2 3)
(1 8 9 2)
)
empty frontAndBack
(
(0 1 5 4)
(1 8 10 5)
(5 10 13 12)
(3 2 6 7)
(2 9 11 6)
(6 11 14 15)
)
);

mergePatchPairs
(
);

// ************************************************** *********************** //

when i try to solve this using simpleFoam I get this following error.


$ simpleFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.7.1 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/


// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting RAS turbulence model kEpsilon


--> FOAM FATAL 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.

From function emptyFvPatchField<Type>::updateCoeffs()
in file fields/fvPatchFields/constraint/empty/emptyFvPatchField.C at line 150.

FOAM exiting


"KINDLY SUGGEST ME TO SOLVE THIS ERROR"

[MartinB]

Hi,

you forgot one side of your domain in your patch definition. Therefor it becomes a "defaultFaces" patch with "empty" condition assigned.

You can see the missing side in paraFoam easily by importing the Mesh Parts one after another.

Good luck

Martin

[skyinventorbt]

Thank You Martin !!

Now the simulations are going on . . . .

Another doubt is whether OpenFOAM does real time simulation ?

[GregorS]

Originally Posted by skyinventorbt

Another doubt is whether OpenFOAM does real time simulation ?

If you use simpleFoam, the answer is no. Because every time step in simpleFoam is really just one iteration. So it doesn't matter whether you set
time step to one second or one milisecond.

But if you use other solvers, like pimpleFoam or pisoFoam, I think the answer is yes.

[skyinventorbt]

I just submitted the 2d plane jet problem which i mentioned above in the first post. The velocity is 10 m/s. Kindly see the snapshot of velocity field.

The velocity field shows discontinuous patches which is unphysical !!

What is the procedure to get physically possible solution for this 2d plane jet problem ??

[lovecraft22]

Can you post a picture of the mesh in that region please? Also, how many iterations have you run?

[skyinventorbt]

About 10,000 iterations !!