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

FLuent simulation of taylor couette flow of concentric cylinder geometry.

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

Like Tree30Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   October 28, 2012, 10:37
Default
  #21
New Member
 
rshbhb
Join Date: Oct 2012
Posts: 26
Rep Power: 4
rshbhb is on a distinguished road
Dear Daniele Sir,

Thank you for the information. I was a bit busy in the previous two days. I'll simulate the system again and will let you know the results or any doubts that I may inccur.

I would like to Thank Mr.Far for joining the conversation and sharing his results and doubts with us.
rshbhb is offline   Reply With Quote

Old   October 28, 2012, 11:50
Default
  #22
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,910
Blog Entries: 6
Rep Power: 38
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
Quote:
Originally Posted by rshbhb View Post
Dear Daniele Sir,

Thank you for the information. I was a bit busy in the previous two days. I'll simulate the system again and will let you know the results or any doubts that I may inccur.

I would like to Thank Mr.Far for joining the conversation and sharing his results and doubts with us.
It was really informative topic and I especially thankful to rshbhb and Daniele for sharing info and subsequent simulations by Daniele. If needed I can also upload ICEM meshes and for theses meshes y+ is less than one as required by SST K-W or standard K-W.
ghost82 likes this.
Far is offline   Reply With Quote

Old   October 28, 2012, 12:00
Default
  #23
Senior Member
 
ghost82's Avatar
 
Daniele
Join Date: Oct 2010
Location: Italy
Posts: 899
Rep Power: 15
ghost82 will become famous soon enough
Quote:
Originally Posted by Far View Post
So what is the Reynolds number for u1= 5 rad/sec, 20 rad/sec ?
In my opinion, the Reynolds number, which is dimensionless is:

Re=(U0*d*density)/mu

U0 cannot be expressed in rad/s, otherwise the Re is not dimensionless, so it must be in m/s, so I think U0 is the peripheral speed of the rotating cylinder.
U0=omega*radius, where omega is the angular velocity in rad/s

d is the characteristic length, in my case=(outer radius - inner radius) and it is expressed in m

density is the density of the fluid (kg/m3)

mu is the dynamic viscosity of the fluid, expressed in Pa*s or, in other words, kg/(s*m)

With the dimensional analysis:

Re=(m/s)*(m)*(kg/m3)*(s*m/kg)

Re is dimensionless.
This is how I calculate the Re

Referring to your data and your cited article I think Re, to be dimensionless is calculated as:
Re=(density*(angular velocity)*radius*radius)/mu

where (angular velocity)*radius is the peripheral velocity; the only difference is the characteristic length

density in kg/m3
angular velocity in rad/s
radius in m (of the rotating cylinder)
mu in Pa*s

So for u1=5 rad/s, Re (your article)=50
for u2=20 rad/s, Re (your article)=200

PS: I think in page 4 of your cited article there's an error: Re is calculated with radius=1 m, but it is not rotating!

PS2: density=1 kg/m3: what type of fluid is it?is it a gas?Is the Taylor-Couette flow valid for gases??
Far likes this.

Last edited by ghost82; October 28, 2012 at 12:31.
ghost82 is offline   Reply With Quote

Old   October 28, 2012, 12:02
Default
  #24
Senior Member
 
ghost82's Avatar
 
Daniele
Join Date: Oct 2010
Location: Italy
Posts: 899
Rep Power: 15
ghost82 will become famous soon enough
Quote:
Originally Posted by Far View Post
If needed I can also upload ICEM meshes and for theses meshes y+ is less than one as required by SST K-W or standard K-W.
Even if I don't use ICEM and I don't know how to use it (till now..) I think it can be of great help!
Far likes this.
ghost82 is offline   Reply With Quote

Old   October 28, 2012, 12:16
Default
  #25
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,910
Blog Entries: 6
Rep Power: 38
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
What if, we use the slip wall instead of translational periodic boundary condition?


Quote:
Originally Posted by ghost82 View Post
Even if I don't use ICEM and I don't know how to use it (till now..) I think it can be of great help!
Full project which i have done uptill now including the research material:
https://dl.dropbox.com/u/68746918/Co...nstability.rar

If only meshes are needed then click on this link :
https://dl.dropbox.com/u/68746918/me...ute_taylor.rar

PS: Dropbox will take some time to upload the files. May be few hours.
ghost82 and rshbhb like this.
Far is offline   Reply With Quote

Old   October 28, 2012, 12:30
Default
  #26
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,910
Blog Entries: 6
Rep Power: 38
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
Quote:
Originally Posted by ghost82 View Post

PS2: density=1 kg/m3: what type of fluid is it?is it a gas?does the Taylor-Couette flow is valid for gase??
I think in the end what matters is the Reynolds number! It does not matter as I believe and it is the essence of non dimensional quantities.
ghost82 likes this.
Far is offline   Reply With Quote

Old   October 28, 2012, 12:37
Default
  #27
New Member
 
rshbhb
Join Date: Oct 2012
Posts: 26
Rep Power: 4
rshbhb is on a distinguished road
Quote:
Originally Posted by Far View Post
I think in the end what matters is the Reynolds number! It does not matter as I believe and it is the essence of non dimensional quantities.
Dear Sirs,

I want to ask you which fluent version are you using?

I am using GAMBIT for geometry and FLuent Solver (Fluent 6.3.26). I want to know that will i get such good images with the versoin of Fluent that I have?

Sir I'll upload the geometry (in GAMBIT) in 15 mins or so could you please verify if it'll be ok for correct simulation. MY Laptop is not very fast (only 3GB RAM).
rshbhb is offline   Reply With Quote

Old   October 28, 2012, 12:44
Default
  #28
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,910
Blog Entries: 6
Rep Power: 38
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
Quote:
Originally Posted by rshbhb View Post
Dear Sirs,

I want to ask you which fluent version are you using?

I am using GAMBIT for geometry and FLuent Solver (Fluent 6.3.26). I want to know that will i get such good images with the versoin of Fluent that I have?

Sir I'll upload the geometry (in GAMBIT) in 15 mins or so could you please verify if it'll be ok for correct simulation. MY Laptop is not very fast (only 3GB RAM).
I am using Ansys Fluent 14.

Yes you can get the same images in Fluent, there is no change in these features. They have just improved flow schemes and included other turbulnece model and wall treatments.

I have also 4GB laptop. But I will be more than happy to look at your mesh. Meanwhile you can compare your mesh with mine or Daniele's. Daniele's opinion will be more valuable in this regard.
Far is offline   Reply With Quote

Old   October 28, 2012, 12:58
Default
  #29
New Member
 
rshbhb
Join Date: Oct 2012
Posts: 26
Rep Power: 4
rshbhb is on a distinguished road
Quote:
Originally Posted by Far View Post
I am using Ansys Fluent 14.

Yes you can get the same images in Fluent, there is no change in these features. They have just improved flow schemes and included other turbulnece model and wall treatments.

I have also 4GB laptop. But I will be more than happy to look at your mesh. Meanwhile you can compare your mesh with mine or Daniele's. Daniele's opinion will be more valuable in this regard.
Dear Sir,

For faster communication i have sent you chat request on your gmail id.
Please accept it.
rshbhb is offline   Reply With Quote

Old   October 28, 2012, 13:11
Default
  #30
New Member
 
rshbhb
Join Date: Oct 2012
Posts: 26
Rep Power: 4
rshbhb is on a distinguished road
Sir,

How do I hard link the faces before meshing in GAMBIT?

I am just getting the option of linking the face meshes.
rshbhb is offline   Reply With Quote

Old   October 29, 2012, 03:24
Default
  #31
Senior Member
 
ghost82's Avatar
 
Daniele
Join Date: Oct 2010
Location: Italy
Posts: 899
Rep Power: 15
ghost82 will become famous soon enough
Quote:
Originally Posted by Far View Post
What if, we use the slip wall instead of translational periodic boundary condition?
The Taylor-Couette flow still exists; in this simulation I cut the domain to create a symmetry plane, so the z coordinate is 0.06 m.
Grid is the same as the previous simulation.
Dimensions are the same as reported in some posts above, Re is 1000 (considering the caracteristic length R2-R1).
I had more troubles to get convergence with this case.

I also have a laptop with win7 64 bit with only 4 gb ram and a workstation with win7 64 bit and 32 gb ram; on both systems ansys 14 is installed.
For these simulations I used the laptop and it took about 1 day to reach a drop in residuals of 1*10^-7.
I think you will have no problem using fluent 6.x version.

Daniele
Attached Images
File Type: jpg taylor60-1-wall.jpg (20.9 KB, 35 views)
File Type: jpg taylor60-2-wall.jpg (60.5 KB, 42 views)
File Type: jpg taylor60-3-wall.jpg (44.0 KB, 39 views)
Far likes this.

Last edited by ghost82; October 29, 2012 at 03:41.
ghost82 is offline   Reply With Quote

Old   October 29, 2012, 03:27
Default
  #32
Senior Member
 
ghost82's Avatar
 
Daniele
Join Date: Oct 2010
Location: Italy
Posts: 899
Rep Power: 15
ghost82 will become famous soon enough
Quote:
Originally Posted by rshbhb View Post
Sir,

How do I hard link the faces before meshing in GAMBIT?

I am just getting the option of linking the face meshes.

Yes, you have to link the faces. See button in the attached image.
Attached Images
File Type: png link.png (1.7 KB, 144 views)
ghost82 is offline   Reply With Quote

Old   October 29, 2012, 04:07
Default
  #33
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,910
Blog Entries: 6
Rep Power: 38
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
Quote:
Originally Posted by ghost82 View Post
I think you will have no problem using fluent 6.x version.

Daniele
But you must be careful with Y+ for old version. As in new version hybrid and scalable wall functions are implemented which are tolerate to users and meshes they create. For older version you must ensure Y+ > 30 for wall function and Y+<1 for integration to wall requirements.

In fluent K-epsilon and K-omega both models are equipped with both : wall function and integration to wall. It should be noted that K-epsilon models are by defaults wall function models but modified (coupled with single equation model in boundary layer) to account for the near wall effects while K-omega models are originally designed for intergration to wall approach.

So be more careful with Y+ for older versions.

Quote:
Dear Sirs,
PS : Since we are Sirs so we will not provide any solution instread we will give you asignments, quizes and then end term exam
ghost82 and rshbhb like this.
Far is offline   Reply With Quote

Old   October 29, 2012, 12:17
Default
  #34
New Member
 
rshbhb
Join Date: Oct 2012
Posts: 26
Rep Power: 4
rshbhb is on a distinguished road
I have now got somewhere close to getting the result. Here is the result of my simulation done yesterday.

I think I've not done the iterations for sufficient no. of times.( Right?)

I had used k-omega SST model for this.

Angular Velcity 6 rad/s.

Also please could you explain what do you mean by 'Y+' (is the residual y values you are talking about or something else and how can I control the Y+ values)

Thank you "SIRs"
Attached Images
File Type: jpg Case 2.JPG (51.7 KB, 36 views)
File Type: jpg Plot of residuals.JPG (31.0 KB, 26 views)
rshbhb is offline   Reply With Quote

Old   October 29, 2012, 12:27
Default
  #35
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,910
Blog Entries: 6
Rep Power: 38
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
Quote:
Originally Posted by rshbhb View Post
I have now got somewhere close to getting the result. Here is the result of my simulation done yesterday.

I think I've not done the iterations for sufficient no. of times.( Right?)

I had used k-omega SST model for this.

Angular Velcity 6 rad/s.

Also please could you explain what do you mean by 'Y+' (is the residual y values you are talking about or something else and how can I control the Y+ values)

Thank you "SIRs"
Yes this is one of source of errors : non converged solution.

Y+ is non dimensional parameters which is due to friction velocity, first mesh point distance from wall and kinematic viscosity. http://www.cfd-online.com/Wiki/Dimen...tance_(y_plus)

You must ensure two things in mesh:

1. First cell height (Expansion ratio should be less than 1.2)
2. Total height of boundary layer. you can get estimate from this link (for flat plate but good enough for initial estimate)
http://www.see.ed.ac.uk/~johnc/teach...dboundary.html
Also try this link: http://www.ae.metu.edu.tr/~ae244/doc...07/node12.html

The mesh, which is shared in post, has y= < 1, so it is good for SST-KW model.
rshbhb likes this.
Far is offline   Reply With Quote

Old   October 29, 2012, 13:09
Default
  #36
Senior Member
 
ghost82's Avatar
 
Daniele
Join Date: Oct 2010
Location: Italy
Posts: 899
Rep Power: 15
ghost82 will become famous soon enough
Quote:
Originally Posted by rshbhb View Post
I have now got somewhere close to getting the result. Here is the result of my simulation done yesterday.

I think I've not done the iterations for sufficient no. of times.( Right?)

I had used k-omega SST model for this.

Angular Velcity 6 rad/s.

Also please could you explain what do you mean by 'Y+' (is the residual y values you are talking about or something else and how can I control the Y+ values)

Thank you "SIRs"
Good, you're on the right way; let the residual go down;
to see if your solution is converged you can monitor some important value during iterations, for example:

integral on the volume of the turbulent kinetic energy; or/and
Area weight surface integral of velocity magnitude on different planes (for example at some planes at z=cost)
etc.

When your monitor value doesn't change during successive iterations your solution is converged.
See my picture in posts above; as you can see in my example convergence is reached for residuals<1*10^-7!!
Fixing only minimum values for residuals is not a good idea to evaluate if solution is converged; you should monitor some other values!

For the y+ subject Far is very expert as you can read on some other threads on this forum.

PS:
you can plot the y+ value; in plot postprocessing you will find "yplus" under "turbulence"; plot yplus values selecting wall surfaces and see if these values respect the limits.
rshbhb likes this.
ghost82 is offline   Reply With Quote

Old   October 29, 2012, 13:11
Default
  #37
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,910
Blog Entries: 6
Rep Power: 38
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
Quote:
integral on the volume of the turbulent kinetic energy
What does this shows?
Far is offline   Reply With Quote

Old   October 29, 2012, 13:13
Default
  #38
New Member
 
rshbhb
Join Date: Oct 2012
Posts: 26
Rep Power: 4
rshbhb is on a distinguished road
In the mesh you all had made, did you use (or when shall i use) the reverse orientation option while hard linking the mesh (in the 'link mesh faces' option).

Do I use it

for 1. the top and bottom periodic faces, OR
2. the vertical periodic side faces, OR
3. should I use this option at all.

Daniele Sir had provided me with this link: Hard-linking faces in Gambit (periodics)

in which there is a mention that we have to take the direction of the arrows same for both the faces but it not coming out to be same whatever I do. (Is it wrongly written?).
Here is the picture for better understanding.


THanks !!
rshbhb is offline   Reply With Quote

Old   October 29, 2012, 13:19
Default
  #39
Senior Member
 
ghost82's Avatar
 
Daniele
Join Date: Oct 2010
Location: Italy
Posts: 899
Rep Power: 15
ghost82 will become famous soon enough
Quote:
Originally Posted by Far View Post
What does this shows?
It computes the integral of kinetic energy on the domain; monitoring some important values is very important to evaluate convergence; in fact in the mixing tank tutorial it is reported:

The default convergence criteria are not sufficient to get the correct flow features in a mixing tank. To judge the convergence, some of the integrated quantities needs to be monitored along with velocity magnitude around the turbine.
In this problem, monitor the velocity magnitude on a surface just above and below the turbine. Also monitor the volume integral of kinetic energy in the tank.
ghost82 is offline   Reply With Quote

Old   October 29, 2012, 13:23
Default
  #40
New Member
 
rshbhb
Join Date: Oct 2012
Posts: 26
Rep Power: 4
rshbhb is on a distinguished road
Here is the image.
Attached Images
File Type: jpg PRoblem Description while meshing.JPG (21.9 KB, 20 views)
rshbhb 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
Simulation of cylinder aligned in the flow direction using k-e model Eliasjal FLUENT 2 June 13, 2012 08:00
Simulation of steam (CO2 and Water vapor mixture) flow through nozzle using Fluent. Jimmy FLUENT 0 March 2, 2011 13:30
Simulation Flow Around cylinder 3D Jwolf CFX 19 November 25, 2009 15:21
Compressible flow simulation using FLUENT arun Main CFD Forum 0 February 16, 2004 16:44
Simulation of the Flow past a circular cylinder using STAR-CD M. S. GUEROUACHE Main CFD Forum 0 October 1, 1998 10:51


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