CFX transition model issue with VATT
Hi everyone,
I have a problem with 2D calculations of a straight bladed vertical axis tidal turbine using the Gamma  theta transition model of CFX.The mesh is divided in 3 parts : a rotating ring containing the 3 blades and 2 stators inside and outside the ring. I have run calculations at several tip speed ratios (TSR) with the fully turbulent SST turbulence model (called SST henceforth) and the SST turbulence model coupled to the gamma  theta transition model (called SSTTM henceforth). I expected differences between the results of the two models but I get big unrealistic differences at high TSR as you can see on the attached picture (averaged Power Coefficients (CP) have been divided by the SST value at TSR = 3 since they have not been published yet). SST calculations seem fairly good since they agree with data obtained by Paraschivoiu at the same solidity. Howerver, I don't understand why the CP keeps increasing after TSR = 3 in the SSTTM calculations. The value at TSR = 6 is even above the Betz limit... I have already checked the convergence of the solution (grid, time step, size of the calculation domain). I use a maximum y+ value of 1 on the blades. I have run calculations on a static hydrofoil with both SST and SSTTM models and they were in very good agreement with experimental data. The turbulence intensity is very low just before the turbine (0.1%) but I don't see how it could have such a big effect on the CP. I have read in the paper of Menter (AIAA Journal, 2009) that one of the limitations of the SSTTM model is that "the transition correlations are formulated nonGalilean invariant". I don't know if it has been corrected in the version I use (CFX v14.0) and this point is not mentioned in the modeling guide or anywhere else in the CFX documentation. If my understanding is good, this should only affect the turbulence intensity calculated as an input of the transition model. I don't understand how it could result in a factor 2 for the CP calculated at TSR = 6. Do you have any idea regarding this problem ? Thanks for your help! 
1 Attachment(s)

Is the model tightly converged?
You would have to look at the detail of the flow over the blades. There is probably something weird going on  maybe the flow separates in one model and does not in the other, maybe the laminar separation bubble of the transition model is doing weird things. And are you modelling it steady state or transient, or something like frozen rotor? 
Thanks for your answer Glenn!
I am modeling it transient. Convergence is good: residuals (rms) are under 10^(5) and the calculation is run during 6 to 10 turbine revolutions so that the power coefficient reaches a constant value at the end of the calculation. A plot of the pressure coefficient shows that the transition model leads to a slightly higher lift and a plot of the friction coefficient shows that the viscous drag is divided by 2 with the transition model. The lift increase and the drag decrease both tend to an increase of the driving force, so the increase of the power coefficient is rather logical. However such a big increase is weird I think... 
What Re is this running at?
Most wind turbines run fairly low Re numbers, resulting in significant laminar flow regions before transition to turbulence. This means a transition model is to be expected to result in significantly lower drag than a full turbulence model. So your report that the transition model is lower drag is expected, and the drag reduction could be significant (like 50%). This probably explains the difference you are seeing. 
In my case, at TSR = 6 the Reynolds (based on the chord and the relative velocity) varies between 1.5 x 10^6 and 2.1 x 10^6 during one revolution.
In the calculation, transition generally takes place at x/c = 0.3 on the suction side and very close to the trailing edge on the pressure side, which explains the decrease of the viscous drag as you say. However, transition location should probably be located upstream of these points, I think, and this behavior may be due to the low tubulence intensity (0.1%). I am currently running the same case with a higher turbulence intensity to assess its effect on the power coefficient. 
All times are GMT 4. The time now is 09:26. 