CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   2D Simulation using simpleFoam (Error) (http://www.cfd-online.com/Forums/openfoam-solving/98382-2d-simulation-using-simplefoam-error.html)

skyinventorbt March 9, 2012 09:42

2D Simulation using simpleFoam (Error)
 
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 March 9, 2012 10:47

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 March 12, 2012 01:32

Thank You Martin !!

Now the simulations are going on . . . .:)

Another doubt is whether OpenFOAM does real time simulation ?

GregorS March 12, 2012 04:30

Quote:

Originally Posted by skyinventorbt (Post 348856)
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 March 13, 2012 23:30

2D plane jet simulation Error in simulation
 
1 Attachment(s)
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.

:mad:The velocity field shows discontinuous patches which is unphysical !!:mad:

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

lovecraft22 March 15, 2012 12:39

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

skyinventorbt March 16, 2012 00:34

1 Attachment(s)
About 10,000 iterations !!


All times are GMT -4. The time now is 08:06.