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

shape optimization problem

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 3, 2017, 05:06
Question shape optimization problem
  #1
New Member
 
JAMES
Join Date: Dec 2016
Posts: 21
Rep Power: 5
yugang881 is on a distinguished road
Dear friends:
I meet 2 problems when I use SU2 to do a BWB shape optimization.
1. The first one is the warning in the log_adjoint.out file, it shows at the start as this:
WARNING: The adjoint solver should use a non-dimensional flow solution.
The flow solution file is also calculated by SU2 using Rans. Isn’t it a non-dimensional flow solution?
2. The second one is the gradient of Cd and Cl are all 0. I change the scale factor form 0.001 to 0.1 and to 1. The results are all the same 0 and shape optimization ended at the first step.
The configureation is as :
//////////////////////////////////////////////////
% -------------------- FREE-FORM DEFORMATION PARAMETERS -----------------------%
%
% Tolerance of the Free-Form Deformation point inversion
FFD_TOLERANCE= 1E-10
%
% Maximum number of iterations in the Free-Form Deformation point inversion
FFD_ITERATIONS= 500
%
% FFD box definition: 3D case (FFD_BoxTag, X1, Y1, Z1, X2, Y2, Z2, X3, Y3, Z3, X4, Y4, Z4,
% X5, Y5, Z5, X6, Y6, Z6, X7, Y7, Z7, X8, Y8, Z8)
% 2D case (FFD_BoxTag, X1, Y1, 0.0, X2, Y2, 0.0, X3, Y3, 0.0, X4, Y4, 0.0,
% 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
FFD_DEFINITION=(0,-2.24858056360336,1.5643188577435723e-13,-4.373899154034036, 47.23476648517159, 1.5643188577435723e-13, -4.373899154034036,40.70375878076216, 32.56444337346593, -4.373899154034036,16.50252511766261, 32.56444337346593, -4.373899154034036,-2.24858056360336, 1.5643188577435723e-13, 2.438575072894233,47.23476648517159, 1.5643188577435723e-13, 2.438575072894233,40.70375878076216, 32.56444337346593, 2.438575072894233,16.50252511766261, 32.56444337346593, 2.438575072894233)
%
% FFD box degree: 3D case (x_degree, y_degree, z_degree)
% 2D case (x_degree, y_degree, 0)
FFD_DEGREE= (5, 4, 1)
%
% Surface continuity at the intersection with the FFD (1ST_DERIVATIVE, 2ND_DERIVATIVE)
FFD_CONTINUITY= 2ND_DERIVATIVE

% ------------------------ GRID DEFORMATION PARAMETERS ------------------------%
%
% Kind of deformation (FFD_SETTING, FFD_CONTROL_POINT_2D, FFD_CAMBER_2D, FFD_THICKNESS_2D,
DV_KIND= FFD_CONTROL_POINT
%
% Marker of the surface in which we are going apply the shape deformation
DV_MARKER= ( WINGBODY )
%
% Parameters of the shape deformation

%
% New value of the shape deformation
DV_VALUE= 0.1
%
% Visualize the deformation (NO, YES)
VISUALIZE_DEFORMATION= NO

% --------------------------- CONVERGENCE PARAMETERS --------------------------%
%
% Convergence criteria (CAUCHY, RESIDUAL)
%
CONV_CRITERIA= CAUCHY
%
% Residual reduction (order of magnitude with respect to the initial value)
RESIDUAL_REDUCTION= 3
%
% Min value of the residual (log10 of the residual)
RESIDUAL_MINVAL= -6
%
% Start convergence criteria at iteration number
STARTCONV_ITER= 1
%
% Number of elements to apply the criteria
CAUCHY_ELEMS= 100
%
% Epsilon to control the series convergence
CAUCHY_EPS= 1E-5
%
% Function to apply the criteria (LIFT, DRAG, NEARFIELD_PRESS, SENS_GEOMETRY,
% SENS_MACH, DELTA_LIFT, DELTA_DRAG)
CAUCHY_FUNC_FLOW= DRAG
CAUCHY_FUNC_ADJFLOW= SENS_GEOMETRY
% --------------------- OPTIMAL SHAPE DESIGN DEFINITION -----------------------%
% Optimization objective function with optional scaling factor
% ex= Objective * Scale
OPT_OBJECTIVE= DRAG * 0.1

% Optimization constraint functions with scaling factors, separated by semicolons
% ex= (Objective = Value ) * Scale, use '>','<','='
OPT_CONSTRAINT= (LIFT > 0.2238) * 0.1

% List of design variables (Design variables are separated by semicolons)
% - FFD_CONTROL_POINT ( 7, Scale | Mark. List | FFD_BoxTag, i_Ind, j_Ind, k_Ind, x_Mov, y_Mov, z_Mov )(from up x=0,1,2,3,4surface to down)
DEFINITION_DV= ( 7, 1.0 | wing | 0, 0, 0, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 1, 0, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 2, 0, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 3, 0, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 4, 0, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 0, 1, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 1, 1, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 2, 1, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 3, 1, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 4, 1, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 0, 2, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 1, 2, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 2, 2, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 3, 2, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 4, 2, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 0, 3, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 1, 3, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 2, 3, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 3, 3, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 4, 3, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 0, 4, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 1, 4, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 2, 4, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 3, 4, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 4, 4, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 0, 0, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 1, 0, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 2, 0, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 3, 0, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 4, 0, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 0, 1, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 1, 1, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 2, 1, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 3, 1, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 4, 1, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 0, 2, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 1, 2, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 2, 2, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 3, 2, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 4, 2, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 0, 3, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 1, 3, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 2, 3, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 3, 3, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 4, 3, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 0, 4, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 1, 4, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 2, 4, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 3, 4, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 4, 4, 0, 0.0, 0.0, 1.0 )
////////////////////////////////////////////////////
The ffd box and shape are also attached.

Could anybody kindly show me how to fix this problem? Many thanks!
Attached Images
File Type: png 541FFDBOX.png (32.3 KB, 24 views)
File Type: png 541FFDBOX1.png (11.5 KB, 18 views)
Attached Files
File Type: txt turb_BWBOPT.txt (13.9 KB, 4 views)
File Type: txt log_Adjoint.txt (126.3 KB, 2 views)
yugang881 is offline   Reply With Quote

Old   December 13, 2017, 03:37
Default
  #2
New Member
 
JAMES
Join Date: Dec 2016
Posts: 21
Rep Power: 5
yugang881 is on a distinguished road
I changed the FFd box and fixed the gradient =0 problem. Now it could caculate the gradient. But it still stopped at the second step because of the SU2_DEF runs wrong. It seems that the SU2_DEF couldn't deform the mesh. And the error message is as:
================================================== =================================
= BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
= PID 9853 RUNNING AT wk-pc
= EXIT CODE: 9
= CLEANING UP REMAINING PROCESSES
= YOU CAN IGNORE THE BELOW CLEANUP MESSAGES
================================================== =================================
YOUR APPLICATION TERMINATED WITH THE EXIT STRING: Killed (signal 9)
This typically refers to a problem with your application.
Please see the FAQ page for debugging suggestions

I checked the FAQ page but cant find the debugging suggestion. Could anybody kindly show what can i do for next step? Thanks very much!
yugang881 is offline   Reply With Quote

Old   December 13, 2017, 08:06
Question
  #3
Member
 
João Lourenço
Join Date: Jan 2017
Posts: 32
Rep Power: 5
JPBLourenco is on a distinguished road
Quote:
Originally Posted by yugang881 View Post
Dear friends:
I meet 2 problems when I use SU2 to do a BWB shape optimization.
1. The first one is the warning in the log_adjoint.out file, it shows at the start as this:
WARNING: The adjoint solver should use a non-dimensional flow solution.
The flow solution file is also calculated by SU2 using Rans. Isn’t it a non-dimensional flow solution?
2. The second one is the gradient of Cd and Cl are all 0. I change the scale factor form 0.001 to 0.1 and to 1. The results are all the same 0 and shape optimization ended at the first step.
The configureation is as :
//////////////////////////////////////////////////
% -------------------- FREE-FORM DEFORMATION PARAMETERS -----------------------%
%
% Tolerance of the Free-Form Deformation point inversion
FFD_TOLERANCE= 1E-10
%
% Maximum number of iterations in the Free-Form Deformation point inversion
FFD_ITERATIONS= 500
%
% FFD box definition: 3D case (FFD_BoxTag, X1, Y1, Z1, X2, Y2, Z2, X3, Y3, Z3, X4, Y4, Z4,
% X5, Y5, Z5, X6, Y6, Z6, X7, Y7, Z7, X8, Y8, Z8)
% 2D case (FFD_BoxTag, X1, Y1, 0.0, X2, Y2, 0.0, X3, Y3, 0.0, X4, Y4, 0.0,
% 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
FFD_DEFINITION=(0,-2.24858056360336,1.5643188577435723e-13,-4.373899154034036, 47.23476648517159, 1.5643188577435723e-13, -4.373899154034036,40.70375878076216, 32.56444337346593, -4.373899154034036,16.50252511766261, 32.56444337346593, -4.373899154034036,-2.24858056360336, 1.5643188577435723e-13, 2.438575072894233,47.23476648517159, 1.5643188577435723e-13, 2.438575072894233,40.70375878076216, 32.56444337346593, 2.438575072894233,16.50252511766261, 32.56444337346593, 2.438575072894233)
%
% FFD box degree: 3D case (x_degree, y_degree, z_degree)
% 2D case (x_degree, y_degree, 0)
FFD_DEGREE= (5, 4, 1)
%
% Surface continuity at the intersection with the FFD (1ST_DERIVATIVE, 2ND_DERIVATIVE)
FFD_CONTINUITY= 2ND_DERIVATIVE

% ------------------------ GRID DEFORMATION PARAMETERS ------------------------%
%
% Kind of deformation (FFD_SETTING, FFD_CONTROL_POINT_2D, FFD_CAMBER_2D, FFD_THICKNESS_2D,
DV_KIND= FFD_CONTROL_POINT
%
% Marker of the surface in which we are going apply the shape deformation
DV_MARKER= ( WINGBODY )
%
% Parameters of the shape deformation

%
% New value of the shape deformation
DV_VALUE= 0.1
%
% Visualize the deformation (NO, YES)
VISUALIZE_DEFORMATION= NO

% --------------------------- CONVERGENCE PARAMETERS --------------------------%
%
% Convergence criteria (CAUCHY, RESIDUAL)
%
CONV_CRITERIA= CAUCHY
%
% Residual reduction (order of magnitude with respect to the initial value)
RESIDUAL_REDUCTION= 3
%
% Min value of the residual (log10 of the residual)
RESIDUAL_MINVAL= -6
%
% Start convergence criteria at iteration number
STARTCONV_ITER= 1
%
% Number of elements to apply the criteria
CAUCHY_ELEMS= 100
%
% Epsilon to control the series convergence
CAUCHY_EPS= 1E-5
%
% Function to apply the criteria (LIFT, DRAG, NEARFIELD_PRESS, SENS_GEOMETRY,
% SENS_MACH, DELTA_LIFT, DELTA_DRAG)
CAUCHY_FUNC_FLOW= DRAG
CAUCHY_FUNC_ADJFLOW= SENS_GEOMETRY
% --------------------- OPTIMAL SHAPE DESIGN DEFINITION -----------------------%
% Optimization objective function with optional scaling factor
% ex= Objective * Scale
OPT_OBJECTIVE= DRAG * 0.1

% Optimization constraint functions with scaling factors, separated by semicolons
% ex= (Objective = Value ) * Scale, use '>','<','='
OPT_CONSTRAINT= (LIFT > 0.2238) * 0.1

% List of design variables (Design variables are separated by semicolons)
% - FFD_CONTROL_POINT ( 7, Scale | Mark. List | FFD_BoxTag, i_Ind, j_Ind, k_Ind, x_Mov, y_Mov, z_Mov )(from up x=0,1,2,3,4surface to down)
DEFINITION_DV= ( 7, 1.0 | wing | 0, 0, 0, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 1, 0, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 2, 0, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 3, 0, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 4, 0, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 0, 1, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 1, 1, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 2, 1, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 3, 1, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 4, 1, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 0, 2, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 1, 2, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 2, 2, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 3, 2, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 4, 2, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 0, 3, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 1, 3, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 2, 3, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 3, 3, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 4, 3, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 0, 4, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 1, 4, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 2, 4, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 3, 4, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 4, 4, 1, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 0, 0, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 1, 0, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 2, 0, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 3, 0, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 4, 0, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 0, 1, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 1, 1, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 2, 1, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 3, 1, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 4, 1, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 0, 2, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 1, 2, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 2, 2, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 3, 2, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 4, 2, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 0, 3, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 1, 3, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 2, 3, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 3, 3, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 4, 3, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 0, 4, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 1, 4, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 2, 4, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 3, 4, 0, 0.0, 0.0, 1.0 ); ( 7, 1.0 | wing | 0, 4, 4, 0, 0.0, 0.0, 1.0 )
////////////////////////////////////////////////////
The ffd box and shape are also attached.

Could anybody kindly show me how to fix this problem? Many thanks!
Hi, yugang881.

I have your first problem too but in an 2D RANS problem: "WARNING: The adjoint solver should use a non-dimensional flow solution.
The flow solution file is also calculated by SU2 using Rans. Isn’t it a non-dimensional flow solution?"

Did you solve this problem? Can you explain me how?

thanks a lot
João Lourenço
JPBLourenco is offline   Reply With Quote

Old   December 14, 2017, 06:26
Default
  #4
New Member
 
JAMES
Join Date: Dec 2016
Posts: 21
Rep Power: 5
yugang881 is on a distinguished road
Hi João Lourenço:
I checked the testcase example "oneram6 optimization" and found that the this warning is also exist. So i think maybe it is just a warning but don't effect the optimization process. This is only my idea, i am also not very sure.
Best regards
yugang881 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
Shape Optimization: How to export optimized geometry? macRC STAR-CCM+ 1 October 16, 2017 00:38
Problem with running tutorial "Constrained Optimal Shape Design of a Fixed Wing" yugang881 SU2 12 February 3, 2017 08:15
Shape Optimization Problem of Mass Flow Rate ggc1991 SU2 Shape Design 1 December 1, 2016 21:24
[ICEM] Problem Projected Mesh Shape DeanM ANSYS Meshing & Geometry 2 December 18, 2014 03:55
Testcase for shape optimization using RANS simualtion diwakaranant SU2 Shape Design 1 April 18, 2013 14:27


All times are GMT -4. The time now is 01:56.