
[Sponsors] 
Transformation of a cartesian velocity gradient volTensorField to cylindrical coord. 

LinkBack  Thread Tools  Display Modes 
June 24, 2009, 16:11 
Transformation of a cartesian velocity gradient volTensorField to cylindrical coord.

#1 
New Member
Alex
Join Date: Jun 2009
Location: Erlangen, Germany
Posts: 7
Rep Power: 8 
Hi Community,
my name is Alex and I'm a PhD student at the University of ErlangenNuremberg. Well, my PhD thesis comprises some numerical simulations, namely the outflow of a nonnewtonian fluid from a cylindroconical vessel. My "problem" now is the following. I need the shear rates that occur in this vessel. As this setup is axisymmetric, the best would be to get the shear rates in cylindrical coordinates. The flow in the cylindrical part (like in a tube) only has a axial velocity component, the angular and radial velocities are zero. Only in the conical part radial components occur. I am mostly interested in d(Uz)/dr, which is the shear rate of interest. Is it possible to transform the cartesian gradient of the velocity into it's corresponding cylindrical counterpart? Does a simple coordinate transformation the job? If yes, how do I do that? I am not a professional programmer, the only thing that helps me is my basic java knowledge. I searched the forum, and the only things I can find are not useful. I would really appreciate your help and I am looking forward to your answers!!! Best regards Alex 

June 25, 2009, 11:19 
Supplement

#2 
New Member
Alex
Join Date: Jun 2009
Location: Erlangen, Germany
Posts: 7
Rep Power: 8 
Hi again,
perhaps I should concretize my question, so it is easier for you to answer. So what i need is the shear rate in this vessel with the outflowing nonnewtonian fluid. Mainly I am interested in d(Uz)/dr which is the L31 of the gradient tensor of the velocity field in cylindrical coordinates. My simulation yields the velocity in cartesian coordinates. For me the most reasonable way would be to transform the velocity in cartesian to cylindrical coordinates and than take the appropriate definition of the gradient of the velocity vector (meaning in cylindrical coordinates). My concern is, if I take the transformed velocity vector (which is one thing I need to know how to do it, or better to say if it is already implemented) and then define a volTensorField L with grad(U) than the code uses the definition of the gradient of U in cartesian coordinates. Or is this fear without any reason because I can tell the code that this vector is in cylindrical coordinates and it takes the right derivation on its own? After solving this topic I will store the field values in a file to view it with paraView, which is already working (with the L31 entry in cartesian coordinates). Than I want to write a script which extrakts the mean shear stress, the standard deviation, the minimum and maximum value and perhaps the median. Are there already functions implemented that calculate the median and the standard deviation? The others I found already, and for the standard deviation it shouldn't be a big problem to compute this, but if it is already implemented... Now I hope you get my point. It would be very kind and nice if you are able to help me! Thanks very much in advance and my best regards!! Alex 

April 25, 2012, 14:23 

#3 
New Member
Caroline Vandame
Join Date: Aug 2010
Posts: 26
Rep Power: 6 
I know this thread is quite old, but any chance you found how to get the velocity vector components in a cylindrical coordinate system from the cartesian coordinates?


April 25, 2012, 18:03 

#4 
Senior Member

What about using the calculator in paraview to do the math?


April 26, 2012, 18:20 

#5 
New Member
Caroline Vandame
Join Date: Aug 2010
Posts: 26
Rep Power: 6 
I'm running OpenFoam on a virtual machine so I'm trying to avoid Paraview. But thanks for the tip, I'll look into it


April 26, 2012, 18:29 

#6 
Senior Member

I'm running it on a virtual machine as well. I created a shared folder so that once the case has run I can access the data using paraview from the host system.
Just so that you know, sometimes, working on a shared folder, can led to some errors due to missing files (that are not missing at all). It0s just a matter of unmounting and remounting the shared folder. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Velocity gradient tensor  desA  Main CFD Forum  36  January 3, 2007 01:31 
FLUENT  velocity gradient  Sue  Main CFD Forum  1  May 8, 2006 15:09 
Calculation of velocity gradient  Chris  Main CFD Forum  8  December 18, 2003 09:35 
Velocity gradient  Stephen  Main CFD Forum  0  April 7, 2003 10:28 
Terrible Mistake In Fluid Dynamics History  Abhi  Main CFD Forum  12  July 8, 2002 09:11 