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/)
-   -   Specify turbulent inlet (http://www.cfd-online.com/Forums/openfoam-solving/60043-specify-turbulent-inlet.html)

chris1980 August 29, 2006 02:49

I want to specify a inlet flow
 
I want to specify a inlet flow with backround turbulence. For this reason I have defined a inlet bc with fixed value for k and epsilon (I am using the k-e epsilon turbulence model). But the following error occurs:

--> FOAM FATAL ERROR : fixedValue is the wrong k patchField type for wall-functions on patch inlet
should be zeroGradient

I know I am using a high-Re turbulence model with wallfunctions and because of that I have to use zeroGradient at the walls.
But in my case I want to define a turbulent inlet.

Does anybody know what is going wrong or what I have to do?

hjasak August 29, 2006 02:55

Yes. You have specified patch
 
Yes. You have specified patch type wall for your inlet, which is not clever. You will find this either in your blockMeshDict or in the boundary file, depending where you got the mesh from.

chris1980 August 29, 2006 03:03

Hrv, I have defined my inle
 
Hrv,

I have defined my inlet patch as type 'patch' not as wall?

chris1980 August 29, 2006 03:09

solved http://www.cfd-online.c
 
solved http://www.cfd-online.com/OpenFOAM_D...part/happy.gif

The problem was caused due an old boundary file in polymesh (blockmesh does not update it?).

hjasak August 29, 2006 03:11

The code says: forAll(p
 
The code says:

forAll(patches, patchi)
{
const fvPatch& curPatch = patches[patchi];

if (isType<wallfvpatch>(curPatch))
{
# include "checkPatchFieldTypes.H"


and in checkPatchFieldTypes.H

if (!isType<zerogradientfvpatchscalarfield>(k_.bounda ryField()[patchi]))
{
FatalErrorIn("wall-function evaluation")
<< k_.boundaryField()[patchi].type()
<< " is the wrong k patchField type for wall-functions on patch "
<< curPatch.name() << nl
<< " should be zeroGradient"
<< exit(FatalError);

}

As you can see, you patch IS a wall.

Hrv


All times are GMT -4. The time now is 10:50.