|
[Sponsors] |
August 27, 2012, 07:14 |
Vortex Shedding on LES
|
#1 |
Senior Member
Julien
Join Date: Jun 2012
Location: France
Posts: 152
Rep Power: 13 |
Hi everybody !
(in fact, this is the continuation of another thread, but with new data. Thus, a new thread ?) This WE was devoted to the calculation of the wake of a simple rectangular cylinder in an homogeneous flow. My purpose is to "tune" OpenFOAM in order to have a correct behaviour concerning the Von Karman Alley, or Vortex Shedding. Here are some characteristics of my run:
A snap shot (arbitrary phase): And zoom near the cylinder: Any advice or comment concerning this run ? The solution seems OK, but the Strouhal number is completely wrong: I found about 0,18, whereas the Eurocode says 0,06 (three time less!!!). What is that ? What did I make wrong to have so huge error? |
|
August 27, 2012, 09:12 |
|
#2 |
Senior Member
A_R
Join Date: Jun 2009
Posts: 122
Rep Power: 16 |
Can I ask you what kind of mesh do you use to create these mesh?
|
|
August 27, 2012, 11:47 |
|
#3 |
Senior Member
Julien
Join Date: Jun 2012
Location: France
Posts: 152
Rep Power: 13 |
blockMesh with very small expanding , and then snappyHexMesh. Does it look strange in your opinion?
|
|
August 28, 2012, 15:19 |
|
#4 |
Senior Member
A_R
Join Date: Jun 2009
Posts: 122
Rep Power: 16 |
the triangle cells may affect your solution
|
|
August 29, 2012, 04:50 |
|
#5 |
Senior Member
Julien
Join Date: Jun 2012
Location: France
Posts: 152
Rep Power: 13 |
I agree. And most of all, some of these triangulars seem in the incorrect way. How to prevent snappyHexMesh from doing this during refinning?
Do you think a bad mesh could explain the error in the Strouhal number? Don't I mysunderstand strongly something else? The wake is 2~3 times faster that it should be! There is a strong bias ! |
|
August 29, 2012, 05:09 |
|
#6 |
Senior Member
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 21 |
It is not a snappyHexMesh issue but a Paraview issue. Deselect 'decompose polyhedra' when loading your foamfile.
|
|
August 29, 2012, 09:45 |
|
#7 |
Senior Member
Julien
Join Date: Jun 2012
Location: France
Posts: 152
Rep Power: 13 |
Thanks Bernhard. So, these triangular have nothing to do with the solution, and thus nothing to do with the tremendous bias on the Strouhal number ...
Nobody has allready modelled correctly the Von Karman Alley with OpenFOAM ? |
|
August 29, 2012, 12:03 |
|
#8 |
Senior Member
Dr. Fabian Schlegel
Join Date: Apr 2009
Location: Dresden, Germany
Posts: 222
Rep Power: 17 |
I never dealed with a square cylinder but first think come into my mind: You should think about your Reynoldsnumber and the regime. May be you have to adjust the model constants of our oneEqnEddy SGS model. It could also be that you are using the wrong SGS model.
I also would refine the mesh in the wake of your square cylinder, not only at the surface. Have you considered to get rid of the wallfunction and resolving the viscous sublayer (y+ < 1)? Fabian |
|
August 29, 2012, 12:52 |
|
#9 |
Senior Member
Julien
Join Date: Jun 2012
Location: France
Posts: 152
Rep Power: 13 |
Fabian,
My Reynolds number is about 22 000. I didn't know the SGS model was dependent on this Re ? Do you have any suggestion instead of "regular" OneEqEddy (saying with default parameters) ? Concerning my mesh, my moderate computing power limits me: refining either the wake or the cylinder's surface would imply a 2D study. Do you think it should be better? |
|
August 29, 2012, 13:13 |
|
#10 |
Senior Member
Dr. Fabian Schlegel
Join Date: Apr 2009
Location: Dresden, Germany
Posts: 222
Rep Power: 17 |
Most of the SGS models are developed for a specific case. As far as I know there exists no SGS model which could be applied to universal cases. Often the coefficients vary if you deal with a low tubulent flow (cylinder at Re=3900) or a high turbulent flow (cylinder 180000) or so. The model is tricky I guess. I would check other papers to see which model the authors use. Probably a ordinary smagorinsky (SM) is not to bad. I would try more than one model and see how the Strouhal number behaves. But especially with a SM model be carfull with the coefficients.
You have a very high Re number. I assume the regimes are comparable to circular cylinder flows. For my opinion, without a refinement in the wake you will never get a right Strouhal number. What are your boundary conditions for the inflow and outflow? You should also ensure that your domain is large enough, to ensure you have no effects from the downstream BC influencing your solution. I have no clue right now how far, but you definitly should check the various papers on cylinder flows for the dimensions. Ok final and most important guess: Take a paper for a circular cylinder, e.g. Wissink, J. & Rodi, W. Numerical study of the near wake of a circular cylinder International Journal of Heat and Fluid Flow, 2008, 29, 1060 - 1070 and try to set up the same case with OF and if you get the right Strouhal number and velocity profiles in the wake. Tune OF as long as you do not fit their results. Than switch back with the same setup to your squared cylinder and see what happens. Fabian P.S: 2D makes no sense in my opinion. This Re number is to high for assuming a mainly 2d flow. If it would be laminar, you can do 2d but not in this highly turbulent case. Please correct me anybody of I am wrong. |
|
August 29, 2012, 13:47 |
|
#11 |
Senior Member
Julien
Join Date: Jun 2012
Location: France
Posts: 152
Rep Power: 13 |
OK, thanks Fabian. I agree with you for everything now (I just don't want to make the same case as others, because I am over confident... and I naively think in this way, I will save time!!!).
'see you in two years when I will have done all what you are saying Sincerely, |
|
August 29, 2012, 14:55 |
|
#12 |
Senior Member
Dr. Fabian Schlegel
Join Date: Apr 2009
Location: Dresden, Germany
Posts: 222
Rep Power: 17 |
I do not think it will take you 2 years. But it would be the right way. The basic setup you have to figure out anyhow. If you do it with your squared cylinder there is much guessing and less knowing. It will coast you very much time to find the right setup. If you calculate a circular cylinder frist and compare it to known results, is less guessing and you get a good setup very quick. Than the only thing you have to change is the geometry (basic mesh should remain the same) and change the inlet velocity. The postprocessing is also established already, because you did need it for the validation and your done. Think about it again :-D
Fabian |
|
September 3, 2012, 05:14 |
|
#13 |
Member
|
First of all: how did you calculate the Strouhal number?
i.e. Which parameter did you base it on? How did you get the principal frequency? Regarding your solution it seems to coherent for being a properly resolved LES. I guess this comes from a too coarse mesh and also setting nuSgs at the inlet to 1e-1 makes no sense (if evolving fluid is air or something with equivalent laminar viscosity). Regards, Cosimo
__________________
Cosimo Bianchini Ergon Research s.r.l. Via Panciatichi, 92 50127 Florence - ITALY Tel: +39 055 0763716 Mob: +39 320 9460153 e-mail: cosimo.bianchini@ergonresearch.it URL: www.ergonresearch.it |
|
September 3, 2012, 06:08 |
|
#14 |
Senior Member
Julien
Join Date: Jun 2012
Location: France
Posts: 152
Rep Power: 13 |
Good Morning Cosimo ! And thanks for your attention.
My Strouhal number is evaluated with the width of the cylinder (dimension perpendicular to the flow, here the shortest side, 0.1 m), the upsteram homogeneous velocity (here 3.3m/s). I put a probe in the wake, and my frequency is evaluated thanks to the Ux velocity time series (see post 378281). 'sure my small eddies are not resolved here. But I don't mind! I want the large ones, specially the vortex shedding. I thought the principle of LES was exactly this: resolve large eddies and model small ones. Am I wrong? In fact, the effective BC concerns k, not nuSGS (calculated from k). And my inlet-k is 0.001. Calculated nuSGS is about 5e-5. Do you think I have to refine my mesh near the cylinder or in the wake? Sincerely |
|
September 3, 2012, 06:46 |
|
#15 | |||
Member
|
Quote:
Quote:
Quote:
You can read this paper by Prof.Davidson to have an overview on methods employed to verify the reliability of your LES. http://www.tfd.chalmers.se/~lada/pos...r_qles2009.pdf Probably both. Indeed using an almost uniform grid in the entire domain is not a very efficient strategy. You should refine more where you aspect to have most of the smallest turbulent structures.
__________________
Cosimo Bianchini Ergon Research s.r.l. Via Panciatichi, 92 50127 Florence - ITALY Tel: +39 055 0763716 Mob: +39 320 9460153 e-mail: cosimo.bianchini@ergonresearch.it URL: www.ergonresearch.it |
||||
September 3, 2012, 07:28 |
|
#16 | |
Senior Member
Julien
Join Date: Jun 2012
Location: France
Posts: 152
Rep Power: 13 |
Thanks a lot. I will read this paper. I am still young in LES problems and I don't "feel" if a LES simulation is OK or not.
Nevertheless, I disagree with you about one point: Quote:
So my huge error in frequency (x3) does not come from this estimation error (x2). Do you think my bad mesh is responsible of multiplying by 3 the vortex shedding frequency ? |
||
February 6, 2013, 12:51 |
|
#17 |
Member
ak
Join Date: May 2011
Posts: 64
Rep Power: 14 |
Hi Julien
Any updates on this? Were you able to get the correct St #? Many thanks, Amit |
|
February 11, 2013, 07:54 |
|
#18 |
Senior Member
Julien
Join Date: Jun 2012
Location: France
Posts: 152
Rep Power: 13 |
Hi everybody!
Excuse me for my long silence: in my language, we say "No news? Good news!". Generally speaking, I write my posts when I have problems. So no post = no problem! Nevertheless, I feel bad thinking of all this kind community dealing only with problems... Thus, I would like to share those good results finally I had with FOAM, concerning the vortex shedding. Basically, I wanted to simulate the vortex-shedding frequency of a rectangular beam, depending on its aspect ratio. This was also to train myself with Foam, LES, meshing,... I worked a lot, with a lot of different parameters: I would not be able to detailed all the steps. But changing things and things, finally I had the correct frequency on the 1:1 case, and then changed the aspect ratio, or certain parameters. Thus, I can say what changes do not work. I cannot explain what to do to have the correct situation, but I know from the correct situation what to do to have a bad case! Basically, I wanted to reproduce this picture from the Eurocode I (NF EN 1991-1-4 Annex E): My case is tuned like this:
Not bad, is it? In fact, the situation is not so clear, because is several cases, there are several "modes" in the wake. Is those cases, there are several peaks in the Fourrier spectrum, with no relations (such as harmonics: octave, fifth,...). Here you can see those additional modes: One of my mystakes has been not to wait enough time: during the transient part, other modes can be wider than the stable mode! I cannot explain these several modes. Nevertheless, I cannot explain neither the high complexity in the Eurocode curve... My guess is that there is actually competition between several modes, at each aspect ratio, but one mode is accentuated: the Eurocode curve should be parts of different curves. But this is a long other story... Now, more basic considerations: by changing some parameters, the calculation fails. Here are those results: Some comments about them:
|
|
February 11, 2013, 08:38 |
|
#19 |
Senior Member
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 21 |
Thanks for sharing!
You are using a uniform inlet velocity (as far as I read here). Did you study the effect of this? How long is your inlet section? |
|
February 11, 2013, 09:09 |
|
#20 |
Senior Member
Julien
Join Date: Jun 2012
Location: France
Posts: 152
Rep Power: 13 |
Yes, I use a uniform inlet. I think it is maybe a bit short, but it is the same as a my reference data (AGARD AR 345 , CMP20 by Lyn, Einav, Rody and Park). Nevertheless, I don't think it affects the Vortex Shedding frequency.
|
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
K-Epsilon for Vortex Shedding | Sham | FLUENT | 33 | March 29, 2017 06:48 |
Cylinder vortex shedding (3D) | Apocolapse | STAR-CCM+ | 2 | April 3, 2014 21:33 |
Vortex shedding behind cylider in cross flow | Muthu | FLUENT | 0 | March 6, 2006 11:29 |
LES simulation of vortex shedding | Hans | Main CFD Forum | 8 | September 25, 2003 07:28 |
LES: vortex shedding past a 2D and 3D circular cylinder | Rocchi Daniele | Main CFD Forum | 7 | April 9, 1999 13:05 |