[lbfhappy]

Hello All,

As I try Optimal Shape Design of a Rotating Airfoil (tutorial 8), it showed it has non-physical points and it diverged. I I entered "python shape_optimization.py -f rotating_NACA0012.cfg -q False" for optimization. I attached configure and mesh file here. And the error is as follow:



802260.070862141.046635147.574381-10000.000000-10000.000000
802270.070862141.045031147.572994-10000.00000010000.000000
802280.070862141.044549147.572567-10000.00000010000.000000
802290.070862141.045506147.573470-10000.00000010000.000000
802300.070862141.046289147.574196-10000.00000010000.000000
802310.070862141.157858147.692308-10000.00000010000.000000
802320.070862141.239897147.771173-10000.00000010000.000000
802330.070862142.549743149.108275-10000.00000010000.000000
802340.070862142.498160149.085855-10000.00000010000.000000
802350.070862142.783924149.366008-10000.00000010000.000000
802360.070862142.383524148.966505-10000.00000010000.000000
802370.070862142.121731148.697318-10000.00000010000.000000
802380.070862141.989612148.562711-10000.00000010000.000000
802390.070862141.917804148.487499-10000.00000010000.000000

---------------------- Local Time Stepping Summary ----------------------
MG level: 0 -> Min. DT: 1.67123e-07. Max. DT: 0.0057421. CFL: 6.
-------------------------------------------------------------------------

----------------------- Residual Evolution Summary ----------------------
log10[Maximum residual]: 143.591.
Maximum residual point 5231, located at (11.3915, 14.0631).
There are 326 non-physical points in the solution.
-------------------------------------------------------------------------

IterTime(s)Res[Rho]Res[RhoE]CLift(Total)CDrag(Total)
802400.070862141.919205148.488830-10000.00000010000.000000
802410.070862141.890851148.461624-10000.00000010000.000000
802420.070862141.907374148.476271-10000.00000010000.000000
802430.070862141.850338148.419919-10000.00000010000.000000
802440.070862141.920849148.493587-10000.00000010000.000000
802450.070862141.347845147.933063-10000.000000-10000.000000
802460.070861141.243377147.867069-10000.000000-10000.000000
802470.070861141.601230147.898575-10000.000000-10000.000000
802480.070861141.281252147.934883-10000.000000-10000.000000
802490.070862141.419211147.958188-10000.000000-10000.000000
802500.070862141.339054147.971665-10000.000000-10000.000000

-------------------------- File Output Summary --------------------------
Writing comma-separated values (CSV) surface files.
Merging connectivities in the Master node.
Merging coordinates in the Master node.
Merging solution in the Master node.
Writing SU2 native restart file.
Writing Tecplot ASCII file volume solution file.
Writing Tecplot ASCII surface solution file.
Writing the forces breakdown file.
-------------------------------------------------------------------------

802510.070864141.364314147.995806-10000.000000-10000.000000
802520.070864141.336538147.943929-10000.000000-10000.000000
802530.070864141.725985148.332103-10000.000000-10000.000000
802540.070864141.350102147.95459110000.000000-10000.000000
802550.070864141.511964148.11914110000.000000-10000.000000
802560.070864141.683121148.30592310000.000000-10000.000000
802570.070864141.874060148.48854910000.000000-10000.000000
802580.070864141.871048148.48667410000.000000-10000.000000
802590.070864141.776553148.38919310000.000000-10000.000000
802600.070864141.660417148.26944410000.000000-10000.000000
802610.070864141.519441148.12146410000.000000-10000.000000
802620.070864141.347365147.94512310000.000000-10000.000000
802630.070864141.133635147.72258810000.000000-10000.000000
802640.070864140.895956147.47376510000.000000-10000.000000
802650.070864140.745786147.19990810000.000000-10000.000000
802660.070864140.505270146.90049210000.000000-10000.000000
802670.070864140.284874146.60089610000.000000-10000.000000
802680.070864140.205566146.38919910000.000000-10000.000000
802690.070864140.411331146.42652710000.000000-10000.000000
802700.070864140.398680146.38738410000.000000-10000.000000
802710.070864140.287615146.36191810000.00000010000.000000
802720.070864142.242658148.30991810000.000000-10000.000000
802730.070864142.135021 148.23228610000.00000010000.000000
802740.070864142.119853148.20461310000.00000010000.000000
802750.070864142.136155148.21911810000.00000010000.000000
802760.070864142.136514148.21933210000.00000010000.000000
802770.070864142.148057148.23146410000.00000010000.000000
802780.070864144.418341150.15343510000.00000010000.000000
802790.070864144.095130150.12161810000.00000010000.000000

---------------------- Local Time Stepping Summary ----------------------
MG level: 0 -> Min. DT: 1.96599e-07. Max. DT: 0.00548642. CFL: 6.
-------------------------------------------------------------------------

----------------------- Residual Evolution Summary ----------------------
log10[Maximum residual]: 146.338.
Maximum residual point 206, located at (12.7485, 15.4103).
There are 195 non-physical points in the solution.
-------------------------------------------------------------------------

IterTime(s)Res[Rho]Res[RhoE]CLift(Total)CDrag(Total)
802800.070864144.480049150.24787910000.000000-10000.000000
802810.070864144.470469150.26627810000.000000-10000.000000
802820.070864144.348874150.17629610000.000000-10000.000000
802830.070864144.461077150.30676610000.000000-10000.000000
802840.070864144.512907150.362175-10000.000000-10000.000000
802850.070864144.560500150.509304-10000.000000-10000.000000
802860.070864144.540206150.519431-10000.000000-10000.000000
802870.070864144.640407150.702012-10000.000000-10000.000000
802880.070864144.792297150.905065-10000.000000-10000.000000
802890.070864144.884054151.01332610000.000000-10000.000000
802900.070864144.719157150.82559610000.000000-10000.000000
802910.070864144.676127150.78552710000.000000-10000.000000
802920.070864144.713505150.81388510000.00000010000.000000
802930.070864144.727893150.82528210000.00000010000.000000
802940.070864144.751866150.84807710000.000000-10000.000000
802950.070864144.787424150.87300810000.000000-10000.000000
802960.070864144.810048150.88856310000.00000010000.000000
802970.070864144.880866150.92778710000.00000010000.000000
802980.070864144.945385150.97423510000.00000010000.000000
802990.070864145.016040151.02283810000.00000010000.000000
803000.070864145.206372151.18454510000.00000010000.000000
803010.070864145.527455151.50660610000.00000010000.000000

!!! Error: SU2 has diverged. Now exiting... !!!

Traceback (most recent call last):
File "shape_optimization.py", line 154, in <module>
main()
File "shape_optimization.py", line 99, in main
options.quiet)
File "shape_optimization.py", line 137, in shape_optimization
SU2.opt.SLSQP(project,x0,xb,its,accu)
File "/home/gentai/Workzone/LBF/SU2-master/bin/SU2/opt/scipy_tools.py", line 128, in scipy_slsqp
epsilon= eps)
File "/home/gentai/Enthought/Canopy_64bit/User/lib/python2.7/site-packages/scipy/optimize/slsqp.py", line 206, in fmin_slsqp
constraints=cons, **opts)
File "/home/gentai/Enthought/Canopy_64bit/User/lib/python2.7/site-packages/scipy/optimize/slsqp.py", line 360, in _minimize_slsqp
fx = func(x)
File "/home/gentai/Enthought/Canopy_64bit/User/lib/python2.7/site-packages/scipy/optimize/optimize.py", line 281, in function_wrapper
return function(*(wrapper_args + args))
File "/home/gentai/Workzone/LBF/SU2-master/bin/SU2/opt/scipy_tools.py", line 144, in obj_f
obj = project.obj_f(x)
File "/home/gentai/Workzone/LBF/SU2-master/bin/SU2/opt/project.py", line 213, in obj_f
return self._eval(konfig, func,dvs)
File "/home/gentai/Workzone/LBF/SU2-master/bin/SU2/opt/project.py", line 182, in _eval
vals = design._eval(func,*args)
File "/home/gentai/Workzone/LBF/SU2-master/bin/SU2/eval/design.py", line 142, in _eval
vals = eval_func(*inputs)
File "/home/gentai/Workzone/LBF/SU2-master/bin/SU2/eval/design.py", line 240, in obj_f
func = su2func(this_obj,config,state)
File "/home/gentai/Workzone/LBF/SU2-master/bin/SU2/eval/functions.py", line 85, in function
aerodynamics( config, state )
File "/home/gentai/Workzone/LBF/SU2-master/bin/SU2/eval/functions.py", line 224, in aerodynamics
info = su2run.direct(config)
File "/home/gentai/Workzone/LBF/SU2-master/bin/SU2/run/direct.py", line 81, in direct
SU2_CFD(konfig)
File "/home/gentai/Workzone/LBF/SU2-master/bin/SU2/run/interface.py", line 110, in CFD
run_command( the_Command )
File "/home/gentai/Workzone/LBF/SU2-master/bin/SU2/run/interface.py", line 268, in run_command
raise exception , message
SU2.DivergenceFailure: Path = /home/gentai/Workzone/LBF/SU2-master/bin/DESIGNS/DSN_049/DIRECT/,
Command = /usr/local/bin/SU2_CFD config_CFD.cfg
SU2 process returned error '2'

[gentai@node09 bin]$

I even met the same error when I try optimization of tutorial 9. But I don't know why? Does anyone met the same problem?

Thank you!