run_MDF_opt_using_subproblem.py

import openmdao.api as om
from components.compute_pitch_angles_using_subproblem import ComputePitchAnglesUsingSubProblem
from components.design_airfoil import DesignAirfoil
from components.compute_modified_power import ComputeModifiedPower


wind_speeds = [4.0, 6.0, 8.0, 10.0]
P_rated = 500.0

prob = om.Problem()
prob.model.add_subsystem(
    "design_airfoil",
    DesignAirfoil(size=len(wind_speeds)),
    promotes=["*"],
)
prob.model.add_subsystem(
    "compute_pitch_angles_using_problem",
    ComputePitchAnglesUsingSubProblem(size=len(wind_speeds), P_rated=P_rated),
    promotes=["*"],
)
prob.model.add_subsystem(
    "compute_modified_power",
    ComputeModifiedPower(size=len(wind_speeds)),
    promotes=["*"],
)

prob.driver = om.ScipyOptimizeDriver()
prob.driver.options["optimizer"] = "SLSQP"
prob.model.approx_totals(method="fd")

prob.model.add_design_var("airfoil_design", lower=6.0, upper=15.0)
prob.model.add_design_var("drag_modifier", lower=6.0, upper=15.0)
prob.model.add_objective("modified_power")

prob.setup()

prob.set_val("wind_speeds", wind_speeds)
prob.run_driver()