
[Sponsors] 
August 28, 2013, 10:11 
LES: Qcriterion visualize/recognize turbulence

#1 
Senior Member
Join Date: Feb 2011
Posts: 140
Rep Power: 7 
Dear all,
I am using the LES model for my stenosis geometry since I found this strange behaviour for the velocity magnitude (and also for the subgrid turbulent viscosity ratio) (see pics) a few diamters behind the stenosis and also read some papers dealing with the apperance of vortices behind a stenosis. But when I plot the Q criterion it is just zero for the whole domain. Does this mean there are no turbulences? What is a criterion for the apperance of a "real" turbulence and how can I visualize this? It would be really nice if somebody could give me ahint! Thanks amillion in advance! Lilly 

August 28, 2013, 10:28 

#2 
Senior Member
Alex
Join Date: Jun 2012
Location: Germany
Posts: 1,281
Rep Power: 22 
When you say that the value of the Qcriterion is zero, do you mean that it is REALLY zero or do you mean that its magnitude throughout the domain is very low compared to the maximum that fluent uses for scaling the countour range?
If it is the first option, did you make sure that the qcriterion is among the result file quantities? If it is option 2, rescale the contour plot to an appropriate range. 

August 29, 2013, 03:16 

#3 
Senior Member
Join Date: Feb 2011
Posts: 140
Rep Power: 7 
Thank you so much, flotus1!
It's really zero. And I did the post processing with fluent itself (so I guess the Q criterion should be among the result quantities, or is this wrong?) Is there another way to visualize turbulence? Or might there be no turbulence if Q is zero, even the the vel magnitude and the viscosity looks so? Or am I doing something wrong at plotting the Q criterion? Thank you for any idea! Lilly 

August 29, 2013, 03:52 

#4 
Senior Member
Alex
Join Date: Jun 2012
Location: Germany
Posts: 1,281
Rep Power: 22 
The normalized Qcriterion is not among the standard quantities that fluent writes to the result file. You could try to add it in the "run calculation" tab under "data file quantities.
But since it is a variable derived from the velocity field, it should be available no matter if it is written to the result file or not, at least when you do the postprocessing in fluent directly. This is odd... Some other identifiers for turbulent flow regions may be the vorticity, helicity, criterion and of course the fluctuations of velocity and pressure. The two contour plots you show in your first post clearly indicate that there are regions of turbulent flow. 

August 29, 2013, 08:53 

#5 
Senior Member
Join Date: Feb 2011
Posts: 140
Rep Power: 7 
Thank you again, flotus1!
I was looking at the results of another sim that I am running at the moment and the Q criterion lies between 9e+6 and 9e6 there (is it normal that it is negative by the way?), so I have no idea what happend to results of the first simulation... Thank you for the description of the other identifiers of turbulence! So a plot of them with nonzero values would indicate turbulence? I made a iso surface plot (Iso Q criterion=100) and colored it by vorticity magnitude. Would this wavy structure at the left of the pic would be an inidicator of turbulence there? Thank you again for any idea! Lilly 

August 29, 2013, 10:24 

#6  
Senior Member
Alex
Join Date: Jun 2012
Location: Germany
Posts: 1,281
Rep Power: 22 
Quote:
Quote:
A similar problem seems to be present in the picture you attatched. You will have to play around with the values of the isosurfaces until the plot looks the way you want it to. The most obvious measure for identifying regions of turbulent flow are the fluctuations of velocity and pressure. 

August 30, 2013, 03:46 

#7  
Senior Member
Join Date: Feb 2011
Posts: 140
Rep Power: 7 
Quote:
I should definitely think before writing : of course the range is 9e+6 to 9e+6, so there are really negative values for the Q criterion (?) I also realized why the Q criterion was zero: A really strange thing, I have to visualize the "Vorticity Magnitude" before visualizing the "Q criterion". Then the "Q criterion" is not zero anymore, I have no idea why this happens... So as a measure of turbulence you would show some plots of the the velocity magnitude at the same time of the period for different periods? If they are not congruent, there are fluctuations in within the velocity field and turbulence is present(?) So a simple eddy at a single time as shown in my first pic is not an evidence of turbulent, yet? Thank you again for any idea! It is really helpful to talk to somebody how is familiar with turbulence modelling, thank you so much! Lilly 

August 30, 2013, 03:56 

#8 
Senior Member

Hi,
i just wanted to add little bit of information regarding Q criterion.. please excuse me if you had already discussed about this.. Q criterion defines the difference between the vorticity magnitude and the strain rate. Therefore sometimes it is better to visualise turbulence throuigh Q criterion instead of vorticity or pressure plots. So basically, if the rotation dominates the strain/shear then the Q criterion is positive and you have active vortices in the flow. you can also manually define the Q criterion inside Fluent through custom defenitions, if i remember correctly. 

September 4, 2013, 05:01 

#9 
Senior Member
Join Date: Feb 2011
Posts: 140
Rep Power: 7 
Thank you for your information, Ananthakrishnan!
So would a negative Q criterion mean that there is no turbulence or would a Qcriterion value of zero mean there is no turbulence present? Thank you for any idea! Lilly 

September 4, 2013, 13:30 

#10 
Senior Member

The Q scalar is defined as:
Q=0.5*(dudx^2+dvdy^2+dwdz^2)dudy*dvdxdudz*dwdxdvdz*dwdy Hence, any parallel flow (e.g., laminar flows in ducts) will have Q = 0. In general, Q = (d^2p/dx^2 + d^2p/dy^2 + d^2p/dz^2)/(2*rho) and will be zero whenever the pressure has just a linear spatial variation, is positive for local pressure minima and negative for local pressure maxima. As you can see: 1) If Q is exactly zero, the flow is necessarily laminar 2) If Q is not zero the flow is NOT necessarily turbulent. Laminar examples with Q not 0 are all the flows with a 2D rotating pattern (just du/dy*dv/dx < 0 and constant, e.g., u=k*y, v=k*x with k = 0.5 times the z component of the vorticity) 3) Turbulent flows will exhibit local maxima and minima for Q whose value are strictly dependent on the Re number of the flow What matters for turbulent flow visualizations are the positive Q isosurfaces. The specific value to pick is problem/visual appeal dependent. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Utilities: post average turbulence fields and create turbulence fields for LES  Hanzo  OpenFOAM Running, Solving & CFD  6  June 20, 2016 05:22 
Default inlet turbulence for LES  N Barton  CFX  11  July 24, 2015 00:34 
Turbulence dampening due to magnetic field in LES and RAS  eelcovv  OpenFOAM  0  June 8, 2010 11:35 
Code release: Flow Transition and Turbulence  Chaoqun Liu  Main CFD Forum  0  September 26, 2008 17:15 
caracteristic turbulence time in LES  skippy  FLUENT  0  September 16, 2004 10:51 