|
[Sponsors] |
Shape optimisation with constraint : NACA0012 |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
New Member
Join Date: Feb 2023
Location: MAROCCO
Posts: 7
Rep Power: 3 ![]() |
Hi everebody,
I'm an absolut beginner in SU2, i have to optimise the shape of NACA0012 Airfoil. The probleme is the following : - Objective : Cd(drag) minimisation. - Constrainte : y>ybasline (with 0<x<1). I've seen on NASA's website than ybasline= +- 0.6*[0.2969*sqrt(x) - 0.1260*x - 0.3516*x2 + 0.2843*x3 - 0.1015*x4] The problem looks very easy but i have just no ideas of what i'have to code. I know than i need to complete the following code : OPT_CONSTRAINT = NONE May you help me please ?! Thank's a lot for your help. OCTOGON |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Senior Member
bigfoot
Join Date: Dec 2011
Location: Netherlands
Posts: 596
Rep Power: 17 ![]() |
Did you already try the tutorials on unconstrained and constrained optimization? They are a good starting point:
https://su2code.github.io/tutorials/home/ What is the use of the design constraint? Is a minimum thickness as constraint not sufficient for you? |
|
![]() |
![]() |
![]() |
![]() |
#3 | |
New Member
Join Date: Feb 2023
Location: MAROCCO
Posts: 7
Rep Power: 3 ![]() |
Quote:
Yes i have already check this two tutorials. For unconstrained optimization problem, i succeed easly, i also succeed for an easy contrained optimization probleme (simply THINKNESS>XXX for example). Ok i'm not an absolut biginner as well BUT a biginner for sur ![]() But for my own problem i have no idea. At first, i would like to "cut" the profil in several small profil like on tuto with "ONERAM6", but i don't know if it's answer well my problem. Maybe it's the only possibility, maybe not i don't know. There is no use, it's an academic project ! |
||
![]() |
![]() |
![]() |
![]() |
#4 |
New Member
Join Date: Feb 2023
Location: MAROCCO
Posts: 7
Rep Power: 3 ![]() |
And the approach with THIKNESS in several x (0, 0.2, 0.4 etc.. , 1) is only coherent for symetric design so i can do better i think
|
|
![]() |
![]() |
![]() |
![]() |
#5 |
Senior Member
bigfoot
Join Date: Dec 2011
Location: Netherlands
Posts: 596
Rep Power: 17 ![]() |
I don't understand your constraint, you want a new airfoil design that is always thicker than the 0012 airfoil, so the 0012 fits completely inside your new design? And the cord length is fixed at 1?
|
|
![]() |
![]() |
![]() |
![]() |
#6 | |
New Member
Join Date: Feb 2023
Location: MAROCCO
Posts: 7
Rep Power: 3 ![]() |
Quote:
And i need to have, for all x (airfoil lenght) y>ybaseline. I have no constrainte on cord i think. So it's true than the constraint previously mentionned is equivalent to : NACA0012 fits completely inside the new design. |
||
![]() |
![]() |
![]() |
![]() |
#7 |
New Member
Join Date: Feb 2023
Location: MAROCCO
Posts: 7
Rep Power: 3 ![]() |
But it's not only thickness because we can imagine a new airfoil design with exactly the same thickness which doesn't respect the constraint y>ybaseline (if the airfoil bends upwards with the same thickness)
|
|
![]() |
![]() |
![]() |
![]() |
#8 |
New Member
Abhijith
Join Date: Nov 2020
Location: United Kingdom
Posts: 19
Rep Power: 5 ![]() |
For NACA 0012 airfoil since the thickness of baseline is 0.12, use the following in the constraint to have a optimized shape with thickness more than or equivalent to baseline thickness:
Code:
OPT_OBJECTIVE= DRAG * 0.001 OPT_CONSTRAINT= (AIRFOIL_THICKNESS = 0.12) * 0.001; (AIRFOIL_THICKNESS > 0.12) * 0.001 It is better to use the full geometry of NACA 0012 for this case. For symmetric case refer this article for insights- https://www.researchgate.net/publica...lems_Using_SU2 Hope it helps Best Abhijith Moni |
|
![]() |
![]() |
![]() |
![]() |
#9 | |
New Member
Join Date: Feb 2023
Location: MAROCCO
Posts: 7
Rep Power: 3 ![]() |
Quote:
i don't understand what it the "Scale" !? |
||
![]() |
![]() |
![]() |
![]() |
#10 |
New Member
Abhijith
Join Date: Nov 2020
Location: United Kingdom
Posts: 19
Rep Power: 5 ![]() |
SLSQP optimizer is used by SU2 for design optimisation in this case. In its initial phase of looking for the local optimum, the optimizer perturbs the design. In order to ensure that this first step is not too big, we often scale the objective function (for example 0.001). The value of this scaling factor needs to be changed for various objective functions; in general, there is no certain thumb rule for choosing this value; rather, a little trial and error is required based on the problem/requirement. Read about Adjoint solver and optimizers available within SU2 documentation for more details.
Best Abhijith Moni |
|
![]() |
![]() |
![]() |
![]() |
#11 |
New Member
Join Date: Feb 2023
Location: MAROCCO
Posts: 7
Rep Power: 3 ![]() |
Thank's.
Ok i have an other idea. What about puting OPT_BOUND_LOWER from -0.1 to 0 ?! That shoud allow to only inscrease the yline. In addition to OPT_CONTRAINT previously viewed. |
|
![]() |
![]() |
![]() |
![]() |
#12 |
New Member
Rafael Gomes Batista
Join Date: Apr 2023
Posts: 12
Rep Power: 3 ![]() |
Do you have any updates about this type of optimization?
Or have you had any success with this optimization? By the way, is this an optimization benchmark from the ADODG group? Thank you. Best Regards, Rafael Batista |
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Shape optimisation not converging | elninoft9 | SU2 Shape Design | 3 | January 25, 2023 09:03 |
3D RANS based shape optimisation mesh deformation | nikki_vlk | SU2 Shape Design | 0 | May 30, 2022 04:36 |
Spline generation for Adjoint Based Shape Optimisation | radioactivity | OpenFOAM | 1 | November 1, 2021 00:11 |
Errors when running Shape Optimization Tutorial 1 - NACA0012 | northfly | SU2 | 7 | February 14, 2019 03:46 |
Segmentation fault in tutorial NACA0012 shape optimization | lazzini | SU2 Shape Design | 3 | March 16, 2018 04:41 |