In-Progress Multi-objective Optimization branch

docs/ftw_overview.rst

@@ -0,0 +1,13 @@
+.. _section-ftw:
+
+====================================
+Formulation Toolbox for WEIS (FTW)
+====================================
+
+
+The Formulation Toolbox for WEIS (FTW) provides a set of tools and methodologies for formulating and analyzing design optimizations.
+
+
+--------------------------
+Introduction
+--------------------------

docs/index.rst

@@ -73,6 +73,15 @@ WEIS Visualization Application
 
    how_weis_viz_works
 
+
+Formulation Toolbox for WEIS
+=============================
+
+.. toctree::
+   :maxdepth: 2
+
+   ftw_overview
+
 Developer Documentation
 =========================
 

examples/00_setup/OpenFAST_models/update_weis_discons.py

@@ -22,3 +22,4 @@
 
     # Make discons    
     update_discons(discon_map)
+

examples/02_run_openfast_cases/analysis_options_filters.yaml

@@ -0,0 +1,42 @@
+general:
+    folder_output: outputs/03_filter_study/6_DISCON_only
+    fname_output: iea15mw
+
+merit_figure: DEL_TwrBsMyt  # Merit figure of the optimization problem. The options are 'AEP' - 'LCOE' - 'Cp' - 'blade_mass' - 'blade_tip_deflection'
+
+design_variables:
+  control:
+    discon:
+      - name: F_FlCornerFreq
+        start: [0.250000, 0.500000]
+        min: [0.001, 0.1]
+        max: [1.0, 3.0]
+
+
+constraints:
+  control:
+    rotor_overspeed:
+      flag: True
+      min: 0.0
+      max: 0.2
+
+driver:
+  optimization:
+    flag: True
+    tol: 1.e-2            # Optimality tolerance
+    max_major_iter: 2     # Maximum number of major design iterations (SNOPT)
+    max_minor_iter: 100   # Maximum number of minor design iterations (SNOPT)
+    max_iter: 2          # Maximum number of iterations (SLSQP)
+    solver: LN_COBYLA         # Optimization solver. Other options are 'SLSQP' - 'CONMIN'
+    step_size: 1.e-3      # Step size for finite differencing
+    form: forward         # Finite differencing mode, either forward or central
+  design_of_experiments:
+    flag: False            # Flag to enable design of experiments
+    run_parallel: False    # Flag to run using parallel processing
+    generator: LatinHypercube    # Type of input generator. (Uniform)
+    num_samples: 8        # number of samples for (Uniform only)
+    criterion: center
+
+recorder:
+    flag: True             # Flag to activate OpenMDAO recorder
+    file_name: log_opt.sql  # Name of OpenMDAO recorder

examples/02_run_openfast_cases/modeling_options_no_tune.yaml

@@ -0,0 +1,116 @@
+General:
+    verbosity: False  # When set to True, the code prints to screen many infos
+    openfast_configuration:
+        use_exe: True
+
+WISDEM:
+    RotorSE:
+        flag: False
+        spar_cap_ss: Spar_Cap_SS
+        spar_cap_ps: Spar_Cap_PS
+        te_ss: TE_reinforcement_SS
+        te_ps: TE_reinforcement_PS
+    TowerSE:
+        flag: False
+    DriveSE:
+        flag: False
+    FloatingSE:
+        flag: False
+    BOS:
+        flag: False
+
+Level3: # Options for WEIS fidelity level 3 = nonlinear time domain
+    flag: True
+    from_openfast: True
+    regulation_trajectory: IEA-15-240-RWT-outputs.yaml
+    openfast_file: IEA-15-240-RWT-UMaineSemi.fst
+    openfast_dir: ../01_aeroelasticse/OpenFAST_models/IEA-15-240-RWT/IEA-15-240-RWT-UMaineSemi
+    AeroDyn:
+        WakeMod: 2
+
+ROSCO:
+    flag: False
+    # tuning_yaml: ../01_aeroelasticse/OpenFAST_models/IEA-15-240-RWT/IEA-15-240-RWT-UMaineSemi/IEA15MW-UMaineSemi.yaml
+    # zeta_pc: 1.23
+    # omega_pc: 0.198
+    # U_pc: [12]
+    # Fl_Mode: 0
+    # DISCON:
+    #     # F_FlCornerFreq: [0.213, 1.]     # filt1
+    #     # F_FlCornerFreq: [0.25, .5]  # filt2
+    #     F_FlCornerFreq: [0.179, .638]   # filt3
+    #     # F_FlCornerFreq: [0.285, .97]    # filt4
+
+
+DLC_driver:
+    metocean_conditions:
+        wind_speed: [1.,3.,5.,7.,9.,11.,13.,15.,17.,19.,21.,23.,25.,27.,29.]
+        wave_height_NSS: [0.84,0.84,0.87,0.99,1.15,1.34,1.58,1.82,2.08,2.34,2.66,2.98,3.28,3.77,3.94]
+        wave_period_NSS: [8.3,8.3,8.3,7.7,7.1,6.3,6.1,6.2,6.2,6.7,7.1,7.1,7.7,7.7,7.7]
+        wave_height_fatigue: [0.84,0.84,0.87,0.99,1.15,1.34,1.58,1.82,2.08,2.34,2.66,2.98,3.28,3.77,3.94]
+        wave_period_fatigue: [8.3,8.3,8.3,7.7,7.1,6.3,6.1,6.2,6.2,6.7,7.1,7.1,7.7,7.7,7.7]
+        wave_height_SSS: [9.7,9.7,9.7,9.7,9.7,9.7,9.7,9.7,9.7,9.7,9.7,9.7,9.7,9.7,9.7]
+        wave_period_SSS: [13.6,13.6,13.6,13.6,13.6,13.6,13.6,13.6,13.6,13.6,13.6,13.6,13.6,13.6,13.6]
+        wave_height1: 5.9
+        wave_period1: 11.2
+        wave_height50: 9.7
+        wave_period50: 13.6
+    DLCs:
+        - DLC: "1.1"
+          n_seeds: 1
+          wind_speed: [12]
+          turbulent_wind:
+              HubHt: 148.73499999999999
+              GridHeight: 297.469 
+              GridWidth: 297.469 
+              RefHt: 148.73499999999999
+          transient_time: 0
+          analysis_time: 1.
+        # - DLC: "1.6"
+        #   n_seeds: 1
+        #   wind_speed: [12]
+        #   transient_time: 0
+        #   analysis_time: 1.
+        # - DLC: "1.3"
+        #   n_seeds: 6
+        # - DLC: "1.4"
+        #   wind_speed: [8.6, 10.6, 12.6]
+        #   n_seeds: 2
+        #   n_azimuth: 3
+        # - DLC: "1.5"
+        # - DLC: "1.6"
+        #   n_seeds: 1
+        #   wind_speed: [12]
+        #   yaw_misalign: [-8,8]
+        #   analysis_time: 10.
+        #   transient_time: 0.
+        # - DLC: "1.6"
+        #   n_seeds: 1
+        #   wind_speed: [12]
+        #   wave_heading: [30]
+        #   analysis_time: 10.
+        #   transient_time: 0.
+        # - DLC: "6.3"
+        #   n_seeds: 6
+        # - DLC: "5.1"
+        #   wind_speed: [12]
+        #   n_seeds: 1
+        #   n_azimuth: 3
+        #   analysis_time: 20.
+        #   shutdown_time: 10.
+        #   transient_time: 0.0
+        # - DLC: "6.5"
+        #   yaw_misalign: [-8.,8.]
+        #   wave_height: [12.57]
+        #   wave_period: [16.0]
+        #   transient_time: 0.0
+        #   wave_gamma: [1.2]
+        #   analysis_time: 10.
+        #   turbulent_wind:
+        #       HubHt: 140
+        #       GridHeight: 220
+        #       GridWidth: 220
+        #       URef: 54.8339
+        #       PLExp: 0.12
+
+

examples/02_run_openfast_cases/weis_driver_filters.py

@@ -0,0 +1,24 @@
+import os
+import time
+import sys
+
+from weis.glue_code.runWEIS     import run_weis
+from openmdao.utils.mpi  import MPI
+
+## File management
+run_dir                 = os.path.dirname( os.path.realpath(__file__) )
+fname_wt_input          = run_dir + os.sep + 'IEA-15-240-RWT.yaml'
+fname_modeling_options  = run_dir + os.sep + 'modeling_options_no_tune.yaml'
+fname_analysis_options  = run_dir + os.sep + 'analysis_options_filters.yaml'
+
+
+tt = time.time()
+wt_opt, modeling_options, opt_options = run_weis(fname_wt_input, fname_modeling_options, fname_analysis_options)
+
+if MPI:
+    rank = MPI.COMM_WORLD.Get_rank()
+else:
+    rank = 0
+if rank == 0:
+    print("Run time: %f"%(time.time()-tt))
+    sys.stdout.flush()