nrel_csm_wrapper.ipynb

{
  "nbformat": 4,
  "nbformat_minor": 0,
  "metadata": {
    "colab": {
      "provenance": [],
      "authorship_tag": "ABX9TyPKRurjLciNlrUd2DpHjZy3",
      "include_colab_link": true
    },
    "kernelspec": {
      "name": "python3",
      "display_name": "Python 3"
    },
    "language_info": {
      "name": "python"
    }
  },
  "cells": [
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "view-in-github",
        "colab_type": "text"
      },
      "source": [
        "<a href=\"https://colab.research.google.com/github/mifm/examples/blob/master/nrel_csm_wrapper.ipynb\" target=\"_parent\"><img src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open In Colab\"/></a>"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": 1,
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "Slxka6QuWMFj",
        "outputId": "54e360dd-0037-4c0f-c586-47b7eefaeb06"
      },
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "Looking in indexes: https://pypi.org/simple, https://us-python.pkg.dev/colab-wheels/public/simple/\n",
            "Collecting git+https://github.com/WISDEM/WISDEM.git\n",
            "  Cloning https://github.com/WISDEM/WISDEM.git to /tmp/pip-req-build-85yn0fu3\n",
            "  Running command git clone --filter=blob:none --quiet https://github.com/WISDEM/WISDEM.git /tmp/pip-req-build-85yn0fu3\n",
            "  Resolved https://github.com/WISDEM/WISDEM.git to commit fb7cddf67d21a12f1a8e04e98833f49a447af7bf\n",
            "  Installing build dependencies ... \u001b[?25l\u001b[?25hdone\n",
            "  Getting requirements to build wheel ... \u001b[?25l\u001b[?25hdone\n",
            "  Installing backend dependencies ... \u001b[?25l\u001b[?25hdone\n",
            "  Preparing metadata (pyproject.toml) ... \u001b[?25l\u001b[?25hdone\n",
            "Building wheels for collected packages: WISDEM\n",
            "  Building wheel for WISDEM (pyproject.toml) ... \u001b[?25l\u001b[?25hdone\n",
            "  Created wheel for WISDEM: filename=WISDEM-3.8-py3-none-any.whl size=4581333 sha256=1cfd4cd9df2e424b6233f1d0ffc7ab5ac1828cc77e82558ab845ca79b2386cd6\n",
            "  Stored in directory: /tmp/pip-ephem-wheel-cache-1ghzodfb/wheels/52/18/6c/8683bc75ae1ca6876395fa1d46d5de2eb35310bd132f8928f3\n",
            "Successfully built WISDEM\n",
            "Installing collected packages: WISDEM\n",
            "Successfully installed WISDEM-3.8\n"
          ]
        }
      ],
      "source": [
        "!pip install git+https://github.com/WISDEM/WISDEM.git --no-deps"
      ]
    },
    {
      "cell_type": "code",
      "source": [
        "!pip install openmdao"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "1h9DOD4OWdoD",
        "outputId": "56ec03d0-8016-4d57-f533-ebf03a274305"
      },
      "execution_count": 2,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "Looking in indexes: https://pypi.org/simple, https://us-python.pkg.dev/colab-wheels/public/simple/\n",
            "Collecting openmdao\n",
            "  Downloading openmdao-3.23.0.tar.gz (5.3 MB)\n",
            "\u001b[2K     \u001b[90m━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━\u001b[0m \u001b[32m5.3/5.3 MB\u001b[0m \u001b[31m10.9 MB/s\u001b[0m eta \u001b[36m0:00:00\u001b[0m\n",
            "\u001b[?25h  Preparing metadata (setup.py) ... \u001b[?25l\u001b[?25hdone\n",
            "Requirement already satisfied: networkx>=2.0 in /usr/local/lib/python3.8/dist-packages (from openmdao) (3.0)\n",
            "Requirement already satisfied: numpy in /usr/local/lib/python3.8/dist-packages (from openmdao) (1.21.6)\n",
            "Requirement already satisfied: scipy in /usr/local/lib/python3.8/dist-packages (from openmdao) (1.7.3)\n",
            "Requirement already satisfied: requests in /usr/local/lib/python3.8/dist-packages (from openmdao) (2.25.1)\n",
            "Requirement already satisfied: packaging in /usr/local/lib/python3.8/dist-packages (from openmdao) (21.3)\n",
            "Requirement already satisfied: pyparsing!=3.0.5,>=2.0.2 in /usr/local/lib/python3.8/dist-packages (from packaging->openmdao) (3.0.9)\n",
            "Requirement already satisfied: certifi>=2017.4.17 in /usr/local/lib/python3.8/dist-packages (from requests->openmdao) (2022.12.7)\n",
            "Requirement already satisfied: chardet<5,>=3.0.2 in /usr/local/lib/python3.8/dist-packages (from requests->openmdao) (4.0.0)\n",
            "Requirement already satisfied: urllib3<1.27,>=1.21.1 in /usr/local/lib/python3.8/dist-packages (from requests->openmdao) (1.24.3)\n",
            "Requirement already satisfied: idna<3,>=2.5 in /usr/local/lib/python3.8/dist-packages (from requests->openmdao) (2.10)\n",
            "Building wheels for collected packages: openmdao\n",
            "  Building wheel for openmdao (setup.py) ... \u001b[?25l\u001b[?25hdone\n",
            "  Created wheel for openmdao: filename=openmdao-3.23.0-py3-none-any.whl size=5806407 sha256=53646e0360ad29be50fdbe31013da6088dfbaebb3ff60deb4cbefa86e7e70015\n",
            "  Stored in directory: /root/.cache/pip/wheels/35/7d/84/03d1663a4f2694741021ba0a18aee54e254bbc89cf35df0ebd\n",
            "Successfully built openmdao\n",
            "Installing collected packages: openmdao\n",
            "\u001b[31mERROR: pip's dependency resolver does not currently take into account all the packages that are installed. This behaviour is the source of the following dependency conflicts.\n",
            "wisdem 3.8 requires dearpygui, which is not installed.\n",
            "wisdem 3.8 requires marmot-agents>=0.2.5, which is not installed.\n",
            "wisdem 3.8 requires nlopt, which is not installed.\n",
            "wisdem 3.8 requires pydoe2, which is not installed.\n",
            "wisdem 3.8 requires python-benedict, which is not installed.\n",
            "wisdem 3.8 requires ruamel.yaml, which is not installed.\n",
            "wisdem 3.8 requires simpy, which is not installed.\u001b[0m\u001b[31m\n",
            "\u001b[0mSuccessfully installed openmdao-3.23.0\n"
          ]
        }
      ]
    },
    {
      "cell_type": "code",
      "source": [
        "!pip install ruamel_yaml"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "pmikAk0JW2Zj",
        "outputId": "ef4a9c21-e8b5-4ed8-ddbb-ff75a58ecfb9"
      },
      "execution_count": 4,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "Looking in indexes: https://pypi.org/simple, https://us-python.pkg.dev/colab-wheels/public/simple/\n",
            "Collecting ruamel_yaml\n",
            "  Downloading ruamel.yaml-0.17.21-py3-none-any.whl (109 kB)\n",
            "\u001b[2K     \u001b[90m━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━\u001b[0m \u001b[32m109.5/109.5 KB\u001b[0m \u001b[31m4.0 MB/s\u001b[0m eta \u001b[36m0:00:00\u001b[0m\n",
            "\u001b[?25hCollecting ruamel.yaml.clib>=0.2.6\n",
            "  Downloading ruamel.yaml.clib-0.2.7-cp38-cp38-manylinux_2_17_x86_64.manylinux2014_x86_64.manylinux_2_24_x86_64.whl (555 kB)\n",
            "\u001b[2K     \u001b[90m━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━\u001b[0m \u001b[32m555.3/555.3 KB\u001b[0m \u001b[31m14.8 MB/s\u001b[0m eta \u001b[36m0:00:00\u001b[0m\n",
            "\u001b[?25hInstalling collected packages: ruamel.yaml.clib, ruamel_yaml\n",
            "\u001b[31mERROR: pip's dependency resolver does not currently take into account all the packages that are installed. This behaviour is the source of the following dependency conflicts.\n",
            "wisdem 3.8 requires dearpygui, which is not installed.\n",
            "wisdem 3.8 requires marmot-agents>=0.2.5, which is not installed.\n",
            "wisdem 3.8 requires nlopt, which is not installed.\n",
            "wisdem 3.8 requires pydoe2, which is not installed.\n",
            "wisdem 3.8 requires python-benedict, which is not installed.\n",
            "wisdem 3.8 requires simpy, which is not installed.\u001b[0m\u001b[31m\n",
            "\u001b[0mSuccessfully installed ruamel.yaml.clib-0.2.7 ruamel_yaml-0.17.21\n"
          ]
        }
      ]
    },
    {
      "cell_type": "code",
      "source": [
        "\n",
        "from wisdem.nrelcsm.nrel_csm_mass_2015 import nrel_csm_2015\n",
        "# 0 ---------- (marker for docs)\n",
        "\n",
        "# 1 ---------- (marker for docs)\n",
        "# OpenMDAO Problem instance\n",
        "prob = om.Problem(reports=False)\n",
        "prob.model = nrel_csm_2015()\n",
        "prob.setup()\n",
        "# 1 ---------- (marker for docs)\n",
        "\n",
        "# 2 ---------- (marker for docs)\n",
        "# Initialize variables for NREL CSM\n",
        "prob[\"machine_rating\"] = 5000.0\n",
        "prob[\"rotor_diameter\"] = 126.0\n",
        "prob[\"turbine_class\"] = 2\n",
        "prob[\"hub_height\"] = 90.0\n",
        "prob[\"blade_number\"] = 3\n",
        "prob[\"blade_has_carbon\"] = False\n",
        "prob[\"max_tip_speed\"] = 80.0\n",
        "prob[\"max_efficiency\"] = 0.90\n",
        "prob[\"main_bearing_number\"] = 2\n",
        "prob[\"crane\"] = True\n",
        "# 2 ---------- (marker for docs)\n",
        "\n",
        "# 3 ---------- (marker for docs)\n",
        "# Evaluate the model\n",
        "prob.run_model()\n",
        "# 3 ---------- (marker for docs)\n",
        "\n",
        "# 4 ---------- (marker for docs)\n",
        "# Print all intermediate inputs and outputs to the screen\n",
        "prob.model.list_inputs(units=True)\n",
        "prob.model.list_outputs(units=True)\n",
        "# 4 ---------- (marker for docs)"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "nkOm_E8bWmLy",
        "outputId": "0fb54ea9-6ff6-4bc0-afa9-4d13c751374d"
      },
      "execution_count": 9,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "231 Input(s) in 'model'\n",
            "\n",
            "varname                                    val                 units      \n",
            "-----------------------------------------  ------------------  -----------\n",
            "nrel_csm_mass\n",
            "  blade\n",
            "    rotor_diameter                         [126.]              m          \n",
            "    blade_mass_coeff                       [0.5]               None       \n",
            "    blade_user_exp                         [2.5]               None       \n",
            "  hub\n",
            "    blade_mass                             [15751.48043042]    kg         \n",
            "    hub_mass_coeff                         [2.3]               None       \n",
            "    hub_mass_intercept                     [1320.]             None       \n",
            "  pitch\n",
            "    blade_mass                             [15751.48043042]    kg         \n",
            "    pitch_bearing_mass_coeff               [0.1295]            None       \n",
            "    pitch_bearing_mass_intercept           [491.31]            None       \n",
            "    bearing_housing_fraction               [0.328]             None       \n",
            "    mass_sys_offset                        [555.]              None       \n",
            "  spinner\n",
            "    rotor_diameter                         [126.]              m          \n",
            "    spinner_mass_coeff                     [15.5]              None       \n",
            "    spinner_mass_intercept                 [-980.]             None       \n",
            "  lss\n",
            "    blade_mass                             [15751.48043042]    kg         \n",
            "    machine_rating                         [5000.]             kW         \n",
            "    lss_mass_coeff                         [13.]               None       \n",
            "    lss_mass_exp                           [0.65]              None       \n",
            "    lss_mass_intercept                     [775.]              None       \n",
            "  bearing\n",
            "    rotor_diameter                         [126.]              m          \n",
            "    bearing_mass_coeff                     [0.0001]            None       \n",
            "    bearing_mass_exp                       [3.5]               None       \n",
            "  torque\n",
            "    rotor_diameter                         [126.]              m          \n",
            "    machine_rating                         [5000.]             kW         \n",
            "    max_tip_speed                          [80.]               m/s        \n",
            "    max_efficiency                         [0.9]               None       \n",
            "  gearbox\n",
            "    rotor_torque                           [4375000.]          N*m        \n",
            "    gearbox_mass_coeff                     [113.]              None       \n",
            "    gearbox_mass_exp                       [0.71]              None       \n",
            "  hss\n",
            "    machine_rating                         [5000.]             kW         \n",
            "    hss_mass_coeff                         [0.19894]           None       \n",
            "  brake\n",
            "    rotor_torque                           [4375000.]          N*m        \n",
            "    brake_mass_coeff                       [0.00122]           None       \n",
            "  generator\n",
            "    machine_rating                         [5000.]             kW         \n",
            "    generator_mass_coeff                   [2300.]             None       \n",
            "    generator_mass_intercept               [3400.]             None       \n",
            "  bedplate\n",
            "    rotor_diameter                         [126.]              m          \n",
            "    bedplate_mass_exp                      [2.2]               None       \n",
            "  yaw\n",
            "    rotor_diameter                         [126.]              m          \n",
            "    yaw_mass_coeff                         [0.0009]            None       \n",
            "    yaw_mass_exp                           [3.314]             None       \n",
            "  hvac\n",
            "    machine_rating                         [5000.]             kW         \n",
            "    hvac_mass_coeff                        [0.08]              None       \n",
            "  cover\n",
            "    machine_rating                         [5000.]             kW         \n",
            "    cover_mass_coeff                       [1.2817]            None       \n",
            "    cover_mass_intercept                   [428.19]            None       \n",
            "  platforms\n",
            "    bedplate_mass                          [41765.26095285]    kg         \n",
            "    platforms_mass_coeff                   [0.125]             None       \n",
            "    crane_weight                           [3000.]             kg         \n",
            "  transformer\n",
            "    machine_rating                         [5000.]             kW         \n",
            "    transformer_mass_coeff                 [1915.]             None       \n",
            "    transformer_mass_intercept             [1910.]             None       \n",
            "  tower\n",
            "    hub_height                             [90.]               m          \n",
            "    tower_mass_coeff                       [19.828]            None       \n",
            "    tower_mass_exp                         [2.0282]            None       \n",
            "  turbine\n",
            "    blade_mass                             [15751.48043042]    kg         \n",
            "    hub_mass                               [37548.40498997]    kg         \n",
            "    pitch_system_mass                      [9334.08947551]     kg         \n",
            "    spinner_mass                           [973.]              kg         \n",
            "    lss_mass                               [20568.96284886]    kg         \n",
            "    main_bearing_mass                      [2245.41649102]     kg         \n",
            "    gearbox_mass                           [43468.32086769]    kg         \n",
            "    hss_mass                               [994.7]             kg         \n",
            "    brake_mass                             [5337.5]            kg         \n",
            "    generator_mass                         [14900.]            kg         \n",
            "    bedplate_mass                          [41765.26095285]    kg         \n",
            "    yaw_mass                               [12329.96247921]    kg         \n",
            "    hvac_mass                              [400.]              kg         \n",
            "    cover_mass                             [6836.69]           kg         \n",
            "    platforms_mass                         [8220.65761911]     kg         \n",
            "    transformer_mass                       [11485.]            kg         \n",
            "    tower_mass                             [182336.48057717]   kg         \n",
            "  blade\n",
            "    turbine_class                          2                   Unavailable\n",
            "    blade_has_carbon                       False               Unavailable\n",
            "  pitch\n",
            "    blade_number                           3                   Unavailable\n",
            "  platforms\n",
            "    crane                                  True                Unavailable\n",
            "  turbine\n",
            "    blade_number                           3                   Unavailable\n",
            "    main_bearing_number                    2                   Unavailable\n",
            "turbine_costs\n",
            "  blade_c\n",
            "    blade_mass                             [15751.48043042]    kg         \n",
            "    blade_mass_cost_coeff                  [14.6]              USD/kg     \n",
            "    blade_cost_external                    [0.]                USD        \n",
            "  hub_c\n",
            "    hub_mass                               [37548.40498997]    kg         \n",
            "    hub_mass_cost_coeff                    [3.9]               USD/kg     \n",
            "  pitch_c\n",
            "    pitch_system_mass                      [9334.08947551]     kg         \n",
            "    pitch_system_mass_cost_coeff           [22.1]              USD/kg     \n",
            "  spinner_c\n",
            "    spinner_mass                           [973.]              kg         \n",
            "    spinner_mass_cost_coeff                [11.1]              USD/kg     \n",
            "  hub_adder\n",
            "    hub_cost                               [146438.77946089]   USD        \n",
            "    hub_mass                               [37548.40498997]    kg         \n",
            "    pitch_system_cost                      [206283.37740871]   USD        \n",
            "    pitch_system_mass                      [9334.08947551]     kg         \n",
            "    spinner_cost                           [10800.3]           USD        \n",
            "    spinner_mass                           [973.]              kg         \n",
            "    hub_assemblyCostMultiplier             [0.]                None       \n",
            "    hub_overheadCostMultiplier             [0.]                None       \n",
            "    hub_profitMultiplier                   [0.]                None       \n",
            "    hub_transportMultiplier                [0.]                None       \n",
            "  rotor_adder\n",
            "    blade_cost                             [229971.61428418]   USD        \n",
            "    blade_mass                             [15751.48043042]    kg         \n",
            "    hub_system_cost                        [363522.45686961]   USD        \n",
            "    hub_system_mass_tcc                    [47855.49446548]    kg         \n",
            "  lss_c\n",
            "    lss_mass                               [20568.96284886]    kg         \n",
            "    lss_mass_cost_coeff                    [11.9]              USD/kg     \n",
            "  bearing_c\n",
            "    main_bearing_mass                      [2245.41649102]     kg         \n",
            "    bearing_mass_cost_coeff                [4.5]               USD/kg     \n",
            "  gearbox_c\n",
            "    gearbox_mass                           [43468.32086769]    kg         \n",
            "    gearbox_mass_cost_coeff                [12.9]              USD/kg     \n",
            "  hss_c\n",
            "    hss_mass                               [994.7]             kg         \n",
            "    hss_mass_cost_coeff                    [6.8]               USD/kg     \n",
            "  brake_c\n",
            "    brake_mass                             [5337.5]            kg         \n",
            "    brake_mass_cost_coeff                  [3.6254]            USD/kg     \n",
            "  generator_c\n",
            "    generator_mass                         [14900.]            kg         \n",
            "    generator_mass_cost_coeff              [12.4]              USD/kg     \n",
            "    generator_cost_external                [0.]                USD        \n",
            "  bedplate_c\n",
            "    bedplate_mass                          [41765.26095285]    kg         \n",
            "    bedplate_mass_cost_coeff               [2.9]               USD/kg     \n",
            "  yaw_c\n",
            "    yaw_mass                               [12329.96247921]    kg         \n",
            "    yaw_mass_cost_coeff                    [8.3]               USD/kg     \n",
            "  hvac_c\n",
            "    hvac_mass                              [400.]              kg         \n",
            "    hvac_mass_cost_coeff                   [124.]              USD/kg     \n",
            "  controls_c\n",
            "    machine_rating                         [5000.]             kW         \n",
            "    controls_machine_rating_cost_coeff     [21.15]             USD/kW     \n",
            "  converter_c\n",
            "    converter_mass                         [0.]                kg         \n",
            "    converter_mass_cost_coeff              [18.8]              USD/kg     \n",
            "  elec_c\n",
            "    machine_rating                         [5000.]             kW         \n",
            "    elec_connec_machine_rating_cost_coeff  [41.85]             USD/kW     \n",
            "  cover_c\n",
            "    cover_mass                             [6836.69]           kg         \n",
            "    cover_mass_cost_coeff                  [5.7]               USD/kg     \n",
            "  platforms_c\n",
            "    platforms_mass                         [8220.65761911]     kg         \n",
            "    platforms_mass_cost_coeff              [17.1]              USD/kg     \n",
            "    crane_cost                             [12000.]            USD        \n",
            "  transformer_c\n",
            "    transformer_mass                       [11485.]            kg         \n",
            "    transformer_mass_cost_coeff            [18.8]              USD/kg     \n",
            "  nacelle_adder\n",
            "    lss_cost                               [244770.65790145]   USD        \n",
            "    lss_mass                               [20568.96284886]    kg         \n",
            "    main_bearing_cost                      [10104.37420958]    USD        \n",
            "    main_bearing_mass                      [2245.41649102]     kg         \n",
            "    gearbox_cost                           [560741.33919321]   USD        \n",
            "    gearbox_mass                           [43468.32086769]    kg         \n",
            "    hss_cost                               [6763.96]           USD        \n",
            "    hss_mass                               [994.7]             kg         \n",
            "    brake_cost                             [19350.5725]        USD        \n",
            "    brake_mass                             [5337.5]            kg         \n",
            "    generator_cost                         [184760.]           USD        \n",
            "    generator_mass                         [14900.]            kg         \n",
            "    bedplate_cost                          [121119.25676326]   USD        \n",
            "    bedplate_mass                          [41765.26095285]    kg         \n",
            "    yaw_system_cost                        [102338.68857747]   USD        \n",
            "    yaw_mass                               [12329.96247921]    kg         \n",
            "    converter_cost                         [0.]                USD        \n",
            "    converter_mass                         [0.]                kg         \n",
            "    hvac_cost                              [49600.]            USD        \n",
            "    hvac_mass                              [400.]              kg         \n",
            "    cover_cost                             [38969.133]         USD        \n",
            "    cover_mass                             [6836.69]           kg         \n",
            "    elec_cost                              [209250.]           USD        \n",
            "    controls_cost                          [105750.]           USD        \n",
            "    platforms_mass                         [8220.65761911]     kg         \n",
            "    platforms_cost                         [101273.24528671]   USD        \n",
            "    transformer_cost                       [215918.]           USD        \n",
            "    transformer_mass                       [11485.]            kg         \n",
            "    nacelle_assemblyCostMultiplier         [0.]                None       \n",
            "    nacelle_overheadCostMultiplier         [0.]                None       \n",
            "    nacelle_profitMultiplier               [0.]                None       \n",
            "    nacelle_transportMultiplier            [0.]                None       \n",
            "  tower_c\n",
            "    tower_mass                             [182336.48057717]   kg         \n",
            "    tower_mass_cost_coeff                  [2.9]               USD/kg     \n",
            "    tower_cost_external                    [0.]                USD        \n",
            "  tower_adder\n",
            "    tower_parts_cost                       [528775.7936738]    USD        \n",
            "    tower_assemblyCostMultiplier           [0.]                None       \n",
            "    tower_overheadCostMultiplier           [0.]                None       \n",
            "    tower_profitMultiplier                 [0.]                None       \n",
            "    tower_transportMultiplier              [0.]                None       \n",
            "  turbine_c\n",
            "    rotor_cost                             [1053437.29972214]  USD        \n",
            "    rotor_mass_tcc                         [95109.93575675]    kg         \n",
            "    nacelle_cost                           [1980813.60164125]  USD        \n",
            "    nacelle_mass_tcc                       [170797.88774975]   kg         \n",
            "    tower_cost                             [528775.7936738]    USD        \n",
            "    tower_mass                             [182336.48057717]   kg         \n",
            "    machine_rating                         [5000.]             kW         \n",
            "    turbine_assemblyCostMultiplier         [0.]                None       \n",
            "    turbine_overheadCostMultiplier         [0.]                None       \n",
            "    turbine_profitMultiplier               [0.]                None       \n",
            "    turbine_transportMultiplier            [0.]                None       \n",
            "  outputs\n",
            "    blade_cost                             [229971.61428418]   USD        \n",
            "    blade_mass                             [15751.48043042]    kg         \n",
            "    hub_cost                               [146438.77946089]   USD        \n",
            "    hub_mass                               [37548.40498997]    kg         \n",
            "    pitch_system_cost                      [206283.37740871]   USD        \n",
            "    pitch_system_mass                      [9334.08947551]     kg         \n",
            "    spinner_cost                           [10800.3]           USD        \n",
            "    spinner_mass                           [973.]              kg         \n",
            "    lss_cost                               [244770.65790145]   USD        \n",
            "    lss_mass                               [20568.96284886]    kg         \n",
            "    main_bearing_cost                      [10104.37420958]    USD        \n",
            "    main_bearing_mass                      [2245.41649102]     kg         \n",
            "    gearbox_cost                           [560741.33919321]   USD        \n",
            "    gearbox_mass                           [43468.32086769]    kg         \n",
            "    hss_cost                               [6763.96]           USD        \n",
            "    hss_mass                               [994.7]             kg         \n",
            "    brake_cost                             [19350.5725]        USD        \n",
            "    brake_mass                             [5337.5]            kg         \n",
            "    generator_cost                         [184760.]           USD        \n",
            "    generator_mass                         [14900.]            kg         \n",
            "    bedplate_cost                          [121119.25676326]   USD        \n",
            "    bedplate_mass                          [41765.26095285]    kg         \n",
            "    yaw_system_cost                        [102338.68857747]   USD        \n",
            "    yaw_mass                               [12329.96247921]    kg         \n",
            "    hvac_cost                              [49600.]            USD        \n",
            "    hvac_mass                              [400.]              kg         \n",
            "    cover_cost                             [38969.133]         USD        \n",
            "    cover_mass                             [6836.69]           kg         \n",
            "    elec_cost                              [209250.]           USD        \n",
            "    controls_cost                          [105750.]           USD        \n",
            "    platforms_cost                         [101273.24528671]   USD        \n",
            "    transformer_cost                       [215918.]           USD        \n",
            "    transformer_mass                       [11485.]            kg         \n",
            "    converter_cost                         [0.]                USD        \n",
            "    converter_mass                         [0.]                kg         \n",
            "    rotor_cost                             [1053437.29972214]  USD        \n",
            "    rotor_mass_tcc                         [95109.93575675]    kg         \n",
            "    nacelle_cost                           [1980813.60164125]  USD        \n",
            "    nacelle_mass_tcc                       [170797.88774975]   kg         \n",
            "    tower_cost                             [528775.7936738]    USD        \n",
            "    tower_mass                             [182336.48057717]   kg         \n",
            "    turbine_cost                           [3563026.69503719]  USD        \n",
            "    turbine_cost_kW                        [712.60533901]      USD/kW     \n",
            "    turbine_mass_tcc                       [448244.30408368]   kg         \n",
            "  rotor_adder\n",
            "    blade_number                           3                   Unavailable\n",
            "  platforms_c\n",
            "    crane                                  True                Unavailable\n",
            "  nacelle_adder\n",
            "    main_bearing_number                    2                   Unavailable\n",
            "\n",
            "\n",
            "53 Explicit Output(s) in 'model'\n",
            "\n",
            "varname                  val                 units \n",
            "-----------------------  ------------------  ------\n",
            "nrel_csm_mass\n",
            "  blade\n",
            "    blade_mass           [15751.48043042]    kg    \n",
            "  hub\n",
            "    hub_mass             [37548.40498997]    kg    \n",
            "  pitch\n",
            "    pitch_system_mass    [9334.08947551]     kg    \n",
            "  spinner\n",
            "    spinner_mass         [973.]              kg    \n",
            "  lss\n",
            "    lss_mass             [20568.96284886]    kg    \n",
            "  bearing\n",
            "    main_bearing_mass    [2245.41649102]     kg    \n",
            "  torque\n",
            "    rated_rpm            [12.1260909]        rpm   \n",
            "    rotor_torque         [4375000.]          N*m   \n",
            "  gearbox\n",
            "    gearbox_mass         [43468.32086769]    kg    \n",
            "  hss\n",
            "    hss_mass             [994.7]             kg    \n",
            "  brake\n",
            "    brake_mass           [5337.5]            kg    \n",
            "  generator\n",
            "    generator_mass       [14900.]            kg    \n",
            "  bedplate\n",
            "    bedplate_mass        [41765.26095285]    kg    \n",
            "  yaw\n",
            "    yaw_mass             [12329.96247921]    kg    \n",
            "  hvac\n",
            "    hvac_mass            [400.]              kg    \n",
            "  cover\n",
            "    cover_mass           [6836.69]           kg    \n",
            "  platforms\n",
            "    platforms_mass       [8220.65761911]     kg    \n",
            "  transformer\n",
            "    transformer_mass     [11485.]            kg    \n",
            "  tower\n",
            "    tower_mass           [182336.48057717]   kg    \n",
            "  turbine\n",
            "    hub_system_mass      [47855.49446548]    kg    \n",
            "    rotor_mass           [95109.93575675]    kg    \n",
            "    nacelle_mass         [170797.88774975]   kg    \n",
            "    turbine_mass         [448244.30408368]   kg    \n",
            "turbine_costs\n",
            "  blade_c\n",
            "    blade_cost           [229971.61428418]   USD   \n",
            "  hub_c\n",
            "    hub_cost             [146438.77946089]   USD   \n",
            "  pitch_c\n",
            "    pitch_system_cost    [206283.37740871]   USD   \n",
            "  spinner_c\n",
            "    spinner_cost         [10800.3]           USD   \n",
            "  hub_adder\n",
            "    hub_system_mass_tcc  [47855.49446548]    kg    \n",
            "    hub_system_cost      [363522.45686961]   USD   \n",
            "  rotor_adder\n",
            "    rotor_cost           [1053437.29972214]  USD   \n",
            "    rotor_mass_tcc       [95109.93575675]    kg    \n",
            "  lss_c\n",
            "    lss_cost             [244770.65790145]   USD   \n",
            "  bearing_c\n",
            "    main_bearing_cost    [10104.37420958]    USD   \n",
            "  gearbox_c\n",
            "    gearbox_cost         [560741.33919321]   USD   \n",
            "  hss_c\n",
            "    hss_cost             [6763.96]           USD   \n",
            "  brake_c\n",
            "    brake_cost           [19350.5725]        USD   \n",
            "  generator_c\n",
            "    generator_cost       [184760.]           USD   \n",
            "  bedplate_c\n",
            "    bedplate_cost        [121119.25676326]   USD   \n",
            "  yaw_c\n",
            "    yaw_system_cost      [102338.68857747]   USD   \n",
            "  hvac_c\n",
            "    hvac_cost            [49600.]            USD   \n",
            "  controls_c\n",
            "    controls_cost        [105750.]           USD   \n",
            "  converter_c\n",
            "    converter_cost       [0.]                USD   \n",
            "  elec_c\n",
            "    elec_cost            [209250.]           USD   \n",
            "  cover_c\n",
            "    cover_cost           [38969.133]         USD   \n",
            "  platforms_c\n",
            "    platforms_cost       [101273.24528671]   USD   \n",
            "  transformer_c\n",
            "    transformer_cost     [215918.]           USD   \n",
            "  nacelle_adder\n",
            "    nacelle_cost         [1980813.60164125]  USD   \n",
            "    nacelle_mass_tcc     [170797.88774975]   kg    \n",
            "  tower_c\n",
            "    tower_parts_cost     [528775.7936738]    USD   \n",
            "  tower_adder\n",
            "    tower_cost           [528775.7936738]    USD   \n",
            "  turbine_c\n",
            "    turbine_mass_tcc     [448244.30408368]   kg    \n",
            "    turbine_cost         [3563026.69503719]  USD   \n",
            "    turbine_cost_kW      [712.60533901]      USD/kW\n",
            "\n",
            "\n",
            "0 Implicit Output(s) in 'model'\n",
            "\n",
            "\n"
          ]
        },
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "[('nrel_csm_mass.blade.blade_mass',\n",
              "  {'val': array([15751.48043042]), 'units': 'kg'}),\n",
              " ('nrel_csm_mass.hub.hub_mass',\n",
              "  {'val': array([37548.40498997]), 'units': 'kg'}),\n",
              " ('nrel_csm_mass.pitch.pitch_system_mass',\n",
              "  {'val': array([9334.08947551]), 'units': 'kg'}),\n",
              " ('nrel_csm_mass.spinner.spinner_mass', {'val': array([973.]), 'units': 'kg'}),\n",
              " ('nrel_csm_mass.lss.lss_mass',\n",
              "  {'val': array([20568.96284886]), 'units': 'kg'}),\n",
              " ('nrel_csm_mass.bearing.main_bearing_mass',\n",
              "  {'val': array([2245.41649102]), 'units': 'kg'}),\n",
              " ('nrel_csm_mass.torque.rated_rpm',\n",
              "  {'val': array([12.1260909]), 'units': 'rpm'}),\n",
              " ('nrel_csm_mass.torque.rotor_torque',\n",
              "  {'val': array([4375000.]), 'units': 'N*m'}),\n",
              " ('nrel_csm_mass.gearbox.gearbox_mass',\n",
              "  {'val': array([43468.32086769]), 'units': 'kg'}),\n",
              " ('nrel_csm_mass.hss.hss_mass', {'val': array([994.7]), 'units': 'kg'}),\n",
              " ('nrel_csm_mass.brake.brake_mass', {'val': array([5337.5]), 'units': 'kg'}),\n",
              " ('nrel_csm_mass.generator.generator_mass',\n",
              "  {'val': array([14900.]), 'units': 'kg'}),\n",
              " ('nrel_csm_mass.bedplate.bedplate_mass',\n",
              "  {'val': array([41765.26095285]), 'units': 'kg'}),\n",
              " ('nrel_csm_mass.yaw.yaw_mass',\n",
              "  {'val': array([12329.96247921]), 'units': 'kg'}),\n",
              " ('nrel_csm_mass.hvac.hvac_mass', {'val': array([400.]), 'units': 'kg'}),\n",
              " ('nrel_csm_mass.cover.cover_mass', {'val': array([6836.69]), 'units': 'kg'}),\n",
              " ('nrel_csm_mass.platforms.platforms_mass',\n",
              "  {'val': array([8220.65761911]), 'units': 'kg'}),\n",
              " ('nrel_csm_mass.transformer.transformer_mass',\n",
              "  {'val': array([11485.]), 'units': 'kg'}),\n",
              " ('nrel_csm_mass.tower.tower_mass',\n",
              "  {'val': array([182336.48057717]), 'units': 'kg'}),\n",
              " ('nrel_csm_mass.turbine.hub_system_mass',\n",
              "  {'val': array([47855.49446548]), 'units': 'kg'}),\n",
              " ('nrel_csm_mass.turbine.rotor_mass',\n",
              "  {'val': array([95109.93575675]), 'units': 'kg'}),\n",
              " ('nrel_csm_mass.turbine.nacelle_mass',\n",
              "  {'val': array([170797.88774975]), 'units': 'kg'}),\n",
              " ('nrel_csm_mass.turbine.turbine_mass',\n",
              "  {'val': array([448244.30408368]), 'units': 'kg'}),\n",
              " ('turbine_costs.blade_c.blade_cost',\n",
              "  {'val': array([229971.61428418]), 'units': 'USD'}),\n",
              " ('turbine_costs.hub_c.hub_cost',\n",
              "  {'val': array([146438.77946089]), 'units': 'USD'}),\n",
              " ('turbine_costs.pitch_c.pitch_system_cost',\n",
              "  {'val': array([206283.37740871]), 'units': 'USD'}),\n",
              " ('turbine_costs.spinner_c.spinner_cost',\n",
              "  {'val': array([10800.3]), 'units': 'USD'}),\n",
              " ('turbine_costs.hub_adder.hub_system_mass_tcc',\n",
              "  {'val': array([47855.49446548]), 'units': 'kg'}),\n",
              " ('turbine_costs.hub_adder.hub_system_cost',\n",
              "  {'val': array([363522.45686961]), 'units': 'USD'}),\n",
              " ('turbine_costs.rotor_adder.rotor_cost',\n",
              "  {'val': array([1053437.29972214]), 'units': 'USD'}),\n",
              " ('turbine_costs.rotor_adder.rotor_mass_tcc',\n",
              "  {'val': array([95109.93575675]), 'units': 'kg'}),\n",
              " ('turbine_costs.lss_c.lss_cost',\n",
              "  {'val': array([244770.65790145]), 'units': 'USD'}),\n",
              " ('turbine_costs.bearing_c.main_bearing_cost',\n",
              "  {'val': array([10104.37420958]), 'units': 'USD'}),\n",
              " ('turbine_costs.gearbox_c.gearbox_cost',\n",
              "  {'val': array([560741.33919321]), 'units': 'USD'}),\n",
              " ('turbine_costs.hss_c.hss_cost', {'val': array([6763.96]), 'units': 'USD'}),\n",
              " ('turbine_costs.brake_c.brake_cost',\n",
              "  {'val': array([19350.5725]), 'units': 'USD'}),\n",
              " ('turbine_costs.generator_c.generator_cost',\n",
              "  {'val': array([184760.]), 'units': 'USD'}),\n",
              " ('turbine_costs.bedplate_c.bedplate_cost',\n",
              "  {'val': array([121119.25676326]), 'units': 'USD'}),\n",
              " ('turbine_costs.yaw_c.yaw_system_cost',\n",
              "  {'val': array([102338.68857747]), 'units': 'USD'}),\n",
              " ('turbine_costs.hvac_c.hvac_cost', {'val': array([49600.]), 'units': 'USD'}),\n",
              " ('turbine_costs.controls_c.controls_cost',\n",
              "  {'val': array([105750.]), 'units': 'USD'}),\n",
              " ('turbine_costs.converter_c.converter_cost',\n",
              "  {'val': array([0.]), 'units': 'USD'}),\n",
              " ('turbine_costs.elec_c.elec_cost', {'val': array([209250.]), 'units': 'USD'}),\n",
              " ('turbine_costs.cover_c.cover_cost',\n",
              "  {'val': array([38969.133]), 'units': 'USD'}),\n",
              " ('turbine_costs.platforms_c.platforms_cost',\n",
              "  {'val': array([101273.24528671]), 'units': 'USD'}),\n",
              " ('turbine_costs.transformer_c.transformer_cost',\n",
              "  {'val': array([215918.]), 'units': 'USD'}),\n",
              " ('turbine_costs.nacelle_adder.nacelle_cost',\n",
              "  {'val': array([1980813.60164125]), 'units': 'USD'}),\n",
              " ('turbine_costs.nacelle_adder.nacelle_mass_tcc',\n",
              "  {'val': array([170797.88774975]), 'units': 'kg'}),\n",
              " ('turbine_costs.tower_c.tower_parts_cost',\n",
              "  {'val': array([528775.7936738]), 'units': 'USD'}),\n",
              " ('turbine_costs.tower_adder.tower_cost',\n",
              "  {'val': array([528775.7936738]), 'units': 'USD'}),\n",
              " ('turbine_costs.turbine_c.turbine_mass_tcc',\n",
              "  {'val': array([448244.30408368]), 'units': 'kg'}),\n",
              " ('turbine_costs.turbine_c.turbine_cost',\n",
              "  {'val': array([3563026.69503719]), 'units': 'USD'}),\n",
              " ('turbine_costs.turbine_c.turbine_cost_kW',\n",
              "  {'val': array([712.60533901]), 'units': 'USD/kW'})]"
            ]
          },
          "metadata": {},
          "execution_count": 9
        }
      ]
    },
    {
      "cell_type": "code",
      "source": [
        "prob['turbine_costs.turbine_c.turbine_cost']"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "iWfAMWaeWsNi",
        "outputId": "92b8032c-2df6-4537-bb9d-1f84da36a6b7"
      },
      "execution_count": 10,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "array([3563026.69503719])"
            ]
          },
          "metadata": {},
          "execution_count": 10
        }
      ]
    },
    {
      "cell_type": "code",
      "source": [
        "from wisdem.nrelcsm.nrel_csm_mass_2015 import nrel_csm_2015\n",
        "def nrelcsm_wrapper(rated_power=5000, rotor_diameter=126):\n",
        "  \n",
        "  # 0 ---------- (marker for docs)\n",
        "\n",
        "  # 1 ---------- (marker for docs)\n",
        "  # OpenMDAO Problem instance\n",
        "  prob = om.Problem(reports=False)\n",
        "  prob.model = nrel_csm_2015()\n",
        "  prob.setup()\n",
        "  # 1 ---------- (marker for docs)\n",
        "\n",
        "  # 2 ---------- (marker for docs)\n",
        "  # Initialize variables for NREL CSM\n",
        "  prob[\"machine_rating\"] = rated_power\n",
        "  prob[\"rotor_diameter\"] = rotor_diameter\n",
        "  prob[\"turbine_class\"] = 2\n",
        "  prob[\"hub_height\"] = 90.0\n",
        "  prob[\"blade_number\"] = 3\n",
        "  prob[\"blade_has_carbon\"] = False\n",
        "  prob[\"max_tip_speed\"] = 80.0\n",
        "  prob[\"max_efficiency\"] = 0.90\n",
        "  prob[\"main_bearing_number\"] = 2\n",
        "  prob[\"crane\"] = True\n",
        "  # 2 ---------- (marker for docs)\n",
        "\n",
        "  # 3 ---------- (marker for docs)\n",
        "  # Evaluate the model\n",
        "  prob.run_model()\n",
        "  # 3 ---------- (marker for docs)\n",
        "\n",
        "  # 4 ---------- (marker for docs)\n",
        "  # Print all intermediate inputs and outputs to the screen\n",
        "   # 4 ---------- (marker for docs)\n",
        "  return prob['turbine_costs.turbine_c.turbine_cost']"
      ],
      "metadata": {
        "id": "8pMX8PcuXFhZ"
      },
      "execution_count": 13,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": [
        "nrelcsm_wrapper(rated_power=5000, rotor_diameter=126)"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "PawQL_J3YG8p",
        "outputId": "6cc7f60f-7883-4bd1-e624-b46506e83423"
      },
      "execution_count": 14,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "array([3563026.69503719])"
            ]
          },
          "metadata": {},
          "execution_count": 14
        }
      ]
    },
    {
      "cell_type": "code",
      "source": [],
      "metadata": {
        "id": "dcrpl6FIYJs5"
      },
      "execution_count": null,
      "outputs": []
    }
  ]
}