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

Need URGENT HELP for convergence problem of a turbine

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 30, 2011, 18:31
Unhappy Need URGENT HELP for convergence problem of a turbine
  #1
New Member
 
Kun Lu
Join Date: Jun 2010
Posts: 4
Rep Power: 16
Pursuor is on a distinguished road
Hello everyone,

I encountered convergence problem when I tried to do simulation of a turbine. I have spent more than two weeks but still can not fix the thing. So hopefully you would help me figure it out or give me some suggestions.

The model is a three-stage turbine. We may call it 'baseline' (Fig. 1 top-left). Then I added endwall contouring on the first rotor (R1), named as 'contoured' (Fig. 1 top-right). For these two case, everything is fine. Usually it takes 300 time steps to get converged. Then I went to 'baseline with purge' (Fig. 1 bottm-left) which has purge flow at upstream of R1. This case also got converged in 500 time steps. After this, I applied purge flow to contoured case 'contoured with purge' (Fig. 1 bottom-right). Then the trouble keeps bothering me. By the way, the turbulence model is SST model.

At the beginning, I used High Resolution for both Advection Scheme and Turbulence Numerics. I got smooth but pretty high residual for Mass (red curve in Fig. 2) and Heat Transfer (Fig. 3), which is much higher (1-2 orders) than what I got from previous cases. This was due to the imbalance within R1. At the same time, the unreasonable result which was not physics was obtained computational domain R1, see temperature distribution in Fig. 4.

Then I tried to enhance the quality of the mesh several times. It never works. So I realized it should not be the problem of mesh, as I used the similar mesh for previous case successfully. I also tried to change the boundary conditions. It did not work either.

Finally, I turned to Timescale. This time I set local timescale factor as 5, using First Order and blender factor 0.75. Although the residual looks better than before (Fig. 5), the residual significantly increased after 500 time steps and became to oscillate. It is quite strange that the high RMS-Mass was found in R3 but not in R1 and thus the temperature distribution in R1 seems reasonable. But I did not see any trend that the residual may drop to less than 1e-6 as previous cases which were all converged. Moreover, after I changed back to physical timescale and high resolution scheme, R1 became my headache again. So my questions are:

1. Does the strong unsteady property of the flow cause this problem?

2. Since the strong unsteady property exists in the 'baseline with purge' too, does the complicated geometry of the endwall contouring lead to instabilities?

3. Could I use the results calculated with local timescale factor as final results?

4. What is the possible solution for this problem if I can not afford to perform transient calculation?

This bothered me a lot, so I eagerly look forward to your solutions or suggestions. Thank you very much.
Attached Images
File Type: jpg cases.jpg (42.4 KB, 50 views)
File Type: jpg RMS-momentum and mass.jpg (91.4 KB, 46 views)
File Type: jpg RMS-Heat Transfer.jpg (94.9 KB, 39 views)
File Type: jpg Temperature.jpg (25.3 KB, 37 views)
File Type: jpg new RMS-momentum and mass.jpg (64.3 KB, 43 views)

Last edited by Pursuor; October 1, 2011 at 05:29.
Pursuor is offline   Reply With Quote

Old   October 1, 2011, 03:28
Default
  #2
Far
Senior Member
 
Sijal
Join Date: Mar 2009
Location: Islamabad
Posts: 4,558
Blog Entries: 6
Rep Power: 54
Far has a spectacular aura aboutFar has a spectacular aura about
Send a message via Skype™ to Far
http://www.cfd-online.com/Forums/cfx...nvergence.html

this might be useful for you
Far is offline   Reply With Quote

Old   October 2, 2011, 02:48
Default
  #3
New Member
 
Kun Lu
Join Date: Jun 2010
Posts: 4
Rep Power: 16
Pursuor is on a distinguished road
Quote:
Originally Posted by Far View Post

Thank you for your reply. I went to the link and reviewed your discussion. I tried to use double precision, but it runs too slowly that I can not afford it. By the way, I am using 10 cpus with total 66GB memory for the running. Is it because the number of the cpus is not enough or the memory is small or I did not set that correctly?

For the mesh, I am afraid that I can not increase the height of the first grid as I use SST turbulence model.
Pursuor is offline   Reply With Quote

Old   October 2, 2011, 04:42
Default
  #4
Far
Senior Member
 
Sijal
Join Date: Mar 2009
Location: Islamabad
Posts: 4,558
Blog Entries: 6
Rep Power: 54
Far has a spectacular aura aboutFar has a spectacular aura about
Send a message via Skype™ to Far
well, I have no clear solution for this problem, this might be due to mesh, flow might be transient or any thing else you name it.

But what I wanted to say is that, use the default time scheme that is auto time scale and set the time scale factor to 0.1 and then run the solution for the 1000-2000 time steps. During this phase you may occasionally have the some divergence, don't worry, ultimately the solution converge to the required tolerance provided that any steady solution exists. (my last reply in the thread SMC models convergence)

PS.

1. With automatic wall treatment you don't need to have the Y+ to be ~1. Y+~ 10 should do the same magic.
2. Can you tell me (or attach the line plot of Y+) the y+ values at the leading edge of rotor 1 (in other words Y+ in the stagnation region)
Far is offline   Reply With Quote

Old   October 3, 2011, 06:52
Default
  #5
Member
 
Join Date: Dec 2009
Posts: 57
Rep Power: 17
Graham81 is on a distinguished road
Pursuor,

Yes the cause may well be unsteady flow patterns, assuming you are now running frozen rotor or mixing plane simulations. The obvious fix is to do a transient simulation instead.

If you dont need the increased accuracy of a transient calculation or just dont have the time, I would suggest setting a physical timescale to a typical residence time as the CFX manual suggests. For turbomachinery: use 1/(2*pi*rpm [min^-1]) where rpm is ofcourse the speed of your machine.

Hope this helps,

Graham
Graham81 is offline   Reply With Quote

Old   October 3, 2011, 18:44
Default
  #6
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,871
Rep Power: 144
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
This looks a FAQ question to me:

http://www.cfd-online.com/Wiki/Ansys...gence_criteria
ghorrocks is offline   Reply With Quote

Reply


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
urgent help needed (rhie-chow interpolation problem) Ardalan Main CFD Forum 2 March 18, 2011 16:22
Problem with interface in a wind turbine simulation be_inspired FLUENT 2 March 8, 2011 00:46
Gas turbine load problem hashimbukhari Main CFD Forum 2 July 24, 2010 04:48
Multiphase Problem (Urgent) Dadang CFX 4 June 21, 2004 08:46
Urgent help on Sand casting problem arunprasad FLUENT 0 April 18, 2004 06:04


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