[balav94]

Hi everyone!

I am currently working on Shape Optimization of automobiles using Adjoint solvers. I have been working with an OpenFOAM based solver called "adjointSensitivityFoam". However, I felt that SU2 has better libraries for this purpose and have started exploring the same.

While I was trying to run the "steady_oneram6" testcase by following the tutorial here (https://github.com/su2code/SU2/wiki/...f-a-Fixed-Wing), specifically when i type "shape_optimization.py -f inv_ONERAM6_adv.cfg", I get the following error.

-------------------------------------------------------------------------
| ___ _ _ ___ |
| / __| | | |_ ) Release 5.0.0 "Raven" |
| \__ \ |_| |/ / |
| |___/\___//___| Aerodynamic Shape Optimization Script |
| |
-------------------------------------------------------------------------
| SU2 Original Developers: Dr. Francisco D. Palacios. |
| Dr. Thomas D. Economon. |
-------------------------------------------------------------------------
| SU2 Developers: |
| - Prof. Juan J. Alonso's group at Stanford University. |
| - Prof. Piero Colonna's group at Delft University of Technology. |
| - Prof. Nicolas R. Gauger's group at Kaiserslautern U. of Technology. |
| - Prof. Alberto Guardone's group at Polytechnic University of Milan. |
| - Prof. Rafael Palacios' group at Imperial College London. |
-------------------------------------------------------------------------
| Copyright (C) 2012-2017 SU2, the open-source CFD code. |
| |
| SU2 is free software; you can redistribute it and/or |
| modify it under the terms of the GNU Lesser General Public |
| License as published by the Free Software Foundation; either |
| version 2.1 of the License, or (at your option) any later version. |
| |
| SU2 is distributed in the hope that it will be useful, |
| but WITHOUT ANY WARRANTY; without even the implied warranty of |
| MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
| Lesser General Public License for more details. |
| |
| You should have received a copy of the GNU Lesser General Public |
| License along with SU2. If not, see <http://www.gnu.org/licenses/>. |
-------------------------------------------------------------------------
Found: mesh_ONERAM6_inv_FFD.su2
Found: solution_flow.dat
Found: solution_adj_cl.dat
Found: solution_adj_cd.dat
New Project: ./
Removing old designs in 10s. Done!

Sequential Least SQuares Programming (SLSQP) parameters:
Number of design variables: 176 ( 176 )
Objective function scaling factor: [0.1]
Maximum number of iterations: 100
Requested accuracy: 1e-11
Initial guess for the independent variable(s): [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, 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, 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, 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, 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, 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, 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, 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, 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, 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, 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, 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, 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, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]
Lower and upper bound for each independent variable: [(-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1)]

Traceback (most recent call last):
File "/usr/local/bin/shape_optimization.py", line 169, in <module>
main()
File "/usr/local/bin/shape_optimization.py", line 104, in main
options.quiet )
File "/usr/local/bin/shape_optimization.py", line 145, in shape_optimization
SU2.opt.SLSQP(project,x0,xb,its,accu)
File "/usr/local/bin/SU2/opt/scipy_tools.py", line 140, in scipy_slsqp
epsilon = eps )
File "/usr/local/lib/python2.7/dist-packages/scipy/optimize/slsqp.py", line 196, in fmin_slsqp
constraints=cons, **opts)
File "/usr/local/lib/python2.7/dist-packages/scipy/optimize/slsqp.py", line 300, in _minimize_slsqp
mieq = sum(map(len, [atleast_1d(c['fun'](x, *c['args'])) for c in cons['ineq']]))
File "/usr/local/bin/SU2/opt/scipy_tools.py", line 464, in con_cieq
cons = project.con_cieq(x)
File "/usr/local/bin/SU2/opt/project.py", line 235, in con_cieq
return self._eval(konfig, func,dvs)
File "/usr/local/bin/SU2/opt/project.py", line 184, in _eval
vals = design._eval(func,*args)
File "/usr/local/bin/SU2/eval/design.py", line 144, in _eval
vals = eval_func(*inputs)
File "/usr/local/bin/SU2/eval/design.py", line 457, in con_cieq
func = su2func(this_con,config,state)
File "/usr/local/bin/SU2/eval/functions.py", line 101, in function
geometry( func_name, config, state )
File "/usr/local/bin/SU2/eval/functions.py", line 488, in geometry
info = su2run.geometry(config)
File "/usr/local/bin/SU2/run/geometry.py", line 102, in geometry
functions = su2io.tools.read_plot(func_filename)
File "/usr/local/bin/SU2/io/tools.py", line 120, in read_plot
line_data = [ float(x.strip()) for x in line_data ]
ValueError: could not convert string to float:

My guess is there are some missing python dependencies. Since I am new, I am finding it difficult to understand the root cause of this error. It would be great if anyone could help me out.

Thanks in advance!

[hlk]

When running the shape optimization simulations, individual log files are stored in the DESIGNS/... subdirectories. Do you seen any more information in these log files?


Based on the date of your post, you may also benefit from downloading the most recent version.

Originally Posted by balav94

Hi everyone!

I am currently working on Shape Optimization of automobiles using Adjoint solvers. I have been working with an OpenFOAM based solver called "adjointSensitivityFoam". However, I felt that SU2 has better libraries for this purpose and have started exploring the same.

While I was trying to run the "steady_oneram6" testcase by following the tutorial here (https://github.com/su2code/SU2/wiki/...f-a-Fixed-Wing), specifically when i type "shape_optimization.py -f inv_ONERAM6_adv.cfg", I get the following error.

-------------------------------------------------------------------------
| ___ _ _ ___ |
| / __| | | |_ ) Release 5.0.0 "Raven" |
| \__ \ |_| |/ / |
| |___/\___//___| Aerodynamic Shape Optimization Script |
| |
-------------------------------------------------------------------------
| SU2 Original Developers: Dr. Francisco D. Palacios. |
| Dr. Thomas D. Economon. |
-------------------------------------------------------------------------
| SU2 Developers: |
| - Prof. Juan J. Alonso's group at Stanford University. |
| - Prof. Piero Colonna's group at Delft University of Technology. |
| - Prof. Nicolas R. Gauger's group at Kaiserslautern U. of Technology. |
| - Prof. Alberto Guardone's group at Polytechnic University of Milan. |
| - Prof. Rafael Palacios' group at Imperial College London. |
-------------------------------------------------------------------------
| Copyright (C) 2012-2017 SU2, the open-source CFD code. |
| |
| SU2 is free software; you can redistribute it and/or |
| modify it under the terms of the GNU Lesser General Public |
| License as published by the Free Software Foundation; either |
| version 2.1 of the License, or (at your option) any later version. |
| |
| SU2 is distributed in the hope that it will be useful, |
| but WITHOUT ANY WARRANTY; without even the implied warranty of |
| MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
| Lesser General Public License for more details. |
| |
| You should have received a copy of the GNU Lesser General Public |
| License along with SU2. If not, see <http://www.gnu.org/licenses/>. |
-------------------------------------------------------------------------
Found: mesh_ONERAM6_inv_FFD.su2
Found: solution_flow.dat
Found: solution_adj_cl.dat
Found: solution_adj_cd.dat
New Project: ./
Removing old designs in 10s. Done!

Sequential Least SQuares Programming (SLSQP) parameters:
Number of design variables: 176 ( 176 )
Objective function scaling factor: [0.1]
Maximum number of iterations: 100
Requested accuracy: 1e-11
Initial guess for the independent variable(s): [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, 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, 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, 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, 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, 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, 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, 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, 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, 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, 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, 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, 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, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]
Lower and upper bound for each independent variable: [(-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1), (-0.1, 0.1)]

Traceback (most recent call last):
File "/usr/local/bin/shape_optimization.py", line 169, in <module>
main()
File "/usr/local/bin/shape_optimization.py", line 104, in main
options.quiet )
File "/usr/local/bin/shape_optimization.py", line 145, in shape_optimization
SU2.opt.SLSQP(project,x0,xb,its,accu)
File "/usr/local/bin/SU2/opt/scipy_tools.py", line 140, in scipy_slsqp
epsilon = eps )
File "/usr/local/lib/python2.7/dist-packages/scipy/optimize/slsqp.py", line 196, in fmin_slsqp
constraints=cons, **opts)
File "/usr/local/lib/python2.7/dist-packages/scipy/optimize/slsqp.py", line 300, in _minimize_slsqp
mieq = sum(map(len, [atleast_1d(c['fun'](x, *c['args'])) for c in cons['ineq']]))
File "/usr/local/bin/SU2/opt/scipy_tools.py", line 464, in con_cieq
cons = project.con_cieq(x)
File "/usr/local/bin/SU2/opt/project.py", line 235, in con_cieq
return self._eval(konfig, func,dvs)
File "/usr/local/bin/SU2/opt/project.py", line 184, in _eval
vals = design._eval(func,*args)
File "/usr/local/bin/SU2/eval/design.py", line 144, in _eval
vals = eval_func(*inputs)
File "/usr/local/bin/SU2/eval/design.py", line 457, in con_cieq
func = su2func(this_con,config,state)
File "/usr/local/bin/SU2/eval/functions.py", line 101, in function
geometry( func_name, config, state )
File "/usr/local/bin/SU2/eval/functions.py", line 488, in geometry
info = su2run.geometry(config)
File "/usr/local/bin/SU2/run/geometry.py", line 102, in geometry
functions = su2io.tools.read_plot(func_filename)
File "/usr/local/bin/SU2/io/tools.py", line 120, in read_plot
line_data = [ float(x.strip()) for x in line_data ]
ValueError: could not convert string to float:

My guess is there are some missing python dependencies. Since I am new, I am finding it difficult to understand the root cause of this error. It would be great if anyone could help me out.

Thanks in advance!