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

Problem predicting transition point.

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 23, 2013, 11:02
Default Problem predicting transition point.
  #1
New Member
 
CFD User
Join Date: Feb 2012
Posts: 18
Rep Power: 6
aditipatel7 is on a distinguished road
Hi there,

Software Used: Fluent, Mesh Generator: Gambit

I am simulating a flow in a turbine cascade. Please find my geometry in the attached picture. My inlet velocity is 1.05m/s and Re = 72000. I am using the SST k-w model. My chord length is 1m. As per the model and other settings, transition happens at 0.2 m into the blade on the suction side i.e. the lower blade in the geometry.

However, transition is expected to be between 0.4 and 0.6. Please suggest some ideas for the above problem.
Attached Images
File Type: png 123.PNG (21.9 KB, 41 views)
aditipatel7 is offline   Reply With Quote

Old   April 1, 2013, 05:25
Default
  #2
Far
Super Moderator
 
Sijal
Join Date: Mar 2009
Location: Islamabad
Posts: 4,258
Blog Entries: 6
Rep Power: 43
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
Hey

I am solving similar problem. Do you need some help?
Far is offline   Reply With Quote

Old   April 2, 2013, 13:36
Default
  #3
New Member
 
CFD User
Join Date: Feb 2012
Posts: 18
Rep Power: 6
aditipatel7 is on a distinguished road
It feels this forum is just dead! No point putting questions up here nymore!
aditipatel7 is offline   Reply With Quote

Old   April 2, 2013, 23:11
Default
  #4
Far
Super Moderator
 
Sijal
Join Date: Mar 2009
Location: Islamabad
Posts: 4,258
Blog Entries: 6
Rep Power: 43
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
What is the problem exactly? You just said you are solving problem involving transition and you are expecting is at 0.2 m location ! 0.2 m from leading edge along the axial chord or on the suction side?

Some more details are required like :

1. How many no of nodes you are using. What is overall mesh? How many nodes you are using in the transition region?

2. What is your Yplus?

3. What are the turbulence parameters at inlet and are the same as per experimental set-up?
Far is offline   Reply With Quote

Old   April 5, 2013, 18:09
Default
  #5
New Member
 
CFD User
Join Date: Feb 2012
Posts: 18
Rep Power: 6
aditipatel7 is on a distinguished road
Thanks for coming forward to help me out. My problem is that I am not able to get the correct transition point as per the journal I am referring to.

For an intensity of 2.8%, my transition point is at 0.6m into the blade on the suction side (lower edge in my geometry). So for a higher intensity i.e. 8.4%, I am expecting the transition point to be between 0.4 and 0.5. I am basing this understanding on Wissink's "Direct Numerical Simulation of flow and heat transfer in a turbine cascade with incoming wakes"-page 231 (Nu Graph). Can you suggest some changes in my setting or give me other ideas?

Please find my geometry and Nu plot in the attached pictures. My inlet velocity is 1.05m/s and Re = 72000. My chord length is 1m.

Software Used: Fluent, Mesh Generator: Gambit

Following are my settings:
1. Refined mesh at the edges to predict accurate boundary layer. y+ is <1
2. I am using the Transition SST-4 equation model with Viscous heating and curvature correction on. Apart from that Model constants are default values
3. I am aiming for an intensity of 8.4% so inlet conditions are set to Intermittancy as 0.64, Intensity 8.4 and length scale as 1 (based on my chord length).

Other settings are "Third-Order Muscle" and "Power law" for Spatial Discretization section in Fluent. My monitors except for continuity are all set to 1e-06.
Attached Images
File Type: jpg New Mesh.JPG (71.4 KB, 10 views)
File Type: jpg Nu graphs for New Mesh.JPG (41.0 KB, 26 views)
aditipatel7 is offline   Reply With Quote

Old   April 5, 2013, 18:16
Default
  #6
New Member
 
CFD User
Join Date: Feb 2012
Posts: 18
Rep Power: 6
aditipatel7 is on a distinguished road
I can send you my case and data files if that works and can help me out quicker.
aditipatel7 is offline   Reply With Quote

Old   April 6, 2013, 00:57
Default
  #7
Far
Super Moderator
 
Sijal
Join Date: Mar 2009
Location: Islamabad
Posts: 4,258
Blog Entries: 6
Rep Power: 43
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
There are many suggestion I can give you to refine your case:

1. It is possible that at higher turbulence intensity there is fully turbulent flow.

2. Don't confuse wakes with turbulence intesnity. Wakes are something trasient flow phenomena and occurs periodically when downstream blade is moving at certain rpm.

3. To have wake effect, you can implement trasient gaussian profile or do the precursor simulation of upstream wakes by just modelling a suitable diameter cylinder.

4. The point where you are expecting transition, mesh should be refined with 20-30 nodes in the vicinity of that location.

5. You should plot mean CP as CP is changing from time step to time step. For this you have turn-on time statistics data sampling after the simulation has achieved periodic convergence.

6. Divide blade into to parts. Which you already have pressure and suction sides. Specify them separately in boundary conditions in gambit. When you plot Cp you will observe there is option of data write. Write data sepratly using the option curve length.

7. Open these files in excel and divide each dimension (curve length) by maximum value. In this way you will normalize curve length between 0 and 1 and now you have suction and pressure side on the same scale otherwise they had different lengths.

7. Did you calculate boundary layer thickness through laminar or turbulent flat plate formula to estimate the boundary layer thickness in your case.

8. Are you using second order schemes for all variables including turbulence.

9. Are you using coupled or piso scheme along with second order implicit time scheme?

PS: Send your reference paper, case and data files.
Far is offline   Reply With Quote

Old   April 7, 2013, 04:09
Default looks like u are solving a similar problem in LES or DNS using other flow analysis S/
  #8
New Member
 
CFD User
Join Date: Feb 2012
Posts: 18
Rep Power: 6
aditipatel7 is on a distinguished road
Thanks for all the suggestions. I have sent the reference paper, case and data files on turboenginner@gmail.com

1.It is possible that at higher turbulence intensity there is fully turbulent flow.
A1. On pressure side I believe it must be fully laminar looking at the velocity contour plot. IS there a method to tell at which point on the blade the flow becomes turbulent??

2. Don't confuse wakes with turbulence intesnity. Wakes are something trasient flow phenomena and occurs periodically when downstream blade is moving at certain rpm.
A2. I understand the difference between wakes and intensity. Its just that I was asked to simulate the cascade model without wakes in 2D RANS and was just asked to match the intensity at the measurement point. i.e 0.2m in front of the blade. (page 215-Wissink's Journal)

3. To have wake effect, you can implement trasient gaussian profile or do the precursor simulation of upstream wakes by just modelling a suitable diameter cylinder.
A3. Thanks for this suggestion, I can mention it in my further work.

4. The point where you are expecting transition, mesh should be refined with 20-30 nodes in the vicinity of that location.
A4. I believe the mesh is well refined at the edges but do I have to refine it in the section between the blades as well? Since it is left coarse in the journal as well.

5. You should plot mean CP as CP is changing from time step to time step. For this you have turn-on time statistics data sampling after the simulation has achieved periodic convergence.
A5. Is periodic convergence applicable for RANS?
If in case of RANS, do u mean plot Cp on Suction and Pressure Side using the dropdown options in Fluent. I get a fish shaped graph but not as exact as Wissink's Journal.

6. Divide blade into to parts. Which you already have pressure and suction sides. Specify them separately in boundary conditions in gambit. When you plot Cp you will observe there is option of data write. Write data sepratly using the option curve length.
A6. Already done that as you can see S/So on my Nu Graph added above in this discussion

7. Open these files in excel and divide each dimension (curve length) by maximum value. In this way you will normalize curve length between 0 and 1 and now you have suction and pressure side on the same scale otherwise they had different lengths.
A7. Already done that.

7. Did you calculate boundary layer thickness through laminar or turbulent flat plate formula to estimate the boundary layer thickness in your case.
A8. Is it using this formula: Delta = 0.382*x/Re_x ^ (1/5)?

8. Are you using second order schemes for all variables including turbulence.
A9. I have tried running basic simulation (Default parameters) using SIMPLEC, PISO and COUPLED while using Second Order or Third Order for Spatial Discretization. However for SST-Transition model.. only SIMPLEC with Third Order Muscle gives reasonable results.

9. Are you using coupled or piso scheme along with second order implicit time scheme?
A10. I haven't seen second order implicit scheme for RANS.
aditipatel7 is offline   Reply With Quote

Old   April 9, 2013, 07:49
Default
  #9
New Member
 
CFD User
Join Date: Feb 2012
Posts: 18
Rep Power: 6
aditipatel7 is on a distinguished road
Nd its dead again! Anyone there with useful suggestions?
aditipatel7 is offline   Reply With Quote

Old   April 9, 2013, 08:11
Default
  #10
Far
Super Moderator
 
Sijal
Join Date: Mar 2009
Location: Islamabad
Posts: 4,258
Blog Entries: 6
Rep Power: 43
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
Improve your mesh.
Far is offline   Reply With Quote

Old   April 9, 2013, 18:15
Default
  #11
New Member
 
CFD User
Join Date: Feb 2012
Posts: 18
Rep Power: 6
aditipatel7 is on a distinguished road
Hi Far,
Have u looked at my email, case and data files yet? I cannot improve my mesh in Gambit anymore coz I am running out of time.

If possible can you answer the questions I have above (That long blue and black reply)

Thanks
aditipatel7 is offline   Reply With Quote

Old   April 13, 2013, 17:39
Default
  #12
New Member
 
CFD User
Join Date: Feb 2012
Posts: 18
Rep Power: 6
aditipatel7 is on a distinguished road
Improved the mesh and the boundary layer, used the same settings as mentioned above or as seen in my case and data files. Still No change in the results.

If possible can I request you to answer some of my doubts raised above.

Any help with be appreciated.

Thanks
aditipatel7 is offline   Reply With Quote

Old   April 14, 2013, 02:05
Default
  #13
Far
Super Moderator
 
Sijal
Join Date: Mar 2009
Location: Islamabad
Posts: 4,258
Blog Entries: 6
Rep Power: 43
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
pls see this thread: Time step for Low pressure turbine simulation
Far is offline   Reply With Quote

Old   April 17, 2013, 07:47
Default
  #14
New Member
 
CFD User
Join Date: Feb 2012
Posts: 18
Rep Power: 6
aditipatel7 is on a distinguished road
Hi Far, if possible can you please email me your case and data files for your study in Transition SST model? I would like to have a look at your settings in detail to understand if I can incorporate some of them into my model.

As I am studying my model in steady state solver, it appears to be quiet hard to get a stable solution.
aditipatel7 is offline   Reply With Quote

Old   April 17, 2013, 09:18
Default
  #15
Far
Super Moderator
 
Sijal
Join Date: Mar 2009
Location: Islamabad
Posts: 4,258
Blog Entries: 6
Rep Power: 43
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
Where to send? Your email address please...
Far is offline   Reply With Quote

Old   April 17, 2013, 09:56
Default
  #16
New Member
 
CFD User
Join Date: Feb 2012
Posts: 18
Rep Power: 6
aditipatel7 is on a distinguished road
its ......

Thanks in advance.

Last edited by aditipatel7; April 17, 2013 at 22:06.
aditipatel7 is offline   Reply With Quote

Old   April 17, 2013, 10:11
Default
  #17
New Member
 
CFD User
Join Date: Feb 2012
Posts: 18
Rep Power: 6
aditipatel7 is on a distinguished road
ok I just have one basic question here...
How can you tell at which point the flow transitioned from laminar to turbulent? I know we can do u+ vs y+ for a flat plat boundary layer profile.
But for an airfoil that is curved...what is the method?

Last edited by aditipatel7; April 17, 2013 at 22:39.
aditipatel7 is offline   Reply With Quote

Old   April 17, 2013, 22:39
Default
  #18
New Member
 
CFD User
Join Date: Feb 2012
Posts: 18
Rep Power: 6
aditipatel7 is on a distinguished road
ok I just have one basic question here...
How can you tell at which point the flow transitioned from laminar to turbulent? I know we can do u+ vs y+ for a flat plat boundary layer profile.
But for an airfoil that is curved...what is the method?
aditipatel7 is offline   Reply With Quote

Old   April 18, 2013, 02:48
Default
  #19
Far
Super Moderator
 
Sijal
Join Date: Mar 2009
Location: Islamabad
Posts: 4,258
Blog Entries: 6
Rep Power: 43
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
Hmmm. there are many methods to do this...

1. Plot intermittency

2. Plot integral parameters (e.g shape factor)

There may be other methods to do this. You need to explorer relevant journal research papers to see how people have decided about the two regions and how they infer transition from other flow variables.
Far is offline   Reply With Quote

Old   August 20, 2014, 04:41
Default Shape factor
  #20
Member
 
venkatesh
Join Date: May 2012
Posts: 93
Rep Power: 6
venkat_aero2007 is on a distinguished road
I am analyzing transition point over NACA 4412 airfoil at low Reynolds number using Transition Kw-SST model. I am using FLUENT for simulation. I dont know how to plot Shape factor from FLUENT result. I know that Shape factor can be calculated using the relation
where H is the shape factor, is the displacement thickness and θ is the momentum thickness.The definition of displacement thickness, δ for incompressible flow is
Where and are the density and velocity in the 'free stream' outside the boundary layer, and is the coordinate normal to the wall.


The definition of Momentum thickness for incompressible flow is

I read in a post that velocity profile can be obtained by creating a line perpendicular to the surface at the point of interest. and with the help of that line velocity variation can be plotted.



I don't know how to plot shape factor in fluent. Can anybody please help me in this regard.
venkat_aero2007 is offline   Reply With Quote

Reply

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
Problem in applying Gamma-Retheta transition model mb.pejvak Main CFD Forum 2 November 13, 2012 11:27
problem in turbulent flow over flat plate with stagnation point mb.pejvak Main CFD Forum 0 September 12, 2012 21:25
Trans. SST Intermittency Factor and Viscosity Ratio eishinsnsayshin FLUENT 3 May 23, 2012 03:02
IFStream read float point data problem liu OpenFOAM Running, Solving & CFD 0 October 24, 2008 12:14
Transition Point shape factors!! please help Simon c FLUENT 0 March 12, 2006 12:00


All times are GMT -4. The time now is 00:47.