
[Sponsors] 
February 25, 2010, 11:16 
Nonlinear k epsilon Shih Model

#1 
Senior Member
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 10 
Hallo Foamers!
Does anybody know a reference where I can find the derivation of the Nonlinear k e Shih model? Cheers 

May 21, 2010, 02:39 

#2 
Senior Member
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 10 
Does really nobody knows a reference, at least one? HELP


May 21, 2010, 03:05 

#3 
Senior Member
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 10 
Does really nobody knows a reference, at least one? HELP


August 24, 2010, 03:28 

#4 
New Member
Olli
Join Date: Mar 2010
Location: Berlin
Posts: 13
Rep Power: 8 
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 

August 25, 2010, 09:29 

#5 
Member
Join Date: Apr 2009
Location: Karlsruhe, Germany
Posts: 96
Rep Power: 9 
Hi,
have you simulated successful using this turbulencemodel? I have problems near walls (U isn't zero orthogonal to the wall in the nearwallpart even in channelflows). Regards Thomas 

August 26, 2010, 04:29 

#6 
New Member
Olli
Join Date: Mar 2010
Location: Berlin
Posts: 13
Rep Power: 8 
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 

August 26, 2010, 05:11 

#7 
Member
Join Date: Apr 2009
Location: Karlsruhe, Germany
Posts: 96
Rep Power: 9 
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 meshdependence. 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 fvSolutiondata, please? Thanks and regards, Thomas 

August 27, 2010, 06:32 

#8 
New Member
Olli
Join Date: Mar 2010
Location: Berlin
Posts: 13
Rep Power: 8 
Hi Thomas,
a y+ of less than 0.4 seems to be a sufficient fine mesh. Nevertheless, are kepsilon 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 meshdependence you should probably do a grid study :( An other way would be to use an much coarser mesh, if boundary layers are not priorranking. 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((1A(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 1e06; relTol 0.01; } U { solver GAMG; preconditioner DILU; tolerance 1e05; relTol 0.1; mergeLevels 1; smoother GaussSeidel; agglomerator faceAreaPair; nCellsInCoarsestLevel 100; tolerance 1e05; relTol 0; } k { solver PBiCG; preconditioner DILU; tolerance 1e05; relTol 0.1; } epsilon { solver PBiCG; preconditioner DILU; tolerance 1e05; relTol 0.1; } R { solver PBiCG; preconditioner DILU; tolerance 1e05; relTol 0.1; } nuTilda { solver PBiCG; preconditioner DILU; tolerance 1e05; 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; } 

August 27, 2010, 11:05 

#9 
Member
Join Date: Apr 2009
Location: Karlsruhe, Germany
Posts: 96
Rep Power: 9 
Hi,
thank you very much. I have simulated using your given condition, but the velocities didn't get better. Have you ever simulated a nonswirling flow using this model? Regards Thomas 

August 30, 2010, 09:37 

#10 
New Member
Olli
Join Date: Mar 2010
Location: Berlin
Posts: 13
Rep Power: 8 
Hi,
no, I never simulated a nonswirling 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 

Thread Tools  
Display Modes  


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 (subgridscale 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 kw model and SST model  Li  CFX  7  June 29, 2007 04:19 
KEpsilon Model  sangit  Main CFD Forum  2  September 9, 2004 13:19 