|
[Sponsors] |
Problem with restart solution in shape_optimization.py |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#21 |
Member
Roberto Pieri
Join Date: Feb 2012
Location: Milan
Posts: 57
Rep Power: 14 ![]() |
Sorry, I didn't read the post on the forum regarding no-slip conditions. Now it works.
Thank you. |
|
![]() |
![]() |
![]() |
![]() |
#22 |
Member
Roberto Pieri
Join Date: Feb 2012
Location: Milan
Posts: 57
Rep Power: 14 ![]() |
I have compiled the 2.0.4 version of the code and I have tested it on my case (NACA0012 turbulent subsonic); direct and adjoint simulations converge, also mesh deformation and optimization framework run without any problem.
When running shape optimization(DRAG minimization with constraint on LIFT), after the first iteration, Cd is not reduced, indeed it is greater than the undeformed case. I post shape optimization configuration for my case: Code:
% --------------------- OPTIMAL SHAPE DESIGN DEFINITION -----------------------% % Objective function (DRAG, LIFT, SIDEFORCE, PRESSURE, MOMENT_X, MOMENT_Y, % MOMENT_Z, EFFICIENCY, CHARGE, SONIC_BOOM) OPT_OBJECTIVE= DRAG*0.01 % % Inequality constraints list separated by comma (DRAG, LIFT, SIDEFORCE, % PRESSURE, MOMENT_X, MOMENT_Y, % MOMENT_Z, EFFICIENCY, SONIC_BOOM) OPT_CONSTRAINT= ( LIFT>1.09 )*0.01 % % Equality constraints list separated by comma (DRAG, LIFT, SIDEFORCE, % PRESSURE, MOMENT_X, MOMENT_Y, % MOMENT_Z, EFFICIENCY, SONIC_BOOM) CONST_EQ= NONE % % Scale equality constraints (separated by comma) CONST_EQ_SCALE= 0.0 % % Value equality constraints list (separated by comma) CONST_EQ_VALUE= 0.0 % % List of design variables (Design variables are separated by semicolons) % From 1 to 99, Geometrycal design variables. % - HICKS_HENNE ( 1, Scale | Mark. List | Lower(0)/Upper(1) side, x_Loc ) % - NACA_4DIGITS ( 4, Scale | Mark. List | 1st digit, 2nd digit, 3rd and 4th digit ) % - ROTATION ( 6, Scale | Mark. List | x_Axis, y_Axis, z_Axis, x_Turn, y_Turn, z_Turn ) % From 100 to 199, Flow solver design variables. % - MACH_NUMBER ( 101, Scale | Markers List ) % - AOA ( 102, Scale | Markers List ) DEFINITION_DV= ( 1, 1.0 | AIRFOIL | 0, 0.05 ); ( 1, 1.0 | AIRFOIL | 0, 0.10 ); ( 1, 1.0 | AIRFOIL | 0, 0.15 ); ( 1, 1.0 | AIRFOIL | 0, 0.20 ); ( 1, 1.0 | AIRFOIL | 0, 0.25 ); ( 1, 1.0 | AIRFOIL | 0, 0.30 ); ( 1, 1.0 | AIRFOIL | 0, 0.35 ); ( 1, 1.0 | AIRFOIL | 0, 0.40 ); ( 1, 1.0 | AIRFOIL | 0, 0.45 ); ( 1, 1.0 | AIRFOIL | 0, 0.50 ); ( 1, 1.0 | AIRFOIL | 0, 0.55 ); ( 1, 1.0 | AIRFOIL | 0, 0.60 ); ( 1, 1.0 | AIRFOIL | 0, 0.65 ); ( 1, 1.0 | AIRFOIL | 0, 0.70 ); ( 1, 1.0 | AIRFOIL | 0, 0.75 ); ( 1, 1.0 | AIRFOIL | 0, 0.80 ); ( 1, 1.0 | AIRFOIL | 0, 0.85 ); ( 1, 1.0 | AIRFOIL | 0, 0.90 ); ( 1, 1.0 | AIRFOIL | 0, 0.95 ); ( 1, 1.0 | AIRFOIL | 1, 0.05 ); ( 1, 1.0 | AIRFOIL | 1, 0.10 ); ( 1, 1.0 | AIRFOIL | 1, 0.15 ); ( 1, 1.0 | AIRFOIL | 1, 0.20 ); ( 1, 1.0 | AIRFOIL | 1, 0.25 ); ( 1, 1.0 | AIRFOIL | 1, 0.30 ); ( 1, 1.0 | AIRFOIL | 1, 0.35 ); ( 1, 1.0 | AIRFOIL | 1, 0.40 ); ( 1, 1.0 | AIRFOIL | 1, 0.45 ); ( 1, 1.0 | AIRFOIL | 1, 0.50 ); ( 1, 1.0 | AIRFOIL | 1, 0.55 ); ( 1, 1.0 | AIRFOIL | 1, 0.60 ); ( 1, 1.0 | AIRFOIL | 1, 0.65 ); ( 1, 1.0 | AIRFOIL | 1, 0.70 ); ( 1, 1.0 | AIRFOIL | 1, 0.75 ); ( 1, 1.0 | AIRFOIL | 1, 0.80 ); ( 1, 1.0 | AIRFOIL | 1, 0.85 ); ( 1, 1.0 | AIRFOIL | 1, 0.90 ); ( 1, 1.0 | AIRFOIL | 1, 0.95 ) Do you have some reference to compare with my shape sensitivities? Thank you very much. |
|
![]() |
![]() |
![]() |
Tags |
optimization, restart solution, su2 |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Problem with SimpleFoam for a solution around an OneraM6 wing | gara1988 | OpenFOAM Running, Solving & CFD | 3 | November 13, 2012 05:47 |
Solution of Lid driven Cavity problem using SIMPLE algorithm | sasank | Main CFD Forum | 2 | May 20, 2012 09:22 |
Problem with restart | henryss | STAR-CD | 2 | October 25, 2009 07:29 |
How to get an accurate CFD solution to aeroheating problem for hypersonic flows? | Leining D | Main CFD Forum | 5 | April 14, 2009 11:26 |
Numerical solution to the rotating disk problem? | johny | Main CFD Forum | 7 | September 5, 2005 05:53 |