
[Sponsors] 
Not converging flow around cylinder with Re=1000000. Help please! 

LinkBack  Thread Tools  Display Modes 
December 3, 2010, 10:51 
Not converging flow around cylinder with Re=1000000. Help please!

#1 
New Member
Witold
Join Date: Nov 2010
Posts: 13
Rep Power: 7 
Dear all,
I am doing my Diplomathesis with ANSYS. I need to simulate a flow around a cylinder as in the Karman vortex street. It is a 3d simulation at a reynolds number of 10^6, which I know is very high. To validate my mesh, I made a smoth mesh with N_x1=42, N_x2=980, N_y=42, N_z=12 and N_diag=42 nodes. The geometry has a hight of 20mm a depth of 10mm and a length of 1000mm. The diameter of the cylinder is 5mm. Please have a look at this image. The mediumfine mesh is the fine one scaled with a factor of 0.5. So is the rough mesh 0.5 times the mediumfine. I have huge problmes because the simulation will not converge. Only the rough net with a "Reynolds Stress" turbulence model converged in few timesteps. The residuals of the finer meshes jitter around at about 10^(2) for hours. Because I do not know, what is relevant for deeper explanation, I linked the result and the outputfile here CLICK. I already tried to modify the "Time Scale" to see, if there is an amortisation by changing it. I made the simulations with the "SST" model (Zylinder_3d_re1000000_stat_g_001), changed it for once in "Reynolds Stress" (Zylinder_3d_re1000000_stat_g_004). The file001 did not converged but had a minimal MAX Residual of 10^(4). File004 converged imidiatly (the cause, it goes a long way straight before, is because I had an "maximum Timescale" activeted. When I deactivated it, it went right down to 10^(6). I would appreciate it very much, if someone would have a look at this. This case must be very simple, so there must be something trivial I do not see! Thank you very much. If there are more files or explanations needed, please ask! Greetings Witold 

December 3, 2010, 21:43 

#2 
Member
Darren Leong
Join Date: Dec 2010
Posts: 63
Rep Power: 6 
Hi Witold,
I tried opening your .res files to no avail. Got a few questions: 1. What version of the ANSYS/fluent are you using? 2. What license do you have; research or teaching? (mesh no. limitation applies here) 3. Simulation in steadystate or transient? < has to be transient XD Suggestion: 1. Grid independence study  you might be able to get away with ~500,000elements, let me know how many you're at. 2. I highly suggest you look into defining your first node height by calculating the y+ at <2 for sst model. You should get your vortex street but it may be overpredicted compared to ke model. Put your growth in range of 1 1.1. The mesh doesn't look fine enough as it is to capture the flow characteristics 3. Reduce your timescale. Have a readup on "Courant Number". 4. Solve in double precision (64bit, more information for better accuracy) 5. Extend/enlarge your fine mesh domain, doesn't look big enough as it is for Re=10^6 (optional) 6. Do a literature review/google. Your problem's been already addressed, there are even tutorials available on it =) dazza 

December 4, 2010, 18:27 

#3 
New Member
Witold
Join Date: Nov 2010
Posts: 13
Rep Power: 7 
Hi dazza,
I do not understand, why the file can't be opened. I am using Ansys 12.1 with ICEM CFX for meshing. All with a research Licence (the one, with a limitation of nodes). As you mentioned I think my biggest fault is, that I didn't make the simulation transient. Of course there is no steady state in this flow. 1: I tried to make an independence study, but since my simulations did not converge, there are noch results. I am limited to 512k nodes, wehn I remember right, so this would be my finest mesh. Then go down by a scale of 0.5 each time. 2: with the formula from the Ansys help I realized a y+ about 1 at the cylinder. Because of the change of the shear around the cylinder y+ changes a bit. I have a growth rate of 1.5 at the wall...I will try a ratio of 1.1 . 3: the Courantnumber was a nice hint. It showed, that my timescale was far to long. 4: double precision is always activated. 5: I had this thought too...but how long? Is there any formula like for the laminar or turbulent flow in straight pipes? 6: I will make some more searches. But most of the paper consider high Renumbers up to 1000. But I continue looking. Thank you very much as far. I will make some simulations and post the results here. Greetings Witold 

December 4, 2010, 20:00 

#4 
Member
Darren Leong
Join Date: Dec 2010
Posts: 63
Rep Power: 6 
Hi Witold,
Glad to know it was of help. Best to shift your direction to transient eventually =) The present setup was alright as a preliaminary result. See if you can plot the wake, it'll give you an idea how far to extend domain. The mesh no. limitation will be a problem depending on whether you have enough to visuallise the vortices smoothly (too coarse > pixelated ). Careful about Re at 10^6. Although it's dimentionless, the reported value differs whether they took the diameter or length as ref. I think for your case, you might barely notice the vortex street as it will be compacted at Re of 1000000. Best wishes, dazza ps. try growth rate of 1.3, 1.1 will likely exceed ur mesh limit, 

December 7, 2010, 14:41 

#5 
New Member
Witold
Join Date: Nov 2010
Posts: 13
Rep Power: 7 
Hi!
My Re ist based on the diameter. I know the vortices will be very compacted. I need to have a look on the flow in it's mixed state. I have a further question. I am not an expert in traniscient flow yet. I have timesteps of 7e8s to achieve a Courant number of about 0.01 . When i want to simulate a 8s flow, this would mean to make 114e6 steps, each with 2 till 5 coefficient loops. With 30s per iteration it will need 108 years There must be a faster way, to achieve the results of a fully developed flow, or am I wrong? I haven't googled this yet, I will immediately, but when someone has an answer right now, I would appreciate it. Thank you! Greets Edit: I already initialized with full velocity 

December 7, 2010, 22:07 

#6 
Member
Darren Leong
Join Date: Dec 2010
Posts: 63
Rep Power: 6 
You'll need to run your simulation in parallel if ur pc has multiple processors or over a cluster. However, i think your timestep's too small (overdefined). Try relaxing the value enough for your solution to just converge.
If you thesis work is on a research basis and not coursework, i would suggest investigating/validating your setup at lower speeds to determine the critical reynolds number first and then proceed to Re 10^6. Forum member, Altomos has done similiar work to yours. Might be helpful msging him on details of his setup. Refer to Help! Emergency!! Karman Vortex using ANSYS He used Strouhal number instead to determine the timestep based on the frequency calculated. Edit: Just to be sure, check your user license prefences for ANSYS. If there's a research license in the list move it to the top (change the value to 1 if 0). I think only the teaching license has the 512k limitation. Last edited by Darren Leong; December 8, 2010 at 00:41. 

Thread Tools  
Display Modes  


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 
flow around a cylinder  pXYZ  Main CFD Forum  14  July 25, 2011 10:05 
Simulation Flow Around cylinder 3D  Jwolf  CFX  19  November 25, 2009 15:21 
Anybody knows smth abt the flow around a cylinder?  Wenxuan  Main CFD Forum  3  March 20, 2007 17:19 
fluid flow fundas  ram  Main CFD Forum  5  June 17, 2000 21:31 