CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Vortexshedding Behind Cylinder at Re 12e05

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 11, 2009, 11:48
Default Hello Foamers, Iam trying t
  #1
Member
 
Maruthamuthu Venkatraman
Join Date: Mar 2009
Location: Norway
Posts: 80
Rep Power: 8
maruthamuthu_venkatraman is on a distinguished road
Hello Foamers,

Iam trying to predict the vortex shedding for Cylinder at Re 1.2E+05 using Low Re KE model Launder and sharma .

Observed Results :

1. Drag coefficient is not at all oscillating in a regular pattern even after 20 Shedding periods. Finally it gives me a constant value of 0.2.

2. Lift coeffcient starts to oscillate in the begining and then it becomes constant with a value of 1.5 after 6 vortex shedding periods.

3. Yplus is less than 1 , All turbulent quantities are set to 1e-10 for cylinder walls. Symmetry BC is used on the top and bottom domain.
Outlet is set to constant with 0 for pressure and all others are set to zero gradient.

4. k and e values are 0.35 and 310 at inlet.

5. Grid is fine enough to capture this frequency.


If anybody has Successfully predicted the shedding frequency at Re 1e+05 , Then shed some light upon this problem.

Residual for last time step.

Time = 74.2625

Courant Number mean: 0.0140374 max: 0.288887
DICPCG: Solving for p, Initial residual = 0.00921648, Final residual = 8.75815e-05, No Iterations 55
time step continuity errors : sum local = 6.25117e-11, global = 8.42844e-13, cumulative = -3.73129e-13
DICPCG: Solving for p, Initial residual = 0.00491658, Final residual = 4.72697e-05, No Iterations 57
time step continuity errors : sum local = 3.37288e-11, global = 8.02879e-13, cumulative = 4.2975e-13
DICPCG: Solving for p, Initial residual = 0.00082418, Final residual = 7.74048e-06, No Iterations 60
time step continuity errors : sum local = 5.51504e-12, global = 1.09791e-12, cumulative = 1.52766e-12
DILUPBiCG: Solving for epsilon, Initial residual = 3.74769e-06, Final residual = 3.74769e-06, No Iterations 0
DILUPBiCG: Solving for k, Initial residual = 5.534e-06, Final residual = 5.534e-06, No Iterations 0


Thanks
maruthamuthu_venkatraman is offline   Reply With Quote

Old   February 12, 2009, 01:59
Default hi I havent done simulations
  #2
Member
 
Sachin Kanetkar
Join Date: Mar 2009
Posts: 57
Rep Power: 8
sachin is on a distinguished road
hi
I havent done simulations in this area but a small doubt in this simulations because your Reynolds no. is too high. Does Low Re k-eps model handle it?
can you try the solutions with high Re k-eps model...just a suggestions

Sachin
sachin is offline   Reply With Quote

Old   February 12, 2009, 05:24
Default Hello Sachin,
  #3
Member
 
Maruthamuthu Venkatraman
Join Date: Mar 2009
Location: Norway
Posts: 80
Rep Power: 8
maruthamuthu_venkatraman is on a distinguished road
Hello Sachin,
Thanks for your Comments. ACtually i need to resolove the turbulent boundary layers compeletely with y+ less than 1. I realize that in OPENFOAM all High RE model are working with wall functions.

I dont know how to resolve turb boundary layers up to wall with High Re models in OpenFoam. If you know any tricks let me know...

Regards
maruthamuthu_venkatraman is offline   Reply With Quote

Old   February 12, 2009, 06:44
Default We have at least 3 low-Re turb
  #4
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
We have at least 3 low-Re turbulence models with integration to the wall - please have a look at what is implemented. I also think somebody spoke about this recently in the Forum.

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   February 12, 2009, 08:08
Default Hello Hrvoje, I have gon
  #5
Member
 
Maruthamuthu Venkatraman
Join Date: Mar 2009
Location: Norway
Posts: 80
Rep Power: 8
maruthamuthu_venkatraman is on a distinguished road
Hello Hrvoje,

I have gonethrough that thread how to incorporate the near wall treatment for low Reke models. That example he gave for Low Re ke model, Launder and sharma. I have followed it as it is.

i changed the default zero gradient for the walls and set the fixed values for k and e in the order of 1e-10.

But if you mean That i should read the thread how to impliment the Near wall treatment for High Re flows then I will go through it.

Shall i follow the same method for Realizable ke model as i do for low Re ke model in the BC settings.

i.e. setting a fixed value (very low ) for all turbulence quantities at walls as posted by someone... Clarify !

Thanks
maruthamuthu_venkatraman is offline   Reply With Quote

Old   February 18, 2009, 10:34
Default Hi, I found that second ord
  #6
Member
 
ms
Join Date: Mar 2009
Location: West London
Posts: 44
Rep Power: 8
anothr_acc is on a distinguished road
Hi,

I found that second order closure turbulence models were required for separation to occur at the correct angle for sub-critical Reynolds number flows around Re=33k. I am investigating OpenFOAM at the moment for exactly this application; I used Fluent, I'm afraid, for previous work, with the Reynolds Stress Model, resolving through the boundary layer, with time steps able to give more than 50 increments per shedding cycle.

BR,
Mark.
anothr_acc is offline   Reply With Quote

Old   May 31, 2010, 06:36
Default
  #7
Member
 
Sven Winkler
Join Date: May 2009
Posts: 70
Rep Power: 8
sven is on a distinguished road
I was also working on vortex shedding around a circular cylinder at re=90k. Even with a second order turbulence model (I used the standard k epsilon model and a especially to vortex shedding modified version of this model) i did not succeed in getting the correct angel of separation. Furthermore, I did not get a constant signal of CD and CL (drag and lift coefficients) over time. There was always a lot of noise in these signals. I would appreciate if you have any suggestions to me or if we could exchange our experiences! Have a good one.

Sven
sven is offline   Reply With Quote

Old   June 4, 2010, 03:00
Default
  #8
Senior Member
 
Ivan Flaminio Cozza
Join Date: Mar 2009
Location: Torino, Piemonte, Italia
Posts: 207
Rep Power: 9
ivan_cozza is on a distinguished road
Send a message via MSN to ivan_cozza
Hi Foamers,
I'm doing some CFD on a 2D cylinder at Re = 1.4x10^5, rhoPisoFoam solver, k-Omega SST turbulence model, with wall function for mut (mutSpalartAllmarasWallFunction) and omega, and fixedValue k = 1x10^12.

I0m able to obtain stable CD and CL for very low values of Co (0.2), the separation angle and St frequency are captured with an error of 20% with respect to experiments, and the average CD value are around 0.2 (from the experiments, this must be order of 1).

Can anyone share some info on how to obtain better results?

Thanks!
ivan_cozza 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
3-d mesh gambit (cylinder within a cylinder) Jason ANSYS Meshing & Geometry 19 June 26, 2013 10:00
cylinder in 2d juliana pesma FLUENT 5 July 19, 2012 04:37
small cylinder and big cylinder Eric CFX 1 January 26, 2007 06:38
Help meeee: Wave cylinder (sine cylinder) Ngoctq FLUENT 0 December 29, 2005 21:46
2d cylinder Achilleas Phoenics 0 August 30, 2000 14:30


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