CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

Q-criterion

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

Like Tree1Likes
  • 1 Post By sbaffini

Reply
 
LinkBack Thread Tools Display Modes
Old   April 29, 2013, 07:32
Default Q-criterion
  #1
New Member
 
Carlo
Join Date: Apr 2013
Posts: 2
Rep Power: 0
imaginor is on a distinguished road
Hi,
could you help me to create a Matlab script for visualization of vortex by Q-criterion? I have the gradient of velocity gained by DNS with software Channelflow. This file because of the construct for becomes 4-D double:
32x33x32x9, where 9 indicates le nine components of the gradient in every nodes.
How can I write the script to display the vortex by the Q-criterion?

Q
= 0.5(Ω^2− S^2),

where

Ω=tr[ΩΩt ]^1/2and S=tr[SSt ]^1/2

with

S
= 1/2 (v + (v)t) and Ω= 1/2(v - (v)t)

but how I use v ?

In Matlab, please.

Thanks a lot

imaginor is offline   Reply With Quote

Old   April 29, 2013, 12:16
Default
  #2
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 534
Blog Entries: 14
Rep Power: 17
sbaffini will become famous soon enough
Let us suppose you have the following matrix:

x(nx,ny,nz): x coordinates of your grid

y(nx,ny,nz): y coordinates of your grid

z(nx,ny,nz): z coordinates of your grid

dudx(nx,ny,nz): x-wise derivative of x-wise velocity component

dudy(nx,ny,nz): y-wise derivative of x-wise velocity component
.
.
.
dwdz(nx,ny,nz): z-wise derivative of z-wise velocity component

This is, more or less, how i would put it down in Matlab:

%MATLAB CODE

iso_q=100; %Pick your number here

%Definition of Q
q=-0.5*(dudx.^2+dvdy.^2+dwdz.^2)-dudy.*dvdx-dudz.*dwdx-dvdz.*dwdy;

%Plotting a Q isosurface, Q=iso_q
figure()
p=patch(isosurface(x,y,z,q,iso_q));
set(p,'FaceColor','red','EdgeColor','none');
daspect([1,1,1])
axis tight
ax = -1; ay = 1; az = 1;
view([ax,ay,az]);
camroll(240)
camlight
lighting gouraud

%END OF MATLAB CODE

hope it works
Far likes this.
sbaffini is offline   Reply With Quote

Old   April 30, 2013, 03:19
Default
  #3
New Member
 
Carlo
Join Date: Apr 2013
Posts: 2
Rep Power: 0
imaginor is on a distinguished road
Thanks a lot...
imaginor is offline   Reply With Quote

Old   April 30, 2013, 09:04
Default
  #4
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,948
Blog Entries: 6
Rep Power: 39
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
sbaffini

Can you please write a blog with all this cool stuff (and also for LES)
Far is offline   Reply With Quote

Old   April 30, 2013, 14:32
Default
  #5
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 534
Blog Entries: 14
Rep Power: 17
sbaffini will become famous soon enough
Dear Far,

i'm glad you appreciate my contributions. I admit that is now a long time i'm thinking about putting all this stuff together in some place that i already identified as a cfd-online blog...

However, as you may imagine, time is never enough as all the material would require some modifications or testing before being directly posted for the general use. In contrast, posting specific answers once in a while is more easy and directly useful for the one who asked.

Nonetheless, as my working group was also considering this opportunity to "share stuff", i promise i will now seriously consider posting on a blog.
sbaffini is offline   Reply With Quote

Old   July 1, 2013, 11:31
Default
  #6
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 534
Blog Entries: 14
Rep Power: 17
sbaffini will become famous soon enough
It took me a while to actually understand what you meant by blog. I'm doing it now
sbaffini is offline   Reply With Quote

Old   July 1, 2013, 13:12
Default
  #7
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,948
Blog Entries: 6
Rep Power: 39
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
I am still trying to understand
Far is offline   Reply With Quote

Old   July 1, 2013, 13:21
Default
  #8
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,948
Blog Entries: 6
Rep Power: 39
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
@sbaffini: You have wonderful skills
Far is offline   Reply With Quote

Old   July 1, 2013, 14:14
Default
  #9
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 534
Blog Entries: 14
Rep Power: 17
sbaffini will become famous soon enough
Still, they are quite difficult to exchange with money

Thanks
sbaffini is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
[ICEM] low edge criterion value problem mraxx ANSYS Meshing & Geometry 6 April 30, 2013 09:03
problem with Min/max rho tH3f0rC3 OpenFOAM 7 February 23, 2013 06:37
Asymptotic stopping criterion is not working rob50904 STAR-CCM+ 2 January 28, 2013 09:25
convergence criterion Dominique FLUENT 5 November 24, 2006 03:36
Convergence criterion Moose Main CFD Forum 5 June 9, 2005 20:39


All times are GMT -4. The time now is 22:29.