CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > Siemens > STAR-CCM+

How to plot `turbulent intensity' in Starccm+

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 6, 2015, 10:04
Default
  #21
Senior Member
 
Matt
Join Date: Aug 2014
Posts: 947
Rep Power: 18
fluid23 is on a distinguished road
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 star-ccm+. 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.
fluid23 is offline   Reply With Quote

Old   May 6, 2015, 10:16
Default
  #22
New Member
 
Join Date: Apr 2015
Posts: 12
Rep Power: 11
ktytagong is on a distinguished road
Thank you very much
ktytagong is offline   Reply With Quote

Old   May 6, 2015, 12:41
Default
  #23
New Member
 
Join Date: Apr 2015
Posts: 12
Rep Power: 11
ktytagong is on a distinguished road
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
ktytagong is offline   Reply With Quote

Old   May 6, 2015, 12:58
Default
  #24
Senior Member
 
Matt
Join Date: Aug 2014
Posts: 947
Rep Power: 18
fluid23 is on a distinguished road
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.
fluid23 is offline   Reply With Quote

Old   May 6, 2015, 13:37
Default
  #25
Senior Member
 
Matt
Join Date: Aug 2014
Posts: 947
Rep Power: 18
fluid23 is on a distinguished road
Use this:

mag($$Velocity)==0 ? 0:sqrt(2/3*${TurbulentKineticEnergy})/mag($${Velocity})
fluid23 is offline   Reply With Quote

Old   May 6, 2015, 13:40
Default
  #26
New Member
 
Join Date: Apr 2015
Posts: 12
Rep Power: 11
ktytagong is on a distinguished road
I will try and update again tomorrow

Thank you very much
ktytagong is offline   Reply With Quote

Old   July 9, 2015, 08:24
Default Hello
  #27
New Member
 
Deutschland
Join Date: May 2015
Posts: 5
Rep Power: 11
sarvesh34 is on a distinguished road
Hello

Can anybody help how to define the code for

(Outlettemperature-Liquidtemperature)/(InletTemperature-LiquidTemperature)

For field Function named Temperature Distribution.

I am New here..

Thanks in advance
sarvesh34 is offline   Reply With Quote

Old   July 9, 2015, 09:42
Default
  #28
Senior Member
 
Matt
Join Date: Aug 2014
Posts: 947
Rep Power: 18
fluid23 is on a distinguished road
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. Right-click 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.
fluid23 is offline   Reply With Quote

Old   July 9, 2015, 11:46
Default
  #29
New Member
 
Deutschland
Join Date: May 2015
Posts: 5
Rep Power: 11
sarvesh34 is on a distinguished road
Thank You very much

I will try it and post if there are any problems
sarvesh34 is offline   Reply With Quote

Old   October 5, 2015, 12:08
Default
  #30
Bin
Member
 
Chan Hiang Bin
Join Date: Apr 2013
Posts: 44
Rep Power: 13
Bin is on a distinguished road
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?
Bin is offline   Reply With Quote

Old   October 5, 2015, 14:37
Default
  #31
Member
 
Join Date: Nov 2014
Posts: 88
Rep Power: 12
hwsv07 is on a distinguished road
Quote:
Originally Posted by Bin View Post
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?
for the cells right next to the wall, y+ should be less than 5 in order for the viscous sublayer to be sufficiently resolved.

of course you can have very small y+, less what you suggested i.e. <1. however, this would mean excessive computation time.
hwsv07 is offline   Reply With Quote

Old   October 6, 2015, 04:56
Default
  #32
Bin
Member
 
Chan Hiang Bin
Join Date: Apr 2013
Posts: 44
Rep Power: 13
Bin is on a distinguished road
Quote:
Originally Posted by hwsv07 View Post
for the cells right next to the wall, y+ should be less than 5 in order for the viscous sublayer to be sufficiently resolved.

of course you can have very small y+, less what you suggested i.e. <1. however, this would mean excessive computation time.
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?
Bin is offline   Reply With Quote

Old   October 6, 2015, 06:48
Default
  #33
Member
 
Join Date: Nov 2014
Posts: 88
Rep Power: 12
hwsv07 is on a distinguished road
Quote:
Originally Posted by Bin View Post
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?
yes.

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.
hwsv07 is offline   Reply With Quote

Old   October 7, 2015, 09:20
Default
  #34
Bin
Member
 
Chan Hiang Bin
Join Date: Apr 2013
Posts: 44
Rep Power: 13
Bin is on a distinguished road
Quote:
Originally Posted by hwsv07 View Post
yes.

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.
I see. Thank for the explanations.
Bin is offline   Reply With Quote

Old   October 15, 2015, 12:08
Default
  #35
Bin
Member
 
Chan Hiang Bin
Join Date: Apr 2013
Posts: 44
Rep Power: 13
Bin is on a distinguished road
Quote:
Originally Posted by MBdonCFD View Post
Use this:

mag($$Velocity)==0 ? 0:sqrt(2/3*${TurbulentKineticEnergy})/mag($${Velocity})
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?
Bin is offline   Reply With Quote

Reply

Tags
graphic, result, starccm+, turbulence intensity


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[swak4Foam] Foam warnings - related to swak4Foam Salam-H 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


All times are GMT -4. The time now is 14:51.