CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > SU2 Shape Design

Problem with restart solution in shape_optimization.py

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 23, 2013, 10:45
Default
  #21
Member
 
Roberto Pieri
Join Date: Feb 2012
Location: Milan
Posts: 57
Rep Power: 4
robyTKD is on a distinguished road
Sorry, I didn't read the post on the forum regarding no-slip conditions. Now it works.

Thank you.
robyTKD is offline   Reply With Quote

Old   May 29, 2013, 10:26
Default
  #22
Member
 
Roberto Pieri
Join Date: Feb 2012
Location: Milan
Posts: 57
Rep Power: 4
robyTKD is on a distinguished road
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 )
I think that gradients are not ok.
Do you have some reference to compare with my shape sensitivities?

Thank you very much.
robyTKD is offline   Reply With Quote

Reply

Tags
optimization, restart solution, su2

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
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 06:47
Solution of Lid driven Cavity problem using SIMPLE algorithm sasank Main CFD Forum 2 May 20, 2012 10:22
Problem with restart henryss STAR-CD 2 October 25, 2009 09:29
How to get an accurate CFD solution to aeroheating problem for hypersonic flows? Leining D Main CFD Forum 5 April 14, 2009 12:26
Numerical solution to the rotating disk problem? johny Main CFD Forum 7 September 5, 2005 06:53


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