
[Sponsors] 
How to plot `turbulent intensity' in Starccm+ 

LinkBack  Thread Tools  Search this Thread  Display Modes 
May 6, 2015, 10:04 

#21 
Senior Member
Matt
Join Date: Aug 2014
Posts: 931
Rep Power: 16 
That depends on the wall treatment you are using, the level of accuracy you want and many other factors. The ranges I gave are very broad and loose. If you want to be as accurate as possible use the near wall model and keep your y+ values as close to 1 as you can. Obviously, in most cases it would require cell by cell modification of the mesh to acheive that which you cannot do in starccm+. However, you should be able to keep it fairly constrained in a straight pipe like this.
I suggest you read up on the wall treatment methods to understand the specifics. The help docs have good discussions of each one. 

May 6, 2015, 10:16 

#22 
New Member
Join Date: Apr 2015
Posts: 12
Rep Power: 9 
Thank you very much


May 6, 2015, 12:41 

#23 
New Member
Join Date: Apr 2015
Posts: 12
Rep Power: 9 
I got a problem again.
According to the field function, I did follow what you said. Apparently, I got 'Turbulent intensity field function: Floating point exception [divide by zero]Turbulent intensity field function: Floating point exception [divide by zero]' Then I could not add the function of turbulent intensity 

May 6, 2015, 12:58 

#24 
Senior Member
Matt
Join Date: Aug 2014
Posts: 931
Rep Power: 16 
The function is added, but it is blowing up because of the zero velocity in the boundary layer. That was an oversight on my part. For now you can just make it:
sqrt(2/3*${TurbulentKineticEnergy})/mag($${Velocity}+0.0001) I am looking to see if I can make it a conditional statemtent to set Intensity to zero if velocity is zero. The equation above should be accurate enough though. You will get high TI anywhere you have zero velocity like the BL, but should be fine in the core flow. That 0.0001 will not impact anything out there. 

May 6, 2015, 13:37 

#25 
Senior Member
Matt
Join Date: Aug 2014
Posts: 931
Rep Power: 16 
Use this:
mag($$Velocity)==0 ? 0:sqrt(2/3*${TurbulentKineticEnergy})/mag($${Velocity}) 

May 6, 2015, 13:40 

#26 
New Member
Join Date: Apr 2015
Posts: 12
Rep Power: 9 
I will try and update again tomorrow
Thank you very much 

July 9, 2015, 08:24 
Hello

#27 
New Member
Deutschland
Join Date: May 2015
Posts: 5
Rep Power: 9 
Hello
Can anybody help how to define the code for (OutlettemperatureLiquidtemperature)/(InletTemperatureLiquidTemperature) For field Function named Temperature Distribution. I am New here.. Thanks in advance 

July 9, 2015, 09:42 

#28 
Senior Member
Matt
Join Date: Aug 2014
Posts: 931
Rep Power: 16 
First, I am assuming you want this to be calculated for every cell not averaged for the whole fluid. If that is wrong let me know...
Create 2 Reports: 1. Mass Flow Averaged Temperature at Inlet (boundary only), call it InletT 2. Mass Flow Averaged Temperature at Outlet (boundary only), call it OutletT Create a Field Function: 1. Rightclick tools > field functions and select new > scalar. 2. Rename the field function and the function name in the properties window to 'Temperature Distribution'. 3. Leave dimensions as 'dimensionless' since you have T/T. 4. Set the definition to ($OutletReport  $Temperature)/($InletTReport$Temperature) It should then appear anywhere you can select scalr field functions (scenes, reports, monitors, etc...). You can treat it just like any other variable. 

July 9, 2015, 11:46 

#29 
New Member
Deutschland
Join Date: May 2015
Posts: 5
Rep Power: 9 
Thank You very much
I will try it and post if there are any problems 

October 5, 2015, 12:08 

#30 
Member
Chan Hiang Bin
Join Date: Apr 2013
Posts: 44
Rep Power: 11 
I have a question. I understand that we need to check the y+ value based on the selected wall treatment. Use low y+ wall treatment as example, the desire y+ should be 1 to 5. y+>5 is definitely not recommended, but how about y+<1?


October 5, 2015, 14:37 

#31  
Member
Join Date: Nov 2014
Posts: 88
Rep Power: 10 
Quote:
of course you can have very small y+, less what you suggested i.e. <1. however, this would mean excessive computation time. 

October 6, 2015, 04:56 

#32 
Member
Chan Hiang Bin
Join Date: Apr 2013
Posts: 44
Rep Power: 11 
This means it won't cause anything bad even the y+ is out of the recommended range (i.e. 1 to 5), just the computation time will be increased?


October 6, 2015, 06:48 

#33  
Member
Join Date: Nov 2014
Posts: 88
Rep Power: 10 
Quote:
consider an unsteady simulation. of course, you can have a timestep size that is infinitely small. but what for? you probably can still can get the same result in a reasonable amount of time at a sufficiently large timestep. 

October 7, 2015, 09:20 

#34 
Member
Chan Hiang Bin
Join Date: Apr 2013
Posts: 44
Rep Power: 11 

October 15, 2015, 12:08 

#35 
Member
Chan Hiang Bin
Join Date: Apr 2013
Posts: 44
Rep Power: 11 
Apart of writing this field function, how about monitoring the variance of velocity [i,j,k]? As far as I know, variance is the second moment of the flow field, means it is (u`)^2; (v`)^2; (w`)^2, which can represent the fluctuaing component. Am I right?


Tags 
graphic, result, starccm+, turbulence intensity 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
[swak4Foam] Foam warnings  related to swak4Foam  SalamH  OpenFOAM Community Contributions  20  August 2, 2015 16:40 
Question on Turbulence Intensity  Eric  FLUENT  1  March 7, 2012 05:30 
About Turbulence Intensity (Pipe flow assimilated)  gRomK13  Main CFD Forum  1  July 10, 2009 04:11 
Loss of Contour Plot colour intensity  Bou  Siemens  3  June 13, 2007 21:19 
graph plot  anindya  Main CFD Forum  2  September 17, 2003 13:00 