# Limiting turbulent viscosity

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

 September 19, 2011, 17:24 Limiting turbulent viscosity #1 Member   Jubayer Join Date: Oct 2009 Location: The University of Western Ontario, London, Ontario Posts: 42 Blog Entries: 1 Rep Power: 7 Hi, I want to limit turbulent viscosity but I don't know how to do that in OpenFOAM. Can someone please help me with this ? Thanks. Jubayer

 September 21, 2011, 14:27 #2 Member   Jubayer Join Date: Oct 2009 Location: The University of Western Ontario, London, Ontario Posts: 42 Blog Entries: 1 Rep Power: 7 Hi, I am using pisoFoam and my fvSchemes is as follows: ddtSchemes { default Euler; } gradSchemes { default cellMDLimited Gauss linear 0.5; grad(p) cellMDLimited Gauss linear 0.5; grad(U) cellMDLimited Gauss linear 0.5; // grad(U) cellLimited Gauss linear 1; } divSchemes { default none; div(phi,U) Gauss limitedLinearV 1; div(phi,k) Gauss limitedLinear 1; div(phi,omega) Gauss limitedLinear 1; div((nuEff*dev(grad(U).T()))) Gauss linear limited 0.5; } laplacianSchemes { default none; laplacian(nuEff,U) Gauss linear limited 0.5; laplacian((1|A(U)),p) Gauss linear limited 0.5; laplacian(DkEff,k) Gauss linear limited 0.5; laplacian(DomegaEff,omega) Gauss linear limited 0.5; } interpolationSchemes { default linear; interpolate(U) linear; } snGradSchemes { default limited 0.5; } fluxRequired { default no; p; } What should I do if I want to limit my nu to specific values like 0 to 1e6? Jubayer

 September 22, 2011, 06:20 #3 Senior Member   Vesselin Krastev Join Date: Jan 2010 Location: University of Tor Vergata, Rome Posts: 361 Rep Power: 10 In order to limit the turbulent viscosity you have to modify the source file related to the turbulence model you are using, by adding a limiter in the eddy (turbulent) viscosity calculation formula. Anyway, I don't know what is your application, but you should be careful in adding an arbitrary limiter in a pre-existing turbulence model. V.

 September 22, 2011, 09:54 #4 Member   Jubayer Join Date: Oct 2009 Location: The University of Western Ontario, London, Ontario Posts: 42 Blog Entries: 1 Rep Power: 7 Thanks Vesselin for your reply. Jubayer

 September 29, 2011, 09:54 #5 Member   Jubayer Join Date: Oct 2009 Location: The University of Western Ontario, London, Ontario Posts: 42 Blog Entries: 1 Rep Power: 7 Hi, To bound nut, I have added this line to the LaunderSharmaKE model, bound(nut_, dimensionedScalar("0", nut_.dimensions(), 10.0)); After compiling and everything, at each time step it showed me the bounding values. But this is not actually bounding the values as I can see that the nut values are going above 10. Can someone please explain me how this bounding actually works? Thanks. Jubayer

 April 6, 2015, 14:18 #6 Member   Charlie Join Date: Dec 2010 Location: 415 Kinross Dr. Newark, DE 19711 Posts: 78 Rep Power: 6 This depends on how to calculate the eddy viscosity, usually the unexpected large eddy viscosity is caused by the deviding a relative small value, for example, in k-Epsilon model, nut = C*k^2/epsilon, to limit the eddy viscosity, an effective way is to limit the smallest epsilon value (but not unphysically large). Charlie

 Tags limit, turbulent viscosity

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post cfdiscool FLUENT 10 June 10, 2015 06:15 shib FLUENT 0 June 22, 2010 12:44 nuimlabib Main CFD Forum 0 August 4, 2009 00:05 varghese FLUENT 2 November 15, 2003 09:56 David Yang FLUENT 3 June 3, 2002 06:13

All times are GMT -4. The time now is 07:34.