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

VIV of a cylinder at Re=120

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 23, 2011, 22:32
Default VIV of a cylinder at Re=120
  #1
New Member
 
ojha.mayank485's Avatar
 
Mayank Ojha
Join Date: May 2011
Posts: 22
Rep Power: 6
ojha.mayank485 is on a distinguished road
Hello,

Am trying to simulate the displacement of a cylinder in a cross flow at Re=120. I first did a mesh and time convergence study on a stationary cylinder and got quite satisfying results but when I used a moving mesh in CFX 12.0, I am unable to get the desired Amplitude. Can anyone help me out as to what is going wrong ????

1. My timesteps are small enough to capture the vortex shedding phenomena.
2. I am using a higher order scheme.
3. No turbulence involved as flow is laminar.
4. Am using a CEL expression for cylinder displacements.

Has anyone done simulations for low Re flow over cylinder????


Thank you very much,

Mayank Ojha
ojha.mayank485 is offline   Reply With Quote

Old   May 24, 2011, 07:48
Default
  #2
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,823
Rep Power: 85
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
Assuming you have accurately modelled the stationary cylinder then it suggests your implementation of the cylinder motion is either incorrect or not accurate. Have you considered doing this with the various built-in rigid body solvers in CFX V13, such as the immersed solid approach and the rigid body solver?
ghorrocks is offline   Reply With Quote

Old   May 24, 2011, 13:44
Default @ Ghorrocks
  #3
New Member
 
ojha.mayank485's Avatar
 
Mayank Ojha
Join Date: May 2011
Posts: 22
Rep Power: 6
ojha.mayank485 is on a distinguished road
Quote:
Originally Posted by ghorrocks View Post
your implementation of the cylinder motion is either incorrect or not accurate.
Does that mean that the ODE [ mass*acc + damper*vel + spring stiffness*disp=force ]
is being modelled wrong ????
When I tried to give a pre-defined forcing value such as force=sin( t / 1[s]), It gives me the O/p what am expecting. I even used a higher order RK4 scheme for it. But it does not work.

Quote:
Originally Posted by ghorrocks View Post
Have you considered doing this with the various built-in rigid body solvers in CFX V13, such as the immersed solid approach and the rigid body solver?
How will using a different solver matter ??? I am having CFX 12.0 and not 13. It will be some time till the University server gets updated with CFX 13

Thank you for your reply.

P.S: I have been breaking my head over this for a long time now. I need to get done with this ASAP. Any kind of suggestions is highly appreciated.


- Mayank Ojha
ojha.mayank485 is offline   Reply With Quote

Old   May 24, 2011, 19:21
Default
  #4
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,823
Rep Power: 85
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
I said either incorrect or not accurate. If you are sure your implementation of the equation is correct then your approach is not accurate. Numerical accuracy is a very different thing.

I suggest V13 as it has several methods of doing exactly what you are doing built-in. As it appears the implementation you have done is not accurate, the different implementation by ANSYS may be accurate.

If your university has up to date TECS/leases then you are entitled to V13. You have paid for software you have not got around to installing - this sounds like waste to me. You could even install CFX V13 on a machine you have access to yourself just to test whether V13 fixes your problem.
ghorrocks is offline   Reply With Quote

Old   June 10, 2011, 14:21
Default
  #5
Senior Member
 
Ugly Kid Joe
Join Date: Aug 2010
Posts: 193
Rep Power: 6
vmlxb6 is on a distinguished road
Quote:
Originally Posted by ghorrocks View Post

I suggest V13 as it has several methods of doing exactly what you are doing built-in. As it appears the implementation you have done is not accurate, the different implementation by ANSYS may be accurate.
Hey Glen, I have got the V13 installed on my m/c. Now I ran a bunch of case on it using the Rigid Body Solver. I tried from tstep=0.01,0.001 & 0.0001. My vortex shedding time period at re=100 is 0.14 [s]. Even if I take 1% of the tstep its 0.0014. But for all three tsteps I do not match the Displacements nor the freq.

Question: Should I try reducing my tsteps ? I think its already way too low. Going even lower doesn't make any sense.

When I ran a turbulence model (SST) just to check my y+, it was 0.06 which is super low. Is it possible to have problems because of very fine mesh.


Quote:
Originally Posted by ghorrocks View Post
such as the immersed solid approach
Documentation say that immersed body solvers should not be used when the simulation requires accurate boundary layer prediction.

Is there any thing that I can do about this ????? Please HELP!!!!!!!
vmlxb6 is offline   Reply With Quote

Old   June 11, 2011, 07:35
Default
  #6
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,823
Rep Power: 85
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
Quote:
Should I try reducing my tsteps ? I think its already way too low.
Why is it too low? This is a common mistake. If the result is approaching the correct value as you decrease timestep then I would keep going smaller.

Quote:
Is it possible to have problems because of very fine mesh.
Yes, especially when you have a very fine mesh which expands to a much larger mesh. This can be helped by using double precision solver.

If you do not like the not-as-accurate boundary layer approach used in immersed solids then use the moving mesh rigid body solver approach. It will be far slower and if the motion is large you will have to be careful to not fold the mesh, but you can retain a good boundary layer mesh.
ghorrocks is offline   Reply With Quote

Old   June 11, 2011, 20:04
Default @Glen
  #7
Senior Member
 
Ugly Kid Joe
Join Date: Aug 2010
Posts: 193
Rep Power: 6
vmlxb6 is on a distinguished road
Hey Glen,

When the Rigid body is defined and we define the Spring constant and external forces should the gravity term be defined along with it ??? The solver theory guide: Ch-9 Rigid body theory (pg 342), It is written that the eq of motion is written as force equal to mass into acc. [mx"=Faero+mg-Kspring(x-xso) + Fext] here the mg force is considered automatically or should we define the gravity term ???

In my case should I be defining the gravity term ???
vmlxb6 is offline   Reply With Quote

Old   June 12, 2011, 09:09
Default
  #8
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,823
Rep Power: 85
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
If gravity is a significant force on the body then yes, include gravity.
ghorrocks 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
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
how to know the cylinder vibration in viv natto FLUENT 1 May 23, 2011 20:41
Specifying boundary faces failes in blockMesh blaise OpenFOAM Meshing & Mesh Conversion 0 May 10, 2010 03:56
Animating VIV around 2D cylinder Dash Swift FLUENT 0 August 19, 2006 20:19


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