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

Negative cl,cd values (rhoSimpleFoam)

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 30, 2023, 06:04
Post Negative cl,cd values (rhoSimpleFoam)
  #1
New Member
 
Rohan
Join Date: May 2023
Posts: 20
Rep Power: 3
Rohanbotta is on a distinguished road
Greetings,

im trying to solve for a finite wing of Naca4415 airfoil data using rhoSimpleFoam. the initial boundary conditions are 0.4M (velocity), and 0 angle of attack. The turbulence model used in k-omega and the flow is compressible. im getting negative coefficients of lift and drag even though i have set the velocity, lift and drag dir based on the coordinate system of the geometry. Tried various things yet nothing works, not sure if im missing something.
kindly help me with this,
regards,
Rohan
Rohanbotta is offline   Reply With Quote

Old   May 31, 2023, 13:48
Question Can anyone help me figure out the issue?
  #2
New Member
 
Rohan
Join Date: May 2023
Posts: 20
Rep Power: 3
Rohanbotta is on a distinguished road
This is the force coeffs code
forces
{
type forceCoeffs;
libs ("libforces.so");
writeControl timeStep;
writeInterval 1;

patches ("wingbottom");
patches ("wingtip");
patches ("wingtop");
rhoInf 1.2;
liftDir (0 1 0);
dragDir (1 0 0);
CofR (0 0 0); // Coefficient of rotation
pitchAxis (0 0 1);
magUInf 136.119;
lRef 1; // Ref length
Aref 1; // Ref area
}
Rohanbotta is offline   Reply With Quote

Old   June 1, 2023, 02:56
Default
  #3
Senior Member
 
Yann
Join Date: Apr 2012
Location: France
Posts: 1,110
Rep Power: 26
Yann will become famous soon enough
Not sure what is going on, but your patches definition is wrong:

Code:
patches ("wingbottom");
patches ("wingtip");
patches ("wingtop");
I'm pretty sure your coefficients are only calculated on wingtop. If you want to include the 3 patches in the coeff calculation it should be all defined at the same place:

Code:
patches (wingbottom wingtip wintop);
Hope this helps,
Yann
Yann is offline   Reply With Quote

Old   June 1, 2023, 03:06
Thumbs up
  #4
New Member
 
Rohan
Join Date: May 2023
Posts: 20
Rep Power: 3
Rohanbotta is on a distinguished road
ohhh, i was thinking lift usually acts on bottom surface of the wing so i need to include wingbottom patch as well. i will try the way you have shown to include all surfaces and check.
Thank you very much,
Rohan
Rohanbotta is offline   Reply With Quote

Old   June 1, 2023, 03:12
Default
  #5
Senior Member
 
Yann
Join Date: Apr 2012
Location: France
Posts: 1,110
Rep Power: 26
Yann will become famous soon enough
Quote:
Originally Posted by Yann View Post
I'm pretty sure your coefficients are only calculated on wingtop.
I think there is a misunderstanding: I mean on your forceCoeffs definition, the last "patches" line override the previous ones so the way your function is defined, it only uses wingtop in the calculation instead of the whole wing (wingtop+wingbottom+wingtip)

Yann
Yann is offline   Reply With Quote

Old   June 1, 2023, 03:26
Question
  #6
New Member
 
Rohan
Join Date: May 2023
Posts: 20
Rep Power: 3
Rohanbotta is on a distinguished road
yes, i made a mistake in defining patches, so the alternative way to include them which you have suggested, [patches (wingbottom wingtip wintop);] will this work for my case?
Thank you,
Rohan
Rohanbotta is offline   Reply With Quote

Old   June 1, 2023, 04:44
Default
  #7
Senior Member
 
Yann
Join Date: Apr 2012
Location: France
Posts: 1,110
Rep Power: 26
Yann will become famous soon enough
Try and see
Check patches names, liftDir and dragDir.

Yann
Yann is offline   Reply With Quote

Old   June 1, 2023, 05:02
Exclamation
  #8
New Member
 
Rohan
Join Date: May 2023
Posts: 20
Rep Power: 3
Rohanbotta is on a distinguished road
i have tried that statement for patches, the negative Cl, Cd is resolved but the Cl, Cd values are quite high like Cl is ranging around 4.83e+01, Cd around 8,71e-01. i have also tried defining only wingtop as a patch but the cl cd values are all zeroes in this case.
Rohanbotta is offline   Reply With Quote

Old   June 1, 2023, 05:21
Default
  #9
Senior Member
 
Yann
Join Date: Apr 2012
Location: France
Posts: 1,110
Rep Power: 26
Yann will become famous soon enough
Does your reference length and area make sense?

About patches, there is no doubt: if you want to compute the coefficient of your profile you need to include the whole profile.
Yann is offline   Reply With Quote

Old   June 1, 2023, 05:26
Post
  #10
New Member
 
Rohan
Join Date: May 2023
Posts: 20
Rep Power: 3
Rohanbotta is on a distinguished road
yes the chord length is 0.046 and the ref area value was calculated based on chord length and span. i will recheck it again anyhow. What units does openfoam consider, for lref and aref, i think there might be a problem there. So i have calculated lref and aref like this (lref = chord length, aref = chord length*span) is this correct or is there any other way to calculate them
Thank you
Rohanbotta is offline   Reply With Quote

Old   June 1, 2023, 05:30
Default
  #11
Senior Member
 
Yann
Join Date: Apr 2012
Location: France
Posts: 1,110
Rep Power: 26
Yann will become famous soon enough
Both lRef and Aref were set to 1 in your first post, I guess you fixed it afterward?
Yann is offline   Reply With Quote

Old   June 1, 2023, 05:34
Post
  #12
New Member
 
Rohan
Join Date: May 2023
Posts: 20
Rep Power: 3
Rohanbotta is on a distinguished road
ohh yes i fixed it later, so lref = chord length and aref = chordlength*span, is this the right way to set the values or is there any other way to find them. and i assumed openfoam takes units as m or m^2 for these values, is this where th problem might be?
Rohanbotta is offline   Reply With Quote

Old   June 1, 2023, 06:01
Default
  #13
Senior Member
 
Yann
Join Date: Apr 2012
Location: France
Posts: 1,110
Rep Power: 26
Yann will become famous soon enough
Sounds good yes. Have you check your simulation results ? Did it converge ? Velocity and pressure fields make sense?
Yann is offline   Reply With Quote

Old   June 1, 2023, 10:13
Default
  #14
New Member
 
Rohan
Join Date: May 2023
Posts: 20
Rep Power: 3
Rohanbotta is on a distinguished road
Yes the velocity fields and pressure contours made sense however I'm not sure about the values, as I said cl and cd values are quite high, so I'm not sure about it, I'll try a different case taking a new wing and calculate lref Aref newly and try
Rohanbotta is offline   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
Preventing negative scalar (UDS) values?? Matthew Brannock FLUENT 5 September 1, 2020 12:32
Negative Drag Values for only some surfaces?? Please help Jordan1997 Main CFD Forum 0 May 2, 2020 06:39
Cannot achieve negative wall shear stress values in CFD Post (ANSYS) SAtkins CFX 3 August 3, 2016 14:27
Negative values for minimum alphas Taataa OpenFOAM Running, Solving & CFD 0 June 28, 2016 04:04
Mass flow rate boundary condition with negative values ashtonJ CFX 3 November 26, 2014 05:21


All times are GMT -4. The time now is 18:39.