CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

Nonlinear k epsilon Shih Model

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   February 25, 2010, 11:16
Default Nonlinear k epsilon Shih Model
  #1
Senior Member
 
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 9
idrama is on a distinguished road
Hallo Foamers!

Does anybody know a reference where I can find the derivation of the Nonlinear k e Shih model?

Cheers
idrama is offline   Reply With Quote

Old   May 21, 2010, 02:39
Default
  #2
Senior Member
 
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 9
idrama is on a distinguished road
Does really nobody knows a reference, at least one? HELP
idrama is offline   Reply With Quote

Old   May 21, 2010, 03:05
Default
  #3
Senior Member
 
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 9
idrama is on a distinguished road
Does really nobody knows a reference, at least one? HELP
idrama is offline   Reply With Quote

Old   August 24, 2010, 03:28
Default
  #4
New Member
 
Olli
Join Date: Mar 2010
Location: Berlin
Posts: 12
Rep Power: 7
Moelfus is on a distinguished road
Hi Claus,

the paper you are searching for can be found at the NASA pages. There you have access to a lot of other papers of Shih:

http://ntrs.nasa.gov/search.jsp?R=35...%257c1%26N%3D0

Best regards,

Olli
Moelfus is offline   Reply With Quote

Old   August 25, 2010, 09:29
Default
  #5
Member
 
Join Date: Apr 2009
Location: Karlsruhe, Germany
Posts: 96
Rep Power: 8
Thomas Baumann is on a distinguished road
Hi,

have you simulated successful using this turbulence-model?
I have problems near walls (U isn't zero orthogonal to the wall in the nearwall-part even in channel-flows).

Regards Thomas
Thomas Baumann is offline   Reply With Quote

Old   August 26, 2010, 04:29
Default
  #6
New Member
 
Olli
Join Date: Mar 2010
Location: Berlin
Posts: 12
Rep Power: 7
Moelfus is on a distinguished road
Hi,

yes I'm using this turbulence model for simulating swirling flows. For this purpose it seems to be perfect. I haven't discovered any problems until know.

Olli
Moelfus is offline   Reply With Quote

Old   August 26, 2010, 05:11
Default
  #7
Member
 
Join Date: Apr 2009
Location: Karlsruhe, Germany
Posts: 96
Rep Power: 8
Thomas Baumann is on a distinguished road
Hi,
I have velocities out of the wall.
I have described this in a posting here in the forum.

Nonlinear turbulence model NonlinearShih

But if I'm using a very fine mesh (y+ < 0.4) the results are better. So the solution has a mesh-dependence.
Maybe it's caused by the fvschemss are the solvers for p,U,k...
I have tried different ones without success.

Could you maybe show your fvschemes and the fvSolution-data, please?

Thanks and regards,

Thomas
Thomas Baumann is offline   Reply With Quote

Old   August 27, 2010, 06:32
Default
  #8
New Member
 
Olli
Join Date: Mar 2010
Location: Berlin
Posts: 12
Rep Power: 7
Moelfus is on a distinguished road
Hi Thomas,

a y+ of less than 0.4 seems to be a sufficient fine mesh. Nevertheless, are k-epsilon models known for their bad boundary layer behaviour. Therefore this model uses an boundary layer model in the near wall region. I can look up how it is set up. If the results have a mesh-dependence you should probably do a grid study :-( An other way would be to use an much coarser mesh, if boundary layers are not prior-ranking.

Olli

fvSchemes:
-----------------------------------------------------------------

ddtSchemes
{
default steadyState;
}

gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
grad(U) Gauss linear;
// grad(U) cellLimited Gauss linear 1;
}

divSchemes
{
default none;
div(phi,U) Gauss linearUpwindV Gauss linear;
div(phi,k) Gauss upwind;
div(phi,omega) Gauss upwind;
div((nuEff*dev(grad(U).T()))) Gauss linear;
div(nonlinearStress) Gauss linear;
div(phi,epsilon) Gauss upwind;
div(phi,R) Gauss upwind;
div(R) Gauss linear;
div(phi,nuTilda) Gauss upwind;
}

laplacianSchemes
{
default Gauss linear corrected;
// default Gauss linear limited 0.5;
// default Gauss linear limited 0.333;
laplacian(nuEff,U) Gauss linear corrected;
laplacian((1|A(U)),p) Gauss linear corrected;
laplacian(DkEff,k) Gauss linear corrected;
laplacian(DepsilonEff,epsilon) Gauss linear corrected;
laplacian(DREff,R) Gauss linear corrected;
laplacian(DnuTildaEff,nuTilda) Gauss linear corrected;

}

interpolationSchemes
{
default linear;
interpolate(U) linear;
}

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p;
}

fvSolution:
-------------------------------------------------------

solvers
{
p
{
solver GAMG;
preconditioner FDIC;
mergeLevels 1;
smoother GaussSeidel;
agglomerator faceAreaPair;
nCellsInCoarsestLevel 100;
tolerance 1e-06;
relTol 0.01;
}

U
{
solver GAMG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
mergeLevels 1;
smoother GaussSeidel;
agglomerator faceAreaPair;
nCellsInCoarsestLevel 100;
tolerance 1e-05;
relTol 0;
}

k
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
}

epsilon
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
}

R
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
}

nuTilda
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
}
}

SIMPLE
{
nNonOrthogonalCorrectors 0;
}

relaxationFactors
{
p 0.2;
U 0.5;
k 0.7;
epsilon 0.7;
R 0.7;
nuTilda 0.7;
}
Moelfus is offline   Reply With Quote

Old   August 27, 2010, 11:05
Default
  #9
Member
 
Join Date: Apr 2009
Location: Karlsruhe, Germany
Posts: 96
Rep Power: 8
Thomas Baumann is on a distinguished road
Hi,

thank you very much.

I have simulated using your given condition, but the velocities didn't get better.
Have you ever simulated a non-swirling flow using this model?

Regards Thomas
Thomas Baumann is offline   Reply With Quote

Old   August 30, 2010, 09:37
Default
  #10
New Member
 
Olli
Join Date: Mar 2010
Location: Berlin
Posts: 12
Rep Power: 7
Moelfus is on a distinguished road
Hi,

no, I never simulated a non-swirling case with this model. I just used the model because all other models either crashed or over predicted the axial velocity. I looked at your other post and really have no idea what is going wrong there.

Olli
Moelfus is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
multiphaseInterFoam for RAS turbulence model chiven OpenFOAM Bugs 7 August 23, 2011 02:52
help for different between les model (subgrid-scale model) liuyuxuan FLUENT 1 October 2, 2009 15:25
2 stage axial turbine model convergence issues sherifkadry CFX 2 September 7, 2009 20:51
simulation results for k-w model and SST model Li CFX 7 June 29, 2007 04:19
K-Epsilon Model sangit Main CFD Forum 2 September 9, 2004 13:19


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