[taxalian]

Hi SU2 users,
Is there already python script exists that automatically takes the converged solution of the current angle of attack as restart solution to compute the next angle of attack and so on ... to complete the lift polar. Instead of always doing this process manually to do the polars.

Thanks and regards.

[dtucker]

Originally Posted by taxalian

Hi SU2 users,
Is there already python script exists that automatically takes the converged solution of the current angle of attack as restart solution to compute the next angle of attack and so on ... to complete the lift polar. Instead of always doing this process manually to do the polars.

Thanks and regards.

I'm willing to bet someone has come up with something a little more elegant than this...but I setup a c-shell script to take care of the leg-work. It would be a prime candidate for a loop call, but this does the trick!

...I'm also not worried about space, so I just copy/paste a run directory to the next run (each AoA is a sub-directory of that) and run the script. Since I'm running on a KRAKEN Cray 5 system, I also need a submit script.

As I say...not terribly sophisticated, and I am currently "tail"ing each output file one by one and transferring the CL and CD values. If you get anything more slick I'd be interested in taking a look!

Good Luck!
Dave

[rktchip]

This is a good example to try with the new SU2 python package! Here's an example python script -

#!/usr/bin/env python 

# imports
import SU2
import numpy as np
import pylab as plt
from copy import deepcopy

# load config, start state
config = SU2.io.Config('inv_NACA0012.cfg')
state  = SU2.io.State()

# prepare config
config.NUMBER_PART = 2
config.EXT_ITER    = 99999
config.RESTART_SOL = 'YES'

# find solution files if they exist
state.find_files(config)

# angles to run
angles = np.linspace(-10.,10.,7)

# start results data
results = SU2.util.bunch()
results.AoA  = angles
results.DRAG = []
results.LIFT = []

# iterate angles
for angle in angles:
    
    # local config and state
    konfig = deepcopy(config)
    ztate  = deepcopy(state)
    
    # set angle of attack
    konfig.AoA = angle
    print 'AoA = ' , konfig.AoA
    
    # run su2
    drag = SU2.eval.func('DRAG',konfig,ztate)
    lift = SU2.eval.func('LIFT',konfig,ztate)
    
    # append results
    results.DRAG.append(drag)
    results.LIFT.append(lift)
    
#: for each angle

# plotting
plt.figure()
plt.plot( results.AoA , results.DRAG )
plt.show()    

# save data
SU2.io.save_data('results.pkl',results)