AeolusAero / Aeolus-AeroSketchPad

Aeolus ASP is an aerodynamic analysis and optimization tool for aircraft and propellers
https://www.aeolus-aero.com
2 stars 0 forks source link

How to get fixed AOA and avoid non-convergence #4

Closed j2l closed 1 year ago

j2l commented 2 years ago

Hello, I'm learning Aeolus basic functions and validate simple profiles with Xfoil, Javafoil and Nasa FoilSim website https://www.grc.nasa.gov/WWW/K-12/airplane/foil3.html : An ellipse 20% thickness with a 2.5m span and a 0.1m chord @ 3.94444 m/s.

FoilSim (2D-3D), Xfoil (2D) and Javafoil (2D) roughly get the same result for Cl and Cd: image image image

Entered Reynolds number is from Aeolus (freshwater). I imported the same profile in Aeolus (first time with a 110 points profile, and the second one with a 48 points profile)

First, I can't tell Aeolus to stick to AoA 6.3: Despite choosing AoA1, 2 and 3 as 6.3, it keeps returning a resulting AoA of 9.997° Even with default -3,5,10° it doesn't get the right AoA around 6.3° but sticks to 9.997° There's a warning saying it has bad convergence despite smoothing discretization, both with 111 and 48 points. I find it pretty slower than for other profiles (10 cores result in 25s) even for a simple 48 points dat profile.

How to only compute one AoA? Any idea why it's not converging and slower than usual? Maybe it's related? Do you have hints to validate this kind of profile?

I attach my Aeolus file and both ellipses dat files ellipses.zip

Thank you very much and have a great week.

Phil

UweSchuster commented 2 years ago

Hi Phil,

Thank you for your message. I am currently travelling and have no access to a computer but here are some notes to help: Aeolus does compute the AOA to provide the target design lift. The 3x AOA you can enter are just the support points for the interpolation. In the Cl,Cd results table you will however find the cl cd values for these support AOAs as well as for the final AOA. For discretization, i propose at least 48 panels in the circumferential direction and an aspect ratio of at least 0.3 for all strips. Hope this helps!

Best regards,

Uwe

On Apr 24, 2022 at 7:58 PM, <Philippe Manzano @.***)> wrote:

Hello, I'm learning Aeolus basic functions and validate simple profiles with Xfoil, Javafoil and Nasa FoilSim website https://www.grc.nasa.gov/WWW/K-12/airplane/foil3.html : An ellipse 20% thickness with a 2.5m span and a 0.1m chord @ 3.94444 m/s.

FoilSim (2D-3D), Xfoil (2D) and Javafoil (2D) roughly get the same result for Cl and Cd: (https://user-images.githubusercontent.com/65325/164988758-cfde1362-737e-4d82-991b-02201c5f4f09.png) (https://user-images.githubusercontent.com/65325/164988773-446afe7f-35d9-46b3-b590-45adebf42b92.png) (https://user-images.githubusercontent.com/65325/164989650-a27b4429-5852-4443-8ebd-5b6cbde292cc.png)

Entered Reynolds number is from Aeolus (freshwater). I imported the same profile in Aeolus (first time with a 110 points profile, and the second one with a 48 points profile)

First, I can't tell Aeolus to stick to AoA 6.3: Despite choosing AoA1, 2 and 3 as 6.3, it keeps returning a resulting AoA of 9.997° Even with default -3,5,10° it doesn't get the right AoA around 6.3° but sticks to 9.997° There's a warning saying it has bad convergence despite smoothing discretization, both with 111 and 48 points. I find it pretty slower than for other profiles (10 cores result in 25s) even for a simple 48 points dat profile.

How to only compute one AoA? Any idea why it's not converging and slower than usual? Maybe it's related? Do you have hints to validate this kind of profile?

I attach my Aeolus file and both ellipses dat files ellipses.zip (https://github.com/AeolusAero/Aeolus-AeroSketchPad/files/8550224/ellipses.zip)

Thank you very much and have a great week.

Phil

— Reply to this email directly, view it on GitHub (https://github.com/AeolusAero/Aeolus-AeroSketchPad/issues/4), or unsubscribe (https://github.com/notifications/unsubscribe-auth/ALJIKHSEM64XSXWFYIHABOLVGWDUPANCNFSM5UGSASBQ). You are receiving this because you are subscribed to this thread.Message ID: @.***>

j2l commented 2 years ago

Thank you very much for replying despite you're traveling!

Here's what I have in the attached file above (BTW, xml extension is not added to saved file by default): image Aspect ratio is far beyond 0.3

For AoA, if I keep default settings: image

Computation and "best" AoA is wrong: image

as it throws wrong results: image

I agree with this message: image

but I don't know how to improve it. Discretization seems fine.

When you have time, if you can test with my attached files so I can fix the missing thing to correctly do tests, thank you :) foil_ellipse_eau.xml.zip

j2l commented 2 years ago

bump :smiling_face_with_tear:

j2l commented 2 years ago

When you have time, if you can test with my attached files so I can fix the missing thing to correctly do tests, thank you :)

Any news?

UweSchuster commented 2 years ago

I am currently busy and will not be able to look into it before August

On Jun 5, 2022 at 3:41 PM, <Philippe Manzano @.***)> wrote:

When you have time, if you can test with my attached files so I can fix the missing thing to correctly do tests, thank you :)

Any news?

— Reply to this email directly, view it on GitHub (https://github.com/AeolusAero/Aeolus-AeroSketchPad/issues/4#issuecomment-1146816919), or unsubscribe (https://github.com/notifications/unsubscribe-auth/ALJIKHQQA6XHR5YJBR4ORULVNS4BRANCNFSM5UGSASBQ). You are receiving this because you commented.Message ID: @.***>

j2l commented 2 years ago

I understand. I'll wait until August.

UweSchuster commented 2 years ago

Hi Phil!

Sorry for the long waiting time. I am now reviewing your files and can reproduce the issue. First, I tweaked your airfoil files by adding the last [1.0 0.0] point to make sure the trailing edge is closed and the airfoil is imported properly.

I think it is because of the Airfoil itself and not the discretization. Changing to NACA0012 works. Today, AeolusASP is based on potential flow which requires airfoils with sharp trailing edges (to apply the Kutta condition). The elliptical airfoil is therefore inherently challenging for potential flow solvers.

UweSchuster commented 2 years ago

Here is a working example, which is not so far off though: Ellipse_Test.zip

reading approximately Cl = 0.74 at 6.3 deg.

image

j2l commented 2 years ago

Thank you @UweSchuster Indeed, I had the same issue (Kutta condition) with another solver, so it makes more sense to me now :smile: Looks like your file was done with 4.22 and latest downloadable version is 4.21.

UweSchuster commented 1 year ago

To summarize here: