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

Increasing instability for flow around as cylinder.

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 31, 2013, 10:31
Default Increasing instability for flow around as cylinder.
  #1
New Member
 
James Grimshaw
Join Date: Jan 2013
Posts: 8
Rep Power: 4
jgrimshaw is on a distinguished road
I am currently attempting to model flow around a stationary cylinder at a Re=200,000 using Fluent 14. However, after running my simulation it appears my solution is too steady as no vortex shedding occurs.
I am currently using the URANS solver coupled with the K-e turbulance model and have tried running the simulation.
I originally used a 1st order upwind scheme and left the Pressure - Velocity coupling as SIMPLE. I then changed this to a 2nd order upwind scheme and used the PISO to do the pressure - velocity coupling.
I have found my wall y+ value and applied the appropriate wall function to my simulation but in doing so i coarsened my mesh.

Does anybody have any suggestions on how to create instability in the flow and hence produce vortex shedding?
jgrimshaw is offline   Reply With Quote

Old   January 31, 2013, 10:51
Default
  #2
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,099
Rep Power: 16
RodriguezFatz will become famous soon enough
How long have you been waiting? My experience is that you have to wait quite a long time until von Karman shows up...
Is your timestep small enough?
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   January 31, 2013, 11:18
Default
  #3
New Member
 
James Grimshaw
Join Date: Jan 2013
Posts: 8
Rep Power: 4
jgrimshaw is on a distinguished road
My time step was originally 0.1s but I then changed this to 0.05s.
jgrimshaw is offline   Reply With Quote

Old   January 31, 2013, 11:24
Default
  #4
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,099
Rep Power: 16
RodriguezFatz will become famous soon enough
What is the radius of your cylinder and the velocity of your fluid?
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   January 31, 2013, 11:25
Default
  #5
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,654
Rep Power: 23
FMDenaro will become famous soon enough
what about your grid resolution? have you used second order accuracy in time?
check also the value in time of the vertcal component in some points at the rear of the cylinder, along the longitudinal direction, and see if it is really steady for several non-dimensional time unit.
FMDenaro is offline   Reply With Quote

Old   January 31, 2013, 11:45
Default
  #6
New Member
 
James Grimshaw
Join Date: Jan 2013
Posts: 8
Rep Power: 4
jgrimshaw is on a distinguished road
My radius is 0.5m and my flow velocity is 2.364 m/s.
jgrimshaw is offline   Reply With Quote

Old   January 31, 2013, 12:39
Default
  #7
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,099
Rep Power: 16
RodriguezFatz will become famous soon enough
Quote:
Originally Posted by jgrimshaw View Post
My radius is 0.5m and my flow velocity is 2.364 m/s.
Ok, that sounds ok for the time step. Now, could you post a picture of your mesh?

How many timesteps have you been waiting before you stopped?
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   February 11, 2013, 09:28
Default
  #8
New Member
 
James Grimshaw
Join Date: Jan 2013
Posts: 8
Rep Power: 4
jgrimshaw is on a distinguished road
Sorry for the slow reply I have been out of the office for a while. I am currently running 600 time steps.

The mesh i am using is posted below, along the axis it has 96 elements with a bias factor of 10. The average element size was 0.328125m which gives a CFL=0.3602.
jgrimshaw is offline   Reply With Quote

Old   February 11, 2013, 09:32
Default
  #9
New Member
 
James Grimshaw
Join Date: Jan 2013
Posts: 8
Rep Power: 4
jgrimshaw is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
what about your grid resolution? have you used second order accuracy in time?
check also the value in time of the vertcal component in some points at the rear of the cylinder, along the longitudinal direction, and see if it is really steady for several non-dimensional time unit.
The time step accuracy is second order, I have been assuming it is steady as on my Cl vs time plot (otherwise known as convergence history) I have an oscillating graph which quickly dampens wafter 20% of the time has been ran.
jgrimshaw is offline   Reply With Quote

Old   February 11, 2013, 09:33
Default
  #10
New Member
 
James Grimshaw
Join Date: Jan 2013
Posts: 8
Rep Power: 4
jgrimshaw is on a distinguished road
Is there any good fluent literature or troubleshooting websites that either of you could recommend.
jgrimshaw is offline   Reply With Quote

Old   February 11, 2013, 09:38
Default
  #11
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,099
Rep Power: 16
RodriguezFatz will become famous soon enough
Quote:
Originally Posted by jgrimshaw View Post
The mesh i am using is posted below
I dont find it...
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   February 11, 2013, 09:44
Default
  #12
New Member
 
James Grimshaw
Join Date: Jan 2013
Posts: 8
Rep Power: 4
jgrimshaw is on a distinguished road
Sorry it is in the attachment of this post
Attached Images
File Type: jpg Mesh 2.jpg (97.1 KB, 18 views)
jgrimshaw is offline   Reply With Quote

Old   February 11, 2013, 09:49
Default
  #13
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,099
Rep Power: 16
RodriguezFatz will become famous soon enough
What is your y+ currently?
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   February 11, 2013, 10:54
Default
  #14
New Member
 
James Grimshaw
Join Date: Jan 2013
Posts: 8
Rep Power: 4
jgrimshaw is on a distinguished road
I do not have a exact value for the y+ around the cylinder, but the values range between 50 - 850. I have attached the plot below.
Attached Images
File Type: jpg Wall y+.jpg (54.5 KB, 13 views)
jgrimshaw is offline   Reply With Quote

Old   February 11, 2013, 10:59
Default
  #15
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,099
Rep Power: 16
RodriguezFatz will become famous soon enough
That means you are using wall functions. Maybe these wall functions restrain eddy shedding? Did you try to resolve the sheath (y+<=1)?

BTW: I don't think you need a turbulent inlet flow to get a von Kármán street. As I understand it, it is an instability that has in the first instance nothing to do with turbulence.
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   February 11, 2013, 11:48
Default
  #16
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,099
Rep Power: 16
RodriguezFatz will become famous soon enough
Are you using air? I quickly rebuild your case and I am getting a nice vortex street right from the start... Maybe you should post your settings...
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   February 11, 2013, 16:33
Default
  #17
Senior Member
 
Martin Hegedus
Join Date: Feb 2011
Posts: 467
Rep Power: 9
Martin Hegedus is on a distinguished road
You can try giving it just a little kick of angle of attack so the problem is not completely symmetric.

That being said, the instability should develop on it's own even if things were "symmetric" but it will take a while. The time to instability is dependent on how accurate the numerical method is. If it was "perfectly" accurate, the problem will theoretically remain symmetric. In other words, it is statically stable but dynamically unstable. However, I'm not sure what your turbulence model will do. It may dampen out the dynamic instability. This is a reason DES is used.
Martin Hegedus is offline   Reply With Quote

Old   February 11, 2013, 16:40
Default
  #18
Senior Member
 
Martin Hegedus
Join Date: Feb 2011
Posts: 467
Rep Power: 9
Martin Hegedus is on a distinguished road
Oh, one thing to watch is lift. If it starts to oscillate from machine zero you're on your way. But, it sounds like this was zero for you? And, just to be clear, you do have two big vortices (symmetric about x axis) behind your cylinder?
Martin Hegedus is offline   Reply With Quote

Old   February 11, 2013, 16:54
Default
  #19
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 374
Rep Power: 10
arjun is on a distinguished road
i had this issue for long time with lattice boltzmann, every time i tried i got exactly like this. No vortex shedding, untill i put the outer boundaries far away.
arjun is offline   Reply With Quote

Old   February 12, 2013, 03:08
Default
  #20
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,099
Rep Power: 16
RodriguezFatz will become famous soon enough
Quote:
Originally Posted by arjun View Post
i had this issue for long time with lattice boltzmann, every time i tried i got exactly like this. No vortex shedding, untill i put the outer boundaries far away.
But do you think his domain looks like it is too small compared to the radius?
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz 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
flow over a cylinder urgent! kevin FLUENT 8 August 11, 2015 13:00
benchmark: flow over a circular cylinder goodegg Main CFD Forum 12 January 22, 2013 12:47
Incorrect Drag and Drag Coefficient for flow over a cylinder ozzythewise Main CFD Forum 8 June 13, 2012 06:24
Specifying boundary faces failes in blockMesh blaise OpenFOAM Meshing & Mesh Conversion 0 May 10, 2010 03:56
rotating cylinder using sliding mesh Tim Daly FLUENT 1 November 10, 2008 00:02


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