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.

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??

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!

What if, we use the slip wall instead of translational periodic boundary condition?


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.

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).

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.

Sir,

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

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

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

Yes, you have to link the faces. See button in the attached image.

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.

PS : Since we are Sirs so we will not provide any solution instread we will give you asignments, quizes and then end term exam:D

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" :D :)

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.

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.

What does this shows?

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: http://www.cfd-online.com/Forums/flu...periodics.html

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 !!

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.

Here is the image.