[elninoft9]

While conducting further analysis into the "on the "Unconstrained shape design of a transonic inviscid airfoil at a cte. AoA" tutorial", for the NACA 0012 airfoil, I am encountering traceback error for the following simulation. I have changed the mach number to 0.91 and the OPT_ITERATIONS to 150. I get this error: Traceback (most recent call last):
File "shape_optimization.py", line 183, in <module>
main()
File "shape_optimization.py", line 92, in main
shape_optimization( options.filename ,
File "shape_optimization.py", line 159, in shape_optimization
SU2.opt.SLSQP(project,x0,xb,its,accu)
File "/usr/local/bin/SU2/opt/scipy_tools.py", line 120, in scipy_slsqp
outputs = fmin_slsqp( x0 = x0 ,
File "/home/suraj/.local/lib/python3.8/site-packages/scipy/optimize/_slsqp_py.py", line 206, in fmin_slsqp
res = _minimize_slsqp(func, x0, args, jac=fprime, bounds=bounds,
File "/home/suraj/.local/lib/python3.8/site-packages/scipy/optimize/_slsqp_py.py", line 428, in _minimize_slsqp
fx = wrapped_fun(x)
File "/home/suraj/.local/lib/python3.8/site-packages/scipy/optimize/_optimize.py", line 277, in eval
return func(x)
File "/home/suraj/.local/lib/python3.8/site-packages/scipy/optimize/_differentiable_functions.py", line 267, in fun
self._update_fun()
File "/home/suraj/.local/lib/python3.8/site-packages/scipy/optimize/_differentiable_functions.py", line 251, in _update_fun
self._update_fun_impl()
File "/home/suraj/.local/lib/python3.8/site-packages/scipy/optimize/_differentiable_functions.py", line 155, in update_fun
self.f = fun_wrapped(self.x)
File "/home/suraj/.local/lib/python3.8/site-packages/scipy/optimize/_differentiable_functions.py", line 137, in fun_wrapped
fx = fun(np.copy(x), *args)
File "/usr/local/bin/SU2/opt/scipy_tools.py", line 383, in obj_f
obj_list = project.obj_f(x)
File "/usr/local/bin/SU2/opt/project.py", line 237, in obj_f
return self._eval(konfig, func,dvs)
File "/usr/local/bin/SU2/opt/project.py", line 206, in _eval
vals = design._eval(func,*args)
File "/usr/local/bin/SU2/eval/design.py", line 147, in _eval
vals = eval_func(*inputs)
File "/usr/local/bin/SU2/eval/design.py", line 244, in obj_f
func += su2func(this_obj,config,state) * sign * scale * global_factor
File "/usr/local/bin/SU2/eval/functions.py", line 92, in function
aerodynamics( config, state )
File "/usr/local/bin/SU2/eval/functions.py", line 274, in aerodynamics
info = su2run.direct(config)
File "/usr/local/bin/SU2/run/direct.py", line 77, in direct
SU2_CFD(konfig)
File "/usr/local/bin/SU2/run/interface.py", line 112, in CFD
run_command( the_Command )
File "/usr/local/bin/SU2/run/interface.py", line 270, in run_command
raise exception(message)
SU2.EvaluationFailure: Path = /home/suraj/Documents/GitHub/FYP/Tutorials/design/15_my_Inviscid_2D_Unconstrained_NACA0012 (copy)/DESIGNS/DSN_015/DIRECT/,
Command = /usr/local/bin/SU2_CFD config_CFD.cfg
SU2 process returned error '1'

Any suggestions on how to get around this? I even tried changing the CFL_REDUCTION_ADJFLOW from 0.8 to 0.9, but to no avail.

[bigfootedrockmidget]

SU2.EvaluationFailure: Path = /home/suraj/Documents/GitHub/FYP/Tutorials/design/15_my_Inviscid_2D_Unconstrained_NACA0012 (copy)/DESIGNS/DSN_015/DIRECT/

It looks like it already did 14 design cycles and fails at the 15th design cycle. So my guess is that the 15th design failed to converge and therefore did not produce any output. So do you now have folders named DSN_001,DSN_002,... containing solution files? You can then have a look at the solution files in the subfolder DSN_00x/DIRECT, which is the CFD solution of the new design at design cycle x.



I don't know how the original test case stops, does it stop after a fixed number of design cycles?

[elninoft9]

Yes there are 15 design files as you mentioned. When I inspected the DSN_015/DIRECT folder, I can see the history_direct.csv file, where the rms(Rho) values are indicated. It begins from about -1.81 and ends of at 11.68, in the 205th iteration. Is there anything that we can observe from this file or these data? I'm not sure why the iterations stop after the 15th design cycle. One thing i noted was that for DSN_015/DIRECT, there were no surface_flow.csv file or surface_flow.vtu files, which is a clear indication of the failure to converge for this design cycle.