CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   k-epsilon and pisoFoam (https://www.cfd-online.com/Forums/openfoam-solving/111696-k-epsilon-pisofoam.html)

 zaynah04 January 13, 2013 08:34

k-epsilon and pisoFoam

hi..
i have a case where i have sent a 2m/s velocity and observe the trail around a cube.
I have successfully done it in Fluent.

i had as boundary condition the inlet velocity, the outflow and the symmetry for my domain then the wall for the cube. i was in turbulent regime and using K-epsilon model with Simple algorithm.

Now i wish to apply it to Openfoam.
so i chose pisoFoam and Ras model

I have changed my U and P Boundary cond accordingly..
My problem is i do not know how to modify the K-epsilon and fv solution and fv scheme accordingly..:(

Any idea from the experienced user on how to modify the above mentioned specifications??

Regards
zaynah

 maHein January 13, 2013 17:12

As a starting point, have a look at one of the tutorial cases.

You should use wall functions for k (kqRWallFunction), epsilon (epsilonWallFunction) and nut (nutUWallFunction) at the walls.

At the inlet you can use the following two boundary conditions to estimate k and epsilon based on the turbulent intensity and the length scale:

turbulentMixingLengthDissipationRateInlet for epsilon

turbulentIntensityKineticEnergyInlet for k

 zaynah04 January 14, 2013 00:13

here is my modified files but i found this error..

Quote:
 keyword s2 is undefined in dictionary "/home/zaynah/Desktop/cavity/0/epsilon::boundaryField" file: /home/zaynah/Desktop/cavity/0/epsilon::boundaryField from line 26 to line 54. From function dictionary::subDict(const word& keyword) const in file db/dictionary/dictionary.C at line 461. FOAM exiting
my epsilon file
Quote:
 dimensions [0 2 -3 0 0 0 0]; internalField uniform 6.768e-6; boundaryField { wall { type epsilonWallFunction; value uniform 0; } p { type fixedValue; value uniform 6.768e-6; s1 { type symmetryPlane; value uniform 0; } s2 { type symmetryPlane; value uniform 0; } v { type fixedValue; value uniform 6.768e-6; } frontAndBack { type empty; } }

 gregor January 17, 2013 08:23

Hi zaynah,you have a missing curly bracket at the defintion of p

 zaynah04 January 20, 2013 04:26

hi Gregor..
in fact you are right and also i had done a mistake on the K file in Boundary condition..
it now work beautifully..
thanks

 All times are GMT -4. The time now is 01:23.