# shape optimization problem

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

November 3, 2017, 04:06
shape optimization problem
#1
New Member

JAMES
Join Date: Dec 2016
Posts: 21
Rep Power: 8
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
% --------------------- 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
 541FFDBOX.png (32.3 KB, 29 views) 541FFDBOX1.png (11.5 KB, 22 views)
Attached Files
 turb_BWBOPT.txt (13.9 KB, 5 views) log_Adjoint.txt (126.3 KB, 3 views)

 December 13, 2017, 02:37 #2 New Member   JAMES Join Date: Dec 2016 Posts: 21 Rep Power: 8 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!

December 13, 2017, 07:06
#3
Member

João Lourenço
Join Date: Jan 2017
Posts: 32
Rep Power: 8
Quote:
 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

 December 14, 2017, 05:26 #4 New Member   JAMES Join Date: Dec 2016 Posts: 21 Rep Power: 8 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