CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

3D vorticity contour

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

LinkBack Thread Tools Display Modes
Old   May 1, 2014, 02:34
Default 3D vorticity contour
New Member
Shahzad Munir
Join Date: Feb 2012
Posts: 14
Rep Power: 5
shahzad_munir is on a distinguished road
Dear All,
I am using following code to calculate vorticity.
vort = zeros(nx,nj,ns); % initialize the matrix for the vorticity
for i = 1:size(u,3)
[dudx,dudy] = gradient(u(:,:,i),h);
[dvdx,dvdy] = gradient(v(:,:,i),h);
vort(:,:,i) = dvdx - dudy;
I can plot it by using contour(x,y,squeeze(vort(:,:,1))) to get a 2D map. As this vort is a 3D matrix this contour will show vorticity plot in z=1 plane and because z=1 to ns (no of files) I will get different contour plots for diff values of z.
I want to know is this right? Because I'm getting different contours of vorticity how will I know which one is truly representing the vorticity in the flow? Or can I use this vort(:,:,i) outside the loop in above routine? In literature authors shows only one vorticity contour map so I am confused a bit. Please guide me to the best.
shahzad_munir is offline   Reply With Quote


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
vorticity contour achyutan OpenFOAM Post-Processing 0 April 29, 2014 08:09
discontinuous region over Vorticity contour jyh3134 FLUENT 0 February 17, 2014 19:37
Vorticity Contour maysmech OpenFOAM 3 November 14, 2013 11:40
Vorticity Contour VSass OpenFOAM Post-Processing 11 July 2, 2013 11:02
vorticity boundary condition bearcharge Main CFD Forum 0 May 14, 2010 11:32

All times are GMT -4. The time now is 01:46.