CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > CFX

Low Reynolds Number SST Model

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

Like Tree4Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 29, 2010, 17:25
Default Low Reynolds Number SST Model
  #1
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 17
Josh is on a distinguished road
Hi all -

I am modeling the T106A low-pressure turbine, shown here:



I am using the SST turbulence model. Here is my pressure distribution graph:



Ignore the orange line (Tu = 1%). Both the red and green lines (my initial results and my results with a quadrupled grid density, respectively) correlate alright with the experimental, but don't capture the laminar separation bubble on the suction surface (the experimental results have a pressure plateau).

At a recent fluids conference, a PhD student told me he obtained very good results (in comparison with Stieger's, found here: http://www-g.eng.cam.ac.uk/whittle/T106/Start.html) when he used the SST model with the low Reynolds feature turned on. I know the low Reynolds number feature exists in Fluent, but in CFX-Pre, I am unsure of how to turn off the automatic wall function.

After digging through the Help file, I have found many topics about the low Reynolds model, such as the required node refinement (at least 15 within the boundary layer for the low Re model, or 10 for the automatic wall function model), but I can't find any info on actually enabling the low Re model.

Does anyone know how to enable the low Re model and/or turn off the automatic wall function?
Josh is offline   Reply With Quote

Old   April 30, 2010, 06:33
Default
  #2
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,696
Rep Power: 143
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
You results are very impressive, very good agreement with experiment, well done.

Rather than change the wall function approach why not turn on the turbulence transition model and try to directly model the laminar to turbulent transition? This will be your best bet to get even better accuracy in my opinion.

But if you REALLY want to change the wall modelling you will find the wall function options under the turbulence model in the domain tab. By default the SST mode has the automatic option enabled.
ghorrocks is online now   Reply With Quote

Old   April 30, 2010, 20:34
Default
  #3
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 17
Josh is on a distinguished road
Thanks for the kind remark, Glenn. Unfortunately, my surface velocity and skin friction results aren't nearly as good.

I'll try what you said and post my new results. Thanks!
Josh is offline   Reply With Quote

Old   May 4, 2010, 12:05
Default
  #4
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 17
Josh is on a distinguished road
When I used the SST model without transitional turbulence, my solution did not converge, but it did oscillate very close to the residual values (and produced the above pressure results, which are pretty close to accurate). I didn't expect it to converge fully, as it is a highly three-dimensional, unsteady problem that I am modeling as two-dimensional and steady. When I tried using the transitional Gamma Theta Model, my residuals oscillated once again, but nowhere near the desired residual value.

A plot of y+ reveals maximum y+ values of 2.5 at the trailing and leading edges, and average y+ values < 1. Although CFX recommends a max y+ value of 1, the documentation said that results should not change significantly for 0.001 < y+ < 5, so a max of 2.5 seems reasonable.

CFX also recommended an expansion ratio of 1.1 with a hexahedral mesh, which I used.

I used the high resolution advection scheme, as recommended.

75-100 nodes in the streamwise direction are recommended. I used more, but that shouldn't drastically affect the results. They found that 30,000 nodes in most turbomachinery cases produces grid-independent results. I used just over 30,000 nodes.

I have read and practiced "obtaining convergence" on CFD Online, in the CFX Help file, and from various other sources. I will try a lower order discretization scheme, but I don't feel that will significantly help my solution converge. All I did was switch from the default SST to the SST with the Gamma Theta Model.

Any recommendations?
Josh is offline   Reply With Quote

Old   May 4, 2010, 18:38
Default
  #5
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,696
Rep Power: 143
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
Quote:
the documentation said that results should not change significantly for 0.001 < y+ < 5, so a max of 2.5 seems reasonable.
This comment is applicable to the default SST model. It is not applicable to SST with turbulence transition. The turbulence transition model requires a finer mesh, but you should do a sensitivity study to find out exactly how fine. You will probably have to lower the max y+ to about 1.

The turbulence transition model is often hard to converge in a steady state simulation. This is because you now have a transition, and for many airfoils you will get a laminar separation bubble associated with it (the transition model is very good at picking up laminar separation bubbles). However these bubbles are often unstable and jiggle about, even though the remainder of the flow is stable. This means you are forced to consider transient simulations.
8cold8hot and AlexRonto like this.
ghorrocks is online now   Reply With Quote

Old   May 4, 2010, 20:49
Default
  #6
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 17
Josh is on a distinguished road
Thanks again, honey (I heard that's what you're getting called these days).

I am trying to capture the laminar separation bubble, and I realize that this is an unsteady, three-dimensional problem in which the LSB likely fluctuates in size and location. However, I read a paper that was able to capture the LSB with steady, 2D simulations, so I thought I'd give it a whirl.
Josh is offline   Reply With Quote

Old   May 5, 2010, 06:44
Default
  #7
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,696
Rep Power: 143
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
Quote:
Thanks again, honey (I heard that's what you're getting called these days).
Yup, that's me

Quote:
However, I read a paper that was able to capture the LSB with steady, 2D simulations, so I thought I'd give it a whirl.
OK, but sounds like a problem to me. You might be really lucky with large timesteps that the turbulence model can stabilise it but you would need lots of luck. Your best bet is just to run transient.
ghorrocks is online now   Reply With Quote

Old   May 5, 2010, 12:24
Default
  #8
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 17
Josh is on a distinguished road
Hey Glenn -

I tried the simulation with my original, unrefined (though higher quality) mesh, SST Gamma Theta Transition Model, and adaptive timesteps (with the initial, minimum, and maximum timesteps based on experimental data). I captured the LSB (teal line):



As for the location/size, we think my geometry is a little off (i.e. the incidence angle is incorrect) and that this is causing added acceleration over the suction surface. We came to this conclusion because:
1) The velocity profiles on the suction surface are too high and show an unexpected spike at about mid-chord. The acceleration should be somewhat constant.
2) My laminar separation bubble is small and delayed, which could be because of added momentum in the boundary layer due to the extra acceleration.

Thanks again, Glenn.
Josh is offline   Reply With Quote

Old   May 7, 2010, 15:31
Default
  #9
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 17
Josh is on a distinguished road
Just an update for your consideration. I changed the incidence angle from 37.7 degrees to 40.5 degrees and obtained much better results. I believe this is because the stagnation point on my geometry was improperly set.

http://img208.imageshack.us/i/cpsteadyunsteadytu1.jpg/

Last edited by wyldckat; September 3, 2015 at 17:32. Reason: disabled embedded images
Josh is offline   Reply With Quote

Old   May 10, 2012, 03:52
Default Boundary condition values
  #10
Member
 
Arun
Join Date: Mar 2012
Location: Vellore,TN,India
Posts: 43
Rep Power: 14
arunintn is on a distinguished road
Hi Josh,
I'm doing the same analysis for the project. I checked and different sources i didn't find a good solution for the boundary value.

I have my value like given below.

Inlet
V=2.334 (for 1.6*10^5 re)
angle of Attack(37.7 deg)
x=2.33 cos(37.7)
y=2.33 sin(37.7)

but the problem is i didn't get the out let pressure value and the angle(do i need to consider angle on outlet pressure?)

i'm doing k and omega sst model

And the most worst case is i didn't get any bubble separation still now. my friends suggested the form. I hope i'll find some answer from you
arunintn is offline   Reply With Quote

Old   May 10, 2012, 11:03
Default
  #11
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 17
Josh is on a distinguished road
Outlet reference pressure is zero if you're specifying an inlet velocity. I don't understand your "outlet angle" question.

You're not seeing bubble separation because you're not using the proper turbulence model. The k-w SST model failed to predict bubble separation (see the above graphs). The SST gamma-Re_theta transition model was what I used to predict the bubble, and an unsteady analysis was required.
mrenergy and AlexRonto like this.
Josh is offline   Reply With Quote

Old   May 10, 2012, 16:55
Default
  #12
Member
 
Arun
Join Date: Mar 2012
Location: Vellore,TN,India
Posts: 43
Rep Power: 14
arunintn is on a distinguished road
thank you for setting me to the right track josh.. I'm doing my analysis in fluent and i'm new for CFD.

i calculated x and y by velocity triangle is this correct?
x=2.33 cos(37.7)
y=2.33 sin(37.7)
as the above post i'll change to 40.5 for my angle of attack.

Last edited by arunintn; May 10, 2012 at 17:29.
arunintn is offline   Reply With Quote

Old   May 10, 2012, 17:39
Default
  #13
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 17
Josh is on a distinguished road
I no longer have access to my simulations, but that looks correct. Post your results once you've obtained them.
Josh is offline   Reply With Quote

Old   May 10, 2012, 17:42
Default
  #14
Member
 
Arun
Join Date: Mar 2012
Location: Vellore,TN,India
Posts: 43
Rep Power: 14
arunintn is on a distinguished road
really i was little relaxed

I'll update my result after i finish my computation.
arunintn is offline   Reply With Quote

Old   May 17, 2012, 09:01
Default
  #15
Member
 
Arun
Join Date: Mar 2012
Location: Vellore,TN,India
Posts: 43
Rep Power: 14
arunintn is on a distinguished road
Hi Josh,
i completed the solution and it is converged but i did't get bubble separation and i have reverse flow in face(some numbers) and on pressure out (some numbers)

i have attached you the jpeg files for reference. I can't figure why i have revers flow. I hope that you can help in this.
Attached Images
File Type: jpg final_part1 vectorfull.jpg (60.7 KB, 119 views)
File Type: jpg final_part2 vector.jpg (59.6 KB, 86 views)
File Type: jpg final_part2 pressure plot.jpg (64.0 KB, 115 views)
File Type: jpg final_part2 y+.jpg (57.5 KB, 101 views)
arunintn is offline   Reply With Quote

Old   May 17, 2012, 18:38
Default
  #16
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,696
Rep Power: 143
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
The reverse flow must be small. I cannot see it on your plots. Are you using Fluent? Any difference when you use CFX?
ghorrocks is online now   Reply With Quote

Old   May 18, 2012, 07:39
Default
  #17
Member
 
Arun
Join Date: Mar 2012
Location: Vellore,TN,India
Posts: 43
Rep Power: 14
arunintn is on a distinguished road
Hi,
I have not used CFX before. I think it should be same.

I showed to my professor and he told it fine but he asked me to make the bubble separation. i tried a little but i can't find it. I need some idea here i need to know where and what change. I'm new to CFD and fluent
arunintn is offline   Reply With Quote

Old   May 18, 2012, 07:54
Default
  #18
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,696
Rep Power: 143
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
This is the CFX forum so we cannot help you specifically with Fluent. I would try this model in CFX as I feel much more comfortable with that solver, but it is your choice. Also note Fluent has lots of different sovlers to try, if you use the coupled solver you will be closer to the CFX solver.

To get laminar separation bubbles you will need the turbulence transition model and a fine, good quality airfoil mesh.
ghorrocks is online now   Reply With Quote

Old   April 8, 2013, 14:17
Default
  #19
New Member
 
Join Date: Jan 2012
Posts: 5
Rep Power: 14
littlek is on a distinguished road
Hi all,

I'm trying to do the same T106 analysis, but am not getting a correct velocity contour plot. I'm using inlet velocity of 8.45 m/s at an angle of 37.7. Just wondering if you guys have any tips for what I'm missing. Thanks very much.
Attached Images
File Type: jpg cas5.jpg (58.9 KB, 59 views)
littlek is offline   Reply With Quote

Old   April 8, 2013, 18:34
Default
  #20
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,696
Rep Power: 143
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
FAQ: http://www.cfd-online.com/Wiki/Ansys..._inaccurate.3F
ghorrocks is online now   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
DecomposePar unequal number of shared faces maka OpenFOAM Pre-Processing 6 August 12, 2010 09:01
Airfoil, LES and Low Reynolds number impecca OpenFOAM 1 July 23, 2010 11:59
Turbulence Model for low Reynolds Number Muhammad Shakaib Main CFD Forum 2 July 3, 2006 15:42
Turbulent Schmidt Number in SST model ? David CFX 0 December 5, 2005 04:43
About low Reynolds number airfoil experiment data. zqnwpu Main CFD Forum 5 December 25, 2004 03:52


All times are GMT -4. The time now is 21:42.