# problem with torque calc

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 May 31, 2009, 06:46 problem with torque calc #1 Member   Hai Yu Join Date: Mar 2009 Location: Harbin Posts: 67 Rep Power: 10 Sponsored Links Dear foamers, I am now calculating the torque of a turbine driven by wind power, with turbDymFoam. at a higher rotating speed, the torque shows a perfect fit with the experiementation data. while, at lower rotating speed, the torque coeff is 20% lower than experimentation data. Can you give some hints for this? cheers.

 May 31, 2009, 08:06 #2 Senior Member   Henrik Rusche Join Date: Mar 2009 Location: Braunschweig, Niedersachsen, Germany Posts: 277 Rep Power: 11 Dear Yuhai, the information you are offering to solve the problem is far from exhaustive. Could this be problem related to the turbulence modelling you use? Henrik

May 31, 2009, 09:41
#3
Member

Hai Yu
Join Date: Mar 2009
Location: Harbin
Posts: 67
Rep Power: 10
Quote:
 Originally Posted by henrik Dear Yuhai, the information you are offering to solve the problem is far from exhaustive. Could this be problem related to the turbulence modelling you use? Henrik
oh, yes, sorry, I am lack of information.
I am now using the standard k-epsilon with values of k and e, that work fine in Fluent, What I mean is that I used same model and k/e value in OpenFOAM and Fluent, but at lower rotating speed, they give a difference of 20%, while quite alike at higher rotating speed.
I have also tried RNG and SST, but the result is unacceptable for my current k/e/w value. For RNG and SST is not conducted in Fluent, I can't make a compare.

hier is some setting in fvScheme:

ddtSchemes
{
default Euler;
}
{
default leastSquares;
}
divSchemes
{
default none;
div(phi,U) Gauss upwind;
div(phi,k) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div(phi,omega) Gauss upwind;
div(phi,R) Gauss upwind;
div(R) Gauss linear;
div(phi,nuTilda) Gauss upwind;
}
laplacianSchemes
{
default none;
laplacian(nuEff,U) Gauss linear corrected;
laplacian((1|A(U)),p) Gauss linear corrected;
laplacian(DkEff,k) Gauss linear corrected;
laplacian(DepsilonEff,epsilon) Gauss linear corrected;
laplacian(DREff,R) Gauss linear corrected;
laplacian(DnuTildaEff,nuTilda) Gauss linear corrected;
laplacian(DomegaEff,omega) Gauss linear corrected;
laplacian((rho*(1|A(U))),p) Gauss linear corrected;
laplacian(alphaEff,h) Gauss linear corrected;
laplacian(rAU,p) Gauss linear corrected;
}
interpolationSchemes
{
default linear;
interpolate(U) linear;
}
{
default corrected;
}
fluxRequired
{
default no;
p;
}

the air flow is from left side to the turbine. wind speed 10m/s.
while I got fine result at angular velocity = 14, but not at av=8.
I am now wondering whether should I calculate different Re and k/e for different angular velocity?

Best Regards.
Attached Images
 Clip_3.jpg (44.4 KB, 16 views)

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Patrick CFX 6 September 19, 2013 05:51 Wouter Fluent UDF and Scheme Programming 6 June 6, 2012 04:43 Gianni FLUENT 0 April 5, 2008 10:33 Se-Hee CFX 2 June 10, 2007 06:29 ParodDav CFX 5 April 29, 2007 19:13