fname_analysis_options = os.path.join(run_dir,
f'analysis_options.{istep}.yaml')
if MPI:
rank = MPI.COMM_WORLD.Get_rank()
else:
rank = 0
if rank == 0:
print('STEP', istep)
## Update analysis options
aopt = load_yaml(os.path.join(run_dir, 'analysis_options.start.yaml'))
aopt['general']['folder_output'] = f'outputs.{istep}'
aopt['general']['fname_output'] = f'NREL-2p5-116-step{istep}'
save_yaml(fname_analysis_options, aopt)
tt = time.time()
# step 1: manually updated turbine rating, rotor size -- NO OPT
# - rated_power: 2500000.0 W
# - rotor_diameter: 116 m
# - hub_height: 80 m
wt_opt, modeling_options, opt_options = run_wisdem(fname_wt_input,
fname_modeling_options,
fname_analysis_options)
if rank == 0:
print('Run time: %f' % (time.time() - tt))
sys.stdout.flush()
print('STEP', istep)
## Update analysis options
aopt = load_yaml(os.path.join(run_dir, 'analysis_options.start.yaml'))
aopt['general']['folder_output'] = f'outputs.{istep}'
aopt['general']['fname_output'] = f'NREL-1p7-103-step{istep}'
# - blade-mass opt constrained by tip deflection
aopt['driver']['optimization']['flag'] = True
aopt['design_variables']['blade']['structure']['spar_cap_ss'][
'flag'] = True
aopt['design_variables']['blade']['structure']['spar_cap_ps'][
'flag'] = True
aopt['constraints']['blade']['tip_deflection']['flag'] = True
aopt['merit_figure'] = 'blade_mass'
save_yaml(fname_analysis_options, aopt)
## Update modeling options
mopt = load_yaml(
os.path.join(run_dir, f'outputs.{istep-1}',
f'NREL-1p7-103-step{istep-1}-modeling.yaml'))
save_yaml(fname_modeling_options, mopt)
if MPI:
MPI.COMM_WORLD.Barrier()
model_changes = {}
tt = time.time()
wt_opt, modeling_options, opt_options = run_wisdem(
