|
[Sponsors] |
How to plot `turbulent intensity' in Starccm+ |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
May 4, 2015, 07:49 |
How to plot `turbulent intensity' in Starccm+
|
#1 |
New Member
Join Date: Apr 2015
Posts: 12
Rep Power: 11 |
I would like to plot graph in Starr cmm+ such as velocity magnitude, axial velocity and turbulent intensity.
I have created the plan section in 'derived part'. Then when I plotted it, there several points (lines) in the graph. What should I do if I would like to get only 1 line in the graph. Last edited by ktytagong; May 21, 2015 at 11:15. |
|
May 4, 2015, 11:48 |
|
#2 |
Member
Join Date: Nov 2014
Posts: 88
Rep Power: 12 |
sort the data - go to the properties | data under Plots and sort the data by their Position (X) to plot a nice line.
for turbulent intensity, you need a time series data velocity vs time. with the data, use excel or matlab to compute the mean and standard deviation to find turbulence intensity (std/mean). |
|
May 4, 2015, 12:03 |
|
#3 |
Senior Member
Matt
Join Date: Aug 2014
Posts: 947
Rep Power: 18 |
Your plots are not visible. Turbulence intensity in an analysis that uses a turbulence model (k-e, k-w, sst, etc...) will be a field function. No need to export and calculate.
Not sure exactly what you are trying to do. You have a derived part (plane) but are plotting points in the plane? Please elaborate. |
|
May 4, 2015, 12:49 |
|
#4 | |
Member
Join Date: Nov 2014
Posts: 88
Rep Power: 12 |
Quote:
selecting the RANS equations with the desired turbulence models does give you some field function relating to turbulence (i see Tke, Viscosity Ratio etc). But it does not report on the turbulence intensity. I think one has to export the data and post-process. |
||
May 4, 2015, 13:01 |
|
#5 |
Senior Member
Matt
Join Date: Aug 2014
Posts: 947
Rep Power: 18 |
I stand (partially) corrected. I assumed intensity was an output, not sure why they don't include it. However, export is definitely not needed.
First, if doing a steady RANS analysis you aren't going to see any turbulent fluctuation a given point in your flow. The assumption of RANS is that turbulence all averages out to zero. To get turbulence intensity simply setup a field function: 1. Right-Click Tools > Field Functions 2. Select New > Scalar 3. Set Function Name to Turbulence_Intensity. 4. Leave as dimensionless. 5. Set Definition to sqrt(2/3*${TurbulentKineticEnergy})/mag($${Velocity}) 7. Rename Field Function to Turbulence Intensity. 8. Now use it like you would any other field function. http://www.cfd-online.com/Wiki/Turbulence_intensity |
|
May 4, 2015, 14:48 |
|
#6 | |
Member
Join Date: Nov 2014
Posts: 88
Rep Power: 12 |
Quote:
But Turbulent Intensity by definition, requires a time series. from the time series, one can find the fluctuation and the mean value and thereby compute TI. Like what you said, if it is a steady case, from the above method, one can get the TI - but it was based on the TKE right at the end of the simulation. It is a steady case because some selected parameter (e.g. velocity) has converged. |
||
May 4, 2015, 14:50 |
|
#7 |
Member
Join Date: Nov 2014
Posts: 88
Rep Power: 12 |
But I have some weird thing going on now. I tried to plot the TKE field scene.
I have a steady RANS case which has converged solution - I can see that I have a vorticity field from my vorticity scene and I was expecting to have an obvious TKE field scene. However, my TKE field scene is showing zero throughout the domain. Any idea what could be wrong? |
|
May 4, 2015, 16:33 |
|
#8 |
Senior Member
Matt
Join Date: Aug 2014
Posts: 947
Rep Power: 18 |
I shouldn't have said steady. turbulent fluctuations will always average to zero in a RANS analysis. You may still see transient behavior but not random. Not sure why you say that uses TKE is from the end of the simulation. Definitely not the case.
When you say TKE is zero, are you sure its not just really really small? If this is the same problem you were working back when I helped before, you probably won't have terribly huge values of TKE. Can you send me a shot or something? |
|
May 5, 2015, 05:10 |
|
#9 |
New Member
Join Date: Apr 2015
Posts: 12
Rep Power: 11 |
Sorry for the late reply.
To obtain the graph in 1 lline, I already figured out what I should do. But for the turbulent intensity, I still do not get your points. Obviously, my simulation is in the turbulent flow and I also selected the method of K-e and bla bla any related to turbulent. In fluent, we can get the intensity graph directly ( as I remember) but in starccm+, there are turbulent kinetic energy, turbulent viscosity , and bla bla And the point is I want to get the turbulent intensity. Do you have any idea what should I do? |
|
May 5, 2015, 10:07 |
|
#10 |
Senior Member
Matt
Join Date: Aug 2014
Posts: 947
Rep Power: 18 |
Still no images. To get a field function (aka variable) that you can plot do exactly as I listed. Its not a defalut field function apparently so you must define it yourself. It's not a big deal. If these step-by-steps aren't clear enough I am not sure I can help you.
"To get turbulence intensity simply setup a field function: 1. Right-Click Tools > Field Functions 2. Select New > Scalar 3. Set Function Name to Turbulence_Intensity. 4. Leave as dimensionless. 5. Set Definition to sqrt(2/3*${TurbulentKineticEnergy})/mag($${Velocity}) 7. Rename Field Function to Turbulence Intensity. 8. Now use it like you would any other field function." |
|
May 5, 2015, 10:29 |
|
#11 |
New Member
Join Date: Apr 2015
Posts: 12
Rep Power: 11 |
I will try to follow your step.
And now I have a problem with the "a request for memory failed this is typically due to insufficient virtual memory" when I am simulating. I am monitoring the memory during the simulation, and the swap free memory is showing '0'. The physical properties of my computer are 7 GB of swap and I am using 64 bit of star cmm+ . Do you have any idea to solve this problem?? Moreover, I am simulating a horizontal pipe (simple simulation) with 12 m length and 0.08 m diameter. I set 0.08 for base size, 40 for number of prism layer . Thank you for your kind |
|
May 5, 2015, 10:58 |
|
#12 |
Senior Member
Matt
Join Date: Aug 2014
Posts: 947
Rep Power: 18 |
That is probably an indication that your mesh is exceeding the capacity of your machine. 7GB isn't a whole lot and will limit you to ~<3.5M poly cells.
40 prism layers is very aggressive too. Sounds like you should reassess your approach for pipe flow. Star-CCM has some cylinder mesh options which will use less memory for a straight pipe. |
|
May 5, 2015, 11:11 |
|
#13 |
New Member
Join Date: Apr 2015
Posts: 12
Rep Power: 11 |
I should reduce the number of prism layer to get the lower number of meshes.
The option of cylinder mesher, I have already done with this one. The result is quite weird. Then I am trying to use another functions to get the comparison result. Anyways, Thank you in advance for helping me |
|
May 5, 2015, 11:14 |
|
#14 |
Senior Member
Matt
Join Date: Aug 2014
Posts: 947
Rep Power: 18 |
Can you post a shot of your mesh or something? It's hard to say, 40 layer may be appropriate but it seems high. The cylinder mesher can be troublesome on occasion. It usualy just means you have something set wrong. Also, make sure you mesh don't mesh in parallel when you have that mesh option activated. It can cause issues.
|
|
May 5, 2015, 11:20 |
|
#15 |
New Member
Join Date: Apr 2015
Posts: 12
Rep Power: 11 |
[/IMG]
I am not sure that you can see the picture or not. |
|
May 5, 2015, 11:21 |
|
#16 |
New Member
Join Date: Apr 2015
Posts: 12
Rep Power: 11 |
delete pic
Last edited by ktytagong; May 7, 2015 at 08:29. Reason: wrong link |
|
May 5, 2015, 17:52 |
|
#17 |
Senior Member
Matt
Join Date: Aug 2014
Posts: 947
Rep Power: 18 |
This mesh is most likely inappropriate for your problem. Have you solved anything yet or are you still trying to mesh? Two rules of thumb for meshing:
1. Target wall y+ values that are between 1-5 (for near wall treatment) 30-60 (for high wall treatment) or keep between 1 and 60 for all wall y+ treatment. (Actually, for all wall y+ avoid 5-30 range if you can. The values are interpolated between near wall and high wall methods in this range). Should be easy enough for a simple pipe flow. 2. Set number of prism layers and prism layer ratio such that your prism cells more-or-less transition into your core cells without a large change in volume. Right now it looks like you jump about 400% between your last prism and the core mesh. Keep that closer to 20%. That being said, for your model I would maybe start with 5-7 prism layers. A thickness that is about 1/3 - 1/2 of what you have now and a prism growth rate of 1.3. Mesh, run and look at your y+ values then re-asses your prism layer settings. |
|
May 5, 2015, 17:56 |
|
#18 |
New Member
Join Date: Apr 2015
Posts: 12
Rep Power: 11 |
No, I am still trying to work with the mesh because when I do the simulation, I always get the memory error.
Tomorrow, I will follow your advise and keep update again. Thank you |
|
May 5, 2015, 17:59 |
|
#19 |
Senior Member
Matt
Join Date: Aug 2014
Posts: 947
Rep Power: 18 |
Memory error is because your mesh is too aggressive. Try the settings I mentioned, that should get you moving at least.
|
|
May 6, 2015, 09:50 |
|
#20 |
New Member
Join Date: Apr 2015
Posts: 12
Rep Power: 11 |
I did follow what you have said.
Everything works well (Simulation can run). One more question how do we know what is y+ desired value? Thank you very much |
|
Tags |
graphic, result, starccm+, turbulence intensity |
|
|
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 |