|
[Sponsors] |
Accuracy Problem with Flow over 2D cylinder at Transcritical Reynolds Number |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
February 25, 2014, 06:30 |
Accuracy Problem with Flow over 2D cylinder at Transcritical Reynolds Number
|
#1 |
New Member
Join Date: Nov 2012
Posts: 16
Rep Power: 14 |
Hi all,
I am trying to perform what I thought would be a fairly straightforward validation case in which I want to predict the Strouhal Number and mean drag coefficient for a 2D circular cylinder at a Reynolds number of 5,000,000 but am getting nowhere near the results I should be, can anyone help?! I am comparing against data given in a paper by Roshko (1961) found here: http://authors.library.caltech.edu/10105/1/ROSjfm61.pdf From this paper I should be achieving a Cd of approx 0.7 and St of approx 0.27 (I know other authors cite anywhere between 0.25 and 0.27 so I'd be happy with anything like this!). My results are way off the mark so any advice would be great! My setup is: Cylinder diameter: 5m Reynolds number: 5,000,000. Domain extents: upstream 10D, downstream 20D, Upper and Lower 10D. SST turb model (runs both with and without trans turbs gamma theta model) Timestep of approx 1/50 expected shedding frequency, my monitor points look really smooth so I don't think there is any issue here although I will do a time-step sensitivity study when I have checked everything in the setup. Mesh: Fully hexahedral with an average yplus oscillating between 0.6 and 0.9. There are approx 125 streamwise nodes along the cylinder surface. Convergence criteria at 1E-5, each iteration is solving within 4 coefficient loops and the RMS Courant Number is around 4.4. Boundary conditions: Inlet: Normal velocity inlet 15.45m/s (Based on air at 25C properties, isothermal flow and required Reynolds Number of 5,000,000). Outlet: Average static pressure 0Pa (across whole outlet) Upper and Lower: Free slip walls. Cylinder end planes: Symmetry. Cylinder: No-slip smooth wall. Default turbulence intensity of 5%. Results: Without transitional turbulence: Cd=approx 0.4, St=approx 0.32 With trans turbulence (gamma theta): Cd=approx 0.23, St= approx 0.4 Results are independent of mesh but no time-step dependency study has been done, I want to have my setup checked first. Any help would be gratefully received, I know there are other posts of this nature but I cant see anything that helps me past where I already am. I can post the CCL for the most recent run if neccessary, Cheers, Chris |
|
February 25, 2014, 18:27 |
|
#2 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,869
Rep Power: 144 |
I would do a time step size and convergence tolerance sensitivity check, but if you say you are converging to 4 coeff loops per iteration then your current time step looks pretty good.
There is an FAQ on accuracy but you appear to have looked at most of the issues discussed on the page. Still it is worth a look, it might remind you of something you missed: http://www.cfd-online.com/Wiki/Ansys..._inaccurate.3F Please post your CCL and images of you mesh. |
|
February 26, 2014, 08:31 |
|
#3 |
New Member
Join Date: Nov 2012
Posts: 16
Rep Power: 14 |
Hi Glenn,
Thanks for looking at this. I've put a link below to the CCL for the run with no transitional turbulence, can put the other up if necessary. Also in the same place are a few mesh images https://www.dropbox.com/sh/kvxj8ziq37rjqpa/EiQ4VU2Ns6 Any help greatly appreciated Chris |
|
February 26, 2014, 17:56 |
|
#4 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,869
Rep Power: 144 |
* Use adaptive time stepping, or at least do a time step sensitivity check on the time step you are using.
* Check that the inlet turbulence you have specified is not diffusing the vorticies. Try with a very low inlet turbulence. * Your y+ is averaged over the whole cylinder. This includes the front section which has a boundary layer, and the rear section which will have a large separation. Getting the average y+ over both these regions does not appear useful. I would think calculating y+ up until separation (of maybe just the front half) would be useful. * Check whether your separation point is about right. It might give you a clue. * You have first order turbulence numerics. You might need second for this. |
|
February 28, 2014, 18:58 |
|
#5 |
New Member
Join Date: Nov 2012
Posts: 16
Rep Power: 14 |
Thanks for the advice, been a bit distracted by another job but will look. At these points on Monday...............
|
|
July 12, 2016, 18:35 |
|
#6 |
New Member
marko milosevic
Join Date: Jun 2016
Posts: 4
Rep Power: 10 |
hi,
this is an older post, but lets give it a try, I am dealing with similar case, did you solve your problem cheers |
|
July 14, 2016, 06:51 |
|
#7 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,869
Rep Power: 144 |
There is a general FAQ on accuracy. Have you read it? http://www.cfd-online.com/Wiki/Ansys..._inaccurate.3F
|
|
Tags |
2d simulation, cylinder, drag coefficients, high reynolds number, vortex shedding |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Foam::error::PrintStack | almir | OpenFOAM Running, Solving & CFD | 92 | May 21, 2024 08:56 |
Flow past a cylinder, Reynolds No: 100, Tracking strouhal Number. | Karan | FLUENT | 11 | November 14, 2016 05:49 |
Flow past rotating cylinder: Problem with ForeCoeffs | raf1111 | OpenFOAM | 1 | December 16, 2013 10:45 |
Can 'shock waves' occur in viscous fluid flows? | diaw | Main CFD Forum | 104 | February 16, 2006 06:44 |
Difficulties in solving a high Reynolds number Flow? | wowakai | Main CFD Forum | 10 | December 29, 1998 14:46 |