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

viscosity ratio in wake of ellipse - physical meaning?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 16, 2013, 08:59
Question viscosity ratio in wake of ellipse - physical meaning?
  #1
Member
 
Join Date: Aug 2011
Posts: 54
Rep Power: 5
MachZero is on a distinguished road
Hello all,

I am trying to model and capture the separation length behind an elliptical cylinder. I have done so reasonably well with k-omega SST but would like to move to an RSM model to get a more certain answer. My grid is high quality in the wake and I am using double precision.

As I activate 1st order RSM, I get a warning that viscosity ratio is greater than 1e5 and then it diverges. I checked my solution and this occurs in the wake. Reading one of ANSYS's response, they state:

Also note that for certain types of flows, for instance atmospheric boundary layer flows, the turbulent viscosity ratio can be as high as 1e+08 or 1e+09. In such cases, the limit can be increased.

Does anyone have any physical idea what range of viscosity ratios are possible in a wake? Should I actually increase this value? I dont have much physical understanding of turbulent viscosity ratio.
MachZero is offline   Reply With Quote

Old   July 16, 2013, 09:30
Default
  #2
Member
 
Reza
Join Date: Mar 2009
Location: Appleton, WI
Posts: 99
Rep Power: 8
triple_r is on a distinguished road
No idea why you are getting turbulence viscosity warnings with RSM, as the whole point of RSM is not having to use the Boussinesq approximation, it might be a post-processing thing. Anyhow, turbulence viscosity value depends on the grid resolution as well as flow conditions, so if the grid is coarse you might get very high values. Have you tried using the SST solution as the initial condition for the RSM run?
triple_r is offline   Reply With Quote

Old   July 16, 2013, 09:53
Default
  #3
Senior Member
 
flotus1's Avatar
 
Alex
Join Date: Jun 2012
Location: Germany
Posts: 1,098
Rep Power: 19
flotus1 will become famous soon enoughflotus1 will become famous soon enough
AFAIK, the turbulent viscosity also depends on the Reynolds number. That is why it can be higher in atmosperic boundary layers.
If the Reynolds number of your flow setup is not exeedingly high, the warning about limited TVR should vanish after some iterations.
If it doesnt, review the boundary conditions and, as triple_r already mentioned, the mesh size and quality in the boundary layer and in the wake region.
flotus1 is offline   Reply With Quote

Old   July 16, 2013, 10:04
Default
  #4
Member
 
Join Date: Aug 2011
Posts: 54
Rep Power: 5
MachZero is on a distinguished road
Thank you for your quick responses!

I tried initially using SST values for initialization, but found that when I switched to RSM, it diverged. Looking over ANSYS notes, they recommend initializing using k-e, then switching to RSM. So that is what I have been doing since. I have been using low UR factors in momentum and rsm, and starting with 1st order before switching to 2nd.

I just checked the solution history. The warning is:

turbulent viscosity limited to viscosity ratio of 1.000000e+05 in 19 cells. I will attach an image of my mesh quality. Maybe it really isnt high enough.

My Re is about 1,000,000.

Thank you both for responses.
MachZero is offline   Reply With Quote

Old   July 16, 2013, 10:14
Default
  #5
Member
 
Join Date: Aug 2011
Posts: 54
Rep Power: 5
MachZero is on a distinguished road
Under closer inspection, I see on the midplane of my "wind tunnel", the TVR is maxing out around 6000, BUT near the symmetry plane (3d simulation of elliptical cylinder using symmetry plane) is where I am getting very high values. I have attached two images that show this (range for picture with mesh is 0-6000, range for image with red blip is 0-85000).

High TVR values are not reducing with iterations, but rather generating in more cells with iterations. I have seen several posts on this forum seeing similar issues , so thank you for your responses! I would love to get this figured out.

midplaneCFD.png

extremeValueCFD.png
MachZero is offline   Reply With Quote

Old   July 16, 2013, 11:33
Default
  #6
Senior Member
 
flotus1's Avatar
 
Alex
Join Date: Jun 2012
Location: Germany
Posts: 1,098
Rep Power: 19
flotus1 will become famous soon enoughflotus1 will become famous soon enough
How about the Y+ values? And the aspect ratio of the cells where the maximum value of TVR occurs? Does this point coincide with the separation point?
flotus1 is offline   Reply With Quote

Old   July 16, 2013, 11:59
Default
  #7
Member
 
Reza
Join Date: Mar 2009
Location: Appleton, WI
Posts: 99
Rep Power: 8
triple_r is on a distinguished road
With regard to initialization, you are absolutely right, depending on the seventh equation for RSM you might be better off with k-e instead of k-o.

As flotus1 mentioned, you might want to check your grid around that blip for aspect ratio. If y+ is low, you can get rectangular elements stretched in the direction parallel to the wall, and usually that is a good way to save on elements as gradients are much higher normal to the wall than parallel to the wall. But, if that point is near the separation point, then you probably need more square-like elements.

From your pictures it looks like you are using a "mapped" mesh, if that is so then try slicing the domain radially right where the blip is, and add gradients to the line meshes on the two sides of the slice so that they become finer near the blip. Let me see if I can show this in a schematic picture :-) I made it in paint, so sorry for the bad quality :-p

triple_r is offline   Reply With Quote

Old   July 16, 2013, 12:02
Default
  #8
Member
 
Reza
Join Date: Mar 2009
Location: Appleton, WI
Posts: 99
Rep Power: 8
triple_r is on a distinguished road
Sorry, the picture didn't get attached, so hopefully here it is.
Attached Images
File Type: jpg grid.jpg (33.2 KB, 11 views)
triple_r is offline   Reply With Quote

Old   July 22, 2013, 16:32
Default
  #9
Member
 
Join Date: Aug 2011
Posts: 54
Rep Power: 5
MachZero is on a distinguished road
First off thank you all for your suggestions. I have taken a step back and taken some time to look at this file with reapect to all your comments.

First thing to note is that yes I have a y plus of 1. Next is that I have found that if I use an fmg initialization and run 1st order rsm immediately it does not diverge. I ran it for 10k iterations so I am pretty sure of that. I switched to second order and it diverged eventually. Note that the regions with ultra high tvr were at the symmetry planes.

An ansys engineer told me that I should use kw sst since it will give a better solution then switch to rsm. He said that when I do that I should freeze the flow and solve turbulence only for a while. I tried this approach with several variations and it still diverged.

Last thought is to build a coarser mesh and use non equilibrium wall functions. Should I be worried about this killing bl calculations and separatipn location calculation since it is not being resolved?
MachZero is offline   Reply With Quote

Reply

Tags
reynolds stress model, separated flow, turbulent viscosity ratio

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
Ratio of eddy viscosity to molecular viscosity : Laminar or turbulent flow? RicochetJ CFX 7 September 9, 2013 07:45
Problem with divergence TDK FLUENT 10 September 8, 2012 01:11
pressure eq. "converges" after few time steps maddalena OpenFOAM Running, Solving & CFD 69 July 21, 2011 07:42
Meaning of Backflow Turbulent Viscosity Ratio? Remy FLUENT 1 February 4, 2008 10:19
Problem of Turbulent Viscosity Ratio Limited David Yang FLUENT 3 June 3, 2002 06:13


All times are GMT -4. The time now is 07:13.