Problem with tutorial "Constrained Optimal Shape Design of a Fixed Wing" version 4.2
 

Hello,

Following a previous thread found in https://www.cfd-online.com/Forums/su...ixed-wing.html , I experience similar issues in executing the tutorial.
Apparently, in the previous thread the problem was related with geometrical constraints (for example: MAX_THICKNESS_SEC1), and limited to version 4.3.0.
However, I experience the same issue when using version 4.2 and lift as a constraint.
The settings on the cfg file are the tutorial default ones.
Specifically:
% Optimization objective function with 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.2864) * 0.1; (MAX_THICKNESS_SEC1 > 0.0570) * 0.1; (MAX_THICKNESS_SEC2 > 0.0513) * 0.1; (MAX_THICKNESS_SEC3 > 0.0457) * 0.1; (MAX_THICKNESS_SEC4 > 0.0399) * 0.1; (MAX_THICKNESS_SEC5 > 0.0343) * 0.1
%

The error I get is:

Traceback (most recent call last):
File "../../../../../apps/SU2/4.2.0/bin/shape_optimization.py", line 155, in <module>
main()
File "../../../../../apps/SU2/4.2.0/bin/shape_optimization.py", line 99, in main
options.quiet )
File "../../../../../apps/SU2/4.2.0/bin/shape_optimization.py", line 138, in shape_optimization
SU2.opt.SLSQP(project,x0,xb,its,accu)
File "/apps/SU2/4.2.0/bin/SU2/opt/scipy_tools.py", line 133, in scipy_slsqp
epsilon = eps )
File "/apps/miniconda/lib/python2.7/site-packages/scipy/optimize/slsqp.py", line 206, in fmin_slsqp
constraints=cons, **opts)
File "/apps/miniconda/lib/python2.7/site-packages/scipy/optimize/slsqp.py", line 311, in _minimize_slsqp
mieq = sum(map(len, [atleast_1d(c['fun'](x, *c['args'])) for c in cons['ineq']]))
File "/apps/SU2/4.2.0/bin/SU2/opt/scipy_tools.py", line 216, in con_cieq
cons = project.con_cieq(x)
File "/apps/SU2/4.2.0/bin/SU2/opt/project.py", line 233, in con_cieq
return self._eval(konfig, func,dvs)
File "/apps/SU2/4.2.0/bin/SU2/opt/project.py", line 182, in _eval
vals = design._eval(func,*args)
File "/apps/SU2/4.2.0/bin/SU2/eval/design.py", line 142, in _eval
vals = eval_func(*inputs)
File "/apps/SU2/4.2.0/bin/SU2/eval/design.py", line 440, in con_cieq
func = su2func(this_con,config,state)
File "/apps/SU2/4.2.0/bin/SU2/eval/functions.py", line 93, in function
geometry( func_name, config, state )
File "/apps/SU2/4.2.0/bin/SU2/eval/functions.py", line 471, in geometry
info = su2run.geometry(config)
File "/apps/SU2/4.2.0/bin/SU2/run/geometry.py", line 93, in geometry
SU2_GEO(konfig)
File "/apps/SU2/4.2.0/bin/SU2/run/interface.py", line 211, in GEO
run_command( the_Command )
File "/apps/SU2/4.2.0/bin/SU2/run/interface.py", line 295, in run_command
raise exception , message
SU2.EvaluationFailure: Path = /home/GRID/s204030/SU2practice/steady_oneram6/DESIGNS/DSN_001/GEOMETRY/,
Command = /apps/SU2/4.2.0/bin/SU2_GEO config_GEO.cfg
SU2 process returned error '1'
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0


Thanks everyone for your help!

Spyros