search

WISDEM / WEIS

Wind Energy with Integrated Servo-controls Toolset
https://weis.readthedocs.io/en/latest/
Apache License 2.0
49 stars 38 forks source link

Warnings when running example 2 #207

Closed rebeccamccabe closed 7 months ago

rebeccamccabe commented 1 year ago

Description

I am new to WEIS and trying to run my first example. I am running 02_control_opt/weis_driver.py since it is on the list of reliable examples as suggested in #198. I get various warnings (related to mpi4py, Rosco components, OpenFAST NaNs, invalid values, AeroDyn nodal output not found, FAST UA and Mach number, and no optimization report) and an LCOE of zero, that makes me think something with my setup/install is wrong. Are these warnings expected?

Steps to reproduce issue

Please provide a minimum working example (MWE) if possible

  1. Clean install of main on WSL following README instructions, with the exception of using numpy=1.23 to get around #202
  2. cd WEIS/examples/02_control opt, python weis_driver.py

Current behavior

(weis-env) rmccabe@EN-MA-AMH2498-F:~/WEIS/examples/02_control_opt$ python weis_driver.py
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/pyoptsparse/pyOpt_MPI.py:68: UserWarning: mpi4py could not be imported. mpi4py is required to use the parallel gradient analysis and parallel objective analysis for non-gradient based optimizers. Continuing using a dummy MPI module from pyOptSparse.
  warnings.warn(warn)
Using weis.aeroelasticse in ROSCO_toolbox...
WARNING: Airfoils provided but no blades/rotor found or RotorSE deactivated
WARNING: Environment provided but no related component found found
Found existing potential model: examples/01_aeroelasticse/OpenFAST_models/IEA-15-240-RWT/IEA-15-240-RWT-UMaineSemi/HydroData/IEA-15-240-RWT-UMaineSemi
    - Trying to use this instead of running PyHAMS.
-----------------------------------------------------------------------------
Loading wind turbine data for NREL's ROSCO tuning and simulation processeses
-----------------------------------------------------------------------------
Loading FAST model: IEA-15-240-RWT-UMaineSemi.fst
Loading rotor performace data from text file: /home/rmccabe/WEIS/examples/01_aeroelasticse/OpenFAST_models/IEA-15-240-RWT/IEA-15-240-RWT/Cp_Ct_Cq.IEA15MW.txt
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/n2_viewer/n2_viewer.py:168: OpenMDAOWarning:All-NaN slice encountered
Driver debug print for iter coord: rank0:NLopt_LN_COBYLA|0
----------------------------------------------------------
Design Vars
{'tune_rosco_ivc.Kp_float': array([-10.]),
 'tune_rosco_ivc.omega_pc': array([0.2, 0.2]),
 'tune_rosco_ivc.ptfm_freq': array([0.2]),
 'tune_rosco_ivc.zeta_pc': array([1., 1.])}

Warning: turbine_costsse_2015.py : TowerCost2015.compute : No tower mass provided.  Assuming $0 for tower cost, LCOE will be artificially low.
RuntimeWarning: /home/rmccabe/WEIS/weis/control/dac.py:734
invalid value encountered in double_scalars-----------------------------------------------------------------------------
   Tuning a reference wind turbine controller using NREL's ROSCO toolbox
-----------------------------------------------------------------------------
Loading rotor performace data from text file: /home/rmccabe/WEIS/examples/01_aeroelasticse/OpenFAST_models/IEA-15-240-RWT/IEA-15-240-RWT/Cp_Ct_Cq.IEA15MW.txt

 **************************************************************************************************
 TurbSim

 Copyright (C) 2023 National Renewable Energy Laboratory
 Copyright (C) 2023 Envision Energy USA LTD

 This program is licensed under Apache License Version 2.0 and comes with ABSOLUTELY NO WARRANTY.
 See the "LICENSE" file distributed with this software for details.
 **************************************************************************************************

  Running TurbSim (from OpenFAST-v1.1)

  Reading the input file "IEA15_NTM_U16.000000_Seed1501552846.0.in".
  Calculating the spectral and transfer function matrices:
     u-component matrices (IEC coherence method)
     v-component matrices (identity coherence)
     w-component matrices (identity coherence)
  Generating time series for all points:
     u-component
     v-component
     w-component
  Computing hub-height statistics
  Writing statistics to summary file
  Generating AeroDyn binary time-series file "IEA15_NTM_U16.000000_Seed1501552846.0.bts"

  Processing complete.  3.5469 CPU seconds used.

  TurbSim terminated normally.

 **************************************************************************************************
 TurbSim

 Copyright (C) 2023 National Renewable Energy Laboratory
 Copyright (C) 2023 Envision Energy USA LTD

 This program is licensed under Apache License Version 2.0 and comes with ABSOLUTELY NO WARRANTY.
 See the "LICENSE" file distributed with this software for details.
 **************************************************************************************************

  Running TurbSim (from OpenFAST-v1.1)

  Reading the input file "IEA15_NTM_U46.789000_Seed488200390.0.in".
  Calculating the spectral and transfer function matrices:
     u-component matrices (IEC coherence method)
     v-component matrices (identity coherence)
     w-component matrices (identity coherence)
  Generating time series for all points:
     u-component
     v-component
     w-component
  Computing hub-height statistics
  Writing statistics to summary file
  Generating AeroDyn binary time-series file "IEA15_NTM_U46.789000_Seed488200390.0.bts"

  Processing complete.  8.5781 CPU seconds used.

  TurbSim terminated normally.

 **************************************************************************************************
 TurbSim

 Copyright (C) 2023 National Renewable Energy Laboratory
 Copyright (C) 2023 Envision Energy USA LTD

 This program is licensed under Apache License Version 2.0 and comes with ABSOLUTELY NO WARRANTY.
 See the "LICENSE" file distributed with this software for details.
 **************************************************************************************************

  Running TurbSim (from OpenFAST-v1.1)

  Reading the input file "IEA15_NTM_U46.789000_Seed1693606511.0.in".
  Calculating the spectral and transfer function matrices:
     u-component matrices (IEC coherence method)
     v-component matrices (identity coherence)
     w-component matrices (identity coherence)
  Generating time series for all points:
     u-component
     v-component
     w-component
  Computing hub-height statistics
  Writing statistics to summary file
  Generating AeroDyn binary time-series file "IEA15_NTM_U46.789000_Seed1693606511.0.bts"

  Processing complete.  8.5625 CPU seconds used.

  TurbSim terminated normally.

Writing rotor performance text file: /home/rmccabe/WEIS/examples/02_control_opt/outputs/02_control_opt/IEA15_0_Cp_Ct_Cq.txt
Writing new controller parameter file parameter file: /home/rmccabe/WEIS/examples/02_control_opt/outputs/02_control_opt/IEA15_0_DISCON.IN.

 **************************************************************************************************
 OpenFAST

 Copyright (C) 2023 National Renewable Energy Laboratory
 Copyright (C) 2023 Envision Energy USA LTD

 This program is licensed under Apache License Version 2.0 and comes with ABSOLUTELY NO WARRANTY.
 See the "LICENSE" file distributed with this software for details.
 **************************************************************************************************

 OpenFAST-v1.1
 Compile Info:
  - Compiler: GCC version 9.4.0
  - Architecture: 64 bit
  - Precision: single
  - OpenMP: No
  - Date: Jan 30 2023
  - Time: 10:03:57
 Execution Info:
  - Date: 01/30/2023
  - Time: 10:35:06-0500

 OpenFAST input file heading:
     Generated with AeroElasticSE FAST driver

 Running ElastoDyn.
 Nodal outputs section of ElastoDyn input file not found or improperly formatted.
 Running AeroDyn.
 AD15 Nodal Outputs: Nodal output section of AeroDyn input file not found or improperly formatted.
 Skipping nodal outputs.
 Warning: Turning off Unsteady Aerodynamics because UA parameters are not included in airfoil
 (airfoil has likely has constant polars). (node 1, blade 1)
 Warning: Turning off Unsteady Aerodynamics because UA parameters are not included in airfoil
 (airfoil has likely has constant polars). (node 2, blade 1)
 Warning: Turning off Unsteady Aerodynamics because UA parameters are not included in airfoil
 (airfoil has likely has constant polars). (node 1, blade 2)
 Warning: Turning off Unsteady Aerodynamics because UA parameters are not included in airfoil
 (airfoil has likely has constant polars). (node 2, blade 2)
 Warning: Turning off Unsteady Aerodynamics because UA parameters are not included in airfoil
 (airfoil has likely has constant polars). (node 1, blade 3)
 Warning: Turning off Unsteady Aerodynamics because UA parameters are not included in airfoil
 (airfoil has likely has constant polars). (node 2, blade 3)
 Running InflowWind.

    Reading a 25x25 grid (297.47 m wide, 4.88281E-04 m to 297.47 m above ground) with a
    characteristic wind speed of 16.149 m/s. This full-field file was generated by TurbSim (from
    OpenFAST-v1.1) on 30-Jan-2023 at 10:34:43.

    Processed 204 time steps of 20-Hz full-field data (period of 10.2 seconds).
 Running HydroDyn.
  WARNING:   The random number generator in use differs from the original code provided by NREL.
  This pRNG uses 8 seeds instead of the 2 in the input file.
  Generating incident wave kinematics and current time history.
  Calculating second order sum frequency wave kinematics.
  Reading in WAMIT output with root name
  "/home/rmccabe/WEIS/examples/01_aeroelasticse/OpenFAST_models/IEA-15-240-RWT/IEA-15-240-RWT-UMain
  eSemi/HydroData/IEA-15-240-RWT-UMaineSemi".
  Computing radiation impulse response functions and wave diffraction forces.
 Running MoorDyn.
   MD_Init: Opening MoorDyn input file:  ./IEA15_0_MoorDyn.dat
   Creating mooring system.   3 fairleads, 3 anchors, 0 connects.
    Finalizing ICs using dynamic relaxation.

   t=60  FairTen 1: 2.43587E+06, 2.44120E+06, 2.43677E+06                                             Fairlead tensions did not converge within TMaxIC=60 seconds.
 Running ServoDyn.
 Running ServoDyn Interface for Bladed Controllers (using GNU Fortran for Linux).
 Using legacy Bladed DLL interface.

 FAST_InitializeAll:HydroDyn_Init:Waves_Init:VariousWaves_Init: The random number generator in use
 differs from the original code provided by NREL. This pRNG uses 8 seeds instead of the 2 in the
 HydroDyn input file.

  Time: 0 of 10.1 seconds.

------------------------------------------------------------------------------
Running ROSCO-v2.6.0
A wind turbine controller framework for public use in the scientific field
Developed in collaboration: National Renewable Energy Laboratory
                            Delft University of Technology, The Netherlands
------------------------------------------------------------------------------
Generator speed:    7.6 RPM, Pitch angle:  15.0 deg, Power:     0.0 kW, Est. wind Speed:  15.9 m/s

 FAST_Solution:CalcOutputs_And_SolveForInputs:SolveOption2:RotCalcOutput:BEMT_CalcOutput(node 3,
 blade 3):UA_CalcOutput:UA_BlendSteady:Temporarily turning off UA due to high angle of attack or
 low relative velocity. This warning will not be repeated though the condition may persist.

 FAST_Solution:CalcOutputs_And_SolveForInputs:SolveOption2:RotCalcOutput:BEMT_CalcOutput(node 50,
 blade 3):UA_CalcOutput:Mach number exceeds 0.3. Theory is invalid. This warning will not be
 repeated though the condition may persist.

Generator speed:    8.1 RPM, Pitch angle:   6.4 deg, Power: 16087.0 kW, Est. wind Speed:  13.1 m/s

  Total Real Time:       68.817 seconds
  Total CPU Time:        22.016 seconds
  Simulation CPU Time:   9.875 seconds
  Simulated Time:        10.1 seconds
  Time Ratio (Sim/CPU):  1.0228

  OpenFAST terminated normally.

Runtime:        IEA15_0.fst = 69.09 s
Writing rotor performance text file: /home/rmccabe/WEIS/examples/02_control_opt/outputs/02_control_opt/IEA15_1_Cp_Ct_Cq.txt
Writing new controller parameter file parameter file: /home/rmccabe/WEIS/examples/02_control_opt/outputs/02_control_opt/IEA15_1_DISCON.IN.

 **************************************************************************************************
 OpenFAST

 Copyright (C) 2023 National Renewable Energy Laboratory
 Copyright (C) 2023 Envision Energy USA LTD

 This program is licensed under Apache License Version 2.0 and comes with ABSOLUTELY NO WARRANTY.
 See the "LICENSE" file distributed with this software for details.
 **************************************************************************************************

 OpenFAST-v1.1
 Compile Info:
  - Compiler: GCC version 9.4.0
  - Architecture: 64 bit
  - Precision: single
  - OpenMP: No
  - Date: Jan 30 2023
  - Time: 10:03:57
 Execution Info:
  - Date: 01/30/2023
  - Time: 10:36:16-0500

 OpenFAST input file heading:
     Generated with AeroElasticSE FAST driver

 Running ElastoDyn.
 Nodal outputs section of ElastoDyn input file not found or improperly formatted.
 Running AeroDyn.
 AD15 Nodal Outputs: Nodal output section of AeroDyn input file not found or improperly formatted.
 Skipping nodal outputs.
 Running InflowWind.

    Reading a 25x25 grid (297.47 m wide, 4.88281E-04 m to 297.47 m above ground) with a
    characteristic wind speed of 47.094 m/s. This full-field file was generated by TurbSim (from
    OpenFAST-v1.1) on 30-Jan-2023 at 10:34:47.

    Processed 204 time steps of 20-Hz full-field data (period of 10.2 seconds).
 Running HydroDyn.
  WARNING:   The random number generator in use differs from the original code provided by NREL.
  This pRNG uses 8 seeds instead of the 2 in the input file.
  Generating incident wave kinematics and current time history.
  Calculating second order sum frequency wave kinematics.
  Reading in WAMIT output with root name
  "/home/rmccabe/WEIS/examples/01_aeroelasticse/OpenFAST_models/IEA-15-240-RWT/IEA-15-240-RWT-UMain
  eSemi/HydroData/IEA-15-240-RWT-UMaineSemi".
  Computing radiation impulse response functions and wave diffraction forces.
 Running MoorDyn.
   MD_Init: Opening MoorDyn input file:  ./IEA15_1_MoorDyn.dat
   Creating mooring system.   3 fairleads, 3 anchors, 0 connects.
    Finalizing ICs using dynamic relaxation.

   t=60  FairTen 1: 2.43587E+06, 2.44120E+06, 2.43677E+06                                             Fairlead tensions did not converge within TMaxIC=60 seconds.
 Running ServoDyn.
 Running ServoDyn Interface for Bladed Controllers (using GNU Fortran for Linux).
 Using legacy Bladed DLL interface.

 FAST_InitializeAll:HydroDyn_Init:Waves_Init:VariousWaves_Init: The random number generator in use
 differs from the original code provided by NREL. This pRNG uses 8 seeds instead of the 2 in the
 HydroDyn input file.

  Time: 0 of 10.1 seconds.

------------------------------------------------------------------------------
Running ROSCO-v2.6.0
A wind turbine controller framework for public use in the scientific field
Developed in collaboration: National Renewable Energy Laboratory
                            Delft University of Technology, The Netherlands
------------------------------------------------------------------------------
Generator speed:    0.0 RPM, Pitch angle:  90.0 deg, Power:     0.0 kW, Est. wind Speed:  47.1 m/s
Generator speed:    0.2 RPM, Pitch angle:  70.0 deg, Power:     0.0 kW, Est. wind Speed: -20.6 m/s

  Total Real Time:       18.884 seconds
  Total CPU Time:        18.844 seconds
  Simulation CPU Time:   9.0625 seconds
  Simulated Time:        10.1 seconds
  Time Ratio (Sim/CPU):  1.1145

  OpenFAST terminated normally.

Runtime:        IEA15_1.fst = 18.92 s
Writing rotor performance text file: /home/rmccabe/WEIS/examples/02_control_opt/outputs/02_control_opt/IEA15_2_Cp_Ct_Cq.txt
Writing new controller parameter file parameter file: /home/rmccabe/WEIS/examples/02_control_opt/outputs/02_control_opt/IEA15_2_DISCON.IN.

 **************************************************************************************************
 OpenFAST

 Copyright (C) 2023 National Renewable Energy Laboratory
 Copyright (C) 2023 Envision Energy USA LTD

 This program is licensed under Apache License Version 2.0 and comes with ABSOLUTELY NO WARRANTY.
 See the "LICENSE" file distributed with this software for details.
 **************************************************************************************************

 OpenFAST-v1.1
 Compile Info:
  - Compiler: GCC version 9.4.0
  - Architecture: 64 bit
  - Precision: single
  - OpenMP: No
  - Date: Jan 30 2023
  - Time: 10:03:57
 Execution Info:
  - Date: 01/30/2023
  - Time: 10:36:36-0500

 OpenFAST input file heading:
     Generated with AeroElasticSE FAST driver

 Running ElastoDyn.
 Nodal outputs section of ElastoDyn input file not found or improperly formatted.
 Running AeroDyn.
 AD15 Nodal Outputs: Nodal output section of AeroDyn input file not found or improperly formatted.
 Skipping nodal outputs.
 Running InflowWind.

    Reading a 25x25 grid (297.47 m wide, 4.88281E-04 m to 297.47 m above ground) with a
    characteristic wind speed of 47.094 m/s. This full-field file was generated by TurbSim (from
    OpenFAST-v1.1) on 30-Jan-2023 at 10:34:55.

    Processed 204 time steps of 20-Hz full-field data (period of 10.2 seconds).
 Running HydroDyn.
  WARNING:   The random number generator in use differs from the original code provided by NREL.
  This pRNG uses 8 seeds instead of the 2 in the input file.
  Generating incident wave kinematics and current time history.
  Calculating second order sum frequency wave kinematics.
  Reading in WAMIT output with root name
  "/home/rmccabe/WEIS/examples/01_aeroelasticse/OpenFAST_models/IEA-15-240-RWT/IEA-15-240-RWT-UMain
  eSemi/HydroData/IEA-15-240-RWT-UMaineSemi".
  Computing radiation impulse response functions and wave diffraction forces.
 Running MoorDyn.
   MD_Init: Opening MoorDyn input file:  ./IEA15_2_MoorDyn.dat
   Creating mooring system.   3 fairleads, 3 anchors, 0 connects.
    Finalizing ICs using dynamic relaxation.

   t=60  FairTen 1: 2.43587E+06, 2.44120E+06, 2.43677E+06                                             Fairlead tensions did not converge within TMaxIC=60 seconds.
 Running ServoDyn.
 Running ServoDyn Interface for Bladed Controllers (using GNU Fortran for Linux).
 Using legacy Bladed DLL interface.

 FAST_InitializeAll:HydroDyn_Init:Waves_Init:VariousWaves_Init: The random number generator in use
 differs from the original code provided by NREL. This pRNG uses 8 seeds instead of the 2 in the
 HydroDyn input file.

  Time: 0 of 10.1 seconds.

------------------------------------------------------------------------------
Running ROSCO-v2.6.0
A wind turbine controller framework for public use in the scientific field
Developed in collaboration: National Renewable Energy Laboratory
                            Delft University of Technology, The Netherlands
------------------------------------------------------------------------------
Generator speed:    0.0 RPM, Pitch angle:  90.0 deg, Power:     0.0 kW, Est. wind Speed:  47.1 m/s
Generator speed:    0.2 RPM, Pitch angle:  70.0 deg, Power:     0.0 kW, Est. wind Speed: -20.6 m/s

  Total Real Time:       19.007 seconds
  Total CPU Time:        18.938 seconds
  Simulation CPU Time:   9.0938 seconds
  Simulated Time:        10.1 seconds
  Time Ratio (Sim/CPU):  1.1107

  OpenFAST terminated normally.

Runtime:        IEA15_2.fst = 19.04 s
WARNING: OpenFAST is run at a single wind speed. AEP cannot be estimated. Using average power instead.
########################################
Objectives
Turbine AEP: 0.0150820187 GWh
Blade Mass:  0.0000000000 kg
LCOE:        0.0000000000 USD/MWh
Tip Defl.:   0.0000000000 m
Pitch PI gain inputs: omega_pc[0] = 0.200, zeta_pc[0] = 1.000
Floating Feedback: Kp_float = -10.000, ptfm_freq = 0.200
IPC Ki1p = 0.000e+00
IPC Ki1p = 0.000e+00
DEL(TwrBsMyt): 224454450.0820847154 Nm
rotor_overspeed: 7.5571429760 %
########################################
Nonlinear constraints
{'aeroelastic.rotor_overspeed': array([0.07557143]),
 'sse_tune.tune_rosco.PC_Ki': array([-0.12552134]),
 'sse_tune.tune_rosco.PC_Kp': array([-1.20562287])}

Linear constraints
None

Objectives
{'aeroelastic.DEL_TwrBsMyt': array([224454.45008208])}

Driver debug print for iter coord: rank0:NLopt_LN_COBYLA|1
----------------------------------------------------------
Design Vars
{'tune_rosco_ivc.Kp_float': array([-10.]),
 'tune_rosco_ivc.omega_pc': array([0.2, 0.2]),
 'tune_rosco_ivc.ptfm_freq': array([0.2]),
 'tune_rosco_ivc.zeta_pc': array([1., 1.])}

Warning: turbine_costsse_2015.py : TowerCost2015.compute : No tower mass provided.  Assuming $0 for tower cost, LCOE will be artificially low.
RuntimeWarning: /home/rmccabe/WEIS/weis/control/dac.py:734
invalid value encountered in double_scalars-----------------------------------------------------------------------------
   Tuning a reference wind turbine controller using NREL's ROSCO toolbox
-----------------------------------------------------------------------------
Loading rotor performace data from text file: /home/rmccabe/WEIS/examples/01_aeroelasticse/OpenFAST_models/IEA-15-240-RWT/IEA-15-240-RWT/Cp_Ct_Cq.IEA15MW.txt
Writing rotor performance text file: /home/rmccabe/WEIS/examples/02_control_opt/outputs/02_control_opt/IEA15_0_Cp_Ct_Cq.txt
Writing new controller parameter file parameter file: /home/rmccabe/WEIS/examples/02_control_opt/outputs/02_control_opt/IEA15_0_DISCON.IN.

 **************************************************************************************************
 OpenFAST

 Copyright (C) 2023 National Renewable Energy Laboratory
 Copyright (C) 2023 Envision Energy USA LTD

 This program is licensed under Apache License Version 2.0 and comes with ABSOLUTELY NO WARRANTY.
 See the "LICENSE" file distributed with this software for details.
 **************************************************************************************************

 OpenFAST-v1.1
 Compile Info:
  - Compiler: GCC version 9.4.0
  - Architecture: 64 bit
  - Precision: single
  - OpenMP: No
  - Date: Jan 30 2023
  - Time: 10:03:57
 Execution Info:
  - Date: 01/30/2023
  - Time: 10:36:56-0500

 OpenFAST input file heading:
     Generated with AeroElasticSE FAST driver

 Running ElastoDyn.
 Nodal outputs section of ElastoDyn input file not found or improperly formatted.
 Running AeroDyn.
 AD15 Nodal Outputs: Nodal output section of AeroDyn input file not found or improperly formatted.
 Skipping nodal outputs.
 Warning: Turning off Unsteady Aerodynamics because UA parameters are not included in airfoil
 (airfoil has likely has constant polars). (node 1, blade 1)
 Warning: Turning off Unsteady Aerodynamics because UA parameters are not included in airfoil
 (airfoil has likely has constant polars). (node 2, blade 1)
 Warning: Turning off Unsteady Aerodynamics because UA parameters are not included in airfoil
 (airfoil has likely has constant polars). (node 1, blade 2)
 Warning: Turning off Unsteady Aerodynamics because UA parameters are not included in airfoil
 (airfoil has likely has constant polars). (node 2, blade 2)
 Warning: Turning off Unsteady Aerodynamics because UA parameters are not included in airfoil
 (airfoil has likely has constant polars). (node 1, blade 3)
 Warning: Turning off Unsteady Aerodynamics because UA parameters are not included in airfoil
 (airfoil has likely has constant polars). (node 2, blade 3)
 Running InflowWind.

    Reading a 25x25 grid (297.47 m wide, 4.88281E-04 m to 297.47 m above ground) with a
    characteristic wind speed of 16.149 m/s. This full-field file was generated by TurbSim (from
    OpenFAST-v1.1) on 30-Jan-2023 at 10:34:43.

    Processed 204 time steps of 20-Hz full-field data (period of 10.2 seconds).
 Running HydroDyn.
  WARNING:   The random number generator in use differs from the original code provided by NREL.
  This pRNG uses 8 seeds instead of the 2 in the input file.
  Generating incident wave kinematics and current time history.
  Calculating second order sum frequency wave kinematics.
  Reading in WAMIT output with root name
  "/home/rmccabe/WEIS/examples/01_aeroelasticse/OpenFAST_models/IEA-15-240-RWT/IEA-15-240-RWT-UMain
  eSemi/HydroData/IEA-15-240-RWT-UMaineSemi".
  Computing radiation impulse response functions and wave diffraction forces.
 Running MoorDyn.
   MD_Init: Opening MoorDyn input file:  ./IEA15_0_MoorDyn.dat
   Creating mooring system.   3 fairleads, 3 anchors, 0 connects.
    Finalizing ICs using dynamic relaxation.

   t=60  FairTen 1: 2.43587E+06, 2.44120E+06, 2.43677E+06                                             Fairlead tensions did not converge within TMaxIC=60 seconds.
 Running ServoDyn.
 Running ServoDyn Interface for Bladed Controllers (using GNU Fortran for Linux).
 Using legacy Bladed DLL interface.

 FAST_InitializeAll:HydroDyn_Init:Waves_Init:VariousWaves_Init: The random number generator in use
 differs from the original code provided by NREL. This pRNG uses 8 seeds instead of the 2 in the
 HydroDyn input file.

  Time: 0 of 10.1 seconds.

------------------------------------------------------------------------------
Running ROSCO-v2.6.0
A wind turbine controller framework for public use in the scientific field
Developed in collaboration: National Renewable Energy Laboratory
                            Delft University of Technology, The Netherlands
------------------------------------------------------------------------------
Generator speed:    7.6 RPM, Pitch angle:  15.0 deg, Power:     0.0 kW, Est. wind Speed:  15.9 m/s

 FAST_Solution:CalcOutputs_And_SolveForInputs:SolveOption2:RotCalcOutput:BEMT_CalcOutput(node 3,
 blade 3):UA_CalcOutput:UA_BlendSteady:Temporarily turning off UA due to high angle of attack or
 low relative velocity. This warning will not be repeated though the condition may persist.

 FAST_Solution:CalcOutputs_And_SolveForInputs:SolveOption2:RotCalcOutput:BEMT_CalcOutput(node 50,
 blade 3):UA_CalcOutput:Mach number exceeds 0.3. Theory is invalid. This warning will not be
 repeated though the condition may persist.

Generator speed:    8.1 RPM, Pitch angle:   6.4 deg, Power: 16087.0 kW, Est. wind Speed:  13.1 m/s

  Total Real Time:       19.748 seconds
  Total CPU Time:        19.719 seconds
  Simulation CPU Time:   9.9531 seconds
  Simulated Time:        10.1 seconds
  Time Ratio (Sim/CPU):  1.0148

  OpenFAST terminated normally.

Runtime:        IEA15_0.fst = 19.78 s
Writing rotor performance text file: /home/rmccabe/WEIS/examples/02_control_opt/outputs/02_control_opt/IEA15_1_Cp_Ct_Cq.txt
Writing new controller parameter file parameter file: /home/rmccabe/WEIS/examples/02_control_opt/outputs/02_control_opt/IEA15_1_DISCON.IN.

 **************************************************************************************************
 OpenFAST

 Copyright (C) 2023 National Renewable Energy Laboratory
 Copyright (C) 2023 Envision Energy USA LTD

 This program is licensed under Apache License Version 2.0 and comes with ABSOLUTELY NO WARRANTY.
 See the "LICENSE" file distributed with this software for details.
 **************************************************************************************************

 OpenFAST-v1.1
 Compile Info:
  - Compiler: GCC version 9.4.0
  - Architecture: 64 bit
  - Precision: single
  - OpenMP: No
  - Date: Jan 30 2023
  - Time: 10:03:57
 Execution Info:
  - Date: 01/30/2023
  - Time: 10:37:17-0500

 OpenFAST input file heading:
     Generated with AeroElasticSE FAST driver

 Running ElastoDyn.
 Nodal outputs section of ElastoDyn input file not found or improperly formatted.
 Running AeroDyn.
 AD15 Nodal Outputs: Nodal output section of AeroDyn input file not found or improperly formatted.
 Skipping nodal outputs.
 Running InflowWind.

    Reading a 25x25 grid (297.47 m wide, 4.88281E-04 m to 297.47 m above ground) with a
    characteristic wind speed of 47.094 m/s. This full-field file was generated by TurbSim (from
    OpenFAST-v1.1) on 30-Jan-2023 at 10:34:47.

    Processed 204 time steps of 20-Hz full-field data (period of 10.2 seconds).
 Running HydroDyn.
  WARNING:   The random number generator in use differs from the original code provided by NREL.
  This pRNG uses 8 seeds instead of the 2 in the input file.
  Generating incident wave kinematics and current time history.
  Calculating second order sum frequency wave kinematics.
  Reading in WAMIT output with root name
  "/home/rmccabe/WEIS/examples/01_aeroelasticse/OpenFAST_models/IEA-15-240-RWT/IEA-15-240-RWT-UMain
  eSemi/HydroData/IEA-15-240-RWT-UMaineSemi".
  Computing radiation impulse response functions and wave diffraction forces.
 Running MoorDyn.
   MD_Init: Opening MoorDyn input file:  ./IEA15_1_MoorDyn.dat
   Creating mooring system.   3 fairleads, 3 anchors, 0 connects.
    Finalizing ICs using dynamic relaxation.

   t=60  FairTen 1: 2.43587E+06, 2.44120E+06, 2.43677E+06                                             Fairlead tensions did not converge within TMaxIC=60 seconds.
 Running ServoDyn.
 Running ServoDyn Interface for Bladed Controllers (using GNU Fortran for Linux).
 Using legacy Bladed DLL interface.

 FAST_InitializeAll:HydroDyn_Init:Waves_Init:VariousWaves_Init: The random number generator in use
 differs from the original code provided by NREL. This pRNG uses 8 seeds instead of the 2 in the
 HydroDyn input file.

  Time: 0 of 10.1 seconds.

------------------------------------------------------------------------------
Running ROSCO-v2.6.0
A wind turbine controller framework for public use in the scientific field
Developed in collaboration: National Renewable Energy Laboratory
                            Delft University of Technology, The Netherlands
------------------------------------------------------------------------------
Generator speed:    0.0 RPM, Pitch angle:  90.0 deg, Power:     0.0 kW, Est. wind Speed:  47.1 m/s
Generator speed:    0.2 RPM, Pitch angle:  70.0 deg, Power:     0.0 kW, Est. wind Speed: -20.6 m/s

  Total Real Time:       18.875 seconds
  Total CPU Time:        18.672 seconds
  Simulation CPU Time:   9.0312 seconds
  Simulated Time:        10.1 seconds
  Time Ratio (Sim/CPU):  1.1183

  OpenFAST terminated normally.

Runtime:        IEA15_1.fst = 18.91 s
Writing rotor performance text file: /home/rmccabe/WEIS/examples/02_control_opt/outputs/02_control_opt/IEA15_2_Cp_Ct_Cq.txt
Writing new controller parameter file parameter file: /home/rmccabe/WEIS/examples/02_control_opt/outputs/02_control_opt/IEA15_2_DISCON.IN.

 **************************************************************************************************
 OpenFAST

 Copyright (C) 2023 National Renewable Energy Laboratory
 Copyright (C) 2023 Envision Energy USA LTD

 This program is licensed under Apache License Version 2.0 and comes with ABSOLUTELY NO WARRANTY.
 See the "LICENSE" file distributed with this software for details.
 **************************************************************************************************

 OpenFAST-v1.1
 Compile Info:
  - Compiler: GCC version 9.4.0
  - Architecture: 64 bit
  - Precision: single
  - OpenMP: No
  - Date: Jan 30 2023
  - Time: 10:03:57
 Execution Info:
  - Date: 01/30/2023
  - Time: 10:37:37-0500

 OpenFAST input file heading:
     Generated with AeroElasticSE FAST driver

 Running ElastoDyn.
 Nodal outputs section of ElastoDyn input file not found or improperly formatted.
 Running AeroDyn.
 AD15 Nodal Outputs: Nodal output section of AeroDyn input file not found or improperly formatted.
 Skipping nodal outputs.
 Running InflowWind.

    Reading a 25x25 grid (297.47 m wide, 4.88281E-04 m to 297.47 m above ground) with a
    characteristic wind speed of 47.094 m/s. This full-field file was generated by TurbSim (from
    OpenFAST-v1.1) on 30-Jan-2023 at 10:34:55.

    Processed 204 time steps of 20-Hz full-field data (period of 10.2 seconds).
 Running HydroDyn.
  WARNING:   The random number generator in use differs from the original code provided by NREL.
  This pRNG uses 8 seeds instead of the 2 in the input file.
  Generating incident wave kinematics and current time history.
  Calculating second order sum frequency wave kinematics.
  Reading in WAMIT output with root name
  "/home/rmccabe/WEIS/examples/01_aeroelasticse/OpenFAST_models/IEA-15-240-RWT/IEA-15-240-RWT-UMain
  eSemi/HydroData/IEA-15-240-RWT-UMaineSemi".
  Computing radiation impulse response functions and wave diffraction forces.
 Running MoorDyn.
   MD_Init: Opening MoorDyn input file:  ./IEA15_2_MoorDyn.dat
   Creating mooring system.   3 fairleads, 3 anchors, 0 connects.
    Finalizing ICs using dynamic relaxation.

   t=60  FairTen 1: 2.43587E+06, 2.44120E+06, 2.43677E+06                                             Fairlead tensions did not converge within TMaxIC=60 seconds.
 Running ServoDyn.
 Running ServoDyn Interface for Bladed Controllers (using GNU Fortran for Linux).
 Using legacy Bladed DLL interface.

 FAST_InitializeAll:HydroDyn_Init:Waves_Init:VariousWaves_Init: The random number generator in use
 differs from the original code provided by NREL. This pRNG uses 8 seeds instead of the 2 in the
 HydroDyn input file.

  Time: 0 of 10.1 seconds.

------------------------------------------------------------------------------
Running ROSCO-v2.6.0
A wind turbine controller framework for public use in the scientific field
Developed in collaboration: National Renewable Energy Laboratory
                            Delft University of Technology, The Netherlands
------------------------------------------------------------------------------
Generator speed:    0.0 RPM, Pitch angle:  90.0 deg, Power:     0.0 kW, Est. wind Speed:  47.1 m/s
Generator speed:    0.2 RPM, Pitch angle:  70.0 deg, Power:     0.0 kW, Est. wind Speed: -20.6 m/s

  Total Real Time:       19.33 seconds
  Total CPU Time:        18.797 seconds
  Simulation CPU Time:   9.0625 seconds
  Simulated Time:        10.1 seconds
  Time Ratio (Sim/CPU):  1.1145

  OpenFAST terminated normally.

Runtime:        IEA15_2.fst = 19.36 s
WARNING: OpenFAST is run at a single wind speed. AEP cannot be estimated. Using average power instead.
########################################
Objectives
Turbine AEP: 0.0150820187 GWh
Blade Mass:  0.0000000000 kg
LCOE:        0.0000000000 USD/MWh
Tip Defl.:   0.0000000000 m
Pitch PI gain inputs: omega_pc[0] = 0.200, zeta_pc[0] = 1.000
Floating Feedback: Kp_float = -10.000, ptfm_freq = 0.200
IPC Ki1p = 0.000e+00
IPC Ki1p = 0.000e+00
DEL(TwrBsMyt): 224454450.0820847154 Nm
rotor_overspeed: 7.5571429760 %
########################################
Nonlinear constraints
{'aeroelastic.rotor_overspeed': array([0.07557143]),
 'sse_tune.tune_rosco.PC_Ki': array([-0.12552134]),
 'sse_tune.tune_rosco.PC_Kp': array([-1.20562287])}

Linear constraints
None

Objectives
{'aeroelastic.DEL_TwrBsMyt': array([224454.45008208])}

Driver debug print for iter coord: rank0:NLopt_LN_COBYLA|2
----------------------------------------------------------
Design Vars
{'tune_rosco_ivc.Kp_float': array([-10.]),
 'tune_rosco_ivc.omega_pc': array([0.3, 0.2]),
 'tune_rosco_ivc.ptfm_freq': array([0.2]),
 'tune_rosco_ivc.zeta_pc': array([1., 1.])}

Warning: turbine_costsse_2015.py : TowerCost2015.compute : No tower mass provided.  Assuming $0 for tower cost, LCOE will be artificially low.
RuntimeWarning: /home/rmccabe/WEIS/weis/control/dac.py:734
invalid value encountered in double_scalars-----------------------------------------------------------------------------
   Tuning a reference wind turbine controller using NREL's ROSCO toolbox
-----------------------------------------------------------------------------
Loading rotor performace data from text file: /home/rmccabe/WEIS/examples/01_aeroelasticse/OpenFAST_models/IEA-15-240-RWT/IEA-15-240-RWT/Cp_Ct_Cq.IEA15MW.txt
Writing rotor performance text file: /home/rmccabe/WEIS/examples/02_control_opt/outputs/02_control_opt/IEA15_0_Cp_Ct_Cq.txt
Writing new controller parameter file parameter file: /home/rmccabe/WEIS/examples/02_control_opt/outputs/02_control_opt/IEA15_0_DISCON.IN.

 **************************************************************************************************
 OpenFAST

 Copyright (C) 2023 National Renewable Energy Laboratory
 Copyright (C) 2023 Envision Energy USA LTD

 This program is licensed under Apache License Version 2.0 and comes with ABSOLUTELY NO WARRANTY.
 See the "LICENSE" file distributed with this software for details.
 **************************************************************************************************

 OpenFAST-v1.1
 Compile Info:
  - Compiler: GCC version 9.4.0
  - Architecture: 64 bit
  - Precision: single
  - OpenMP: No
  - Date: Jan 30 2023
  - Time: 10:03:57
 Execution Info:
  - Date: 01/30/2023
  - Time: 10:37:59-0500

 OpenFAST input file heading:
     Generated with AeroElasticSE FAST driver

 Running ElastoDyn.
 Nodal outputs section of ElastoDyn input file not found or improperly formatted.
 Running AeroDyn.
 AD15 Nodal Outputs: Nodal output section of AeroDyn input file not found or improperly formatted.
 Skipping nodal outputs.
 Warning: Turning off Unsteady Aerodynamics because UA parameters are not included in airfoil
 (airfoil has likely has constant polars). (node 1, blade 1)
 Warning: Turning off Unsteady Aerodynamics because UA parameters are not included in airfoil
 (airfoil has likely has constant polars). (node 2, blade 1)
 Warning: Turning off Unsteady Aerodynamics because UA parameters are not included in airfoil
 (airfoil has likely has constant polars). (node 1, blade 2)
 Warning: Turning off Unsteady Aerodynamics because UA parameters are not included in airfoil
 (airfoil has likely has constant polars). (node 2, blade 2)
 Warning: Turning off Unsteady Aerodynamics because UA parameters are not included in airfoil
 (airfoil has likely has constant polars). (node 1, blade 3)
 Warning: Turning off Unsteady Aerodynamics because UA parameters are not included in airfoil
 (airfoil has likely has constant polars). (node 2, blade 3)
 Running InflowWind.

    Reading a 25x25 grid (297.47 m wide, 4.88281E-04 m to 297.47 m above ground) with a
    characteristic wind speed of 16.149 m/s. This full-field file was generated by TurbSim (from
    OpenFAST-v1.1) on 30-Jan-2023 at 10:34:43.

    Processed 204 time steps of 20-Hz full-field data (period of 10.2 seconds).
 Running HydroDyn.
  WARNING:   The random number generator in use differs from the original code provided by NREL.
  This pRNG uses 8 seeds instead of the 2 in the input file.
  Generating incident wave kinematics and current time history.
  Calculating second order sum frequency wave kinematics.
  Reading in WAMIT output with root name
  "/home/rmccabe/WEIS/examples/01_aeroelasticse/OpenFAST_models/IEA-15-240-RWT/IEA-15-240-RWT-UMain
  eSemi/HydroData/IEA-15-240-RWT-UMaineSemi".
  Computing radiation impulse response functions and wave diffraction forces.
 Running MoorDyn.
   MD_Init: Opening MoorDyn input file:  ./IEA15_0_MoorDyn.dat
   Creating mooring system.   3 fairleads, 3 anchors, 0 connects.
    Finalizing ICs using dynamic relaxation.

   t=60  FairTen 1: 2.43587E+06, 2.44120E+06, 2.43677E+06                                             Fairlead tensions did not converge within TMaxIC=60 seconds.
 Running ServoDyn.
 Running ServoDyn Interface for Bladed Controllers (using GNU Fortran for Linux).
 Using legacy Bladed DLL interface.

 FAST_InitializeAll:HydroDyn_Init:Waves_Init:VariousWaves_Init: The random number generator in use
 differs from the original code provided by NREL. This pRNG uses 8 seeds instead of the 2 in the
 HydroDyn input file.

  Time: 0 of 10.1 seconds.

------------------------------------------------------------------------------
Running ROSCO-v2.6.0
A wind turbine controller framework for public use in the scientific field
Developed in collaboration: National Renewable Energy Laboratory
                            Delft University of Technology, The Netherlands
------------------------------------------------------------------------------
Generator speed:    7.6 RPM, Pitch angle:  15.0 deg, Power:     0.0 kW, Est. wind Speed:  15.9 m/s

 FAST_Solution:CalcOutputs_And_SolveForInputs:SolveOption2:RotCalcOutput:BEMT_CalcOutput(node 3,
 blade 3):UA_CalcOutput:UA_BlendSteady:Temporarily turning off UA due to high angle of attack or
 low relative velocity. This warning will not be repeated though the condition may persist.

 FAST_Solution:CalcOutputs_And_SolveForInputs:SolveOption2:RotCalcOutput:BEMT_CalcOutput(node 50,
 blade 3):UA_CalcOutput:Mach number exceeds 0.3. Theory is invalid. This warning will not be
 repeated though the condition may persist.

Generator speed:    8.1 RPM, Pitch angle:   6.4 deg, Power: 16087.0 kW, Est. wind Speed:  13.1 m/s

  Total Real Time:       20.541 seconds
  Total CPU Time:        20.125 seconds
  Simulation CPU Time:   10.234 seconds
  Simulated Time:        10.1 seconds
  Time Ratio (Sim/CPU):  0.98687

  OpenFAST terminated normally.

Runtime:        IEA15_0.fst = 20.57 s
Writing rotor performance text file: /home/rmccabe/WEIS/examples/02_control_opt/outputs/02_control_opt/IEA15_1_Cp_Ct_Cq.txt
Writing new controller parameter file parameter file: /home/rmccabe/WEIS/examples/02_control_opt/outputs/02_control_opt/IEA15_1_DISCON.IN.

 **************************************************************************************************
 OpenFAST

 Copyright (C) 2023 National Renewable Energy Laboratory
 Copyright (C) 2023 Envision Energy USA LTD

 This program is licensed under Apache License Version 2.0 and comes with ABSOLUTELY NO WARRANTY.
 See the "LICENSE" file distributed with this software for details.
 **************************************************************************************************

 OpenFAST-v1.1
 Compile Info:
  - Compiler: GCC version 9.4.0
  - Architecture: 64 bit
  - Precision: single
  - OpenMP: No
  - Date: Jan 30 2023
  - Time: 10:03:57
 Execution Info:
  - Date: 01/30/2023
  - Time: 10:38:20-0500

 OpenFAST input file heading:
     Generated with AeroElasticSE FAST driver

 Running ElastoDyn.
 Nodal outputs section of ElastoDyn input file not found or improperly formatted.
 Running AeroDyn.
 AD15 Nodal Outputs: Nodal output section of AeroDyn input file not found or improperly formatted.
 Skipping nodal outputs.
 Running InflowWind.

    Reading a 25x25 grid (297.47 m wide, 4.88281E-04 m to 297.47 m above ground) with a
    characteristic wind speed of 47.094 m/s. This full-field file was generated by TurbSim (from
    OpenFAST-v1.1) on 30-Jan-2023 at 10:34:47.

    Processed 204 time steps of 20-Hz full-field data (period of 10.2 seconds).
 Running HydroDyn.
  WARNING:   The random number generator in use differs from the original code provided by NREL.
  This pRNG uses 8 seeds instead of the 2 in the input file.
  Generating incident wave kinematics and current time history.
  Calculating second order sum frequency wave kinematics.
  Reading in WAMIT output with root name
  "/home/rmccabe/WEIS/examples/01_aeroelasticse/OpenFAST_models/IEA-15-240-RWT/IEA-15-240-RWT-UMain
  eSemi/HydroData/IEA-15-240-RWT-UMaineSemi".
  Computing radiation impulse response functions and wave diffraction forces.
 Running MoorDyn.
   MD_Init: Opening MoorDyn input file:  ./IEA15_1_MoorDyn.dat
   Creating mooring system.   3 fairleads, 3 anchors, 0 connects.
    Finalizing ICs using dynamic relaxation.

   t=60  FairTen 1: 2.43587E+06, 2.44120E+06, 2.43677E+06                                             Fairlead tensions did not converge within TMaxIC=60 seconds.
 Running ServoDyn.
 Running ServoDyn Interface for Bladed Controllers (using GNU Fortran for Linux).
 Using legacy Bladed DLL interface.

 FAST_InitializeAll:HydroDyn_Init:Waves_Init:VariousWaves_Init: The random number generator in use
 differs from the original code provided by NREL. This pRNG uses 8 seeds instead of the 2 in the
 HydroDyn input file.

  Time: 0 of 10.1 seconds.

------------------------------------------------------------------------------
Running ROSCO-v2.6.0
A wind turbine controller framework for public use in the scientific field
Developed in collaboration: National Renewable Energy Laboratory
                            Delft University of Technology, The Netherlands
------------------------------------------------------------------------------
Generator speed:    0.0 RPM, Pitch angle:  90.0 deg, Power:     0.0 kW, Est. wind Speed:  47.1 m/s
Generator speed:    0.2 RPM, Pitch angle:  70.0 deg, Power:     0.0 kW, Est. wind Speed: -20.6 m/s

  Total Real Time:       19.566 seconds
  Total CPU Time:        18.703 seconds
  Simulation CPU Time:   9.0938 seconds
  Simulated Time:        10.1 seconds
  Time Ratio (Sim/CPU):  1.1107

  OpenFAST terminated normally.

Runtime:        IEA15_1.fst = 19.60 s
Writing rotor performance text file: /home/rmccabe/WEIS/examples/02_control_opt/outputs/02_control_opt/IEA15_2_Cp_Ct_Cq.txt
Writing new controller parameter file parameter file: /home/rmccabe/WEIS/examples/02_control_opt/outputs/02_control_opt/IEA15_2_DISCON.IN.

 **************************************************************************************************
 OpenFAST

 Copyright (C) 2023 National Renewable Energy Laboratory
 Copyright (C) 2023 Envision Energy USA LTD

 This program is licensed under Apache License Version 2.0 and comes with ABSOLUTELY NO WARRANTY.
 See the "LICENSE" file distributed with this software for details.
 **************************************************************************************************

 OpenFAST-v1.1
 Compile Info:
  - Compiler: GCC version 9.4.0
  - Architecture: 64 bit
  - Precision: single
  - OpenMP: No
  - Date: Jan 30 2023
  - Time: 10:03:57
 Execution Info:
  - Date: 01/30/2023
  - Time: 10:38:41-0500

 OpenFAST input file heading:
     Generated with AeroElasticSE FAST driver

 Running ElastoDyn.
 Nodal outputs section of ElastoDyn input file not found or improperly formatted.
 Running AeroDyn.
 AD15 Nodal Outputs: Nodal output section of AeroDyn input file not found or improperly formatted.
 Skipping nodal outputs.
 Running InflowWind.

    Reading a 25x25 grid (297.47 m wide, 4.88281E-04 m to 297.47 m above ground) with a
    characteristic wind speed of 47.094 m/s. This full-field file was generated by TurbSim (from
    OpenFAST-v1.1) on 30-Jan-2023 at 10:34:55.

    Processed 204 time steps of 20-Hz full-field data (period of 10.2 seconds).
 Running HydroDyn.
  WARNING:   The random number generator in use differs from the original code provided by NREL.
  This pRNG uses 8 seeds instead of the 2 in the input file.
  Generating incident wave kinematics and current time history.
  Calculating second order sum frequency wave kinematics.
  Reading in WAMIT output with root name
  "/home/rmccabe/WEIS/examples/01_aeroelasticse/OpenFAST_models/IEA-15-240-RWT/IEA-15-240-RWT-UMain
  eSemi/HydroData/IEA-15-240-RWT-UMaineSemi".
  Computing radiation impulse response functions and wave diffraction forces.
 Running MoorDyn.
   MD_Init: Opening MoorDyn input file:  ./IEA15_2_MoorDyn.dat
   Creating mooring system.   3 fairleads, 3 anchors, 0 connects.
    Finalizing ICs using dynamic relaxation.

   t=60  FairTen 1: 2.43587E+06, 2.44120E+06, 2.43677E+06                                             Fairlead tensions did not converge within TMaxIC=60 seconds.
 Running ServoDyn.
 Running ServoDyn Interface for Bladed Controllers (using GNU Fortran for Linux).
 Using legacy Bladed DLL interface.

 FAST_InitializeAll:HydroDyn_Init:Waves_Init:VariousWaves_Init: The random number generator in use
 differs from the original code provided by NREL. This pRNG uses 8 seeds instead of the 2 in the
 HydroDyn input file.

  Time: 0 of 10.1 seconds.

------------------------------------------------------------------------------
Running ROSCO-v2.6.0
A wind turbine controller framework for public use in the scientific field
Developed in collaboration: National Renewable Energy Laboratory
                            Delft University of Technology, The Netherlands
------------------------------------------------------------------------------
Generator speed:    0.0 RPM, Pitch angle:  90.0 deg, Power:     0.0 kW, Est. wind Speed:  47.1 m/s
Generator speed:    0.2 RPM, Pitch angle:  70.0 deg, Power:     0.0 kW, Est. wind Speed: -20.6 m/s

  Total Real Time:       19.507 seconds
  Total CPU Time:        18.844 seconds
  Simulation CPU Time:   9.0625 seconds
  Simulated Time:        10.1 seconds
  Time Ratio (Sim/CPU):  1.1145

  OpenFAST terminated normally.

Runtime:        IEA15_2.fst = 19.54 s
WARNING: OpenFAST is run at a single wind speed. AEP cannot be estimated. Using average power instead.
########################################
Objectives
Turbine AEP: 0.0150820187 GWh
Blade Mass:  0.0000000000 kg
LCOE:        0.0000000000 USD/MWh
Tip Defl.:   0.0000000000 m
Pitch PI gain inputs: omega_pc[0] = 0.300, zeta_pc[0] = 1.000
Floating Feedback: Kp_float = -10.000, ptfm_freq = 0.200
IPC Ki1p = 0.000e+00
IPC Ki1p = 0.000e+00
DEL(TwrBsMyt): 224454450.0820847154 Nm
rotor_overspeed: 7.5571429760 %
########################################
Nonlinear constraints
{'aeroelastic.rotor_overspeed': array([0.07557143]),
 'sse_tune.tune_rosco.PC_Ki': array([-0.28242301]),
 'sse_tune.tune_rosco.PC_Kp': array([-1.83322957])}

Linear constraints
None

Objectives
{'aeroelastic.DEL_TwrBsMyt': array([224454.45008208])}

/home/rmccabe/anaconda3/envs/weis-env/lib/python3.10/site-packages/openmdao/visualization/opt_report/opt_report.py:119: DriverWarning:The optimizer report is not applicable for Driver type 'NLoptDriver', which does not support optimization
Run time: 265.874547

Expected behavior

Runs without warnings with a nonzero LCOE result.

Code versions

List versions only if relevant

dzalkind commented 1 year ago

Hi,

Despite the abundance of warnings, this example ran to completion properly.

LCOE is expected to be zero in this example because most of the WISDEM system engineering modules are disabled, so the cost of blades, tower, etc. are not calculated. The 02_examples use the OpenFAST model only, rather than generating one using WISDEM models.

The only warnings that I would be concerned about are the ones about MPI and the N2 diagrams. If you want to parallelize the runs using MPI, you might need to make sure mpi4py is installed properly. I'm not sure about the N2 diagrams, though I suspect your version of openmdao might need to be depreciated. If you don't want N2 diagrams at this point, you can ignore them.

I hope this helps.

Best, Dan

rebeccamccabe commented 1 year ago

Hi Dan, Thanks so much for your quick response! I will look into mpi4py installation and openmdao versions. Is there a place I should look to understand more clearly what each example is doing? I didn't see much in the documentation. Thanks!

dzalkind commented 7 months ago

Hi Rebecca,

We have recently kicked off a new phase of funded development, and added documentation to this example here: https://github.com/WISDEM/WEIS/tree/main/examples/02_control_opt

I apologize for the long delay in response. We are actively developing WEIS again, and should be more responsive in the future.

Best, Dan