[yugang881]

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!

[yugang881]

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!

[JPBLourenco]

Originally Posted by yugang881

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

[yugang881]

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