def direct ( config ):
""" info = SU2.run.direct(config)
Runs an adjoint analysis with:
SU2.run.decomp()
SU2.run.CFD()
SU2.run.merge()
Assumptions:
Redundant decomposition if config.DECOMPOSED == True
Does not rename restart filename to solution filename
Adds 'direct' suffix to convergence filename
Outputs:
info - SU2 State with keys:
FUNCTIONS
HISTORY.DIRECT
FILES.DIRECT
Updates:
config.DECOMPOSED
config.MATH_PROBLEM
Executes in:
./
"""
# local copy
konfig = copy.deepcopy(config)
# decompose
su2decomp(konfig)
# setup direct problem
konfig['MATH_PROBLEM'] = 'DIRECT'
konfig['CONV_FILENAME'] = konfig['CONV_FILENAME'] + '_direct'
# Run Solution
SU2_CFD(konfig)
# merge
konfig['SOLUTION_FLOW_FILENAME'] = konfig['RESTART_FLOW_FILENAME']
su2merge(konfig)
# filenames
plot_format = konfig['OUTPUT_FORMAT']
plot_extension = su2io.get_extension(plot_format)
history_filename = konfig['CONV_FILENAME'] + plot_extension
special_cases = su2io.get_specialCases(konfig)
# get history and objectives
history = su2io.read_history( history_filename )
aerodynamics = su2io.read_aerodynamics( history_filename , special_cases )
# update super config
config.update({ 'DECOMPOSED' : konfig['DECOMPOSED'] ,
'MATH_PROBLEM' : konfig['MATH_PROBLEM'] })
# info out
info = su2io.State()
info.FUNCTIONS.update( aerodynamics )
info.FILES.DIRECT = konfig['RESTART_FLOW_FILENAME']
if 'EQUIV_AREA' in special_cases:
info.FILES.WEIGHT_NF = 'WeightNF.dat'
info.HISTORY.DIRECT = history
return info
def adjoint(config):
""" info = SU2.run.adjoint(config)
Runs an adjoint analysis with:
SU2.run.decomp()
SU2.run.CFD()
SU2.run.merge()
Assumptions:
Does not run Gradient Projection
Does not rename restart filename to solution filename
Adds 'adjoint' suffix to convergence filename
Outputs:
info - SU2 State with keys:
HISTORY.ADJOINT_NAME
FILES.ADJOINT_NAME
Updates:
config.MATH_PROBLEM
Executes in:
./
"""
# local copy
konfig = copy.deepcopy(config)
# setup problem
if konfig.get('GRADIENT_METHOD',
'CONTINUOUS_ADJOINT') == 'DISCRETE_ADJOINT':
konfig['MATH_PROBLEM'] = 'DISCRETE_ADJOINT'
else:
konfig['MATH_PROBLEM'] = 'CONTINUOUS_ADJOINT'
konfig['CONV_FILENAME'] = konfig['CONV_FILENAME'] + '_adjoint'
# Run Solution
SU2_CFD(konfig)
# merge
konfig['SOLUTION_ADJ_FILENAME'] = konfig['RESTART_ADJ_FILENAME']
su2merge(konfig)
# filenames
plot_format = konfig['OUTPUT_FORMAT']
plot_extension = su2io.get_extension(plot_format)
history_filename = konfig['CONV_FILENAME'] + plot_extension
special_cases = su2io.get_specialCases(konfig)
# get history
history = su2io.read_history(history_filename)
# update super config
config.update({
'MATH_PROBLEM': konfig['MATH_PROBLEM'],
'OBJECTIVE_FUNCTION': konfig['OBJECTIVE_FUNCTION']
})
# files out
objective = konfig['OBJECTIVE_FUNCTION']
adj_title = 'ADJOINT_' + objective
suffix = su2io.get_adjointSuffix(objective)
restart_name = konfig['RESTART_FLOW_FILENAME']
restart_name = su2io.add_suffix(restart_name, suffix)
# info out
info = su2io.State()
info.FILES[adj_title] = restart_name
info.HISTORY[adj_title] = history
return info
def adjoint( config ):
""" info = SU2.run.adjoint(config)
Runs an adjoint analysis with:
SU2.run.decomp()
SU2.run.CFD()
SU2.run.merge()
Assumptions:
Does not run Gradient Projection
Does not rename restart filename to solution filename
Adds 'adjoint' suffix to convergence filename
Outputs:
info - SU2 State with keys:
HISTORY.ADJOINT_NAME
FILES.ADJOINT_NAME
Updates:
config.MATH_PROBLEM
Executes in:
./
"""
# local copy
konfig = copy.deepcopy(config)
# setup problem
if konfig.get('GRADIENT_METHOD', 'CONTINUOUS_ADJOINT') == 'DISCRETE_ADJOINT':
konfig['MATH_PROBLEM'] = 'DISCRETE_ADJOINT'
else:
konfig['MATH_PROBLEM'] = 'CONTINUOUS_ADJOINT'
konfig['CONV_FILENAME'] = konfig['CONV_FILENAME'] + '_adjoint'
# Run Solution
SU2_CFD(konfig)
# merge
konfig['SOLUTION_ADJ_FILENAME'] = konfig['RESTART_ADJ_FILENAME']
su2merge(konfig)
# filenames
plot_format = konfig['OUTPUT_FORMAT']
plot_extension = su2io.get_extension(plot_format)
history_filename = konfig['CONV_FILENAME'] + plot_extension
special_cases = su2io.get_specialCases(konfig)
# get history
history = su2io.read_history( history_filename )
# update super config
config.update({ 'MATH_PROBLEM' : konfig['MATH_PROBLEM'] ,
'OBJECTIVE_FUNCTION' : konfig['OBJECTIVE_FUNCTION'] })
# files out
objective = konfig['OBJECTIVE_FUNCTION']
adj_title = 'ADJOINT_' + objective
suffix = su2io.get_adjointSuffix(objective)
restart_name = konfig['RESTART_FLOW_FILENAME']
restart_name = su2io.add_suffix(restart_name,suffix)
# info out
info = su2io.State()
info.FILES[adj_title] = restart_name
info.HISTORY[adj_title] = history
return info
def direct ( config ):
""" info = SU2.run.direct(config)
Runs an adjoint analysis with:
SU2.run.decomp()
SU2.run.CFD()
SU2.run.merge()
Assumptions:
Does not rename restart filename to solution filename
Adds 'direct' suffix to convergence filename
Outputs:
info - SU2 State with keys:
FUNCTIONS
HISTORY.DIRECT
FILES.DIRECT
Updates:
config.MATH_PROBLEM
Executes in:
./
"""
# local copy
konfig = copy.deepcopy(config)
# setup direct problem
konfig['MATH_PROBLEM'] = 'DIRECT'
konfig['CONV_FILENAME'] = konfig['CONV_FILENAME'] + '_direct'
direct_diff = konfig.get('DIRECT_DIFF','NO') == "YES"
# Run Solution
SU2_CFD(konfig)
# multizone cases
multizone_cases = su2io.get_multizone(konfig)
# merge
konfig['SOLUTION_FLOW_FILENAME'] = konfig['RESTART_FLOW_FILENAME']
if 'FLUID_STRUCTURE_INTERACTION' in multizone_cases:
konfig['SOLUTION_STRUCTURE_FILENAME'] = konfig['RESTART_STRUCTURE_FILENAME']
su2merge(konfig)
# filenames
plot_format = konfig['OUTPUT_FORMAT']
plot_extension = su2io.get_extension(plot_format)
history_filename = konfig['CONV_FILENAME'] + plot_extension
special_cases = su2io.get_specialCases(konfig)
# averaging final iterations
final_avg = config.get('ITER_AVERAGE_OBJ',0)
# get history and objectives
history = su2io.read_history( history_filename )
aerodynamics = su2io.read_aerodynamics( history_filename , special_cases, final_avg )
# update super config
config.update({ 'MATH_PROBLEM' : konfig['MATH_PROBLEM'] })
# info out
info = su2io.State()
info.FUNCTIONS.update( aerodynamics )
info.FILES.DIRECT = konfig['RESTART_FLOW_FILENAME']
if 'EQUIV_AREA' in special_cases:
info.FILES.WEIGHT_NF = 'WeightNF.dat'
if 'INV_DESIGN_CP' in special_cases:
info.FILES.TARGET_CP = 'TargetCp.dat'
if 'INV_DESIGN_HEATFLUX' in special_cases:
info.FILES.TARGET_HEATFLUX = 'TargetHeatFlux.dat'
info.HISTORY.DIRECT = history
return info
def direct(config):
""" info = SU2.run.direct(config)
Runs an adjoint analysis with:
SU2.run.decomp()
SU2.run.CFD()
SU2.run.merge()
Assumptions:
Redundant decomposition if config.DECOMPOSED == True
Does not rename restart filename to solution filename
Adds 'direct' suffix to convergence filename
Outputs:
info - SU2 State with keys:
FUNCTIONS
HISTORY.DIRECT
FILES.DIRECT
Updates:
config.DECOMPOSED
config.MATH_PROBLEM
Executes in:
./
"""
# local copy
konfig = copy.deepcopy(config)
# decompose
su2decomp(konfig)
# setup direct problem
konfig['MATH_PROBLEM'] = 'DIRECT'
konfig['CONV_FILENAME'] = konfig['CONV_FILENAME'] + '_direct'
# Run Solution
SU2_CFD(konfig)
# merge
konfig['SOLUTION_FLOW_FILENAME'] = konfig['RESTART_FLOW_FILENAME']
su2merge(konfig)
# filenames
plot_format = konfig['OUTPUT_FORMAT']
plot_extension = su2io.get_extension(plot_format)
history_filename = konfig['CONV_FILENAME'] + plot_extension
special_cases = su2io.get_specialCases(konfig)
# averaging final iterations
final_avg = config.get('ITER_AVERAGE_OBJ', 0)
# get history and objectives
history = su2io.read_history(history_filename)
aerodynamics = su2io.read_aerodynamics(history_filename, special_cases,
final_avg)
# update super config
config.update({
'DECOMPOSED': konfig['DECOMPOSED'],
'MATH_PROBLEM': konfig['MATH_PROBLEM']
})
# info out
info = su2io.State()
info.FUNCTIONS.update(aerodynamics)
info.FILES.DIRECT = konfig['RESTART_FLOW_FILENAME']
if 'EQUIV_AREA' in special_cases:
info.FILES.WEIGHT_NF = 'WeightNF.dat'
if 'INV_DESIGN_CP' in special_cases:
info.FILES.TARGET_CP = 'TargetCp.dat'
if 'INV_DESIGN_HEATFLUX' in special_cases:
info.FILES.TARGET_HEATFLUX = 'TargetHeatFlux.dat'
info.HISTORY.DIRECT = history
return info
def adjoint(config):
""" info = SU2.run.adjoint(config)
Runs an adjoint analysis with:
SU2.run.decomp()
SU2.run.CFD()
SU2.run.merge()
Assumptions:
Redundant decomposition if config.DECOMPOSED == True
Does not run Gradient Projection
Does not rename restart filename to solution filename
Adds 'adjoint' suffix to convergence filename
Outputs:
info - SU2 State with keys:
HISTORY.ADJOINT_NAME
FILES.ADJOINT_NAME
Updates:
config.DECOMPOSED
config.MATH_PROBLEM
Executes in:
./
"""
# local copy
konfig = copy.deepcopy(config)
# decompose
su2decomp(konfig)
# setup problem
konfig["MATH_PROBLEM"] = "ADJOINT"
konfig["CONV_FILENAME"] = konfig["CONV_FILENAME"] + "_adjoint"
# Run Solution
SU2_CFD(konfig)
# merge
konfig["SOLUTION_ADJ_FILENAME"] = konfig["RESTART_ADJ_FILENAME"]
su2merge(konfig)
# filenames
plot_format = konfig["OUTPUT_FORMAT"]
plot_extension = su2io.get_extension(plot_format)
history_filename = konfig["CONV_FILENAME"] + plot_extension
special_cases = su2io.get_specialCases(konfig)
# get history
history = su2io.read_history(history_filename)
# update super config
config.update(
{
"DECOMPOSED": konfig["DECOMPOSED"],
"MATH_PROBLEM": konfig["MATH_PROBLEM"],
"ADJ_OBJFUNC": konfig["ADJ_OBJFUNC"],
}
)
# files out
objective = konfig["ADJ_OBJFUNC"]
adj_title = "ADJOINT_" + objective
suffix = su2io.get_adjointSuffix(objective)
restart_name = konfig["RESTART_FLOW_FILENAME"]
restart_name = su2io.add_suffix(restart_name, suffix)
# info out
info = su2io.State()
info.FILES[adj_title] = restart_name
info.HISTORY[adj_title] = history
return info